What is the source of baseflow in agriculturally fragmented catchments? Complex groundwater/surface‐water interactions in three tributary catchments of the Wabash River,Indiana, USA

Some conceptual models suggest that baseflow in agriculturally fragmented watersheds may contain little, if any, groundwater. This has critical implications for stream quality and ecosystem functioning. Here, we (a) identify the sources and flowpaths contributing to baseflow using ²²²Rn and ⁸⁷Sr/⁸⁶Sr and (b) quantify mean apparent ages of groundwater and baseflow using multiple isotopic tracers (CFC, SF₆, ³⁶Cl, and ³H) in 4 small (0.08 to 0.64 km²) tributary catchments to the Wabash River in Indiana, USA. ²²²Rn activities and ⁸⁷Sr/⁸⁶Sr ratios indicate that baseflow in 3 catchments is sourced primarily from groundwater; baseflow in the fourth is dominated by a source similar to agricultural run‐off. CFC‐12 data indicate that springs in 1 catchment are discharging significant proportions of water that recharged between 1974 (42 ± 2 years) and 1961 (55 ± 2 years). Those same springs have ³⁶Cl/Cl ratios between 1,381.08 ± 29.37 (×10^?15) and 1,530.64 ± 27.65 (×10^?15) indicating that a substantial proportion of the discharge likely recharged between 1975 (41 years) and 1950 (66 years). Groundwater samples collected from streambed mini‐piezometers in a separate catchment have CFC‐12 concentrations indicating that a large proportion of the recharge occurred between 1948 (68 ± 2 years) and 1950 (66 ± 2 years). Repeat sampling conducted in September 2015 after above‐average summer rainfall did not show significant decreases in mean apparent age. The relatively old ages observed in 3 of the catchments can be explained by geological complexities that are likely present in all 4 catchments, but overwhelmed by flow from the shallow phreatic aquifer in the fourth catchment.     

Experimental determination of REE partition coefficients between zircon,garnet and melt: a key to understanding high‐T crustal processes

The partitioning of rare earth elements (REE) between zircon, garnet and silicate melt was determined using synthetic compositions designed to represent partial melts formed in the lower crust during anatexis. The experiments, performed using internally heated gas pressure vessels at 7 kbar and 900–1000 °C, represent equilibrium partitioning of the middle to heavy REE between zircon and garnet during high‐grade metamorphism in the mid to lower crust. The D_REE (zircon/garnet) values show a clear partitioning signature close to unity from Gd to Lu. Because the light REE have low concentrations in both minerals, values are calculated from strain modelling of the middle to heavy REE experimental data; these results show that zircon is favoured over garnet by up to two orders of magnitude. The resulting general concave‐up shape to the partitioning pattern across the REE reflects the preferential incorporation of middle REE into garnet, with D_Gd (zircon/garnet) ranging from 0.7 to 1.1, D_Ho (zircon/garnet) from 0.4 to 0.7 and D_Lu (zircon/garnet) from 0.6 to 1.3. There is no significant temperature dependence in the zircon–garnet REE partitioning at 7 kbar and 900–1000 °C, suggesting that these values can be applied to the interpretation of zircon–garnet equilibrium and timing relationships in the ultrahigh‐T metamorphism of low‐Ca pelitic and aluminous granulites.     

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.     

High‐grade metamorphism and partial melting in Archean composite grey gneiss complexes

Much of the exposed Archean crust is composed of composite gneiss which includes a large proportion of intermediate to tonalitic material. These gneiss terranes were typically metamorphosed to amphibolite to granulite facies conditions, with evidence for substantial partial melting at higher grade. Recently published activity–composition (a?x) models for partial melting of metabasic to intermediate compositions allows calculation of the stable metamorphic minerals, melt production and melt composition in such rocks for the first time. Calculated P?T pseudosections are presented for six bulk rock compositions taken from the literature, comprising two metabasic compositions, two intermediate/dioritic compositions and two tonalitic compositions. This range of bulk compositions captures much of the diversity of rock types found in Archean banded gneiss terranes, enabling us to present an overview of metamorphism and partial melting in such terranes. If such rocks are fluid saturated at the solidus, they first begin to melt in the upper amphibolite facies. However, at such conditions, very little (< 5%) melt is produced and this melt is granitic in composition for all rocks. The production of greater proportions of melt requires temperatures ～800–850 °C and is associated with the first appearance of orthopyroxene at pressures below 8–9 kbar or with the appearance and growth of garnet at higher pressures. The temperature at which orthopyroxene appears varies little with composition providing a robust estimate of the amphibolite–granulite facies boundary. Across this boundary, melt production is coincident with the breakdown of hornblende and/or biotite. Melts produced at granulite facies range from tonalite–trondhjemite–granodiorite for the metabasic protoliths, granodiorite to granite for the intermediate protoliths and granite for the tonalitic protoliths. Under fluid‐absent conditions the melt fertility of the different protoliths is largely controlled by the relative proportions of hornblende and quartz at high grade, with the intermediate compositions being the most fertile. The least fertile rocks are the most leucocratic tonalites due to their relatively small proportions of hydrous mafic phases such as hornblende or biotite. In the metabasic rocks, melt production becomes limited by the complete consumption of quartz to higher temperatures. The use of phase equilibrium forward‐modelling provides a thermodynamic framework for understanding melt production, melt loss and intracrustal differentiation during the Archean.     

Zircon U–Pb ages and Hf isotope compositions of migmatites from the North Qinling terrane and their geological implications

Migmatites are predominant in the North Qinling (NQ) orogen, but their formation ages are poorly constrained. This paper presents a combined study of cathodoluminescence imaging, U–Pb age, trace element and Hf isotopes of zircon in migmatites from the NQ unit. In the migmatites, most zircon grains occur as new, homogeneous crystals, while some are present as overgrowth rims around inherited cores. Morphological and trace element features suggest that the zircon crystals are metamorphic and formed during partial melting. The inherited cores have oscillatory zoning and yield U–Pb ages of c. 900 Ma, representing their protolith ages. The early Neoproterozoic protoliths probably formed in an active continental margin, being a response to the assembly of the supercontinent Rodinia. The migmatite zircon yields Hf model ages of 1911 ± 20 to 990 ± 22 Ma, indicating that the protoliths were derived from reworking of Palaeoproterozoic to Neoproterozoic crustal materials. The anatexis zircon yields formation ages ranging from 455 ± 5 to 420 ± 4 Ma, with a peak at c. 435 Ma. Combined with previous results, we suggest that the migmatization of the NQ terrane occurred at c. 455–400 Ma. The migmatization was c. 50 Ma later than the c. 490 Ma ultra‐high‐P (UHP) metamorphism, indicating that they occurred in two independent tectonic events. By contrast, the migmatization was coeval with the granulite facies metamorphism and the granitic magmatism in the NQ unit, which collectively argue for their formation due to the northward subduction of the Shangdan Ocean. UHP rocks were distributed mainly along the northern margin and occasionally in the inner part of the NQ unit, indicating that they were exhumed along the northern edge and detached from the basement by the subsequent migmatization process.     

New mineral activity–composition relations for thermodynamic calculations in metapelitic systems

New activity–composition (a–x) relations for minerals commonly occurring in metapelites are presented for use with the internally consistent thermodynamic dataset of Holland & Powell ( 2011 , Journal of Metamorphic Geology, 29 , 333–383). The a–x relations include a broader consideration of Fe₂O₃ in minerals, changes to the formalism of several phases and order–disorder in all ferromagnesian minerals where Fe–Mg mixing occurs on multiple sites. The a–x relations for chlorite, biotite, garnet, chloritoid, staurolite, cordierite, orthopyroxene, muscovite, paragonite and margarite have been substantially reparameterized using the approach outlined in the companion paper in this issue. For the first time, the entire set of a–x relations for the common ferromagnesian minerals in metapelitic rocks is parameterized simultaneously, with attention paid to ensuring that they can be used together to calculate phase diagrams of geologically appropriate topology. The a–x relations developed are for use in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCKFMASHTO) system for both subsolidus and suprasolidus conditions. Petrogenetic grids in KFMASH and KFMASHTO are similar in topology to those produced with earlier end‐member datasets and a–x relations, but with some notable differences. In particular, in subsolidus equilibria, the FeO/(FeO + MgO) of garnet is now greater than in coexisting staurolite, bringing a number of key staurolite‐bearing equilibria into better agreement with inferences from field and petrographic observations. Furthermore, the addition of Fe³⁺ and Ti to a number of silicate phases allows more plausible equilibria to be calculated in relevant systems. Pseudosections calculated with the new a–x relations are also topologically similar to equivalent diagrams using earlier a–x relations, although with many low variance fields shifting in P–T space to somewhat lower pressure conditions.     

Numerical simulation on mining-induced water inrushes related to geologic structures using a damage-based hydromechanical model

A large number of statistics indicate that water inrush has a direct relationship with geological structures such as fault and karst collapse columns. Understanding the mechanism of water inrushes controlled by geologic structures is of vital importance for adopting effective measures to prevent their occurrence. The work begins with formulization of a damage-based hydromechanical model based on elastic damage theory. Next, the model is numerically implemented with finite element method by employing a finite element package called COMSOL Multiphysics, and is also validated against some existing experimental observations. Finally, the model is used to simulate the mining-induced groundwater inrushes when the effect of faults and karst collapse columns is considered in the numerical simulation, and some suggestive conclusions for preventing water inrushes and optimizing underground mining operations are drawn.     

Using monazite and zircon petrochronology to constrain the P–T–t evolution of the middle crust in the Bhutan Himalaya

The growth and dissolution behaviour of accessory phases (and especially those of geochronological interest) in metamorphosed pelites depends on, among others, the bulk composition, the prograde metamorphic evolution and the cooling path. Monazite and zircon are arguably the most commonly used geochronometers for dating felsic metamorphic rocks, yet crystal growth mechanisms as a function of rock composition, pressure and temperature are still incompletely understood. Ages of different growth zones in zircon and monazite in a garnet‐bearing anatectic metapelite from the Greater Himalayan Sequence in NW Bhutan were investigated via a combination of thermodynamic modelling, microtextural data and interpretation of trace‐element chemical ‘fingerprint’ indicators in order to link them to the metamorphic stage at which they crystallized. Differences in the trace‐element composition (HREE, Y, Eu_N/Eu*_N) of different phases were used to track the growth/dissolution of major (e.g. plagioclase, garnet) and accessory phases (e.g. monazite, zircon, xenotime, allanite). Taken together, these data constrain multiple pressure–temperature–time (P–T–t) points from low temperature (<550 °C) to upper amphibolite facies (partial melting, >700 °C) conditions. The results suggest that the metapelite experienced a cryptic early metamorphic stage at c. 38 Ma at <550 °C, ≥0.85 GPa during which plagioclase was probably absent. This was followed by a prolonged high‐T, medium‐pressure (~600 °C, 0.55 GPa) evolution at 35–29 Ma during which the garnet grew, and subsequent partial melting at >690 °C and >18 Ma. Our data confirm that both geochronometers can crystallize independently at different times along the same P–T path and that neither monazite nor zircon necessarily provides timing constraints on ‘peak’ metamorphism. Therefore, collecting monazite and zircon ages as well as major and trace‐element data from major and accessory phases in the same sample is essential for reconstructing the most coherent metamorphic P–T–t evolution and thus for robustly constraining the rates and timescales of metamorphic cycles.     

Constraints on the timing and conditions of high‐grade metamorphism,charnockite formation and fluid–rock interaction in the Trivandrum Block,southern India

Incipient charnockites have been widely used as evidence for the infiltration of CO₂‐rich fluids driving dehydration of the lower crust. Rocks exposed at Kakkod quarry in the Trivandrum Block of southern India allow for a thorough investigation of the metamorphic evolution by preserving not only orthopyroxene‐bearing charnockite patches in a host garnet–biotite felsic gneiss, but also layers of garnet–sillimanite metapelite gneiss. Thermodynamic phase equilibria modelling of all three bulk compositions indicates consistent peak‐metamorphic conditions of 830–925 °C and 6–9 kbar with retrograde evolution involving suprasolidus decompression at high temperature. These models suggest that orthopyroxene was most likely stabilized close to the metamorphic peak as a result of small compositional heterogeneities in the host garnet–biotite gneiss. There is insufficient evidence to determine whether the heterogeneities were inherited from the protolith or introduced during syn‐metamorphic fluid flow. U–Pb geochronology of monazite and zircon from all three rock types constrains the peak of metamorphism and orthopyroxene growth to have occurred between the onset of high‐grade metamorphism at c. 590 Ma and the onset of melt crystallization at c. 540 Ma. The majority of metamorphic zircon growth occurred during protracted melt crystallization between c. 540 and 510 Ma. Melt crystallization was followed by the influx of aqueous, alkali‐rich fluids likely derived from melts crystallizing at depth. This late fluid flow led to retrogression of orthopyroxene, the observed outcrop pattern and to the textural and isotopic modification of monazite grains at c. 525–490 Ma.