Variable power‐law scaling of hillslope Hortonian rainfall–runoff processes

Hydrological studies focused on Hortonian rainfall–run‐off scaling have found that the run‐off depth generally declines with the plot length in power‐law scaling. Both the power‐law proportional coefficient and the scaling exponent show great variability for specific conditions, but why and how they vary remain unclear. In the present study, the scaling of hillslope Hortonian rainfall–run‐off processes is investigated for different rainfall, soil infiltration, and hillslope surface characteristics using the physically based cell‐based rainfall‐infiltration‐run‐off model. The results show that both temporally intermittent and steady rainfalls can result in prominent power‐law scaling at the initial stage of run‐off generation. Then, the magnitude of the power‐law scaling decreases gradually due to the decreasing run‐on effect. The power‐law scaling is most sensitive to the rainfall and soil infiltration parameters. When the ratio of rainfall to infiltration exceeds a critical value, the magnitude of the power‐law scaling tends to decrease notably. For different intermittent rainfall patterns, the power‐law exponent varies in the range of ?1.0 to ?0.113, which shows an approximately logarithmic increasing trend for the proportional coefficient as a function of the run‐off coefficient. The scaling is also sensitive to the surface roughness, soil sealing, slope angle, and hillslope geometry because these factors control the run‐off routing and run‐on infiltration processes. These results provide insights into the variable scaling of the Hortonian rainfall–run‐off process, which are expected to benefit modelling of large‐scale hydrological and ecological processes.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Spatially distributed long‐term hydrologic simulation using a continuous SCS CN method‐based hybrid hydrologic model

A continuous Soil Conservation Service (SCS) curve number (CN) method that considers time‐varied SCS CN values was developed based on the original SCS CN method with a revised soil moisture accounting approach to estimate run‐off depth for long‐term discontinuous storm events. The method was applied to spatially distributed long‐term hydrologic simulation of rainfall‐run‐off flow with an underlying assumption for its spatial variability using a geographic information systems‐based spatially distributed Clark's unit hydrograph method (Distributed‐Clark; hybrid hydrologic model), which is a simple few parameter run‐off routing method for input of spatiotemporally varied run‐off depth, incorporating conditional unit hydrograph adoption for different run‐off precipitation depth‐based direct run‐off flow convolution. Case studies of spatially distributed long‐term (total of 6 years) hydrologic simulation for four river basins using daily NEXRAD quantitative precipitation estimations demonstrate overall performances of Nash–Sutcliffe efficiency (E_NS) 0.62, coefficient of determination (R²) 0.64, and percent bias 0.33% in direct run‐off and E_NS 0.71, R² 0.72, and percent bias 0.15% in total streamflow for model result comparison against observed streamflow. These results show better fit (improvement in E_NS of 42.0% and R² of 33.3% for total streamflow) than the same model using spatially averaged gauged rainfall. Incorporation of logic for conditional initial abstraction in a continuous SCS CN method, which can accommodate initial run‐off loss amounts based on previous rainfall, slightly enhances model simulation performance; both E_NS and R² increased by 1.4% for total streamflow in a 4‐year calibration period. A continuous SCS CN method‐based hybrid hydrologic model presented in this study is, therefore, potentially significant to improved implementation of long‐term hydrologic applications for spatially distributed rainfall‐run‐off generation and routing, as a relatively simple hydrologic modelling approach for the use of more reliable gridded types of quantitative precipitation estimations.     

Simulating pan‐Arctic runoff with a macro‐scale terrestrial water balance model

A terrestrial hydrological model, developed to simulate the high‐latitude water cycle, is described, along with comparisons with observed data across the pan‐Arctic drainage basin. Gridded fields of plant rooting depth, soil characteristics (texture, organic content), vegetation, and daily time series of precipitation and air temperature provide the primary inputs used to derive simulated runoff at a grid resolution of 25 km across the pan‐Arctic. The pan‐Arctic water balance model (P/WBM) includes a simple scheme for simulating daily changes in soil frozen and liquid water amounts, with the thaw–freeze model (TFM) driven by air temperature, modelled soil moisture content, and physiographic data. Climate time series (precipitation and air temperature) are from the National Centers for Environmental Prediction (NCEP) reanalysis project for the period 1980–2001. P/WBM‐generated maximum summer active‐layer thickness estimates differ from a set of observed data by an average of 12 cm at 27 sites in Alaska, with many of the differences within the variability (1σ) seen in field samples. Simulated long‐term annual runoffs are in the range 100 to 400 mm year^?1. The highest runoffs are found across northeastern Canada, southern Alaska, and Norway, and lower estimates are noted along the highest latitudes of the terrestrial Arctic in North America and Asia. Good agreement exists between simulated and observed long‐term seasonal (winter, spring, summer–fall) runoff to the ten Arctic sea basins (r = 0·84). Model water budgets are most sensitive to changes in precipitation and air temperature, whereas less affect is noted when other model parameters are altered. Increasing daily precipitation by 25% amplifies annual runoff by 50 to 80% for the largest Arctic drainage basins. Ignoring soil ice by eliminating the TFM sub‐model leads to runoffs that are 7 to 27% lower than the control run. The results of these model sensitivity experiments, along with other uncertainties in both observed validation data and model inputs, emphasize the need to develop improved spatial data sets of key geophysical quantities (particularly climate time series) to estimate terrestrial Arctic hydrological budgets better. Copyright © 2003 John Wiley & Sons, Ltd.     

First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach

Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 John Wiley & Sons, Ltd.     

A mass‐balance approach to estimate in‐stream processes in a large river

A mass‐balance approach was used to estimate in‐stream processes related to inorganic nitrogen species (NH₄⁺, NO₂^? and NO₃^?) in a large river characterized by highly variable hydrological conditions, the Garonne River (south‐west France). Studies were conducted in two consecutive reaches of 30 km located downstream of the Toulouse agglomeration (population 760 000, seventh order), impacted by modification of discharge regime and high nitrogen concentrations. The mass‐balance was calculated by two methods: the first is based on a variable residence time (VRT) simulated by a one‐dimensional (1‐D) hydraulic model; the second is a based on a calculation using constant residence time (CRT) evaluated according to hydrographic peaks. In the context of the study, removal of dissolved inorganic nitrogen (DIN) for a reach of 30 km is underestimated by 11% with the CRT method. In sub‐reaches, the discrepancy between the two methods led to a 50% overestimation of DIN removal in the upper reach (13 km) and a 43% underestimation in the lower reach (17 km) using the CRT method. The study highlights the importance of residence time determination when using modelling approaches in the assessment of whole stream processes in short‐duration mass‐balance for a large river under variable hydrological conditions. Copyright © 2007 John Wiley & Sons, Ltd.     

A transport‐distance approach to scaling erosion rates: 3. Evaluating scaling characteristics of Mahleran

In the two previous papers of this series, we demonstrated how a novel approach to erosion modelling (Mahleran – Model for Assessing Hillslope‐Landscape Erosion, Runoff And Nutrients) provided distinct advantages in terms of process representation and explicit scaling characteristics when compared with existing models. A first evaluation furthermore demonstrated the ability of the model to reproduce spatial and temporal patterns of erosion and their particle‐size characteristics on a large rainfall‐simulation plot. In this paper, we carry out a more detailed evaluation of the model using monitored erosion events on plots of different size. The evaluation uses four plots of 21·01, 115·94, 56·84 and 302·19 m², with lengths of 4·12, 14·48, 18·95 and 27·78 m, respectively, on similar soils to the rainfall‐simulation plot, for which runoff and erosion were monitored under natural rainfall. Although the model produces the correct ranking of the magnitude of erosion events, it performs less well in reproducing the absolute values and particle‐size distributions of the eroded sediment. The implications of these results are evaluated in terms of requirements for process understanding and data for parameterization of improved soil‐erosion models. We suggest that there are major weaknesses in the current understanding and data underpinning existing models. Consequently, a more holistic re‐evaluation is required that produces functional relationships for different processes that are mutually consistent, and that have appropriate parameterization data to support their use in a wide range of environmental conditions. Copyright © 2008 John Wiley & Sons, Ltd.     

Modelling soil solute release and transport in run‐off on a loessial slope with and without surface stones

Stone covers on loessial slopes can increase the time of infiltration by slowing the velocity of the overland flow, which reduces the transport of solutes, but few mechanistic models have been tested under water‐scouring conditions. We carried out field experiments to test a previously proposed, physically based model of water and solute transport. The area of soil infiltration was calculated from the uncovered surface area, and Richards' equation and the kinematic wave equation were used to describe water infiltration and flow along slopes with stone covers. The transport of chemicals into the run‐off from the surface soil, presumably by diffusion, and their movement in the soil profile could be described by the convection–diffusion equations of the model. The simulated and measured data correlated well. The stones on the soil surface reduced the area available for infiltration but increased the Manning coefficient, eventually leading to increased water infiltration and decreased solute loss with run‐off. Our results indicated that the traditional model of water movement and solute migration could be used to simulate water transport and solute migration for stone‐covered soil on loessial slopes.     

Understanding changes in annual runoff following land use changes: a systematic data‐based approach

A simple modelling framework for assessing the response of ungauged catchments to land use change in South‐Western Australia is presented. The framework uses knowledge of transpiration losses from native vegetation and pasture and then partitions the ‘excess’ water (resulting from reduced transpiration after land use change) between runoff and deep storage. The simple partitioning is achieved by using soft information (satellite imagery, previous mapping and field assessment) to delimit the spread of the permanently saturated area close to the stream. Runoff is then assumed to increase in proportion to the saturated area, with the residual difference going to deep storage. The model parameters to describe the annual water yield are obtained a priori and no calibration to streamflow is required. We tested the model using gauged records over 25 years from paired catchment experiments in South‐Western Australia. Very good estimates of runoff were obtained from high rainfall (>1100 mm yr⁻¹) catchments (R² > 0·9) and for low rainfall (<900 mm yr⁻¹) catchments after clearing (R² = 0·96) but results were poorer (R² = 0·55) for an uncleared low rainfall catchment, although overall balances were reasonable. In the drier uncleared catchments, the within‐year distributions of rainfall may exert a substantial influence on runoff response that is not completely captured by the presented model. Copyright © 2005 John Wiley & Sons, Ltd.     

Estimating source regions for snowmelt runoff in a Rocky Mountain basin: tests of a data‐based conceptual modeling approach

In many mountain basins, river discharge measurements are located far away from runoff source areas. This study tests whether a basic snowmelt runoff conceptual model can be used to estimate relative contributions of different elevation zones to basin‐scale discharge in the Cache la Poudre, a snowmelt‐dominated Rocky Mountain river. Model tests evaluate scenarios that vary model configuration, input variables, and parameter values to determine how these factors affect discharge simulation and the distribution of runoff generation with elevation. Results show that the model simulates basin discharge well (NSCE and R >0.90) when input precipitation and temperature are distributed with different lapse rates, with a rain‐snow threshold parameter between 0 and 3.3 °C, and with a melt rate parameter between 2 and 4 mm °C^?1 d^?1 because these variables and parameters can have compensating interactions with each other and with the runoff coefficient parameter. Only the hydrograph recession parameter can be uniquely defined with this model structure. These non‐unique model scenarios with different configurations, input variables, and parameter values all indicate that the majority of basin discharge comes from elevations above 2900 m, or less than 25% of the basin total area, with a steep increase in runoff generation above 2600 m. However, the simulations produce unrealistically low runoff ratios for elevations above 3000 m, highlighting the need for additional measurements of snow and discharge at under‐sampled elevations to evaluate the accuracy of simulated snow and runoff patterns. Copyright © 2013 John Wiley & Sons, Ltd.     

Run‐off formation in a humid,temperate headwater catchment using a combined hydrological,hydrochemical and isotopic approach (Jizera Mountains,Czech Republic)

The formation of baseflow and stormflow was examined in the 1.18 km² part of the headwater catchment Uhlí?ská, Jizera Mountains, Czech Republic, over the period 2007–2011, by means of run‐off data and environmental tracers ¹⁸O and SiO₂. The baseflow, computed using the digital filter approach BFLOW, contributes 67% to total streamflow and has a mean residence time of 12.3 months. It is formed by groundwater discharge from the valley deluviofluvial granitic sediments, in combination with soil water in weathered layers on hillslopes during rainfall and snowmelt periods. The prevailing source of the groundwater is the infiltration of snowmelt water. Analysis of 20 run‐off events and their hysteretic patterns demonstrated that the stormflow water has a residence time of about 4 months and is generated by preferential flow on hillslopes combined by soil matrix drainage. Because of slower flow in the soil matrix, the enrichment of pore water in SiO₂ is more pronounced. The stormflow and snowmelt water flowing via preferential pathways of upslope minerals soils pushes the pre‐event groundwater through the pathways in wetlands to the stream, and the wetland can be therefore considered as groundwater supplied. This mechanism has been found to be typical for the groundwater‐supplied headwater catchments of the Jizera Mountains and can be also assumed in other mountainous headwaters of the granitic massif in Central Europe. The main methodological contribution of this study are the residence time calculations stratified by baseflow and event flow, identifying run‐off components of different travel times to streams and linking them with geochemical run‐off sources. This achievement was possible because of a comprehensive dataset on hydrology, stable isotopes and silica hydrochemistry in all relevant run‐off generation components. This concept indicates that a possible long‐term change in snowmelt may affect the run‐off regime of headwater catchments to climate or land‐use changes. Copyright © 2013 John Wiley & Sons, Ltd.     

Removal of water‐surface reflection effects with a temporal minimum filter for UAV‐based shallow‐water photogrammetry

The recent development of structure‐from‐motion (SfM) and multi‐view stereo (MVS) photogrammetry techniques has enabled semi‐automatic high‐resolution bathymetry using aerial images taken by consumer‐grade digital cameras mounted on unmanned aerial vehicles (UAVs). However, the applicability of these techniques is sometimes limited by sun and sky reflections at the water surface, which render the point‐cloud density and accuracy insufficient. In this research, we present a new imaging technique to suppress the effect of these water‐surface reflections. In this technique, we order a drone to take a short video instead of a still picture at each waypoint. We then apply a temporal minimum filter to the video. This filter extracts the smallest RGB values in all the video frames for each pixel, and composes an image with greatly reduced reflection effects. To assess the performance of this technique, we applied it at three small shallow‐water sites. Specifically, we evaluated the effect of the technique on the point cloud density and the accuracy and precision of the photogrammetry. The results showed that the proposed technique achieved a far denser point cloud than the case in which a randomly chosen frame was used for each waypoint, and also showed better overall accuracy and precision in estimating water‐bottom elevation. The effectiveness of this new technique should depend on the surface wave state and sky radiance distribution, and this dependence, as well as the applicability to large areas, should be investigated in future research. Copyright © 2018 John Wiley & Sons, Ltd.     

Identification of factors influencing hydrologic model performance using a top-down approach in a large number of U.S. catchments

Investigating the performance that can be achieved with different hydrological models across catchments with varying characteristics is a requirement for identifying an adequate model for any catchment, gauged or ungauged, just based on information about its climate and catchment properties. As parameter uncertainty increases with the number of model parameters, it is important not only to identify a model achieving good results but also to aim at the simplest model still able to provide acceptable results. The main objective of this study is to identify the climate and catchment properties determining the minimal required complexity of a hydrological model. As previous studies indicate that the required model complexity varies with the temporal scale, the study considers the performance at the daily, monthly, and annual timescales. In agreement with previous studies, the results show that catchments located in arid areas tend to be more difficult to model. They therefore require more complex models for achieving an acceptable performance. For determining which other factors influence model performance, an analysis was carried out for four catchment groups (snowy, arid, and eastern and western catchments). The results show that the baseflow and aridity indices are the most consistent predictors of model performance across catchment groups and timescales. Both properties are negatively correlated with model performance. Other relevant predictors are the fraction of snow in the annual precipitation (negative correlation with model performance), soil depth (negative correlation with model performance), and some other soil properties. It was observed that the sign of the correlation between the catchment characteristics and model performance varies between clusters in some cases, stressing the difficulties encountered in large sample analyses. Regarding the impact of the timescale, the study confirmed previous results indicating that more complex models are needed for shorter timescales.     

Season‐based rainfall–runoff modelling using the probability‐distributed model (PDM) for large basins in southeastern Brazil

Southeastern Brazil is characterized by seasonal rainfall variability. This can have a great social, economic, and environmental impact due to both excessive and deficient water availability. During 2014 and 2015, the region experienced one of the most severe droughts since 1960. The resulting water crisis has seriously affected water supply to the metropolitan region of São Paulo and hydroelectric power generation throughout the entire country. This research considered the upstream basins of the southeastern Brazilian reservoirs Cantareira (2,279 km²; water supply) and Emborcação (29,076 km²), Três Marias (51,576 km²), Furnas (52,197 km²), and Mascarenhas (71,649 km²; hydropower) for hydrological modelling. It made the first attempt at configuring a season‐based probability‐distributed model (PDM‐CEMADEN) for simulating different hydrological processes during wet and dry seasons. The model successfully reproduced the intra‐annual and interannual variability of the upstream inflows during 1985–2015. The performance of the model was very satisfactory not only during the wet, dry, and transitional seasons separately but also during the whole period. The best performance was obtained for the upstream basin of Furnas, as it had the highest quality daily precipitation and potential evapotranspiration data. The Nash–Sutcliffe efficiency and logarithmic Nash–Sutcliffe efficiency were 0.92 and 0.93 for the calibration period 1984–2001, 0.87 and 0.88 for the validation period 2001–2010, and 0.93 and 0.90 for the validation period 2010–2015, respectively. Results indicated that during the wet season, the upstream basins have a larger capacity and variation of soil water storage, a larger soil water conductivity, and quicker surface water flow than during the dry season. The added complexity of configuring a season‐based PDM‐CEMADEN relative to the traditional model is well justified by its capacity to better reproduce initial conditions for hydrological forecasting and prediction. The PDM‐CEMADEN is a simple, efficient, and easy‐to‐use model, and it will facilitate early decision making and implement adaptation measures relating to disaster prevention for reservoirs with large‐sized upstream basins.     

用稳定同位素方法估算大型浅水湖泊蒸发量——以太湖为例

湖泊蒸发量的准确估算对于水文学、气象学和湖泊学等研究有重要的意义.基于2013-2015年太湖水量收支资料、气象观测数据和稳定同位素观测资料,采用稳定同位素质量守恒模型、水量平衡法和Priestley-Taylor模型估算太湖蒸发量,分析太湖蒸发量的季节变化和年际变化特征,并以Priestley-Taylor模型结果为参考值,评价水量平衡法和同位素质量守恒方程的计算精度.结果表明:5-9月太湖蒸发量较高,冬季最低.2013-2015年太湖年总蒸发量分别为1069、894和935 mm,蒸发量的年际变化受到天气条件的影响.2013年12月2014年11月期间,用Priestley-Taylor模型计算的湖泊蒸发量为885 mm;同位素质量守恒模型的估算结果较一致,为893 mm;而水量平衡方程的估算结果明显偏高,为1247 mm.     

Characterization of total suspended solid dynamics in a large shallow lake using long‐term daily satellite images

Large, shallow‐water lakes located on floodplains play an important role in creating highly productive ecosystems and are prone to high concentrations of suspended solids due to sediment resuspension. In this study, the aim was to determine the dominant processes governing the total suspended solid (TSS) concentration at the water surface in Tonle Sap Lake, Cambodia, which is a large, shallow‐water lake. Satellite remotely sensed daily reflectance data from 2003 to 2017 were used. Seasonal changes in TSS concentration indicated that bottom sediment resuspension during dry seasons was mostly caused by wind and the TSS concentration was closely correlated with the water depth of the lake. The TSS concentration during flood periods was controlled by both wind and inflow currents from the Tonle Sap River. Additionally, we confirmed that surface/subsurface flow with a low TSS concentration from forests on the floodplain lowered the TSS concentration year round, except during August and September. This fact implied that the floodplain forest area decrease may increase the lake TSS concentration. An analysis of the long‐term changes in TSS indicated that a decrease in the water level during flood periods resulted in the high TSS concentrations observed during the subsequent dry periods. Therefore, climate change and water resource development, which are likely to cause water level reductions in the Mekong River during flood periods, may increase the TSS concentration in Tonle Sap Lake, particularly during the dry season.     

Application of the WEPP model to determine sources of run‐off and sediment in a forested watershed

This study investigates critical run‐off and sediment production sources in a forested Kasilian watershed located in northern Iran. The Water Erosion Prediction Project (WEPP) watershed model was set up to simulate the run‐off and sediment yields. WEPP was calibrated and validated against measured rainfall–run‐off–sediment data. Results showed that simulated run‐off and sediment yields of the watershed were in agreement with the measured data for the calibration and validation periods. While low and medium values of run‐off and sediment yields were adequately simulated by the WEPP model, high run‐off and sediment yield values were underestimated. Performance of the model was evaluated as very good and satisfactory during the calibration and validation stages, respectively. Total soil erosion and sediment load of the study watershed during the study period were determined to be 10 108 t yr^?1 and 8735 t yr^?1, respectively. The northern areas of the watershed with dry farming were identified as the critical erosion prone zones. To prioritize the subwatersheds based on their contribution to the run‐off and sediment production at the watershed's main outlet, unit response approach (URA) was applied. In this regard, subwatersheds close to the main outlet were found to have the highest contribution to sediment yield of the whole watershed. Results indicated that depending on the objective of land and water conservation practices, particularly, for controlling sediment yield at the main outlet, critical areas for implementing the best management practices may be identified through conjunctive application of WEPP and URA. Copyright © 2014 John Wiley & Sons, Ltd.     

A tale of two isotopes: differences in hydrograph separation for a runoff event when using δD versus δ18O

It is often assumed that stable water isotopes (δD and δ¹⁸O) provide redundant information for a given sample of water. In this note we illustrate that the choice of isotope used may influence the resultant hydrograph separation. This is especially true in light of the spatial and temporal variability in the isotopic composition of rainfall water at the catchment scale. We present several possible hydrograph separations based on both δD and δ¹⁸O observed in rainfall for a single runoff event occurring in the southwest USA. This study demonstrates the potential of using both stable water isotopes by showing that δD and δ¹⁸O may provide unique information for catchment hydrologists. We also report on the utility of new technology capable of simultaneous measurements of both δD and δ¹⁸O using off‐axis integrated cavity output spectroscopy (OA‐ICOS) methods. This may be of interest to catchment hydrologists seeking to incorporate this type of equipment into their laboratory. Copyright © 2009 John Wiley & Sons, Ltd.     

Variable water‐saturated areas and streamflow generation in the small Ringelbach catchment (Vosges Mountains,France): the master recession curve as an equilibrium curve for interactions between atmosphere,surface and ground waters

In the small Ringelbach research catchment, where studies on the water cycle components in a granitic mountainous environment have been conducted since 1976, the water‐saturated areas that are hydraulically connected to the outlet play a major role in the streamflow generation, as it is here that complex interactions between atmosphere, surface and ground waters take place. During baseflow recession periods, which may last several months between two groundwater recharge events, the atmospheric inputs of water and energy on these contributing areas only explain the streamflow fluctuations observed around the master recession curve, which defines the groundwater contribution: fluctuating above it in the case of precipitation input on these areas, below it in the case of evaporation output from these areas. Streamflow may therefore largely deviate from the master recession curve in the case of long, hot, dry spells. Detailed mapping has shown that their variable extent is well related to baseflow by a loglinear curve. On the other hand, a synthetic master recession curve, well fitted by a second‐order hyperbolic function, has been obtained from numerous pure recession periods. Both based on these two curves, a simple procedure and a simple model have been used to (i) validate the hypothesis that the connected saturated areas are the only permanent variable contributing areas and (ii) simulate the daily streamflow volumes over long baseflow recession periods by a water balance of the aquifer below these areas only. The storm runoff ratio for small to moderate rainfall events is indeed corresponding to the catchment saturated fraction at that time. The volume of daily streamflow oscillations is indeed corresponding to the evaporation at the potential rate from the saturated areas only. In both cases, streamflow naturally tends towards the master recession curve after the end of any atmospheric perturbation. Introducing these findings into TOPMODEL led to significantly improved simulation results during baseflow recession periods. The master recession curve may therefore be considered as a dynamic equilibrium curve. Together with the relationship between saturated extent and baseflow, it provides the main characteristics necessary to understand and model the interactions at this complex interface and the resulting daily streamflow variations during baseflow recession periods in this type of catchment. Copyright © 2016 John Wiley & Sons, Ltd.