Identification of glacial meltwater runoff in a karstic environment and its implication for present and future water availability

Glaciers all over the world are expected to continue to retreat due to the global warming throughout the 21st century. Consequently, future seasonal water availability might become scarce once glacier areas have declined below a certain threshold affecting future water management strategies. Particular attention should be paid to glaciers located in a karstic environment, as parts of the meltwater can be drained by underlying karst systems, making it difficult to assess water availability. In this study tracer experiments, karst modeling and glacier melt modeling are combined in order to identify flow paths in a high alpine, glacierized, karstic environment (Glacier de la Plaine Morte, Switzerland) and to investigate current and predict future downstream water availability. Flow paths through the karst underground were determined with natural and fluorescent tracers. Subsequently, geologic information and the findings from tracer experiments were assembled in a karst model. Finally, glacier melt projections driven with a climate scenario were performed to discuss future water availability in the area surrounding the glacier. The results suggest that during late summer glacier meltwater is rapidly drained through well-developed channels at the glacier bottom to the north of the glacier, while during low flow season meltwater enters into the karst and is drained to the south. Climate change projections with the glacier melt model reveal that by the end of the century glacier melt will be significantly reduced in the summer, jeopardizing water availability in glacier-fed karst springs.

Quantitative high-resolution warm season rainfall recorded in varved sediments of Lake Oeschinen, northern Swiss Alps: calibration and validation AD 1901–2008

High-resolution, well-calibrated records of lake sediments are critically important for quantitative climate reconstructions, but they remain a methodological and analytical challenge. While several comprehensive paleotemperature reconstructions have been developed across Europe, only a few quantitative high-resolution studies exist for precipitation. Here we present a calibration and verification study of lithoclastic sediment proxies from proglacial Lake Oeschinen (46°30′N, 7°44′E, 1,580 m a.s.l., north–west Swiss Alps) that are sensitive to rainfall for the period AD 1901–2008. We collected two sediment cores, one in 2007 and another in 2011. The sediments are characterized by two facies: (A) mm-laminated clastic varves and (B) turbidites. The annual character of the laminae couplets was confirmed by radiometric dating (210Pb, 137Cs) and independent flood-layer chronomarkers. Individual varves consist of a dark sand-size spring-summer layer enriched in siliciclastic minerals and a lighter clay-size calcite-rich winter layer. Three subtypes of varves are distinguished: Type I with a 1–1.5 mm fining upward sequence; Type II with a distinct fine-sand base up to 3 mm thick; and Type III containing multiple internal microlaminae caused by individual summer rainstorm deposits. Delta-fan surface samples and sediment trap data fingerprint different sediment source areas and transport processes from the watershed and confirm the instant response of sediment flux to rainfall and erosion. Based on a highly accurate, precise and reproducible chronology, we demonstrate that sediment accumulation (varve thickness) is a quantitative predictor for cumulative boreal alpine spring (May–June) and spring/summer (May–August) rainfall (rMJ = 0.71, rMJJA = 0.60, p < 0.01). Bootstrap-based verification of the calibration model reveals a root mean squared error of prediction (RMSEPMJ = 32.7 mm, RMSEPMJJA = 57.8 mm) which is on the order of 10–13 % of mean MJ and MJJA cumulative precipitation, respectively. These results highlight the potential of the Lake Oeschinen sediments for high-resolution reconstructions of past rainfall conditions in the northern Swiss Alps, central and eastern France and south-west Germany.

Patterns and controls of sediment production, transfer and yield in the Illgraben

Quantification of the volumes of sediment removed by rock–slope failure and debris flows and identification of their coupling and controls are pertinent to understanding mountain basin sediment yield and landscape evolution. This study captures a multi-decadal period of hillslope erosion and channel change following an extreme rock avalanche in 1961 in the Illgraben, a catchment prone to debris flows in the Swiss Alps. We analyzed photogrammetrically-derived datasets of hillslope and channel erosion and deposition along with climatic and seismic variables for a 43 year period from 1963 to 2005. Based on these analyses we identify and discuss (1) patterns of hillslope production, channel transfer and catchment sediment yield, (2) their dominant interactions with climatic and seismic variables, and (3) the nature of hillslope–channel coupling and implications for sediment yield and landscape evolution in this mountain basin. Our results show an increase in the mean hillslope erosion rate in the 1980s from 0.24 ± 0.01 m yr− 1 to 0.42 ± 0.03 m yr− 1 that coincided with a significant increase in air temperature and decrease in snow cover depth and duration, which we presume led to an increase in the exposure of the slopes to thermal weathering processes. The combination of highly fractured slopes close to the threshold angle for failure, and multiple potential triggering mechanisms, means that it is difficult to identify an individual control on slope failure. On the other hand, the rate of channel change was strongly related to variables influencing runoff. A period of particularly high channel erosion rate of 0.74 ± 0.02 m yr− 1 (1992–1998) coincided with an increase in the frequency and magnitude of intense rainfall events. Hillslope erosion exceeded channel erosion on average, indicative of a downslope-directed coupling relationship between hillslope and channel, and demonstrating the first order control of rock–slope failure on catchment sediment yield and landscape evolution.

Characterisation of flash floods in small ungauged mountain basins of Central Spain using an integrated approach

One of the main problems of flood hazard assessment in ungauged or poorly gauged basins is the lack of runoff data. In an attempt to overcome this problem we have combined archival records, dendrogeomorphic time series and instrumental data (daily rainfall and discharge) from four ungauged and poorly gauged mountain basins in Central Spain with the aim of reconstructing and compiling information on 41 flash flood events since the end of the 19th century. Estimation of historical discharge and the incorporation of uncertainty for the at-site and regional flood frequency analysis were performed with an empirical rainfall–runoff assessment as well as stochastic and Bayesian Markov Chain Monte Carlo (MCMC) approaches. Results for each of the ungauged basins include flood frequency, severity, seasonality and triggers (synoptic meteorological situations). The reconstructed data series clearly demonstrates how uncertainty can be reduced by including historical information, but also points to the considerable influence of different approaches on quantile estimation. This uncertainty should be taken into account when these data are used for flood risk management.

Rainfall input, throughfall and stemflow of nutrients in a central African rain forest dominated by ectomycorrhizal trees

Incident rainfall is a major source of nutrient input to a forest ecosystem and the consequent throughfall and stemflow contribute to nutrient cycling. These rain-based fluxes were measured over 12 mo in two forest types in Korup National Park, Cameroon, one with low (LEM) and one with high (HEM) ectomycorrhizal abundances of trees. Throughfall was 96.6 and 92.4% of the incident annual rainfall (5370 mm) in LEM and HEM forests respectively; stemflow was correspondingly 1.5 and 2.2%. Architectural analysis showed that ln(funneling ratio) declined linearly with increasing ln(basal area) of trees. Mean annual inputs of N, P, K, Mg and Ca in incident rainfall were 1.50, 1.07, 7.77, 5.25 and 9.27 kg ha(-1), and total rain-based inputs to the forest floor were 5.0, 3.2, 123.4, 14.4 and 37.7 kg ha-1 respectively. The value for K is high for tropical forests and that for N is low. Nitrogen showed a significantly lower loading of throughfall and stemflow in HEM than in LEM forest, this being associated in the HEM forest with a greater abundance of epiphytic bryophytes which may absorb more N. Incident rainfall provided c. 35% of the gross input of P to the forest floor (i. e., rain-based plus small litter inputs), a surprisingly high contribution given the sandy P-poor soils. At the start of the wet season leaching of K from the canopy was particularly high. Calcium in the rain was also highest at this time, most likely due to washing off of dry-deposited Harmattan dusts. It is proposed that throughfall has an important `priming' function in the rapid decomposition of litter and mineralization of P at the start of the wet season. The contribution of P inputted from the atmosphere appears to be significant when compared to the rates of P mineralization from leaf litter.

Dynamics of the last four glacial terminations recorded in Lake Van, Turkey

A well-dated suite of Lake Van climate-proxy data covering the last 360 ka documents environmental changes over 4 glacial/interglacial cycles in Eastern Anatolia, Turkey. The picture of cold and dry glacials and warm and wet interglacials emerging from pollen, organic carbon, authigenic carbonate content, elemental profiling by XRF and lithological analyses is inconsistent with classical interpretation of ox- ygen isotopic composition of carbonates pointing to a more complex pattern in Lake Van region. Detailed analysis of glacial terminations allows for the constraining of a depositional model explaining different patterns observed in all the proxies. We hypothesize that variations in relative contribution of rainfall, snowmelt and glacier meltwater recharging the basin have a very important role for all sedimentary processes in Lake Van. Lake level of glacial Lake Van, predominantly fed by snowmelt, was low, the water column was oxic, and carbonates precipitating in the epilimnion recorded the light isotopic signature of inflow. During terminations, increasing rainfall and significant supply of mountain glaciers' meltwater contributed to lake level rise. Increased rainfall enhanced density gradients in the water column, and hindered mixing leading to development of bottom-water anoxia. Carbonates precipitating during terminations show large fluctuations in their isotopic composition. Full interglacial conditions in Lake Van are characterized by high or slowly falling lake level. Rainfall and snowmelt feed the lake but due to re-established mixing, the isotopic composition of authigenic carbonates is heavier and closer to that of evaporation-influenced lake water than that of runoff representing snowmelt and atmospheric precipitation.

Retrospective analysis of a nonforecasted rain-on-snow flood in the Alps – a matter of model limitations or unpredictable nature?

A rain-on-snow flood occurred in the Bernese Alps, Switzerland, on 10 October 2011, and caused significant damage. As the flood peak was unpredicted by the flood forecast system, questions were raised concerning the causes and the predictability of the event. Here, we aimed to reconstruct the anatomy of this rain-on-snow flood in the Lötschen Valley (160 km2) by analyzing meteorological data from the synoptic to the local scale and by reproducing the flood peak with the hydrological model WaSiM-ETH (Water Flow and Balance Simulation Model). This in order to gain process understanding and to evaluate the predictability. The atmospheric drivers of this rain-on-snow flood were (i) sustained snowfall followed by (ii) the passage of an atmospheric river bringing warm and moist air towards the Alps. As a result, intensive rainfall (average of 100 mm day-1) was accompanied by a temperature increase that shifted the 0° line from 1500 to 3200 m a.s.l. (meters above sea level) in 24 h with a maximum increase of 9 K in 9 h. The south-facing slope of the valley received significantly more precipitation than the north-facing slope, leading to flooding only in tributaries along the south-facing slope. We hypothesized that the reason for this very local rainfall distribution was a cavity circulation combined with a seeder-feeder-cloud system enhancing local rainfall and snowmelt along the south-facing slope. By applying and considerably recalibrating the standard hydrological model setup, we proved that both latent and sensible heat fluxes were needed to reconstruct the snow cover dynamic, and that locally high-precipitation sums (160 mm in 12 h) were required to produce the estimated flood peak. However, to reproduce the rapid runoff responses during the event, we conceptually represent likely lateral flow dynamics within the snow cover causing the model to react "oversensitively" to meltwater. Driving the optimized model with COSMO (Consortium for Small-scale Modeling)-2 forecast data, we still failed to simulate the flood because COSMO-2 forecast data underestimated both the local precipitation peak and the temperature increase. Thus we conclude that this rain-on-snow flood was, in general, predictable, but requires a special hydrological model setup and extensive and locally precise meteorological input data. Although, this data quality may not be achieved with forecast data, an additional model with a specific rain-on-snow configuration can provide useful information when rain-on-snow events are likely to occur.

Ecosystem functions are resistant to extreme changes to rainfall regimes in a mesotrophic grassland

Major changes to rainfall regimes are predicted for the future but the effect of such changes on terrestrial ecosystem function is largely unknown. We created a rainfall manipulation experiment to investigate the effects of extreme changes in rainfall regimes on ecosystem functioning in a grassland system. We applied two rainfall regimes; a prolonged drought treatment (30 % reduction over spring and summer) and drought/downpour treatment (long periods of no rainfall interspersed with downpours), with an ambient control. Both rainfall manipulations included increased winter rainfall. We measured plant community composition, CO2 fluxes and soil nutrient availability. Plant species richness and cover were lower in the drought/downpour treatment, and showed little recovery after the treatment ceased. Ecosystem processes were less affected, possibly due to winter rainfall additions buffering reduced summer rainfall, which saw relatively small soil moisture changes. However, soil extractable P and ecosystem respiration were significantly higher in rainfall change treatments than in the control. This grassland appears fairly resistant, in the short term, to even the more extreme rainfall changes that are predicted for the region, although prolonged study is needed to measure longer-term impacts. Differences in ecosystem responses between the two treatments emphasise the variety of ecosystem responses to changes in both the size and frequency of rainfall events. Given that model predictions are inconsistent there is therefore a need to assess ecosystem function under a range of potential climate change scenarios.

Combined effect of soil bund with biological soil and water conservation measures in the northwestern Ethiopian highlands

Excessive runoff and soil erosion in the upper Blue Nile Basin poses a threat that has attracted the attention of the Ethiopian government because of the serious on-site effects in addition to downstream effects, such as the siltation of water harvesting structures and reservoirs. The objective of the study was to evaluate and recommend effective biophysical soil and water conservation measure(s) in the Debre Mewi watershed, about 30 km south of the Lake Tana. Six conservation measures were evaluated for their effects on runoff, soil loss, and forage yield using runoff plots. There was a significant difference between treatments for both runoff and soil loss. The four-year average annual soil loss in the different plots ranged from 26 to 71 t ha−1, and total runoff ranged from 180 to 302 mm, while annual rainfall varied between 854 mm in 2008 and 1247 mm in 2011. Soil bund combined with elephant grass had the lowest runoff and soil loss as compared to the other treatments, whereas the untreated control plot had the highest for both parameters. As an additional benefit, 2.8 and 0.7 t ha−1 year−1 of dried forage was obtained from elephant and local grasses, respectively. Furthermore, it was found that soil bund combined with Tephrosia increased soil organic matter by 13% compared to the control plot. Soil bund efficiency was significantly enhanced by combining them with biological measures and improved farmers’ perception of soil and water conservation measures.

Plant and arthropod community sensitivity to rainfall manipulation but not nitrogen enrichment in a successional grassland ecosystem

Grasslands provide many ecosystem services including carbon storage, biodiversity preservation and livestock forage production. These ecosystem services will change in the future in response to multiple global environmental changes, including climate change and increased nitrogen inputs. We conducted an experimental study over 3 years in a mesotrophic grassland ecosystem in southern England. We aimed to expose plots to rainfall manipulation that simulated IPCC 4th Assessment projections for 2100 (+15 % winter rainfall and −30 % summer rainfall) or ambient climate, achieving +15 % winter rainfall and −39 % summer rainfall in rainfall-manipulated plots. Nitrogen (40 kg ha−1 year−1) was also added to half of the experimental plots in factorial combination. Plant species composition and above ground biomass were not affected by rainfall in the first 2 years and the plant community did not respond to nitrogen enrichment throughout the experiment. In the third year, above-ground plant biomass declined in rainfall-manipulated plots, driven by a decline in the abundances of grass species characteristic of moist soils. Declining plant biomass was also associated with changes to arthropod communities, with lower abundances of plant-feeding Auchenorrhyncha and carnivorous Araneae indicating multi-trophic responses to rainfall manipulation. Plant and arthropod community composition and plant biomass responses to rainfall manipulation were not modified by nitrogen enrichment, which was not expected, but may have resulted from prior nitrogen saturation and/or phosphorus limitation. Overall, our study demonstrates that climate change may in future influence plant productivity and induce multi-trophic responses in grasslands.