ALBERT

Description: The objective of this study is to incorporate a time-dependent SCS CN method (SMA_CN) in Soil and Water Assessment Tool (SWAT) and compare its performance with the existing CN method in SWAT by simulating the hydrology of two agricultural watersheds in Indiana, United States. Results show that fusion of the SMA_CN method causes decrease in runoff volume and increase in profile soil moisture content, associated with larger groundwater contribution to the streamflow. In addition, the higher amount of moisture in the soil profile slightly elevates the actual evapotranspiration. The SMA-based SWAT configuration consistently produces improved goodness of fit scores and less uncertain outputs with respect to streamflow during both calibration and validation. The SMA_CN method exhibits better match with the observed data for all flow regimes, thereby addressing issues related to peak and low flow predictions by SWAT in many past studies. Comparison of the calibrated model outputs with field-scale soil moisture observations reveal that the SMA overhauling enables SWAT to represent soil moisture condition more accurately, with better response to the incident rainfall dynamics. While the results from the modification of the SCS method in SWAT are promising, more studies including watersheds with various physical and climatic settings are needed to validate the proposed approach. This article is protected by copyright. All rights reserved.

Description: We adapted Newton's Law of Cooling to model downstream water temperature change in response to stream-adjacent forest harvest on small and medium streams (average 327 ha in size) throughout the Oregon Coast Range, USA. The model requires measured stream gradient, width, depth and upstream control reach temperatures as inputs and contains two free parameters which were determined by fitting the model to measured stream temperature data. This model reproduces the measured downstream temperature responses to within 0.4 C ° for 15 of the 16 streams studied and provides insight into the physical sources of site-to-site variation among those responses. We also use the model to examine how the pre-to-post harvest change in daily maximum stream temperature depends on distance from the harvest reach. The model suggests that the pre-to-post harvest temperature change approximately 300  m downstream of the harvest will range from roughly 82% to less than 1% of that temperature change which occurred within the harvest reach, depending primarily on the downstream width, depth, and gradient. Using study-averaged values for these channel characteristics the model suggests that for a stream representative of those in the study, the temperature change approximately 300  m downstream of the harvest will be 56% of the temperature change which occurred within the harvest reach. This adapted Newton's Law of Cooling procedure represents a highly practical means for predicting stream temperature behavior downstream of timber harvests relative to conventional heat budget approaches, and is informative of the dominant processes affecting stream temperature. This article is protected by copyright. All rights reserved.

Description: Digital elevation models (DEMs) that are used in hydrological applications must be processed to remove sinks, mainly topographic depressions. Flow enforcement techniques include filling methods, which raise elevations within depressions, breaching, which carves channels through blockages, and hybrid methods. Despite previous research demonstrating the large impact to DEMs and subsequent analyses of depression filling, it is common practice apply this technique to flow enforcement. This is partly due to the greater efficiency of depression filling tools compared to breaching counterparts, which often limits breaching to applications of small- to moderate-sized DEMs. A new hybrid flow enforcement algorithm is presented in this study. The method can be run in complete breaching, selective breaching (either breached or filled), or constrained breaching (partial breaching) modes, allowing for greater flexibility in how practitioners enforce continuous flow paths. Algorithm performance was tested with DEMs of varying topography, spatial extents, and resolution. The sites included three moderate sized DEMs (52,000,000 to 190,000,000 cells) and three massive DEMs of the Iberian Peninsula, and the Amazon and Nile River basins, the largest containing nearly one billion cells. In complete breaching mode, the new algorithm required 87% of the time needed by a filling method to process the test DEMs, while the selective breaching and constrained breaching modes, operating with maximum breach depth constraints, increased run times by 8% and 27% respectively. Therefore, the new algorithm offers comparable performance to filling and the ability to process massive topographic data sets, while giving practitioners greater flexibility and lowering DEM impact. This article is protected by copyright. All rights reserved.

Description: Mountain snowpacks provide most of the annual discharge of western U.S. rivers, but the future of water resources in the western U.S. is tenuous, as climatic changes have resulted in earlier spring melts that have exacerbated summer droughts. Compounding changes to the physical environment are biotic disturbances including the mountain pine beetle (MPB), which has decimated millions of acres of western North American forests. At the watershed scale, MPB disturbance increases the peak hydrograph, and at the stand scale the ‘gray’ phase of MPB canopy disturbance decreases canopy snow interception, increases snow albedo, increases net shortwave radiation and decreases net longwave radiation versus the ‘red’ phase. Fewer studies have been conducted on the red phase of MPB disturbance, and in the mixed coniferous stands that may follow MPB-damaged forests. We measured the energy balance of four snowpacks representing different stages of MPB damage, management, and recovery: a lodgepole pine stand, a MPB-infested stand in the red phase, a mixed coniferous stand (representing one successional trajectory), and a clearcut (representing reactive management) in the Tenderfoot Creek Experimental Forest in Montana, USA. Net longwave radiation was lower in the MPB-infested stand despite higher basal area and plant area index of the other forests, suggesting that the dessicated needles serve as a less effective thermal buffer against longwave radiative losses. Eddy covariance observations of sensible and latent heat flux indicate that they are of similar but opposite magnitude, on the order of 20 MJ m −2 during the melt period. Further analyses reveal that net turbulent energy fluxes were near zero due to the temperature and atmospheric vapor pressure encountered during the melt period. Future research should place snow science in the context of forest succession and management, and address important uncertainties regarding the timing and magnitude of needlefall events. This article is protected by copyright. All rights reserved.

Description: As a result of climate change/variation and its aggravation by human activities over the past several decades, the hydrological conditions in the middle Yellow River in China have dramatically changed, which has led to a sharp decrease of streamflow and the drying up of certain tributaries. This paper simulated and analysed the impact of sediment-trapping dams (STDs, a type of large-sized check dam used to prevent sediment from entering the Yellow River main stem) on hydrological processes, and the study area was located in the 3,246 km 2 Huangfuchuan (HFC) River basin. Changes in the hydrological processes were analysed, and periods of natural and disturbed states were defined. Subsequently, the number and distribution of the STDs were determined based on data collected from statistical reports and identified from remote sensing images, and the topological relationships between the STDs and high-resolution river reaches were established. A hydrological model, the Digital Yellow River Integrated Model, was used to simulate the STD impact on the hydrological processes, and the maximum STD impact was evaluated through a comparison between the simulation results with and without the STDs, which revealed that the interception effect of the STDs contributed to the decrease of the streamflow by approximately 39%. This paper also analysed the relationship between the spatial distribution of the STDs and rainfall in the HFC River basin and revealed that future soil and water conservation measures should focus on areas with a higher average annual rainfall and higher number of rainstorm hours. This article is protected by copyright. All rights reserved.

Description: Distributed, continuous hydrologic models promote better understanding of hydrology and enable integrated hydrologic analyses by providing a more detailed picture of water transport processes across the varying landscape. However, such models are not widely used in routine modeling practices, due in part to the extensive data input requirements, computational demands, and complexity of routing algorithms. We developed a two-dimensional continuous hydrologic model, HYSTAR, using a time-area method within a grid-based spatial data model with the goal of providing an alternative way to simulate spatiotemporally varied watershed-scale hydrologic processes. The model calculates the direct runoff hydrograph by coupling a time-area routing scheme with a dynamic rainfall excess sub-model implemented here using a modified curve number method with an hourly time step, explicitly considering downstream ‘reinfiltration’ of routed surface runoff. Soil moisture content is determined at each time interval based on a water balance equation, and overland and channel runoff is routed on time-area maps, representing spatial variation in hydraulic characteristics for each time interval in a storm event. Simulating runoff hydrographs does not depend on unit hydrograph theory or on solution of the Saint Venant equation, yet retains the simplicity of a unit hydrograph approach and the capability of explicitly simulating two-dimensional flow routing. The model provided acceptable performance in predicting daily and monthly runoff for a 6-year period for a watershed in Virginia (USA) using readily available geographic information about the watershed landscape. Spatial and temporal variability in simulated effective runoff depth and time area maps dynamically show the areas of the watershed contributing to the direct runoff hydrograph at the outlet over time, consistent with the variable source area overland flow generation mechanism. The model offers a way to simulate watershed processes and runoff hydrographs using the time-area method, providing a simple, efficient, and sound framework that explicitly represents mechanisms of spatially and temporally varied hydrologic processes. This article is protected by copyright. All rights reserved.

Description: Given the importance of groundwater temperature to the biogeochemical health of aquatic ecosystems, a floodplain study was implemented to improve understanding of rural land use impacts on shallow groundwater (SGW) temperature. Study sites included a historic agricultural field (Ag) and bottomland hardwood forest (BHF), each with nine piezometers in an 80 × 80 m grid. Piezometers were equipped with pressure transducers to monitor SGW temperature and level at 30 minute intervals during the 2011, 2012, 2013, and 2014 water years. The study is one of the first to utilize long-term, continuous, automated, in situ monitoring to investigate rural land use impacts on shallow groundwater temperatures. Average SGW temperature during the study period was 11.1 and 11.2 °C at the Ag and BHF sites, respectively. However, temperature range at the Ag site was 72% greater than at the BHF site. Results indicate a greater responsiveness to seasonal climate fluctuations in Ag site SGW temperature related to absence of forest canopy. Patterns of intra-site groundwater temperature differences at both study sites illustrate the influence of stream-aquifer thermal conduction and occasional baseflow reversals. Considering similar surface soil temperature amplitudes and low average groundwater flow values at both sites, results suggest that contrasting rates of plant water use, groundwater recharge, and subsurface hydraulic conductivity are likely mechanistic causes for the observed SGW temperature differences. Results highlight the long-term impact of forest removal on subsurface hydrology and groundwater temperature regime. This article is protected by copyright. All rights reserved.

Description: Growing demand on groundwater resources and the semi-arid climate in the North China Plain (NCP) highlight the need for improved understanding of connections between regional climate change and groundwater recharge. Hydrologic time series of precipitation and groundwater levels were analyzed in three representative geographical zones throughout the NCP for the period of 1960-2008 using trend analysis and spectral analysis methods. A significant change point around 1975 is followed by a long term decline trend in precipitation time series, which coincides with the Pacific Decadal Oscillation (PDO) positive phase. However, the magnitudes of groundwater levels variability due to heavily pumping overwhelm the low-frequency signal of groundwater levels. Nonlinear trends that related to long-term climatic variability and anthropogenic activities are removed by using the Singular Spectrum Analysis (SSA) method. Spectral analyses of the detrended residuals demonstrate significant short-term oscillations at the frequencies of 2–7 years, which have strong correlations with the El Niño-Southern Oscillation (ENSO) modes. This study contributes to improved understanding of dynamic relationship between groundwater and climate variability modes in the NCP, and demonstrates the importance of reliable detrending methods for groundwater levels that are affected greatly by pumping. This article is protected by copyright. All rights reserved.

Description: On September 3, 1998, a glacial lake outburst flood (GLOF) that originated from Tam Pokhari occurred in the Hinku valley of the eastern Nepal Himalaya. This study analyzes the lake's geomorphic and hydrologic conditions prior to the outburst, and evaluates the conditions that could contribute to a future flood through photogrammetric techniques. We processed high-resolution Corona KH-4A (2.7 m) and ALOS PRISM (2.5 m) stereo-images taken before and after the GLOF event, and produced detailed topographic maps (2-m contour interval) and DEMs (5 m × 5 m). We (re-) constructed lake water surfaces before (4410 ± 5 m) and after (4356 ± 5 m) the outburst, and reliably estimated the lake water surface lowering (54 ± 5 m) and the water volume released (19.5 ± 2.2 × 10 6  m 3 ) from the lake, showing good agreement with the results obtained from ground-based measurements. The most relevant conditions that may have influenced the catastrophic drainage of Tam Pokhari in 1998 include the presence of: i) a narrow (75 ± 6 m), steep (up to 50°) and high (120 ± 5 m) moraine dam; ii) high lake level (8 ± 5 m of freeboard); and iii) a steep overhanging glacier (〉40°). The lake outburst substantially altered the immediate area, creating a low and wide (〉500 m) outwash plain below the lake, a wide lake outlet channel (~50 m) and a gentle channel slope (~3–5°). Our new data suggest that the likelihood of a future lake outburst is low. Our results demonstrate that the datasets produced by photogrammetric techniques provide an excellent representation of micro-landform features on moraine dams, lake water surfaces and the changes in both over time, thereby allowing highly accurate pre- and post-GLOF (volumetric) change analysis of glacial lakes. Furthermore, it enables precise measurement of several predictive variables of GLOFs that can be useful for identifying potentially dangerous glacial lakes or prioritizing them for detailed field investigations. This article is protected by copyright. All rights reserved.

Description: Although the importance to account for microrelief in the calculation of specific yields for shallow groundwater systems is well recognized, the microrelief influence is often treated very simplified, which can cause considerable errors. We provide a general one-dimensional expression that correctly represents the effect of a microrelief on the total specific yield that is composed of the soil and surface specific yield. The one-dimensional expression can be applied for different soil hydraulic parameterizations and soil surface elevation frequency distributions. Applying different van Genuchten parameters and a simple linear microrelief model, we demonstrate that the specific yield is influenced by the microrelief not only when surface storage directly contributes to specific yield by (partial) inundation but also when water levels are lower than the minimum surface elevation. Compared to a simplified representation of the soil specific yield, in which a mean soil surface is assumed for the calculation of soil specific yield, the correct representation can lead to lower as well as higher soil specific yields depending on the specific interaction of the soil water retention characteristics and the microrelief. The new equation can be used to obtain more accurate evapotranspiration estimates from water level fluctuations and to account for the effect of microtopographic subgrid variability on simulated water levels of spatially-distributed hydrological models. This article is protected by copyright. All rights reserved.

Description: Globally, Dissolved Inorganic Carbon (DIC) accounts for more than half the annual flux of carbon exported from terrestrial ecosystems via rivers. Here we assess the relative influences of biogeochemical and hydrological processes on DIC fluxes exported from a tropical river catchment characterized by distinct land cover, climate and geology transition from the wet tropical mountains to the low lying savanna plains. Processes controlling changes in river DIC were investigated using dissolved organic carbon (DOC), particulate organic carbon (POC) and DIC concentrations and stable isotope ratios of DIC (δ 13 C DIC ) at two time scales; seasonal and diel. The recently developed Isotopic Continuous Dissolved Inorganic Carbon Analyser (ISO-CADICA) was used to measure diel DIC concentration and δ 13 C DIC changes at a 15 minute temporal resolution. Results highlight the predominance of biologically mediated processes (photosynthesis and respiration) controlling diel changes in DIC. These resulted in DIC concentrations varying between 3.55-3.82 mg/L, and δ 13 C DIC values ranging from -19.7 ± 0.31 to -17.1 ± 0.08 ‰. In contrast, at the seasonal scale we observe wet season DIC variations predominantly from mixing processes, and dry season DIC variations due to both mixing processes and biological processes. The observed wet season increases in DIC concentrations (by 6.81 mg/L) and δ 13 C DIC values of river water (by 5.4 ‰) largely result from proportional increases in subsurface inflows from the savanna plains (C 4 vegetation) region relative to inflows from the rainforest (C 3 vegetation) highlands. The high DIC river load during the wet season results in the transfer of 97% of the annual river carbon load. Therefore, in this gaining river there are significant seasonal variations in both the hydrological and carbon cycles, and there is evidence of substantial coupling between the carbon cycles of the terrestrial and the fluvial environments. Recent identification of a substantial savanna carbon sink in wetter years in the recent past does not take into account the possibility of a substantial, rapid, lateral flux of carbon to rivers and back to the atmosphere. This article is protected by copyright. All rights reserved.

Description: Various remote-sensing methods are available to estimate soil moisture, but few address the fine spatial resolutions (e.g., 30 m grid cells) and root-zone depth requirements of agricultural and other similar applications. One approach that has been previously proposed to estimate fine-resolution soil moisture is to first estimate the evaporative fraction from an energy balance that is inferred from optical and thermal remote-sensing images (e.g., using the ReSET algorithm) and then estimate soil moisture through an empirical relationship to evaporative fraction. A similar approach has also been proposed to estimate the degree of saturation. The primary objective of this study is to evaluate these methods for estimating soil moisture and degree of saturation, particularly for a semiarid grassland with relatively dry conditions. Soil moisture was monitored at twenty-eight field locations in southeastern Colorado with herbaceous vegetation during the summer months of three years. In-situ soil moisture and degree of saturation observations are compared with estimates calculated from Landsat imagery using the ReSET algorithm. The in-situ observations suggest that the empirical relationships with evaporative fraction that have been proposed in previous studies typically provide overestimates of soil moisture and degree of saturation in this region. However, calibrated functions produce estimates with an accuracy that may be adequate for various applications. The estimates produced by this approach are more reliable for degree of saturation than for soil moisture, and the method is more successful at identifying temporal variability than spatial variability in degree of saturation for this region. This article is protected by copyright. All rights reserved.

Description: Groundwater movements in volcanic mountains and their effects on streamflow discharge and representative elementary area (REA) have remained largely unclear. We surveyed the discharge and chemical composition of spring and stream water in two catchments: the Hontani river (NR) catchment (6.6 km 2 ) and the Hosotani river (SR) catchment (4.0 km 2 ) at the southern part of Daisen volcano, Japan. Daisen volcano is a young volcano (17 × 10 3  years) at an early stage of erosion. Our study indicated that deep groundwater that moved through thick lava and pyroclastic flows and that could not be explained by shallow movements controlled by surface topography contributed dominantly to streamflow at larger catchment areas. At the NR catchment, the deep groundwater contribution clearly increased at a catchment boundary defined by an area of 3.0 km 2 and an elevation of 800 m. At the SR catchment, the contribution deep groundwater to the stream also increased suddenly at a boundary threshold of 2.0 km 2 and 700 m. Beyond these thresholds, the contributions of deep bedrock groundwater remained constant, indicating that the REA is between 2 and 3 km 2 at the observed area. These results indicate that the hydrological conditions of base flow were controlled mainly by the deep bedrock groundwater that moved through thick lava and pyroclastic flows in the undissected volcanic body of the upper part of the catchment. Our study demonstrates that deep and long groundwater movements via a deep bedrock layer including thick deposits of volcanic materials at the two catchments on Daisen volcano strongly determined streamflow discharge instead of the mixing of small-scale hydrological conditions. This article is protected by copyright. All rights reserved.

Description: We propose a novel technique for improving a long-term multi-step-ahead streamflow forecasts. A model based on wavelet decomposition and a multivariate Bayesian machine learning approach is developed for forecasting the streamflow three, six, nine and twelve months ahead simultaneously. The inputs of the model utilize only the past monthly streamflow records. They are decomposed into components formulated in terms of wavelet multiresolution analysis. It is shown that the model accuracy can be increased by using the wavelet boundary rule introduced in this study. A simulation study is performed to evaluate the effects of different wavelet boundary rules using synthetic and real streamflow data from the Yellowstone River in the Uinta Basin in Utah. The model based on the combination of wavelet and Bayesian machine learning regression techniques is compared to the wavelet and artificial neural networks based model. The robustness of the models is evaluated. This article is protected by copyright. All rights reserved.

Description: The combined use of water erosion models and geographic information systems (GIS) has facilitated soil loss estimation at the watershed scale. Tools such as the Geo-spatial interface for the Water Erosion Prediction Project (GeoWEPP) model provide a convenient spatially-distributed soil loss estimate but require discretization to identify hillslopes and channels. In GeoWEPP, the TOpographic PArameteriZation (TOPAZ) model is used as an automated procedure to extract a watershed boundary, hillslopes and channels from a digital elevation model (DEM). Previous studies in small watersheds have shown that the size of the hillslopes and the channel distribution affect the model estimates, but in large watersheds the effects on the soil loss estimates have yet to be tested. Therefore, the objective of this study was to evaluate the effect of discretization on the hillslope sediment yield estimates using GeoWEPP in two large watersheds (〉10 km 2 ). The watersheds were selected and discretized varying the TOPAZ parameters (CSA – Critical Source Area, and MSCL – Minimum Source Channel Length) in a 30-m resolution DEM. The drainage networks built with TOPAZ were compared among each other using the drainage density index. The results showed that the discretization affected hillslope sediment yield estimates and their spatial distribution more than total runoff. The drainage density index and the hillslope sediment yield were proportional but inversely related, thus soil loss estimates were highly affected by the spatial discretization. As a result of this analysis, a method to choose the CSA and MSCL values that generates the greatest fraction of hillslopes having profile lengths less than 200 m was developed. This slope length condition is particularly crucial when using the WEPP and GeoWEPP models, in order for them to produce realistic estimates of sheet and rill erosion. Finally, and as a result of this analysis, a more reliable method was developed for selecting the TOPAZ channel network parameters (CSA and MSCL).

Description: Linkages between the controls on surface storage and catchment streamflow response were examined in a wetland dominated basin in the Canadian Prairie Pothole region. Snowmelt, surface storage, water table elevation, atmospheric fluxes, and streamflow were monitored during spring snowmelt and summer in a 1 km 2 sub-catchment containing a semi-permanent pond complex connected via an intermittent stream. Snow accumulation in the basin in spring of the 2013 study year was the largest in the 24-year record. Rainfall totals in 2013 were close to the long term average, though June was an anomalously wet month. The water budget of the pond complex indicates that there was a significant subsurface contribution to surface storage. Activation of an effective transmission zone occurred between uplands and the stream network where the water table was located near the ground surface, which allowed significant lateral movement of subsurface water into the stream network. This was also important for maintaining and re-establishing surface connectivity and streamflow during rainfall events. The observed period of surface-water connectivity was one of the longest on record in the catchment due to unusually wet conditions; nevertheless, the results of this study have implications for how contributing area and runoff should be considered in monitoring and modelling studies in the region, as inclusion of more frequent and varied runoff processes will be essential to understanding changing streamflow regimes. This article is protected by copyright. All rights reserved.

Description: Stream-subsurface exchange plays a significant role in the fate and transport of contaminants in streams. It has been modeled explicitly by considering fundamental processes such as hydraulic exchange, colloid filtration, and contaminant interactions with streambed sediments and colloids. The models have been successfully applied to simulate the transport of inorganic metals and nutrients. In this study, laboratory experiments were conducted in a recirculating flume to investigate the exchange of a hydrophobic organic contaminant (HOC), p,p′ -DDE, between a stream and a quartz sand bed. A previously developed process-based multiphase exchange model was modified by accounting for the p,p′ -DDE kinetic adsorption to and desorption from the bed sediments/colloids and was applied to interpret the experimental results. Model input parameters were obtained by conducting independent small-scale batch experiments. Results indicate that the immobilization of p,p′ -DDE in the quartz sand bed can occur under representative natural stream conditions. The observed p,p′ -DDE exchange was successfully simulated by the process-based model. The model sensitivity analysis results show that the exchange of p,p′ -DDE can be sensitive to either the sediment sorption/desorption parameters or colloidal parameters depending on the experimental conditions tested. For the experimental conditions employed here, the effect of colloids on contaminant transport is expected to be minimal and the stream-subsurface exchange of p,p′ -DDE is dominated by the interaction of p,p′ -DDE with bed sediment. The work presented here contributes to a better mechanistic understanding of the complex transport process that HOCs undergo in natural streams, and to the development of reliable, predictive models for the assessment of impacted streams. This article is protected by copyright. All rights reserved.

Description: The physical and hydrological conditions in extracted peatlands often act as barriers to the regeneration of the keystone peat-forming genus Sphagnum . Although previous work has suggested that Sphagnum mosses regenerating on cutover peat surfaces quickly become vulnerable to water stress as the thickness of the regenerated layer increases, uncertainties regarding the storage and transmission properties of this layer and how these might evolve over time have made this assertion difficult to evaluate. This study investigates the hydrophysical properties and hydrological behaviour of regenerating Sphagnum layers ranging from 3-43 years in age using both field and laboratory methods. The 〉40 year old regenerated layers had significantly (p 〈 0.001) higher bulk density and retention capacity in the 5 cm thick basal layer directly overlying the cutover peat than the newer (〈10 year old) regenerated layers. Capillarity was a much stronger control on surficial water content (θ) than precipitation, which was poorly retained in the Sphagnum canopy, suggesting that regulation of water table position is an effective method of controlling θ as a means of optimizing productivity. In general, the θ sustained at a given water table position decreased as regenerated layer thickness increased. Analysis of water table position relative to the former cutover peat surface in different areas of the site suggests that the soil water dynamics of the 〉40 year old regenerated layers may be becoming increasingly similar to those of a natural bog peatland. This article is protected by copyright. All rights reserved.

Description: Nutrient concentrations in stream water, rainfall, throughfall, stem flow, surface flow and ground water were compared before, during, and after strip thinning (intensive 50%) in plantation forested watersheds in Tochigi, Japan. Influences were evaluated comparing four thinning-applied and two reference basins for one year before, six months during and one year after the thinning. Results show that this strip thinning significantly increased DTP, TP and DOC (DTP: 0.01 mg l -1 , TP: 0.04 mg l -1 , DOC: 0.53 mg l -1 ) during the thinning period and DTN and TN (DTN: 0.34 mg l -1 , TN: 0.46 mg l -1 ) after the thinning in stream waters relative to the unthinned basins. The increased phosphorus during thinning indicated ground disturbances by the strip thinning, with a concomitant increase in DOC. Changes in biotic and abiotic conditions resulted in increased nitrogen after the thinning, particularly in the dissolved pool. Changes in hydrological processes due to thinning, e.g. a change in flow distributions (less high nutrient stem flow and more low nutrient throughfall) and an increase in water discharge in stream water, possibly weakened the direct influences of thinning on nutrient concentrations. This article is protected by copyright. All rights reserved.

Description: Soil water repellency can impact soil hydrology, overland flow generation and associated soil losses. However, current hydrological models do not take it into account, which creates a challenge in repellency-prone regions. This work focused on the adaptation for soil water repellency of a daily water balance model. Repellency is estimated from soil moisture content using site-specific empirical relations, and used to limit maximum soil moisture. This model was developed and tested using c. 2 years of data from one long-unburned and two recently burned eucalypt plantations in northern Portugal, all of which showed strong seasonal soil water repellency cycles. Results indicated important improvements for the burned plantations, with the Nash-Sutcliffe efficiency increasing from -0.55 and -0.49 to 0.55 and 0.65. For the unburned site, model performance was already good without the modification and efficiency only improved slightly from 0.71 to 0.74, mostly due to the better simulation of delayed soil wetting after dry periods. Results suggested that even a simple approach to simulate soil water repellency can markedly improve the performance of hydrological models in eucalypt forests, especially after fire. This article is protected by copyright. All rights reserved.

Description: This study describes the use of linearly modulated optically stimulated luminescence (LM-OSL) to distinguish surface-soil derived sediments from those derived from channel bank erosion. LM-OSL signals from quartz extracted from fifteen surface-soil and five channel bank samples were analysed and compared to signals from samples collected from two downstream river sites. Discriminant analysis showed that the detrapping probabilities of fast, first slow and second slow components of the LM-OSL signal can be used to differentiate between the samples collected from the channel bank and surface-soil sources. We show that for each of these source end members these components are all normally distributed. These distributions are then used to estimate the relative contribution of surface-soil derived and channel bank derived sediment to the river bed sediments. The results indicate that channel bank derived sediments dominate the sediment sources at both sites, with 90.1 ± 3% and 91.9 ± 1.9% contributions. These results are in agreement with a previous study which used measurements of 137 Cs and 210 Pb ex fallout radionuclides to estimate the relative contribution from these two sources. This result shows that LM-OSL may be a useful method, at least in the studied catchment, to estimate the relative contribution of surface soil and channel erosion to river sediments. However, further research in different settings is required to test the difference of OSL signals in distinguishing these sediment sources. And if generally acceptable, this technique may provide an alternative to the use of fallout radionuclides for source tracing. This article is protected by copyright. All rights reserved.

Description: The vadose zone plays a crucial role in the water cycle for storing water, providing water to vegetation and transporting solutes or degrading contaminants. Earth scientists have long acknowledged the importance of the vadose zone and numerous methods have been developed to better understand and predict hydrological processes within this “critical zone”. For several decades, stable isotopes ( 18 O and 2 H) of pore water are used as environmental tracers to gain insights into vadose zone water movement and other processes. To determine the pore water stable isotopic composition various sampling procedures have been developed. We present the procedure and the accompanied advantages and drawbacks of each method. We further discuss possible opportunities and limitations regarding the scale of interest and the pore space that is sampled. The methodological review reveals that the choice of the sampling method is crucial for the interpretation of pore water stable isotopes in the vadose zone, but a thorough comparison between the different methods is yet missing. Spiking experiments, where water of known isotopic composition is added to oven-dried soil, have shown to be questionable, since the extracted water is usually depleted compared to the standard water. A comparative study analyzing soil samples with the recently developed direct water-vapor equilibration method and the widely used cryogenic extraction shows deviations, which can only be partly explained, but discloses the need for a more thorough experimental comparative study. Especially promising are developments of continuous isotope measurements based on laser-based spectrometry that will open up new opportunities of analyzing pore water isotopes with higher temporal and spatial resolutions, revealing new insights into hydrological processes across various temporal and spatial scales. This article is protected by copyright. All rights reserved.

Description: In the Vietnamese Mekong Delta (VMD), water levels at some stations have increased. However, the factors that cause this rise in the VMD have not been identified. We considered four factors that may have contributed to the water level rise: (1) increased runoff from upstream, (2) sea-level rise, (3) land subsidence, and (4) decrease in flood mitigation function due to construction of high dykes. We analyzed daily maximum and minimum water levels, and mean daily water levels from 24 monitoring stations from 1987 to 2006. Using daily and annual water level differences, we classified the delta into two groups; one is dominated by flows from upstream, while the other is tide-dominated. We then tested the trends of annual maximum and minimum water levels using the Mann-Kendall test, and identified the slope of the trend using the method of Sen. The areas of dyke construction were estimated using the Enhanced Vegetation Index (EVI) and Land Surface Water Index (LSWI) from the Moderate Resolution Imaging Spectroradiometer (MODIS) data. Results show (1) river inflow has little impact on rising water levels in the VMD, (2) the influence of high dykes on water level rise could not be quantified in this study, (3) both maximum and minimum water levels significantly increased in the tide-dominated area. Trend of annual minimum water level can be considered as the sum sea-level rise and land subsidence. Therefore, we attribute 6.05 mm year −1 (80%) to land subsidence and 1.42 mm year −1 (20%) to sea level rise, indicating inundations have been severe in the VMD, caused primarily by land subsidence. This article is protected by copyright. All rights reserved.

Description: Groundwater is the principal water resource in semi-arid and arid environments. Therefore, quantitative estimates of its replenishment rate are important for managing groundwater systems. In dry regions, karst outcrops often show enhanced recharge rates compared to other surface and sub-surface conditions. Areas with exposed karst features like sinkholes or open shafts allow point recharge, even from single rainfall events. Using the example of the As Sulb plateau in Saudi Arabia, this study introduces a cost-effective and robust method for recharge monitoring and modelling in karst outcrops. The measurement of discharge of a representative small catchment (4.0·10 4  m 2 ) into a sinkhole, and hence the direct recharge into the aquifer, was carried out with a time lapse camera. During the monitoring period of two rainy seasons (autumn 2012 to spring 2014) four recharge events were recorded. Afterwards, recharge data as well as proxy data about the drying of the sediment cover are used to set up a conceptual water balance model. The model was run for 17 years (1971 to 1986 and 2012 to 2014). Simulation results show highly variable seasonal recharge-precipitation ratios between zero and 0.27. In addition to the amount of seasonal precipitation, this ratio is influenced by the interannual distribution of rainfall events. Overall, an average annual groundwater recharge for the doline (sinkhole) catchment on As Sulb plateau of 5.1 mm a −1 has estimated for the simulation period. This article is protected by copyright. All rights reserved.

Description: This paper presents an analytical model for describing the tidal effects in a two-dimensional leaky confined aquifer system in an estuarine delta where ocean and river meet. This system has an unconfined aquifer on top, a confined aquifer on the bottom with an aquitard in between the two. The unconfined and confined aquifers interact with each other through leakage. It was assumed that the aquitard storage was negligible, and that the leakage was linearly proportional to the head difference between the unconfined and confined aquifers. This model's solution was based on the separation of variables method. Two existing solutions that deal with the head fluctuation in one-dimensional or two-dimensional leaky confined aquifers are shown as special cases in the present solution. Based on this new solution, the dynamic effect of the water table's fluctuations can be clearly explored, as well as the influence of leakage on the behavior of fluctuations in groundwater levels in the leaky aquifer system. This article is protected by copyright. All rights reserved.

Description: A portable Wavelength Scanned-Cavity Ring-Down Spectrometer (Picarro L2120) fitted with a diffusion sampler (DS-CRDS) was used for the first time to continuously measure δ 18 O and δ 2 H of stream water. The experiment took place during a storm event in a wet tropical agricultural catchment in north-eastern Australia. At a temporal resolution of one minute, the DS-CRDS measured 2160 δ 18 O and δ 2 H values continuously over a period of 36 hours with a precision of ± 0.08 and 0.5 ‰ for δ 18 O and δ 2 H, respectively. Four main advantages in using high temporal resolution stream δ 18 O and δ 2 H data during a storm event are highlighted from this study. Firstly, they enabled us to separate components of the hydrograph, which was not possible using high temporal resolution electrical conductivity data that represented changes in solute transfers during the storm event rather than physical hydrological processes. The results from the hydrograph separation confirm fast groundwater contribution to the stream, with the first 5 hours of increases in stream discharge comprising over 70% pre-event water. Secondly, the high temporal resolution stream δ 18 O and δ 2 H data allowed us to detect a short-lived reversal in stream isotopic values (δ 18 O increase by 0.4 ‰ over 9 minutes), which was observed immediately after the heavy rainfall period. Thirdly, δ 18 O values were used to calculate a time lag of 20 minutes between the physical and chemical stream responses during the storm event. Finally, the hydrograph separation highlights the role of event waters in the runoff transfers of herbicides and nutrients from this heavily cultivated catchment to the Great Barrier Reef. This article is protected by copyright. All rights reserved.

Description: Quantifying the linkages between vegetation disturbance by fire and the changes in hydrologic processes leading to post-fire erosional response remains a challenge. We measured the influence of fire severity, defined as vegetation disturbance (using a satellite derived vegetation disturbance index, VDI), landscape features that organize hydrologic flow pathways (relief and elongation ratios), and pre-fire vegetation type on the probability of the occurrence of post-fire gully rejuvenation (GR). We combined field surveys across 270 burned low-order catchments (112 occurrences of GR) and geospatial analysis to generate a probabilistic model through logistic regression. VDI alone discriminated well between catchments where GR did and did not occur (area under curve = 0.78, model accuracy = 0.72). The strong effect of vegetation disturbance on GR suggests that vegetation exerts a primary influence on the occurrence of infiltration excess runoff and post-fire erosion and that major gully erosion will not occur until fire consumes above ground biomass. Other topographic and local factors also influenced GR response, including catchment elongation, percent pre-fire shrub, mid-slope riparian vegetation, armored headwaters, firehose effects, and concentration of severe burn in source areas. These factors highlight the need to consider vegetation effects in concert with local topography and site conditions to understand the propensity for flow accumulation leading to GR. We present a process-based conceptual hydrologic model where vegetation loss from fire decreases rainfall attenuation and surface roughness, leading to accelerated flow accumulation and erosion; these effects are also influenced by interactions between fire severity and landscape structure. The VDI metric provides a continuous measure of vegetation disturbance and, when placed in a hydrologic context, may improve quantitative analysis of burned-area susceptibility to erosive rainfall, hazard prediction, ecological effects of fire, landform evolution, and sensitivity to climate change. This article is protected by copyright. All rights reserved.

Description: There is considerable interest in naturalizing flow regime on managed rivers to slow the spread of saltcedar ( Tamarix ramosissima ) invasion in southwestern United States or to preserve riparian forests dominated by saltcedar and other species in northwestern China. However, little is known about the responses of established saltcedar in water sources to frequent intra-annual fluctuation of water table resulting from this new, more dynamic flow regime. This study investigates how saltcedar at a riparian site in the middle reaches of the Heihe River, northwest China responds in water sources use to intra-annual water table fluctuations. Stable oxygen isotope was employed to determine accurate depth at which saltcedar obtains its water supply, and soil moisture monitoring was used to determine sources of plant-available soil water. We found that the primary zone of water uptake by saltcedar was stable at 25-60 cm depth, but the water sources used by saltcedar switched between groundwater and soil moisture with the water table fluctuations. Saltcedar derived its water from groundwater when water table was at depth less than 60 cm, but switched to soil moisture at 25-60 cm depth when water table declined. It is supposed that the well-developed clay layer at 60–80 cm depth constrained lateral roots of saltcedar to the soil layers above 60 cm, while the fine-textured soils at this site, which was periodically re-saturated by rising groundwater before the stored soil moisture had become depleted, provided an important water reservoir for saltcedar when groundwater dropped below the primary zone of fine roots. The root distribution of saltcedar may also be related to local groundwater history. The quick decline in water table in the early 1980s when the riparian saltcedar had established may strand its roots in the shallow unsaturated zone. We suggested that raising the water table periodically instead of maintaining it invariably above the rooting depth could sustain desired facultative phreatophytes while maximizing water deliveries. This article is protected by copyright. All rights reserved.

Description: Secondary circulation is the component of three-dimensional (3-D) flow in river channels perpendicular to the primary flow direction. Secondary circulation calculated from ADCP transects is sensitive to the calculation method, and is affected by the transect angle relative to the mean flow direction and variations in the flow direction along a transect. To quantify bounds on transect alignment relative to river flow for field data collection and examine tidal time-scale variability in secondary circulation, the 3-D hydrodynamic model UnTRIM was applied to simulate the hydrodynamics in the lower reach of the Sacramento River (CA, USA). Secondary circulation was calculated using the Rozovskii and the zero net discharge methods on repeated transects extracted from the model results in regions of both relatively uniform and complex flows. When the depth-averaged flow direction along a transect varied by more than about 5 degrees, occurring when the transect was as little as 10 to 20 degrees out of normal to the mean flow direction, the Rozovskii method produced more realistic secondary circulation than the zero net discharge method. Analysis indicated that ADCP transects should be within 20 degrees of perpendicular to the mean flow direction when calculating secondary circulation. Secondary circulation strength around two tidally influenced bends generally increased with increasing flow and broke down near slack water. However, the strength of the secondary circulation was not only a function of the flow magnitude, but also depended on the direction of the water flow, and the transect location relative to the river curvature, which varied with the tidal flow direction. This article is protected by copyright. All rights reserved.

Description: Flooding associated with tropical storms can cause extreme perturbations in riverine and coastal ecosystems. Measuring isotope variability of tropical storm events can help investigate the impacts of flooding. We measured the water isotope composition (δD and δ 18 O) of rain and associated floodwater collected during two storms and subsequent major and minor flooding events in the subtropical coast of eastern Australia. Compared to baseline regional rainfall isotope values of -15.0 ± 1.9‰ for δD and -3.3 ± 0.2‰ for δ 18 O, floodwater had lower values with -33.8 ± 2.5‰ δD and -5.1 ± 0.4‰ δ 18 O for the major flood and -29.4 ± 1.0‰ δD and -4.6 ± 0.1‰ δ 18 O for the minor flood. The low isotope composition of the floodwater was associated with the transport of large quantities of suspended sediments, with sediment loads 30 to 70 times larger than during base flow conditions. Floods carried up to 35% of the annual phosphorus and up to 208% of the currently calculated average annual nitrogen load of the Brisbane River. The dramatic changes caused by a rapid increase in discharge from 2 to 2,015 m 3 s -1 over two days in the major flood would have major consequences in riverine and coastal ecosystems of the region. These changes could potentially be traced using the isotope composition of the floodwaters. This article is protected by copyright. All rights reserved.

Description: Large river floods are a key water source for many lakes in fluvial periglacial settings. Where permeable sediments occur, the distribution of permafrost may play an important role in the routing of floodwaters across a floodplain. This relationship is explored for lakes in the discontinuous permafrost of Yukon Flats, interior Alaska, using an analysis that integrates satellite-derived gradients in water surface elevation, knowledge of hydrogeology, and hydrologic modeling. We observed gradients in water surface elevation between neighboring lakes ranging from 0.001 to 0.004. These high gradients, despite a ubiquitous layer of continuous shallow gravel across the flats, are consistent with limited groundwater flow across lake basins resulting from the presence of permafrost. Permafrost impedes the propagation of floodwaters in the shallow subsurface and constrains transmission to “fill-and-spill” over topographic depressions (surface sills), as we observed for the Twelvemile-Buddy Lake pair following a May 2013 ice-jam flood on the Yukon River. Model results indicate that permafrost table deepening of 1–11 m in gravel, depending on watershed geometry and subsurface properties, could shift important routing of floodwater to lakes from overland flow (fill-and-spill) to shallow groundwater flow (“fill-and-seep”). Such a shift is possible in the next several hundred years of ground surface warming, and may bring about more synchronous water level changes between neighboring lakes following large flood events. This relationship offers a potentially useful tool, well-suited to remote sensing, for identifying long-term changes in shallow groundwater flow resulting from thawing of permafrost. This article is protected by copyright. All rights reserved.

Description: Imperviousness, considered as a critical indicator of the hydrologic impacts of urbanization, has gained increasing attention both in the research field and in practice. However, the effectiveness of imperviousness on rainfall-runoff dynamics has not been fully determined in a fine spatiotemporal scale. In this study, 69 drainage subareas 〈 1 ha of a typical residential catchment in Beijing were selected to evaluate the hydrologic impacts of imperviousness, under a typical storm event with a three-year return period. Two metrics, total impervious area (TIA) and effective impervious area (EIA), were identified to represent the impervious characteristics of the selected subareas. Three runoff variables, total runoff depth (TR), peak runoff depth (PR), and lag time (LT), were simulated by using a validated hydrologic model. Regression analyses were developed to explore the quantitative associations between imperviousness and runoff variables. Then, three scenarios were established to test the applicability of the results in considering the different infiltration conditions. Our results showed that runoff variables are significantly related to imperviousness. However, the hydrologic performances of TIA and EIA were scale-dependent. Specifically, with finer spatial scale and the condition heavy rainfall, TIA rather than EIA was found to contribute more to TR and PR. EIA tended to have a greater impact on LT and showed a negative relationship. Moreover, the relative significance of TIA and EIA was maintained under the different infiltration conditions. These findings may provide potential implications for landscape and drainage design in urban areas, which help to mitigate the runoff risk. This article is protected by copyright. All rights reserved.

Description: The presence or absence of snow and its seasonal variability have profound effects on the climate, energy balance, and ecology in high-latitude and high-altitude regions of the world. While the timing, amount and duration of snow and its water content is of fundamental importance for several natural processes, it is also a convincing indicator of climate change in seasonally snow covered regions. The reduction in the amount and duration of the snow pack that already may have begun to happen across large geographic areas risks affecting several of the important ecosystem services that snow provides. Despite its fundamental role, snow and especially its water content is one of the least monitored climate variables in many seasonally snow-covered regions. Therefore we add complimentary long-term snow depth data (1980–2015) and snow water equivalent (SWE) data (1986–2015) to an already well established field research infrastructure – the Krycklan Catchment Study in northern Sweden. By combining digitized snow data to the already available long-term daily precipitation, runoff, and water quality information the ambition is to provide a more holistic database. While the primary purpose of this data release is to make the data readily available to any interested user, it also shows that the length of seasonal winter snow cover has shortened, especially because of its later arrival in the autumn, at the study site over the 35 year monitoring period. This article is protected by copyright. All rights reserved.

Description: Tree-ring-based reconstructions of paleo-hydrology have proved useful for better understanding the irregularities and extent of past climate changes, and therefore, for more effective water resources management. Despite considerable advances in the field, there still exist challenges that introduce significant uncertainties into paleo-reconstructions. This study outlines these challenges and address them by developing two themes: (1) the effect of temporal scaling on the strength of the relationship between the hydrologic variables, streamflow in this study, and tree growth rates, and (2) the reconstruction uncertainty of streamflow due to the dissimilarity or inconsistency in the pool of tree-ring chronologies (predictors in reconstruction) in a basin. Based on the insight gained, a methodology is developed to move beyond only relying on the annual hydrology-growth correlations, and to utilize additional information embedded in the annual time series at longer time scales (e.g, multi-year to decadal time scales). This methodology also generates an ensemble of streamflow reconstructions to formally account for uncertainty in the pool of chronology sites. The major headwater tributaries of the Saskatchewan River basin, the main source of surface water in the Canadian Prairie Provinces, are used as the case study. It is shown that the developed methodology explains the variance of streamflows to a larger extent than the conventional approach, and better preserves the persistence and variability of streamflows across time scales (Hurst-type behaviour). The resulting ensemble of paleo-hydrologic time series is able to more credibly pinpoint the timing and extent of past dry and wet periods, and provides a dynamic range of uncertainty in reconstruction. This range varies with time over the course of the reconstruction period, indicating that the utility of tree-ring chronologies for paleo-reconstruction differs for different time periods over the past several centuries in the history of the region. The proposed ensemble approach provides a credible range of multiple-century-long water availability scenarios that can be used for vulnerability assessment of the existing water infrastructure and improving water resources management. This article is protected by copyright. All rights reserved.

Description: Discharge in mountain streams may be a mixture of snowmelt, water from surface runoff, and deep return flow through valley bottom alluvia. We used δ 18 O and δ 2 H, solute concentrations, and 222 Rn to determine water sources of a headwater stream located at the McDonald Creek watershed, Glacier National Park, USA during summer recession flow period. We analyzed minimal water isotope ranges of -17.6 to -16.5‰ and -133 to -121‰ for δ 18 O and δ 2 H, respectively, potentially due to dominance of snow-derived water in the stream. Likewise, solute concentrations measured in the stream through the watershed showed minimal variation with little indication of sub-surface water input into the stream. However, we observed 222 Rn activities in the stream that ranged from 39 to 2,646 Bq/m 3 with the highest value measured in middle of the watershed associated with channel constriction corresponding to changes in local orientation of underlying rocks. Downstream from this point, 222 Rn activity decreased from 581 to 117 Bq/m 3 in a series of punctuated steps associated with small rapids and waterfalls that we hypothesized to cause radon degassing with a maximum predicted loss of 427 Bq/m 3 along a 400 m distance. Based on mass balance calculations using 222 Rn activity values, streamflow, and channel characteristics, we estimated that groundwater contributed between 0.3 and 29% of total flow. Overall, we estimated a 5.9% of groundwater contribution integrated for stream reach measured at McDonald Creek during recession flow period. Finally, a lower mean hyporheic flux of 14 m 3 /day was estimated compared to the groundwater flux of 70,710 m 3 /day. These assessments highlight the potential for radon as a conservative tracer that can be used to estimate sub-surface water contribution in mountain streams within a complex geologic setting. This article is protected by copyright. All rights reserved.

Description: Vegetative filter strips (VFS) can effectively trap sediment in overland flow, but little information is available on its performance in controlling high-concentration sediment and the runoff hydraulics in VFS. Flume experiments were conducted to investigate the sediment deposition, hydraulics of overland flow and their relationships in simulating VFS under a great range of sediment concentrations with four levels of vegetation cover (bare slope and 4, 11 and 17%) and two flow rates (15 and 30 L min −1 ). Sediment concentrations varied from 30 to 400 kg m −3 and slope gradient was 9°. Both the deposited sediment load and deposition efficiency in VFS increased as the vegetation cover increased. Sediment concentration had a positive effect on the deposited load but no effect on deposition efficiency. A lower flow rate corresponded to greater deposition efficiency but had little effect on deposited load. Flow velocities decreased as vegetation cover increased. Sediment concentration had a negative effect on the mean velocity but no effect on surface velocity. Hydraulic resistance increased as the vegetation cover and sediment concentration increased. Sediment deposition efficiency had a much more pronounced relationship with overland flow hydraulics compared with deposited load, especially with the mean flow velocity, and there was a power relationship between them. Flow regime also affected the sediment deposition efficiency, and the efficiency was much higher under subcritical than supercritical flow. The results will be useful for the design of VFS and the control of sediment flowing into rivers in areas with serious soil erosion.

Description: Recent changes in land use and new industrial developments in river basins around the world have raised serious concerns regarding future water quality. Assessment of regional nutrient export from river basins are needed in order to identify main nutrient sources and hydrological variables involved. The Red Deer River catchment in Alberta, Canada has been chosen as a test case to assess total phosphorus export regionally. Agricultural and livestock activities in the catchment have increased manure production and the supply of fertilizer to crops. Oil and gas exploitation has also increased the risk of surface water and groundwater contamination. The rapid population growth has not only lead to increases in water consumption and wastewater discharges, but also to further construction of transportation infrastructure and the expansion of new urban developments. This has imposed hydraulic controls on waterways, affecting the catchment hydrology and changing the dynamics of sediment and nutrient transport. River ecosystems are not exempt from the negative effects of water quality deterioration. Downstream from the city of Red Deer, the physiology and reproduction habits of native fish species have changed due to high riverine productivity. Although improvements have been made to surface water quality standards by Alberta Environment, further research is needed in order to identify major nutrient sources and quantify nutrient export. The SPAtially Referenced Regression On Watershed (SPARROW) model has therefore been used in this study of the Red Deer River catchment to assess regional water quality, in order to describe the spatial and temporal patterns of the processes that affect water quality. The model is suitable for interpreting monitoring data sets that suffer from network sparseness, bias and basin heterogeneity. Ultimately, the model could provide improved information to environmental agencies to guide future water quality management practices and policies.

Description: The White method is a simple but the most frequently applied approach to estimate groundwater evapotranspiration ( ET g ) from groundwater level diurnal signals. Due to a lack of direct measurements of ET g , it is difficult to evaluate the performance of the White method, particularly in field environments with variable groundwater fluctuations. A two-year field observation in a hyper-arid riparian tamarisk ( Tamarix spp.) stand with deep groundwater depth in the lower Tarim River basin of China was conducted to measure the surface evapotranspiration ( ET s ) and groundwater table. The performance of the White method and the influences of the variable groundwater fluctuations on the determinations of the specific yield ( S y ) and recharge rate of groundwater ( r ) in the White method were investigated. The results showed that the readily available S y determined by Loheide's method was feasible, but must be finely determined based on the soil textures in the layers in which the groundwater level fluctuated. A constant S y value for a defined porous medium could be assumed regardless of the discharge or recharge processes of groundwater. The time span of 0 A.M. to 6 A.M. for r estimation for the White method worked best than other time spans. A two-day moving average of r values further enhance ET g estimation . Slight effects of environmental or anthropogenic disturbances on the diurnal fluctuations of groundwater level did not influence the ET g estimations by the White method. Our results provide valuable references to the application of the White method for estimating daily ET g in desert riparian forests with deep groundwater depth. This article is protected by copyright. All rights reserved.

Description: This study analyzed age-related water use dynamics across three temperate conifer forest plantations (aged 11-, 39-, and 74-years old, as of 2013, henceforth referred to as TP02, TP74, and TP39, where the last two digits represent the year of planting) in southern Ontario, Canada from 2008 to 2013. Eddy covariance-measured mean evapotranspiration over the growing season (April-October) was 438 ± 19, 392 ± 19 and 323 ± 25 mm at TP39, TP74 and TP02 respectively. Daytime bulk surface conductance was highest and most variable at the TP39 site (8.5 ± 4.0 mm s −1 ), followed by the TP74 (7.0 ± 2.8 mm s −1 ) and TP02 (5.4 ± 2.5 mm s −1 ) sites. Evapotranspiration at all the forests was sensitive to air temperature and also tended to decrease with increasing atmospheric dryness. The youngest forest's evapotranspiration was most conservative, which led to an increase in water use efficiency throughout the study period, in particular during drought events. The oldest forest was the least restrictive in its water use during drought, which led to lower water use efficiency during such events as compared to the younger forests. The oldest forest was thinned in early 2012, where about 1/3 of trees were commercially harvested. No significant change in evapotranspiration or water use efficiency was observed at this site following thinning, however daytime bulk surface conductance declined. Our results suggested that changes in stand structure with forest ageing, such as reduction in stem density and increase in sapwood area, were responsible for differences in soil water demand during drought and non-drought periods, leading to differences in forest water use. Hence, forest age, due to its structural implications, is an important control on the stand-level water use efficiency and forests’ response to drought events. Our study suggested that younger forests may be best suited to maximize growth and carbon uptake efficiency under rising air temperatures and increasing precipitation variability as predicted by climate models for eastern North America. This article is protected by copyright. All rights reserved.

Description: Physically-based distributed and semi-distributed hydrological models have become some of the primary tools for water resources studies and management over the past decades owing to increased computational capabilities and advances in data measurements. Representation of the existing heterogeneity in nature still remains one of the main challenges in these models, and is accomplished primarily via watershed discretization, subdividing watersheds into hydrologically similar land parcels. Discretization decisions in distributed modeling studies are often ad hoc and determined with little or no quantitative analysis to support these decisions. In this work, we present a quantitative methodology for assessing alternative watershed discretization schemes in terms of their corresponding model performance in ungauged basins. The effect of the computational time spent for calibrating each scheme (calibration budget) is considered as part of the assessment. Here, these schemes differ in how they represent landscape attributes, and range from a simple lumped scheme to more complex ones by adding spatial land cover and then soil information. The methodology was demonstrated using the Modélisation Environmentale–Surface et Hydrologie (MESH) model as applied to the Nottawasaga river basin in Ontario, Canada. Results reveal that model performance in ungauged basins depends upon the location of the validation sub-basin (i.e. upstream or downstream) with respect to the calibration sub-basins. Also, using a more complex scheme did not necessarily lead to improved performance in validation, when constrained by calibration budget. Therefore, the calibration budget also should be considered as a factor in the assessment process. This methodology was also implemented using a shorter sub-period for calibration, which leads to substantial computational saving. Results of the sub-period test were promising and consistent particularly when sufficient budget is spent to calibrate the model. Other strategies utilized for reducing the computational burden of the proposed analyses are also discussed in this study. This article is protected by copyright. All rights reserved.

Description: Floods are the most frequent natural disaster, causing more loss of life and property than any other in the United States. Floods also strongly influence the structure and function of watersheds, stream channels, and aquatic ecosystems. The Pacific Northwest is particularly vulnerable to climatically-driven changes in flood frequency and magnitude, because snowpacks that strongly influence flood generation are near the freezing point and thus sensitive to small changes in temperature. To improve predictions of future flooding potential and inform strategies to adapt to these changes, we mapped the sensitivity of landscapes to changes in peak flows due to climate warming across Oregon and Washington. We first developed principal component based models for predicting peak flows across a range of recurrence intervals (2-, 10-, 25-, 50-, and 100-years) based on historical instantaneous peak flow data from 1000 gaged watersheds in Oregon and Washington. Key predictors of peak flows included drainage area and principal component scores for climate, land cover, soil, and topographic metrics. We then used these regression models to predict future peak flows by perturbing the climate variables based on future climate projections (2020s, 2040s, and 2080s) for the A1B emission scenario. For each recurrence interval, peak flow sensitivities were computed as the ratio of future to current peak flow magnitudes. Our analysis suggests that temperature induced changes in snowpack dynamics will result in large (〉30-40%) increases in peak flow magnitude in some areas, principally the Cascades, Olympics, and Blue Mountains and parts of the western edge of the Rocky Mountains. Flood generation processes in lower elevation areas are less likely to be affected, but some of these areas may be impacted by floodwaters from upstream. These results can assist land, water, and infrastructure managers in identifying watersheds and resources that are particularly vulnerable to increased peak flows and developing plans to increase their resilience. This article is protected by copyright. All rights reserved.

Description: Numerous studies investigated the influence of abiotic (meteorological conditions) and biotic factors (tree characteristics) on stemflow generation. Though these studies identified the variables that influence stemflow volumes in simply structured forests, the combination of tree characteristics that allows a robust prediction of stemflow volumes in species-rich forests is not well known. Many hydrological applications, however, require at least a rough estimate of stemflow volumes based on the characteristics of a forest stand. The need for robust predictions of stemflow motivated us to investigate the relationships between tree characteristics and stemflow volumes in a species-rich tropical forest located in central Panama. Based on a sampling setup consisting of 10 rainfall collectors, 300 throughfall samplers, and 60 stemflow collectors and cumulated data comprising 26 rain events, we derive three main findings. First, stemflow represents a minor hydrological component in the studied 1 ha forest patch (1.0 % of cumulated rainfall). Second, in the studied species-rich forest, single tree characteristics are only weakly related to stemflow volumes. The influence of multiple tree parameters (e.g. crown diameter, presence of large epiphytes, and inclination of branches) and the dependencies among these parameters require a multivariate approach to understand the generation of stemflow. Third, predicting stemflow in species-rich forests based on tree parameters is a difficult task. Although our best model can capture the variation in stemflow to some degree, a critical validation reveals that the model cannot provide robust predictions of stemflow. A reanalysis of data from previous studies in species-rich forests corroborates this finding. Based on these results and considering that for most hydrological applications stemflow is only one parameter among others to estimate, we advocate using the base model, i.e. the mean of the stemflow data, to quantify stemflow volumes for a given study area. Studies in species-rich forests that wish to obtain predictions of stemflow based on tree parameters probably need to conduct a much more extensive sampling than currently implemented by most studies. This article is protected by copyright. All rights reserved.

Description: This study analyzes the stable isotopic compositions of hydrogen and oxygen ( δ 2 H, δ 18 O) in montane meteoric waters including precipitation and stream water of central Taiwan to identify hydrological processes in montane catchments. Results of precipitation demonstrate that monsoon and altitude effects are two principal processes affecting δ and deuterium excess ( d E ) values of inland precipitation in central Taiwan. Furthermore, slope and intercept values of summer and winter local Meteoric Water Line (MWL) are modified by secondary evaporation effects such as moisture recycling and raindrop evaporation. Additionally, stream water's results indicate that differences in δ values among stream waters reflect isotopic altitude effect whereby lower values are more evident in stream water originating from high-elevation catchments than low-elevation catchments. Comparison of the isotopic results between precipitation and stream water indicates that summer precipitation containing recycled moisture is the most important water source for the studied stream waters and indicates that catchment effect and base-flow contribution are the two major hydrological processes affecting mountain stream hydrology. The hydrological processes identified by the isotopic study re-stress the important role of forests in mountain hydrology. This article is protected by copyright. All rights reserved.

Description: Groundwater ridging is the rapid rise of a shallow water table during a rainfall event, in an environment where, in the pre-event period, the capillary fringe extends to the ground surface. Groundwater ridging is widely cited to account for the observed significant appearance of pre-event water in a stream stormflow hydrograph. Various hypotheses have been advanced to explain the groundwater ridging mechanism; and most recently, from a field study site in South Africa, an energy hypothesis was proposed, which explains that groundwater ridging water table rise is as a result of rapid introduction and transmission of additional pressure head into the capillary fringe from an intense rainfall at the ground surface. However, there is need for further analysis and evidence from other field study sites to confirm and support this newly proposed energy hypothesis. The objectives of this paper are, therefore, to: review previous observations on groundwater ridging, from other study sites, in order to deduce evidence of the newly proposed energy hypothesis; present and evaluate a one-dimensional diffusion mathematical model that can simulate groundwater ridging water table rise, based on the newly proposed energy hypothesis; and evaluate the importance of a capillary fringe in streamflow generation. Analysis of previous observations from other study sites generally indicated that the rate of groundwater ridging water table rise is directly related to the rainfall intensity, hence, confirming and agreeing with the newly proposed energy hypothesis. Additionally, theoretical results by the mathematical model agreed fairly well with the field results observed under natural rainfall, confirming that the rapidly rainfall-induced energy is diffusively transmitted downwards through pore water, elevating the pressure head at every depth. The results in this study also support the concept of a three-end member stream stormflow hydrograph, and contribute to the explanation of how catchments can store water for long periods, but then release it rapidly during storm events. This article is protected by copyright. All rights reserved.

Description: The most popular practice for analyzing nonstationarity of flood series is to use a fixed single-type probability distribution incorporated with the time-varying moments. However, the type of probability distribution could be both complex because of distinct flood populations and time-varying under changing environments. To allow the investigation of this complex nature, the time-varying two-component mixture distributions (TTMD) method is proposed in this study by considering the time variations of not only the moments of its component distributions but also the weighting coefficients. Having identified the existence of mixed flood populations based on circular statistics, the proposed TTMD was applied to model the annual maximum flood series (AMFS) of two stations in the Weihe River basin (WRB), with the model parameters calibrated by the meta-heuristic maximum likelihood (MHML) method. The performance of TTMD was evaluated by different diagnostic plots and indexes and compared with stationary single-type distributions, stationary mixture distributions and time-varying single-type distributions. The results highlighted the advantages of TTMD with physically-based covariates for both stations. Besides, the optimal TTMD models were considered to be capable of settling the issue of nonstationarity and capturing the mixed flood populations satisfactorily. This article is protected by copyright. All rights reserved.

Description: While the effects of land use change in urban areas have been widely examined, the combined effects of climate and land use change on the quality of urban and urbanizing streams have received much less attention. We describe a modeling framework that is applicable to the evaluation of potential changes in urban water quality and associated hydrologic changes in response to ongoing climate and landscape alteration. The grid-based spatially distributed model, DHSVM-WQ, is an outgrowth of the Distributed Hydrology-Soil-Vegetation Model (DHSVM) that incorporates modules for assessing hydrology and water quality in urbanized watersheds at a high spatial and temporal resolution. DHSVM-WQ simulates surface runoff quality and in-stream processes that control the transport of nonpoint-source (NPS) pollutants into urban streams. We configure DHSVM-WQ for three partially urbanized catchments in the Puget Sound region to evaluate the water quality responses to current conditions and projected changes in climate and/or land use over the next century. Here we focus on total suspended solids (TSS) and total phosphorus (TP) from nonpoint sources (runoff), as well as stream temperature. The projection of future land use is characterized by a combination of densification in existing urban or partially urban areas, and expansion of the urban footprint. The climate change scenarios consist of individual and concurrent changes in temperature and precipitation. Future precipitation is projected to increase in winter and decrease in summer, while future temperature is projected to increase throughout the year. Our results show that urbanization has a much greater effect than climate change on both the magnitude and seasonal variability of streamflow, TSS and TP loads largely due to substantially increased streamflow, and particularly winter flow peaks. Water temperature is more sensitive to climate warming scenarios than to urbanization and precipitation changes. Future urbanization and climate change together are predicted to significantly increase annual mean streamflow (up to 55%), water temperature (up to 1.9 °C), TSS load (up to 182%), and TP load (up to 74%). This article is protected by copyright. All rights reserved.

Description: The resolution of a digital elevation model (DEM) is a crucial factor in watershed hydrologic and environmental modeling. DEM resolution can cause significant variability in the representation of surface topography, which further affects quantification of hydrologic connectivity and simulation of hydrologic processes. The objective of this study is to examine the effects of DEM resolution on (1) surface microtopographic characteristics, (2) hydrologic connectivity, and (3) the spatial and temporal variations of hydrologic processes. A puddle-to-puddle (P2P) modeling system was utilized for surface delineation and modeling of the P2P overland flow dynamics, surface runoff, infiltration, and unsaturated flow for nine DEM resolution scenarios of a field plot surface. Comparisons of the nine modeling scenarios demonstrated that coarser DEM resolutions tended to eliminate topographic features, reduce surface depression storage, and strengthen hydrologic connectivity and surface runoff. We found that reduction in maximum depression storage and maximum ponding area was as high as 97.56% and 76.36%, respectively, as the DEM grid size increased from 2 cm to 80 cm. The paired t-test and fractal analysis demonstrated the existence of a threshold DEM resolution (10 cm for the field plot), within which the DEM-based hydrologic modeling was effective and acceptable. The effects of DEM resolution were further evaluated for a larger surface in the Prairie Pothole Region (PPR) subjected to observed rainfall events. It was found that simulations based on coarser resolution DEMs (〉10 m) tended to overestimate ponded areas and underestimate runoff discharge peaks. The simulated peak discharge from the PPR surface reduced by approximately 50% as the DEM resolution changed from 2 m to 90 m. Fractal analysis results elucidated scale dependency of hydrologic and topographic processes. In particular, scale analysis highlighted a unique constant-threshold-power relationship between DEM scale and topographic and hydrologic parameters/variables. Not only does this finding allow one to identify threshold DEM, but also further develop functional relationships for scaling to achieve valid topographic characterization as well as effective and efficient hydrologic modeling. This article is protected by copyright. All rights reserved.

Description: This paper presents a methodology for hydrograph separation in mountain watersheds, which aims at identifying flow sources among ungauged headwater sub-catchments through a combination of observed streamflow and data on natural tracers including isotope and dissolved solids. Daily summer and bi-daily spring season water samples obtained at the outlet of the Juncal River Basin in the Andes of Central Chile were analyzed for all major ions as well as stable water isotopes, δ 18 O and δD. Additionally, various samples from rain, snow, surface streams and exfiltrating subsurface water (springs), were sampled throughout the catchment. A principal component analysis (PCA) was performed in order to address cross-correlation in the tracer dataset, reduce the dimensionality of the problem and to uncover patterns of variability. Potential sources were identified in a two-component U-space that explains 94% of the observed tracer variability at the catchment outlet. Hydrograph separation was performed through an Informative - Bayesian model. Our results indicate that the Juncal Norte Glacier headwater sub-catchment contributed at least 50% of summer flows at the Juncal River Basin outlet during the 2011–12 water year (a hydrologically dry period in the Region), even though it accounts for only 27% of the basin area. Our study confirms the value of combining solute and isotope information for estimating source contributions in complex hydrologic systems, and provides insights regarding experimental design in high-elevation semi-arid catchments. The findings of this study can be useful for evaluating modeling studies of the hydrological consequences of the rapid decrease in glacier cover observed in this region, by providing insights into the origin of river water in basins with little hydrometeorological information.

Description: The average flow of Silver Springs, one of the largest magnitude springs in central Florida, declined 32% from 2000 to 2012. The average groundwater head in the springshed declined 0.14 m, and the spring pool altitude increased 0.24 m during the same period. This paper presents a novel explanation of the spring flow recession curve for Silver Springs using Torricelli model, which uses the groundwater head at a sentinel well, the spring pool altitude, and the net recharge to groundwater. The effective springshed area and net recharge (defined as recharge minus groundwater pumping and evapotranspiration) were estimated based on the observed recession slopes for spring flow, groundwater head, and spring pool altitude. The results show that the effective springshed area continuously declined since 1989, and the net recharge declined since the 1970s with a significant drop in 2002. Subsequent to 2002, the net recharge increased modestly but not to the levels prior to the 1990s. The reduction in net recharge was caused by changes in hydroclimatic conditions including precipitation and temperature, along with groundwater withdrawals, which contributed to the declined spring flow.

Description: Discharge simulation from snow-dominated catchments seems to be an easy task. Any spatially-explicit precipitation-runoff model coupled to a temperature-index snow model generally yields simulations that mimic well the observed daily discharges. The robustness of such models is, however, questionable: in presence of strong annual discharge cycles, small model residuals do not guarantee high explanatory power of the underlying model. This paper proposes a methodology for snow hydrological model identification within a limits-of-acceptability framework, where acceptable model simulations are the ones that reproduce a set of signatures within an a priori specified range. The signatures proposed here namely include the relationship between the air temperature regime and the discharge regime, a new snow hydrology signature that can be readily transferred to other Alpine settings. The discriminatory power of all analyzed signatures is assessed with a new measure of their discriminatory power in the model prediction domain. The value of the proposed snow hydrology signatures and of the limits-of-acceptability approach is demonstrated for the Dischma river in Switzerland, whose discharge shows a strong temporal variability of hydrologic forcing conditions over the last 30 years. The signature-based model identification for this case study leads to the surprising conclusion that the observed discharge data contains a multi-year period that cannot be reproduced with the model at hand. This model-data mismatch might well result from a yet to be identified problem with the discharge observations, which would have been difficult to detect in a classical residual-based model identification approach. Overall, the detailed results for this case study underline the robustness of the limits-of-acceptability approach in the presence of error-prone observations if it is applied in combination with relatively robust signatures. Future work will show whether snow hydrology signatures and their limits-of-acceptability can be regionalized to ungauged catchments, which would make this model selection approach particularly powerful for Alpine environments. This article is protected by copyright. All rights reserved.

Description: Changes in land use and land cover are major drivers of hydrological alteration in the tropical Andes. However, quantifying their impacts is fraught with difficulties because of the extreme diversity in meteorological boundary conditions, which contrasts strongly with the lack of knowledge about local hydrological processes. Although local studies have reduced data scarcity in certain regions, the complexity of the tropical Andes poses a big challenge to regional hydrological prediction. This study analyses data generated from a participatory monitoring network of 25 headwater catchments covering three of the major Andean biomes ( páramo , jalca , and puna ), and link their hydrological responses to main types of human interventions (cultivation, afforestation and grazing). A paired catchment setup was implemented to evaluate the impacts of change using a “trading space-for-time” approach. Catchments were selected based on regional representativeness and contrasting land use types. Precipitation and discharge have been monitored and analysed at high temporal resolution for a time period between 1 and 5 years. The observed catchment responses clearly reflect the extraordinarily wide spectrum of hydrological processes of the tropical Andes. They range from perennially humid páramos in Ecuador and northern Peru with extremely large specific discharge and baseflows, to highly seasonal, flashy catchments in the drier punas of southern Peru and Bolivia. The impacts of land use are similarly diverse and their magnitudes are a function of catchment properties, original and replacement vegetation, and management type. Cultivation and afforestation consistently affect the entire range of discharges, particularly low flows. The impacts of grazing are more variable, but have the largest effect on the catchment hydrological regulation. Overall, anthropogenic interventions result in increased streamflow variability and significant reductions in catchment regulation capacity and water yield, irrespective of the hydrological properties of the original biome. This article is protected by copyright. All rights reserved.

Description: Temporal streamflow variability in an inland hydrologic station and temporal trends and frequency changes at three weather stations in a semiarid river basin located in Loess Plateau, China were detected using linear regression, Mann-Kendall (MK) analysis, and wavelet transform methods. Double Cumulative Curve and Ordered Clustering were used to identify the hydrological periods of upper Sang-kan (USK) basin between 1957 and 2012. The results indicate that: (1) precipitation in the USK basin over the study period did not show any trend while the temperature showed a significant increase; (2) streamflow flowing out of the USK basin indicated a significant decrease; (3) two distinct hydrological periods; the “natural period” from 1957 to 1984 and the “human impact period” from 1985 to 2012 were present; and (4) the contribution of climate change and human activities to reduce the streamflow was 36.9% and 63.1%, respectively. The results indicate that human activities may be contributing to a decrease in streamflow in the USK basin. This article is protected by copyright. All rights reserved.

Description: Hydrochemistry methods were used to decipher the weathering and geochemical processes controlling solute acquisition of river waters in the dry season in the middle Loess Plateau (MLP), one of the most severely eroded areas and turbid riverine systems in the world. River waters were neutral to slightly alkaline with pH varying from 7.6 to 9.6. The total dissolved solids (TDS) decreased from northwest to southeast with a mean value of 804 mg/L, much higher than the global average and other large rivers in China. Ternary diagram showed that river waters were dominated by Na + , HCO 3 - , and Cl - with the main water-type of HCO 3 - -Cl - -Na + . Saturation index (SI) values, Mg 2+ , Ca 2+ , and HCO 3 - analyses indicated the preferential Ca 2+ removal by calcite precipitation. Gibbs plots and stoichiometry plots indicated that the dissolved solutes were mainly derived from rock weathering with minor anthropogenic and atmospheric inputs. Samples in the northwestern basin are also influenced by evaporation. A forward model of mass budget calculation showed that, owing to high soluble characteristics, evaporite dissolution was a major feature of river waters and contributed 41% to the total dissolved cations on average, while carbonate and silicate weathering contributed 28%, and 25% on average, respectively. Beside evaporite dissolution, cation exchange is also responsible for the high concentrations of Na + in river water. Spatial variations showed that evaporite dissolution and silicate weathering were higher in the northern basin, whereas carbonate weathering was higher in the southern basin. Different from most rivers in the world, the physical erosion rates (varying from 117.7 to 4116.6 t/km 2 y) are much higher than the chemical weathering rates (varying from 3.54 to 6.76 t/km 2 y) in the MLP due to the loose structure of loess and poor vegetation in the basin. In the future, studies on comparison of water geochemistry in different seasons and on influence of different types of land use and soil salinization on water geochemistry, denudation rates, and water quality should be strengthened in the MLP. These results shed some lights on processes responsible for modern loess weathering and also indicated the importance of time-series sampling strategy for river water chemistry.

Description: The Tagus River Basin is an ultimately important water source for hydropower production, urban and agricultural water supply in Spain and Portugal. Growing electricity and water supply demands, over-regulation of the river and construction of new dams, as well as large inter-basin and intra-basin water transfers aggravated by strong natural variability of climate in the catchment have already imposed significant pressures on the river. The substantial reduction of discharge is observed already now, and projected climatic change is expected to alter the water budget of the catchment further. In this study we address the effects of projected climate change on the water resources availability in the Tagus River Basin, and influence of potential changes on hydropower generation of the three important reservoirs in the basin. The catchment scale, process-based eco-hydrological model Soil and Water Integrated Model (SWIM) was set up, calibrated and validated for the entire Tagus River Basin, taking into account fifteen large reservoirs in the catchment. The future climate projections were selected from those generated within the Inter-Sectoral Impact Model Intercomparison Project. They include five bias-corrected climatic datasets for the region, obtained from Global Circulation Models runs under two emissions scenario – moderate and extreme ones, and covered the whole century. The results show a strong agreement among model runs in projecting substantial decrease of discharge of the Tagus River discharge and, consequently, a strong decrease in hydropower production under both future climate scenarios.

Description: Simulations from hydrological models are affected by potentially large uncertainties stemming from various sources, including model parameters and observational uncertainty in the input/output data. Understanding the relative importance of such sources of uncertainty is essential to support model calibration, validation and diagnostic evaluation, and to prioritise efforts for uncertainty reduction. It can also support the identification of ‘disinformative data’ whose values are the consequence of measurement errors or inadequate observations. Sensitivity Analysis (SA) provides the theoretical framework and the numerical tools to quantify the relative contribution of different sources of uncertainty to the variability of the model outputs. In traditional applications of GSA, model outputs are aggregations of the full set of a simulated variable. For example, many GSA applications use a performance metric (e.g. the root mean squared error) as model output that aggregates the distances of a simulated time series to available observations. This aggregation of propagated uncertainties prior to GSA may lead to a significant loss of information and may cover up local behaviour that could be of great interest. Time-varying sensitivity analysis (TVSA), where the aggregation and SA are repeated at different time-steps, is a viable option to reduce this loss of information. In this work, we use TVSA to address two questions: [1] Can we distinguish between the relative importance of parameter uncertainty versus data uncertainty in time? [2] Do these influences change in catchments with different characteristics? To our knowledge, the results present one of the first quantitative investigation on the relative importance of parameter and data uncertainty across time. We find that the approach is capable of separating influential periods across data and parameter uncertainties, while also highlighting significant differences between the catchments analysed.

Description: In this paper we assess the performance of the catchment model SIMCAT, to predict nitrate and soluble reactive phosphorus (SRP) concentrations against four monitoring regimes with different spatial and temporal sampling frequencies. The GLUE uncertainty framework is used, along with a general sensitivity analysis to understand relative parameter sensitivity. Improvements to model calibration are explored by introducing more detailed process representation using the INCA water quality model, driven by the E-HYPE hydrological model. The results show how targeted sampling of headwater watercourses upstream of point discharges is essential for calibrating diffuse loads, and can exert a strong influence on the whole-catchment model performance. Further downstream, if the point discharges and loads are accurately represented, then the improvement in the catchment scale model performance is relatively small as more calibration points are added or frequency is increased. The higher order, dynamic model INCA-P, which incorporates sediment and biotic interaction, resulted in improved whole-catchment performance over SIMCAT, although there are still large epistemic uncertainties from land phase export coefficients and runoff. However, the very large sampling errors in routine monitoring make it difficult to invest confidence in the modelling, especially since we know phosphorous transport to be very episodic and driven by high flow conditions for which there are few samples. The environmental modelling community seems to have been stuck in this position for some time, and whilst it is useful to use an uncertainty framework to highlight these issues, it has not widely been adopted, perhaps because there is no clear mechanism to allow uncertainties to influence investment decisions. This raises the question as to whether it might better place a cost on uncertainty, and use this to drive more data collection or improved models, before making investment decisions concerning for example mitigation strategies. This article is protected by copyright. All rights reserved.

Description: This study was undertaken to evaluate the effects of climatic variability on inter-annual variations in each component of evapotranspiration ( ET ) and the total ET in a temperate coniferous forest in Japan. We conducted eddy covariance flux and meteorological measurements for 7 years and parameterized a one-dimensional multi-layer biosphere-atmosphere model (Kosugi et al., 2006) that partitions ET to transpiration ( Tr ), wet-canopy evaporation ( E wet ), and soil evaporation ( E soil ). The model was validated with the observed flux data. Using the model, the components of ET were estimated for the 7 years. Annual precipitation, ET , Tr , E wet , and E soil over the 7 years were 1536 ± 334 mm, 752 ± 29 mm, 425 ± 37 mm, 219 ± 34 mm, and 108 ± 10 mm, respectively. The maximum inter-annual fluctuation of observed ET was 64 mm with a coefficient of variance (CV) of 2.7%, in contrast to relatively large year-to-year variations in annual rainfall (CV = 20.1%). Tr was related to the vapor pressure deficit, incoming radiation, and air temperature with relatively small inter-annual variations (CV = 8.2%). E soil (CV = 8.6%) was related mainly to the vapor pressure deficit. E wet was related to precipitation with large inter-annual variations (CV = 14.3%) because of the variability in precipitation. The variations in E wet were counterbalanced by the variations in Tr and E soil , producing the small inter-annual variations in total ET . This article is protected by copyright. All rights reserved.

Description: Backflow, the temporary reversal of discharge at the outlet of a lake, is an important mechanism controlling flow and transport in many connected river-lake systems. This study used statistical methods to examine long-term variations and primary causal factors of backflow from the Yangtze River to a laterally connected, large floodplain lake (Poyang Lake, China). Additionally, the effects of backflow on the lake hydrology were explored using a physically based hydrodynamic model and a particle-tracking model. Although backflow into Poyang Lake occurs frequently, with an average of 16 backflow events per year, and varies greatly in magnitude between years, statistical analysis indicates that both the frequency and magnitude of backflow reduced significantly during 2001-2010 relative to the previous period of 1960-2000. The ratio of Poyang Lake catchment inflows to Yangtze River discharge can be used as an indication of the daily occurrence of backflow, which is most likely to occur during periods when this ratio is lower than 5%. Statistical analysis also indicates that the Yangtze River discharge is the main controlling factor of backflow during July-October, rather than catchment inflows to the lake. Hydrodynamic modeling reveals that, in general, backflow disturbs the normal northward water flow direction in Poyang Lake, and transports mass ~20 km southward into the lake. The effects of backflow on flow direction, water velocities and water levels propagate to virtually its upstream extremity. The current study represents a first attempt to explore backflow and causal factors for a highly dynamic lake-floodplain system. An improved understanding of Poyang Lake backflow is critical for guiding future strategies to manage the lake, its water quality and ecosystem value, given proposals to modify the lake-river connectivity. This article is protected by copyright. All rights reserved.

Description: This paper reports the first results on δ 18 O and δ 2 H analysis of precipitations, cave drip waters, and groundwaters from sites in Mallorca (Balearic Islands, western Mediterranean), a key region for paleoclimate studies. Understanding the isotopic variability and the sources of moisture in modern climate systems is required to develop speleothem isotope-based climate reconstructions. The stable isotopic composition of precipitation was analyzed in samples collected between March 2012 and March 2013. The values are in the range reported by GNIP Palma station. Based on these results, the local meteoric water line δ 2 H = 7.9 (±0.3) δ 18 O + 10.8 (±2.5) was derived, with slightly lower slope than GMWL. The results help tracking two main sources of air masses affecting the study sites: rain events with the highest δ 18 O values (〉 –5 ‰) originate over the Mediterranean Sea, whereas the more depleted samples (〈 –8 ‰) are sourced in the North Atlantic region. The back trajectory analysis and deuterium excess values, ranging from 0.4 to 18.4 ‰, further support our findings. To assess the isotopic variation across the island, water samples from eight caves were collected. The δ 18 O values range between –6.9 and –1.6 ‰. With one exception (Artà), the isotopic composition of waters in caves located along the coast (Drac, Vallgornera, Cala Varques, Tancada, and Son Sant Martí) indicates Mediterranean-sourced moisture masses. By contrast, the drip water δ 18 O values for inland caves (Campanet, ses Rates Pinyades) or developed under a thick (〉50 m) limestone cap (Artà) exhibit more negative values. A well-homogenized aquifer supplied by rainwaters of both origins is clearly indicated by groundwater δ 18 O values, which show to be within 2.4 ‰ of the unweighted arithmetic mean of –7.4 ‰. Although limited, the isotopic data presented here constitute the baseline for future studies using speleothem δ 18 O records for western Mediterranean paleoclimate reconstructions. This article is protected by copyright. All rights reserved.

Description: Extremely high precipitation occurs in the Southern Alps of New Zealand, associated with both orographic enhancement and synoptic-scale weather processes. In this study, we test the hypothesis that atmospheric rivers (ARs) are a key driver of floods in the Southern Alps of New Zealand. Vertically-integrated water vapour and horizontal water vapour transport, and atmospheric circulation, are investigated concurrently with major floods on the Waitaki River (a major South Island river). Analysis of the largest eight annual maximum floods between 1979 and 2012 indicates that all are associated with ARs. Geopotential height fields reveal that these ARs are located in slow eastward moving extra-tropical cyclones, with high pressure to the northeast of New Zealand. The confirmation of ARs as a contributor to Waitaki flooding indicates the need for their further exploration to better understand South Island hydrometeorological extremes. This article is protected by copyright. All rights reserved.

Description: Land cover has been increasingly recognized as an important factor affecting hydrologic processes at the basin and regional level. Therefore, improved understanding of how land cover change affects hydrologic systems is needed for better management of water resources. The objective of this study is to investigate the effects of land cover change on the duration and severity of high and low flows by using the Soil Water Assessment Tool (SWAT) model, Bayesian Model Averaging (BMA) and copulas. Two basins dominated by different land cover in the Ohio River basin are used as study area in this study. Two historic land covers from the 1950s and 1990s are considered as input to the SWAT model, thereby investigating the hydrologic high and low flow response of different land cover conditions of these two basins. The relationships between the duration and severity of both low and high flow are defined by applying the copula method; changes in the frequency of the duration and severity are investigated. The results show that land cover changes affect both the duration and severity of both high and low flows. An increase in forest area leads to a decrease in the duration and severity during both high and low flows, but its impact is highest during extreme flows. The results also show that the land cover changes have had significant influences on changes in the joint return periods of duration and severity of low and high flows. While this study sheds light on the role of land cover change on severity and duration of high and low flow conditions, more studies using various land cover conditions and climate types are required in order to draw more reliable conclusions in future. This article is protected by copyright. All rights reserved.

Description: The mountain headwater Bow River at Banff, Alberta, Canada was subject to a large flood in June 2013, over which considerable debate has ensued regarding its probability of occurrence. It is therefore instructive to consider what information long term streamflow discharge records provide about environmental change in the Upper Bow River basin above Banff. Though protected as part of Banff National Park, since 1885, the basin has experienced considerable climate and land cover changes, each of which has the potential to impact observations, and hence the interpretations of flood probability. The Bow River at Banff hydrometric station is one of Canada's longest operating reference hydrological basin network stations and so has great value for assessing changes in flow regime over time. Furthermore, the station measures a river that provides an extremely important water supply for Calgary and irrigation district downstream and so is of great interest for assessing regional water security. These records were examined for changes in several flood attributes and to determine whether flow changes may have been related to landscape change within the basin as caused by forest fires, conversion from grasslands to forest with fire suppression, and regional climate variations and/or trends. Floods in the Upper Bow River are generated by both snowmelt and rain-on-snow (ROS) events, the latter type which include floods events generated by spatially and temporally large storms such as occurred in 2013. The two types of floods also have different frequency characteristics. Snowmelt and ROS flood attributes were not correlated significantly with any climate index or with burned area except that snowmelt event duration correlated negatively to the Pacific Decadal Oscillation. While there is a significant negative trend in all floods over the past 100 years, when separated based on generating process, neither snowmelt floods nor large ROS floods associated with mesoscale storms show any trends over time. Despite extensive changes to the landscape of the basin and in within the climate system, the flood regime remains unchanged, something identified at smaller scales in the region but never at larger scales.

Description: A theoretical, dimensionless rainfall-runoff model was used to simulate the discharge of Wulongdong spring in western Hubei Province, South China. The single parameter (time constant τ) in the model is easy to obtain by fitting the recession rate of the observed hydrographs. The model was scaled by simply matching the total annual flow volume of the model to the observed value. Annual distribution of actual evapotranspiration was embedded in the model input to calculate the accumulated deficit of soil moisture before each rain event. Hourly precipitation input data performed better than daily data, defining τ of 0.85 days, and returning a Nash-Sutcliffe efficiency of 0.89 and the RMSE of 0.07. This model offers an effective way to simulate the discharge of karst springs that respond sensitively to rainfall events. The model parameters of a successful simulation can be used to estimate the recharge area and indicate the intrinsic response time of the basin.

Description: The change of hydrological regimes may cause impacts on human and natural system. Therefore, investigation of hydrologic alteration induced by climate change is essential for preparing timely proper adaptation to the changes. This study employed 24 climate projections from the Coupled Model Intercomparison Project Phase 5 (CMIP5) under Representative Concentration Pathway (RCP) 4.5 scenario. The climate projections were downscaled at a station-spacing for seven Korean catchments by a statistical downscaling method that preserves a long-term trend in climate projections. Using an ensemble of future hydrologic projections simulated by three conceptual rainfall-runoff models (GR4J, IHACRES, and Sacramento models), we calculated Hydrologic Alteration Factors (HAFs) to investigate degrees of variations in Indicators of Hydrologic Alteration (IHAs) derived from the hydrologic projections. The results showed that the seven catchments had similar trend in terms of the HAFs for the 24 IHAs. Given that more frequent severe floods and droughts were projected over Korean catchments, sound water supply strategies are definitely required to adapt to the alteration of streamflow. A wide range of HAFs between rainfall-runoff models for each catchment was detected by large variations in the magnitude of HAFs with the hydrologic models and the difference could be the hydrologic prediction uncertainty. There were no-consistent tendency in the order of HAFs between the hydrologic models. In addition, we found that the alterations of hydrologic regimes by climate change are smaller as the size of catchment is larger.

Description: The lower stretch of the Vistula is the most ice-jammed river section on the North European Plain. Since 1982, the structure of hanging dams has been studied by means of a mechanical non-core sampler. In this article, a selected of field research results of the hanging dams structure and the degree of filling of the cross-section with ice obtained during surface ice jam events in the years 1995-2014 are presented, along with an explanation of their causes. Surface ice jams occurred during spring snowmelt surges and ice breakups and also during freeze-up and ice covered periods. Their main cause was changes in the river flow and were also affected by anthropogenic sources. A characteristic feature of the analysed cross-sections was the considerable share of the underhanging ice dam's firm accumulation with ice floes, when the cross-section would be filled with ice in excess of 70%. In most cases, due to low river discharge, there was no substantial flooding damage. This article is protected by copyright. All rights reserved.

Description: Uruguay has stimulated the development of its forest sector since the promulgation of Forest Law N° 15 939 in December of 1987. Nevertheless, the substitution of natural grasslands with forest plantations for industrial use has raised concerns regarding hydrological processes of groundwater recharge and water consumption involving evapotranspiration. The purpose of this study is to assess the effects of this substitution approach on water resources. Input data were collected from two small experimental watersheds of roughly 100-200 hectares located in western Uruguay. The watersheds are characterized by Eucalyptus Globulus ssp. Maidenni and natural grasslands for cattle use. Total rainfall, stream discharge, rainfall redistribution, soil water content and groundwater level data were collected. Groundwater recharge was estimated from water table fluctuations and from groundwater contributions to base flows. Seasonal and annual water budgets were computed from October of 2006 to September of 2014 to evaluate changes in the hydrological processes. The data show a decrease in annual specific discharge of roughly 17% for mean hydrological years and no conclusive effects on annual groundwater recharge in the forested watershed relative to the reference pasture watershed. Reduced annual specific discharge is equivalent to the mean annual interception. The computed actual annual evapotranspiration is consistent with international catchment measurements. Reduction rates vary seasonally and according to accumulated rainfall and its temporary distribution. The degree of specific discharge decline is particularly high for drier autumns and winters (32 to 28%) when the corresponding rainfall varies from 275 to 400 mm. These results are of relevance for water resources management efforts, as water uses downstream can be affected. These findings, based on a study period dominated by anomalous wet springs and summers and by dry autumns and winters, oppose earlier results based on 34 years of rainfall and discharge data drawn from Uruguayan large basins. Copyright © 2016 John Wiley & Sons, Ltd.

Description: River water temperature is a key physical variable controlling several chemical, biological and ecological processes. Its reliable prediction is a main issue in many environmental applications, which however is hampered by data scarcity, when using data-demanding deterministic models, and modeling limitations, when using simpler statistical models. In this work we test a suite of models belonging to air2stream family (Toffolon and Piccolroaz, 2015), which are characterized by a hybrid formulation that combines a physical derivation of the key equation with a stochastic calibration of parameters. The air2stream models rely solely on air temperature and streamflow, and are of similar complexity as standard statistical models. The performances of the different versions of air2stream in predicting river water temperature are compared with those of the most common statistical models typically used in the literature. To this aim, a dataset of 38 Swiss rivers is used, which includes rivers classified into four different categories according to their hydrological characteristics: low-land natural rivers, lake outlets, snow-fed rivers, and regulated rivers. The results of the analysis provide practical indications regarding the type of model that is most suitable to simulate river water temperature across different time scales (from daily to seasonal) and for different hydrological regimes. A model intercomparison exercise suggests that the family of air2stream hybrid models generally outperforms statistical models, while cross-validation conducted over a 30-year period indicates that they can be suitably adopted for long-term analyses.

Description: As large, high-severity forest fires increase and snowpacks become more vulnerable to climate change across the western US, it is important to understand post-fire disturbance impacts on snow hydrology. Here, we examine, quantify, parameterize, model, and assess the post-fire radiative forcing effects on snow to improve hydrologic modeling of snow-dominated watersheds having experienced severe forest fires. Following a 2011 high-severity forest fire in the Oregon Cascades, we measured snow albedo, monitored snow and micrometeorological conditions, sampled snow surface debris, and modeled snowpack energy and mass balance in adjacent burned and unburned forest sites. For three winters following the fire, charred debris in the burned forest reduced snow albedo, accelerated snow albedo decay, and increased snowmelt rates thereby advancing the date of snow disappearance compared to the unburned forest. We demonstrate a new parameterization of post-fire snow albedo as a function of days-since-snowfall and net snowpack energy balance using an empirically-based exponential decay function. Incorporating our new post-fire snow albedo decay parameterization in a spatially-distributed energy and mass balance snow model, we show significantly improved predictions of snow cover duration and spatial variability of snow water equivalent across the burned forest, particularly during the late snowmelt period. Field measurements, snow model results, and remote sensing data demonstrate that charred forests increase the radiative forcing to snow and advance the timing of snow disappearance for several years following fire. This article is protected by copyright. All rights reserved.

Description: A spatially distributed representation of basin hydrology and transport processes in hydrologic models facilitates the identification of critical source areas and the placement of management and conservation measures. Floodplains are critical landscape features that differ from neighboring uplands in terms of their hydrological processes and functions. Accordingly, an important step in watershed modeling is the representation of floodplain and upland areas within a watershed. The aim of this study is (1) to evaluate four floodplain-upland delineation methods that use readily available topographic data (topographic wetness index, slope position, uniform flood stage, and variable flood stage) with regard to their suitability for hydrological models and (2) to introduce an evaluation scheme for the delineated landscape units. The methods are tested in three U.S. watersheds ranging in size from 334 to 629 km 2 with different climatic, hydrological and geomorphological characteristics. Evaluation of the landscape delineation methods includes visual comparisons, error matrices (i.e. cross-tabulations of delineated versus reference data), and geometric accuracy metrics. Reference data was obtained from SSURGO (Soil Survey Geographic database) and FEMA (Federal Emergency Management Agency) flood maps. Results suggest that the slope position and the variable flood stage method work very well in all three watersheds. Overall percentages of floodplain and upland areas allocated correctly were obtained by comparing delineated and reference data. Values range from 83 to 93 % for the slope position and from 80 to 95 % for the variable flood stage method. Future studies will incorporate these two floodplain-upland delineation methods into the subwatershed-based hydrologic model SWAT to improve the representation of hydrological processes within floodplain and upland areas. This article is protected by copyright. All rights reserved.

Description: There has been a great deal of research interest regarding changes in flow path/runoff source with increases in catchment area. However, there have been very few quantitative studies taking subscale variability and convergence of flow path/runoff source into account, especially in relation to headwater catchments. This study was performed to elucidate how the contributions and discharge rates of subsurface water (water in the soil layer) and groundwater (water in fractured bedrock) aggregate and change with catchment area increase, and to elucidate whether the spatial variability of the discharge rate of groundwater determines the spatial variability of stream discharge or groundwater contribution. The study area was a 5-km 2 forested headwater catchment in Japan. We measured stream discharge at 113 points and water chemistry at 159 points under base flow conditions. End-member mixing analysis (EMMA) was used to separate stream water into subsurface water and groundwater. The contributions of both subsurface water and groundwater had large variability below 1 km 2 . The contribution of subsurface water decreased markedly, while that of groundwater increased markedly, with increases in catchment area. The specific discharge of subsurface water showed a large degree of variability and decreased with catchment area below 0.1 km 2 , becoming almost constant above 0.1 km 2 . The specific discharge of groundwater showed large variability below 1 km 2 and increased with catchment area. These results indicated that the variabilities of stream discharge and groundwater contribution corresponded well with the variability of the discharge rate of groundwater. However, below 0.1 km 2 , it was necessary to consider variations in the discharge rates of both subsurface water and groundwater.

Description: In arctic and sub-arctic environments, mercury (Hg), and more specifically toxic methylmercury (MeHg) is of growing concern to local communities because of its accumulation in fish. In these regions, there is particular interest in the potential mobilisation of atmospherically-deposited Hg sequestered in permafrost that is thawing at unprecedented rates. Permafrost thaw and the resulting ground surface subsidence transforms forested peat plateaus into treeless and permafrost-free thermokarst wetlands where inorganic Hg released from the thawed permafrost and draining from the surrounding peat plateaus may be transformed to MeHg. This study begins to characterise the spatial distribution of MeHg in a peat plateau-thermokarst wetland complex, a feature that prevails throughout the wetland-dominated southern margin of thawing discontinuous permafrost in Canada's Northwest Territories. We measured pore water total Hg, MeHg, dissolved organic matter (DOM) characteristics and general water chemistry parameters to evaluate the role of permafrost thaw on the pattern of water chemistry. A gradient in vegetation composition, water chemistry and DOM characteristics followed a toposequence from the ombrotrophic bogs near the crest of the complex to the poor fens at its downslope margins. We found that pore waters in poor fens contained elevated levels of MeHg and the water draining from these features had dissolved MeHg concentrations 4.5 to 14.5 times higher than the water draining from the bogs. It was determined through analysis of historical aerial images that the poor fens in the toposequence had formed relatively recently (early 1970s) as a result of permafrost thaw. Differences between the fens and bogs are likely to be due to their differences in groundwater function and this suggests that permafrost thaw in this landscape can result in hotspots for Hg methylation that are hydrologically connected to downstream ecosystems.

Description: The objective of this work was to build a prognostic water flow model and potentially toxic elements (lead, cadmium, zinc) transport model in the unsaturated zone. Research was conducted in the catchment area of Kosnica regional wellfield, where the unsaturated zone is characterised by Fluvisol. Lower sorption capacities were determined in the first horizons for all three potentially toxic elements. Correlation coefficient of the measured and simulated values of tracer concentration is 0.58 for the AC horizon and 0.84 for the 2C/C1 horizon. Based on calibrated water flow and transport parameters, a prognostic water flow model and potentially toxic elements (lead, cadmium, zinc) transport model in the unsaturated zone was built. In case of an accidental spill of potentially toxic elements with concentrations of 1000 mg/L, the risk of contamination of the aquifer is present.

Description: Headwater storage-discharge remains one of the least understood processes, and there is renewed interest in the S-Q relation. How well can the S-Q relation be interpreted mechanistically using geometric factors? In this paper, the hillslope storage Boussinesq (hsB) and hillslope storage kinematic wave (hsKW) equation were adopted to guide the theoretical derivations. Analytical solutions were derived based on the hsKW equation for nine idealized hillslope aquifers, which were subdivided into two groups, i.e., hillslope aquifers with exponential hillslope width function (C1) and hillslope aquifers with Gaussian hillslope width function (C2). We found that analytical expressions of the S-Q relation can be derived for C1 hillslope aquifers. For more compound hillslope aquifers, i.e., C2, no explicit S-Q relation can be obtained. The whole subsurface recession after a rainstorm is simulated by applying the initial saturation condition. We found that the simulated S-Q processes can be characterized by a two-phase recession, i.e., quick and slow recession. The time (t b ) at the dividing point of the quick and slow recessions depends on the geometric factors, such as the plan and profile curvature. In the quick recession for C1, many of the S-Q curves can be described as linear or quasi-linear functions, which indicate that linear reservoir models can be applied approximately for recession simulations. However, during the slow recession phase of C1 and during the whole recession of C2, the S-Q relations are highly non-linear. Finally, we compared the hsKW and hsB models for simulating subsurface water recession after a rainstorm event in a real-world headwater catchment (G5) in China. Through comparison of the recession slope curves, we found that the simulated results of the models employing the Gaussian hillslope width function match the observed hydrograph. The results indicated that appropriate organization of the hillslope geometric factors enhances our ability to make storage-discharge predictions. This article is protected by copyright. All rights reserved.

Description: Spatial and temporal variation in wet canopy conditions following precipitation events can influence processes such as transpiration and photosynthesis, which can be further enhanced as upper canopy leaves dry more rapidly than the understory following each event. As part of a larger study aimed at improving land-surface modeling of evapotranspiration processes in wet tropical forests, we compared transpiration among trees with exposed and shaded crowns under both wet and dry canopy conditions in central Costa Rica, which has an average 4200 mm annual rainfall. Transpiration was estimated for 5 months using 43 sap flux sensors in 8 dominant, 10 midstory, and 8 suppressed trees in a mature forest stand surrounding a 40-m tower equipped with micrometeorological sensors. Dominant trees were 13% of the plot's trees and contributed around 76% to total transpiration at this site, whereas midstory and suppressed trees contributed 18% and 5%, respectively. After accounting for vapor pressure deficit and solar radiation, leaf wetness was a significant driver of sap flux, reducing it by as much as 28%. Under dry conditions, sap flux rates ( J s ) of dominant trees were similar to midstory trees and were almost double that of suppressed trees. On wet days, all trees had similarly low J s . As expected, semi-dry conditions (dry upper canopy) led to higher J s in dominant trees than midstory, which had wetter leaves, but semi-dry conditions only reduced total stand transpiration slightly and did not change the relative proportion of transpiration from dominant and midstory. Therefore, models that better capture forest stand wet-dry canopy dynamics and individual tree water use strategies, are needed to improve accuracy of predictions of water recycling over tropical forests.

Description: Snow and glacial melt processes are an important part of the Himalayan water balance. Correct quantification of melt runoff processes is necessary to understand the region's vulnerability to climate change. This paper describes in detail an application of conceptual GR4J hydrological model in the Tamor catchment in Eastern Nepal using typical elevation band and degree-day factor approaches to model Himalayan snow and glacial melt processes. The model aims to provide a simple model that meets most water planning applications. The paper contributes a model conceptualization (GR4JSG) that enables coarse evaluation of modeled snow extents against remotely sensed MODIS snow extent. Novel aspects include the glacial store in GR4JSG and examination of how the parameters controlling snow and glacial stores correlate with existing parameters of GR4J. The model is calibrated using a Bayesian Monte Carlo Markov Chain method against observed streamflow for one glaciated catchment with reliable data. Evaluation of the modelled streamflow with observed streamflow gave Nash Sutcliffe Efficiency of 0.88 and Percent Bias of 〈4%. Comparison of the modelled snow extents with MODIS gave R 2 of 0.46, with calibration against streamflow only. The contribution of melt runoff to total discharge from the catchment is 14-16% across different experiments. The model is highly sensitive to rainfall and temperature data, which suffer from known problems and biases, for example due to stations being located predominantly in valleys and at lower elevations. Testing of the model in other Himalayan catchments may reveal additional limitations. This article is protected by copyright. All rights reserved.

Description: The impact of climate change on the behaviour of intensity-duration-frequency (IDF) curves is critical to the estimation of design storms, and thus to the safe design of drainage infrastructure. The present study develops a regional time trend methodology that detects the impact of climate change on extreme precipitation from 1960 to 2010. The regional time trend linear regression method is fitted to different durations of annual maximum precipitation intensities derived from multiple sites in Ontario, Canada. The results show the relationship between climate change and increased extreme precipitation in this province. The regional trend analysis demonstrates, under nonstationary conditions arising from climate change, that the intensity of extreme precipitation increased decennially between 1.25 per cent for the 30 min storm and 1.82 per cent for the 24-hour storm. A comparison of the results with a regional Mann-Kendall test validates the found regional time-trend results. The results are employed to extrapolate the IDF curves temporally and spatially for future decades across the province. The results of the regional time trend assessment helps with the establishment of new safety margins for infrastructure design in Ontario.

Description: Changes in potential evapotranspiration and surface runoff can have profound implications for hydrological processes in arid and semiarid regions. In this study, we investigated the response of hydrological processes to climate change in Upper Heihe River Basin (UHRB) in Northwest China for the period from 1981 to 2010. We used agronomic, climatic and hydrological data to drive the SWAT (Soil and Water Assessment Tool) model for changes in potential evapotranspiration (ET 0 ) and surface runoff and the driving factors in the study area. The results showed that increasing autumn temperature increased snow melt, resulting in increased surface runoff, especially in September and October. The spatial distribution of annual runoff was different from that of seasonal runoff, with the highest runoff in Yeniugou River, followed by Babaohe River and then the tributaries in the northern of the basin. There was no evaporation paradox at annual and seasonal time scales, and annual ET 0 was driven mainly by wind speed. ET 0 was driven by relative humidity in spring, sunshine hour duration in autumn, and both sunshine hour duration and relative humility in summer. Surface runoff was controlled by temperature in spring and winter and by precipitation in summer (flood season). Although surface runoff increased in autumn with increasing temperature, it depended on rainfall in September and on temperature in October and November. This article is protected by copyright. All rights reserved.

Description: As a consequence of the remote location of the Andean páramo is knowledge on their hydrologic functioning limited, notwithstanding this alpine tundra ecosystem act as water towers for a large fraction of the society. Given the harsh environmental conditions in this region is year-round monitoring cumbersome, and it would be beneficially if the monitoring needed for the understanding of the rainfall-runoff response could be limited in time. To identify the hydrological response and the effect of temporal monitoring a nested (n = 7) hydrological monitoring network was set up in the Zhurucay catchment (7.53 km 2 ), south Ecuador. The research questions were: (1) can event sampling provide similar information in comparison to continuous monitoring, and (2) if so, how many events are needed to achieve a similar degree of information? A subset of 34 rainfall runoff events was compared to monthly values derived from a continuous monitoring scheme from December 2010 to November 2013. Land cover and physiographic characteristics were correlated with eleven hydrological indices. Results show that despite some distinct differences between event and continuous sampling, both datasets reveal similar information; more in particular the monitoring of a single event in the rainy season provides the same information as continuous monitoring, while during the dry season 10 events ought to be monitored. This article is protected by copyright. All rights reserved.

Description: We apply the process-based, distributed TOPKAPI-ETH glacio-hydrological model to glacierized catchment (19% glacierized) in the semiarid Andes of central Chile. The semiarid Andes provides vital freshwater resources to valleys in Chile and Argentina, but only few glacio-hydrological modelling studies have been conducted and its dominant hydrological processes remain poorly understood. The catchment contains two debris-free glaciers reaching down to 3900 m asl (Bello and Yeso glaciers) and one debris-covered avalanche-fed glacier reaching to 3200 m asl (Piramide Glacier). Our main objective is to compare the mass balance and runoff contributions of both glacier types under current climatic conditions. We use a unique dataset of field measurements collected over two ablation seasons combined with the distributed TOPKAPI-ETH model that includes physically-oriented parameterizations of snow and ice ablation, gravitational distribution of snow, snow albedo evolution and the ablation of debris-covered ice. Model outputs indicate that while the mass balance of Bello and Yeso glaciers is mostly explained by temperature gradients, the Piramide Glacier mass balance is governed by debris thickness and avalanches and has a clear non-linear profile with elevation as a result. Despite the thermal insulation effect of the debris cover, the mass balance and contribution to runoff from debris-free and debris-covered glaciers is similar in magnitude, mainly due to elevation differences. However, runoff contributions are distinct in time and seasonality with ice melt starting approximately four weeks earlier from the debris-covered glacier, what is of relevance for water resources management. At the catchment scale, snowmelt is the dominant contributor to runoff during both years. However, during the driest year of our simulations, ice melt contributes 42 ± 8% and 67 ± 6% of the annual and summer runoff, respectively. Sensitivity analyses show that runoff is most sensitive to temperature and precipitation gradients, melt factors and debris cover thickness. This article is protected by copyright. All rights reserved.

Description: Nine mid- to high-latitude headwater catchments —part of the North-Watch (Northern Watershed Ecosystem Response to Climate Change) program — were used to analyze threshold response to rainfall and snowmelt-driven events, and link the different responses to the catchment characteristics of the nine sites. The North-Watch data include daily time-series of various lengths of multiple variables such as air temperature, precipitation and discharge. Rainfall and meltwater inputs were differentiated using a degree-day snowmelt approach. Distinct hydrological events were identified, and precipitation-runoff response curves were visually assessed. Results showed that eight of nine catchments showed runoff initiation thresholds and effective precipitation input thresholds. For rainfall-triggered events, catchment hydroclimatic and physical characteristics (e.g., mean annual air temperature, median flow path distance to the stream, median sub-catchment area) were strong predictors of threshold strength. For snowmelt-driven events, however, thresholds and their governing factors controlling precipitation-runoff response were difficult to identify. The variability in catchments responses to snowmelt was not fully explained by runoff initiation thresholds and input magnitude thresholds. The quantification of input intensity thresholds (e.g., snow melting and permafrost thawing rates) is likely required for an adequate characterization of nonlinear spring runoff generation in such northern environments. This article is protected by copyright. All rights reserved.

Description: Rivers in the Mediterranean region often exhibit an intermittent character. An understanding and classification of the flow regimes of these rivers is needed, since flow patterns control both physicochemical and biological processes. This paper reports an attempt to classify flow regimes in Mediterranean rivers based on hydrological variables extracted from discharge time series. Long-term discharge records from 60 rivers within the Mediterranean region were analysed in order to classify the streams into different flow regime groups. Hydrological indices (HIs) were derived for each stream and principal component analysis (PCA) then applied to these indices to identify subsets of HIs describing the major sources of variations, whilst simultaneously minimising redundancy. PCA was performed for two groups of streams (perennial and temporary) and for all streams combined. The results show that whereas perennial streams are mainly described by high flow indices, temporary streams are described by duration, variability and predictability indices. Agglomerative cluster analysis based on hydrological indices identified six groups of rivers classified according to differences in intermittency and variability. A methodology allowing such a classification for ungauged catchments was also tested. Broad-scale catchment characteristics based on digital elevation, climate, soil and land-use data were derived for each long-term station where these data were available. By using stepwise multiple regression analysis, statistically significant relationships were fitted linking the three selected hydrological variables (mean annual number of zero flow days, predictability and flashiness) to the catchment characteristics. The method provides a means of simplifying the complexity of river systems and is thus useful for river basin management. This article is protected by copyright. All rights reserved.

Description: This study uses stable isotopes and major ions to examine the seasonal evolution of penitentes on the surface of Tapado Glacier, in the Norte Chico region of the Chilean Andes. A snow pit was sampled in November 2011 and penitentes were sampled during the summer (December 2011 and January 2012). The major ion load of the winter snowpack is dominated by Ca 2+ (60%), SO 4 2− (16%) and NO 3 − (13%) and there is little influence from marine air masses at the site, with most SO 4 2− , Mg 2+ , Ca 2+ , and Na + , derived from non-sea salt sources. During the early ablation season we observe increases in stable isotope ratios and major ion concentrations (particularly lithic ions Na + , Mg 2+ , and Ca 2+ ) in the upper reaches of penitentes, which is attributed to sublimation and the aeolian deposition of dust particles. In the late-summer, melt replaces sublimation as the dominant ablation process, and results in smoothing of the stable isotope profile and the elution of major ions within the penitente snow and ice matrix. This article is protected by copyright. All rights reserved.

Description: Mathematical modelling is a well-accepted framework to evaluate the effects of wetlands on stream flow and watershed hydrology in general. Although the integration of wetland modules into a distributed hydrological model represents a cost-effective way to make this assessment, the added value brought by landscape-specific modules to a model's ability to replicate basic hydrograph characteristics remains unclear. The objectives of this paper were to: (i) present the adaptation of PHYSITEL (a GIS) to parameterize isolated and riparian wetlands; (ii) describe the integration of specific isolated wetland and riparian wetland modules into HYDROTEL, a distributed hydrological model; and (iii) evaluate the performance of the updated modelling platform with respect to the capacity of replicating various hydrograph characteristics. To achieve this, two sets of simulations were performed (with and without wetland modules) and the added-value was assessed at three river segments of the Becancour River watershed, Quebec, Canada, using six general goodness-of-fit indicators (GOFIs) and fourteen water flow criteria (WFC). A sensitivity analysis of the wetland module parameters was performed to characterize their impact on stream flows of the modelled watershed. Results of this study indicate that: (i) integration of specific wetland modules can slightly increase the capacity of HYDROTEL to replicate basic hydrograph characteristics and (ii) the updated modelling platform allows for the explicit assessment of the impact of wetlands ( e.g ., typology, location) on watershed hydrology. This article is protected by copyright. All rights reserved.