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.     

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).     

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.     

Zircon U–Pb age, trace element and Hf isotope composition of Kongling terrane in the Yangtze Craton: refining the timing of Palaeoproterozoic high-grade metamorphism

U–Pb age, trace element and Hf isotope compositions of zircon were analysed for a metasedimentary rock and two amphibolites from the Kongling terrane in the northern part of the Yangtze Craton. The zircon shows distinct morphological and chemical characteristics. Most zircon in an amphibolite shows oscillatory zoning, high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, high formation temperature, high trace element contents, clear negative Eu anomaly, as well as HREE-enriched patterns, suggesting that it is igneous. The zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2857 ± 8 Ma, representing the age of the magmatic protolith. The zircon in the other two samples is metamorphic. It has low Th/U ratios, low trace element concentrations, variable HREE contents (33.8 ≥ Lu_N≥2213; 14.7 ≤ Lu_N/Sm_N ≤ 354) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.000030–0.001168). The data indicate that the zircon formed in the presence of garnet and under upper amphibolite facies conditions. The metamorphic zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2010 ± 13 Ma. These results combined with previously obtained Palaeoproterozoic metamorphic ages suggest a c. 2.0 Ga Palaeoproterozoic collisional event in the Yangtze Craton, which may result from the assembly of the supercontinent Columbia. The zircon in two samples yields weighted mean two-stage Hf model ( T _DM2) ages of 3217 ± 110 and 2943 ± 50 Ma, respectively, indicating that their protoliths were mainly derived from Archean crust.     

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.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

Green Lake Landslide and other giant and very large postglacial landslides in Fiordland,New Zealand

Green Lake Landslide is an ancient giant rock slide in gneiss and granodiorite located in the deeply glaciated Fiordland region of New Zealand. The landslide covers an area of 45 km² and has a volume of about 27 km³. It is believed to be New Zealand's largest landslide, and possibly the largest landslide of its type on Earth. It is one of 39 known very large (10⁶–10⁷ m³) and giant (≥10⁸ m³) postglacial landslides in Fiordland discussed in the paper. Green Lake Landslide resulted in the collapse of a 9 km segment of the southern Hunter Mountains. Slide debris moved up to 2.5 km laterally and 700 m vertically, and formed a landslide dam about 800 m high, impounding a lake about 11 km long that was eventually infilled with sediments. Geomorphic evidence supported by radiocarbon dating indicates that Green Lake Landslide probably occurred 12 000–13 000 years ago, near the end of the last (Otira) glaciation. The landslide is described, and its geomorphic significance, age, failure mechanism, cause, and relevance in the region are discussed, in relation to other large landslides and recent earthquake-induced landslides in Fiordland. The slope failure occurred on a low-angle fault zone undercut by glacial erosion, and was probably triggered by strong shaking (MM IX–X) associated with a large (≥ M 7.5–8) earthquake, on the Alpine Fault c. 80 km to the northwest. Geology was a major factor that controlled the style and size of Green Lake landslide, and in that respect it is significantly different from most other gigantic landslides. Future large earthquakes on the Alpine Fault in Fiordland are likely to trigger more very large and giant landslides across the region, causing ground damage and devastation on a scale that has not occurred during the last 160 years, with potentially disastrous effects on towns, tourist centres, roads, and infrastructure. The probability of such an event occurring within the next 50 years may be as high as 45%.     

Isotope fractionation of neon during stepheating extraction?: a comment on ‘Re‐interpretation of the existence of a primitive plume under Australia based on neon isotope fractionation during step heating’ by Gautheron and Moreira (2003)

Using a Rayleigh distillation fractionation model, we calculate that the maximum isotope fractionation potentially achievable is less than 5% during the early stages of gas release from a sample. Our calculation corrects the erroneous conclusions of Gautheron and Moreira (2003), who re‐interpreted the plume‐like neon isotopic compositions found in metasomatic apatite from a south‐eastern Australian xenolith (Matsumoto et al., 1997) to be the result of Rayleigh‐type isotope fractionation of originally MORB‐type neon during stepheating gas extraction. We stress that the modelling of neon isotopic fractionation by Gautheron and Moreira (2003) is incorrect, and that the finding of a plume‐like neon isotopic composition in the apatite by Matsumoto et al. (1997) remains a quite valid and robust conclusion.