Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

Elm bark beetle in Holocene peat deposits and the northwest European elm decline

The elm decline of 5000 ¹⁴C yr ago has been the most widely discussed phenomenon in post‐glacial vegetation history. This pan‐European reduction of elm populations, echoed in the decimation of elmwoods in Europe during the twentieth century, has attracted a series of interrelated hypotheses involving climate change, human activity, disease and soil deterioration. The elm bark beetle (Scolytus scolytus L.) is an essential component of disease explanations. We present evidence for the presence of the beetle over a prolonged period (ca. 7950–4910 yr BP [8800–5660 cal. yr BP]) from a lowland raised mire deposit in northeast Scotland, with its final appearance at this site, and the first and only appearance in another mire of a single scolytid find, around the time of the elm decline. The subfossil S. scolytus finds are not only the first from Scotland, but they also represent the most comprehensive sequence of finds anywhere. Copyright © 2004 John Wiley & Sons, Ltd.     

A reanalysis of the light curve of RT Persei with the W-D code

The UBV photometric observations of RT Per, from Sanwal and Chaubey (1981), were analyzed by the Wilson and Devinney code (1971). The light curves include reflection effects that for the first time has been suggested by Dugan (1911). RT Per has a semi-detached configuration where the lower-mass component is in contact with its respective Roche surface. The higher-mass component very nearly fills its Roche lobe. It has the characteristic of an Algol type system. The absolute dimensions for the primary and secondary of this system were calculated from its spectral types and by combining the photometric solution with inferred component radial velocities (Lu, 1990).     

Urban geology: Glasgow's hidden industrial heritage

The growth of Glasgow from a small town into a heavily industrialised conurbation depended greatly on its local geological riches. Extensive mining and quarrying of a range of minerals took place from the 18th century onwards. The early underground stoop and room (pillar and stall) workings, and the backfilled quarries, together with variably consolidated natural superficial deposits, have bequeathed to the city a heritage of unwanted engineering problems which cannot easily be quantified. Recent work by the British Geological Survey (BGS) illustrates the scale of both geological and man-made problems.     

Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

Age and early metamorphic history of the Sanbagawa belt: Lu–Hf and P–T constraints from the Western Iratsu eclogite

Two distinct age estimates for eclogite-facies metamorphism in the Sanbagawa belt have been proposed: (i) c.  120–110 Ma based on a zircon SHRIMP age for the Western Iratsu unit and (ii) c.  88–89 Ma based on a garnet–omphacite Lu–Hf isochron age from the Seba and Kotsu eclogite units. Despite the contrasting estimates of formation ages, petrological studies suggest the formation conditions of the Western Iratsu unit are indistinguishable from those of the other two units—all ∼20 kbar and 600–650 °C. Studies of the associated geological structures suggest the Seba and Western Iratsu units are parts of a larger semi-continuous eclogite unit. A combination of geochronological and petrological studies for the Western Iratsu eclogite offers a resolution to this discrepancy in age estimates. New Lu–Hf dating for the Western Iratsu eclogite yields an age of 115.9 ± 0.5 Ma that is compatible with the zircon SHRIMP age. However, petrological studies show that there was significant garnet growth in the Western Iratsu eclogite before eclogite facies metamorphism, and the early core growth is associated with a strong concentration of Lu. Pre-eclogite facies garnet (Grt1) includes epidote–amphibolite facies parageneses equilibrated at 550–650 °C and ∼10 kbar, and this is overgrown by prograde eclogite facies garnet (Grt2). The Lu–Hf age of c.  116 Ma is strongly skewed to the isotopic composition of Grt1 and is interpreted to reflect the age of the pre-eclogite phase. The considerable time gap ( c.  27 Myr) between the two Lu–Hf ages suggests they may be related to separate tectonic events or distinct phases in the evolution of the Sanbagawa subduction zone.