Simultaneous identification of a single pollution point-source location and contamination time under known flow field conditions

A theoretical framework is presented that allows direct identification of a single point-source pollution location and time in heterogeneous multidimensional systems under known flow field conditions. Based on the concept of the transfer function theory, it is shown that an observed pollution plume contains all the necessary information to predict the concentration at the unknown pollution source when a reversed flow field transport simulation is performed. This target concentration C₀ is obtained from a quadratic integral of the observed pollution plume itself. Backwards simulation of the pollution plume leads to shrinkage of the C₀-contour due to dispersion. When the C₀-contour reduces to a singular point, i.e. becomes a concentration maximum, the position of the pollution source is identified and the backward simulation time indicates the time elapsed since the contaminant release. The theoretical basis of the method is first developed for the ideal case that the pollution plume is entirely known and is illustrated using a synthetic heterogeneous 2D example where all the hydro-dispersive parameters are known. The same example is then used to illustrate the procedure for a more realistic case, i.e. where only few observation points exist.     

Direct simulation of solute recycling in irrigated areas

Solute recycling from irrigation can be described as the process that occurs when the salt load that is extracted from irrigation wells and distributed on the fields is returned to the groundwater below irrigated surfaces by deep percolation. Unless the salt load leaves the system by means of drains or surface runoff, transfer to the groundwater will take place, sooner or later. This can lead to solute accumulation and thus to groundwater degradation, particularly in areas where extraction rates exceed infiltration rates (semi-arid and arid regions). Thus, considerable errors can occur in a predictive solute mass budget if the recycling process is not accounted for in the calculation. A method is proposed which allows direct simulation of solute recycling. The transient solute response at an extraction well is shown to be a superposition of solute mass flux contributions from n recycling cycles and is described as a function of the travel time distribution between a recycling point and a well. This leads to an expression for a transient ‘recycling source’ term in the advection–dispersion equation, which generates the effect of solute recycling. At long times, the ‘recycling source’ is a function of the local capture probability of the irrigation well and the solute mass flux captured by the well from the boundaries. The predicted concentration distribution at steady state reflects the maximum spatial concentration distribution in response to solute recycling and can thus be considered as the solute recycling potential or vulnerability of the entire domain for a given hydraulic setting and exploitation scheme. Simulation of the solute recycling potential is computationally undemanding and can therefore, for instance, be used for optimisation purposes. Also, the proposed method allows transient simulation of solute recycling with any standard flow and transport code.     

Sequence and style of major post-nappe structures, Simplon—Pennine Alps

In the Upper Pennine nappe complex of the Simplon—Pennine Alps (Switzerland and Italy), at least three phases of major post-nappe folding (in places associated with thrusting) can be distinguished. These are superimposed on an earlier-formed, partly chaotic, complex of tectonic units, including the Bernhard and Monte Rosa continental flakes and the Zermatt—Saas and Antrona ophiolite complexes. The earliest post-nappe folds were essentially isoclinal throughout the whole region and were accompanied by a strong schistosity which is the main foliation in most areas. Later, two successive phases of back-folding led to the present overall structure. Both phases typically show rapid variations in style from open folds lacking axial planar schistosity to very tight structures with complete foliation transposition. This has been demonstrated by systematically mapping the major axial traces over the whole region. Successively removing the major structures in reverse order shows that the ophiolite complexes were originally part of a continuous unit marking an important suture between the Bernhard and Monte Rosa nappes.     

Rubidium‐strontium geochronology of the encounter bay granite and adjacent Metasedimentary rocks,South Australia

Rubidium‐strontium and strontium isotope data for eight whole‐rock samples of granite varieties from the Encounter Bay area, South Australia, yield an isochron age of 487 ± 37 m.y. Two specimens of albitised granite, formed as a result of late‐stage metasomatic alteration of original megacrystic granite, conform to this isochron. These data support a genetic relation between granites and late‐stage metasomatic alteration as suspected from field, petrographical and geochemical studies. Eight samples from contiguous Kanmantoo Group metasedimentary rocks have an isochron age of 487 ± 60 m.y. Thus this metamorphic event is coincident with emplacement of the Encounter Bay Granite.

The initial Sr⁸⁷Sr⁸⁶ ratio for the Encounter Bay Granite (0.719) is significantly higher than initial ratios for the Palmer (0.709) and Anabama (0.705) Granites from the same region and can be attributed to either remobilisation or incorporation of strontium from older crustal material in the intrusion. The apparent initial Sr⁸⁷/Sr⁸⁶ ratio for the Kanmantoo Group metasedimentary rocks (0.722) can not be distinguished from that for the Encounter Bay Granite within the analytical uncertainties. Compatability of ages and high initial Sr⁸⁷Sr⁸⁶ ratios suggest that the granites formed by remobilisation of associated crustal rock.     

Petrology and mineralogy of ‘laterites’ in southern and eastern Australia and southern Africa

Field observations of ‘laterites’ in southern and eastern Australia and in southern Africa reveal a variety of ferruginous horizons and crusts referred to herein as ferricretes. Their geomorphic and stratigraphic relationships with bedrock, sediments and soils indicate formation throughout long intervals of geological time in landscapes which are also characterised by zones of bleached and iron-mottled materials. There does not appear to be a genetic relationship between the ferricretes and the weathered zones in the sense of the so-called ‘laterite profile’. Many of the ferricretes form part of existing soil profiles.

Petrographic studies of a variety of ferricretes have identified three broad categories: (a) ferruginised bedrock; (b) Fe-impregnated and -indurated sediments, including sands, clays and organic sediments; and (c) ferricretes of complex sedimentary and pedogenic origin. Type-(a) and -(b) ferricretes characteristically have simple fabrics, often with single-generation, secondary Fe-oxides. Type-(c) ferricretes have complex fabrics, with many generations of hematite, goethite and in some variants, gibbsite, in the matrix and in ferruginous clasts and pisoliths. Maghemite is a common constituent of the pisoliths. The characteristics of ironstone gravelly duplex soils, which are common in the contemporary landscapes, provide the framework for a model involving multiple stages in the development of these ferricretes.

The origins of the various secondary oxide minerals in ferricretes are assessed on the basis of knowledge about the formation of these minerals in pedogenic environments. Examples are given of the intricate patterns of distribution of the minerals in thin section from which definitive data may be obtained on environmental conditions for integration with field-based geomorphic studies.     

Modelling Discharge Rates and Ground Settlement Induced by Tunnel Excavation

Interception of aquifers by tunnel excavation results in water inflow and leads to drawdown of the water table which may induce ground settlement. In this work, analytical and numerical models are presented which specifically address these groundwater related processes in tunnel excavation. These developed models are compared and their performance as predictive tools is evaluated. Firstly, the water inflow in deep tunnels is treated. It is shown that introducing a reduction factor accounting for the effect of effective stress on hydrodynamic parameters avoids overestimation. This effect can be considered in numerical models using effective stress-dependent parameters. Then, quantification of ground settlement is addressed by a transient analytical solution. These solutions are then successfully applied to the data obtained during the excavation of the La Praz exploratory tunnel in the Western Alps (France), validating their usefulness as predictive tools.     

Pre‐ to syn‐tectonic emplacement of early Palaeozoic granites in southeastern South Australia

Stratigraphic and structural observations indicate that the Encounter Bay Granites concordantly intruded the youngest formations of the Kanmantoo Group in the Mount Lofty Ranges metamorphic belt prior to the culmination of the first phase of folding and associated schistosity development recorded during the early Palaeozoic Delamerian Orogeny. Metamorphic textures in the metasediments of the Kanmantoo Group suggest that cordierite crystallized locally near the granites prior to and during the F ₁ folding, whereas andalusite crystallized on a regional scale during the F ₁ folding and in the post‐F ₁ and pre‐F ₂ static phase.

Rb‐Sr isotope data for total‐rock, feldspar, and muscovite samples of the meta‐sediment‐contaminated border facies and the uncontaminated inner facies of the Encounter Bay Granites indicate that the granites were emplaced between 515 ± 8 m.y. and 506 ± 6 m.y. ago in the Late Cambrian epoch. Rb‐Sr and K‐Ar data for biotite from the granites record variable radiogenic Sr loss until about 469 m.y. ago and comparatively uniform radiogenic Ar loss until 460–475 m.y. ago. Rb‐Sr data for Kanmantoo Group metasediments and a metamorphic pegmatite indicate crystallization ages between 459–463 m.y. ago. Thus the regional andalusite‐grade temperatures and pressures, which appear responsible for the leakage of radiogenic Sr and Ar from biotite in the granites and the redistribution of Rb and Sr in the metasediments, seem to have persisted for some 50 m.y. after emplacement of the granites until the Early Ordovician epoch. There is evidence for further leakage of Sr and Ar from biotite in deformed granites from the margins of the intrusion more than 50 m.y. afterwards in the Late Silurian or Early Devonian, possibly during the F ₂ folding.

Geological observations and radiometric data for other granitic rocks in southeastern South Australia, including the Palmer Granite, are consistent with this structural and metamorphic history of the Encounter Bay region.     

Origin of aragonite in a continental saline seep,north Queensland

Aragonite, low‐magnesian calcite, gypsum and halite were identified by X‐ray diffraction and electron microbeam techniques in mineral precipitates near a salt seep 50 km southwest of Charters Towers in north Queensland. The chemistry of water from the creek and from the groundwater at the salt seep shows that Mg:Ca ratios are greater than or equal to 1.5 throughout the year. The formation of halite and gypsum is due to evaporative concentration of the water at the seep and that of the carbonates, in particular aragonite, is probably due to a combination of evaporation and photosynthetic activity by diatoms.