A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

Retrograde fluids in the Archean Shawmere anorthosite, Kapuskasing Structural Zone, Ontario, Canada

The Archean Shawmere anorthosite lies within the granulite facies portion of the Kapuskasing Structural Zone (KSZ), Ontario, and is crosscut by numerous linear alteration veins containing calcite + quartz ± dolomite ± zoisite ± clinozoisite ± margarite ±paragonite ± chlorite. These veins roughly parallel the trend of the Ivanhoe Lake Cataclastic Zone. Equilibria involving clinozoisite + margarite + quartz ± calcite ± plagioclase show that the vein minerals were stable at T < 600 °C, XCO₂ < 0.4 at P ≈ 6 kbar. The stabilities of margarite and paragonite in equilibrium with quartz are also consistent with T < 600 °C and XCO₂ < 0.4 at 6 kbar. Additional assemblages consisting of calcite + clinochlore + quartz + talc + margarite indicate T < 500 °C with XCO₂ > 0.9. Thus, vein formation, while clearly retrograde, spanned a range of temperatures, and fluid compositions evolved from H₂O-rich to CO₂-rich. The calcite in the retrograde veins has δ¹⁸O values that range from 8.4 to 11.2‰ (average = +9.7 ± 0.9‰) and δ¹³C values that range from −3.9 to −1.6‰ (average = −3.1 ± 0.6‰). These values indicate that the fluids from which calcite precipitated underwent extensive exchange with the anorthosite and other crustal lithologies. The fluids may have been initially derived either from devolatilization of metamorphic rocks or crystallization of igneous rocks in the adjacent Abitibi subprovince. Vein quartz contains CO₂-rich fluid inclusions (final melting T = −57.0 to −58.7 °C) that range in size from 5 to 17 μm. Measured homogenization temperatures (T h) range from −44.0 to 14.5 °C, however for most inclusions (46 of S1), T h = −44.0 to −21.1 °C (ρCO₂ ≈ 1.13 to 1.05 g/cm³). At 400 to 600 °C, these densities correspond to pressures of 3.5 to 7 kbar, which is the best estimate of pressures of vein formation. It has been argued that some high density CO₂-rich fluid inclusions found in the KSZ were formed during peak metamorphism and thus document the presence of a CO₂-rich fluid during peak granulite facies metamorphism (Rudnick et al. 1984). The association of high density CO₂-rich fluid inclusions with clearly retrograde veins documents the formation of similar composition and density inclusions after the peak of metamorphism. Thus, the coincidence of entrapment pressures calculated from fluid inclusion density measurements with peak metamorphic pressures alone should not be considered strong evidence for peak metamorphic inclusion entrapment. All fluid inclusion results are consistent with an initially semi-isobaric retrograde P–T path. Received: 2 April 1996 / Accepted: 15 November 1996     

Origin and evolution of 2 Ga old carbonatite complex (Ihouhaouene,Ahaggar, Algeria): Nd and Sr isotopic evidence

Associated syenitic rocks and carbonatites from Ihouhaouene, Algeria, have been investigated for their Sr and Nd isotope and trace element geochemistry. A zircon U-Pb emplacement age (1994 +22 -17 Ma) has been obtained from the carbonatites. The REE characteristics, among which the significant and approximately constant negative Eu anomaly and the evolution of (La/Ce)N and (Yb/Lu)N ratios which increase sympathetically with total REE abundances, are considered to be of purely magmatic origin. They are used to constrain the genetic links between syenites and carbonatites. Sr and Nd isotopes suggest a similar source for carbonatites and syenites, which is enriched compared to a chondritic reference: Nd(T)=-6.4 to -8.6 and ISr(T)=0.7097. These features are interpreted as evidence of contamination of a mantle-derived magmatic precursor by continental crust, occurring in a magmatic chamber.     

Buried and Submerged Greek Archaeological Coastal Structures and Artifacts as Gauges to Measure Late Holocene Seafloor Subsidence off Calabria,Italy

This synthesis integrates recently acquired archaeological and geological data with earlier documented observations to shed light on the subsidence of ancient Greek coastal facilities in southern Italy. These are now positioned between former shorelines and inner shelf sectors at five Calabrian margin localities. Submergence of coastal to inner shelf facilities has resulted in part from sea‐level rise by about 2 m associated with glacio‐hydro‐isostatic factors since archaic to classic Greek time. This phenomenon alone, however, does not explain the wide variation of measured subsidence rates from site‐to‐site. The marked lowering of coastal site substrates by seismo‐tectonic activity (including extensional fault motion), stratal readjustments at depth, and compaction of underlying sediment sequences is significant. Four of the subsided facilities are positioned near emerged Calabrian areas where prevailing Holocene average annual land uplift rates range to ˜1.0 mm/yr; at the fifth, near Hipponion, terrains have risen by nearly 2 mm/yr. In marked contrast, submerged and/or buried structures record the following late Holocene long‐term average rates of coastal margin subsidence: Sybaris‐Thuri on the Taranto Gulf margin (˜0.5–1.0 mm/yr); Hipponion‐Vibo Valentia along the Tyrrhenian coast (˜0.8 to ˜3.2 mm/yr); and Locri‐Epizefiri, Kaulonia, and Capo Colonna on Calabria's Ionian margin (˜1.6, ˜1.6–2.4, and ˜4.0 mm/yr, respectively). © 2012 Wiley Periodicals, Inc.     

Subduction-related mafic to felsic magmatism in the Malayer–Boroujerd plutonic complex,western Iran

The Malayer–Boroujerd plutonic complex (MBPC) in western Iran, consists of a portion of a magmatic arc built by the northeast verging subduction of the Neo-Tethys plate beneath the Central Iranian Microcontinent (CIMC). Middle Jurassic-aged felsic magmatic activity in MBPC is manifested by I-type and S-type granites. The mafic rocks include gabbroic intrusions and dykes and intermediate rocks are dioritic dykes and minor intrusions, as well as mafic microgranular enclaves (MMEs). MBPC Jurassic-aged rocks exhibit arc-like geochemical signatures, as they are LILE- and LREE-enriched and HFSE- and HREE-depleted and display negative Nb–Ta anomalies. The gabbro dykes and intrusions originated from metasomatically enriched garnet-spinel lherzolite [Degree of melting (f_mel) ~ 15%] and exhibit negative Nd and positive to slightly negative ε_Hf(T) (+ 3.0 to ? 1.6). The data reveal that evolution of Middle Jurassic magmatism occurred in two stages: (1) deep mantle-crust interplay zone and (2) the shallow level upper crustal magma chamber. The geochemical and isotopic data, as well as trace element modeling, indicate the parent magma for the MBPC S-type granites are products of upper crustal greywacke (f_mel: 0.2), while I-type granites formed by partial melting of amphibolitic lower crust (f_mel: 0.25) and mixing with upper crustal greywacke melt in a shallow level magma chamber [Degree of mixing (f_mix): 0.3]. Mixing between andesitic melt leaving behind a refractory dense cumulates during partial crystallization of mantle-derived magma and lower crustal partial melt most likely produced MMEs (f_mix: 0.2). However, enriched and moderately variable ε_Nd(T) (? 3.21 to ? 4.33) and high (⁸⁷Sr/⁸⁶Sr)_i (0.7085–0.7092) in dioritic intrusions indicate that these magmas are likely experienced assimilation of upper crustal materials. The interpretations of magmatic activity in the MBPC is consistent with the role considered for mantle-derived magma as heat and mass supplier for initiation and evolution of magmatism in continental arc setting, elsewhere.     

Geometry and kinematics of brittle deformation in the Central Cameroon Shear Zone (Kékem area): Implication for gold exploration within the Central Africa Fold Belt in Cameroon

The Central Africa Fold Belt (CAFB) is a collision belt endowed with gold deposits in Eastern Cameroon area mined for about 50 years. However, favorable areas for gold exploration are poorly known. This paper presents (1) the kinematics of the brittle deformation in the Kékem area in the SW portion of the Central Cameroon Shear Zone and (2) constraints gold mineralization events with respect to the collisional evolution of the CAFB. The authors interpret that the conjugate ENE to E and NNW to NW trending lineament corresponds to the synthetic (R) and the antithetic (R’) shears, which accompanied the dextral slip along the NE to ENE striking shear. The latter coincides with the last 570–552 Ma D₃ dextral simple shear-dominated transpression, which is parallel to the Bétaré Oya shear zone hosting gold deposits. Gold mineralizations, which mainly occurred during the last dextral shearing, are disseminated within quartz veins associated to Riedel’s previous structures reactivated due to late collisional activities of the CAFB as brittle deformation. Gold mineralizations occurred mainly during the 570–552 Ma D₃ event. The reactivation, which might be due to dextral simple shear during mylonitzation, plausibly remobilized the early gold deposits hosted in syn-compressional rocks and/or possibly focused deep-sourced fluid mixed with those released by dehydration. Therefore, the Central Cameroon Shear Zone where Kékem is located, and which shows similar petrographical and structural features to those controling Batouri gold district, is a target area for gold exploration in Cameroon.     

Tectonic evolution of the Lachlan Fold Belt, southeastern Australia: constraints from coupled numerical models of crustal deformation and surface erosion driven by subduction of the underlying mantle

We have used a coupled thermo-mechanical finite-element (FE) model of crustal deformation driven by mantle/oceanic subduction to demonstrate that the tectonic evolution of the Lachlan Fold Belt (LFB) during the Mid-Palaeozoic (Late Ordovician to Early Carboniferous) can be linked to continuous subduction along a single subduction zone. This contrasts with most models proposed to date which assume that separate subduction zones were active beneath the western, central and eastern sections of the Lachlan Orogen. We demonstrate how the existing data on the structural, volcanic and erosional evolution of the Lachlan Fold Belt can be accounted for by our model. We focus particularly on the timing of fault movement in the various sectors of the orogen. We demonstrate that the presence of the weak basal decollement on which most of the Lachlan Fold Belt is constructed effectively decouples crustal structures from those in the underlying mantle. The patterns of faulting in the upper crust appears therefore to be controlled by lateral strength contrasts inherited from previous orogenic events rather than the location of one or several subduction zones. The model also predicts that the uplift and deep exhumation of the Wagga-Omeo Metamorphic Belt (WOMB) is associated with the advection of this terrane above the subduction point and is the only tectonic event that gives us direct constraints on the location of the subduction zone. We also discuss the implications of our model for the nature of the basement underlying the present-day orogen.     

System parameters of soil foundations for time domain dynamic analysis

Soil-structure interaction analysis is usually carried out in the frequency domain, because the compliance functions of the half-space are known only in the frequency domain. Since non-linear analysis cannot be carried out in the frequency domain, a system with frequency independent parameters is used to represent the half-space soil medium so that a nonlinear analysis in the time domain becomes possible. The objective of this paper is to propose a system with lumped parameters, which are independent of frequency, to represent the half-space soil medium. The proposed frequency independent system consists of a number of real discrete structure elements; thus the existing dynamic analysis programs may be adoptable with little modification. In this paper, the parameters are found by minimizing the sum of the squares of deviations between the steady-state responses of the theoretical half-space model and those of the lumped parameter system over a specified frequency range. Once the parameters have been found, the lumped parameter system can be used in practical applications for time domain dynamic analysis of either linear or non-linear structures. In comparison with the dynamic response of the theoretical half-space model, the lumped parameter system yields satisfactory results.