Local models for exploratory analysis of hydrological extremes

Trend analysis is widely used for detecting changes in hydrological data. Parametric methods for this employ pre-specified models and associated tests to assess significance, whereas non-parametric methods generally apply rank tests to the data. Neither approach is suitable for exploratory analysis, because parametric models impose a particular, perhaps unsuitable, form of trend, while testing may confirm that trend is present but does not describe its form. This paper describes semi-parametric approaches to trend analysis using local likelihood fitting of annual maximum and partial duration series and illustrates their application to the exploratory analysis of changes in extremes in sea level and river flow data. Bootstrap methods are used to quantify the variability of estimates.

Long-term variability and environmental preferences of calycophoran siphonophores in the Bay of Villefranche (North-Western Mediterranean)

Long-term variability of the main calycophoran siphonophores was investigated between 1974 and 1999 in a coastal station in the north-western Mediterranean. The data were collected at weekly frequency using a macroplankton net (680 μm mesh size) adapted to quantitatively sample delicate gelatinous plankton. A 3-year collection (1967–1969) of siphonophores from offshore waters using the same methodology showed that the patterns of variability observed inshore were representative of siphonophores’ changes at a regional scale. The aims of the study were: (i) to investigate the patterns of variability that characterised the dominant calycophoran species and assemblages; (ii) to identify the environmental optima that were associated with a significant increase in the dominant siphonophore species and (iii) to verify the influence of hydroclimatic variability on long-term changes of siphonophores. Our results showed that during nearly 3 decades the standing stock of calycophoran siphonophores did not show any significant change, with the annual maximum usually recorded in spring as a result of high densities of the dominant species Lensia subtilis, Muggiaea kochi and Muggiaea atlantica. Nevertheless, major changes in community composition occurred within the calycophoran population. Since the middle 1980s, M. kochi, once the most dominant species, started to decrease allowing other species, the congeneric M. atlantica and Chelophyes appendiculata, to increasingly dominate in spring and summer–autumn, respectively. The comparison of environmental and biotic long-term trends suggests that the decrease of M. kochi was triggered by hydrological changes that occurred in the north-western Mediterranean under the forcing of large-scale climate oscillations. Salinity, water stratification and water temperature were the main hydroclimatic factors associated with a significant increase of siphonophores, different species showing different environmental preferences.

Temporal variability in total, micro- and nano-phytoplankton primary production at a coastal site in the Western English Channel

Primary productivity and subsequent carbon cycling in the coastal zone have a significant impact on the global carbon budget. It is currently unclear how anthropogenic activity could alter these budgets but long term coastal time series of hydrological, biogeochemical and biological measurements represent a key means to better understand past drivers, and hence to predicting future seasonal and inter-annual variability in carbon fixation in coastal ecosystems. An 8-year time series of primary production from 2003 to 2010, estimated using a recently developed absorption-based algorithm, was used to determine the nature and extent of change in primary production at a coastal station (L4) in the Western English Channel (WEC). Analysis of the seasonal and inter-annual variability in production demonstrated that on average, nano- and pico-phytoplankton account for 48% of the total carbon fixation and micro-phytoplankton for 52%. A recent decline in the primary production of nano- and pico-phytoplankton from 2005 to 2010 was observed, corresponding with a decrease in winter nutrient concentrations and a decrease in the biomass of Phaeocystis sp. Micro-phytoplankton primary production (PPM) remained relatively constant over the time series and was enhanced in summer during periods of high precipitation. Increases in sea surface temperature, and decreases in wind speeds and salinity were associated with later spring maxima in PPM. Together these trends indicate that predicted increases in temperature and decrease in wind speeds in future would drive later spring production whilst predicted increases in precipitation would also continue these blooms throughout the summer at this site.

Long-term hydroclimatic change and interannual variability in water sources, Apex River (Iqaluit), Baffin Island, Nunavut

The eastern Canadian Arctic is home to Canada’s largest Indigenous population, which depends on local freshwater sources for drinking water. However, small watersheds have rarely been analyzed for long-term hydrologic response to changing climate. This study aims to address this issue by examining the Apex River, a small watershed with a long hydroclimatic record, near Iqaluit, Nunavut. Particular emphasis was placed on the long-term changes in climate and river discharge, and the seasonal variability of water sources between two snapshots in time, 1983 and 2013. Long-term hydrological data were obtained from gauge station 10UH002, operated by Environment and Climate Change Canada, and long-term meteorological data were acquired from Environment Canada–operated stations near Iqaluit Airport. Breakpoint analysis suggested that long-term mean annual surface air temperatures have increased since 1994. In contrast, no long-term total precipitation or annual discharge changes were observed. However, river flow initiation and cessation analyses of the Apex River flow season indicates that flow extended into the autumn since the 2000s. The 2013 flow season lasted 44 days longer than the 1983 flow season. Systematic river sampling was undertaken throughout the 2013 thaw season to determine contributing proportions of event (snowmelt or rainfall) and pre-event (baseflow) water to river runoff. Results from the stable isotope hydrograph separation for 2013 were compared to findings for 1983. Snow was the main source of water to the river during the snowmelt period in 1983 and 2013, however baseflow was still an important contributor. Although there was high similarity of water sources early in the season in 1983 and 2013, the two years differed during the autumn. In 2013 there was a high rainfall runoff response that was not present in 1983, suggesting high release of late-season sub-surface water storage and an increased sensitivity to late-season rainfall events in 2013. This research provides insights into the hydrologic response of the Apex River to long-term climatic change, and highlights the need for high-quality precipitation and discharge data for effective long-term hydrological assessment.

Mid- to late-Holocene climate variability and anthropogenic impacts: multi-proxy evidence from Lake Bliden, Denmark

We conducted multi-proxy geochemical analyses (including measurements of organic carbon, nitrogen and sulphur stable isotope composition, and carbonate carbon and oxygen isotope composition) on a 13.5 m sediment core from Lake Bliden, Denmark, which provide a record of shifting hydrological conditions for the past 6,700 years. The early part of the stratigraphic record (6,700-5,740 cal year BP) was wet, based on delta O-18(carb) and lithology, and corresponds to the Holocene Thermal Maximum. Shifts in primarily delta O-18(carb) indicate dry conditions prevailed from 5,740 to 2,800 cal year BP, although this was interrupted by very wet conditions from 5,300 to 5,150, 4,300 to 4,050 and 3,700 to 3,450 cal year BP. The timing of the latter two moist intervals is consistent with other Scandinavian paleoclimatic records. Dry conditions at Lake Bliden between 3,450 and 2,800 cal year BP is consistent with other paleolimnological records from southern Sweden but contrasts with records in central Sweden, possibly suggesting a more northerly trajectory of prevailing westerlies carrying moisture from the North Atlantic at this time. Overall, fluctuating moisture conditions at Lake Bliden appear to be strongly linked to changing sea surface temperatures in the Greenland, Iceland and Norwegian seas. After 2,800 cal year BP, sedimentology, magnetic susceptibility, delta C-13(ORG), delta C-13(carb) and delta O-18(carb) indicate a major reduction on water level, which caused the depositional setting at the coring site to shift from the profundal to littoral zone. The Roman Warm Period (2,200-1,500 cal year BP) appears dry based on enriched delta O-18(carb) values. Possible effects of human disturbance in the watershed after 820 cal year BP complicate attempts to interpret the stratigraphic record although tentative interpretation of the delta O-18(carb), magnetic susceptibility, delta C-13(ORG), delta C-13(carb) and delta O-18(carb) records suggest that the Medieval Warm Period was dry and the Little Ice Age was wet.

Interannual and interdecadal oscillations in hydrological variables: Sources and modeling of the persistence in the Elbe River Basin

Wenn man die Existenz von physikalischen Mechanismen ignoriert, die für die Struktur hydrologischer Zeitreihen verantwortlich sind, kann das zu falschen Schlussfolgerungen bzgl. des Vorhandenseins möglicher Gedächtnis (memory) -Effekte, d.h. von Persistenz, führen. Die hier vorgelegte Doktorarbeit spürt der niedrigfrequenten klimatischen Variabilität innerhalb den hydrologischen Zyklus nach und bietet auf dieser "Reise" neue Einsichten in die Transformation der charakteristischen Eigenschaften von Zeitreihen mit einem Langzeitgedächtnis. Diese Studie vereint statistische Methoden der Zeitreihenanalyse mit empirisch-basierten Modelltechniken, um operative Modelle zu entwickeln, die in der Lage sind (1) die Dynamik des Abflusses zu modellieren, (2) sein zukünftiges Verhalten zu prognostizieren und (3) die Abflusszeitreihen an unbeobachteten Stellen abzuschätzen. Als solches präsentiert die hier vorgelegte Dissertation eine ausführliche Untersuchung zu den Ursachen der niedrigfrequenten Variabilität von hydrologischen Zeitreihen im deutschen Teil des Elbe-Einzugsgebietes, den Folgen dieser Variabilität und den physikalisch basierten Reaktionen von Oberflächen- und Grundwassermodellen auf die niedrigfrequenten Niederschlags-Eingangsganglinien. Die Doktorarbeit gliedert sich wie folgt: In Kapitel 1 wird als Hintergrundinformation das Hurst Phänomen beschrieben und ein kurzer Rückblick auf diesbezügliche Studien gegeben. Das Kapitel 2 diskutiert den Einfluss der Präsenz von niedrigfrequenten periodischen Zeitreihen auf die Zuverlässigkeit verschiedener Hurst-Parameter-Schätztechniken. Kapitel 3 korreliert die niedrigfrequente Niederschlagsvariabilität mit dem Index der Nord-Atlantischen Ozillations (NAO). Kapitel 4-6 sind auf den deutschen Teil des Elbe-Einzugsgebietes fokussiert. So werden in Kapitel 4 die niedrigfrequenten Variabilitäten der unterschiedlichen hydro-meteorologischen Parameter untersucht und es werden Modelle beschrieben, die die Dynamik dieser Niedrigfrequenzen und deren zukünftiges Verhalten simulieren. Kapitel 5 diskutiert die mögliche Anwendung der Ergebnisse für die charakteristische Skalen und die Verfahren der Analyse der zeitlichen Variabilität auf praktische Fragestellungen im Wasserbau sowie auf die zeitliche Bestimmung des Gebiets-Abflusses an unbeobachteten Stellen. Kapitel 6 verfolgt die Spur der Niedrigfrequenzzyklen im Niederschlag durch die einzelnen Komponenten des hydrologischen Zyklus, nämlich dem Direktabfluss, dem Basisabfluss, der Grundwasserströmung und dem Gebiets-Abfluss durch empirische Modellierung. Die Schlussfolgerungen werden im Kapitel 7 präsentiert. In einem Anhang werden technische Einzelheiten zu den verwendeten statistischen Methoden und die entwickelten Software-Tools beschrieben.

Variability in clear-sky longwave radiative cooling of the atmosphere

The longwave radiative cooling of the clear-sky atmosphere (Q(LWc)) is a crucial component of the global hydrological cycle and is composed of the clear-sky outgoing longwave radiation to space (OLRc) and the net downward minus upward clear-sky longwave radiation to the surface (SNLc). Estimates of QLWc from reanalyses and observations are presented for the period 1979-2004. Compared to other reanalyses data sets, the European Centre for Medium-range Weather Forecasts 40-year reanalysis (ERA40) produces the largest Q(LWc) over the tropical oceans (217 W m(-2)), explained by the least negative SNLc. On the basis of comparisons with data derived from satellite measurements, ERA40 provides the most realistic QLWc climatology over the tropical oceans but exhibits a spurious interannual variability for column integrated water vapor (CWV) and SNLc. Interannual monthly anomalies of QLWc are broadly consistent between data sets with large increases during the warm El Nino events. Since relative humidity ( RH) errors applying throughout the troposphere result in compensating effects on the cooling to space and to the surface, they exert only a marginal effect on QLWc. An observed increase in CWV with surface temperature of 3 kg m(-2) K-1 over the tropical oceans is important in explaining a positive relationship between QLWc and surface temperature, in particular over ascending regimes; over tropical ocean descending regions this relationship ranges from 3.6 to 4.6 +/- 0.4 W m(-2) K-1 for the data sets considered, consistent with idealized sensitivity tests in which tropospheric warming is applied and RH is held constant and implying an increase in precipitation with warming.

Fine sediment delivery and transfer in lowland catchments: modelling suspended sediment concentrations in response to hydrological forcing

Fine sediment delivery to and storage in stream channel reaches can disrupt aquatic habitats, impact river hydromorphology, and transfer adsorbed nutrients and pollutants from catchment slopes to the fluvial system. This paper presents a modelling toot for simulating the time-dependent response of the fine sediment system in catchments, using an integrated approach that incorporates both land phase and in-stream processes of sediment generation, storage and transfer. The performance of the model is demonstrated by applying it to simulate in-stream suspended sediment concentrations in two lowland catchments in southern England, the Enborne and the Lambourn, which exhibit contrasting hydrological and sediment responses due to differences in substrate permeability. The sediment model performs well in the Enborne catchment, where direct runoff events are frequent and peak suspended sediment concentrations can exceed 600 mg l(-1). The general trends in the in-stream concentrations in the Lambourn catchment are also reproduced by the model, although the observed concentrations are low (rarely exceeding 50 mg l(-1)) and the background variability in the concentrations is not fully characterized by the model. Direct runoff events are rare in this highly permeable catchment, resulting in a weak coupling between the sediment delivery system and the catchment hydrology. The generic performance of the model is also assessed using a generalized sensitivity analysis based on the parameter bounds identified in the catchment applications. Results indicate that the hydrological parameters contributing to the sediment response include those controlling (1) the partitioning of runoff between surface and soil zone flows and (2) the fractional loss of direct runoff volume prior to channel delivery. The principal sediment processes controlling model behaviour in the simulations are the transport capacity of direct runoff and the in-stream generation, storage and release of the fine sediment fraction. The in-stream processes appear to be important in maintaining the suspended sediment concentrations during low flows in the River Enborne and throughout much of the year in the River Lambourn. Copyright (c) 2007 John Wiley & Sons, Ltd.

Simulating current global river runoff with a global hydrological model: model revisions, validation, and sensitivity analysis

Global hydrological models (GHMs) model the land surface hydrologic dynamics of continental-scale river basins. Here we describe one such GHM, the Macro-scale - Probability-Distributed Moisture model.09 (Mac-PDM.09). The model has undergone a number of revisions since it was last applied in the hydrological literature. This paper serves to provide a detailed description of the latest version of the model. The main revisions include the following: (1) the ability for the model to be run for n repetitions, which provides more robust estimates of extreme hydrological behaviour, (2) the ability of the model to use a gridded field of coefficient of variation (CV) of daily rainfall for the stochastic disaggregation of monthly precipitation to daily precipitation, and (3) the model can now be forced with daily input climate data as well as monthly input climate data. We demonstrate the effects that each of these three revisions has on simulated runoff relative to before the revisions were applied. Importantly, we show that when Mac-PDM.09 is forced with monthly input data, it results in a negative runoff bias relative to when daily forcings are applied, for regions of the globe where the day-to-day variability in relative humidity is high. The runoff bias can be up to - 80% for a small selection of catchments but the absolute magnitude of the bias may be small. As such, we recommend future applications of Mac-PDM.09 that use monthly climate forcings acknowledge the bias as a limitation of the model. The performance of Mac-PDM.09 is evaluated by validating simulated runoff against observed runoff for 50 catchments. We also present a sensitivity analysis that demonstrates that simulated runoff is considerably more sensitive to method of PE calculation than to perturbations in soil moisture and field capacity parameters.

Decadal-scale rainfall variability in Ethiopia recorded in an annually laminated, Holocene-age, stalagmite

An annually laminated, uranium-series dated, Holocene stalagmite from southeast Ethiopia has been analysed for growth rate and δ13C and δ18O variations at annual to biennial resolution, in order to provide the first long duration proxy record of decadal-scale rainfall variability in this climatically sensitive region. Our study site (10°N) is climatically influenced by both summer (June—August) and spring (March—May) rainfall caused by the annual movement of the Inter-Tropical Convergence Zone (ITCZ) and modulated by large-scale anomalies in the atmospheric circulation and in ocean temperatures. Here we show that stalagmite growth, episodic throughout the last 7800 years, demonstrates decadal-scale (8—25 yr) variability in both growth rate and δ 18O. A hydrological model was employed and indicates that this decadal variability is due to variations in the relative amounts of rainfall in the two rain seasons. Our record, unique in its combination of length (a total of ~1000 years), annual chronology and high resolution δ18O, shows for the first time that such decadal-scale variability in rainfall in this region has occurred through the Holocene, which implies persistent decadal-scale variability for the large-scale atmospheric and oceanic driving factors.

A molluscan perspective on hydrological cycle dynamics in northwestern Europe

Shell aragonite δ18O values of unionid freshwater mussels are applied as a proxy for past river discharges in the rivers Rhine and Meuse, using a set of nine shells from selected climatic intervals during the late Holocene. A single Meuse shell derives from the Subboreal and its δ18O values are similar to modern values. The Rhine specimens represent the Subboreal, the Roman Warm Period and the Medieval Warm Period (MWP). These shells also show averages and ranges of aragonite δ18O values similar to modern specimens. This indicates that environmental conditions such as Rhine river dynamics, Alpine meltwater input and drought severity during these intervals were similar to the 20th century. These shells do not record subtle centennial to millennial climatic variation due to their relatively short lifespan and the large inter-annual and intra-seasonal variation in environmental conditions. However, they are very suitable for studying seasonal to decadal scale climate variability. The two shells with the longest lifespan appear to show decadal scale variability in reconstructed water δ18O values during the MWP, possibly forced by the North Atlantic Oscillation (NAO), which is the dominant mode of variability influencing precipitation regimes over Europe.

Uncertainty in the relationship between climate forcing and hydrological response in UK catchments

This paper assesses the relationship between amount of climate forcing – as indexed by global mean temperature change – and hydrological response in a sample of UK catchments. It constructs climate scenarios representing different changes in global mean temperature from an ensemble of 21 climate models assessed in the IPCC AR4. The results show a considerable range in impact between the 21 climate models, with – for example - change in summer runoff at a 2oC increase in global mean temperature varying between -40% and +20%. There is evidence of clustering in the results, particularly in projected changes in summer runoff and indicators of low flows, implying that the ensemble mean is not an appropriate generalised indicator of impact, and that the standard deviation of responses does not adequately characterise uncertainty. The uncertainty in hydrological impact is therefore best characterised by considering the shape of the distribution of responses across multiple climate scenarios. For some climate model patterns, and some catchments, there is also evidence that linear climate change forcings produce non-linear hydrological impacts. For most variables and catchments, the effects of climate change are apparent above the effects of natural multi-decadal variability with an increase in global mean temperature above 1oC, but there are differences between catchments. Based on the scenarios represented in the ensemble, the effect of climate change in northern upland catchments will be seen soonest in indicators of high flows, but in southern catchments effects will be apparent soonest in measures of summer and low flows. The uncertainty in response between different climate model patterns is considerably greater than the range due to uncertainty in hydrological model parameterisation.

Hydrological cycle in the Danube basin in present-day and XXII century simulations by IPCCAR4 global climate models

We present an intercomparison and verification analysis of 20 GCMs (Global Circulation Models) included in the 4th IPCC assessment report regarding their representation of the hydrological cycle on the Danube river basin for 1961–2000 and for the 2161–2200 SRESA1B scenario runs. The basin-scale properties of the hydrological cycle are computed by spatially integrating the precipitation, evaporation, and runoff fields using the Voronoi-Thiessen tessellation formalism. The span of the model- simulated mean annual water balances is of the same order of magnitude of the observed Danube discharge of the Delta; the true value is within the range simulated by the models. Some land components seem to have deficiencies since there are cases of violation of water conservation when annual means are considered. The overall performance and the degree of agreement of the GCMs are comparable to those of the RCMs (Regional Climate Models) analyzed in a previous work, in spite of the much higher resolution and common nesting of the RCMs. The reanalyses are shown to feature several inconsistencies and cannot be used as a verification benchmark for the hydrological cycle in the Danubian region. In the scenario runs, for basically all models the water balance decreases, whereas its interannual variability increases. Changes in the strength of the hydrological cycle are not consistent among models: it is confirmed that capturing the impact of climate change on the hydrological cycle is not an easy task over land areas. Moreover, in several cases we find that qualitatively different behaviors emerge among the models: the ensemble mean does not represent any sort of average model, and often it falls between the models’ clusters.

Physically consistent responses of the global atmospheric hydrological cycle in models and observations

Robust and physically understandable responses of the global atmospheric water cycle to a warming climate are presented. By considering interannual responses to changes in surface temperature (T), observations and AMIP5 simulations agree on an increase in column integrated water vapor at the rate 7 %/K (in line with the Clausius­Clapeyron equation) and of precipitation at the rate 2-­3 %/K (in line with energetic constraints). Using simple and complex climate models, we demonstrate that radiative forcing by greenhouse gases is currently suppressing global precipitation (P) at ~ -0.15 %/decade. Along with natural variability, this can explain why observed trends in global P over the period 1988-2008 are close to zero. Regional responses in the global water cycle are strongly constrained by changes in moisture fluxes. Model simulations show an increased moisture flux into the tropical wet region at 900 hPa and an enhanced outflow (of smaller magnitude) at around 600 hPa with warming. Moisture transport explains an increase in P in the wet tropical regions and small or negative changes in the dry regions of the subtropics in CMIP5 simulations of a warming climate. For AMIP5 simulations and satellite observations, the heaviest 5-day rainfall totals increase in intensity at ~15 %/K over the ocean with reductions at all percentiles over land. The climate change response in CMIP5 simulations shows consistent increases in P over ocean and land for the highest intensities, close to the Clausius-Clapeyron scaling of 7 %/K, while P declines for the lowest percentiles, indicating that interannual variability over land may not be a good proxy for climate change. The local changes in precipitation and its extremes are highly dependent upon small shifts in the large-scale atmospheric circulation and regional feedbacks.

In situ preservation of wetland heritage: hydrological and chemical change in the burial environment of the Somerset Levels, UK.

In Situ preservation is a core strategy for the conservation and management of waterlogged remains at wetland sites. Inorganic and organic remains can, however, quickly become degraded, or lost entirely, as a result of chemical or hydrological changes. Monitoring is therefore crucial in identifying baseline data for a site, the extent of spatial and or temporal variability, and in evaluating the potential impacts of these variables on current and future In Situ preservation potential. Since August 2009, monthly monitoring has taken place at the internationally important Iron Age site of Glastonbury Lake Village in the Somerset Levels, UK. A spatial, stratigraphic, and analytical approach to the analysis of sediment horizons and monitoring of groundwater chemistry, redox potential, water table depth and soil moisture (using TDR) was used to characterize the site. Significant spatial and temporal variability has been identified, with results from water-table monitoring and some initial chemical analysis from Glastonbury presented here. It appears that during dry periods parts of this site are at risk from desiccation. Analysis of the chemical data, in addition to integrating the results from the other parameters, is ongoing, with the aim of clarifying the risk to the entire site.