ALBERT

Beschreibung: Abstract Reliable estimation of the volume and timing of snowmelt runoff is vital for water supply and flood forecasting in snow‐dominated regions. Snowmelt is often simulated using temperature‐index (TI) models due to their applicability in data‐sparse environments. Previous research has shown that a modified‐TI model, which uses a radiation‐derived proxy temperature instead of air temperature as its surrogate for available energy, can produce more accurate snow covered area (SCA) maps than a traditional TI model. However, it is unclear whether the improved SCA maps are associated with improved snow water equivalent (SWE) estimation across the watershed or improved snowmelt‐derived streamflow simulation. This paper evaluates whether a modified‐TI model produces better streamflow estimates than a TI model when they are used within a fully‐distributed hydrologic model. It further evaluates the performance of the two models when they are calibrated using either point SWE measurements or SCA maps. The Senator Beck Basin in Colorado is used as the study site because its surface is largely bedrock, which reduces the role of infiltration and emphasizes the role of the SWE pattern on streamflow generation. Streamflow is simulated using both models for six years. The modified‐TI model produces more accurate streamflow estimates (including flow volume and peak flow rate) than the TI model, likely because the modified‐TI model better reproduces the SWE pattern across the watershed. Both models also produce better performance when calibrated with SCA maps instead of point SWE data, likely because the SCA maps better constrain the space‐time pattern of SWE.

Beschreibung: Abstract Groundwater transit time is an essential hydrologic metric for groundwater resources management. However, especially in tropical environments studies on the transit time distribution (TTD) of groundwater infiltration and its corresponding mean transit time (mTT) have been extremely limited due to data sparsity. In this study, we primarily use stable isotopes to examine the TTDs and their mTTs of both vertical and horizontal infiltration at a riverbank infiltration area in the Vietnamese Mekong Delta (VMD), representative of the tropical climate in Asian Monsoon regions. Precipitation, river water, groundwater, and local ponding surface water were sampled for three to nine years and analyzed for stable isotopes (δ18O and δ2H), providing a unique data set of stable isotope records for a tropical region. We quantified the contribution that the two sources contributed to the local shallow groundwater by a novel concept of two‐component lumped parameter models (LPMs) that are solved using δ18O records. The study illustrates that two‐component LPMs, in conjunction with hydrological and isotopic measurements, are able to identify subsurface flow conditions and water mixing at riverbank infiltration systems. However, the predictive skill and the reliability of the models decrease for locations farther from the river, where recharge by precipitation dominates, and a low‐permeable aquitard layer above the highly permeable aquifer is present. This specific setting impairs the identifiability of model parameters. For river infiltration short mTTs (〈40 weeks) were determined for sites closer to the river (〈200 m), whereas for the precipitation infiltration the mTTs were longer (〉80 weeks) and independent of the distance to the river. The results not only enhance the understanding of the groundwater recharge dynamics in the VMD but also suggest that the highly complex mechanisms of surface‐groundwater interaction can be conceptualized by exploiting two‐component LPMs in general. The model concept could thus be a powerful tool for better understanding both the hydrological functioning of mixing processes and the movement of different water components in riverbank infiltration systems.

Beschreibung: Abstract Physically‐based models are useful frameworks for testing intervention strategies designed to reduce elevated sediment loads in agricultural catchments. Evaluating the success of these strategies depends on model accuracy, generally established by a calibration and evaluation process. In this contribution, the physically‐based SHETRAN model was assessed in two similar UK agricultural catchments. The model was calibrated on the Blackwater catchment (18 km2) and evaluated in the adjacent Kit Brook catchment (22 km2) using 4‐years of 15‐minute discharge and suspended sediment flux data. Model sensitivity to changes in single and multiple combinations of parameters as well as sensitivity to changes in Digital Elevation Model (DEM) resolution were assessed. Model flow performance was reasonably accurate; with a Nash‐Sutcliffe efficiency coefficient (NSE) of 0.78 in Blackwater and 0.60 in Kit Brook. In terms of event prediction, the mean of the absolute percentage of difference (μAbsdiff) between measured and simulated flow volume (Qv), peak discharge (Qp), sediment yield (Sy) and peak sediment flux (Sp) showed larger values in Kit Brook (48% [Qv], 66% [Qp], 298% [Sy], 438% [Sp]) compared to the Blackwater catchment (30% [Qv], 41% [Qp], 106% [Sy], 86% [Sp]). Results indicate that SHETRAN can produce reasonable flow prediction but performs less well in estimation of sediment flux, despite reasonably similar hydro‐sedimentary behaviour between catchments. The sensitivity index showed flow volume sensitive to saturated hydraulic conductivity and peak discharge to the Strickler coefficient; sediment yield was sensitive to the overland flow erodibility coefficient and peak sediment flux to raindrop/leaf soil erodibility coefficient. The multi‐parameter sensitivity analysis showed that different combinations of parameters produced similar model responses. Model sensitivity to grid resolution presented similar flow volumes for different DEM resolutions, whereas event peak and duration (for both flow and sediment flux) were highly sensitive to changes in grid size.

Beschreibung: Abstract High‐frequency stable isotope data are useful for validating atmospheric moisture circulation models and provide improved understanding of the mechanisms controlling isotopic compositions in tropical rainfall. Here we present a near‐continuous 6‐month record of O‐ and H‐isotope compositions in both water vapour and daily rainfall from Northeast Australia measured by laser spectroscopy. The data set spans both Wet and Dry Seasons to help address a significant data and knowledge gap in the southern hemisphere tropics. We interpret the isotopic records for water vapour and rainfall in the context of contemporaneous meteorological observations. Surface air moisture provided near‐continuous tracking of the links between isotopic variations and meteorological events on local to regional spatial scales. Power spectrum analysis of the isotopic variation showed a range of significant periodicities, from hourly to monthly scales and cross‐wavelet analysis identified significant regions of common power for hourly‐averaged water vapour isotopic composition and relative humidity, wind direction and solar radiation. Relative humidity had the greatest sub‐diurnal influence on isotopic composition. On longer timescales (weeks to months) isotope variability was strongly correlated with both wind direction and relative humidity. The high‐frequency records showed diurnal isotopic variations in O‐ and H‐isotope compositions due to local dew formation and, for deuterium excess, as a result of evapotranspiration. Several significant negative isotope anomalies on a daily scale were associated with the activity of regional mesoscale convective systems and the occurrence of two tropical cyclones. Calculated air parcel back‐trajectories identified the predominant moisture transport paths from the Southwest Pacific Ocean while moisture transport from northerly directions occurred mainly during the Wet Season monsoonal air flow. Water vapour isotope compositions reflected the same meteorological events as recorded in rainfall isotopes but provided much more detailed and continuous information on atmospheric moisture cycling than the intermittent isotopic record provided by rainfall. Improved global coverage of stable isotope data for atmospheric water vapour is likely to improve simulations of future changes to climate drivers of the hydrological cycle.

Beschreibung: Hydrological Processes, Volume 33, Issue 19, Page 2499-2501, 15 September 2019.

Beschreibung: Abstract Monitoring and estimation of snow depth in alpine catchments is needed for a proper assessment of management alternatives for water supply in these water resources systems. The distribution of snowpack thickness is usually approached by using field data that come from snow samples collected at a given number of locations that constitute the monitoring network. Optimal design of this network is required to obtain the best possible estimates. Assuming that there is an existing monitoring network, its optimization may imply the selection of an optimal network as a subset of the existing one (if there are not funds to maintain them) or enlarging the existing network by one or more stations (optimal augmentation problem). We propose an optimization procedure that minimizes the total variance in the estimate of snowpack thickness. The novelty of this work is to treat, for the first time, the problem of snow observation network optimization for an entire mountain range rather than for small catchments as done in previous studies. Taking into account the reduced data available, which is a common problem in many mountain ranges, the importance of a proper design of these observation networks is even larger. Snowpack thickness is estimated by combining regression models to approach the effect of the explanatory variables and kriging techniques to consider the influence of the stakes location. We solve the optimization problems under different hypotheses, studying the impacts of augmentation and reduction, both, one by one and in pairs. We also analyse the sensitivity of results to non‐snow measurements deduced from satellite information. Finally, we design a new optimal network by combining the reduction and augmentation methods. The methodology has been applied to the Sierra Nevada mountain range (southern Spain), where very limited resources are employed to monitor snowfall and where an optimal snow network design could prove critical. An optimal snow observation network is defined by relocating some observation points. It would reduce the estimation variance by around 600 cm2 (15%).

Beschreibung: Abstract The effects of root systems on soil detachment by overland flow are closely related to vegetation types. The objective of this study was to quantify the effects of two gramineous roots (Paspalum mandiocanum with shallow roots and Pennisetum giganteum with deep roots) on soil detachment capacity, rill erodibility and critical shear stress on alluvial fans of benggang in southeast China. A 4 m long and 0.12 m wide flume was used. Slope steepness ranged from 9% to 27%, and unit flow discharge ranged from 1.39×10‐3 to 4.19×10‐3 m2 s‐1. The mean detachment capacities of Paspalum mandiocanum and Pennisetum giganteum lands were 18% and 38% lower than that of bare land, respectively, and the effects of root on reducing soil detachment were mainly reflected in the 0‐5 cm soil layer. The most important factors in characterizing soil detachment capacity were root length density and soil cohesion, and soil detachment capacity of the two grass lands could be estimated using flow shear stress, soil cohesion, and root length density (NSE=0.90). With the increase in soil depth, rill erodibility increased, while shear stress decreased. The mean rill erodibilities of Paspalum mandiocanum and Pennisetum giganteum lands were 81% and 61% as much as that of bare land, respectively. Additionally, rill erodibilities of the two grass lands could be estimated as an exponential function by root length density and soil cohesion (NSE=0.88). The mean critical shear stress of Paspalum mandiocanum and Pennisetum giganteum lands were 1.29 and 1.39 times that of bare land, respectively, and it could be estimated with a linear function by root length density (NSE=0.76). This study demonstrated that planting of the two grasses Paspalum mandiocanum and Pennisetum giganteum could effectively reduce soil detachment and enhance soil resistance to erosion on alluvial fans, with the deep roots of Pennisetum giganteum being more effective than the shallow roots of Paspalum mandiocanum. The results are helpful for understanding the influencing mechanism of root systems on soil detachment process.

Beschreibung: Abstract The aim of this study was to investigate how the spatial distribution of grass influenced runoff and erosion from a hillslope with loess and cinnamon soils in the rocky area of Northern China. We set up a trial to test the two soils with different treatments, including bare soil (BS), grass strips on the upper (UGS) and lower (DGS) parts of the slope, grass cover over the entire slope (GS), and a grass carpet on the lower part of the slope (GC), under simulated rainfall conditions. The results showed that the runoff coefficients for the loess and cinnamon soils decreased by between 4% and 20% and by between 2% and 37%, respectively, when covered with grass. Grass spatial distribution had little effect on the runoff, but more effect on erosion than vegetation coverage degree. The most effective location of grass cover for decreasing hillslope erosion was at the foot, and the high efficiency was mainly due to controlling of rill formation and sediment deposition. The soil loss from GS, DGS, and GC on the loess and cinnamon soils was between 77% and 93% less and 55% and 80% less, respectively, compared to the loss from BS. However, the soil characteristics had little effect on soil erosion for well‐vegetated slopes. The results highlight the importance of vegetation re‐establishment at the foot of hillslope in controlling soil erosion.

Beschreibung: The manuscript analyses links between flood spatial arrangement and soil water balance in a plain watershed, and, with that purpose, landscape metrics are calculated in maps obtain by remote sensed data in different hydrological scenarios. In contrast with previous works that investigate connectivity, we applied landscape metrics focusing on flooding pattern, their spatial and temporal variability, and their relationship with soil water balance. In addition, the analysis of patterns allows highlighting the internal heterogeneity that plain landscapes usually exhibit. Abstract In areas with very mild relief, water drains in a disordered way due to the lack of a developed drainage network, as it occurs in extremely flat sedimentary regions like the Argentine Pampas. The study analysed the flood spatial arrangements in 2014 by calculating landscape metrics and relating them to soil water balance. The study area is located at Del Azul creek lower basin (Pampa Ecoregion, Argentina). Daily soil water balances were obtained, and seven landscape metrics were calculated in 15 windows in five LandSat images, all along 2014, to explore the relationship between hydrological scenarios and spatial pattern summarized with principal component analysis. Water excess concentrated in winter (June and August); deficits were in late spring and summer (January and November), whereas the beginning of autumn (March) was an intermediate situation. Principal component 1 (44.7%) reflected area and shape metrics and correlated positively with water table level; principal component 2 (32.3%) summarized aggregation ones and was negatively associated with accumulated water excesses or deficits in previous 30 days and useful reserve. Both exhibited possible threshold‐driven behaviour. Internal heterogeneity between NW and SE zones within the study area coincided with the existence of ancient alluvial fans. The results highlight the peculiarities of the flood spatial patterns in regions with very mild relief, where landforms usually determine water flows.

Beschreibung: Abstract This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.

Beschreibung: Scaling issues in snow hydrology persist due to limitations in instrumentation and inability to measure physical properties and processes at spatiotemporal scales required for analysis. Snow depth and water equivalent (SWE) across scale estimated using time‐lapse photos, transects, and model grids (Canadian Meteorological Centre depth, GlobSnow SWE) were found to represent different physical processes and have substantially different statistical moments. Findings have implications for understanding limitations of distributing snowpack measurements, data assimilation, and validation of remotely sensed estimates. Abstract This study investigates scaling issues by evaluating snow processes and quantifying bias in snowpack properties across scale in a northern Great Lakes–St. Lawrence forest. Snow depth and density were measured along transects stratified by land cover over the 2015/2016 and 2016/2017 winters. Daily snow depth was measured using a time‐lapse (TL) camera at each transect. Semivariogram analysis of the transect data was conducted, and no autocorrelation was found, indicating little spatial structure along the transects. Pairwise differences in snow depth and snow water equivalent (SWE) between land covers were calculated and compared across scales. Differences in snowpack between forested sites at the TL points corresponded to differences in canopy cover, but this relationship was not evident at the transect scale, indicating a difference in observed process across scale. TL and transect estimates had substantial bias, but consistency in error was observed, which indicates that scaling coefficients may be derived to improve point scale estimates. TL and transect measurements were upscaled to estimate grid scale means. Upscaled estimates were compared and found to be consistent, indicating that appropriately stratified point scale measurements can be used to approximate a grid scale mean when transect data are not available. These findings are important in remote regions such as the study area, where frequent transect data may be difficult to obtain. TL, transect, and upscaled means were compared with modelled depth and SWE. Model comparisons with TL and transect data indicated that bias was dependent on land cover, measurement scale, and seasonality. Modelled means compared well with upscaled estimates, but model SWE was underestimated during spring melt. These findings highlight the importance of understanding the spatial representativeness of in situ measurements and the processes those measurements represent when validating gridded snow products or assimilating data into models.

Beschreibung: Abstract Forest canopies present irregular surfaces that alter both the quantity and spatiotemporal variability of precipitation inputs. The drop size distribution (DSD) of rainfall varies with rainfall event characteristics and is altered substantially by the forest stand properties. Yet, the influence of two major European tree species, European beech (Fagus sylvatica L.) and Norway spruce (Picea abies (L.) H. Karst), on throughfall DSD is largely unknown. In order to assess the impact of these two species with differing canopy structures on throughfall DSD, two optical disdrometers, one above and one below the canopy of each European beech and Norway spruce, measured DSD of both incident rainfall and throughfall over two months at a 10 second resolution. Fractions of different throughfall categories were analysed for single precipitation events of different intensities. While penetrating the canopies, clear shifts in drop size and temporal distributions of incoming rainfall were observed. Beech and spruce, however, had different DSD, behaved differently in their effect on diameter volume percentiles as well as width of drop spectrum. The maximum drop sizes under beech were higher than under spruce. The mean ± standard deviation of the median volume drops size (D50) over all rain events was 2.7 ± 0.28 mm for beech and 0.80 ± 0.04 mm for spruce, respectively. In general, there was a high DSD variability within events indicating varying amounts of the different throughfall fractions. These findings help to better understand the effects of different tree species on rainfall partitioning processes and small scale variations in subcanopy rainfall inputs, thereby demonstrating the need for further research in high resolution spatial and temporal properties of rainfall and throughfall.

Beschreibung: A recently theoretically deduced rill flow resistance equation, based on a power‐velocity profile, was tested experimentally on plots of varying slopes (ranging from 9% to 26%) in which mobile and fixed bed rills were incised. The measurements carried out in both conditions confirmed that the Darcy–Weisbach friction factor can be accurately estimated using the proposed theoretical approach. For the investigated conditions, the effect of sediment transport on the flow resistance law is negligible respect to the grain roughness effect. Abstract Rills caused by run‐off concentration on erodible hillslopes have very irregular profiles and cross‐section shapes. Rill erosion directly depends on the hydraulics of flow in the rills, which may differ greatly from hydraulics of flow in larger and regular channels. In this paper, a recently theoretically deduced rill flow resistance equation, based on a power–velocity profile, was tested experimentally on plots of varying slopes (ranging from 9% to 26%) in which mobile and fixed bed rills were incised. Initially, measurements of flow velocity, water depth, cross‐section area, wetted perimeter, and bed slope, carried out in 320 reaches of mobile bed rills and in 165 reaches of fixed rills, were used for calibrating the theoretical flow resistance equation. Then the relationship between the velocity profile parameter Γ, the channel slope, and the flow Froude number was separately calibrated for the mobile bed rills and for the fixed ones. The measurements carried out in both conditions (fixed and mobile bed rills) confirmed that the Darcy–Weisbach friction factor can be accurately estimated using the proposed theoretical approach. For mobile bed rills, the data were supportive of the slope independence hypothesis of velocity, due to the feedback mechanism, stated by Govers. The feedback mechanism was able to produce quasicritical flow conditions. For fixed bed rills, obtained by fixing the rill channel, by a glue, at the end of the experimental run with a mobile bed rill, the slope independence of the flow velocity measurements was also detected. Therefore, an experimental run carried out by a rill bed fixed after modelling flow action is useful to detect the feedback mechanism. Finally, the analysis showed that, for the investigated conditions, the effect of sediment transport on the flow resistance law can be considered negligible respect to the grain roughness effect.

Beschreibung: Analyses of stem water stable isotope composition were used to identify water taken up by plants. The Populus trees and Shepherdia and Symphoricarpos shrubs had contrasting functional rooting depths as illustrated by the different seasonal patterns of change in deuterium excess values, which indicated greater shallow soil water use by the shrub species. Abstract Riparian cottonwood forests in dry regions of western North America do not typically receive sufficient growing season precipitation to completely support their relatively high transpiration requirements. Water used in transpiration by riparian ecosystems must include alluvial groundwater or water stored in the potentially large reservoir of the unsaturated soil zone. We used the stable oxygen and hydrogen isotope composition of stem xylem water to evaluate water sources used by the dominant riparian cottonwood (Populus spp.) trees and shrubs (Shepherdia argentea and Symphoricarpos occidentalis) in Lethbridge, Alberta, during 3 years of contrasting environmental conditions. Cottonwoods did not exclusively take up alluvial groundwater but made extensive use of water sourced from the unsaturated soil zone. The oxygen and hydrogen isotope compositions of cottonwood stem water did not strongly overlap with those of alluvial groundwater, which were closely associated with the local meteoric water line. Instead, cottonwood stem water δ18O and δ2H values were located below the local meteoric water line, forming a line with a low slope that was indicative of water exposed to evaporative enrichment of heavy isotopes. In addition, cottonwood xylem water isotope compositions had negative values of deuterium excess (d‐excess) and line‐conditioned (deuterium) excess (lc‐excess), both of which provided evidence that water taken up by the cottonwoods had been exposed to fractionation during evaporation. The shrub species had lower values of d‐excess and lc‐excess than had the cottonwood trees due to shallower rooting depths, and the d‐excess values declined during the growing season, as shallow soil water that was taken up by the plants was exposed to increasing, cumulative evaporative enrichment. The apparent differences in functional rooting pattern between cottonwoods and the shrub species, strongly influenced the ratio of net photosynthesis to stomatal conductance (intrinsic water‐use efficiency), as shown by variation among species in the δ13C values of leaf tissue.

Beschreibung: Hydrological Processes, Volume 33, Issue 17, Page 2263-2265, 15 August 2019.

Beschreibung: Hydrological Processes, Volume 0, Issue ja, -Not available-.

Beschreibung: Abstract Meltwater from glaciers is not only a stable source of water but also affects downstream streamflow dynamics. One of these dynamics is the interannual variability of streamflow. Glaciers can moderate streamflow variability, because the runoff in the glacierised part, driven by temperature, correlates negatively with the runoff in the non‐glacierised part of a catchment, driven by precipitation, thereby counterbalancing each other. This is also called the glacier compensation effect (GCE) and the effect is assumed to depend on relative glacier cover. Previous studies found a convex relationship between streamflow variability and glacier cover of different glacierised catchments, with lowest streamflow variability at a certain optimum glacier cover. In this study we aim to revisit these previously found curves to find out if a universal relationship between interannual streamflow variability and glacier cover exists, which could potentially be used in a space‐for‐time substitution analysis. Moreover, we test the hypothesis that the dominant climate drivers (here precipitation and temperature) switch around the suggested optimum of the curve. First a set of virtual nested catchments, with the same absolute glacier area but varying non‐glacierised area, were modelled to isolate the effect of glacier cover on streamflow variability. The modelled relationship was then compared to a multi‐catchment dataset of gauged glacierised catchments in the European Alps. In a third step, changes of the GCE curve over time were analysed. Model results showed a convex relationship and the optimum in the simulated curve aligned with a switch in the dominant climate driver. However, the multi‐catchment data and the time change analyses did not suggest the existence of a universal convex relationship. Overall, we conclude that GCE is complex due to entangled controls and changes over time in glacierised catchments. Therefore, care should be taken to use a GCE curve for estimating and/or predicting interannual streamflow variability in glacierised catchments.

Beschreibung: Abstract Base flows are important for tropical regions with pronounced dry seasons which are facing increasing water demands. Base flow generation, however, is one of the most challenging hydrological processes to characterize in the tropics. In many years during the May‐December wet season in the Panama Canal Watershed (PCW), base flows in rivers abruptly increase. This increase persists until the start of the December‐April dry season. Understanding this unusual base flow jump (BFJ) behavior is critical to improve water provisioning in seasonal tropics, especially during droughts and extended dry seasons. This study developed an integrated approach combining piecewise regression on cumulative average base flow and sensitivity analysis to calculate the timing and magnitude of BFJ. Rainfall, forest cover, mean land surface slope, catchment area and estimated subsurface storage were tested as predictors for the occurrence and magnitude of the BFJs in seven sub‐catchments of the PCW. Sensitivity analysis on correlated predictors allowed ranking of predictor contributions due to isolated and cross correlation effects. Correlations between observed BFJs and BFJs predicted by watershed and rainfall‐related predictors were 0.92 and 0.65 for BFJ timing and magnitude, respectively. Forest cover was the second most significant predictor after cumulative rainfall for jump magnitude, owing to larger subsurface storage and groundwater recharge in forests than pastures. Catchments in the mountainous eastern PCW always generated larger jumps due to their higher rainfall and greater forest cover than the western PCW catchments. The cross‐correlations between predictors contributed to more than 50% of the jump variances. The results demonstrate the importance of rainfall gradient and catchment characteristics in affecting the sudden and sustained BFJs, which can help inform land management decisions intended to enhance water supplies in tropics. This study underscores the need for more research to further understand the hydrological processes involved in the BFJ phenomenon, including better BFJ models and field characterizations, to help improve tropical ecosystem services under a changing environment.

Beschreibung: Melting seasonal ground ice reduces potential evapotranspiration and may be a mechanism for peatland persistence in the Western Boreal Plain Abstract Peatlands in the Western Boreal Plains act as important water sources in the landscape. Their persistence, despite potential evapotranspiration (PET) often exceeding annual precipitation, is attributed to various water storage mechanisms. One storage element that has been understudied is seasonal ground ice (SGI). This study characterized spring SGI conditions and explored its impacts on available energy, actual evapotranspiration, water table, and near surface soil moisture in a western boreal plains peatland. The majority of SGI melt took place over May 2017. Microtopography had limited impact on melt rates due to wet conditions. SGI melt released 139mm in ice water equivalent (IWE) within the top 30cm of the peat, and weak significant relationships with water table and surface moisture suggest that SGI could be important for maintaining vegetation transpiration during dry springs. Melting SGI decreased available energy causing small reductions in PET (〈10mm over the melt period) and appeared to reduce actual evapotranspiration variability but not mean rates, likely due to slow melt rates. This suggests that melting SGI supplies water, allowing evapotranspiration to occur at near potential rates, but reduces the overall rate at which evapotranspiration could occur (PET). The role of SGI may help peatlands in headwater catchments act as a conveyor of water to downstream landscapes during the spring while acting as a supply of water for the peatland. Future work should investigate SGI influences on evapotranspiration under differing peatland types, wet and dry spring conditions, and if the spatial variability of SGI melt leads to spatial variability in evapotranspiration.

Beschreibung: Abstract In urban areas, the presence of impervious surfaces limits natural drainage and routes water to stormwater infrastructure with finite capacity, making these areas especially prone to flooding. Though large floods are responsible for endangering lives and causing extensive damage, there is growing evidence that more frequent floods with shallow water depths, termed nuisance flooding, can have a high cumulative cost and many direct and indirect damages. To determine whether locations of nuisance flooding may be linked to topography, we took a parsimonious, spatially distributed approach to explore whether high topographic index values co‐occur with citizen‐reported nuisance flooding. We obtained nuisance flood reports from the municipal data service 311 for several watersheds in New York City and Baltimore, USA. Our analysis tested two topographic indices (TI)—topographic wetness index (TWI) and sink depth—both calculated from high‐resolution (~1 m) digital elevation models. Generally, our findings suggest that not all but many locations of reported flooding tend to coincide with deep sinks or large TWI. However, nuisance flooding reports most commonly coincided with deep sinks and high TWI when using a maximum, instead of coincident, TI value extracted around each reported location of flooding, an approach we used due to the uncertainty in location accuracy of flooding reports. Overall, our results show promise for application of topographic indices, typically applied in more natural settings, as indicators of nuisance flooding areas in urbanized environments. Although limitations to this approach exist, the application of TIs and crowd‐sourced reporting in tandem could provide a useful starting point for mapping flood‐prone areas in many cities with technologically adept community members and ample geospatial data.

Beschreibung: 1 There was a poor correlation between the flow intensity indices of velocity category and nondimensional sediment transport capacity.2 Ignoring the partition phenomenon of the relation curves, stream power can be used to predict sediment transport capacity, with a coefficient of determination of 0.85.3 An empirical formula for predicting sediment transport with a coefficient of determination of 0.90 was established by multiple regression analysis based on the general flow intensity index. Abstract Soil erosion is a major contributor to land degradation in the Loess Plateau in China. To clarify the sediment transport capacity of overland flow influenced by hydraulic parameters, such as shear stress, sand shear stress (hydraulic gradient partition method and hydraulic radius partition method), mean flow velocity, Froude number, stream power, and unit stream power, indoor experiments with eight‐unit‐width flow discharges from 0.0667 × 10−3 to 0.3333 × 10−3 m2·s−1, six slope gradients from 3.49 to 20.79%, and two kinds of sand soils (d50 = 0.17 and 0.53 mm) were systematically investigated. A nondimensional method was adopted in data processing. Results showed that there was a partition phenomenon of relation curves because of the different median grain diameters. The correlation between the nondimensional stream power and nondimensional sediment transport capacity was the highest, followed by the correlation between the nondimensional unit stream power and nondimensional sediment transport capacity. However, there was a poor correlation between the flow intensity indices of velocity category and nondimensional sediment transport capacity. Nondimensional stream power, nondimensional unit stream power, and nondimensional shear stress could predict sediment transport capacity well. Ignoring the partition phenomenon of the relation curves, stream power could be used to predict sediment transport capacity, with a coefficient of determination of .85. Furthermore, a general flow intensity index was obtained to predict sediment transport capacity of overland flow. Finally, an empirical formula for predicting sediment transport capacity with a coefficient of determination of .90 was established by multiple regression analyses based on the general flow intensity index. During the analysis between measured sediment transport capacities in present study and predicted values based on Zhang model, Mahmoodabadi model, and Wu model, it was found that these three models could not accurately predict sediment transport capacities of this study because different models are estimated on the basis of different experimental conditions.

Beschreibung: Abstract A few relatively large reservoirs, hundreds of small reservoirs and numerous farm dams were built in the upper Gan River Basin, China. The operation of such a reservoir network can serve as a significant source of variability in the local hydrological regime and should be included in research to better understand the interaction between multiple hydrological processes and watershed management. In this study, a reservoir network module that included reservoirs of multiple sizes was developed and fully integrated into a coupled land surface and distributed hydrologic model, CLHMS, for a detailed description of the hydrological impact of a reservoir network. A generalized release scheme was employed to determine the outflow of both large and small reservoirs. The integrated model was then evaluated against observations and reanalysis data, which indicate that the model can reasonably reconstruct the reservoir operation, streamflow and other hydrological variables. Results quantitatively demonstrate that a reservoir network can result in an increased streamflow in dry seasons, a decreased streamflow in wet seasons, a generally larger groundwater discharge, higher groundwater level, a slightly damper soil condition and a larger amount of evapotranspiration at the basin level. With the integrated model, it is feasible to achieve more sustainable watershed planning and management.

Beschreibung: Abstract Seasonal snowpacks in marginal snow environments are typically warm and nearly isothermal, exhibiting high inter‐ and intra‐annual variability. Measurements of snow depth and snow water equivalent were made across a small subalpine catchment in the Australian Alps over two snow seasons in order to investigate the extent and implications of snowpack spatial variability in this marginal setting. The distribution and dynamics of the snowpack were found to be influenced by upwind terrain, vegetation, solar radiation and slope. The role of upwind vegetation was quantified using a novel parameter based on gridded vegetation height. The elevation range of the catchment was relatively modest (185 m), and elevation impacted distribution but not dynamics. Two characteristic features of marginal snowpack behaviour are presented. Firstly, the evolution of the snowpack is described in terms of a relatively unstable accumulation state and a highly stable ablation state, as revealed by temporal variations in the mean and standard deviation of snow water equivalent. Secondly, the validity of partitioning the snow season into distinct accumulation and ablation phases is shown to be compromised in such a setting. Snow at the most marginal locations may undergo complete melt several times during a season and, even where snow cover is more persistent, ablation processes begin to have an effect on the distribution of the snowpack early in the season. Our results are consistent with previous research showing that individual point measurements are unable to fully represent the variability in the snowpack across a catchment, and we show that recognising and addressing this variability is particularly important for studies in marginal snow environments.

Beschreibung: Abstract Spatiotemporal heterogeneity in soil water content is recognized as a common phenomenon, but heterogeneity in the hydrogen and oxygen isotope composition of soil water, which can reveal processes of water cycling within soils, has not been well studied. New advances are being driven by measurement approaches allowing sampling with high density in both space and time. Using in situ soil water vapor probe techniques, combined with conventional soil and plant water vacuum distillation extraction, we monitored the hydrogen and oxygen stable isotopic composition of soil and plant waters at paired sites dominated by grasses and Gambel's oak (Quercus gambelli) within a semi‐arid montane ecosystem over the course of a growing season. We found that sites spaced only 20 m apart had profoundly different soil water isotopic and volumetric conditions. We document patterns of depth‐ and time‐explicit variation in soil water isotopic conditions at these sites, and consider mechanisms for the observed heterogeneity. We found that soil water content and isotopic variability was damped under Quercus gambelli, perhaps due in part to hydraulic redistribution of deep soil water or groundwater by Quercus gambelli in these soils relative to the grass‐dominated site. We also found some support for H isotope discrimination effects during water uptake by Quercus gambelli. In this ecosystem, the soil water content was higher than that at the neighboring grass site, and thus 25% more water was available for transpiration by Quercus gambelli compared to the grass site. This work highlights the role of plants in governing soil water variation and demonstrates that they can also strongly influence the isotope ratios of soil water. The resulting fine‐scale heterogeneity has implications for the use of isotope tracers to study soil hydrology and evaporation and transpiration fluxes to improve understanding of water cycling through the soil‐plant‐atmosphere continuum.

Beschreibung: Abstract While we know that rainfall interception (the rain caught, stored, and evaporated from aboveground vegetative surfaces and ground litter) is affected by rain and throughfall drop size, what was unknown until now is the relative proportion of each throughfall type (free throughfall, splash throughfall, canopy drip) beneath coniferous and broadleaved trees. Based on a multi‐national dataset of 〉 120 million throughfall drops, we found that the type, number, and volume of throughfall drops are different between coniferous and broadleaved tree species, leaf states, and timing within rain events. Compared to leafed broadleaved trees, conifers had a lower percentage of canopy drip (51% vs. 69% with respect to total throughfall volume) and slightly smaller diameter splash throughfall and canopy drip. Canopy drip from leafless broadleaved trees consisted of fewer and smaller diameter drops (D50_DR, fiftieth cumulative drop volume percentile for canopy drip, of 2.24 mm) than leafed broadleaved trees (D50_DR of 4.32 mm). Canopy drip was much larger in diameter under woody drip points (D50_DR of 5.92 mm) than leafed broadleaved trees. Based on throughfall volume, the percentage of canopy drip was significantly different between conifers, leafed broadleaved trees, leafless broadleaved trees, and woody surface drip points (p ranged from 〈 0.001 to 0.005). These findings are partly attributable to differences in canopy structure and plant surface characteristics between plant functional types and canopy state (leaf, leafless), among other factors. Hence, our results demonstrating the importance of drop‐size dependent partitioning between coniferous and broadleaved tree species could be useful to those requiring more detailed information on throughfall fluxes to the forest floor.

Beschreibung: Abstract High‐elevation tropical grassland systems, called Páramo, provide essential ecosystem services such as water storage and supply for surrounding and lowland areas. Páramo systems are threatened by climate and land use changes. Rainfall generation processes and moisture transport pathways influencing precipitation in the Páramo are poorly understood but needed to estimate the impact of these changes, particularly during El Niño conditions which largely affect hydrometeorological conditions in tropical regions. To fill this knowledge gap, we present a stable isotope analysis of rainfall samples collected on a daily to weekly basis between January 2015 and May 2016 during the strongest El Niño event recorded in history (2014‐2016) in two Páramo regions of Central America (Chirripó, Costa Rica) and the northern Andes (Cajas, south Ecuador). Isotopic compositions were used to identify how rainfall generation processes (convective and orographic) change seasonally at each study site. HYSPLIT air mass back trajectory analysis was used to identify preferential moisture transport pathways to each Páramo site. Our results show the strong influence of northeast trade winds to transport moisture from the Caribbean Sea to Chirripó and the South American low‐level jet to transport moisture from the Amazon forest to Cajas. These moisture contributions were also related to the formation of convective rainfall associated with the passage of the Intertropical Convergence Zone over Costa Rica and Ecuador during the wetter seasons and to orographic precipitation during the transition and drier seasons. Our findings provide essential baseline information for further research applications of water stable isotopes as tracers of rainfall generation processes and transport in the Páramo and other montane ecosystems in the tropics.

Beschreibung: Abstract Seasonal snow cover in mountainous regions will affect local climate and hydrology. In this study, we assessed the role of altitude in determining the relative importance of temperature and precipitation in snow cover variability in the Central Tianshan Mountains. The results show that: (1) in the study area, temperature has a greater influence on snow cover than precipitation during most of the time period studied and in most altitudes. (2) In the high‐elevation area, there is a threshold altitude of 3900±400 m, below which temperature is negatively while precipitation is positively correlated to snow cover, above which the situation is the opposite. Besides, this threshold altitude decreases from snow accumulated period to snow stable period and then increases from snowmelt period to snow‐free period. (3) Below 2000 m, there is another threshold altitude of 1400±100 m during the snow stable period, below (above) which precipitation (temperature) is the main driver of snow cover.

Beschreibung: Hydrological Processes, Volume 33, Issue 6, Page 889-891, 15 March 2019.

Beschreibung: Abstract Low flow events can cause significant impacts to river ecosystems and water‐use sectors; as such it is important to understand their variability and drivers. In this study, we characterize the variability and timing of annual total frequency of low streamflow days across a range of headwater streams within the continental United States (US). To quantify this, we use a metric that counts the annual number of low flow days below a given threshold, defined as the Cumulative Dry days Occurrence (CDO). First, we identify three large clusters of streamgauge locations using a Partitioning Around Medoids (PAM) clustering algorithm. In terms of timing, results reveal that for most clusters, the majority of low streamflow days occur from the middle of summer until early fall, though several locations in Central and Western US also experience low flow days in cold seasons. Further, we aim to identify the regional climate and larger‐scale drivers for these low streamflow days. Regionally, we find that precipitation deficits largely associate with low streamflow days in the western US, while within the central and eastern US clusters, high temperature indicators are also linked to low streamflow days. In terms of larger‐scale, we examine sea surface temperature (SST) anomalies, finding that extreme dry years exhibit a high degree of co‐occurrence with different patterns of warmer SST anomalies across the Pacific and northern Atlantic Oceans. The linkages identified with regional climate and SSTs offer promise towards regional prediction of changing conditions of low streamflow events.

Beschreibung: Abstract Plant transpiration depends on environmental conditions, and soil water availability is its primary control under water deficit conditions. In this study, we improve a simplified process‐based model proposed by Buckley, Turnbull & Adams (2012) (hereafter "BTA") by including soil water potential (ψsoil) to explicitly represent the dependence of plant transpiration on root‐zone moisture conditions. The improved model is denoted as the BTA‐ψ model. We assessed the performance of the BTA and BTA‐ψ models in a subtropical monsoon climate and a Mediterranean climate with different levels of water stress. The BTA model performed reasonably in estimating daily and hourly transpiration under sufficient water conditions, but it failed during dry periods. Overall, the BTA‐ψ model provided a significant improvement for estimating transpiration under a wide range of soil moisture conditions. Although both models could estimate transpiration (sap flow) at night, BTA‐ψ was superior to BTA in this regard. Species differences in the calibrated parameters of both models were consistent with leaf‐level photosynthetic measurements on each species, as expected given the physiological basis of these parameters. By combining a simplified representation of physiological regulation with reasonable performance across a range of soil moisture conditions, the BTA‐ψ model provides a useful alternative to strictly empirical models for modeling transpiration.

Beschreibung: Abstract Characterization of spatial and temporal variability of stable isotopes (δ18O & δ2H) of surface waters is essential to interpret hydrological processes and establish modern isotope‐elevation gradients across mountainous terrains. Here, we present stable isotope data for river waters across Kyrgyzstan. River water isotopes exhibit substantial spatial heterogeneity among different watersheds in Kyrgyzstan. Higher river water isotope values were found mainly in the Issyk–Kul Lake watershed, whereas waters in the Son–Kul Lake watershed display lower values. Results show a close δ18O–δ2H relation between river water and the LMWL (Local Meteoric Water Line), implying that river water experiences little evaporative enrichment. River water from the high elevation regions (e.g., Naryn and Son–Kul Lake watershed) had the most negative isotope values, implying that river water is dominated by snowmelt. Higher d‐excess (average d=13.9‰) in river water probably represents the isotopic signature of combined contributions from direct precipitation and glacier melt in stream discharge across Kyrgyzstan. A significant relationship between river water δ18O and elevation was observed with a vertical lapse rate of 0.13 ‰/100 m. These findings provide crucial information about hydrological processes across Kyrgyzstan and contribute to a better understanding of the paleoclimate/elevation reconstruction of this region.

Beschreibung: Abstract Evaporative flux is a key component of hydrological budgets. Water loss through evapotranspiration reduces volumes available for runoff. The transition from liquid to water vapour on open water surfaces requires heat. Consequently, evaporation act as a cooling mechanism during summer. Both river discharge and water temperature simulations are thus influenced by the methods used to model evaporation. In this paper, the impact of evapotranspiration estimation methods on simulated discharge is assessed using a semi‐distributed model on two Canadian watersheds. The impact of evaporation estimation methods on water temperature simulations is also evaluated. Finally, the validity of using the same formulation to simulate both of these processes is verified. Five well known evapotranspiration models and five evaporation models with different wind functions were tested. Results show a large disparity (18‐22% of mean annual total evapotranspiration) among the evapotranspiration methods, leading to important differences in simulated discharge (3‐25% of observed discharge). Larger differences results from evaporation estimation methods with mean annual divergences of 34‐48%. This translates into a difference in mean summer water temperature of 1‐15%. Results also show that the choice of model parameter has less influence than the choice of evapotranspiration method in discharge simulations. However, the parameter values influence thermal simulations in the same order of magnitude as the choice of evaporation estimation method. Overall, the results of this study suggest that evapotranspiration and open water evaporation should be represented separately in a hydrological modelling framework, especially when water temperature simulations are required.

Beschreibung: Abstract Dissolved organic carbon (DOC) originating in peatlands can be mineralized to carbon dioxide (CO2) and methane (CH4), two potent greenhouse gases. Knowledge of the dynamics of DOC export via runoff is needed for a more robust quantification of C cycling in peatland ecosystems, a prerequisite for realistic predictions of future climate change. We studied dispersion pathways of DOC in a mountain‐top peat bog in the Czech Republic (Central Europe), using a dual isotope approach. While δ13CDOC values made it possible to link exported DOC with its within‐bog source, δ18OH2O values of precipitation and runoff helped to understand runoff generation. Our two‐year DOC‐H2O isotope monitoring was complemented by a laboratory peat incubation study generating an experimental time‐series of δ13CDOC values. DOC concentrations in runoff during high‐flow periods were 20‐30 mg L‐1. The top 2 cm of the peat profile, composed of decaying green moss, contained isotopically lighter C than deeper peat, and this isotopically light C was present in runoff in high‐flow periods. In contrast, baseflow contained only 2‐10 mg DOC L‐1, and its more variable C isotope composition intermittently fingerprinted deeper peat. DOC in runoff occasionally contained isotopically extremely light C whose source in solid peat substrate was not identified. Pre‐event water made up on average 60 % of the water runoff flux, while direct precipitation contributed 40 %. Runoff response to precipitation was relatively fast. A highly leached horizon was identified in shallow catotelm. This peat layer was likely affected by a lateral influx of precipitation. Within 36 days of laboratory incubation, isotopically heavy DOC that had been initially released from the peat was replaced by isotopically lighter DOC, whose δ13C values converged to the solid substrate and natural runoff. We suggest that δ13C systematics can be useful in identification of vertically stratified within‐bog DOC sources for peatland runoff.

Beschreibung: Abstract Reservoir tillage (RT) improves the soil rainwater harvesting capacity and reduces soil erosion on cropland, but there is some debate regarding its effectiveness. The objective of this study was to further verify the effect of RT on soil erosion and explore the reasons for this effect by analysing microrelief changes during rainfall. Rainfall intensities of 60, 90 and 120 mm/h and three slope degrees (5°, 15° and 25°, representing gentle, medium and steep slopes) were considered. A smooth surface (SS) served as the control. The microrelief changes were determined based on digital elevation models (DEMs), which were measured using a laser scanner with a 2‐cm grid before and after rainfall events. The results showed that compared to the values for the SS, RT reduced both the runoff and sediment by approximately 10‐20% on the gentle slope; on the medium slope, although RT also reduced the runoff in the 90 and 120 mm/h intensity rainfall events, the sediment increased by 158.90% and 246.08%; on the steep slope, the sediment increased by 92.33 to 296.47%. Overall, when the runoff control benefit of RT was lower than 5%, there was no sediment control benefit. RT was effective at controlling soil loss on the gentle slopes but was not effective on the medium and steep slopes. This is because the surface depressions created by RT were filled in with sediment that eroded from the upslopes, and the surface microrelief became smoother, which then caused greater soil and water loss than that on an SS at the later rainfall stage.

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

Beschreibung: Abstract Catchments consist of distinct landforms that affect the storage and release of subsurface water. Certain landforms may be the main contributors to streamflow during extended dry periods and these may vary for different catchments in a given region. We present a unique dataset from snapshot field campaigns during low‐flow conditions in eleven catchments across Switzerland to illustrate this. The catchments differed in size (10 to 110 km2), varied from predominantly agricultural lowlands to Alpine areas, and covered a range of physical characteristics. During each snapshot campaign, we jointly measured streamflow and collected water samples for the analysis of major ions and stable water isotopes. For every sampling location (basin), we determined several landscape characteristics from national geo‐datasets, including drainage area, elevation, slope, flowpath length, dominant land use, and geological and geomorphological characteristics, such as the lithology and fraction of Quaternary deposits. The results demonstrate very large spatial variability in specific low‐flow discharge and water chemistry: neighboring sampling locations could differ significantly in their specific discharge, isotopic composition and ion concentrations, indicating that different sources contribute to streamflow during extended dry periods. However, none of the landscape characteristics that we analyzed could explain the spatial variability in specific discharge or stream water chemistry in multiple catchments. This suggests that local features determine the spatial differences in discharge and water chemistry during low‐flow conditions and that this variability cannot be assessed a priori from available geodata and statistical relations to landscape characteristics. The results furthermore suggest that measurements at the catchment outlet during low‐flow conditions do not reflect the heterogeneity of the different source areas in the catchment that contribute to streamflow.

Beschreibung: Abstract In headwater catchments, streamflow recedes between periods of rainfall at a predictable rate generally defined by a power‐law relationship relating streamflow decay to streamflow. Research over the last four decades has applied this relationship to predictions of water resource availability as well as estimations of basin‐wide physiographic characteristics and ecohydrologic conditions. However, the interaction of biophysical processes giving rise to the form of these power‐law relationships remain poorly understood, and recent investigations into the variability of streamflow recession characteristics between discrete events have alternatively suggested evapotranspiration, water table elevation, and stream network contraction as dominant factors, without consensus. To assess potential temporal variability and interactions in the mechanism(s) driving streamflow recession, we combine long‐term observational data from a headwater stream in the southern Appalachian Mountains with state and flux conditions from a process‐based ecohydrologic model. Streamflow recession characteristics are non‐unique, and vary systematically with seasonal fluctuations in both rates of transpiration and watershed wetness conditions, such that transpiration dominates recession signals in the early growing season and diminishes in effect as the water table elevation progressively drops below and decouples with the root zone with topographic position. As a result of this decoupling, there exists a seasonal hysteretic relationship between streamflow decay and both evapotranspiration and watershed wetness conditions. Results indicate that for portions of the year, forest transpiration may actively compete with subsurface drainage for the same water resource that supplies streamflow, though for extended time periods these processes exploit distinct water stores. Our analysis raises concerns about the efficacy of assessing humid headwater systems using traditional recession analysis, with recession curve parameters treated as static features of the watershed, and we provide novel alternatives for evaluating interacting biological and geophysical drivers of streamflow recession.

Beschreibung: Abstract The study on the hydraulic properties of coastal aquifers has significant implications both in hydrological sciences and environmental engineering. Although many analytical solutions are available, most of them are based on the same basic assumption that assumes aquifers extend landward semi‐infinitely which does not necessarily reflect the reality. In this study, the general solutions for a leaky confined coastal aquifer have been developed that consider both finitely landward constant‐head and no‐flow boundaries. The newly developed solutions were then used to examine theoretically the joint effects of leakage and aquifer length on hydraulic head fluctuations within the leaky confined aquifer, and the validity of using the simplified solution that assumes the aquifer is semi‐infinite. The results illustrated that the use of the simplified solution may cause significant errors, depending on joint effects of leakage and aquifer length. A dimensionless characteristic parameter was then proposed as an index for judging the applicability of the simplified solution. In addition, practical application of the general solution for the constant‐head inland boundary was used to characterize the hydraulic properties of a leaky confined aquifer using the data collected from a field site at the Seine River estuary, France, and the versatility of the general solution was further justified.

Beschreibung: Abstract The stable isotopes of water (δ2H and δ18O) are useful conservative tracers for tracking the movement of water in soil. But while the tracking of water infiltrating through the soil profile and its movement as runoff and groundwater recharge are well developed, water movement through the soil can also include evaporative fractionation. Soil water fractionation factors have, until now, been largely empirical. Unlike open water evaporation where temperature, humidity and vapor pressure gradient define fractionation, soil water evaporation includes fractionation by soil matrix effects. These effects are still poorly characterized. Here we present preliminary results from a simple laboratory experiment with four soil admixtures with grain sizes that range from sand to silt and clay. Our results show that soil tension seems to control the isotope fractionation of resident soil water. The relationship between soil tension and equilibrium fractionation appears to be independent of soil texture and appears well supported by thermodynamic theory. While these results are preliminary, they suggest that future work should go after soil tension effects as a possible explanatory factor of soil water and water vapor fractionation.

Beschreibung: Abstract Surface water flooding (SWF) is a recurrent hazard that affects lives and livelihoods. Climate change is projected to change the frequency of extreme rainfall events that can lead to SWF. Increasingly, data from Regional Climate Models (RCMs) are being used to investigate the potential water‐related impacts of climate change; such assessments often focus on broad‐scale fluvial flooding and the use of coarse resolution (〉12km) RCMs. However, high‐resolution (〈4km) convection‐permitting RCMs are now becoming available that allow impact assessments of more localised SWF to be made. At the same time, there has been an increasing demand for more robust and timely real‐time forecast and alert information on SWF. In the UK, a real‐time SWF Hazard Impact Model framework has been developed. The system uses 1km gridded surface runoff estimates from a hydrological model to simulate the SWF hazard. These are linked to detailed inundation model outputs through an Impact Library to assess impacts on property, people, transport and infrastructure for four severity levels. Here, a set of high‐resolution (1.5km and 12km) RCM data has been used as input to a grid‐based hydrological model over southern Britain to simulate Current (1996‐2009) and Future (~2100s; RCP8.5) surface runoff. Counts of threshold‐exceedance for surface runoff and precipitation (at 1‐, 3‐ and 6‐hour durations) are analysed. Results show that the percentage increases in surface runoff extremes, are less than those of precipitation extremes. The higher‐resolution RCM simulates the largest percentage increases, which occur in winter, and the winter exceedance counts are greater than summer exceedance counts. For property impacts the largest percentage increases are also in winter however, it is the 12km RCM output that leads to the largest percentage increase in impacts. The added‐value of high‐resolution climate model data for hydrological modelling is from capturing the more intense convective storms in surface runoff estimates.

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

Beschreibung: Abstract The paper presents oxygen and hydrogen isotopes of 284 precipitation event samples systematically collected in Irkutsk, in the Baikal region (southeast Siberia), between June 2011 and April 2017. This is the first high‐resolution dataset of stable isotopes of precipitation from this poorly studied region of continental Asia, which has a high potential for isotope‐based paleoclimate research. The dataset revealed distinct seasonal variations: relatively high δ18O (up to –4‰) and δD (up to –40‰) values characterise summer air masses, while lighter isotope composition (–41‰ for δ18O and –322‰ for δD) is characteristic of winter precipitation. Our results show that air temperature mainly affects the isotope composition of precipitation, while no significant correlations were obtained for precipitation amount and relative humidity. A new temperature dependence was established for weighted mean monthly precipitation: +0.50‰/°C (r2 = 0.83; p 〈 0.01; n = 55) for δ18O and +3.8‰/°C (r2 = 0.83, p 〈 0.01; n = 55) for δD. Secondary fractionation processes (e.g. contribution of recycled moisture) were identified mainly in summer from low d excess. Backward trajectories assessed with the HYSPLIT model indicate that precipitation with the lowest mean δ18O and δD values reaches Irkutsk in winter related to moisture transport from the Arctic. Precipitation originating from the west/southwest with the heaviest mean isotope composition reaches Irkutsk in summer, thus representing moisture transport across Eurasia. Generally, moisture transport from the west i.e. the Atlantic Ocean predominates throughout the year. A comparison of our new isotope dataset with simulation results using the ECHAM5‐wiso climate model reveals a good agreement of variations in δ18O (r2 = 0.87; p 〈 0.01; n = 55) and air temperature (r2 = 0.99; p 〈 0.01; n = 71). However, the ECHAM5‐wiso model fails to capture observed variations in d excess (r2 = 0.14; p 〈 0.01; n = 55). This disagreement can be partly explained by a model deficit of capturing regional hydrological processes associated with secondary moisture supply in summer.

Beschreibung: Abstract Mapping of Groundwater‐Dependant Ecosystems (GDEs) relies largely on assumption‐laden evaporation models and few global, direct, and real‐time monitoring techniques exist. We propose a new Synthetic Aperture Radar imagery‐derived index, SARGDE, to identify and monitor these ecosystems across Australia. The index captures vegetation reliance on groundwater during dry periods by estimating the relative stability of foliage and branch structure from the Vertical/Horizontal cross‐polarized band and InSAR coherence. SARGDE is tested over two contrasting study sites in Australia. To build and verify the index, a total of 90 Sentinel‐1 Interferometric Wide images are processed. GDE response to the SAR signal is explored using a non‐linear dimension reduction algorithm. Relevant statistical parameters are derived from data‐cubes and combined to form the index. As the index relies on a one‐year time‐series of globally, freely available, and cloud‐insensitive SAR imagery, SARGDE offers unprecedented capabilities for large‐scale, annual monitoring of GDEs. Such monitoring will aid reconciliation of human and ecosystem groundwater needs by acting as a systematic monitoring tool, helping policy makers to assure ecosystem sustainability where impacts related to mining, agriculture, or climate change may occur.

Beschreibung: Hydrological Processes, EarlyView.

Beschreibung: An integrated hydrologic model is used to simulate watershed hydrodynamics following land cover changes due to a wildfire. Differences between present‐day and postwildfire groundwater pressure show nonlinear increases and decreases that are not spatially limited to burn scar areas. Abstract In recent years, wildfires in the western United States have occurred with increasing frequency and scale. Climate change scenarios in California predict prolonged periods of droughts with even greater potential for conditions amenable to wildfires. The Sierra Nevada Mountains provide 70% of water resources in California, yet how wildfires will impact watershed‐scale hydrology is highly uncertain. In this work, we assess the impacts of wildfires perturbations on watershed hydrodynamics using a physically based integrated hydrologic model in a high‐performance‐computing framework. A representative Californian watershed, the Cosumnes River, is used to demonstrate how postwildfire conditions impact the water and energy balance. Results from the high‐resolution model show counterintuitive feedbacks that occur following a wildfire and allow us to identify the regions most sensitive to wildfires conditions, as well as the hydrologic processes that are most affected. For example, whereas evapotranspiration generally decreases in the postfire simulations, some regions experience an increase due to changes in surface water run‐off patterns in and near burn scars. Postfire conditions also yield greater winter snowpack and subsequently greater summer run‐off as well as groundwater storage in the postfire simulations. Comparisons between dry and wet water years show that climate is the main factor controlling the timing at which some hydrologic processes occur (such as snow accumulation) whereas postwildfire changes to other metrics (such as streamflow) show seasonally dependent impacts primarily due to the timing of snowmelt, illustrative of the integrative nature of hydrologic processes across the Sierra Nevada‐Central Valley interface.

Beschreibung: Abstract Rainfall simulators can enhance our understanding of the hydrologic processes affecting the total runoff to urban drainage systems. This knowledge can be used to improve urban drainage designs. In this study, a rainfall simulator is developed to simulate rainfall on urban green surfaces. The rainfall simulator is controlled by a microcomputer programmed to replicate the temporal variations in rainfall intensity of both historical and synthetic rainfall events with constant rainfall intensity on an area of one square metre. The performance of the rainfall simulator is tested under laboratory conditions with regard to spatial uniformity of the rainfall, the kinetic energy of the raindrops, and the ability to replicate historical and synthetic rainfall events with temporally varying intensity. The rainfall simulator is applied in the field to evaluate its functionality under field conditions and the influence of wind on simulated rainfall. Finally, a field study is carried out on the relationship between runoff, soil volumetric water content, and surface slope. Performance and field tests show that the simulated rainfall has a uniform spatial distribution while the kinetic energy of the raindrops is slightly higher than that of other comparable rainfall simulators. The rainfall simulator performs best in low wind speed conditions. The simulator performs well in replicating historical and synthetic rainfall events by matching both intensity variations and accumulated rainfall depth. The field study shows good correlation between rainfall, runoff, infiltration, soil water content, and surface slope.

Beschreibung: Hydrological Processes, Volume 33, Issue 18, Page 2381-2383, 30 August 2019.

Beschreibung: Municipalities may alter their storm water management focus depending on the most relevant processes (left); the analytical framework developed in this study can be used with measured soil properties to estimate the propensity of urban versus predeveloped reference soil profiles towards saturation‐excess overland flow (SEOF) or infiltration‐excess overland flow (IEOF), with 11 cities in the United States analysed (right). Abstract Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.

Beschreibung: Abstract Permeability reduction of infiltration media due to suspended solid (SS) clogging is the bane of groundwater artificial recharge (GAR). To overcome the clogging problem and advance the understanding of the process‐based spatial‐temporal evolution of SS clogging, a 1D laboratory column simulation was carried out, followed by numerical modelling of the experimental data in this study. It was found that clogging caused a reduction in the hydraulic conductivity (K) in the upper layer at the beginning and extended deeper to approximately 50 cm, and no reduction in K was detected below 52 cm throughout the experimental period of 129 h. The most clogged layer spanned from the surface to a depth of 11 cm, and the middle 11‐52 cm was characterized by a slight decrease in K. The clogging rates of the different layers decreased with the depth, which was based on data analysis, with the largest value of 0.038 h‐1 in the upper 1 cm. The overall K began to decrease from the surface layer and was increasingly affected by clogging with time. A mathematical model was established to simulate the SS clogging process evolution based on considerations of the attachments and detachments of particles. Then, the model was applied to perform several scenario analyses after calibration and validation using the data obtained in the experiment. The simulation results indicated that the SS concentration was much more sensitive than the groundwater depth (GD) below the land surface, and 10 days of constant recharge is recommended as the disposal cycle of the clogged layer under the given conditions.

Beschreibung: Abstract During the last decades, increasing exports of both dissolved organic carbon (DOC) and iron were observed from peat catchments in North America and Europe with potential consequences for water quality of streamwater and carbon storages of soils. As mobilisation and transport processes of DOC and iron in peat catchments are only partly understood, the purpose of this study was to elucidate these processes in an intensively monitored and studied system. Specifically, it was hypothesised that dissimilatory iron reduction in riparian peatland soils mobilises DOC initially adsorbed to iron minerals. During stormflow conditions, both DOC and iron will be transported into the stream network. Ferrous iron may be reoxidised at redox interfaces on its way to the stream and subsequently ferric iron could be transported together with DOC as complexes. To test these hypotheses, generalised additive models (GAM) were applied to 14 years of weekly time series of discharge and concentrations of selected solutes measured in a German headwater stream called Lehstenbach. This stream drains a 4.19 km2 forested mountain catchment, one third of which is covered by riparian peatland soils. We interpreted results of different types of GAM in the way that (a) iron reduction drove the mobilisation of DOC from peatland soils and that (b) both iron and DOC were transported as complexes after their joint mobilisation to and within the steam. It was speculated that low nitrate availability in the uppermost wetland soil layer, particularly during the growing season, promoted iron reduction and thus the mobilisation of DOC. However, the influence of nitrate on the DOC mobilisation remains relatively uncertain. This influence could be further investigated using methods similar to the GAM analysis conducted here for other catchments with long‐term data as well as detailed measurements of the relevant species in riparian wetland soils and the adjacent stream network.

Beschreibung: Hydrological Processes, Volume 33, Issue 16, Page 2157-2159, 30 July 2019.

Beschreibung: Abstract Synchronously and accurately estimating the flood discharges and dynamic changes in the fluid density is essential for hydraulic analysis and forecasting of flash floods, as well as for risk assessment. However, such information is rare for steep mountain catchments, especially in regions that are hotspots for earthquakes. Therefore, six hydrological monitoring sites were established in mainstream and tributaries of the 78.3 km2 Longxi River catchment, an affected region of the Wenchuan earthquake region in China. Direct real‐time monitoring equipment was installed to measure the flow depths, velocities and fluid total pressures of the flood hydrographs. Based on field measurements, real‐time mean cross‐sectional velocities during the flood hydrographs could be derived from easily obtainable parameters: cross‐sectional maximum velocities and the calibrated dimensionless parameter Kh. Real‐time discharges were determined based on this non‐contact method to establish the effective rating curves of this mountainous stream, ranging from 1.46 m3/s to 386.34 m3/s with the RMSEs of ≤10.22 m3/s. Compared to the traditional point‐velocity method and empirical Manning's formula, the proposed non‐contact method was reliable and safe for monitoring whole flood hydrographs. Additionally, the real‐time fluid density during the flood hydrographs was calculated based on the direct monitoring parameters for fluid total pressures and water depths. During the flood hydrograph, transient flow behaviour with higher fluid density generally occurred downstream during the flood peak periods when the flow was in the supercritical flow regime. The observed behaviour greatly increased the threat of damage to infrastructure and human life near the river. Thus, it is important to accurately estimate flood discharge and identify for fluid densities so that people at risk from an impending flash flood are given reliable, advanced warning.

Beschreibung: Hydrological Processes, Volume 33, Issue 15, Page 2045-2047, 15 July 2019.

Beschreibung: A numerical model coupling flow conditions with biochemical reactions is developed. A method for delineating the zones of nitrification and denitrification is proposed. Genetic programming is used to optimize simulation and realize real‐time forecasts. Abstract The hyporheic zone (HZ) plays a vital role in the stream ecosystem. Reactions in the HZ such as denitrification and nitrification have been examined in previous studies. However, no numerical model has yet been developed that can accurately simulate nitrogen concentration changes in the HZ, because the zones for the two reactions can change throughout the reactions. This study proposes a method of evaluating the nitrogen removal rate in the HZ through numerical modelling. First, a basic two‐dimensional numerical model coupling flow conditions with biochemical reactions is proposed to consider both nitrification and denitrification. The zones for different reactions are determined under the assumption that related environmental variables (i.e., the dissolved oxygen) will not change throughout the reactions. Next, to examine changes in environmental variables throughout the reactions, an improved model is proposed, and a method is developed for delineating the boundary between nitrification and denitrification zones and identifying a transition zone where either reaction might take place. However, more information about biochemical reactions in the HZ is required to use the improved model. To overcome this shortcoming, a new model that couples the basic model and genetic programming (GP) is proposed to optimize the simulation results of the basic model and allow for real‐time forecasting. The results show that the basic model obtains acceptable simulation results for nitrate concentration distribution in the HZ. The improved model performs better than the basic model, but the model coupling the basic model with GP performs best. In addition, the function of the HZ in nitrogen removal is examined through a case study of four scenarios, leading to the conclusion that the HZ has a higher nitrogen removal rate when water quality is neither too poor nor too good. Overall, this study enhances our understanding of the HZ and can benefit the restoration and management of HZs and streams in the face of the continual degradation caused by human activity.

Beschreibung: Abstract An automated disk infiltrometer was developed to improve the measurements of soil hydraulic properties (saturated hydraulic conductivity and sorptivity) of soils affected by wildfire. Guideline are given for interpreting curves showing cumulative infiltration as a function of time measured by the autodisk. The autodisk was used to measure the variability of these soil hydraulic properties in three different sample sets: (1) a reference soil consisting of a non‐repellent, uniform, fine sand; (2) soils with the same soil‐textural classification derived from the same bedrock geology but having different initial burn severities; and (3) soils from different bedrock geology but having the same burn severity. The autodisk infiltrometer had greater sampling rates and volume resolution when compared to the visual mini‐disk infiltrometer from previous studies. There was no statistical difference in the mean values measured using the autodisk and visual mini‐disk, but the variability of the autodisk measurements was significantly less than the visual mini‐disk for a given set of samples. The greatest variability of soil hydraulic properties in reference samples with uniform particle size was attributed to different pore geometries (coefficient of variation, COV = 0.28‐0.34). Unburned field samples (same soil type) with heterogeneous particle sizes had greater variability (COV=0.57‐0.78) than the reference samples. However, this basic variability decreased or remained constant in these field samples as burn severity increased. Additional sources of variability (COV=0.53‐1.99) were attributed to multiple‐layers resulting from ash or sediment deposition. Results indicate that resolving differences in soil hydraulic properties from different sites requires more than the common 10 random samples because of the multiple sources of variability.

Beschreibung: Abstract Planning soil conservation strategies requires predictive techniques at event scale because a large percentage of soil loss over a long‐time period is due to relatively few large storms. Considering runoff is expected to improve soil loss predictions and allows relation of the process‐oriented approach with the empirical one, furthermore, the effects of detachment and transport on soil erosion processes can be distinguished by a runoff component. In this paper, the empirical model USLE‐MB (USLE‐M based), including a rainfall‐runoff erosivity factor in which the event rainfall erosivity index EI30 of the Universal Soil Loss Equation (USLE) multiplies the runoff coefficient QR raised to an exponent b1 〉 1 is tested by the measurements carried out for the Masse (10 plots) and Sparacia (22 plots) experimental stations in Italy. For the Masse experimental station, an exponent b1 〉 1 was also estimated by tests carried out by a nozzle‐type rainfall simulator. For each experimental site in fallow conditions, the effect of the sample size of the plot soil loss measurements on the estimate of the b1 coefficient was also studied by the extraction of a fixed number N of randomly obtained pairs of the normalized soil loss and runoff coefficient. The analysis showed that the variability of b1 with N is low and that 350 pairs are sufficient to obtain a stable estimate of b1. A total of 1,262 soil loss data were used to parameterize the model both locally and considering the two sites simultaneously. The b1 exponent varied between the two sites (1.298–1.520), but using a common exponent (1.386) was possible. Using a common b1 exponent for the two experimental areas increases the practical interest for the model and allows the estimation of a baseline component of the soil erodibility factor, which is representative of the at‐site soil intrinsic and quasi‐static properties. Development of a single USLE‐MB model appears possible, and sampling other sites is advisable to develop a single USLE‐MB model for general use.

Beschreibung: Hydrological Processes, Volume 33, Issue 13, Page 1781-1783, 30 June 2019.

Beschreibung: Abstract The root‐zone moisture replenishment mechanisms are key unknowns required to understand soil hydrological processes and water sources used by plants. Temporal patterns of root‐zone moisture replenishment reflect wetting events that contribute to plant growth and survival and to catchment water yield. In this study, stable oxygen and hydrogen isotopes of twigs and throughfall were continuously monitored to characterize the seasonal variations of the root‐zone moisture replenishment in a native vegetated catchment under Mediterranean climate in South Australia. The two studied hillslopes (the north‐facing slope [NFS] and the south‐facing slope [SFS]) had different environmental conditions with opposite aspects. The twig and throughfall samples were collected every ~20 days over 1 year on both hillslopes. The root‐zone moisture replenishment, defined as percentage of newly replenished root‐zone moisture as a complement to antecedent moisture for plant use, calculated by an isotope balance model, was about zero (±25% for the NFS and ± 15% for the SFS) at the end of the wet season (October), increased to almost 100% (±26% for the NFS and ± 29% for the SFS) after the dry season (April and May), then decreased close to zero (±24% for the NFS and ± 28% for the SFS) in the middle of the following wet season (August). This seasonal pattern of root‐zone moisture replenishment suggests that the very first rainfall events of the wet season were significant for soil moisture replenishment and supported the plants over wet and subsequent dry seasons, and that NFS completed replenishment over a longer time than SFS in the wet season and depleted the root zone moisture quicker in the dry season. The stable oxygen isotope composition of the intraevent samples and twigs further confirms that rain water in the late wet season contributed little to root‐zone moisture. This study highlights the significant role of the very first rain events in the early wet season for ecosystem and provides insights to understanding ecohydrological separation, catchment water yield, and vegetation response to climate changes.

Beschreibung: Abstract Preferential flowpaths transport phosphorus (P) to agricultural tile drains. However, if and to what extent this may vary with soil texture, moisture conditions, and P placement is poorly understood. This study investigated (a) interactions between soil texture, antecedent moisture conditions, and the relative contributions of matrix and preferential flow and (b) associated P distributions through the soil profile when fertilizers were applied to the surface or subsurface. Brilliant blue dye was used to stain subsurface flowpaths in clay and silt loam plots during simulated rainfall events under wet and dry conditions. Fertilizer P was applied to the surface or via subsurface placement to plots of different soil texture and moisture condition. Photographs of dye stains were analysed to classify the flow patterns as matrix dominated or macropore dominated, and soils within plots were analysed for their water‐extractable P (WEP) content. Preferential flow occurred under all soil texture and moisture conditions. Dye penetrated deeper into clay soils via macropores and had lower interaction with the soil matrix, compared with silt loam soil. Moisture conditions influenced preferential flowpaths in clay, with dry clay having deeper infiltration (92 ± 7.6 cm) and less dye–matrix interaction than wet clay (77 ± 4.7 cm). Depth of staining did not differ between wet (56 ± 7.2 cm) and dry (50 ± 6.6 cm) silt loam, nor did dominant flowpaths. WEP distribution in the top 10 cm of the soil profile differed with fertilizer placement, but no differences in soil WEP were observed at depth. These results demonstrate that large rainfall events following drought conditions in clay soil may be prone to rapid P transport to tile drains due to increased preferential flow, whereas flow in silt loams is less affected by antecedent moisture. Subsurface placement of fertilizer may minimize the risk of subsurface P transport, particularily in clay.

Beschreibung: Abstract Flood irrigation is globally one of the most used irrigation methods. Typically, not all water that is applied during flood irrigation is consumed by plants or lost to evaporation. Return flow, the portion of applied water from flood irrigation that returns back to streams either via surface or subsurface flow, can constitute a large part of the water balance. Few studies have addressed the connection between vertical and lateral subsurface flows and its potential role in determining return flow pathways due to the difficulty in observing and quantifying these processes at plot or field scale. We employed a novel approach, combining induced polarization, time‐lapse electrical resistivity tomography, and time‐lapse borehole nuclear magnetic resonance, to identify flow paths and quantify changes in soil hydrological conditions under nonuniform application of flood irrigation water. We developed and tested a new method to track the wetting front in the subsurface using the full range of inverted resistivity values. Antecedent soil moisture conditions did not play an important role in preferential flow path activation. More importantly, boundaries between lithological zones in the soil profile were observed to control preferential flow pathways with subsurface run‐off occurring at these boundaries when saturation occurred. Using the new method to analyse time‐lapse resistivity measurements, we were able to track the wetting front and identify subsurface flow paths. Both uniform infiltration and preferential lateral flows were observed. Combining three geophysical methods, we documented the influence of lithology on subsurface flow processes. This study highlights the importance of characterizing the subsurface when the objective is to identify and quantify subsurface return flow pathways under flood irrigation.

Beschreibung: Abstract Evapotranspiration (ET) is an essential component of the hydrological cycle and plays a critical role in water resource management. However, ET is often overlooked in order to transform rainfall to runoff for better streamflow simulation. Hydrological models are commonly used to estimate areal actual evapotranspiration (AET) after calibration against observed discharge. However, classical approaches are often inadequate to appropriately simulate other hydrologic components. Hence, it is important to introduce natural heterogeneity to enhance hydrological processes and reduce water balance errors. In this study, the effectiveness of introducing a constant crop coefficient (Kc), flux tower‐based Kc, and leaf area index (LAI) to three hydrological models (Three‐Parametric Hydrologic Model [TPHM], Génie Rural à 4 paramètres Journalier [GR4J], and Catchment hydrologic cycle Assessment Tool [CAT]) is assessed for the simulation of daily streamflow and AET in a mountainous mixed forest watershed (8.54 km2) in South Korea. The results show that the streamflow simulations after introduction of Kc and LAI to the original models are quite similar. However, the effectiveness of the AET estimation was significantly enhanced after introduction of the flux tower‐based Kc and LAI. The information criterion computed to compare the models reveals that the flux tower‐based Kc‐simulated AET demonstrated better agreement with the observed AET. The Pearson's correlation coefficients (R2) of the TPHM (8%), GR4J (55%), and CAT (55%) models also showed improvements that were greater than the constant based Kc simulation. Similarly, the limitations of the three models with respect to capturing seasonal variation as well as high and low flows were enhanced after the introduction of the flux tower‐based Kc, which adequately reproduced hydrological processes with minimum water balance errors and bias. A regression analysis revealed the potential of estimating Kc as a linear function of LAI (R2 = 0.84). The results of this study indicate that introduction of Kc and LAI to the conceptual rainfall–runoff models can be considered an effective approach to reduce water balance errors and uncertainties in hydrological models and improve the reliability of climate change studies and water resource management.

Beschreibung: Abstract We present the results of a 3‐year monitoring programme of the stable isotope composition of lake water and precipitation at Taozi Lake, in the East Asian monsoon region of China. Our aims were to reveal the spatiotemporal pattern of variation of stable isotopes in a small closed‐basin lake and to quantitatively determine the impacts of precipitation and evaporation on the stable isotope composition of lake water under a humid monsoon climate. In the time domain, the stable oxygen isotopic ratio of the lake water (δ18OL) exhibited substantial seasonal and interannual variations, but the isotope variations between different precipitation events substantially exceeded seasonal and interannual variations. Compared with the stable isotopes in precipitation, δ18OL was substantially positive and dL was negative. In the space domains, the lake water was homogeneously mixed. Indicated by statistic analyses, precipitation plays a dominant role in dynamic of the lake stable isotope during precipitation events of relatively large magnitude, whereas the effect of evaporation is dominant during smaller precipitation events. Results advance our understanding of the stable isotope change rule in the process of lake water evaporation, and it is helpful to identify the climatic significance recorded in stable isotopic compositions of lake bottom sediments.

Beschreibung: Abstract This paper describes surface hydrometeorological and spectral datasets collected from two tower sites located in the University of Melbourne's Dookie experimental farm, Victoria, Australia. The datasets were collected from different vegetation types including wheat, canola, and grazed pasture during the 2012, 2013, and 2014 cropping seasons. The dataset includes 30‐min frequency latent and sensible heat flux measurements and layer‐average soil moisture data at profile depths of 0–5, 0–30, 30–60, 60–90, and 90–120 cm. Air temperature, wind speed, wind direction, relative humidity, precipitation, and incoming and outgoing longwave and shortwave radiation data were also collected from two locations in the study area. The dataset described in this paper is available online.

Beschreibung: Abstract Integrated watershed models can be used to calculate streamflow generation in snow‐dominated mountainous catchments. Parameterization of water flow is often complicated by the lack of information on subsurface hydraulic properties. In this study, bulk density optimization was used to determine hydraulic parameters for the upper and lower regolith in the GEOtop model. The methodology was tested in two small catchments in the Dry Creek Watershed in Idaho and the Libby Creek Watershed in Wyoming. Modelling efficiencies for profile‐average soil–water content for the two catchments were between 0.52 and 0.64. Modelling efficiencies for stream discharge (cumulative stream discharge) were 0.45 (0.91) and 0.54 (0.94) for the Idaho and Wyoming catchments, respectively. The calculated hydraulic properties suggest that lateral flow across the upper–lower regolith interface is an important driver of streamflow in both the Idaho and Wyoming watersheds. The overall calibration procedure is computationally efficient because only two bulk density values are optimized. The two‐parameter calibration procedure was complicated by uncertainty in hydraulic conductivity anisotropy. Different upper regolith hydraulic conductivity anisotropy factors had to be tested in order to describe streamflow in both catchments.

Beschreibung: Abstract The prediction of snowmelt in mountainous forests strongly depends on the accurate description of sensible and latent heat turbulent fluxes. Uncertainty about the within‐canopy wind conditions especially poses a challenge, with relatively few studies examining both above‐ and below‐canopy turbulent fluxes. In this study, turbulent flux predictions from a state‐of‐the‐art watershed model GEOtop were verified against eddy covariance data from one above‐canopy tower and two below‐canopy towers in a snow‐dominated coniferous forest in south‐eastern Wyoming. The model was applied in one‐dimensional vertical mode using field‐observed vegetation parameters and laboratory‐measured soil water retention data. The model was calibrated by identifying optimum values for the canopy fraction and the within‐canopy eddy decay coefficient using the brute‐force method. Above‐canopy sensible heat flux at the Glacier Lakes Ecosystem Experiments Site was predicted reasonably well (r2 = .851). The prediction of above‐canopy latent heat flux was weaker (r2 = .426). For latent heat flux, errors in 30‐min values offset each other when fluxes were aggregated over time, resulting in realistic mean diurnal trends. Below‐canopy turbulent flux at two sites in the Libby Creek Experimental Watershed were predicted with variable success with r2 = .031–.146 for sensible heat flux and r2 = .445–.581 for latent heat flux. Modelled below‐canopy sensible heat flux was too low due to the underestimation of daytime ground surface temperature, because of not enough solar radiation reaching the soil surface. This study suggests that future work on GEOtop and related models should include better parameterizations of the ground surface energy balance to more reliably predict snowmelt and streamflow from mountainous forests.

Beschreibung: Abstract Daily river inflow time series are highly valuable for water resources and water environment management of large lakes. However, the availability of continuous inflow data for large lakes is still relatively limited, especially for large lakes situated within humid plain regions with tens or even hundreds of tributaries. In this study, we choose the fifth largest freshwater Lake Chaohu in China as our study area to introduce a new approach to reconstruct historical daily inflows at ungauged subcatchments of large lakes. This approach makes use of water level, lake surface rainfall, evaporation from the lake, and catchment rainfall observations. Rainfall‐runoff relationship at a reference catchment was analyzed to select rainfall input and estimate runoff coefficient firstly, the runoff coefficient was then transferred to ungauged subcatchments to initially estimate daily inflows. Runoff coefficient was scaled to adjust daily inflows at ungauged subcatchments according to water balance of the lake. This approach was evaluated using sparsely measured inflows at eight subcatchments of Lake Chaohu and compared to the commonly used Drainage Area Ratio (DAR) method. Results suggest that the inflow time series reconstructed from this approach consistent well to corresponding observations, with mean R2 and Nash‐Sutcliffe efficiency (NSE) values of 0.69 and 0.6, respectively. This approach outperforms DAR method in terms of mean R2 and NSE values. Accuracy of this approach holds well when the number of water‐level station being used decreased from four to one.

Beschreibung: Abstract This study analyzes the changes in sediment transport regimes in the middle Yellow River basin (MYRB) using sediment rating parameters. Daily streamflow and suspended sediment concentration (SSC) data were collected at 35 hydrological stations from the 1950s to 2016, which can be divided into three periods based on the type and intensity of human activities: the base stage before 1970, the restraining stage from 1971 to 1989, and the restoration stage after 2002. Data within each period were fitted by log‐linear sediment rating curves (SRC) and the sediment rating parameters were utilized to analyze the spatial and temporal variations in sediment transport regimes. The results show that sediment rating parameters are indicative of sediment transport regimes. In the base stage and the restraining stage, the hydrological stations can be categorized into four groups based on their locations on the rating parameter plot. The stations with small drainage basins were characterized by the highest sediment transport regime, followed by those located in the coarse‐particle zone, the loess zone and the mountainous/forest zone. In the restoration stage, the difference in sediment transport regimes between different geomorphic zones became less distinguishable than in previous stages. During the transition from the base stage to the restraining stage, sediment rating parameters showed no significant changes in sediment transport regimes in all four geomorphic groups. During the transition from the restraining stage to the restoration stage, significant changes were observed in the coarse‐particle zone and the mountain/forest zone, indicating that the revegetation program and large reservoirs imposed a stronger influence on sediment transport regimes in these two zones than in the rest of the MYRB. This study provides theoretical support for evaluating sediment transport regimes with sediment rating parameters.

Beschreibung: Hydrological Processes, EarlyView.

Beschreibung: Abstract In this study, we investigate the impact of the spatial variability of daily precipitation on hydrological projections based on a comparative assessment of streamflow simulations driven by a global climate model (GCM) and two regional climate models (RCMs). A total of 12 different climate input datasets, that is, the raw and bias‐corrected GCM and raw and bias‐corrected two RCMs for the reference and future periods, are fed to a semidistributed hydrological model to assess whether the bias correction using quantile mapping and dynamical downscaling using RCMs can improve streamflow simulation in the Han River basin, Korea. A statistical analysis of the daily precipitation demonstrates that the precipitation simulated by the GCM fails to capture the large variability of the observed daily precipitation, in which the spatial autocorrelation decreases sharply within a relatively short distance. However, the spatial variability of precipitation simulated by the two RCMs shows better agreement with the observations. After applying bias correction to the raw GCM and raw RCMs outputs, only a slight change is observed in the spatial variability, whereas an improvement is observed in the precipitation intensity. Intensified precipitation but with the same spatial variability of the raw output from the bias‐corrected GCM does not improve the heterogeneous runoff distributions, which in turn regulate unrealistically high peak downstream streamflow. GCM‐simulated precipitation with a large bias correction that is necessary to compensate for the poor performance in present climate simulation appears to distort streamflow patterns in the future projection, which leads to misleading projections of climate change impacts on hydrological extremes.

Beschreibung: Abstract Collecting a representative time‐integrated sample of fluvial fine‐grained suspended sediment (〈63 μm) is an important requirement for the understanding of environmental, geomorphological, and hydrological processes operating within watersheds. This study (a) characterized the hydrodynamic behaviour of a commonly used time‐integrated fine sediment sampler (TIFSS) using an acoustic Doppler velocimeter (ADV) in controlled laboratory conditions and (b) measured the mass collection efficiency (MCE) of the sampler by an acoustic Doppler current profiler under field conditions. The laboratory results indicated that the hydrodynamic evaluations associated with the original development of the TIFSS involved an underestimation of the inlet flow velocity of the sampler that results in a significant overestimation of the theoretical MCE. The ADV data illustrated that the ratio of the inlet flow velocity of the sampler to the ambient velocity was 87% and consequently, it can be assumed that a representative sample of the ambient fine suspended particles entered into the sampler. The field results showed that the particle size distribution of the sediment collected by the TIFSS was statistically similar to that for the ambient sediment in the Red River, Manitoba, Canada. The MCE of the TIFSS in the field trials appeared to be as low as 10%. Collecting a representative sample in the field was consistent with the previous findings that the TIFSS is a suitable sampler for the collection of a representative sample of sufficient mass (e.g., 〉1 g) for the investigation of the properties of fluvial fine‐grained suspended sediment. Hydrodynamic evaluation of the TIFSS under a wider range of hydraulic conditions is suggested to assess the performance of the sampler during high run‐off events.

Beschreibung: Abstract Hydrological threshold behaviour has been observed across hillslopes and catchments with varying characteristics. Few studies, however, have evaluated rainfall–run‐off response in areas dominated by agricultural land use and artificial subsurface drainage. Hydrograph analysis was used to identify distinct hydrological events over a 9‐year period and examine rainfall characteristics, dynamic water storage, and surface and subsurface run‐off generation in a drained and farmed closed depression in north‐eastern Indiana, USA. Results showed that both surface flow and subsurface tile flow displayed a threshold relationship with the sum of rainfall amount and soil moisture deficit (SMD). Neither surface flow nor subsurface tile flow was observed unless rainfall amount exceeded the SMD. Timing of subsurface tile flow relative to soil moisture response on the shoulder slope of the depression indicated that the formation and drainage of perched water tables on depression hillslopes were likely the main mechanism that produced subsurface connectivity. Surface flow generation was delayed compared with subsurface tile flow during rainfall events due to differences in soil water storage along depression hillslopes and run‐off generation mechanisms. These findings highlight the substantial impact of subsurface tile drainage on the hydrology of closed depressions; the bottom of the depression, the wettest area prior to drainage installation, becomes the driest part of the depression after installation of subsurface drainage. Rapid connectivity of localized subsurface saturation zones during rainfall events is also greatly enhanced because of subsurface drainage. Thus, less fill is required to generate substantial spill. Understanding hydrologic processes in drained and farmed closed depressions is a critical first step in developing improved water and nutrient management strategies in this landscape.

Beschreibung: Abstract The Guarani Aquifer System (GAS) has been studied since the 1970's, a time frame that coincides with the advent of isotopic techniques in Brazil. The GAS isotope data from many studies are organized in different phases: (i) the advent of isotope techniques; (ii) consolidation and new applications; (iii) isotope assessments and hydrochemistry evolution and (iv) a roadmap to a new conceptual model. The reasons behind the phases, their methodological approaches and impacts on the regional flow conceptual models are examined. Starting with local δ2H and δ18O assessments of values for water fingerprinting and estimates of recharge paleoclimate scenarios, studies evolved to more integrated approaches based on multiple tracers. Stable isotope application techniques were consolidated during the 1980's, when new dating approaches dealing with radiogenic and heavy isotopes were introduced. Through the execution of an international transboundary project, the GAS was studied and extensively sampled for isotopes. These results have triggered wider application of isotope techniques, reflecting also world research trends. Presently, hydrochemical evolution models along flow lines from recharge to discharge areas, across large scale tectonic features within the entire sedimentary basin, are being combined with residence time estimates at GAS outcrop areas and deep confined units. In a complex system it is normal that many, and even contradictory hypotheses are proposed, but isotope techniques provide a unique chance to test them. Stable isotope assessments are still needed near recharge areas and they can be combined with groundwater classical dating procedures, complemented by newer techniques (3H‐3He, CFCs, SF6). Recent noble gas sampling and world pioneer analytical efforts focused on the confined units in the GAS will certainly led to new findings on the overall GAS circulation. The objective of this article is to discuss how isotope information can contribute to the evolution of conceptual groundwater flow models for regional continental aquifers, such as the GAS.

Beschreibung: Abstract The aim of this paper is to quantify peakflow attenuation and/or amplification in a river, investigating lateral flow from the intermediate catchment during floods. This is a challenge for the study of the hydrological response of permeable/intermittent streams, and our contribution refers to a modelling framework based on the inverse problem for the diffusive wave model applied in a karst catchment. Knowing the upstream and downstream hydrographs on a reach between two stations, we can model the lateral one, given information on the hydrological processes involved in the intermediate catchment. The model is applied to 33 flood events in the karst reach of the Iton River in French Normandy where peakflow attenuation is observed. The monitored zone consists of a succession of losing and gaining reaches controlled by strong surface‐water/groundwater (SW/GW) interactions. Our results show that despite a high baseflow increase in the reach, peakflow is attenuated. Model application shows that the intensity of lateral outflow for the flood component is linked to upstream discharge. A combination of river loss and overbank flow for highest floods is proposed for explaining the relationships. Our approach differentiates the role of outflow (river loss and overbank flow) and that of wave diffusion on peakflow attenuation. Based on several sets of model parameterization, diffusion is the main attenuation process for most cases, despite high river losses of up to several m3/s (half of peakflow for some parameterization strategies). Finally, this framework gives new insight into the SW/GW interactions during floods in karst basins, and more globally in basins characterized by disconnected river‐aquifer systems.

Beschreibung: Abstract Groundwater often accounts for a substantial fraction of flood hydrographs, but the processes responsible for this have been unclear. However, many aquifers have preferential flow and this explains how aquifers can be so responsive. In bedrock aquifers, weathering enhances the connectivity and apertures along the most efficient flow paths and hence enhances the permeability. This results in celerities and velocities of the preferential flow in these dual‐porosity aquifers that are two to three orders of magnitude higher than if the aquifers behaved as single‐porosity media. The celerities have been determined from artificial and natural flood pulses, from tidal lags, and from pumping tests. Preferential‐flow velocities have been calculated from tests using applied tracers. Celerities in bedrock aquifers are typically one to two orders of magnitude faster than velocities. The ubiquitous preferential flow in aquifers provides an additional explanation, besides groundwater ridging, for the rapid release of groundwater to streams during storm events.

Beschreibung: Abstract Soil water is an important limiting factor for restoring alpine meadows on the northern Tibetan Plateau. Field studies of soil‐water content (SWC), however, are rare due to the harsh environment, especially in a mesoscale alpine‐meadow ecosystem. The objective of this study was to assess the spatial variability of SWC and the temporal variation of the spatial variability in a typical alpine meadow using a geostatistical approach. SWC was measured using a neutron probe to a depth of 50 cm at 113 locations on 22 sampling occasions in a 33.5 hm2 alpine meadow during the 2015 and 2016 growing seasons. Mean SWC in the study plot for the two growing seasons was 18.7, 14.0, 13.9, 14.3 and 14.8% for depths of 10, 20, 30, 40 and 50 cm, respectively, and SWC was significantly larger at 10 cm than at other depths. SWC was negatively correlated with its spatial variability, and the spatial variability was higher when SWC was lower. Thirty‐three sampling locations in this study plot met the requirement of accuracy of the central limit theorem. A Gaussian model was the best fit for SWC semivariance at depths of 10, 20 and 30 cm, and the spatial structural ratio was between 0.997 and 1, indicating a strong spatial dependence of SWC. The sill and range fluctuated temporally, and the nugget and spatial structural ratio did not generally vary with time. The sill was significantly positively correlated with SWC and was initially stable and then tend to increase with SWC. The nugget, range and spatial structure ratio, however, were not correlated with SWC. These results contribute to our understanding of SWC spatial distribution and variation in alpine meadows and provide basic empirical SWC data for mesoscale model simulations, optimizing sampling strategies and managing meadows on the Tibetan Plateau.

Beschreibung: Hydrological Processes, EarlyView.

Beschreibung: There is great interest in modelling the export of nitrogen (N) and phosphorus (P) from agricultural fields because of ongoing challenges of eutrophication. However, the use of existing hydrochemistry models can be problematic in cold regions because models frequently employ incomplete or conceptually incorrect representations of the dominant cold regions hydrological processes and are over‐parameterised, often with insufficient data for validation. Here, a process‐based N model, WINTRA, which is coupled to a physically based cold regions hydrological model, was expanded to simulate P and account for overwinter soil nutrient biochemical cycling. An inverse modelling approach, using this model with consideration of parameter equifinality, was applied to an intensively monitored agricultural basin in Manitoba, Canada, to help identify the main climate, soil and anthropogenic controls on nutrient export. Consistent with observations, the model results suggest that snow water equivalent, melt rate, snow cover depletion rate and contributing area for runoff generation determine the opportunity time and surface area for runoff‐soil interaction. These physical controls have not been addressed in existing models. Results also show that the time lag between the start of snowmelt and the arrival of peak nutrient concentration in runoff increased with decreasing antecedent soil moisture content, highlighting potential implications of frozen soils on runoff processes and hydrochemistry. The simulations showed TDP concentration peaks generally arriving earlier than NO3, but also decreasing faster afterwards, which suggests a significant contribution of plant residue TDP to early snowmelt runoff. Antecedent fall tillage and fertiliser application increased TDP concentrations in spring snowmelt runoff but did not consistently affect NO3 runoff. In this case, the antecedent soil moisture content seemed to have had a dominant effect on overwinter soil N biogeochemical processes such as mineralization, which are often ignored in models. This work demonstrates both the need for better representation of cold regions processes in hydrochemical models and the model improvements that are possible if these are included.

Beschreibung: Abstract We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels.

Beschreibung: Abstract Biocrust effects on soil infiltration have attracted increasing attention in dryland ecosystems, but their seasonal variations in infiltrability have not yet been well understood. On the Chinese Loess Plateau, soil infiltrability indicated by saturated hydraulic conductivity (Ks) of biocrusts and bare soil, both on aeolian sand and loess soil, was determined by disc infiltrometer in late spring (SPR), midsummer (SUM), and early fall (FAL). Then their correlations with soil biological and physiochemical properties and water repellency index (RI) were analysed. The results showed that the biocrusts significantly decreased Ks both on sand during SPR, SUM, and FAL (by 43%, 66%, and 35%, respectively; P 〈 .05) and on loess (by 42%, 92%, and 10%, respectively; P 〈.05). As compared with the bare soil, the decreased Ks in the biocrusted surfaces was mostly attributed to the microorganism biomass and also to the increasing content of fine particles and organic matter. Most importantly, both the biocrusts and bare soil exhibited significant (F ≥ 11.89, P ≤ .003) seasonal variations in Ks, but their patterns were quite different. Specifically, the Ks of bare soil gradually decreased from SPR to SUM (32% and 42% for sand and loess, respectively) and FAL (29% and 39%); the Ks of biocrusts also decreased from SPR to SUM (59% and 92%) but then increased in FAL (36% and 588%). Whereas the seasonal variations in Ks of the biocrusts were closely correlated with the seasonal variations in RI, the RI values were not high enough to point at hydrophobicity. Instead of that, the seasonal variations of Ks were principally explained by the changes in the crust biomass and possibly by the microbial exopolysaccharides. We conclude that the biocrusts significantly decreased soil infiltrability and exhibited a different seasonal variation pattern, which should be carefully considered in future analyses of hydropedological processes.

Beschreibung: Abstract Parametrization of transformation and transport processes of dissolved organic carbon (DOC) in soils is challenging especially under variable hydrological conditions. In this study, DOC concentrations in stormflow were analyzed with a physically‐based modeling approach. A one‐dimensional dual‐continuum vertical flow and transport model was applied to simulate subsurface processes in a macroporous forest hillslope soil over a period of 4.5 years. Microbially mediated transformations of DOC were assumed to depend primarily on soil moisture and soil temperature. Two conceptually different descriptions of the sorption of DOC to soil were examined with equilibrium and kinetic approaches. In order to quantify the uncertainties associated with the model parameterization, Monte Carlo analyses in conjunction with Latin hypercube sampling was performed. Despite the complexity of microbial transformations, the simulated temporal patterns of DOC concentrations in stormflow showed similar behavior to those reflected in the observed DOC fluxes. Due to preferential flow, the hillslope DOC export (5.0 ± 0.5 g C m−2 y−1) was higher than the amounts usually reported in the literature. Overall DOC transport in hillslope scenarios could be described appropriately using the equilibrium sorption assumption. The performed analyses showed that the inclusion of the kinetic description of DOC sorption only slightly improved the predictions of the DOC hillslope export. Moreover, influences of seasonal hydro‐climatological conditions on hillslope export of DOC could be observed. Reduced DOC transport during an extreme warm and dry summer was described with lower accuracy, thus indicating the difficulties in representation DOC transformations under dry conditions.

Beschreibung: Abstract As the largest hydroelectric dam in the world, the Three Gorges Dam (TGD) has raised wide concerns over the environmental and ecological impacts since its dramatic effect on the downstream flow regime of the Yangtze River. Since 2003, the TGD has progressed from the initial operation period to the full operation period, with different effects on the downstream flow regime over each period. Although the upstream inflow change (USIC) of the TGD is a possibly additional factor affecting the downstream flow regime, this has drawn little attention. This study aims to quantify the individual contributions of the TGD and the USIC to the changes of the downstream flow regime over different operation periods of the dam. Using the Muskingum routing model and the Xin'anjiang rainfall‐runoff model, we reconstruct the discharge unregulated by the TGD for the post‐TGD period from 2003 to 2015. On this basis, the effects of the TGD and the USIC on the downstream flow regime are quantitatively assessed. Benchmarked against the flow regime during the pre‐TGD period from 1955 to 2002, it is found that the TGD and the USIC play considerable and comparable roles in affecting the downstream flow regime during the whole post‐TGD period from 2003 to 2015. Furthermore, the TGD appears to have a limited effect on the downstream flow regime during the initial operation period from 2003 to 2008 relative to the USIC. In contrast, during the full operation period from 2009 to 2015, the TGD plays a dominant role in changing the downstream flow regime, although the effect of the USIC cannot be neglected. The findings of this study are helpful to understand the exact impacts of the TGD on the downstream flow regime, thereby facilitating the development of a rational strategy for operating the dam.

Beschreibung: Abstract Human‐accelerated climate change is quickly leading to glacier‐free mountains, with consequences for the ecology and hydrology of alpine river systems. Water origin (i.e. glacier, snowmelt, precipitation, groundwater) is a key control on multiple facets of alpine stream ecosystems, since it drives the physico‐chemical template of the habitat in which ecological communities reside and interact, and ecosystem processes occur. Accordingly, distinct alpine stream types and associated communities have been identified. However, unlike streams fed by glaciers (i.e. kryal), groundwater (i.e. krenal), and snowmelt‐precipitation (i.e. rhithral), those fed by rock glaciers are still poorly documented. We characterized the physical and chemical features of these streams and investigated the influence of rock glaciers on the habitat template of alpine river networks. We analysed two subcatchments in a deglaciating area of the Central European Alps, where rock glacier‐fed, groundwater‐fed, and glacier‐fed streams are all present. We monitored the spatial, seasonal, and diel variability of physical conditions (i.e. water temperature, turbidity, channel stability, discharge) and chemical variables (electrical conductivity, major ions and trace element concentrations) during the snowmelt, glacier ablation, and flow recession periods of two consecutive years. We observed distinct physical and chemical conditions and seasonal responses for the different stream types. Rock glacial streams were characterized by very low and constant water temperatures, stable channels, clear waters, and high concentrations of ions and trace elements that increased as summer progressed. Furthermore, one rock glacier strongly influenced the habitat template of downstream waters due to high solute export, especially in late summer under increased permafrost thaw. Given their unique set of environmental conditions, we suggest that streams fed by thawing rock glaciers are distinct river habitats that differ from those normally classified for alpine streams. Rock glaciers may become increasingly important in shaping the hydroecology of alpine river systems under continued deglaciation.

Beschreibung: Abstract For most of the year, a dry‐bed desert wash is void of water flow. Intensive rain events, however, could trigger significant flash floods that bring about highly complicated hydro‐ and morphodynamics processes within a desert stream. We present a fully coupled three‐phase flow model of air, water, and sediment to simulate numerically the propagation of a flash flood in a field‐scale fluvial desert stream, the so‐called Tex Wash located in the Mojave Desert, California, United States. The turbulent flow of the flash flood is computed using the three‐dimensional unsteady Reynolds‐averaged Navier–Stokes (URANS) equations closed with the shear stress transport (SST) k − ω model. The free surface of the flash flood at the interface of air and water phases is computed with the level‐set method, which enables instantaneous tracking of the water surface as the flash flood propagates over the dry bed of the desert stream. The evolution of the desert fluvial stream's morphology, due to the action of the propagating flash flood on the mobile bed, is calculated using a Eulerian morphodynamics model based on the curvilinear immersed boundary method. The capabilities of the proposed numerical framework are demonstrated by applying it to simulate a flash flood event in a 0.65‐km‐long reach of the Tex Wash, the intricate channel morphology of which is obtained using light imaging detection and ranging (LiDAR) technology. The simulated region of the stream includes a number of bridge foundations. The simulation results of the model for the flash flood event revealed the formation of a highly complex flow field and scour patterns within the stream. Moreover, our simulation results showed that most scour processes take place during the steady phase of the flash flood, i.e. after the flash flood fills the stream. The transient phase of the flash flood is rather short and contributes to a very limited amount of erosion within the desert stream.

Beschreibung: Abstract Stable water isotope ratios are measured as a tracer of environmental processes in materials such as leaves, soils, and lakes. Water in these archives may experience evaporation, which increases the abundance of heavy isotopologues in a manner proportional to the gradients in humidity and isotope ratio between the evaporating water and the surrounding atmosphere. Until recently, however, the isotope ratio of the atmosphere has been difficult to measure, and measurements remain scarce. As a result, several assumptions have been adopted to estimate isotope ratios in atmospheric water vapor. Perhaps the most commonly employed assumption in terrestrial environments is that water vapor is in isotopic equilibrium with precipitation. We evaluate this assumption using an 8‐member ensemble of general circulation model (GCM) simulations that include explicit calculation of isotope ratios in precipitation and vapor. We find that across the model ensemble water vapor is typically less depleted in heavy isotopologues than expected if it were in equilibrium with annual precipitation. Atmospheric vapor likely possesses higher‐than‐expected isotope ratios because precipitation isotope ratios are determined by atmospheric conditions that favor condensation, which do not reflect atmospheric mixing and advection processes outside of precipitation events. The effect of this deviation on theoretical estimates of evaporation isotope ratios scales with relative humidity. As a result, the equilibrium assumption gives relatively accurate estimates of the isotope ratios of evaporation in low latitudes, but performs increasingly poorly at increasing latitudes. Future studies of evaporative water pools should include measurements of atmospheric isotope ratios or constrain potential bias with isotope‐enabled GCM simulations.

Beschreibung: Abstract Temporal variation of runoff chemistry and its seasonal controls relating to chemical weathering processes and drainage system evolution were examined at Urumqi Glacier No.1 (UG1) in Xinjiang, China over a full melt season. The dominant ions in meltwater runoff are HCO3‐, Ca2+ and SO42‐; and Fe, Sr and Al are dominant elements. Concentrations of major ions and some elements show periodic variations with seasons and negatively correlate with discharge, while other elements (e.g. Al, Ni, Cu, Zn, Cd, and Pb) show a random change, providing insights to the hydrological and physicochemical controls. HCO3‐ and Ca2+ are primarily derived from calcite, SO42‐ and Fe mainly come from pyrite, and Sr and Al principally originate from silicate. Hydrochemical fluxes of solutes exhibit strong seasonality but are positively related to discharge, suggesting an increasing release of solutes during higher flow conditions. Solute yields, CDR and CWI observed at UG1 are higher than those at most basins worldwide. This suggests that chemical weathering in central Asia may be stronger than at other glacial basins with similar specific discharge. Concentrations of some elements (e.g. Fe, Al, As, Pb, and Zn) are close to or exceed the guidelines for drinking water standards in meltwater‐fed rivers. These rivers may face future challenges of water quality degradation, and relationships between changing flow and water quality conditions should be established soon, given that development of channelized flow is expected to be earlier over a melt season in a warming climate.

Beschreibung: Abstract Solute concentrations in streamflow typically vary systematically with stream discharge, and the resulting concentration‐discharge relationships are important signatures of catchment biogeochemical processes. Solutes derived from mineral weathering often exhibit decreasing concentrations with increasing flows, suggesting dilution of a kinetically limited weathering flux by a variable flux of water. However, previous work showed that concentration‐discharge relationships of weathering‐derived solutes in 59 headwater catchments were much weaker than this simple dilution model would predict. Instead, catchments behaved as chemostats, with rates of solute production and/or mobilization that were nearly proportional to water fluxes, on both event and inter‐annual time scales. Here we re‐examine these findings using data for a wider range of solutes from 2186 catchments, ranging from ~10 to 〉1,000,000 km2 in drainage area, and spanning a wide range of lithologic and climatic settings. Concentration‐discharge relationships among this much larger set of larger catchments are broadly consistent with the previously described chemostatic behavior, at least on event and inter‐annual time scales for weathering‐derived solutes. Among these same catchments, however, site‐to‐site variations in mean concentrations of weathering‐derived solutes exhibit strong negative correlations with long‐term average precipitation and discharge, reflecting strong climatic control on long‐term leaching of the critical zone. We use multiple regression of site characteristics including discharge to identify potential controls on long‐term mean concentrations, and find that lithologic and land cover controls are significant predictors for many analytes. The picture that emerges is one in which, on event and inter‐annual time scales, weathering‐derived stream solute concentrations are chemostatically buffered by groundwater storage and fast chemical reactions, but each catchment's chemostatic "set point" reflects site‐to‐site differences in climatically driven evolution of the critical zone. In contrast to these weathering products, some nutrients and particulates are often near‐chemostatic across all timescales, and their long‐term mean concentrations correlate more strongly with land use than climatic characteristics.

Beschreibung: Abstract For an erosion event (October 2016) occurred at the Sparacia experimental area (Southern Italy), both terrestrial and low‐altitude aerial surveys were carried out by consumer grade camera and quadcopter (low‐cost unmanned aerial vehicle [UAV]) to measure rill erosion on two plots with steepness of 22% and 26%. Applying the structure from motion (SfM) technique, the three‐dimensional digital terrain models (3D‐DTMs) and the quasi three‐dimensional models (2.5D‐digital elevation model [DEM]) were obtained by the two surveys. Furthermore, 3D‐DTM and DEM were built using the available aerial photographs (166) and adding 40 terrestrial photographs. For the first time, the convergence index was applied to high‐resolution rill data for extracting the rill network, and a subsequent separation into contributing and non‐contributing rills was carried out. The comparison among the three surveys (terrestrial, UAV, and UAV + terrestrial) was developed using two morphometric parameters of the rill network (drainage density and drainage frequency). Moreover, using as reference the weight of sediment stored on the tanks located downstream of the plots, the reliability of soil loss measurement by 3D models was tested. For both contributing and non‐contributing rills, the morphometric parameters were higher for the terrestrial than for UAV and UAV + terrestrial surveys. For both plots, SfM always provided reliable soil loss measurements, which were affected by errors ranging from −8% to 13%. Although the applied technique used a low‐cost UAV and a consumer grade camera, the obtained results demonstrated that a reliable estimate of rill erosion can be obtained in an area of interest.

Beschreibung: Abstract Characterization of stable isotope compositions (δ2H and δ18O) of surface water and groundwater in a catchment is critical for refining moisture sources and establishing modern isotope–elevation relationships for paleoelevation reconstructions. There is no consensus on the moisture sources of precipitation in the Yellow River source region during summer season. This study presents δ2H and δ18O data from 111 water samples collected from tributaries, mainstream, lakes, and groundwater across the Yellow River source region during summertime. Measured δ18O values of the tributaries range from −13.5‰ to −5.8‰ with an average of −11.0‰. Measured δ18O values of the groundwater samples range from −12.7‰ to −10.5‰ with an average of −11.9‰. The δ18O data of tributary waters display a northward increase of 1.66‰ per degree latitude. The δ18O data and d‐excess values imply that moisture sources of the Yellow River source region during summertime are mainly from the mixing of the Indian Summer Monsoon and the Westerlies, local water recycling, and subcloud evaporation. Analysis of tributary δ18O data from the Yellow River source region and streamwater and precipitation δ18O data from its surrounding areas leads to a best‐fit second‐order polynomial relationship between δ18O and elevation over a 4,600 m elevation range. A δ18O elevation gradient of −1.6‰/km is also established using these data, and the gradient is in consistence with the δ18O elevation gradient of north and eastern plateau. Such relationships can be used for paleoelevation reconstructions in the Yellow River source region.

Beschreibung: Abstract The flux of fluvial carbon from the terrestrial biosphere to the world's oceans is known to be an important component of the global carbon cycle, but within this pathway, the flux and return of carbon to the river network via sewage effluent has not been quantified. In this study, monitoring data from 2000 to 2016 for the dissolved organic carbon (DOC) concentration, biochemical oxygen demand, and chemical oxygen demand of the final effluent of sewage treatment works from across England were examined to assess the amount of DOC contributing to national‐scale fluvial fluxes of carbon. The study shows that the median concentration of DOC in final effluent was 9.4 compared with 4.8 mg C/L for all surface waters for the United Kingdom over the study period and that the DOC in final effluent significantly declined over the study period from 11.0 to 6.4 mg C/L. Rivers receiving sewage effluent showed a significant, on average 19%, increase in DOC concentration downstream of sewage discharges. At the scale of the United Kingdom, the flux of DOC in final effluent was 31 ktonnes C/year with a per capita export of 0.55 kg C/year and compared with an average annual flux of DOC from the United Kingdom of 859 ktonnes C/year, that is, only 3.6% of national‐scale flux. The lability of this DOC was limited, with only 7.4% loss of final effluent DOC concentration over in‐stream residence times of up to 5 days. The direct decline in DOC concentration from sewage treatment works was not large enough on its own to explain the declines observed in DOC concentration in U.K. rivers at their tidal limit.

Beschreibung: Abstract Native Nothofagus forests in the midlatitude region of the Andes Cordillera are notorious biodiversity hot spots, uniquely situated in the Southern Hemisphere such that they develop in snow‐dominated reaches of this mountain range. Spanning a smaller surface area than similar ecosystems, where forests and snow coexist in the Northern Hemisphere, the interaction between vegetation and snow processes in this ecotone has received lesser attention. We present the first systematic study of snow–vegetation interactions in the Nothofagus forests of the Southern Andes, focusing on how the interplay between interception and climate determines patterns of snow water equivalent (SWE) variability. The Valle Hermoso experimental catchment, located in the Nevados de Chillán vicinity, was fitted with eight snow depth sensors that provided continuous measurements at varying elevations, aspect, and forest cover. Also, manual measurements of snow properties were obtained during snow surveys conducted during end of winter and spring seasons for 3 years, between 2015 and 2017. Each year was characterized by distinct climatological conditions, with 2016 representing one of the driest winters on record in this region. Distance to canopy, leaf area index, and total gap area were measured at each observational site. A regression model was built on the basis of statistical analysis of local parameters to model snow interception in this kind of forest. We find that interception implied a 23.2% reduction in snow accumulation in forested sites compared with clearings. The interception in these deciduous trees represents, on average, 23.6% of total annual snowfall, reaching a maximum measured interception value of 13.8‐mm SWE for all snowfall events analysed in this research.

Beschreibung: Abstract Cushion plant dominated peatlands are key ecosystems in tropical alpine regions of the Andes in South America. The cushion plants have formed peat bodies over thousands of years that fill many valley bottoms, and the forage produced by the plants is critical for native and nonnative domesticated mammals. The sources and flow paths of water supporting these peatlands remain largely unknown. Some studies have suggested that glacier meltwater streams support some peatlands, and that the ongoing loss of glaciers and their meltwaters could lead to the loss or diminishment of peatlands. We analysed the hydrologic regime of 10 peatlands in four mountain regions of Bolivia and Peru using groundwater monitoring. Groundwater levels in peatlands were relatively stable and within 20 cm of the ground surface during the rainy season, and many sites had water tables 40–90 cm below the ground surface in the dry season. Topographic and groundwater elevations in the peatlands demonstrated that the water source of all 10 peatlands was hillslope groundwater flowing from lateral moraines, talus, colluvium, or bedrock aquifers into the peatlands. There was little to no input from streams, whether derived from glacier melt or other sources, and glacier melt could not have recharged the hillslope aquifers supporting peatlands. We measured the stable water isotopes in water samples taken during different seasons, distributed throughout the catchments, and the values are consistent with this interpretation. Our findings indicate that peatlands in the study region are recharged by hillslope groundwater discharge rather than stream water and may not be as vulnerable to glacial decline as other studies have indicated. However, both glaciers and peatlands are susceptible to changing thermal and precipitation regimes that could affect the persistence of peatlands.

Beschreibung: Abstract Headwater streams (HSs) are generally naturally prone to flow intermittence. These intermittent rivers and ephemeral streams have recently seen a marked increase in interest, especially to assess the impact of drying on aquatic ecosystems. The two objectives of this work are (a) to identify the main drivers of flow intermittence dynamics in HS and (b) to reconstruct local daily drying dynamics. Discrete flow states—“flowing” versus “drying”—are modelled as functions of covariates that include information on climate, hydrology, groundwater, and basin descriptors. Three classifiers to estimate flow states using covariates are tested on four contrasted regions in France: (a) a linear classifier with regularization (LASSO for least absolute shrinkage and selection operator) and two non‐linear non‐parametric classifiers, (b) a one‐hidden‐layer feedforward artificial neural network (ANN) classifier, and (c) a random forest (RF) classifier. The three classifiers are compared with a benchmark classifier (BC) that simply estimates dominant flow state for each month based on observations (without using covariates). The performance assessment over the period 2012–2016 carried out by cross‐validation shows that the three classifiers for flow state based on covariates outperformed the BC. This demonstrates the predictive power of the covariates. ANN is the classifier that globally achieves the best performance to predict the daily drying dynamics whereas both RF and LASSO tend to underestimate the proportion of drying states. The covariates are ranked in terms of relevance for each classifier. The monthly proportion of drying states provided by the discrete observation network has a major importance for the three classifiers ANN, LASSO, and RF. This may reflect the proclivity of a site to flow intermittence. ANN gives higher importance to climatic and hydrological covariates and its non‐linearity allows a greater degree of freedom.

Beschreibung: Abstract Movement of soil moisture associated with tree root‐water uptake is ecologically important but technically challenging to measure. Here, the self‐potential (SP) method, a passive electrical geophysical method, is used to characterize water flow in situ. Unlike tensiometers, which use a measurement of state (i.e., matric pressure) at two locations to infer fluid flow, the SP method directly measures signals generated by water movement. We collected SP measurements in a two‐dimensional array at the base of a Douglas‐fir tree (Pseudotsuga menziesii) in the H.J. Andrews Experimental Forest in western Oregon over 5 months to provide insight on the propagation of transpiration signals into the subsurface under variable soil moisture. During dry conditions, SP data appear to show downward unsaturated flow, whereas nearby tensiometer data appear to suggest upward flow during this period. After the trees enter dormancy in the fall, precipitation‐induced vertical flow dominates in the SP and tensiometer data. Diel variations in SP data correspond to periods of tree transpiration. Changes in volumetric water content occurring from soil moisture movement during transpiration are not large enough to appear in volumetric water content data. Fluid flow and electrokinetic coupling (i.e., electrical potential distribution) were simulated using COMSOL Multiphysics to explore the system controls on field data. The coupled model, which included a root‐water uptake term, reproduced components of both the long‐term and diel variations in SP measurements, thus indicating that SP has potential to provide spatially and temporally dense measurements of transpiration‐induced changes in water flow. This manuscript presents the first SP measurements focusing on the movement of soil moisture in response to tree transpiration.

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

Beschreibung: Abstract With the intensification of climate change, its impact on runoff variations cannot be ignored. The main purpose of this study is to analyze the nonstationarity of runoff frequency adjusted for future climate change in the Luanhe River basin, China, and quantify the different sources of uncertainties in nonstationary runoff frequency analysis. The advantage of our method is the combination of generalized additive models in location, scale and shape (GAMLSS) and downscaling models. The nonstationary GAMLSS models were established for the nonstationary frequency analysis of runoff (1961–2010) by using the observed precipitation as a covariate, which is closely related to runoff and contributes significantly to its nonstationarity. To consider the nonstationary effects of future climate change on future runoff variations, the downscaled precipitation series in the future (2011–2080) from the general circulation models (GCMs) were substituted into the selected nonstationary model to calculate the statistical parameters and runoff frequency in the future. A variance decomposition method was applied to quantify the impacts of different sources of uncertainty on the nonstationary runoff frequency analysis. The results show that the impacts of uncertainty in the GCMs, scenarios and statistical parameters of the GAMLSS model increase with increasing runoff magnitude. In addition, GCMs and GAMLSS model parameters have the main impacts on runoff uncertainty, accounting for 14% and 83% of the total uncertainty sources, respectively. Conversely, the interactions and scenarios make limited contributions, accounting for 2% and 1%, respectively. Further analysis shows that the sources of uncertainty in the statistical parameters of the nonstationary model mainly result from the fluctuations in the precipitation sequence. This result indicates the necessity of considering the precipitation sequence as a covariate for runoff frequency analysis in the future.

Beschreibung: Abstract Fog phenomena and their associated meteorological variables were continuously monitored during four years in an evergreen laurisilva cloud forest of the Anaga Massif Biosphere Reserve (Tenerife, Canary Islands), in order to establish its current dynamics. Fog was more frequent during night through early morning and in the afternoon, and particularly from May until September, coincidental with a frequent immersion of the 1025 m a.s.l. experimental site in the cloud layer of wind‐driven stratocumulus. The concomitant meteorological conditions during different fog regimes, characterized according to visibility (Ω) ranges, were compared with those when fog was absent. The presence of fog was associated with a significant reduction in global solar radiation, Rg, increased wind speed and lower and more stable ambient temperatures. The foggy vs. fog‐free hourly medians of Rg were found to be linearly related, while the proportion of median Rg reduction due to fog varied logarithmically with Ω. However, foggy vs. fog‐free extreme values of the hourly Rg distributions departed from such a linear trend. By contrast, hourly temperatures during foggy versus fog‐free periods behaved linearly for most of the Ω range, except for very dense fog, Ω 〈 100 m. Transpiration of the canopy, intermittently wetted due to interception of both rain and fog water droplets, was determined by quantifying the water balance at leaf scale with a mathematical model for the two representative hypostomatous species present at the site: the arboreal shrub Erica platycodon, with needle‐like leaves, and the laurophyll tree Myrica faya. Both tree transpiration and evaporation of the intercepted fog water were predictively higher during summer. By contrast, transpiration was reduced during February, in agreement with a one year period of sap velocity measurements, and was not appreciably affected by soil moisture content. The consequences of an anticipated downward shift of the stratocumulus cloud layer and of various projected Representative Concentration Pathways (RCPs) scenarios in the Macaronesian area were simulated, yielding in all cases a significant rise in transpiration for both species. Particularly the simulated RCPs scenarios implied 29‐73% increments in transpiration from the actual values. Since fog is concomitant with lower temperatures and vapour pressure deficit, the modification of its current distribution as a consequence of climate change may have a direct effect on such associated meteorological variables, and therefore a meaningful impact in the water relations of the laurel cloud forests.

Beschreibung: Abstract A synthetic storm generator – Dynamic Moving Storm (DMS) – is developed in this study to represent spatio‐temporal variabilities of rainfall and storm movement in synthetic storms. Using an urban watershed as the testbed, the authors investigate the hydrologic responses to the DMS parameters and their interactions. In order to reveal the complex nature of rainfall‐runoff processes, previously simplified assumptions are relaxed in this study regarding 1) temporal variability of rainfall intensity and 2) time‐invariant flow velocity in channel routing. The results of this study demonstrate the significant contribution of storm moving velocity to the variation of peak discharge based on a global sensitivity analysis (GSA). Furthermore, a pairwise sensitivity analysis (PSA) is conducted to elucidate not only the patterns in individual contributions from parameters to hydrologic responses but also their interactions with storm moving velocity. The intricacies of peak discharges resulting from sensitivity analyses are then dissected into independent hydrologic metrics, i.e. runoff volume and standard deviation of runoff timings, for deeper insights. It is confirmed that peak discharge is increased when storms travel downstream along the main channel at the speed that corresponds to a temporal superposition of runoff. Spatial concentration of catchment rainfall is found to be a critical linkage through which characteristics of moving storms affect peak discharges. In addition, altering peak timing of rainfall intensity in conjunction with storm movement results in varied storm core locations in the channel network, which further changes the flow attenuation effects from channel routing. For future directions, the DMS generator will be embedded in a stochastic modeling framework and applied in rainfall/flow frequency analysis.

Beschreibung: Abstract Although it is well known that the vast majority of the time only a portion of any watershed contributes run‐off to the outlet, this extent is rarely documented. Also, the power law form of the streamflow and contributing area (Q‐Ac) relationship has been known for a half century, but it is uncommon for it to be quantified, and time series of contributing area extensive enough to calculate its frequency distribution are almost non‐existent. Data from the Canadian Prairies, where there are extensive estimates of contributing area during the median annual flood, imply that the power law coefficient for any Q‐Ac curve is a function of flow magnitude and return period. These data also suggest that regional flood frequency curves are a construct of Q‐Ac curves from individual basins. This paper will discuss research that attempted to reproduce the Q‐Ac curves for the La Salle River Watershed with a semidistributed numerical hydrological model, MESH‐PDMROF. The model simulated streamflow reasonably well (Nash Sutcliffe values = 0.62) compared with published examples of comparable models applied in the region. Estimates of the coefficient and exponent of the Q‐Ac power law function ranged from 0.08–0.14 and 0.9–1.12, respectively. These exponent values were lower than those of regional flood frequency curves and support the theory that regional flood frequency curves are a construct of Q‐Ac curves. Simulations of the area contributing to the median annual flood were lower (0.3) than those derived from independent topographic analysis (0.9) described in earlier literature though there is uncertainty in both these estimates. This uncertainty was extended across the flood frequency distribution and may be too large to definitively verify the study hypothesis.

Beschreibung: Abstract This work proposes two modelling frameworks for diagnosing temporal variations in nonlinear rating curves that describe suspended sediment–discharge relationships. A variant of the weighted regression on time, discharge, and season model is proposed and is compared against dynamic nonlinear modelling, a newly developed nonlinear time series filter based on sequential Monte Carlo sampling. Both approaches estimate a time series of rating curve parameters, with uncertainty, that can be used to diagnose variability in the sediment–discharge relationship over time. We evaluate the models with a variety of synthetic scenarios to highlight their ability to estimate signals of known rating curve change. Results reveal important bias‐variance trade‐offs unique to each approach, and in general, suggest that dynamic nonlinear modelling is better suited for rapid rating curve changes, whereas the weighted regression on time, discharge, and season variant more precisely estimates slow change. The techniques are then applied in two case studies in the Upper Hudson and Mohawk Rivers in New York. We conclude with a discussion of the implications of dynamic rating curves for the management of water quality in riverine and estuary systems.

Beschreibung: Abstract Recent studies have highlighted the agronomic and environmental importance of phosphorus (P) movement through the soil profile. Thus, faced with challenges such as high‐profile cases of P enrichment of surface water, better understanding of nutrient movement through soil is needed to better manage agricultural fertilizers and manures and their contribution to water quality degradation. In particular, field‐scale research is especially needed in soils with preferential flow transport pathways. Thus, we collected nitrogen (N) and P transport data in run‐off and seepage (lateral subsurface return flow) from 13 field‐ and farm‐scale watersheds on Vertisols in Central Texas for a 14‐year period. For 2004–2017, seepage accounted for ~20% of the total surface flow, and nutrient concentrations were generally similar in run‐off and seepage. As surface run‐off contributed ~80% of the flow, it follows that median annual N and P loads in run‐off were significantly greater than in seepage for every watershed. N loads in both run‐off and seepage flow from cultivated land were an order of magnitude greater than in native prairie and improved pasture, and the highest run‐off and seepage P loads both occurred on cultivated land with organic fertilizer sources. Increasing watershed scale (size) did not to produce consistent patterns in N or P loss in run‐off or seepage. Land use and watershed scale produced significant differences in seepage volume but did not affect run‐off volumes or total surface flow/rainfall. Although less significant in terms of total offsite flux, nutrient movement in vadose zones has important agronomic and environmental implications as considerable N and P are transported through and within the root zone and eventually offsite. And in terms of P, this contradicts the traditionally held scientific viewpoint that P movement through the vadose zone is unimportant agronomically and environmentally.