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.     

Soil cracking effects on hydrological and erosive processes: a study case in Mediterranean cultivated vertisols

Shrink–swell soils, such as those in a Mediterranean climate regime, can cause changes in terms of hydrological and erosive responses because of the changing soil water storage conditions. Only a limited number of long‐term studies have focused on the impacts on both hydrological and erosive responses and their interactions in an agricultural environment. In this context, this study aims to document the dynamics of cracks, runoff and soil erosion within a small Mediterranean cultivated catchment and to quantify the influence of crack processes on the water and sediment supplied to a reservoir located at the catchment outlet using water and sediment measurements at a cultivated field outlet as baseline. Detailed monitoring of the presence of topsoil cracks was conducted within the Kamech catchment (ORE OMERE, Tunisia), and runoff and suspended sediment loads were continuously measured over a long period of time (2005–2012) at the outlets of a field (1.32 ha) and a catchment (263 ha). Analysis of the data showed that topsoil cracks were open approximately half of the year and that the rainfall regime and water table level conditions locally control the seasonal cracking dynamics. Topsoil cracks appeared to seriously affect the generation of runoff and sediment concentrations and, consequently, sediment yields, with similar dynamics observed at the field and catchment outlets. A similar time lag in the seasonality between water and sediment delivery was observed at these two scales: although the runoff rates were globally low during the presence of topsoil cracks, most sediment transport occurred during this period associated with very high sediment concentrations. This study underlines the importance of a good prediction of runoff during the presence of cracks for reservoir siltation considerations. In this context, the prediction of cracking effects on runoff and soil erosion is a key factor for the development of effective soil and water management strategies and downstream reservoir preservation. Copyright © 2016 John Wiley & Sons, Ltd.     

Extreme value modelling of daily areal rainfall over Mediterranean catchments in a changing climate

Heavy rainfall events during the fall season are causing extended damages in Mediterranean catchments. A peaks‐over‐threshold model is developed for the extreme daily areal rainfall occurrence and magnitude in fall over six catchments in Southern France. The main driver of the heavy rainfall events observed in this region is the humidity flux (FHUM) from the Mediterranean Sea. Reanalysis data are used to compute the daily FHUM during the period 1958–2008, to be included as a covariate in the model parameters. Results indicate that the introduction of FHUM as a covariate can improve the modelling of extreme areal precipitation. The seasonal average of FHUM can improve the modelling of the seasonal occurrences of heavy rainfall events, whereas daily FHUM values can improve the modelling of the events magnitudes. In addition, an ensemble of simulations produced by five different general circulation models are considered to compute FHUM in future climate with the emission scenario A1B and hence to evaluate the effect of climate change on the heavy rainfall distribution in the selected catchments. This ensemble of climate models allows the evaluation of the uncertainties in climate projections. By comparison to the reference period 1960–1990, all models project an amplification of the mean seasonal FHUM from the Mediterranean Sea for the projection period 2070–2099, on average by +22%. This increase in FHUM leads to an increase in the number of heavy rainfall events, from an average of 2.55 events during the fall season in present climate to 3.57 events projected for the period 2070–2099. However, the projected changes have limited effects on the magnitude of extreme events, with only a 5% increase in the median of the 100‐year quantiles. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrogen emissions from Erebus volcano, Antarctica

The continuous measurement of molecular hydrogen (H₂) emissions from passively degassing volcanoes has recently been made possible using a new generation of low-cost electrochemical sensors. We have used such sensors to measure H₂, along with SO₂, H₂O and CO₂, in the gas and aerosol plume emitted from the phonolite lava lake at Erebus volcano, Antarctica. The measurements were made at the crater rim between December 2010 and January 2011. Combined with measurements of the long-term SO₂ emission rate for Erebus, they indicate a characteristic H₂ flux of 0.03?kg s^–1 (2.8?Mg? day^–1). The observed H₂ content in the plume is consistent with previous estimates of redox conditions in the lava lake inferred from mineral compositions and the observed CO₂/CO ratio in the gas plume (～0.9 log units below the quartz–fayalite–magnetite buffer). These measurements suggest that H₂ does not combust at the surface of the lake, and that H₂ is kinetically inert in the gas/aerosol plume, retaining the signature of the high-temperature chemical equilibrium reached in the lava lake. We also observe a cyclical variation in the H₂/SO₂ ratio with a period of ～10?min. These cycles correspond to oscillatory patterns of surface motion of the lava lake that have been interpreted as signs of a pulsatory magma supply at the top of the magmatic conduit.     

Mechanisms of orthopyroxene dissolution in silica-undersaturated melts at 1 atmosphere and implications for the origin of silica-rich glass in mantle xenoliths

Experiments dissolving orthopyroxene (En₉₃) in a variety of Si-undersaturated alkaline melts at 1 atmosphere and variable f _O2 demonstrate that orthopyroxene dissolves to form olivine, Si-rich melt and clinopyroxene. These phases form a texturally and chemically distinct boundary layer around the partly dissolved orthopyroxene crystals. The occurrence of clinopyroxene in the boundary layer is due to inward diffusion of Ca from the solvent melt to the boundary layer causing clinopyroxene saturation. Compositional profiles through the solvent and the boundary layer for a number of experiments demonstrate rapid diffusion of cations across the boundary layer – solvent interface. SiO₂ diffuses outward from the boundary layer whereas CaO and Al₂O₃ diffuse toward the Si-enriched boundary layer melt. The rate of Al diffusion is slower under reducing conditions compared to the rates in experiments performed in air. Concentrations of FeO and MgO in the boundary layer and solvent are approximately equal indicating rapid diffusion and attainment of equilibrium despite ongoing crystallisation of clinopyroxene within the boundary layer. The behaviour of Na₂O and K₂O is strongly affected by f _O2. Under reducing conditions Na₂O and K₂O concentrations are approximately equal in the boundary layer and solvent indicating normal diffusion down the concentration gradient and attainment of equilibrium. Under oxidising conditions, K₂O and to a lesser extent Na₂O, have compositional profiles indicative of uphill diffusion likely due to their preference for more polymerised Si- and Al-rich melts. Under reduced conditions Al-enrichment in the boundary layer melt is not as extreme and uphill diffusion did not occur. The composition of the solvent melt after the experiments indicates that it was contaminated by the boundary layer by convective mixing due to the onset of hydrodynamic instabilities brought on by density and viscosity contrasts between the two melts. Despite using a wide variety of solvent melt compositions we find that the boundary layer melts converge toward a common composition at high SiO₂ contents. The composition of glass generated by orthopyroxene dissolution at 1 atmosphere is similar in many respects to Si-rich glass found in many orthopyroxene-rich mantle xenoliths that have been attributed to high pressure in situ processes including mantle metasomatism. The results of this study suggest that at least some Si-rich melts are likely to have formed by dissolution of xenolith orthopyroxene at low pressure possibly by their Si-undersaturated host magmas. Received: 30 August 1996 / Accepted: 15 April 1998     

Evidence of an early alteration process driven by magmatic fluid in Mururoa volcano

Evidence for a deuteric alteration process induced by a magmatic fluid has been found in the feeder zone of the Mururoa volcano (French Polynesia). Within the dikes, where basaltic glass does not show any evidence of pervasive alteration, vesicles are filled with dioctahedral smectites and calcite, while olivine phenocrysts are replaced by dioctahedral smectites, ankerite and calcite.The ¹³C signature of carbonates, the carbon and H₂O content of the whole rocks and their impoverishment in deuterium are compatible with the presence of magmatic CO₂ during the crystallization of intruding lavas and exclude contamination by seawater. Mass balance calculations on selected thin sections photographs of partly filled up vesicles and replaced olivine crystals, constrain, assuming a closed system interaction, the chemical composition of the initial fluid and the respective amounts of the initial solid phases involved in the alteration process. Thermodynamic modelings using the EQ3/6 software package correctly predict the mineralogic, chemical and isotopic exchanges accompanying alteration, thus validating the closed system assumption. The model which allows prediction of the influence of CO₂ on the alteration products, shows that, above a 0.25 CO₂ mole fraction in the initial fluid, the alteration is entirely controlled by the chemical composition of the initial solid phases. The presence of CO₂ implies the precipitation of dioctahedral smectites and carbonates instead of the magnesian smectites commonly observed in CO₂-free systems.The Mururoa feeder zone shows alteration features typical of a closed system interaction between the basaltic rock and a magmatic fluid in which seawater did not take part.     

The snowball Earth aftermath: Exploring the limits of continental weathering processes

Carbonates capping Neoproterozoic glacial deposits contain peculiar sedimentological features and geochemical anomalies ascribed to extraordinary environmental conditions in the snowball Earth aftermath. It is commonly assumed that post-snowball climate dominated by CO₂ partial pressures several hundred times greater than modern levels, would be characterized by extreme temperatures, a vigorous hydrological cycle, and associated high continental weathering rates. However, the climate in the aftermath of a global glaciation has never been rigorously modelled. Here, we use a hierarchy of numerical models, from an atmospheric general circulation model to a mechanistic model describing continental weathering processes, to explore characteristics of the Earth system during the supergreenhouse climate following a snowball glaciation. These models suggest that the hydrological cycle intensifies only moderately in response to the elevated greenhouse. Indeed, constraints imposed by the surface energy budget sharply limit global mean evaporation once the temperature has warmed sufficiently that the evaporation approaches the total absorbed solar radiation. Even at 400 times the present day pressure of atmospheric CO₂, continental runoff is only 1.2 times the modern runoff. Under these conditions and accounting for the grinding of the continental surface by the ice sheet during the snowball event, the simulated maximum discharge of dissolved elements from continental weathering into the ocean is approximately 10 times greater than the modern flux. Consequently, it takes millions of years for the silicate weathering cycle to reduce post-snowball CO₂ levels to background Neoproterozoic levels. Regarding the origin of the cap dolostones, we show that continental weathering alone does not supply enough cations during the snowball melting phase to account for their observed volume.     

Comparison of LiDAR waveform processing methods for very shallow water bathymetry using Raman,near‐infrared and green signals

Airborne light detection and ranging (LiDAR) bathymetry appears to be a useful technology for bed topography mapping of non‐navigable areas, offering high data density and a high acquisition rate. However, few studies have focused on continental waters, in particular, on very shallow waters (<2 m) where it is difficult to extract the surface and bottom positions that are typically mixed in the green LiDAR signal. This paper proposes two new processing methods for depth extraction based on the use of different LiDAR signals [green, near‐infrared (NIR), Raman] of the SHOALS‐1000T sensor. They have been tested on a very shallow coastal area (Golfe du Morbihan, France) as an analogy to very shallow rivers. The first method is based on a combination of mathematical and heuristic methods using the green and the NIR LiDAR signals to cross validate the information delivered by each signal. The second method extracts water depths from the Raman signal using statistical methods such as principal components analysis (PCA) and classification and regression tree (CART) analysis. The obtained results are then compared to the reference depths, and the performances of the different methods, as well as their advantages/disadvantages are evaluated. The green/NIR method supplies 42% more points compared to the operator process, with an equivalent mean error (?4·2 cm verusu ?4·5 cm) and a smaller standard deviation (25·3 cm verusu 33·5 cm). The Raman processing method provides very scattered results (standard deviation of 40·3 cm) with the lowest mean error (?3·1 cm) and 40% more points. The minimum detectable depth is also improved by the two presented methods, being around 1 m for the green/NIR approach and 0·5 m for the statistical approach, compared to 1·5 m for the data processed by the operator. Despite its ability to measure other parameters like water temperature, the Raman method needed a large amount of reference data to provide reliable depth measurements, as opposed to the green/NIR method. Copyright © 2010 John Wiley & Sons, Ltd.