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.     

Integration of aquifer geology,groundwater flow and arsenic distribution in deltaic aquifers – A unifying concept

Groundwater arsenic (As) presents a public health risk of great magnitude in densely populated Asian delta regions, most acutely in the Bengal Basin (West Bengal, India and Bangladesh). Research has focused on the sources, mobilisation, and heterogeneity of groundwater As, but a consistent explanation of As distribution from local to basin scale remains elusive. We show for the Bengal Aquifer System that the numerous, discontinuous silt‐clay layers together with surface topography impose a hierarchical pattern of groundwater flow, which constrains As penetration into the aquifer and controls its redistribution towards discharge zones, where it is re‐sequestered to solid phases. This is particularly so for the discrete periods of As release to groundwater in the shallow subsurface associated with sea level high‐stand conditions of Quaternary inter‐glacial periods. We propose a hypothesis concerning groundwater flow ( S ilt‐clay layers I mpose H ierarchical groundwater flow patterns constraining A rsenic progression [SIHA]), which links consensus views on the As source and history of sedimentation in the basin to the variety of spatial and depth distributions of groundwater As reported in the literature. SIHA reconciles apparent inconsistencies between independent, in some cases contrasting, field observations. We infer that lithological and topographic controls on groundwater flow, inherent to SIHA, apply more generally to deltaic aquifers elsewhere. The analysis suggests that groundwater As may persist in the aquifers of Asian deltas over thousands of years, but in certain regions, particularly at deeper levels, As will not exceed low background concentrations unless groundwater flow systems are short‐circuited by excessive pumping.     

A GLUE‐based uncertainty assessment framework for tritium‐inferred transit time estimations under baseflow conditions

The last decade has seen major technical and scientific improvements in the study of water transfer time through catchments. Nevertheless, it has been argued that most of these developments used conservative tracers that may disregard the oldest component of water transfer, which often has transit times greater than 5 years. Indeed, although the analytical reproducibility of tracers limits the detection of the older flow components associated with the most dampened seasonal fluctuations, this is very rarely taken into account in modelling applications. Tritium is the only environmental tracer at hand to investigate transfer times in the 5‐ to 50‐year range in surface waters, as dissolved gases are not suitable due to the degassing process. Water dating with tritium has often been difficult because of the complex history of its atmospheric concentration, but its current stabilization together with recent analytical improvements open promising perspectives. In this context, the innovative contribution of this study lies in the development of a generalized likelihood uncertainty estimation‐based approach for analysing the uncertainties associated with the modelling of transit time due to both parameter identification and tracer analytical precision issues. A coupled resampling procedure allows assessment of the statistical significance of the transfer time differences found in diverse waters. This approach was developed for tritium and the exponential‐piston model but can be implemented for virtually any tracer and model. Stream baseflow, spring and shallow aquifer waters from the Vallcebre research catchments, analysed for tritium in different years with different analytical precisions, were investigated by using this approach and taking into account other sources of uncertainty. The results showed three groups of waters of different mean transit times, with all the stream baseflow and spring waters older than the 5‐year threshold needing tritium. Low sensitivity of the results to the model structure was also demonstrated. Dual solutions were found for the waters sampled in 2013, but these results may be disambiguated when additional analyses will be made in a few years. Copyright © 2016 John Wiley & Sons, Ltd.     

The June 2013 Alberta catastrophic flooding event – part 2: fine‐scale precipitation and associated features

Data obtained from a variety of sources including the Canadian Lightning Detection Network, weather radars, weather stations and operational numerical weather model analyses were used to address the evolution of precipitation during the June 2013 southern Alberta flood. The event was linked to a mid‐level closed low pressure system to the west of the region and a surface low pressure region initially to its south. This configuration brought warm, moist unstable air into the region that led to dramatic, organized convection with an abundance of lightning and some hail. Such conditions occurred in the southern parts of the region whereas the northern parts were devoid of lightning. Initially, precipitation rates were high (extreme 15‐min rainfall rates up to 102 mm h^?1 were measured) but decreased to lower values as the precipitation shifted to long‐lived stratiform conditions. Both the convective and stratiform precipitation components were affected by the topography. Similar flooding events, such as June 2002, have occurred over this region although the 2002 event was colder and precipitation was not associated with substantial convection over southwest Alberta. Copyright © 2016 Her Majesty the Queen in Right of Canada. Hydrological Processes. © John Wiley & Sons, Ltd.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

Characterization of dissolved organic matter across the Subtropical Convergence off the South Island,New Zealand

Dissolved organic matter (DOM) was investigated along a gradient across the Subtropical Convergence (STC) off the South Island, New Zealand. Ultrahigh resolution electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS), excitation emission matrix fluorescence (EEM) spectroscopy, and molecular lignin analysis techniques were used to study this DOM. The analysis revealed a group of compounds found only in the coastal DOM samples, which were also characterized by an elevated terrestrial DOM fluorescence pattern and elevated lignin content if compared to the offshore samples. This group exhibited a high degree of carbon unsaturation, as evident from high double bond equivalence minus oxygen values (DBE-O > 9) and maximum fluorescence intensity. Sulfur-containing molecular formulae for summer DOM samples were much more abundant across the entire transect of the STC compared to winter and exhibited distinctly different hydrogen:carbon and oxygen:carbon elemental ratios, suggesting a possible correlation between biological activity and sulfur compounds. The molecular formulae common to all STC samples were compared with those calculated for DOM extracted from freshwater collected from a stream discharging into Doubtful Sound (West Coast, New Zealand). ESI-FT-ICR-MS measurements undertaken in negative electrospray ionization (ESI?) mode indicated that 30% of the molecular formulae were present in both types of DOM, while in positive ESI mode (ESI+) over 90% of the formulae were present in all samples. Hence, a significant portion of the molecular formulae assigned to the solid-phase extractable DOM pool (SPE-DOM) appear to be identical in oceanic and freshwater samples.     

Reach‐scale contributions of road‐surface sediment to the Honna River,Haida Gwaii,BC

Gravel road surfaces can be a major source of fine sediment to streams, yet their contribution to channel reach sediment balances remains poorly documented. To quantify the input of road surface material and to compare this input with natural sediment sources at the reach scale, suspended sediment dynamics was examined and a 16‐month sediment balance was developed for a ~35 channel‐width (approx. 425 m) reach of the Honna River, a medium‐size, road‐affected stream located in coastal British Columbia. Of the 105 ± 33 t of suspended material passing through the reach, 18 ± 6% was attributed to the road surface. The high availability of sediment on the road surface appears to limit hysteresis in road run‐off. During rainstorms that increase streamflow, road surface material composed 0.5–15% of sediment inputs during relatively dry conditions from April to the end of September and 5–70% through wetter conditions from October to the end of March, but our data do not show evidence of major sediment accumulation on the riverbed in the reach. A comparison of modelled sediment production on the road surface with observed yields from drainage channels suggests that (1) during low intensity rainfall, ditches and drainage channels may trap sediment from road run‐off, which is subsequently released during events of greater intensity, and/or (2) production models do not effectively describe processes, such as deposition or erosion of sediment in ditches, which control sediment transport and delivery. Our findings further emphasize the risk of unpaved roads in polluting river systems and highlight the continued need for careful road design and location away from sensitive aquatic environments. Copyright © 2016 John Wiley & Sons, Ltd.     

Climate change and future flows of Rocky Mountain rivers: converging forecasts from empirical trend projection and down‐scaled global circulation modelling

To predict future river flows, empirical trend projection (ETP) analyses and extends historic trends, while hydroclimatic modelling (HCM) incorporates regional downscaling from global circulation model (GCM) outputs. We applied both approaches to the extensively allocated Oldman River Basin that drains the North American Rocky Mountains and provides an international focus for water sharing. For ETP, we analysed monthly discharges from 1912 to 2008 with non‐parametric regression, and extrapolated changes to 2055. For modelling, we refined the physical models MTCLIM and SNOPAC to provide water inputs into RIVRQ (river discharge), a model that assesses the streamflow regime as involving dynamic peaks superimposed on stable baseflow. After parameterization with 1960–1989 data, we assessed climate forecasts from six GCMs: CGCM1‐A, HadCM3, NCAR‐CCM3, ECHAM4 and 5 and GCM2. Modelling reasonably reconstructed monthly hydrographs (R² about 0·7), and averaging over three decades closely reconstructed the monthly pattern (R² = 0·94). When applied to the GCM forecasts, the model predicted that summer flows would decline considerably, while winter and early spring flows would increase, producing a slight decline in the annual discharge (?3%, 2005–2055). The ETP predicted similarly decreased summer flows but slight change in winter flows and greater annual flow reduction (?9%). The partial convergence of the seasonal flow projections increases confidence in a composite analysis and we thus predict further declines in summer (about ? 15%) and annual flows (about ? 5%). This composite projection indicates a more modest change than had been anticipated based on earlier GCM analyses or trend projections that considered only three or four decades. For other river basins, we recommend the utilization of ETP based on the longest available streamflow records, and HCM with multiple GCMs. The degree of correspondence from these two independent approaches would provide a basis for assessing the confidence in projections for future river flows and surface water supplies. Copyright © 2010 John Wiley & Sons, Ltd.     

Progradational and backstepping shoreface deposits in the Ladinian to Early Norian Snadd Formation of the Barents Sea

Shoreface sandstone deposits within the Early Carnian part of the Snadd Formation of the Norwegian Barents Sea can be traced for hundreds of kilometres in the depositional strike direction and for tens of kilometres in the depositional‐dip direction. This study uses three‐dimensional seismic attribute mapping and two‐dimensional regional seismic profiles to visualize the seismic facies of these shoreface deposits and to map their internal stratigraphic architecture at a regional scale. The shoreface deposits are generally elongate but show variable width from north‐east to south‐west, which corresponds to a sediment source in the northern part of the basin and a southward decrease in longshore sediment transport. The Snadd Formation presents an example of how large‐scale progradational shoreface deposits develop. The linear nature of its shoreface deposits contrasts with more irregular, cuspate wave‐dominated deltaic shorelines that contain river outlets, and instead implies longshore drift as the main sediment source. In map view, discrete sets of linear features bounded by truncation surfaces scale directly to beach ridge sets in modern counterparts. The shoreface deposits studied here are characteristic in terms of scale and basin‐wide continuity, and offer insight into the contrast between shallow marine deposition under stable Triassic Greenhouse and fluctuating Holocene Icehouse climates. Findings presented herein are also important for hydrocarbon exploration in the Barents Sea, because they describe a hitherto poorly understood reservoir play in the Triassic interval, wherein the most prominent reservoir plays have so far been considered to be found in channelized deposits in net‐progradational delta‐plain strata that form the topsets to shelf‐edge clinoforms. The documented presence of widespread wave‐dominated shoreface deposits also has implications for how the relative importance of different sedimentary processes is considered within the basin during this period.     

Anisotropic mechanical properties of zircon and the effect of radiation damage

This study provides new insights into the relationship between radiation-dose-dependent structural damage due to natural U and Th impurities and the anisotropic mechanical properties (Poisson’s ratio, elastic modulus and hardness) of zircon. Natural zircon samples from Sri Lanka (see Muarakami et al. in Am Mineral 76:1510–1532, 1991) and synthetic samples, covering a dose range of zero up to 6.8 × 10¹⁸ α-decays/g, have been studied by nanoindentation. Measurements along the [100] crystallographic direction and calculations, based on elastic stiffness constants determined by Özkan (J Appl Phys 47:4772–4779, 1976), revealed a general radiation-induced decrease in stiffness (~54 %) and hardness (~48 %) and an increase in the Poisson’s ratio (~54 %) with increasing dose. Additional indentations on selected samples along the [001] allowed one to follow the amorphization process to the point that the mechanical properties are isotropic. This work shows that the radiation-dose-dependent changes of the mechanical properties of zircon can be directly correlated with the amorphous fraction as determined by previous investigations with local and global probes (Ríos et al. in J Phys Condens Matter 12:2401–2412, 2000a; Farnan and Salje in J Appl Phys 89:2084–2090, 2001; Zhang and Salje in J Phys Condens Matter 13:3057–3071, 2001). The excellent agreement, revealed by the different methods, indicates a large influence of structural and even local phenomena on the macroscopic mechanical properties. Therefore, this study indicates the importance of acquiring better knowledge about the mechanical long-term stability of radiation-damaged materials.