Refine
Year of publication
Document Type
- Article (32)
- Postprint (11)
- Doctoral Thesis (1)
- Review (1)
Language
- English (45)
Is part of the Bibliography
- yes (45)
Keywords
- Nordeste (4)
- Precipitation (4)
- time-series (4)
- variability (4)
- Connectivity (3)
- Erosion (3)
- seasonality (3)
- snow (3)
- Climate (2)
- Forecasting Framework (2)
Saturated hydraulic conductivity (K-s) is an important soil characteristic affecting soil water storage, runoff generation and erosion processes. In some areas where high-intensity rainfall coincides with low K-s values at shallow soil depths, frequent overland flow entails dense drainage networks. Consequently, linear structures such as flowlines alternate with inter-flowline areas. So far, investigations of the spatial variability of K-s mainly relied on isotropic covariance models which are unsuitable to reveal patterns resulting from linear structures. In the present study, we applied two sampling approaches so as to adequately characterize K-s spatial variability in a tropical forest catchment that features a high density of flowlines: A classical nested sampling survey and a purposive sampling strategy adapted to the presence of flowlines. The nested sampling approach revealed the dominance of small-scale variability, which is in line with previous findings. Our purposive sampling, however, detected a strong spatial gradient: surface K-s increased substantially as a function of distance to flowline; 10 m off flowlines, values were similar to the spatial mean of K-s. This deterministic trend can be included as a fixed effect in a linear mixed modeling framework to obtain realistic spatial fields of K-s. In a next step we used probability maps based on those fields and prevailing rainfall intensities to assess the hydrological relevance of the detected pattern. This approach suggests a particularly good agreement between the probability statements of K-s exceedance and observed overland flow occurrence during wet stages of the rainy season.
Forests seem to represent low-erosion systems, according to most, but not all, studies of suspended-sediment yield. We surmised that this impression reflects an accidental bias in the selection of monitoring sites towards those with prevailing vertical hydrological flowpaths, rather than a tight causal link between vegetation cover and erosion alone. To evaluate this conjecture, we monitored, over a 2-year period, a 3.3 ha old-growth rainforest catchment prone to frequent and widespread overland flow. We sampled stream flow at two and overland flow at three sites in a nested arrangement on a within-event basis, and monitored the spatial and temporal frequency of overland flow. Suspended-sediment concentrations were modeled with Random Forest and Quantile Regression Forest to be able to estimate the annual yields for the 2 years, which amounted to 1 t ha(-1) and 2 t ha(-1) in a year with below-average and with average precipitation, respectively. These estimates place our monitoring site near the high end of reported suspended-sediment yields and lend credence to the notion that low yields reflect primarily the dominance of vertical flowpaths and not necessarily and exclusively the kind of vegetative cover. Undisturbed forest and surface erosion are certainly no contradiction in terms even in the absence of mass movements.
Soils in various places of the Panama Canal Watershed feature a low saturated hydraulic conductivity (K-s) at shallow depth, which promotes overland-flow generation and associated flashy catchment responses. In undisturbed forests of these areas, overland flow is concentrated in flow lines that extend the channel network and provide hydrological connectivity between hillslopes and streams. To understand the dynamics of overland-flow connectivity, as well as the impact of connectivity on catchment response, we studied an undisturbed headwater catchment by monitoring overland-flow occurrence in all flow lines and discharge, suspended sediment, and total phosphorus at the catchment outlet. We find that connectivity is strongly influenced by seasonal variation in antecedent wetness and can develop even under light rainfall conditions. Connectivity increased rapidly as rainfall frequency increased, eventually leading to full connectivity and surficial drainage of entire hillslopes. Connectivity was nonlinearly related to catchment response. However, additional information on factors such as overland-flow volume would be required to constrain relationships between connectivity, stormflow, and the export of suspended sediment and phosphorus. The effort to monitor those factors would be substantial, so we advocate applying the established links between rain event characteristics, drainage network expansion by flow lines, and catchment response for predictive modeling and catchment classification in forests of the Panama Canal Watershed and in similar regions elsewhere.
An effective strategy for combining variance- and distribution-based global sensitivity analysis
(2020)
We present a new strategy for performing global sensitivity analysis capable to estimate main and interaction effects from a generic sampling design. The new strategy is based on a meaningful combination of varianceand distribution-based approaches. The strategy is tested on four analytic functions and on a hydrological model. Results show that the analysis is consistent with the state-of-the-art Saltelli/Jansen formula but to better quantify the interaction effect between the input factors when the output distribution is skewed. Moreover, the estimation of the sensitivity indices is much more robust requiring a smaller number of simulations runs. Specific settings and alternative methods that can be integrated in the new strategy are also discussed. Overall, the strategy is considered as a new simple and effective tool for performing global sensitivity analysis that can be easily integrated in any environmental modelling framework.
Cosmic-ray neutron sensing (CRNS) is a powerful technique for retrieving representative estimates of soil water content at a horizontal scale of hectometres (the “field scale”) and depths of tens of centimetres (“the root zone”). This study demonstrates the potential of the CRNS technique to obtain spatio-temporal patterns of soil moisture beyond the integrated volume from isolated CRNS footprints. We use data from an observational campaign carried out between May and July 2019 that featured a dense network of more than 20 neutron detectors with partly overlapping footprints in an area that exhibits pronounced soil moisture gradients within one square kilometre. The present study is the first to combine these observations in order to represent the heterogeneity of soil water content at the sub-footprint scale as well as between the CRNS stations. First, we apply a state-of-the-art procedure to correct the observed neutron count rates for static effects (heterogeneity in space, e.g. soil organic matter) and dynamic effects (heterogeneity in time, e.g. barometric pressure). Based on the homogenized neutron data, we investigate the robustness of a calibration approach that uses a single calibration parameter across all CRNS stations. Finally, we benchmark two different interpolation techniques for obtaining spatio-temporal representations of soil moisture: first, ordinary Kriging with a fixed range; second, spatial interpolation complemented by geophysical inversion (“constrained interpolation”). To that end, we optimize the parameters of a geostatistical interpolation model so that the error in the forward-simulated neutron count rates is minimized, and suggest a heuristic forward operator to make the optimization problem computationally feasible. Comparison with independent measurements from a cluster of soil moisture sensors (SoilNet) shows that the constrained interpolation approach is superior for representing horizontal soil moisture gradients at the hectometre scale. The study demonstrates how a CRNS network can be used to generate coherent, consistent, and continuous soil moisture patterns that could be used to validate hydrological models or remote sensing products.
Modelers can improve a model by addressing the causes for the model errors (data errors and structural errors). This leads to implementing model enhancements (MEs), for example, meteorological data based on more monitoring stations, improved calibration data, and/or modifications in process formulations. However, deciding on which MEs to implement remains a matter of expert knowledge. After implementing multiple MEs, any improvement in model performance is not easily attributed, especially when considering different objectives or aspects of this improvement (e.g., better dynamics vs. reduced bias). We present an approach for comparing the effect of multiple MEs based on real observations and considering multiple objectives (MMEMO). A stepwise selection approach and structured plots help to address the multidimensionality of the problem. Tailored analyses allow a differentiated view on the effect of MEs and their interactions. MMEMO is applied to a case study employing the mesoscale hydro-sedimentological model WASA-SED for the Mediterranean-mountainous Isabena catchment, northeast Spain. The investigated seven MEs show diverse effects: some MEs (e.g., rainfall data) cause improvements for most objectives, while other MEs (e.g., land use data) only affect a few objectives or even decrease model performance. Interaction of MEs was observed for roughly half of the MEs, confirming the need to address them in the analysis. Calibration and increasing the temporal resolution showed by far stronger impact than any of the other MEs. The proposed framework can be adopted in other studies to analyze the effect of MEs and, thus, facilitate the identification and implementation of the most promising MEs for comparable cases.
A comprehensive hydro-sedimentological dataset for the Isabena catchment, northeastern (NE) Spain, for the period 2010-2018 is presented to analyse water and sediment fluxes in a Mediterranean mesoscale catchment. The dataset includes rainfall data from 12 rain gauges distributed within the study area complemented by meteorological data of 12 official meteo-stations. It comprises discharge data derived from water stage measurements as well as suspended sediment concentrations (SSCs) at six gauging stations of the River Isabena and its sub-catchments. Soil spectroscopic data from 351 suspended sediment samples and 152 soil samples were collected to characterize sediment source regions and sediment properties via fingerprinting analyses. The Isabena catchment (445 km(2)) is located in the southern central Pyrenees ranging from 450 m to 2720 m a.s.l.; together with a pronounced topography, this leads to distinct temperature and precipitation gradients. The River Isabena shows marked discharge variations and high sediment yields causing severe siltation problems in the downstream Barasona Reservoir. The main sediment source is badland areas located on Eocene marls that are well connected to the river network. The dataset features a comprehensive set of variables in a high spatial and temporal resolution suitable for the advanced process understanding of water and sediment fluxes, their origin and connectivity and sediment budgeting and for the evaluation and further development of hydro-sedimentological models in Mediterranean mesoscale mountainous catchments.
Cosmic-ray neutron sensing (CRNS) is a powerful technique for retrieving representative estimates of soil water content at a horizontal scale of hectometres (the “field scale”) and depths of tens of centimetres (“the root zone”). This study demonstrates the potential of the CRNS technique to obtain spatio-temporal patterns of soil moisture beyond the integrated volume from isolated CRNS footprints. We use data from an observational campaign carried out between May and July 2019 that featured a dense network of more than 20 neutron detectors with partly overlapping footprints in an area that exhibits pronounced soil moisture gradients within one square kilometre. The present study is the first to combine these observations in order to represent the heterogeneity of soil water content at the sub-footprint scale as well as between the CRNS stations. First, we apply a state-of-the-art procedure to correct the observed neutron count rates for static effects (heterogeneity in space, e.g. soil organic matter) and dynamic effects (heterogeneity in time, e.g. barometric pressure). Based on the homogenized neutron data, we investigate the robustness of a calibration approach that uses a single calibration parameter across all CRNS stations. Finally, we benchmark two different interpolation techniques for obtaining spatio-temporal representations of soil moisture: first, ordinary Kriging with a fixed range; second, spatial interpolation complemented by geophysical inversion (“constrained interpolation”). To that end, we optimize the parameters of a geostatistical interpolation model so that the error in the forward-simulated neutron count rates is minimized, and suggest a heuristic forward operator to make the optimization problem computationally feasible. Comparison with independent measurements from a cluster of soil moisture sensors (SoilNet) shows that the constrained interpolation approach is superior for representing horizontal soil moisture gradients at the hectometre scale. The study demonstrates how a CRNS network can be used to generate coherent, consistent, and continuous soil moisture patterns that could be used to validate hydrological models or remote sensing products.
The characteristics of a landscape pose essential factors for hydrological processes. Therefore, an adequate representation of the landscape of a catchment in hydrological models is vital. However, many of such models exist differing, amongst others, in spatial concept and discretisation. The latter constitutes an essential pre-processing step, for which many different algorithms along with numerous software implementations exist. In that context, existing solutions are often model specific, commercial, or depend on commercial back-end software, and allow only a limited or no workflow automation at all. Consequently, a new package for the scientific software and scripting environment R, called lumpR, was developed. lumpR employs an algorithm for hillslope-based landscape discretisation directed to large-scale application via a hierarchical multi-scale approach. The package addresses existing limitations as it is free and open source, easily extendible to other hydrological models, and the workflow can be fully automated. Moreover, it is user-friendly as the direct coupling to a GIS allows for immediate visual inspection and manual adjustment. Sufficient control is furthermore retained via parameter specification and the option to include expert knowledge. Conversely, completely automatic operation also allows for extensive analysis of aspects related to landscape discretisation. In a case study, the application of the package is presented. A sensitivity analysis of the most important discretisation parameters demonstrates its efficient workflow automation. Considering multiple streamflow metrics, the employed model proved reasonably robust to the discretisation parameters. However, parameters determining the sizes of subbasins and hillslopes proved to be more important than the others, including the number of representative hillslopes, the number of attributes employed for the lumping algorithm, and the number of sub-discretisations of the representative hillslopes.
The semiarid northeast of Brazil is one of the most densely populated dryland regions in the world and recurrently affected by severe droughts. Thus, reliable seasonal forecasts of streamflow and reservoir storage are of high value for water managers. Such forecasts can be generated by applying either hydrological models representing underlying processes or statistical relationships exploiting correlations among meteorological and hydrological variables. This work evaluates and compares the performances of seasonal reservoir storage forecasts derived by a process-based hydrological model and a statistical approach.
Driven by observations, both models achieve similar simulation accuracies. In a hindcast experiment, however, the accuracy of estimating regional reservoir storages was considerably lower using the process-based hydrological model, whereas the resolution and reliability of drought event predictions were similar by both approaches. Further investigations regarding the deficiencies of the process-based model revealed a significant influence of antecedent wetness conditions and a higher sensitivity of model prediction performance to rainfall forecast quality.
Within the scope of this study, the statistical model proved to be the more straightforward approach for predictions of reservoir level and drought events at regionally and monthly aggregated scales. However, for forecasts at finer scales of space and time or for the investigation of underlying processes, the costly initialisation and application of a process-based model can be worthwhile. Furthermore, the application of innovative data products, such as remote sensing data, and operational model correction methods, like data assimilation, may allow for an enhanced exploitation of the advanced capabilities of process-based hydrological models.