Recent progress on combining geomorphological and geochronological data with ice sheet modelling,demonstrated using the last British–Irish Ice Sheet

Palaeo-ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo-ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model-based approaches. To achieve this goal, model-data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice-flow direction and flow-switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice-free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo-ice sheet reconstructions with greater explicative and ultimately predictive power.     

A probabilistic approach for computing water retention of particulate systems from statistics of grain size and tessellated pore network

The paper offers an analytical determination of the hydraulic properties of an unsaturated soil with reference to its retention curve, which describes the relationship between the volumetric water content and capillarity through matric suction. The analysis combines a particulate approach focused on the physics at the pore scale, including microstructural aspects, with a probabilistic approach where the void space and grain size are considered as random variables. In the end, the soil water characteristic curve of an unsaturated granular medium along a drying path can be derived analytically based on the sole information of particle size distribution. The analysis hinges on the tessellation of a wet granular system into an assemblage of tetrahedral unit cells revealing a pore network upon which capillary physics are computed with respect to pore throat invasion by a non-wetting fluid with evolving pendular capillary bridges. The crux of the paper is to pass from particle size probability distribution to a matching void space distribution to eventually reveal key information such as void cell and solid volume statistics. Making reasonable statistically based assumptions to render calculations tractable, the water retention curve can be readily constructed. Model predictions compare quite favourably with experimental data available for actual soils, especially in the high saturation range. Having a sound scientific basis, the model can be made amenable to address a variety of soils with a wider range of particle sizes.     

A hydropedological approach to simulate streamflow and soil water contents with SWAT+

Reflecting internal catchment hydrological processes in hydrological models is important for accurate predictions of the impact of climate and land-use change on water resources. Characterizing these processes is however difficult and expensive due to their dynamic nature and spatio-temporal variability. Hydropedology is a relatively new discipline focusing on the synergistic integration of hydrology, soil physics and pedology. Hydropedological interpretations of soils and soil distribution can be used to characterize key hydrological processes, especially in areas with no or limited hydrometric measurements. Here we applied a hydropedological approach to reflect flowpaths through detailed routing in SWAT+ for a 157 ha catchment (Weatherley) in South Africa. We compared the hydropedological approach and a standard (no routing) approach against measured streamflow (two weirs) and soil water contents (13 locations). The catchment was treated as ‘ungauged’ and the model was not calibrated against hydrometric measurements in order to determine the direct contribution of hydropedology on modelling efficiency. Streamflow was predicted well without calibration (NSE > 0.8; R² > 0.82) for both approaches at both weirs. The standard approach yielded slightly better streamflow predictions. The hydropedological approach resulted in considerable improvements in the simulation of soil water contents (R² increased from 0.40 to 0.49 and PBIAS decreased from 40% to 20%). The routing capacity of SWAT+ as employed in the hydropedological approach reduced the underestimation of wetland water regimes drastically and resulted in a more accurate representation of the dominant hydrological processes in this catchment. We concluded that hydropedology can be a valuable source of ‘soft data’ to reflect internal catchment structure and processes and, potentially, for realistic calibrations in other studies, especially those conducted in areas with limited hydrometric measurements.     

Hydrologic and water isotope characterization of a regulated Canadian Shield river basin

Stable water isotope surveys have increasingly been integrated into river basins studies, but fewer have used them to evaluate impact of hydropower regulation. This study applies hydrologic and water isotope survey approaches to a Canadian Shield river basin with both regulated and natural flows. Historical streamflow records were used to evaluate the influence of three hydroelectric reservoirs and unregulated portions of the basin on downstream flows and changes in water level management implemented after an extreme flood year (1979). In 2013, water isotope surveys of surface and source waters (e.g., rainfall, groundwater, snowmelt) were conducted to examine spatial and temporal variation in contributions to river flow. Seasonal changes in relative groundwater contribution were assessed using a water‐isotope mass balance approach. Within the basin, two regulated reservoirs exhibited inverted hydrographs with augmented winter flows, whereas a third exhibited a hydrograph dominated by spring snowmelt. In 2013, spatial variation in rain‐on‐snow and air temperatures resulted in a critical lag in snowmelt initiation in the southern and northern portions of the basin resulting in a dispersed, double peak spring hydrograph, contrasting with 1979 when a combination of rain‐on‐snow and coincident snowmelt led to the highest flood on record. Although eastern basin reservoirs become seasonally enriched in δ¹⁸O and δ²H values, unregulated western basin flows remain less variable due to groundwater driven baseflow with increasing influence downstream. Combined analysis of historical streamflow (e.g., flood of 1979, drought of 2010) and the 2013 water isotope surveys illustrate extreme meteorological conditions that current management activities are unable to prevent. In this study, the influence of evaporative fractionation on large surface water reservoirs provides important evidence of streamflow partitioning, illustrating the value of stable water isotope tracers for study of larger catchments.     

Traces of volcanic ash from the Mediterranean,Iceland and North America in a Holocene record from south Wales,UK

A tephra record is presented for a sediment core from Llyn Llech Owain, south Wales, spanning the early- to mid-Holocene. Seven cryptotephra deposits are discovered with three thought to correlate with known eruptions and the remaining four considered to represent previously undocumented events. One deposit is suggested to correlate with the ~6.9 cal ka bp Lairg A tephra from Iceland, whereas more distant sources are proposed as the origin for two of the tephra deposits. A peak of colourless shards in early-Holocene sediments is thought to tentatively correlate with the ~9.6 cal ka bp Fondi di Baia tephra (Campi Flegrei) and a second cryptotephra is tentatively correlated with the ~3.6 cal ka bp Aniakchak (CFE) II tephra (Alaska). The Fondi di Baia tephra has never been recorded beyond proximal sites and its discovery in south Wales significantly extends the geographical distribution of ash from this eruption. The remaining four cryptotephra deposits are yet to be correlated with known eruptions, demonstrating that our current understanding of widespread tephra deposits is incomplete. This new tephra record highlights the potential for sites at more southerly and westerly locations in northwest Europe to act as repositories for ash from several volcanic regions.     

Micromechanical approach to swelling behavior of capillary-porous media with coupled physics

A detailed multiscale analysis is presented of the swelling phenomenon in unsaturated clay-rich materials in the linear regime through homogenization. Herein, the structural complexity of the material is formulated as a three-scale, triple porosity medium within which microstructural information is transmitted across the various scales, leading ultimately to an enriched stress-deformation relation at the macroscopic scale. As a side note, such derived relationship leads to a tensorial stress partitioning that is reminiscent of a Terzaghi-like effective stress measure. Otherwise, a major result that stands out from previous works is the explicit expression of swelling stress and capillary stress in terms of micromechanical interactions at the very fine scale down to the clay platelet level, along with capillary stress emerging due to interactions between fluid phases at the different scales, including surface tension, pore size, and morphology. More importantly, the swelling stress is correlated with the disjoining forces due to electrochemical effects of charged ions on clay minerals and van der Waals forces at the nanoscale. The resulting analytical expressions also elucidate the role of the various physics in the deformational behavior of clayey material. Finally, the capability of the proposed formulation in capturing salient behaviors of unsaturated expansive clays is illustrated through some numerical examples.     

Dissolved nitrous oxide (N2O) dynamics in agricultural field drains and headwater streams in an intensive arable catchment

Indirect nitrous oxide (N₂O) emissions produced by nitrogen (N) leaching into surface water and groundwater bodies are poorly understood in comparison to direct N₂O emissions from soils. In this study, dissolved N₂O concentrations were measured weekly in both lowland headwater streams and subsurface agricultural field drain discharges over a 2‐year period (2013–2015) in an intensive arable catchment, Norfolk, UK. All field drain and stream water samples were found to have dissolved N₂O concentrations higher than the water–air equilibrium concentration, illustrating that all sites were acting as a net source of N₂O emissions to the atmosphere. Soil texture was found to significantly influence field drain N₂O dynamics, with mean concentrations from drains in clay loam soils (5.3 μg N L^?1) being greater than drains in sandy loam soils (4.0 μg N L^?1). Soil texture also impacted upon the relationships between field drain N₂O concentrations and other water quality parameters (pH, flow rate, and nitrate (NO₃) and nitrite (NO₂) concentrations), highlighting possible differences in N₂O production mechanisms in different soil types. Catchment antecedent moisture conditions influenced the storm event mobilisation of N₂O in both field drains and streams, with the greatest concentration increases recorded during precipitation events preceded by prolonged wet conditions. N₂O concentrations also varied seasonally, with the lowest mean concentrations typically occurring during the summer months (JJA). Nitrogen fertiliser application rates and different soil inversion regimes were found to have no effect on dissolved N₂O concentrations, whereas higher N₂O concentrations recorded in field drains under a winter cover crop compared to fallow fields revealed cover crops are an ineffective greenhouse gas emission mitigation strategy. Overall, this study highlights the complex interactions governing the dynamics of dissolved N₂O concentrations in field drains and headwater streams in a lowland intensive agricultural catchment.