δCa evolution in a carbonate aquifer and its bearing on the equilibrium isotope fractionation factor for calcite

To improve understanding of Ca isotope transport during water-rock interaction on the continents, we measured dissolved δ⁴⁴Ca values along a 236 km flow path in the Madison aquifer, South Dakota, where fluids have chemically evolved according to dolomite and anhydrite dissolution, calcite precipitation, and Ca-for-Na ion-exchange over a timescale spanning ~ 15 kyr. We used a reactive transport model employing rate data constrained from major ion mass-balances to evaluate the extent to which calcite precipitation and ion-exchange fractionate Ca isotopes. Elevated δ⁴⁴Ca values during the initial and final stages of water transport possibly result from calcite precipitation under supersaturated conditions and Ca-for-Na ion-exchange, respectively. However, for the bulk of the flow path, δ⁴⁴Ca values evolve by mixing between anhydrite and dolomite dissolution, with no fractionation during calcite precipitation under saturated conditions. We attribute the absence of Ca isotope fractionation to the long timescale of water-rock interaction and slow rate of calcite precipitation, which have enabled fluids to chemically and isotopically equilibrate with calcite. We therefore conclude that the equilibrium Ca isotope fractionation factor between calcite and water (Δ_cal–w) is very close to zero. To the extent that the Madison aquifer typifies other groundwater systems where calcite slowly precipitates from solutions at or near chemical equilibrium, this study suggests that groundwater contributions to δ⁴⁴Ca variability on the continents can be modeled according to simple mixing theory without invoking isotope discrimination.     

Channel flow competence and sediment transport in upland streams in southeast Australia

Bedload transport data from planebed and step‐pool reach types are used to determine grain size transport thresholds for selected upland streams in southeast Australia. Morphological differences between the reach types allow the effects of frictional losses from bedforms, microtopography and bed packing to be incorporated into the dimensionless critical shear stress value. Local sediment transport data are also included in a regime model and applied to mountain streams, to investigate whether empirical data improve the delineation of reach types on the basis of dimensionless discharge per unit width (q*) and dimensionless bedload transport (q_b*). Instrumented planebed and step‐pool sites are not competent to transport surface median grains (D_50s) at bankfull discharge (Q_bf). Application of a locally parametrized entrainment equation to the full range of reach types in the study area indicates that the majority of cascades, cascade‐pools, step‐pools and planebeds are also not competent at Q_bf and require a 10 year recurrence interval flood to mobilize their D_50s. Consequently, the hydraulic parameters of the regime diagram, which assume equilibrium conditions at bankfull, are ill suited to these streams and provide a poor basis of channel delineation. Modifying the diagram to better reflect the dominant transported bedload size (equivalent to the D₁₆ of surface sediment) made only slight improvements to reach delineation and had greatest effect on the morphologies with smaller surface grain sizes such as forced pool‐riffles and planebeds. Likewise, the Corey shape factor was incorporated into the regime diagram as an objective method for adjusting a base dimensionless critical shear stress (τ*_c50b) to account for lithologically controlled grain shape on bed packing and entrainment. However, it too provided only minor adjustments to reach type delineation. Copyright © 2007 John Wiley & Sons, Ltd.     

Conceptual modelling to assess the influence of hydro-climatic variability on runoff processes in data scarce semi-arid Andean catchments

Using hydro-meteorological time series of 50 years and in situ measurements, the dominant runoff processes in perennial Andean headwater catchments in Chile were determined using the hydrological model HBV light. First, cluster analysis was used to identify dry, wet and intermediate years. From these, sub-periods were identified with contrasting seasonal climatic influences on streamflow. By calibrating the model across different periods, impacts on model performance, parameter sensitivity and identifiability were investigated, providing insights into differences in hydrological processes. The modelling approach suggested that, independently of a dry or wet period of calibration, the streamflow response is mostly consistent with flux from groundwater storage, while only a small fraction comes from direct routing of snowmelt. The variation of model parameters, such as the groundwater rate coefficient, was found to be consistent with differing recharge in wet and dry years. The resulting snowmelt–groundwater model is a realistic hypothesis of the hydrological operation of such complex, data scarce and semi-arid Andean catchments. This model may also be a useful tool for predictions of seasonal water availability and a basis for further field studies.     

Using stable isotopes to estimate travel times in a data‐sparse Arctic catchment: Challenges and possible solutions

Use of isotopes to quantify the temporal dynamics of the transformation of precipitation into run‐off has revealed fundamental new insights into catchment flow paths and mixing processes that influence biogeochemical transport. However, catchments underlain by permafrost have received little attention in isotope‐based studies, despite their global importance in terms of rapid environmental change. These high‐latitude regions offer limited access for data collection during critical periods (e.g., early phases of snowmelt). Additionally, spatio‐temporal variable freeze–thaw cycles, together with the development of an active layer, have a time variant influence on catchment hydrology. All of these characteristics make the application of traditional transit time estimation approaches challenging. We describe an isotope‐based study undertaken to provide a preliminary assessment of travel times at Siksik Creek in the western Canadian Arctic. We adopted a model–data fusion approach to estimate the volumes and isotopic characteristics of snowpack and meltwater. Using samples collected in the spring/summer, we characterize the isotopic composition of summer rainfall, melt from snow, soil water, and stream water. In addition, soil moisture dynamics and the temporal evolution of the active layer profile were monitored. First approximations of transit times were estimated for soil and streamwater compositions using lumped convolution integral models and temporally variable inputs including snowmelt, ice thaw, and summer rainfall. Comparing transit time estimates using a variety of inputs revealed that transit time was best estimated using all available inflows (i.e., snowmelt, soil ice thaw, and rainfall). Early spring transit times were short, dominated by snowmelt and soil ice thaw and limited catchment storage when soils are predominantly frozen. However, significant and increasing mixing with water in the active layer during the summer resulted in more damped steam water variation and longer mean travel times (~1.5 years). The study has also highlighted key data needs to better constrain travel time estimates in permafrost catchments.     

Glacial curvilineations found along the southern sector of the Laurentide Ice sheet and a hypothesis of formation involving subglacial slope failure in tunnel valleys and subglacial lakes

First discovered in Poland, glacial curvilineations (GCLs) are enigmatic landforms comprising parallel sets of sinuous ridges and troughs of metres amplitude and around 150 m wavelength, found within kilometres‐wide valleys interpreted as being produced by meltwater flowing subglacially. Their morphological and sedimentary characteristics and association with tunnel valleys has been described for some prominent Polish examples. From these observations the existing hypothesis is that they form as a consequence of erosion by longitudinal vortices that develop in subglacial floods. Here we report, for the first time, GCLs found along the southern sector of the Laurentide Ice Sheet in three northern states of the USA. Using mapping and topographic analysis from high resolution digital elevation models we report observations on their morphological properties and landform associations. We find aspects of their context and morphology difficult to explain using the existing hypothesis. We instead suggest that these glacial curvilineations are produced by subglacial bank and slope failures that locally widen tunnel valleys, or that occur near subglacial lake shorelines. Further investigation is required to test this hypothesis and to ascertain the mechanisms of proposed mass movements, which may have occurred by rotational or translational slope failure or by creep deformation. Our preferred mechanism is that such movements occurred where subglacial water was emplaced over previously perma‐frozen ground. Under such circumstances, sediment blocks thawed by the water may then easily glide over a frozen décollement at low slope angles; analogous to subaerial active‐layer glides in permafrost environments. Permafrost spring sapping may have provided lines of weakness for slope failure. If the requirement for permafrost is found to hold, then GCLs may become an important indicator of the palaeo‐distribution of permafrost. Copyright © 2017 John Wiley & Sons, Ltd.     

Improving bank erosion modelling at catchment scale by incorporating temporal and spatial variability

Bank erosion can contribute a significant portion of the sediment budget within temperate catchments, yet few catchment scale models include an explicit representation of bank erosion processes. Furthermore, representation is often simplistic resulting in an inability to capture realistic spatial and temporal variability in simulated bank erosion. In this study, the sediment component of the catchment scale model SHETRAN is developed to incorporate key factors influencing the spatio‐temporal rate of bank erosion, due to the effects of channel sinuosity and channel bank vegetation. The model is applied to the Eden catchment, north‐west England, and validated using data derived from a GIS methodology. The developed model simulates magnitudes of total catchment annual bank erosion (617–4063 t y^‐1) within the range of observed values (211–4426 t yr^‐1). In addition, the model provides both greater inter‐annual and spatial variability of bank eroded sediment generation when compared with the basic model, and indicates a potential 61% increase of bank eroded sediment as a result of temporal flood clustering. The approach developed within this study can be used within a number of distributed hydrologic models and has general applicability to temperate catchments, yet further development of model representation of bank erosion processes is required. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Using stable water isotopes to identify spatio-temporal controls on groundwater recharge in two contrasting East African aquifer systems

Understanding the spatio-temporal variability in groundwater recharge is a prerequisite to sustainable management of aquifers. Spatial analysis of groundwater stable isotopes uncovered predominant controls on groundwater recharge in the Nairobi aquifer system (NAS) and the South Coast aquifer (SC), two exemplar East African aquifers relied upon by 7 million people. We analysed 368 samples for stable isotopes and basic physico-chemical parameters. The NAS groundwater isotopes are controlled by precipitation orographic effects and enriched recharge from impounded lakes/wetlands; the SC isotopes are correlated with water-table depth influencing evapotranspiration. Global Network of Isotopes in Precipitation (GNIP) data revealed groundwater recharge during months of heavy rains in the NAS, whilst the SC experiences spatio-temporally diffuse recharge. Inferred “isoscapes” show: in the NAS, (1) direct, rapid recharge favoured by faults, well-drained soils and ample rainfall in uplands, (2) delayed recharge from impounded lakes and wetlands in mid-lands, and (3) focused, event-based recharge in floodplains; and in the SC, diffuse recharge complicated by significant water-table evapotranspiration processes.     

Field Study of Soil Vapor Extraction for Reducing Off-Site Vapor Intrusion

Soil vapor extraction (SVE) is effective for removing volatile organic compound (VOC) mass from the vadose zone and reducing the potential for vapor intrusion (VI) into overlying and surrounding buildings. However, the relationship between residual mass in the subsurface and VI is complex. Through a series of alternating extraction (SVE on) and rebound (SVE off) periods, this field study explored the relationship and aspects of SVE applicable to VI mitigation in a commercial/light-industrial setting. The primary objective was to determine if SVE could provide VI mitigation over a wide area encompassing multiple buildings, city streets, and subsurface utilities and eliminate the need for individual subslab depressurization systems. We determined that SVE effectively mitigates offsite VI by intercepting or diluting contaminant vapors that would otherwise enter buildings through foundation slabs. Data indicate a measurable (5 Pa) influence of SVE on subslab/indoor pressure differential may occur but is not essential for effective VI mitigation. Indoor air quality improvements were evident in buildings 100 to 200 feet away from SVE including those without a measurable reversal of differential pressure across the slab or substantial reductions in subslab VOC concentration. These cases also demonstrated mitigation effects across a four-lane avenue with subsurface utilities. These findings suggest that SVE affects distant VI entry points with little observable impact on differential pressures and without relying on subslab VOC concentration reductions.