Space–time modelling of catchment scale drought characteristics

Drought may affect all components of the water cycle and covers commonly a large part of the catchment area. This paper examines drought propagation at the catchment scale using spatially aggregated drought characteristics and illustrates the importance of catchment processes in modifying the drought signal in both time and space. Analysis is conducted using monthly time series covering the period 1961–1997 for the Pang catchment, UK. The time series include observed rainfall and groundwater recharge, head and discharge simulated by physically-based soil water and groundwater models. Drought events derived separately for each unit area and variable are combined to yield catchment scale drought characteristics. The study reveals relatively large differences in the spatial and temporal characteristics of drought for the different variables. Meteorological droughts cover frequently the whole catchment; and they are more numerous and last for a short time (1–2 months). In comparison, droughts in recharge and hydraulic head cover typically a smaller area and last longer (4–5 months). Hydraulic head and groundwater discharge exhibit similar drought characteristics, which can be expected in a groundwater fed catchment. Deficit volume is considered a robust measure of the severity of a drought event over the catchment area for all variables; whereas, duration is less sensitive, particular for rainfall. Spatial variability in drought characteristics for groundwater recharge, head and discharge are primarily controlled by catchment properties. It is recommended not to use drought area separately as a measure of drought severity at the catchment scale, rather it should be used in combination with other drought characteristics like duration and deficit volume.     

Estimating mean residence time of karst groundwater in mountainous catchments of Western Himalaya,India

The mean residence time (MRT) of karst groundwater in three mountainous catchments of the Western Himalaya was estimated using multiple approaches: the tritium method, the sine wave model and tracer tests. Water samples were collected from precipitation, glacier melt, streams and karst springs for δ²H and tritium analysis during 2012 and 2013. High tritium values were observed in winter precipitation and low values in summer precipitation. The variation of tritium in karst springs was similar to that of the streams, whereas glacier melt showed lower tritium values. The MRT of cold karst springs was shorter than that of warm karst springs. The tracer breakthrough curves (TBC) retrieved for different springs suggested a short travel time for groundwater and possibly conduit flow. Deterioration of water quality and variation in flux magnitude are the two main practical consequences of the short travel time of karst groundwater in the region.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR K. Heal     

Estimation of mean residence times of subsurface waters using seasonal variation in deuterium excess in a small headwater catchment in Japan

We measured deuterium excess (d = δD ? 8δ¹⁸O) in throughfall, groundwater, soil water, spring water, and stream water for 3 years in a small headwater catchment (Matsuzawa, 0·68 ha) in the Kiryu Experimental Watershed in Japan. The d value represents a kinetic effect produced when water evaporates. The d value of the throughfall showed a sinusoidal change (amplitude: 6·9‰ relative to Vienna standard mean ocean water (V‐SMOW)) derived from seasonal changes in the source of water vapour. The amplitude of this sinusoidal change was attenuated to 1·3–6·9‰ V‐SMOW in soil water, groundwater, spring water, and stream water. It is thought that these attenuations derive from hydrodynamic transport processes in the subsurface and mixing processes at an outflow point (stream or spring) or a well. The mean residence time (MRT) of water was estimated from d value variations using an exponential‐piston flow model and a dispersion model. MRTs for soil water were 0–5 months and were not necessarily proportional to the depth. This may imply the existence of bypass flow in the soil. Groundwater in the hillslope zone had short residence times, similar to those of the soil water. For groundwater in the saturated zone near the spring outflow point, the MRTs differed between shallow and deeper groundwater; shallow groundwater had a shorter residence time (5–8 months) than deeper groundwater (more than 9 months). The MRT of stream water (8–9 months) was between that of shallow groundwater near the spring and deeper groundwater near the spring. The seasonal variation in the d value of precipitation arises from changes in isotopic water vapour composition associated with seasonal activity of the Asian monsoon mechanism. The d value is probably an effective tracer for estimating the MRT of subsurface water not only in Japan, but also in other East Asian countries influenced by the Asian monsoon. Copyright © 2006 John Wiley & Sons, Ltd.     

Stable isotope tracers as diagnostic tools in upscaling flow path understanding and residence time estimates in a mountainous mesoscale catchment

δ¹⁸O measurements of precipitation and stream waters were used as a natural tracer to investigate hydrological pathways and residence times in the River Feshie, a complex mesoscale (231 km²) catchment in the Cairngorm Mountains of Scotland. Precipitation δ¹⁸O exhibited strong seasonal variation over the 2001–02 hydrological year, ranging from −6·9‰ in the summer, to −12·0‰ during winter snowfalls (mean δ¹⁸O −9·59‰). Although damped, this seasonality was reflected in stream water outputs at seven sampling sites in the catchment, allowing δ¹⁸O variations to be used to infer hydrological source areas. Thus, stream water δ¹⁸O was generally controlled by a seasonally variable storm flow end member, mixing with groundwater of more constant isotopic composition. Periodic regression analysis allowed the differences in this mixing process between monitoring subcatchments to be assessed more quantitatively to provide a preliminary estimate of mean stream water residence time. This demonstrated the importance of responsive hydrological pathways associated with peat and shallow alpine soils in the headwater subcatchments in producing seasonally variable runoff with short mean residence times (33–113 days). In contrast, other tributaries with more freely draining soils and larger groundwater storage in shallow aquifers provided more effective mixing of variable precipitation inputs, resulting in longer residence time estimates (178–445 days). The mean residence time of runoff leaving the Feshie catchment reflected an integration of these contrasting influences (110–200 days). These insights from δ¹⁸O measurements extend the hydrological understanding of the Feshie catchment gained from other hydrochemical tracers, and demonstrate the utility of isotope tracers in investigating hydrological processes at the mesoscale. Copyright © 2005 John Wiley & Sons, Ltd.     

Assessment of the adequacy of the regional relationships between catchment attributes and catchment response dynamics,calibrated by a multi‐objective approach

This study aimed to evaluate the effectiveness of the regionalization method on the basis of a combination of a parsimonious model structure and a multi‐objective calibration technique. For this study, 12 gauged catchments in the Republic of Korea were used. The parsimonious model structure, requiring minimal input data, was used to avoid adverse effects arising from model complexity, over‐parameterization and data requirements. The IHACRES rainfall‐runoff model was applied to represent the dynamic response characteristics of catchments in Korea. A multi‐objective approach was adopted to reduce the predictive uncertainty arising from the calibration of a rainfall‐runoff model, by increasing the amount of information retrieved from the available data. The regional relationships (or models) between the model parameters and the catchment attributes were established via a multiple regression approach, incorporating correlation analysis and stepwise regression on linear and logarithmic scales. The impacts of the parameters, calibrated by the multi‐objective approach, on the adequacy of regional relationships were assessed by comparison with impacts obtained by the single‐objective approach. The regional relationships were well defined, despite limited available data. The drainage area, the effective soil depth, the mean catchment slope and the catchment gradient appeared to be the main factors for describing the hydrologic response characteristics in the areas studied. The overall model performance of the regional models based on the multi‐objective approach was good, producing reasonable results for high and low flows and for the overall water balance, simultaneously. The regional models based on the single‐objective approach yielded accurate predictions in high flows but showed limited predictive capability for low flows and the overall water balance. This was due to the optimal model parameter estimates when using a single‐objective measure. The parameters calibrated by the single‐objective approach decreased the predictability of the regional models. Copyright © 2013 John Wiley & Sons, Ltd.     

Direct relationships of discharges to the mean and standard deviation of the intervals between their exceedences

Abstract

This article paves a way for assessing flood risk by the use of two-parameter distributions, for the intervals between threshold exceedences rather than by the traditional exponential distribution. In a case study, the apparent properties of intervals between exceedences of runoff events differ from those anticipated for exponentially distributed series. A procedure is proposed to relate two statistical parameters of the intervals to threshold discharges. It considers partial duration series (PDS) with thresholds equal to all high enough observed discharges. To avoid unnecessary assumptions on the behaviour of those parameters and effects of dependence between parameters for different PDS, a non-parametric trend-free pre-whitened scheme is applied. It leads to power-law relationships between a discharge and the mean and standard deviation of the intervals between its exceedences. Predicted mean inter-exceedence intervals, for the highest observed discharges at the stations, are closer to the observational periods than those predicted by GEV distributions fitted to AMS, and by GP distributions to fitted PDS. In the present case, the latter predictions are longer than the observational periods whereas some of the predicted mean inter-exceedences are shorter than the corresponding observational periods and some others are longer.

Citation Ben-Zvi, A. & Azmon, B. (2010) Direct relationships of discharges to the mean and standard deviation of the intervals between their exceedences. Hydrol. Sci J. 55(4), 565–577.     

Chloride in precipitation and streamwater for the upland catchment of river severn,mid-wales; some consequences for hydrochemical models

Variations in the concentration of Cl in rainfall and stream runoff are presented for two catchments in the Hafren forest of mid-Wales, Great Britain. Despite the large fluctuations in rainfall concentrations, Cl in the streamwater remains relatively constant. Using the two-reservoir Birkenes model, an attempt was made to simulate observed Cl in streamwater. The original model was unable to reproduce the observations and several modifications are suggested to provide better simulations. The resulting model is not the only one capable of reproducing the observations; other hydrochemical models will most probably also achieve this although emphasis will in each case be placed on different aspects. In this paper, it is suggested that the stochastic properties of water movement and chemical processes can account for the streamwater chemistry responses observed. On the catchment scale these processes will lead to an apparently deterministic behaviour that may well be described by simple relationships.     

The boundary layer of the residence time field

The residence time of a tracer in a control domain is usually computed by releasing tracer parcels and registering the time when each of these tracer parcels cross the boundary of the control domain. In this Lagrangian procedure, the particles are discarded or omitted as soon as they leave the control domain. In a Eulerian approach, the same approach can be implemented by integrating forward in time the advection–diffusion equation for a tracer. So far, the conditions to be applied at the boundary of the control domain were uncertain. We show here that it is necessary to prescribe that the tracer concentration vanishes at the boundary of the control domain to ensure the compatibility between the Lagrangian and Eulerian approaches. When we use the Constituent oriented Age and Residence time Theory (CART), this amounts to solving the differential equation for the residence time with boundary conditions forcing the residence time to vanish at the open boundaries of the control domain. Such boundary conditions are likely to induce the development of boundary layers (at outflow boundaries for the tracer concentration and at inflow boundaries for the residence time). The thickness of these boundary layers is of the order of the ratio of the diffusivity to the velocity. They can however be partly smoothed by tidal and other oscillating flows.     

Nitrate removal and young stream water fractions at the catchment scale

Despite extensive research on nitrate export and removal, nutrient contamination remains a major threat to water bodies worldwide. At the local scale, nitrate removal is governed by biogeochemical conditions that vary in space and time, making integration to entire landscapes critical. Water transit times have often been used to describe solute transport, but the relation between water age and nitrate removal at the catchment scale is still poorly understood. We test the hypothesis that nitrate removal peaks when the fraction of young water in discharge is at its minimum, because nitrate removal occurs mostly under dry conditions where deeper, older groundwater dominates streamflow. We tested this hypothesis by exploring a detailed water quality record from the Kervidy–Naizin catchment (FR) and comparing the dynamics of nitrate to those of a conservative solute (chloride). We find that estimates of nitrate removal are consistent with previous estimates at the site and they show a good (inverse) correlation with the fraction of streamflow that is younger than 2.5 months. However, this young water fraction cannot be used to predict nitrate removal in the winter–spring period, when no removal is observed regardless of streamflow age. While this leads us to reject our hypothesis during the winter period, it also suggests that water age distributions and their correlation with nitrate removal can possibly reveal distinct sources of stream water at different hydrologic regimes and relevant biogeochemical reactions.     

Hydrological flow paths at the allt A' mharcaidh catchment: An analysis of plot and catchment scale observations

An integrated programme of hydrological monitoring at the 10 km² Allt a' Mharcaidh catchment in north-east Scotland has been based on observations at plot, hillslope and catchment scale. The resonse of the principal soil types has been characterized from a combination of throughflow and three-dimensional tensiometer data at plot scale, and plot sequences have been used to investigate hillslope scale effects. Seep emergence is associated with downslope drainage and local topographic convergence; in parallel preferential pathways generate a highly dynamic throughflow response. Catchment and subcatchment hydrographs mirror the twin dynamic observed at hillslope scale, and a unified hypothesis of response is presented which is consistent with all scales of observations.     

Completing the Neogene geological time scale between 8.5 and 12.5 Ma

A much improved high-resolution magnetostratigraphy is presented for the interval between 8.5 and 12.5 Ma in the deep marine astronomically tuned Monte dei Corvi Beach section in northern Italy. This section contains the formally designated boundary between the Middle and Upper Miocene defined by the Tortonian GSSP. The natural remanent magnetization is carried by an unusual low-temperature component of primary origin. According to interpretations of isothermal remanent magnetization (IRM) acquisition curves, the palaeomagnetic signal is most probably carried by iron sulphides, such as greigite or pyrrhotite. The resultant magnetostratigraphy is reliable and can be straightforwardly calibrated to the Astronomically Tuned Neogene Time Scale (ATNTS2004) showing that the section ranges from Chron C5An.2n up to C4.2r-1. The inferred correlation of the Tortonian GSSP to the older part of Chron C5r.2n is confirmed, thereby guaranteeing global correlation potential. In addition, the Monte dei Corvi Beach section has been re-tuned to the new numerical solution La2004 for the interval between 8.5 and 12.9 Ma. This exercise shows that the characteristic cycle patterns are not consistent with a single solution with fixed values for dynamical ellipticity and tidal dissipation. For the younger part of the section, sedimentary cycles show an excellent fit with intricate details of precession/obliquity interference in the nominal La2004_(1,1) solution. However, this solution produces serious mismatches between the sedimentary cycles and the precession/obliquity interference in the interval older than ∼ 10.2 Ma. For this interval changes in tidal dissipation value (while keeping dynamical ellipticity constant at its present-day value) are necessary to maintain a perfect fit. A detailed visual comparison reveals a very good to excellent fit with La2004_(1,1.2) and, to lesser extent, with La2004_(1,0.5). However, La2004_(1,1.2) is the preferred tuning target for the moment in the interval older than 10.2 Ma. The astronomical (re)tuning on precessional scale to La2004_(1,1) yielded slightly different astronomical ages for each basic cycle and, hence, also for the calcareous plankton events and magnetic reversal boundaries with uncertainties ranging within a few thousand years. The age of the Tortonian GSSP turns into 11.625 Ma instead of the previously published 11.608 Ma. Our ages should replace existing ages in the Astronomically Tuned Neogene Time Scale (ATNTS2004).     

Coupled modeling of rainfall-induced floods and sediment transport at the catchment scale

Rainfall-induced floods may trigger intense sediment transport on erodible catchments, especially on the Loess Plateau in China, which in turn modifies the floods. However, the role of sediment transport in modifying floods has to date remained poorly understood. Concurrently, traditional hydrodynamic models for rainfall-induced floods typically ignore sediment transport, which may lead to inaccurate results for highly erodible catchments. Here, a two-dimensional (2D) coupled shallow water hydro-sediment-morphodynamic (SHSM) model, based on the Finite Volume Method on unstructured meshes and parallel computing, is proposed and applied to simulate rainfall-induced floods in the Zhidan catchment on the Loess Plateau, Shaanxi Province, China. For six historical floods of return periods up to 2 years, the numerical results compare well with observations of discharge hydrographs at the catchment outlet. The computed runoff-sediment yield relation is quantitatively reasonable as compared with other catchments under similar geographical conditions. It is revealed that neglecting sediment transport leads to underestimation of peak discharge of the flood by 14%–45%, whilst its effect on the timing of the peak discharge varies for different flood events. For 18 design floods with return periods of 10–500 years, sediment transport may lead to higher peak discharge by around 9%–15%. The temporal pattern of concentrated rainfall in a short period may lead to a larger exponent value of the power function for the runoff-sediment yield relation. The current finding leads us to propose that incorporating sediment transport in rainfall-induced flood modeling is warranted. The SHSM model is applicable to flood and sediment modeling at the catchment scale in support of risk management and water and soil management.     

A catchment scale evaluation of the SIBERIA and CAESAR landscape evolution models

Landscape evolution models provide a way to determine erosion rates and landscape stability over times scales from tens to thousands of years. The SIBERIA and CAESAR landscape evolution models both have the capability to simulate catchment–wide erosion and deposition over these time scales. They are both cellular, operate over a digital elevation model of the landscape, and represent fluvial and slope processes. However, they were initially developed to solve research questions at different time and space scales and subsequently the perspective, detail and process representation vary considerably between the models. Notably, CAESAR simulates individual events with a greater emphasis on fluvial processes whereas SIBERIA averages erosion rates across annual time scales. This paper describes how both models are applied to Tin Camp Creek, Northern Territory, Australia, where soil erosion rates have been closely monitored over the last 10 years. Results simulating 10 000 years of erosion are similar, yet also pick up subtle differences that indicate the relative strengths and weaknesses of the two models. The results from both the SIBERIA and CAESAR models compare well with independent field data determined for the site over different time scales. Representative hillslope cross‐sections are very similar between the models. Geomorphologically there was little difference between the modelled catchments after 1000 years but significant differences were revealed at longer simulation times. Importantly, both models show that they are sensitive to input parameters and that hydrology and erosion parameter derivation has long‐term implications for sediment transport prediction. Therefore selection of input parameters is critical. This study also provides a good example of how different models may be better suited to different applications or research questions. Copyright © 2010 John Wiley & Sons, Ltd and Commonwealth of Australia     

鄱阳湖换水周期与示踪剂传输时间特征的数值模拟

在复杂湖泊水动力环境作用下,换水周期和传输时间变化直接影响着污染物的迁移和转化.本文运用数值模拟方法,定量研究了季节水情动态下鄱阳湖换水周期和示踪剂传输时间的空间分布.结果表明,不同季节下鄱阳湖换水周期均具有较高的空间异质性,贯穿整个湖区的主河道换水周期约10 d,大多湖湾区的换水周期则长达300多天.尽管不同季节下换水周期空间分布格局几乎相似,但受鄱阳湖水动力场的季节变化影响,夏、秋季的换水周期要明显大于春、冬季.基于换水周期频率分布曲线的统计表明,80%的鄱阳湖区的换水周期约30 d,其余湖区换水周期为几十天至几百天,表明鄱阳湖应该更加确切地描述为一个快速换水和慢速换水同时共存的湖泊系统.鄱阳湖示踪剂传输时间介于4~32 d,夏、秋季的传输时间(11~32 d)约为春、冬季(4~8 d)的4倍,主要与鄱阳湖季节性水情特征及示踪剂的迁移路径有关.本文所获取的换水周期和示踪剂传输时间的时空分布信息可为今后鄱阳湖水质、水环境和生态系统管理和维护等方面提供重要科学参考.     

An analysis of mean transition time between flood and drought in the large river basin

 A soil moisture balance equation over large spatial regions is studied at seasonal and annual time scales for the Arkansas river basin. Interaction and feedback effects between land-surface and atmospheric moisture are studied in the parameterization for this basin. Due to the interaction between the land-surface and atmosphere at large scales, the surface hydrology of large land areas is susceptible to two distinct stable modes in the long-term probability density function: a dry and a wet state. In the soil moisture balance equation, stochastic fluctuations lead to separate preferred statistical stable states with transitions between these stable states induced by environmental fluctuations. On the basis of historical data, the soil moisture balance equation is calibrated for the Arkansas river basin. The transition times between the stable modes in the model are studied based on the stochastic representation of the physical processes and the calibrated model parameters. This study has implications for prediction of the transition times between stable modes or residence times, that is, the time the system spends in a given stable mode, since this would be equivalent to predicting the duration of droughts or wet conditions.     

Hydrometeorological relationships in a glacierized catchment in the Canadian high arctic

Runoff from a small glacierized catchment in the Canadian high Arctic was monitored throughout one melt season. The stream discharge record is one aspect of a larger project involving glacier mass balance, superimposed ice formation and local climate on a glacier in the Sawtooth Range, Ellesmere Island, Northwest Territories, Canada. To better understand the main factors influencing the production of runoff on the glacier during the period of main summer melt, regression analyses were performed relating daily air temperature, shortwave incoming and net radiation, absorptivity and wind speed to daily glacier discharge. Air temperature at the glacier meteorological station on rain-free days is the element with the greatest correlation with runoff (r² = 0.57; n = 34). A multiple regression of discharge with air temperature, shortwave incoming radiation, net radiation hours and wind speed achieved the best fit (r² = 0.84; n = 34). Rain events (> 10mmd^?1) can dominate daily discharge when they occur during the period of ice melt, creating more runoff per unit area than can be produced by melt alone, and significantly reduce the accuracy of runoff predictions.     

Meso‐scale catchment sediment budgets: combining field surveys and modeling in the Dragonja catchment,southwest Slovenia

In this paper, we present a methodology to construct a sediment budget for meso‐scale catchments. We combine extensive field surveys and expert knowledge of the catchment with a sediment delivery model. The meso‐scale Mediterranean drainage basin of the Dragonja (91 km²), southwest Slovenia, was chosen as case study area. During the field surveys, sheet wash was observed on sloping agricultural fields during numerous rainfall events, which was found to be the main source of sediment. With the sediment yield model WATEM/SEDEM the estimated net erosion on the hillslopes 4·1 t ha^–1 y^–1 (91% of inputs). The second source, bank erosion (4·2%; 0·25 t ha^–1 y^–1) was monitored during several years with erosion pins and photogrammetric techniques. The last source, channel incision, was derived from geomorphological mapping and lichenomery and provided 3·8% (0·17 t ha^–1 y^–1) of the sediment input. The river transports its suspended sediment mainly during high‐flow events (sampled with automated water samplers). About 27% (1·2 t ha^–1 y^–1) of the sediment delivered to the channel is deposited on floodplains and low terraces downstream (estimated with geomorphological mapping, coring and cesium‐137 measurements). The sediment transported as bedload disintegrates during transport to the outlet due to the softness of the bedrock material. As a result, the river carries no bedload when it reaches the sea. The results imply a build‐up of sediment in the valleys catchment. However, extreme flood events may flush large amounts of sediment stored in the lower parts of the system. Geomorphological evidence exists in the catchment that such high magnitude, low frequency events have happened in the past. Copyright © 2009 John Wiley & Sons, Ltd.     

Development and application of a catchment scale diffuse nitrate modelling tool

Water pollution from diffuse sources is a problem of increasing concern. Efforts to control diffuse pollution have been confined mainly to agricultural land and forests. Little attention has been paid to sources of diffuse pollution from urban areas. A diffuse nitrate modelling tool (DNMT) has been developed for modelling the fate of nitrate in urban areas. This tool works at the catchment scale and has a modular structure that consists of three components: the hydrological module, the nitrogen cycle module and the nitrate transport module. The hydrological model describes the possible flow pathways. The nitrogen cycle model accounts for the mass balance of nitrate and calculates the amount of nitrate for potential loss. The nitrate transport module simulates the movement of nitrate within and from the soil to the receiving water. This paper demonstrates the development of the tool and its application in the White Cart Water catchment. This implementation of the tool shows that it has a good capability for simulating the fate of nitrate in urban catchments. Copyright © 2005 John Wiley & Sons, Ltd.