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.     

A critical evaluation of multi‐proxy dendroclimatology in northern Finland

Twentieth‐century summer (July–August) temperatures in northern Finland are reconstructed using ring widths, maximum density and stable carbon isotope ratios (δ¹³C) of Scots pine tree rings, and using combinations of these proxies. Verification is based on the coefficient of determination (r²), reduction of error (RE) and coefficient of efficiency (CE) statistics. Of the individual proxies, δ¹³C performs best, followed by maximum density. Combining δ¹³C and maximum density strengthens the climate signal but adding ring widths leads to little improvement. Blue intensity, an inexpensive alternative to X‐ray densitometry, is shown to perform similarly. Multi‐proxy reconstruction of summer temperatures from a single site produces strong correlations with gridded climate data over most of northern Fennoscandia. Since relatively few trees are required (<15) the approach could be applied to long sub‐fossil chronologies where replication may be episodically low. Copyright © 2010 John Wiley & Sons, Ltd.     

Ressi experimental catchment: Ecohydrological research in the Italian pre-Alps

Ressi is a small (2.4 ha) forested catchment located in the Italian pre-Alps. The site became an experimental catchment to investigate the water fluxes in the soil–plant–atmosphere continuum and the impact of vegetation on runoff generation in 2012. The elevation of the catchment ranges from 598 to 721 m a.s.l. and the climate is humid temperate. The bedrock consists of rhyolites and dacites; the soil is a Cambisol. The catchment is covered by a dense forest, dominated by beech, chestnut, maple, and hazel trees. The field set up includes measurements of the rainfall in an open area, streamflow at the outlet, soil moisture at various depths and locations, and depth to water table in six piezometers at a 5- or 10-min interval. Samples of precipitation, stream water, shallow groundwater and soil water are collected monthly for tracer analysis (stable isotopes (²H and ¹⁸O), electrical conductivity and major ions), and during selected rainfall–runoff events to determine the contribution of the various sources to runoff. Since 2017, soil and plant water samples have been collected to determine the sources of tree transpiration. Data collected in the period 2012–2016 are publicly available. Data collection is ongoing, and the data set is expected to be updated on an annual basis to include the most recent measurements.     

Dynamic storage: a potential metric of inter‐basin differences in storage properties

The potential for dynamic storage to serve as a metric of basin behaviour was assessed using data from five drainage basins with headwaters on the thick sand and gravel deposits of the Oak Ridges Moraine in southern Ontario, Canada. Dynamic storage was directly correlated with the ratio of variability of δ²H in streamflow relative to that in precipitation. This ratio has previously been shown to be inversely related to basin mean transit time (MTT), suggesting an inverse relationship between dynamic storage and MTT for the study basins. Dynamic storage was also directly correlated with interannual variability in stream runoff, baseflow and baseflow:runoff ratio, implying that basins with smaller dynamic storage have less interannual variability in their streamflow regimes. These preliminary results suggest that dynamic storage may serve as a readily derived and useful metric of basin behaviour for inter‐basin comparisons. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

Zircon in amphibolites from Naxos,Aegean Sea,Greece: origin,significance and tectonic setting

We report U–Pb zircon ages of c. 700–550 Ma, 262–220 Ma, 47–38 Ma and 15–14 Ma from amphibolites on Naxos Island in the Aegean extensional province of Greece. The zircon has complex internal structures. Based on cathodoluminescence response, zoning and crosscutting relationships a minimum of four zircon growth stages are identified: inherited core, magmatic core, inner metamorphic (?) rim and an outer metamorphic rim. Trace element compositions of the amphibolites suggest igneous differentiation and crustal assimilation. Zircon solubility as a function of saturation temperatures, Zr content and melt composition indicates that the zircon did not originally crystallize in the mafic bodies but was inherited from felsic precursor rocks, and subsequently assimilated into the mafic intrusives during emplacement. Zircon inheritance is corroborated by the complex, xenocrystic nature of the zircon in one sample. Ages of c. 700–550 Ma and 262–220 Ma are assigned to inherited zircon. Available geochemical data suggest that the 15–14 Ma metamorphic rims grew in situ in the amphibolites, corresponding to a high‐grade metamorphic event at this time. However, the geochemical data cannot conclusively establish if the c. 40 Ma zircon rims also grew in situ, or whether they were inherited along with the xenocrystic cores. Two scenarios for emplacement of the mafic intrusives are discussed: (i) Intrusion during late‐Triassic to Jurassic ocean basin development of the Aegean realm, in which case the 40 Ma zircon rims would have grown in situ, and (ii) emplacement in the Miocene as a result mafic underplating during large‐scale extension. In this case, only the 15–14 Ma metamorphic outer rims would have formed in situ in the amphibolitic host rocks.     

A long‐term study of stable isotopes as tracers of processes governing water flow and quality in a lowland river basin: the upper Thames,UK

A long‐term study of O, H and C stable isotopes has been undertaken on river waters across the 7000‐km² upper Thames lowland river basin in the southern UK. During the period, flow conditions ranged from drought to flood. A 10‐year monthly record (2003–2012) of the main River Thames showed a maximum variation of 3‰ (δ¹⁸O) and 20‰ (δ²H), although interannual average values varied little around a mean of –6.5‰ (δ¹⁸O) and –44‰ (δ²H). A δ²H/δ¹⁸O slope of 5.3 suggested a degree of evaporative enrichment, consistent with derivation from local rainfall with a weighted mean of –7.2‰ (δ¹⁸O) and –48‰ (δ²H) for the period. A tendency towards isotopic depletion of the river with increasing flow rate was noted, but at very high flows (>100 m³/s), a reversion to the mean was interpreted as the displacement of bank storage by rising groundwater levels (corroborated by measurements of specific electrical conductivity). A shorter quarterly study (October 2011–April 2013) of isotope variations in 15 tributaries with varying geology revealed different responses to evaporation, with a well‐correlated inverse relationship between Δ¹⁸O and baseflow index for most of the rivers. A comparison with aquifer waters in the basin showed that even at low flow, rivers rarely consist solely of isotopically unmodified groundwater. Long‐term monitoring (2003–2007) of carbon stable isotopes in dissolved inorganic carbon (DIC) in the Thames revealed a complex interplay between respiration, photosynthesis and evasion, but with a mean interannual δ¹³C‐DIC value of –14.8 ± 0.5‰, exchange with atmospheric carbon could be ruled out. Quarterly monitoring of the tributaries (October 2011–April 2013) indicated that in addition to the aforementioned factors, river flow variations and catchment characteristics were likely to affect δ¹³C‐DIC. Comparison with basin groundwaters of different alkalinity and δ¹³C‐DIC values showed that the origin of river baseflow is usually obscured. The findings show that long‐term monitoring of environmental tracers can help to improve the understanding of how lowland river catchments function. Copyright © NERC 2015. Hydrological Processes © 2015 John Wiley & Sons, Ltd.     

Hammond Hill Research Catchment: Supporting hydrologic investigations of rooting zone and vegetation water dynamics under climate change

The Hammond Hill Research Catchment (HH) is a small (120 ha), temperate, second order tributary to Six Mile Creek, Cayuga Lake, and the Great Lakes (42.42°, −76.32°). The HH has been monitored since January 2017 for the purpose of understanding how recent infiltration mixes with antecedent soil water on hillslope forest floors and the spatial and temporal patterns of Root Water Uptake (RWU) by temperate northeastern US tree species (eastern hemlock [Tsuga canadensis], American beech [Fagus grandifolia], and sugar maple [Acer saccharum]). These data are informing us about the hydrologic consequences of anticipated tree species composition change and supporting the development of more refined ecohydrological models. The glaciated catchment is underlain by a shallow confining siltstone layer (1–1.5 m depth) and densely covered with an approximately 60 year old regrowth mixed species forest of hemlock, beech, and other deciduous tree species common to the northeastern US. Current datasets from the HH include precipitation snow water equivalent, discharge, and associated isotopic water compositions, δ²H & δ¹⁸O. Measurements of (top 10 cm) soil water content, as well as bulk soil water and hemlock and beech xylem isotopic compositions are made at several locations across a topographic wetness gradient. The near-term role of the HH is to support an understanding of the environmental and ecological drivers of plant RWU competition. All data from the HH are publicly available.