A novel approach for estimating karst groundwater recharge in mountainous regions and its application in Switzerland

A pragmatic and simple approach for estimating the groundwater recharge of karst aquifers in mountainous regions by extrapolation of the hydrological regimes of gauged and well‐documented systems is presented. Specific discharge rates are derived using annual precipitation and spring measurements by taking into account catchment size and elevation, which are assumed to be the dominant factors. Reference sites with high data reliability are used for calibration and regional extrapolation. This is performed with normalized values employing spatial precipitation deviations and correlation with the elevation of the catchment areas. A tiered step procedure provides minimum and maximum normalized gradients for the relationship between recharge quantity and elevation for karst regions. The normalized recharge can therefore be obtained and extrapolated for any location using the spatial precipitation variability to provide an estimate of annual groundwater recharge. The approach was applied to Switzerland (approximately 7500 km² of karst terrain situated between 200 and over 4000 m a.s.l.) using annual precipitation data from meteorological stations for the years 2000 to 2011. Results show that the average recharge rates of different Swiss karst domains range from 20 to 46 L/km²s, which corresponds to an infiltration ratio between 0.6 and 0.9 of total precipitation. Despite uncertainties inherent in the approach, these results provide a benchmark for renewable karst groundwater resources in Switzerland of about 8.4 km³/year. The approach can be applied to any other mountainous karst region, that is, where a clear relationship between elevation, precipitation and recharge can be assumed. Copyright © 2015 John Wiley & Sons, Ltd.     

A gradual transition to plate tectonics on Earth between 3.2 to 2.7 billion years ago

Zircon crystals precipitated from granitoid magmas contain a robust record of the age and chemistry of continental magmatism spanning some 4.375 Ga of Earth history, a record that charts initiation of plate tectonics. However, constraining when exactly plate tectonics began to dominate crustal growth processes is challenging as the geochemical signatures of individual rocks may reflect local subduction processes rather than global plate tectonics. Here we apply counting statistics to a global database of coupled U–Pb and Hf isotope analyses on magmatic zircon grains from continental igneous and sedimentary rocks to quantify changes in the compositions of their source rocks. The analysis reveals a globally significant change in the sources of granitoid magmas between 3.2 and 2.7 Ga. These secular changes in zircon chemistry are driven by a coupling of the deep (depleted mantle) and shallow (crustal) Earth reservoirs, consistent with a geodynamic regime dominated by Wilson cycle style plate tectonics.     

Glacier outburst floods reconstructed from lake sediments and their implications for Holocene variations of the plateau glacier Folgefonna in western Norway

A sudden release of large volumes of water during a glacier outburst flood (GLOF) is a major hazard worldwide. Here, we identify the sedimentary signature of glacial and non‐glacial processes, including GLOFs, based on lacustrine sediments from the distal glacier‐fed Lake Buarvatnet in western Norway. Historically documented GLOFs in 2002 CE and during the 1980s CE are identified in the ²¹⁰Pb‐ and ¹⁴C‐dated sediments. These events have the same sedimentary signature as 12 earlier events throughout the Holocene interpreted to represent previous GLOFs in the catchment. The GLOFs are interpreted to have occurred during periods when the glacier extent was similar to the modern positions, and the events are thus used to pinpoint past positions of the glacier terminus and, hence, the equilibrium line altitudes (ELAs). The results indicate that the glacier Svartenutbreen, located at the eastern part of Folgefonna, had a similar size in 2002 CE as c. 8200–8300 cal. a BP, corresponding to the 8.2 ka event in the North Atlantic region. The regrowth of Sørfonna after the Holocene Thermal Optimum occurred at c. 6900 cal. a BP and Svartenutbreen was at modern size and extent in the periods c. 6400, c. 5450, c. 4850, c. 3850, c. 3550 and c. 1650 cal. a BP. Since 1650 cal. a BP, we infer that the glacier was larger than the 2002 CE glacier extent until 1910 CE when a GLOF occurred. Svartenutbreen has been retreating since 1910 CE, which led to the ice damming of the two historical GLOFs in the 1980s and 2002 CE separated by a glacier advance in the 1990s CE. The findings are discussed and compared to other regional glacier reconstructions in Norway, and emphasize the value of identifying and utilizing GLOFs as an indicator of past ELA variability.     

A plasticity constitutive model for unsaturated,anisotropic, nonexpansive soils

The paper describes and evaluates an incremental plasticity constitutive model for unsaturated, anisotropic, nonexpansive soils (CMUA). It is based on the modified Cam-Clay (MCC) model for saturated soils and enhances it by introducing anisotropy (via rotation of the MCC yield surface) and an unsaturated compressibility framework describing a double dependence of compressibility on suction and on the degree of saturation of macroporosity. As the anisotropic and unsaturated features can be activated independently, the model is downwards compatible with the MCC model. The CMUA model can simulate effectively: the dependence of compressibility on the level of developed anisotropy, uniqueness of critical state independent of the initial anisotropy, an evolving compressibility during constant suction compression, and a maximum of collapse. The model uses Bishop's average skeleton stress as its first constitutive variable, favouring its numerical implementation in commercial numerical analysis codes (eg, finite element codes) and a unified treatment of saturated and unsaturated material states.     

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.     

Evidence of ameliorated Middle Weichselian climate and sub‐arctic environment in the western Baltic region: coring lake sediments at Klintholm,Møn,Denmark

Coring through glaciotectonically stacked Quaternary sediments situated below sea level on the island of Møn, Denmark, recovered a succession of interstadial sediments of Middle Weichselian age. Plant and animal remains including insects found in laminated sand and mud indicate deposition in a lake surrounded by dwarf shrubs, herbs, mosses and rare trees. The insect fauna indicates a mean July temperature of 8–12 °C, suggesting an arctic to sub‐arctic environment, while winter temperatures around ?8 to ?22 °C suggest periglacial conditions with permafrost. Luminescence dating of sediment samples gave ages from 48–29 ka, and radiocarbon dating indicates deposition of plant fragments between 45 and 36 ka BP. The fossil assemblage from Møn shows close resemblance to those from other sites with similar ages found in the vicinity of the western Baltic Basin.     

Understanding the relative impacts of natural processes and human activities on the hydrology of the Central Rift Valley lakes,East Africa

Significant changes have been observed in the hydrology of Central Rift Valley (CRV) lakes in Ethiopia, East Africa as a result of both natural processes and human activities during the past three decades. This study applied an integrated approach (remote sensing, hydrologic modelling, and statistical analysis) to understand the relative effects of natural processes and human activities over a sparsely gauged CRV basin. Lake storage estimates were calculated from a hydrologic model constructed without inputs from human impacts such as water abstraction and compared with satellite‐based (observed) lake storage measurements to characterize the magnitude of human‐induced impacts. A non‐parametric Mann–Kendall test was used to detect the presence of climatic trends (e.g. a decreasing or increasing trends in precipitation), while the Standard Precipitation Index (SPI) analysis was used to assess the long‐term, inter‐annual climate variability within the basin. Results indicate human activities (e.g. abstraction) significantly contributed to the changes in the hydrology of the lakes, while no statistically significant climatic trend was seen in the basin, however inter‐annual natural climate variability, extreme dryness, and prolonged drought has negatively affected the lakes. The relative contributions of natural and human‐induced impacts on the lakes were quantified and evaluated by comparing hydrographs of the CRV lakes. Lake Abiyata has lost ~6.5 m in total lake height between 1985 and 2006, 70% (~4.5 m) of the loss has been attributed to human‐induced causes, whereas the remaining 30% is related to natural climate variability. The relative impact analysis utilized in this study could potentially be used to better plan and create effective water‐management practices in the basin and demonstrates the utility of this integrated methodology for similar studies assessing the relative natural and human‐induced impacts on lakes in data sparse areas. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Diatom-inferred depth models in 8 Canadian boreal lakes: inferred changes in the benthic:planktonic depth boundary and implications for assessment of past droughts

Assessment of past drought in boreal regions, a region predicted to be strongly affected by climate warming, can provide insights into future availability of water. However, limited instrumental data and paleoclimatic data are available for this assessment. To address this lack of data in the boreal region of northwest Ontario, a regional study of lakes in the Winnipeg River Drainage Basin was initiated. Diatom-inferred (D-I) depth models were developed based on surface samples collected along a depth gradient within 8 small boreal lakes. Weighted-averaging and modern analog approaches provided robust within-lake depth models for each of the study lakes, with bootstrapped r² values ranging from 0.90 to 0.98, and root-mean-squared-errors of prediction (RMSEP) between 1.1 and 2.5 m. Large differences in the estimated depth optima for three representative, but common diatom species across our 8 study lakes suggested that within-lake calibration datasets are more appropriate for inferring past drought based in depth models than a regional multi-lake calibration dataset, and that light and related variables are controlling factors governing the maximum depth of benthic taxa. A down-core application of the D-I depth models on a near-shore core from Meekin Lake, retrieved near the present-day ecotone between the benthic and planktonic diatom assemblages indicated highly similar trends in inferred depth (r = 0.96). The models have significant correlations with other metrics of changes in depth including diatom species richness (r = 0.74–0.78) and evenness (r = 0.76–0.8), thereby allowing a check on the strength and direction of the depth inferences down-core. Near-shore cores located near the benthic:planktonic transition is a sensitive region that can provide estimates of past droughts in lakes where such inferences have been difficult to estimate.