ALBERT

Description: Recent climate change and vegetation greening have important implications for global terrestrial hydrological cycles and other ecosystem functions, raising concerns about the watershed water supply capacity for large water diversion projects. To address this emerging concern, we built a hybrid model based on the Coupled Carbon and Water (CCW) and Water Supply Stress Index (WaSSI) models and conducted a case study on the upper Han River basin (UHRB) in Central China that serves as the water source area to the middle route of the South-to-North Water Diversion Project (SNWDP). Significant vegetation greening occurred in the UHRB during 2001–2018, largely driven by the widespread afforestation in the region, with the normalized difference vegetation index increasing at a rate of 0.5±0.1 % yr−1 (p

Description: Determining the spatiotemporal variability in the annual streamflow volume plays a relevant role in hydrology with regard to improving and implementing sustainable and resilient policies and practices of water resource management. This study investigates annual streamflow volume trends in a newly assembled, consolidated, and validated data set of daily mean river flow records from more than 3000 stations which cover near-natural basins in more than 40 countries across Europe. Although the data set contains streamflow time series from 1900 to 2013 in some stations, the statistical analyses were carried out by including observations from 1950 to 2013 in order to have a consistent and reliable data set over the continent. Trends were detected by calculating the slope of the Theil–Sen line over the annual anomalies of streamflow volume. The results show that annual streamflow volume trends have emerged at European scale, with a marked negative tendency in Mediterranean regions, with about -1×103 m3/(km2 yr−2), and a generally positive trend in northern ones, with about 0.5×103 m3/(km−2 yr−2). The annual streamflow volume trend patterns appear to be in agreement with the continental-scale meteorological observations in response to climate change drivers. In the Mediterranean area, the decline of annual streamflow volumes started in 1965, and since the early 1980s, volumes have consistently been lower than the 1950–2013 average. The spatiotemporal annual streamflow volume patterns observed in this work can help to contextualize short-term trends and regional studies already available in the scientific literature, as well as to provide a valid benchmark for further accurate quantitative analysis of annual streamflow volumes.

Description: The effect of the spatial resolution of digital terrain models (DTMs) on topography and soil erosion modelling is well documented for low resolutions. Nowadays, the availability of high spatial resolution DTMs from unmanned aerial vehicles (UAVs) opens new horizons for detailed assessment of soil erosion with hydrological models, but the effects of DTM resolution on model outputs at this scale have not been systematically tested. This study combines plot-scale soil erosion measurements, UAV-derived DTMs, and spatially explicit soil erosion modelling to select an appropriate spatial resolution based on allowable loss of information. During 39 precipitation events, sediment and soil samples were collected on five bounded and unbounded plots and four land covers (forest, fallow, maize, and eroded bare land). Additional soil samples were collected across a 220 ha watershed to generate soil maps. Precipitation was collected by two rain gauges and vegetation was mapped. A total of two UAV campaigns over the watershed resulted in a 0.60 m spatial-resolution DTM used for resampling to 1, 2, 4, 8, and 15 m and a multispectral orthomosaic to generate a land cover map. The OpenLISEM model was calibrated at plot level at 1 m resolution and then extended to the watershed level at the different DTM resolutions. Resampling the 1 m DTM to lower resolutions resulted in an overall reduction in slope. This reduction was driven by migration of pixels from higher to lower slope values; its magnitude was proportional to resolution. At the watershed outlet, 1 and 2 m resolution models exhibited the largest hydrograph and sedigraph peaks, total runoff, and soil loss; they proportionally decreased with resolution. Sedigraphs were more sensitive than hydrographs to spatial resolution, particularly at the highest resolutions. The highest-resolution models exhibited a wider range of predicted soil loss due to their larger number of pixels and steeper slopes. The proposed evaluation method was shown to be appropriate and transferable for soil erosion modelling studies, indicating that 4 m resolution (

Description: This study employs a stochastic hydrologic modeling framework to evaluate the sensitivity of flood frequency analyses to different components of the hydrologic modeling chain. The major components of the stochastic hydrologic modeling chain, including model structure, model parameter estimation, initial conditions, and precipitation inputs were examined across return periods from 2 to 100 000 years at two watersheds representing different hydroclimates across the western USA. A total of 10 hydrologic model structures were configured, calibrated, and run within the Framework for Understanding Structural Errors (FUSE) modular modeling framework for each of the two watersheds. Model parameters and initial conditions were derived from long-term calibrated simulations using a 100 member historical meteorology ensemble. A stochastic event-based hydrologic modeling workflow was developed using the calibrated models in which millions of flood event simulations were performed for each basin. The analysis of variance method was then used to quantify the relative contributions of model structure, model parameters, initial conditions, and precipitation inputs to flood magnitudes for different return periods. Results demonstrate that different components of the modeling chain have different sensitivities for different return periods. Precipitation inputs contribute most to the variance of rare floods, while initial conditions are most influential for more frequent events. However, the hydrological model structure and structure–parameter interactions together play an equally important role in specific cases, depending on the basin characteristics and type of flood metric of interest. This study highlights the importance of critically assessing model underpinnings, understanding flood generation processes, and selecting appropriate hydrological models that are consistent with our understanding of flood generation processes.

Description: Many recent studies have sought to characterize variations of the annual maximum flood discharge series over time and across space in Europe, including some that have elucidated different process controls on different statistical properties of these series. To further support these studies, we conduct a pan-European assessment of process controls on key properties of this series, including the mean annual flood (MAF) and coefficients of variation (CV) and skewness (CS) of flood discharges. These annual maximum flood discharge series consist of instantaneous peaks and daily means observed in 2370 catchments in Europe without strong human modifications covering the period 1960–2010. We explore how the estimated moments MAF, CV and CS vary due to catchment size, climate and other controls across Europe, where their averages are 0.17 m3 s−1 km−2, 0.52 and 1.28, respectively. The results indicate that MAF is largest along the Atlantic coast, in the high-rainfall areas of the Mediterranean coast and in mountainous regions, while it is smallest in the sheltered parts of the East European Plain. The CV is largest in southern and eastern Europe, while it is smallest in the regions subject to strong Atlantic influence. The pattern of the CS is similar, albeit more erratic, in line with the greater sampling variability of CS. In the Mediterranean, MAF, CV and CS decrease strongly with catchment area, suggesting that floods in small catchments are relatively very large, while in eastern Europe this dependence is much weaker, mainly due to more synchronized timing of snowmelt over large areas. The process controls on the flood moments in five predetermined hydroclimatic regions are identified through correlation and multiple linear regression analyses with a range of covariates, and the interpretation is aided by a seasonality analysis. Precipitation-related covariates are found to be the main controls of the spatial patterns of MAF in most of Europe except for regions in which snowmelt contributes to MAF, where air temperature is more important. The Aridity Index is, by far, the most important control on the spatial pattern of CV in all of Europe. Overall, the findings suggest that, at the continental scale, climate variables dominate over land surface characteristics, such as land use and soil type, in controlling the spatial patterns of flood moments. Finally, to provide a performance baseline for more local studies, we assess the estimation accuracy of regional multiple linear regression models for estimating flood moments in ungauged basins.

Description: Collection efficiency transfer functions that compensate for wind-induced collection loss are presented and evaluated for unshielded precipitation gauges. Three novel transfer functions with wind speed and precipitation fall velocity dependence are developed, including a function from computational fluid dynamics modelling (CFD), an experimental fall velocity threshold function (HE1), and an experimental linear fall velocity dependence function (HE2). These functions are evaluated alongside universal (KUniversal) and climate-specific (KCARE) transfer functions with wind speed and temperature dependence. Transfer function performance is assessed using 30 min precipitation event accumulations reported by unshielded and shielded Geonor T-200B3 precipitation gauges over two winter seasons. The latter gauge was installed in a Double Fence Automated Reference (DFAR) configuration. Estimates of fall velocity were provided by the Precipitation Occurrence Sensor System (POSS). The CFD function reduced the RMSE (0.08 mm) relative to KUniversal (0.20 mm), KCARE (0.13 mm), and the unadjusted measurements (0.24 mm), with a bias error of 0.011 mm. The HE1 function provided a RMSE of 0.09 mm and bias error of 0.006 mm, capturing the collection efficiency trends for rain and snow well. The HE2 function better captured the overall collection efficiency, including mixed precipitation, resulting in a RMSE of 0.07 mm and bias error of 0.006 mm. These functions are assessed across solid and liquid hydrometeor types and for temperatures between −22 and 19 ∘C. The results demonstrate that transfer functions incorporating hydrometeor fall velocity can dramatically reduce the uncertainty of adjusted precipitation measurements relative to functions based on temperature.

Description: Carbon-14 (14C) is routinely used to determine mean residence times (MRTs) of groundwater. 14C-based MRT calculations typically assume that the unsaturated zone is in equilibrium with the atmosphere, controlling the input 14C activity. However, multiple studies have shown that unsaturated zone 14C activities are lower than atmospheric values. Despite the availability of unsaturated zone 14C data, no attempt has been made to generalise initial 14C activities with depth to the water table. We utilise measurements of unsaturated zone 14C activities from 13 studies to produce a 14C–depth relationship to estimate initial 14C activities. The technique only requires the depth to the water table at the time of sampling or an estimate of depth to water in the recharge zone to determine the input 14C activity, making it straightforward to apply. Applying this new relationship to two Australian datasets (113 14C measurements in groundwater) shows that MRT estimates were up to 9250 years younger when the 14C–depth correction was applied relative to conventional MRTs. These findings may have important implications for groundwater samples that suggest the mixing of young and old waters and the determination of the relative proportions of young and waters, whereby the estimated fraction of older water may be much younger than previously assumed. Owing to the simplicity of the application of the technique, this approach can be easily incorporated into existing correction schemes to assess the sensitivity of unsaturated zone 14C to MRTs derived from 14C data.

Description: Temporal (serial) clustering of extreme precipitation events on sub-seasonal timescales is a type of compound event. It can cause large precipitation accumulations and lead to floods. We present a novel, count-based procedure to identify episodes of sub-seasonal clustering of extreme precipitation. We introduce two metrics to characterise the prevalence of sub-seasonal clustering episodes and their contribution to large precipitation accumulations. The procedure does not require the investigated variable (here precipitation) to satisfy any specific statistical properties. Applying this procedure to daily precipitation from the ERA5 reanalysis data set, we identify regions where sub-seasonal clustering occurs frequently and contributes substantially to large precipitation accumulations. The regions are the east and northeast of the Asian continent (northeast of China, North and South Korea, Siberia and east of Mongolia), central Canada and south of California, Afghanistan, Pakistan, the southwest of the Iberian Peninsula, and the north of Argentina and south of Bolivia. Our method is robust with respect to the parameters used to define the extreme events (the percentile threshold and the run length) and the length of the sub-seasonal time window (here 2–4 weeks). This procedure could also be used to identify temporal clustering of other variables (e.g. heat waves) and can be applied on different timescales (sub-seasonal to decadal). The code is available at the listed GitHub repository.

Description: While the availability and affordability of unmanned aerial systems (UASs) has led to the rapid development of remote sensing applications in hydrology and hydrometry, uncertainties related to such measurements must be quantified and mitigated. The physical instability of the UAS platform inevitably induces motion in the acquired videos and can have a significant impact on the accuracy of camera-based measurements, such as velocimetry. A common practice in data preprocessing is compensation of platform-induced motion by means of digital image stabilisation (DIS) methods, which use the visual information from the captured videos – in the form of static features – to first estimate and then compensate for such motion. Most existing stabilisation approaches rely either on customised tools developed in-house, based on different algorithms, or on general purpose commercial software. Intercomparison of different stabilisation tools for UAS remote sensing purposes that could serve as a basis for selecting a particular tool in given conditions has not been found in the literature. In this paper, we have attempted to summarise and describe several freely available DIS tools applicable to UAS velocimetry. A total of seven tools – six aimed specifically at velocimetry and one general purpose software – were investigated in terms of their (1) stabilisation accuracy in various conditions, (2) robustness, (3) computational complexity, and (4) user experience, using three case study videos with different flight and ground conditions. In an attempt to adequately quantify the accuracy of the stabilisation using different tools, we have also presented a comparison metric based on root mean squared differences (RMSDs) of inter-frame pixel intensities for selected static features. The most apparent differences between the investigated tools have been found with regards to the method for identifying static features in videos, i.e. manual selection of features or automatic. State-of-the-art methods which rely on automatic selection of features require fewer user-provided parameters and are able to select a significantly higher number of potentially static features (by several orders of magnitude) when compared to the methods which require manual identification of such features. This allows the former to achieve a higher stabilisation accuracy, but manual feature selection methods have demonstrated lower computational complexity and better robustness in complex field conditions. While this paper does not intend to identify the optimal stabilisation tool for UAS-based velocimetry purposes, it does aim to shed light on details of implementation, which can help engineers and researchers choose the tool suitable for their needs and specific field conditions. Additionally, the RMSD comparison metric presented in this paper can be used in order to measure the velocity estimation uncertainty induced by UAS motion.

Description: Sorptivity is a parameter of primary importance in the study of unsaturated flow in soils. This hydraulic parameter is required to model water infiltration into vertical soil profiles. Sorptivity can be directly estimated from the soil hydraulic functions (water retention and hydraulic conductivity curves), using the integral formulation of Parlange (1975). However, calculating sorptivity in this manner requires the prior determination of the soil hydraulic diffusivity and its numerical integration between initial and final saturation degrees, which may be difficult in some situations (e.g., coarse soil with diffusivity functions that are quasi-infinite close to saturation). In this paper, we present a procedure to compute sorptivity using a scaling parameter, cp, that corresponds to the sorptivity of a unit soil (i.e., unit values for all parameters and zero residual water content) that is utterly dry at the initial state and saturated at the final state. The cp parameter was computed numerically and analytically for five hydraulic models: delta (i.e., Green and Ampt), Brooks and Corey, van Genuchten–Mualem, van Genuchten–Burdine, and Kosugi. Based on the results, we proposed brand new analytical expressions for some of the models and validated previous formulations for the other models. We also tabulated the output values so that they can easily be used to determine the actual sorptivity value for any case. At the same time, our numerical results showed that the relation between cp and the hydraulic shape parameters strongly depends on the chosen model. These results highlight the need for careful selection of the proper model for the description of the water retention and hydraulic conductivity functions when estimating sorptivity.