A revised evolutionary history of armadillos (Dasypus) in North America based on ancient mitochondrial DNA

The large, beautiful armadillo, Dasypus bellus, first appeared in North America about 2.5 million years ago, and was extinct across its southeastern US range by 11 thousand years ago (ka). Within the last 150 years, the much smaller nine‐banded armadillo, D. novemcinctus, has expanded rapidly out of Mexico and colonized much of the former range of the beautiful armadillo. The high degree of morphological similarity between these two species has led to speculation that they might be a single, highly adaptable species with phenotypical responses and geographical range fluctuations resulting from environmental changes. If this is correct, then the biology and tolerance limits for D. novemcinctus could be directly applied to the Pleistocene species, D. bellus. To investigate this, we isolated ancient mitochondrial DNA from late Pleistocene‐age specimens of Dasypus from Missouri and Florida. We identified two genetically distinct mitochondrial lineages, which most likely correspond to D. bellus (Missouri) and D. novemcinctus (Florida). Surprisingly, both lineages were isolated from large specimens that were identified previously as D. bellus. Our results suggest that D. novemcinctus, which is currently classified as an invasive species, was already present in central Florida around 10 ka, significantly earlier than previously believed.     

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.     

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.     

Linkages between spatio‐temporal patterns of environmental factors and distribution of plant assemblages across a boreal peatland complex

Here we examine the arrangement of plant species across an oligotrophic bog/poor fen peatland complex in the North American boreal plain and the relationships of these species to their physical and chemical environment. A semi‐uniform spatial sampling approach was utilized to describe the species assemblages, pore‐water chemistry and physical condition of 100 plots throughout a single peatland complex. Regardless of sharing the same ground cover of Sphagnum mosses, the remaining species separated into four distinct assemblages, each with unique indicators. These species groups along with associated chemical and physical factors are organized into four ecosites: bog, margin (edge) and two poor fen ecosites. The plant assemblages of this peatland have a complex relationship with numerous gradients, both physical and chemical, including depth to water table, shade, pH, nutrient and base cation. Rather than being homogenous across the landscape, most environmental variables exhibit distinct spatial patterns and do so in relationship to the plant assemblages, forming spatially distinct ecosites across the complex. Base cation concentrations play a smaller role than previously thought in differentiating these ecosites, and in addition to shade and depth to water table, nitrogen in the form of dissolved organic nitrogen was highly related to the placement of these ecosites. Many significant chemical factors appear related to evaporative water loss within the peatland complex, and these chemical factors are used to differentiate the ecosites. However, the mediation of evaporative water loss is due largely to self‐generated responses of the plant assemblages related to shade through plant morphology and peat acrotelm development related to depth to water table. We conclude that plant species and associated environmental gradients act together to form spatially distinct ecosites. The distribution of these ecosites within this large, environmentally complex peatland is largely controlled by differing self‐generated responses along the hydrotopographical gradient of differential water loss.     

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.