The role of heterogeneous strain in the development and preservation of a polymetamorphic record in high-P granulites, western Canadian Shield

Mafic rocks in the Chipman domain of the Athabasca granulite terrane, western Canadian Shield, provide the first well‐documented record of two distinct high‐P granulite facies events in the same domain in this region. Textural relations and the results of petrological modelling (NCFMASHT system) of mafic granulites are interpreted in terms of a three‐stage tectonometamorphic history. Stage 1 involved development of the assemblage Grt + Cpx + Qtz ± Pl (M₁) from a primary Opx‐bearing igneous precursor at conditions of 1.3 GPa, 850–900 °C. Field and microstructural observations suggest that M₁ developed synchronously with an early S₁ gneissic fabric. Stage 2 is characterized by heterogeneous deformation (D₂) and synkinematic partial retrogression of the peak assemblage to an amphibole‐bearing assemblage (M₂). Stage 3 involved a third phase of deformation and a return to granulite facies conditions marked by the prograde breakdown of amphibole (Amph₂) to produce matrix garnet (Grt_3a) and the coronitic assemblage Cpx_3b + Opx_3b + Ilm_3b + Pl_3b (M_3b) at 1.0 GPa, 800–900 °C. M₁ and M_3b are correlated with 2.55 and 1.9 Ga metamorphic generations of zircon, respectively, which were dated in a separate study. Heterogeneous strain played a crucial role in both the development and preservation of these rare examples of multiple granulite facies events within single samples. Without this fortuitous set of circumstances, the apparent reaction history could have incorrectly led to an interpretation involving a single‐cycle high‐grade event. The detailed P–T–t–D history constructed for these rocks provides the best evidence to date that much of the east Lake Athabasca region experienced long‐term lower crustal residence from 2.55 to 1.9 Ga, and thus the region represents a rare window into the reactivation and ultimate stabilization processes of cratonic lithosphere.     

Bedrock, soils, and hillslope hydrology in the Central Texas Hill Country, USA: implications on environmental management in a carbonate-rock terrain

The Hill Country of Central Texas, USA, is undergoing rapid socioeconomic development, but environmental management of this region is hampered by misconceptions about local bedrock, soils, terrain, and hydrologic processes. The Hill Country is underlain mostly by Glen Rose Limestone (Lower Cretaceous) and exhibits a stepped terrain, which has been incorrectly attributed to alternating hard and soft bedrock strata. Other characteristics mistakenly attributed to this landscape include thin soils with scant water-retention capabilities, and rapid runoff as the dominant hydrologic process. This report presents new findings: unweathered bedrock is well indurated, but interbeds exhibit variable weathering rates. Recessive slopes (“risers”) on this stepped terrain result from rapid deterioration of strata having generally heterogeneous depositional fabrics (bioturbation and irregular clay partings) in contrast to ledge-forming strata having homogeneous fabrics. A stony regolith is thus formed beneath risers, providing porous and permeable ground that retards runoff and promotes the formation of moderately deep to deep (two-tiered) regolith/soil zones. These surficial materials on local steep slopes compose important natural environmental buffers; they support diverse biota and enhanced geochemical cycling of nutrients; they also exhibit significant water retention and enhanced erosion abatement. Proper land management demands recognition of these attributes in the siting, design, and construction of facilities.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

Low-oxygen and chemical kinetic constraints on the geochemical niche of neutrophilic iron(II) oxidizing microorganisms

Neutrophilic iron oxidizing bacteria (FeOB) must actively compete with rapid abiotic processes governing Fe(II) oxidation and as a result have adapted to primarily inhabit low-O₂ environments where they can more successfully compete with abiotic Fe(II) oxidation. The spatial distribution of these microorganisms can be observed through the chemical gradients they affect, as measured using in situ voltammetric analysis for dissolved Fe(II), Fe(III), O₂, and FeS_(aq). Field and laboratory determination of the chemical environments inhabited by the FeOB were coupled with detailed kinetic competition studies for abiotic and biotic oxidation processes using a pure culture of FeOB to quantify the geochemical niche these organisms inhabit. In gradient culture tubes, the maximum oxygen levels, which were associated with growth bands of Sideroxydans lithotrophicus (ES-1, a novel FeOB), were 15-50 μM. Kinetic measurements made on S. lithotrophicus compared biotic/abiotic (killed control) Fe oxidation rates. The biotic rate can be a significant and measurable fraction of the total Fe oxidation rate below O₂ concentrations of approximately 50 μM, but biotic Fe(II) oxidation (via the biotic/abiotic rate comparison) becomes difficult to detect at higher O₂ levels. These results are further supported by observations of conditions supporting FeOB communities in field settings. Variablity in cell densities and cellular activity as well as variations in hydrous ferrous oxide mineral quantities significantly affect the laboratory kinetic rates. The microbial habitat (or geochemical niche) where FeOB are active is thus largely controlled by the competition between abiotic and biotic kinetics, which are dependent on Fe(II) concentration, P_O2, temperature and pH in addition to the surface area of hydrous ferric oxide minerals and the cell density/activity of FeOB. Additional field and lab culture observations suggest a potentially important role for the iron-sulfide aqueous molecular cluster, FeS_(aq), in the overall cycling of iron associated with the environments these microorganisms inhabit.     

Geochemistry and mineralogy of arsenic in (natural) anaerobic groundwaters

Here new data from field bioremediation experiments and geochemical modeling are reported to illustrate the principal geochemical behavior of As in anaerobic groundwaters. In the field bioremediation experiments, groundwater in Holocene alluvial aquifers in Bangladesh was amended with labile water-soluble organic C (molasses) and MgSO₄ to stimulate metabolism of indigenous SO₄-reducing bacteria (SRB). In the USA, the groundwater was contaminated by Zn, Cd and SO₄, and contained <10 μg/L As under oxidized conditions, and a mixture of sucrose and methanol were injected to stimulate SRB metabolism. In Bangladesh, groundwater was under moderately reducing conditions and contained ∼10 mg/L Fe and ∼100 μg/L As. In the USA experiment, groundwater rapidly became anaerobic, and dissolved Fe and As increased dramatically (As > 1000 μg/L) under geochemical conditions consistent with bacterial Fe-reducing conditions. With time, groundwater became more reducing and biogenic SO₄ reduction began, and Cd and Zn were virtually completely removed due to precipitation of sphalerite (ZnS) and other metal sulfide mineral(s). Following precipitation of chalcophile elements Zn and Cd, the concentrations of Fe and As both began to decrease in groundwater, presumably due to formation of As-bearing FeS/FeS₂. By the end of the six-month experiment, dissolved As had returned to below background levels. In the initial Bangladesh experiment, As decreased to virtually zero once biogenic SO₄ reduction commenced but increased to pre-experiment level once SO₄ reduction ended. In the ongoing experiment, both SO₄ and Fe(II) were amended to groundwater to evaluate if FeS/FeS₂ formation causes longer-lived As removal. Because As-bearing pyrite is the common product of SRB metabolism in Holocene alluvial aquifers in both the USA and Southeast Asia, it was endeavored to derive thermodynamic data for arsenian pyrite to better predict geochemical processes in naturally reducing groundwaters. Including the new data for arsenian pyrite into Geochemist’s Workbench, its stability field completely dominates in reducing Eh–pH space and “displaces” other As-sulfides (orpiment, realgar) that have been implied to be important in previous modeling exercises and reported in rare field conditions.     

In situ strength measurements on natural upper-mantle minerals

Using in situ strength measurements at pressures up to 10 GPa and at room temperature, 400, 600, and 700°C, we examined rheological properties of olivine, orthopyroxene, and chromian-spinel contained in a mantle-derived xenolith. Mineral strengths were estimated using widths of X-ray diffraction peaks as a function of pressure, temperature, and time. Differential stresses of all minerals increase with increasing pressure, but they decrease with increasing temperature because of elastic strain on compression and stress relaxation during heating. During compression at room temperature, all minerals deform plastically at differential stress of 4–6 GPa. During subsequent heating, thermally induced yielding is observed in olivine at 600°C. Neither orthopyroxene nor spinel shows complete stress relaxation, but both retain some stress even at 700°C. The strength of the minerals decreases in the order of chromian-spinel ≈ orthopyroxene > olivine for these conditions. This order of strength is consistent with the residual pressure of fluid inclusions in mantle xenoliths.     

Methodology of hydrogeological characterization of deep carbonate aquifers as potential reservoirs of groundwater. Case of study: the Jurassic aquifer of <Emphasis Type="Italic">El Maestrazgo</Emphasis> (Castellón,Spain)

A methodology for the characterization of deep carbonate aquifers has been developed and applied to El Maestrazgo Jurassic aquifer in Castellón, Spain. Characterization of these aquifer formations, located at more than 300 m deep, consisted of a previous phase of compilation, analysis and synthesis of the existing information about the area, followed by a coordinated combination of different speciality studies: geology, stratigraphy, structural analysis, hydrogeology, hydrochemistry, geophysics and remote sensing. Geological studies included geological mapping, definition of stratigraphical units and facies and structural analysis. The aim of the hydrogeology study was to define aquifer formations, recharge area, aquifer points inventory and groundwater flow directions for the establishment of piezometric and water quality observation nets. Special techniques were applied, like thermal infrared aerial images and the evaluation of submarine groundwater discharge by means of natural radium isotopes. Hydrochemical techniques, including majority elements characterization and stable isotopes (¹⁸O, ²H and ³H) determination, allowed classifying hydrochemical facies and establishing a renewal pattern for water within the system. Geophysics was useful in determining the aquifer geometry, the features of the basement and the petrophysical characteristics of the geological formations. Preliminary results show an important tectonic complexity and the possibilities for groundwater uses in the area of study.     

Spatial Variability of Rock Depth in Bangalore Using Geostatistical, Neural Network and Support Vector Machine Models

Geospatial technology is increasing in demand for many applications in geosciences. Spatial variability of the bed/hard rock is vital for many applications in geotechnical and earthquake engineering problems such as design of deep foundations, site amplification, ground response studies, liquefaction, microzonation etc. In this paper, reduced level of rock at Bangalore, India is arrived from the 652 boreholes data in the area covering 220 km². In the context of prediction of reduced level of rock in the subsurface of Bangalore and to study the spatial variability of the rock depth, Geostatistical model based on Ordinary Kriging technique, Artificial Neural Network (ANN) and Support Vector Machine (SVM) models have been developed. In Ordinary Kriging, the knowledge of the semi-variogram of the reduced level of rock from 652 points in Bangalore is used to predict the reduced level of rock at any point in the subsurface of the Bangalore, where field measurements are not available. A new type of cross-validation analysis developed proves the robustness of the Ordinary Kriging model. ANN model based on multi layer perceptrons (MLPs) that are trained with Levenberg–Marquardt backpropagation algorithm has been adopted to train the model with 90% of the data available. The SVM is a novel type of learning machine based on statistical learning theory, uses regression technique by introducing loss function has been used to predict the reduced level of rock from a large set of data. In this study, a comparative study of three numerical models to predict reduced level of rock has been presented and discussed.     

SOME EXPERIMENTS CONCERNING THE PRIMARY SEISMIC PULSE*

Recordings were made with three types of detector of the primary compressional (P) and shear (S) wave pulses generated by explosions in boreholes. Charge weights varied from 0.08 kg to 9.5 kg and detector distances varied from about 3 m to about 80 m. Scaling by the simple factor W^1/3 where W is the charge weight, enabled observations from different sized charges to be fitted to a single expression. Experiments were carried out in the Bunter sandstone and the London clay and both fluid and solid tamping were used. This variation in tamping had no significant effect on the P-waves but it may have affected the generation of SV-waves. In both media the P-wave energy carried at 30 m from the shot by frequencies less than 100 Hz decreased rapidly with depth and was usually 1–2 % of the available chemical energy for a shot depth of 15 m. The S-wave energy was much less than this, but was highly directional. The P-wave pulse had the appearance of a damped sinusoid in very good agreement with the predictions of the ‘equivalent radiator’ hypothesis. However, the surface of this radiator should be identified not with the blown cavity but with the surface at which the tensile stresses associated with the stress wave become less than the tensile strength of the rock. The predominant frequency for a 1 kg charge at a depth of 15 m was 24 Hz in the clay and 52 Hz in the sandstone. In these and similar media, therefore, an effort should be made to keep individual charges less than 1 kg in reflection shooting and less than 10 kg in refraction shooting. The value of Q was about 50 in clay and about 25 in the sandstone. These estimates are rather uncertain because of the small distances over which the pulses were observed. The Z-transforms of the sampled pulses indicated that they were all of minimum phase, or very near to it.