Proglacial groundwater storage dynamics under climate change and glacier retreat

Proglacial aquifers are an important water store in glacierised mountain catchments that supplement meltwater-fed river flows and support freshwater ecosystems. Climate change and glacier retreat will perturb water storage in these aquifers, yet the climate-glacier-groundwater response cascade has rarely been studied and remains poorly understood. This study implements an integrated modelling approach that combines distributed glacio-hydrological and groundwater models with climate change projections to evaluate the evolution of groundwater storage dynamics and surface-groundwater exchanges in a temperate, glacierised catchment in Iceland. Focused infiltration along the meltwater-fed Virkisá River channel is found to be an important source of groundwater recharge and is projected to provide 14%–20% of total groundwater recharge by the 2080s. The simulations highlight a mechanism by which glacier retreat could inhibit river recharge in the future due to the loss of diurnal melt cycling in the runoff hydrograph. However, the evolution of proglacial groundwater level dynamics show considerable resilience to changes in river recharge and, instead, are driven by changes in the magnitude and seasonal timing of diffuse recharge from year-round rainfall. The majority of scenarios simulate an overall reduction in groundwater levels with a maximum 30-day average groundwater level reduction of 1 m. The simulations replicate observational studies of baseflow to the river, where up to 15% of the 30-day average river flow comes from groundwater outside of the melt season. This is forecast to reduce to 3%–8% by the 2080s due to increased contributions from rainfall and meltwater runoff. During the melt season, groundwater will continue to contribute 1%–3% of river flow despite significant reductions in meltwater runoff inputs. Therefore it is concluded that, in the proglacial region, groundwater will continue to provide only limited buffering of river flows as the glacier retreats.     

Seismic receiver coupling to the seafloor

The presence of geophysical receivers on the seafloor changes the local wave field due to the receiver seafloor interaction. The resulting PP- and PS-wave distortion of the wave field is often referred to as receiver coupling to the seafloor and can make data processing challenging and sometimes impossible. This paper provides an overview of the mathematical approaches to describe receiver coupling, how to estimate receiver coupling and what the difficulties and possible solutions are. The first section shows how the mathematical approach developed from a simple model considering only the vertical receiver component to include all three receiver components and their interactions with the seafloor. In the second section, I show how receiver coupling can be measured and how it can be improved using mathematical and data-driven approaches.     

Fluvial dynamics and 14C-10Be disequilibrium on the Bolivian Altiplano

Determining sediment transfer times is key to understanding source-to-sink dynamics and the transmission of environmental signals through the fluvial system. Previous work on the Bolivian Altiplano applied the in situ cosmogenic ¹⁴C-¹⁰Be-chronometer to river sands and proposed sediment storage times of ~10–20 kyr in four catchments southeast of Lake Titicaca. However, the fidelity of those results hinges upon isotopic steady-state within sediment supplied from the source area. With the aim of independently quantifying sediment storage times and testing the ¹⁴C-¹⁰Be steady-state assumption, we dated sediment storage units within one of the previously investigated catchments using radiocarbon dating, cosmogenic ¹⁰Be-²⁶Al isochron burial dating, and ¹⁰Be-²⁶Al depth-profile dating. Palaeosurfaces appear to preserve remnants of a former fluvial system, which has undergone drainage reversal, reduction in catchment area, and local isostatic uplift since ~2.8 Ma. From alluvium mantling the palaeosurfaces we gained a deposition age of ~580 ka, while lower down fluvial terraces yielded ≤34 ka, and floodplains ~3–1 ka. Owing to restricted channel connectivity with the terraces and palaeosurfaces, the main source of channel sediment is via reworking of the late Holocene floodplain. Yet modelling a set of feasible scenarios reveals that floodplain storage and burial depth are incompatible with the ¹⁴C-¹⁰Be disequilibrium measured in the channel. Instead we propose that the ¹⁴C-¹⁰Be offset results from: (i) non-uniform erosion whereby deep gullies supply hillslope-derived debris; and/or (ii) holocene landscape transience associated with climate or human impact. The reliability of the ¹⁴C-¹⁰Be chronometer vitally depends upon careful evaluation of sources of isotopic disequilibrium in a wide range of depositional and erosional landforms in the landscape. © 2018 John Wiley & Sons, Ltd.     

Investigating groundwater flow components in an Alpine relict rock glacier (Austria) using a numerical model

Relict rock glaciers are complex hydrogeological systems that might act as relevant groundwater storages; therefore, the discharge behavior of these alpine landforms needs to be better understood. Hydrogeological and geophysical investigations at a relict rock glacier in the Niedere Tauern Range (Austria) reveal a slow and fast flow component that appear to be related to the heterogeneous structure of the aquifer. A numerical groundwater flow model was used to indicate the influence of important internal structures such as layering, preferential flow paths and aquifer-base topography. Discharge dynamics can be reproduced reasonably by both introducing layers of strongly different hydraulic conductivities or by a network of highly conductive channels within a low-conductivity zone. Moreover, the topography of the aquifer base influences the discharge dynamics, which can be observed particularly in simply structured aquifers. Hydraulic conductivity differences of three orders of magnitude are required to account for the observed discharge behavior: a highly conductive layer and/or channel network controlling the fast and flashy spring responses to recharge events, as opposed to less conductive sediment accumulations sustaining the long-term base flow. The results show that the hydraulic behavior of this relict rock glacier and likely that of others can be adequately represented by two aquifer components. However, the attempt to characterize the two components by inverse modeling results in ambiguity of internal structures when solely discharge data are available.     

Modeling planform evolution of a mud‐dominated meandering river: Quinn River,Nevada, USA

A depth‐averaged linearized meander evolution model was calibrated and tested using the field data collected at the Quinn River in the Black Rock Desert, Nevada. Two approaches used to test the model were: (1) simulating meander evolution and comparing the results with the observed 38 year migration pattern; and (2) fitting the model parameters to present bank asymmetry (the ratio of the maximum bank gradients on opposite sides of the channel). The data required as input were collected in the field during a high flow in May 2011 and from aerial photographs and LiDAR data. Both approaches yielded similar results for the best fit parameter values. The bank asymmetry analysis showed that the bank asymmetry and the velocity perturbation have high correlation at close to zero spatial lag while the maximum correlation between the bank asymmetry and maximum bend curvature is offset by about 25 m. The model sufficiently replicated 38 years of channel migration, with a few locations significantly under‐ or over‐predicted. Inadequacies of the flow model and/or variation in bank properties unaccounted for are most likely the causes for these discrepancies. Flow through the Quinn River was also simulated by a more general 3D model. The downstream pattern of near‐bank shear stresses simulated by the 3D model is nearly identical to those resulting from the linearized flow model. Topographic profiles across interior bends are essentially invariant over a wide range of migration rates, suggesting that the traditional formulation that cut bank erosion processes govern migration rates is appropriate for the Quinn River. Copyright © 2014 John Wiley & Sons, Ltd.     

Chemical evolution in the high arsenic groundwater of the Huhhot basin (Inner Mongolia,PR China) and its difference from the western Bengal basin (India)

Elevated As concentrations in groundwater of the Huhhot basin (HB), Inner Mongolia, China, and the western Bengal basin (WBB), India, have been known for decades. However, few studies have been performed to comprehend the processes controlling overall groundwater chemistry in the HB. In this study, the controls on solute chemistry in the HB have been interpreted and compared with the well-studied WBB, which has a very different climate, physiography, lithology, and aquifer characteristics than the HB. In general, there are marked differences in solute chemistry between HB and WBB groundwaters. Stable isotopic signatures indicate meteoric recharge in the HB in a colder climate, distant from the source of moisture, in comparison to the warm, humid WBB. The major-ion composition of the moderately reducing HB groundwater is dominated by a mixed-ion (Ca–Na–HCO₃–Cl) hydrochemical facies with an evolutionary trend along the regional hydraulic gradient. Molar ratios and thermodynamic calculations show that HB groundwater has not been affected by cation exchange, but is dominated by weathering of feldspars (allitization) and equilibrium with gibbsite and anorthite. Mineral weathering and mobilization of As could occur as recharging water flows through fractured, argillaceous, metamorphic or volcanic rocks in the adjoining mountain-front areas, and deposits solutes near the center of the basin. In contrast, WBB groundwater is Ca–HCO₃-dominated, indicative of calcite weathering, with some cation exchange and silicate weathering (monosiallitization).     

Assessing the contribution of off‐fault deformation to slip‐rate estimates within the Taupo Rift,New Zealand,using 3‐D ground‐penetrating radar surveying and trenching

We demonstrate that conventional palaeoseismic trenching and mapping techniques that do not account for the effects of off‐fault deformation can significantly underestimate a fault’s slip rate. Using combined interpretations of 3‐D ground‐penetrating radar (GPR) and palaeoseismic trench data, we show that drag folding and hangingwall and footwall horizontal‐axis rotations have accommodated up to 41% of total extension across a normal fault within the Taupo Rift, New Zealand, over the past 24.6 ± 1.0 cal. ka BP. Our results may explain why geologically determined fault‐slip rates for the central and southern Taupo Rift are anomalously low when compared with geodetic estimates. We suggest that a combination of GPR surveying and palaeoseismic trenching may help resolve differences between geodetically and geologically determined strain rates observed across active extensional regimes worldwide.     

Changes in Jupiter’s zonal velocity between 1979 and 2008

We show that the peak velocity of Jupiter’s visible-cloud-level zonal winds near 24°N (planetographic) increased from 2000 to 2008. This increase was the only change in the zonal velocity from 2000 to 2008 for latitudes between ±70° that was statistically significant and not obviously associated with visible weather. We present the first automated retrieval of fast (∼130 m s⁻¹) zonal velocities at 8°N planetographic latitude, and show that some previous retrievals incorrectly found slower zonal winds because the eastward drift of the dark projections (associated with 5-μm hot spots) “fooled” the retrieval algorithms.We determined the zonal velocity in 2000 from Cassini images from NASA’s Planetary Data System using a global method similar to previous longitude-shifting correlation methods used by others, and a new local method based on the longitudinal average of the two-dimensional velocity field. We obtained global velocities from images acquired in May 2008 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). Longer-term variability of the zonal winds is based on comparisons with published velocities based on 1979 Voyager 2 and 1995-1998 HST images. Fluctuations in the zonal wind speeds on the order of 10 m s⁻¹ on timescales ranging from weeks to months were found in the 1979 Voyager 2 and the 1995-1998 HST velocities. In data separated by 10 h, we find that the east-west velocity uncertainty due to longitudinal fluctuations are nearly 10 m s⁻¹, so velocity fluctuations of 10 m s⁻¹ may occur on timescales that are even smaller than 10 h. Fluctuations across such a wide range of timescales limit the accuracy of zonal wind measurements. The concept of an average zonal velocity may be ill-posed, and defining a “temporal mean” zonal velocity as the average of several zonal velocity fields spanning months or years may not be physically meaningful.At 8°N, we use our global method to find peak zonal velocities of ∼110 m s⁻¹ in 2000 and ∼130 m s⁻¹ in 2008. Zonal velocities from 2000 Cassini data produced by our local and global methods agree everywhere, except in the vicinity of 8°N. There, the local algorithm shows that the east-west velocity has large variations in longitude; vast regions exceed ∼140 m s⁻¹. Our global algorithm, and all of the velocity-extraction algorithms used in previously-published studies, found the east-west drift velocities of the visible dark projections, rather than the true zonal velocity at the visible-cloud level. Therefore, the apparent increase in zonal winds between 2000 and 2008 at 8°N is not a true change in zonal velocity.At 7.3°N, the Galileo probe found zonal velocities of 170 m s⁻¹ at the 3-bar level. If the true zonal velocity at the visible-cloud level at this latitude is ∼140 m s⁻¹ rather than ∼105 m s⁻¹, then the vertical zonal wind shear is much less than the currently accepted value.     

Land‐use impacts and water quality targets in the intensive dairying catchment of the Toenepi Stream,New Zealand

Water quality monitoring in Toenepi Stream, New Zealand, started in 1995 in a study of dairy farming influences on lowland stream quality and has continued since then with brief interruptions. Surveys have provided information about changes in farm and soil management practices as they relate to environmental sustainability. Although average water quality in Toenepi Stream has changed little during 1995–2004, there have been some notable improvements. Water clarity measured by black disc has improved from 0.6m to 1.5m, and median ammonia‐N and nitrate‐N concentrations have declined by 70% and 57%, respectively. The frequency and magnitude of extreme concentrations have declined—most notably for nitrogen (N) forms, which also had decreased mean values. Specific yields for suspended solids (SS) and phosphorus (P) forms in 2002–04 were 47–67% of 1995–97 values, mainly because of lower water yields. Reduced specific yields for N forms in 2002–04 (34–37% of 1995–97 yields) were also attributable to lower mean concentrations in stream water. Faecal bacteria concentrations have not abated and are on average 2–3 times recommended guideline values for contact recreation. Fewer dairy farms and an increased proportion irrigating dairyshed effluent to land, rather than discharging it to the stream via two‐pond systems, were likely causes of improvement in water quality. Water quality targets were developed for Toenepi Stream to achieve contact recreation criteria for the Piako River (downstream) and for intrinsic habitat values for Toenepi Stream. A range of mitigation measures has been formulated to meet these targets, but substantial uptake of sustainable farming practices is needed to improve water quality in Toenepi Stream.     

Influence of sea surface wind wave turbulence upon wind-induced circulation, tide–surge interaction and bed stress

A three-dimensional finite volume unstructured mesh model of the west coast of Britain, with high resolution in the coastal regions, is used to investigate the role of wind wave turbulence and wind and tide forced currents in producing maximum bed stress in the eastern Irish Sea. The spatial distribution of the maximum bed stress, which is important in sediment transport problems, is determined, together with how it is modified by the direction of wind forced currents, tide–surge interaction and a surface source of wind wave turbulence associated with wave breaking. Initial calculations show that to first order the distribution of maximum bed stress is determined by the tide. However, since maximum sediment transport occurs at times of episodic events, such as storm surges, their effects upon maximum bed stresses are examined for the case of strong northerly, southerly and westerly wind forcing. Calculations show that due to tide–surge interaction both the tidal distribution and the surge are modified by non-linear effects. Consequently, the magnitude and spatial distribution of maximum bed stress during major wind events depends upon wind direction. In addition calculations show that a surface source of turbulence due to wind wave breaking in shallow water can influence the maximum bed stress. In turn, this influences the wind forced flow and hence the movement of suspended sediment. Calculations of the spatial variability of maximum bed stress indicate the level of measurements required for model validation.