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)
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.
To support scientifically sound water management in dryland environments a modelling system has been developed for the quantitative assessment of water and sediment fluxes in catchments, transport in the river system, and retention in reservoirs. The spatial scale of interest is the mesoscale because this is the scale most relevant for management of water and land resources.
This modelling system comprises process-oriented hydrological components tailored for dryland characteristics coupled with components comprising hillslope erosion, sediment transport and reservoir deposition processes. The spatial discretization is hierarchically designed according to a multi-scale concept to account for particular relevant process scales. The non-linear and partly intermittent run-off generation and sediment dynamics are dealt with by accounting for connectivity phenomena at the intersections of landscape compartments. The modelling system has been developed by means of data from nested research catchments in NE-Spain and in NE-Brazil.
In the semi-arid NE of Brazil sediment retention along the topography is the main process for sediment retention at all scales, i.e. the sediment delivery is transport limited. This kind of deposition retains roughly 50 to 60 % of eroded sediment, maintaining a similar deposition proportion in all spatial scales investigated. On the other hand, the sediment retained in reservoirs is clearly related to the scale, increasing with catchment area. With increasing area, there are more reservoirs, increasing the possibility of deposition. Furthermore, the area increase also promotes an increase in flow volume, favouring the construction of larger reservoirs, which generally overflow less frequently and retain higher sediment fractions. The second example comprises a highly dynamic Mediterranean catchment in NE-Spain with nested sub-catchments and reveals the full dynamics of hydrological, erosion and deposition features. The run-off modelling performed well with only some overestimation during low-flow periods due to the neglect of water losses along the river. The simulated peaks in sediment flux are reproduced well, while low-flow sediment transport is less well captured, due to the disregard of sediment remobilization in the riverbed during low flow.
This combined observation and modelling study deepened the understanding of hydro-sedimentological systems characterized by flashy run-off generation and by erosion and sediment transport pulses through the different landscape compartments. The connectivity between the different landscape compartments plays a very relevant role, regarding both the total mass of water and sediment transport and the transport time through the catchment.
A set of seasonal drought forecast models was assessed and verified for the Jaguaribe River in semiarid northeastern Brazil. Meteorological seasonal forecasts were provided by the operational forecasting system used at FUNCEME (Ceará’s research foundation for meteorology)and by the European Centre for Medium-Range Weather Forecasts (ECMWF). Three downscaling approaches (empirical quantile mapping, extended downscaling and weather pattern classification) were tested and combined with the models in hindcast mode for the period 1981 to 2014. The forecast issue time was January and the forecast period was January to June. Hydrological drought indices were obtained by fitting a multivariate linear regression to observations. In short, it was possible to obtain forecasts for (a) monthly precipitation,(b) meteorological drought indices, and (c) hydrological drought indices. The skill of the forecasting systems was evaluated with regard to root mean square error (RMSE), the Brier skill score (BSS) and the relative operating characteristic skill score (ROCSS). The tested forecasting products showed similar performance in the analyzed metrics. Forecasts of monthly precipitation had little or no skill considering RMSE and mostly no skill with BSS. A similar picture was seen when forecasting meteorological drought indices: low skill regarding RMSE and BSS and significant skill when discriminating hit rate and false alarm rate given by the ROCSS (forecasting drought events of, e.g., SPEI1 showed a ROCSS of around 0.5). Regarding the temporal variation of the forecast skill of the meteorological indices, it was greatest for April, when compared to the remaining months of the rainy season, while the skill of reservoir volume forecasts decreased with lead time. This work showed that a multi-model ensemble can forecast drought events of timescales relevant to water managers in northeastern Brazil with skill. But no or little skill could be found in the forecasts of monthly precipitation or drought indices of lower scales, like SPI1. Both this work and those here revisited showed that major steps forward are needed in forecasting the rainy season in northeastern Brazil.
A set of seasonal drought forecast models was assessed and verified for the Jaguaribe River in semiarid northeastern Brazil. Meteorological seasonal forecasts were provided by the operational forecasting system used at FUNCEME (Ceará’s research foundation for meteorology)and by the European Centre for Medium-Range Weather Forecasts (ECMWF). Three downscaling approaches (empirical quantile mapping, extended downscaling and weather pattern classification) were tested and combined with the models in hindcast mode for the period 1981 to 2014. The forecast issue time was January and the forecast period was January to June. Hydrological drought indices were obtained by fitting a multivariate linear regression to observations. In short, it was possible to obtain forecasts for (a) monthly precipitation,(b) meteorological drought indices, and (c) hydrological drought indices. The skill of the forecasting systems was evaluated with regard to root mean square error (RMSE), the Brier skill score (BSS) and the relative operating characteristic skill score (ROCSS). The tested forecasting products showed similar performance in the analyzed metrics. Forecasts of monthly precipitation had little or no skill considering RMSE and mostly no skill with BSS. A similar picture was seen when forecasting meteorological drought indices: low skill regarding RMSE and BSS and significant skill when discriminating hit rate and false alarm rate given by the ROCSS (forecasting drought events of, e.g., SPEI1 showed a ROCSS of around 0.5). Regarding the temporal variation of the forecast skill of the meteorological indices, it was greatest for April, when compared to the remaining months of the rainy season, while the skill of reservoir volume forecasts decreased with lead time. This work showed that a multi-model ensemble can forecast drought events of timescales relevant to water managers in northeastern Brazil with skill. But no or little skill could be found in the forecasts of monthly precipitation or drought indices of lower scales, like SPI1. Both this work and those here revisited showed that major steps forward are needed in forecasting the rainy season in northeastern Brazil.
Water shortage is a serious threat for many societies worldwide. In drylands, water management measures like the construction of reservoirs are affected by eroded sediments transported in the rivers. Thus, the capability of assessing water and sediment fluxes at the river basin scale is of vital importance to support management decisions and policy making. This subject was addressed by the DFG-funded SESAM-project (Sediment Export from large Semi-Arid catchments: Measurements and Modelling). As a part of this project, this thesis focuses on (1) the development and implementation of an erosion module for a meso-scale catchment model, (2) the development of upscaling and generalization methods for the parameterization of such model, (3) the execution of measurements to obtain data required for the modelling and (4) the application of the model to different study areas and its evaluation. The research was carried out in two meso-scale dryland catchments in NE-Spain: Ribera Salada (200 km²) and Isábena (450 km²). Adressing objective 1, WASA-SED, a spatially semi-distributed model for water and sediment transport at the meso-scale was developed. The model simulates runoff and erosion processes at the hillslope scale, transport processes of suspended and bedload fluxes in the river reaches, and retention and remobilisation processes of sediments in reservoirs. This thesis introduces the model concept, presents current model applications and discusses its capabilities and limitations. Modelling at larger scales faces the dilemma of describing relevant processes while maintaining a manageable demand for input data and computation time. WASA-SED addresses this challenge by employing an innovative catena-based upscaling approach: the landscape is represented by characteristic toposequences. For deriving these toposequences with regard to multiple attributes (eg. topography, soils, vegetation) the LUMP-algorithm (Landscape Unit Mapping Program) was developed and related to objective 2. It incorporates an algorithm to retrieve representative catenas and their attributes, based on a Digital Elevation Model and supplemental spatial data. These catenas are classified to provide the discretization for the WASA-SED model. For objective 3, water and sediment fluxes were monitored at the catchment outlet of the Isábena and some of its sub-catchments. For sediment yield estimation, the intermittent measurements of suspended sediment concentration (SSC) had to be interpolated. This thesis presents a comparison of traditional sediment rating curves (SRCs), generalized linear models (GLMs) and non-parametric regression using Random Forests (RF) and Quantile Regression Forests (QRF). The observed SSCs are highly variable and range over six orders of magnitude. For these data, traditional SRCs performed poorly, as did GLMs, despite including other relevant process variables (e.g. rainfall intensities, discharge characteristics). RF and QRF proved to be very robust and performed favourably for reproducing sediment dynamics. QRF additionally excels in providing estimates on the accuracy of the predictions. Subsequent analysis showed that most of the sediment was exported during intense storms of late summer. Later floods yielded successively less sediment. Comparing sediment generation to yield at the outlet suggested considerable storage effects within the river channel. Addressing objective 4, the WASA-SED model was parameterized for the two study areas in NE Spain and applied with different foci. For Ribera Salada, the uncalibrated model yielded reasonable results for runoff and sediment. It provided quantitative measures of the change in runoff and sediment yield for different land-uses. Additional land management scenarios were presented and compared to impacts caused by climate change projections. In contrast, the application for the Isábena focussed on exploring the full potential of the model's predictive capabilities. The calibrated model achieved an acceptable performance for the validation period in terms of water and sediment fluxes. The inadequate representation of the lower sub-catchments inflicted considerable reductions on model performance, while results for the headwater catchments showed good agreement despite stark contrasts in sediment yield. In summary, the application of WASA-SED to three catchments proved the model framework to be a practicable multi-scale approach. It successfully links the hillslope to the catchment scale and integrates the three components hillslope, river and reservoir in one model. Thus, it provides a feasible approach for tackling issues of water and sediment yield at the meso-scale. The crucial role of processes like transmission losses and sediment storage in the river has been identified. Further advances can be expected when the representation of connectivity of water and sediment fluxes (intra-hillslope, hillslope-river, intra-river) is refined and input data improves.
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 Isábena 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 Isábena 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 Isábena 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 Isábena 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.
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 non-invasive tool for measuring hydrogen pools such as soil moisture, snow or vegetation. The intrinsic integration over a radial hectare-scale footprint is a clear advantage for averaging out small-scale heterogeneity, but on the other hand the data may become hard to interpret in complex terrain with patchy land use.
This study presents a directional shielding approach to prevent neutrons from certain angles from being counted while counting neutrons entering the detector from other angles and explores its potential to gain a sharper horizontal view on the surrounding soil moisture distribution.
Using the Monte Carlo code URANOS (Ultra Rapid Neutron-Only Simulation), we modelled the effect of additional polyethylene shields on the horizontal field of view and assessed its impact on the epithermal count rate, propagated uncertainties and aggregation time.
The results demonstrate that directional CRNS measurements are strongly dominated by isotropic neutron transport, which dilutes the signal of the targeted direction especially from the far field. For typical count rates of customary CRNS stations, directional shielding of half-spaces could not lead to acceptable precision at a daily time resolution. However, the mere statistical distinction of two rates should be feasible.
Cosmic-ray neutron sensing (CRNS) is a non-invasive tool for measuring hydrogen pools such as soil moisture, snow or vegetation. The intrinsic integration over a radial hectare-scale footprint is a clear advantage for averaging out small-scale heterogeneity, but on the other hand the data may become hard to interpret in complex terrain with patchy land use.
This study presents a directional shielding approach to prevent neutrons from certain angles from being counted while counting neutrons entering the detector from other angles and explores its potential to gain a sharper horizontal view on the surrounding soil moisture distribution.
Using the Monte Carlo code URANOS (Ultra Rapid Neutron-Only Simulation), we modelled the effect of additional polyethylene shields on the horizontal field of view and assessed its impact on the epithermal count rate, propagated uncertainties and aggregation time.
The results demonstrate that directional CRNS measurements are strongly dominated by isotropic neutron transport, which dilutes the signal of the targeted direction especially from the far field. For typical count rates of customary CRNS stations, directional shielding of half-spaces could not lead to acceptable precision at a daily time resolution. However, the mere statistical distinction of two rates should be feasible.
The Isabena catchment (445 km(2)), Spain, features highly diverse spatial heterogeneity in land use, lithology and rainfall. Consequently, the relative contribution in terms of water and sediment yield varies immensely between its subcatchments, and also temporally. This study presents the synthesis of similar to 2.5 years of monitoring rainfall, discharge and suspended sediment concentration (SSC) in the five main subcatchments of the Isabena and its outlet.
Continuous discharge at the subcatchment outlets, nine tipping bucket rainfall and automatic SSC samplers (complemented by manual samples), were collected from June 2011 until November 2013. The water stage records were converted to discharge using a rating curve derived with Bayesian regression. For reconstructing sediment yields, the data from the intermittent SSC sampling needed to be interpolated. We employed non-parametric multivariate regression (Quantile Regression Forests, QRF) using the discharge and rainfall data plus different aggregation levels of these as ancillary predictors. The subsequent Monte Carlo simulations allowed the determination of monthly sediment yields and their uncertainty.
The Isabena catchment shows high erosion dynamics with great variability in space and time, with stark contrasts even between adjacent subcatchments. The natural conditions make water and sediment monitoring and instrumentation very challenging; the measurement of discharge is particularly prone to considerable uncertainties. The QRF method employed for reconstructing sedigraphs and monthly yields proved well suited for the task.
This study focuses on the prediction of event-based runoff coefficients (an important descriptor of flood events) for nested catchments up to an area of 50?km(2) in the Eastern Ore Mountains. The four main objectives of the study are (i) the prediction of runoff coefficients with the statistical method of generalized linear models, (ii) the comparison of the results of the linear models with estimates of a distributed conceptual model, (iii) the comparison of the dynamics of observed soil moisture and simulated saturation deficit of the hydrological model and (iv) the analysis of the relationship between runoff coefficient and observed and simulated wetness. Different predictor variables were selected to describe the runoff coefficient and were differentiated into variables describing the catchment’s antecedent wetness and meteorological forcing. The best statistical model was estimated in a stepwise approach on the basis of hierarchical partitioning, an exhaustive search algorithm and model validation with jackknifing. We then applied the rainfall runoff model WaSiM ETH to predict the runoff processes for the two larger catchments. Locally measured small-scale soil moisture (acquired at a scale of four to five magnitudes smaller than the catchment) was identified as one of the key predictor variables for the estimation of the runoff coefficient with the general linear model. It was found that the relationship betweenobserved and simulated (using WaSiM ETH) wetness is strongly hysteretic. The runoff coefficients derived from the rainfall runoff simulations systematically underestimate the observed values. Copyright (C) 2012 John Wiley & Sons, Ltd.
The 933 km(2) Bengue catchment in northeastern Brazil is characterized by distinct rainy and dry seasons. Precipitation is stored in variously sized reservoirs, which is essential for the local population. In this study, we used TerraSAR-X SM(HH) data for an one-year monitoring of seasonal changes in the reservoir areas from July 2011 to July 2012. The monitoring was based on acquisitions in the ascending pass direction, complemented by occasional descending-pass images. To detect water surface areas, a histogram analysis followed by a global threshold classification was performed, and the results were validated using in situ GPS data. Distinguishing between small reservoirs and similar looking dark areas was difficult. Therefore, we tested several approaches for identifying misclassified areas. An analysis of the surface area dynamics of the reservoirs indicated high spatial and temporal heterogeneities and a large decrease in the total water surface area of the reservoirs in the catchment by approximately 30% within one year.
Cosmic-ray neutron sensing (CRNS) has become an effective method to measure soil moisture at a horizontal scale of hundreds of metres and a depth of decimetres. Recent studies proposed operating CRNS in a network with overlapping footprints in order to cover root-zone water dynamics at the small catchment scale and, at the same time, to represent spatial heterogeneity. In a joint field campaign from September to November 2020 (JFC-2020), five German research institutions deployed 15 CRNS sensors in the 0.4 km2 Wüstebach catchment (Eifel mountains, Germany). The catchment is dominantly forested (but includes a substantial fraction of open vegetation) and features a topographically distinct catchment boundary. In addition to the dense CRNS coverage, the campaign featured a unique combination of additional instruments and techniques: hydro-gravimetry (to detect water storage dynamics also below the root zone); ground-based and, for the first time, airborne CRNS roving; an extensive wireless soil sensor network, supplemented by manual measurements; and six weighable lysimeters. Together with comprehensive data from the long-term local research infrastructure, the published data set (available at https://doi.org/10.23728/b2share.756ca0485800474e9dc7f5949c63b872; Heistermann et al., 2022) will be a valuable asset in various research contexts: to advance the retrieval of landscape water storage from CRNS, wireless soil sensor networks, or hydrogravimetry; to identify scale-specific combinations of sensors and methods to represent soil moisture variability; to improve the understanding and simulation of land–atmosphere exchange as well as hydrological and hydrogeological processes at the hillslope and the catchment scale; and to support the retrieval of soil water content from airborne and spaceborne remote sensing platforms.
Cosmic-ray neutron sensing (CRNS) has become an effective method to measure soil moisture at a horizontal scale of hundreds of metres and a depth of decimetres. Recent studies proposed operating CRNS in a network with overlapping footprints in order to cover root-zone water dynamics at the small catchment scale and, at the same time, to represent spatial heterogeneity. In a joint field campaign from September to November 2020 (JFC-2020), five German research institutions deployed 15 CRNS sensors in the 0.4 km2 Wüstebach catchment (Eifel mountains, Germany). The catchment is dominantly forested (but includes a substantial fraction of open vegetation) and features a topographically distinct catchment boundary. In addition to the dense CRNS coverage, the campaign featured a unique combination of additional instruments and techniques: hydro-gravimetry (to detect water storage dynamics also below the root zone); ground-based and, for the first time, airborne CRNS roving; an extensive wireless soil sensor network, supplemented by manual measurements; and six weighable lysimeters. Together with comprehensive data from the long-term local research infrastructure, the published data set (available at https://doi.org/10.23728/b2share.756ca0485800474e9dc7f5949c63b872; Heistermann et al., 2022) will be a valuable asset in various research contexts: to advance the retrieval of landscape water storage from CRNS, wireless soil sensor networks, or hydrogravimetry; to identify scale-specific combinations of sensors and methods to represent soil moisture variability; to improve the understanding and simulation of land–atmosphere exchange as well as hydrological and hydrogeological processes at the hillslope and the catchment scale; and to support the retrieval of soil water content from airborne and spaceborne remote sensing platforms.
In a study from 2008, Lariviere and colleagues showed, for the field of natural sciences and engineering, that the median age of cited references is increasing over time. This result was considered counterintuitive: with the advent of electronic search engines, online journal issues and open access publications, one could have expected that cited literature is becoming younger. That study has motivated us to take a closer look at the changes in the age distribution of references that have been cited in water resources journals since 1965. Not only could we confirm the findings of Lariviere and colleagues. We were also able to show that the aging is mainly happening in the oldest 10-25% of an average reference list. This is consistent with our analysis of top-cited papers in the field of water resources. Rankings based on total citations since 1965 consistently show the dominance of old literature, including text books and research papers in equal shares. For most top-cited old-timers, citations are still growing exponentially. There is strong evidence that most citations are attracted by publications that introduced methods which meanwhile belong to the standard toolset of researchers and practitioners in the field of water resources. Although we think that this trend should not be overinterpreted as a sign of stagnancy, there might be cause for concern regarding how authors select their references. We question the increasing citation of textbook knowledge as it holds the risk that reference lists become overcrowded, and that the readability of papers deteriorates.
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.
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 results of streamflow trend studies are often characterized by mostly insignificant trends and inexplicable spatial patterns. In our study region, Western Austria, this applies especially for trends of annually averaged runoff. However, analysing the altitudinal aspect, we found that there is a trend gradient from higher-altitude to lower-altitude stations, i.e. a pattern of mostly positive annual trends at higher stations and negative ones at lower stations. At midaltitudes, the trends are mostly insignificant. Here we hypothesize that the streamflow trends are caused by the following two main processes: on the one hand, melting glaciers produce excess runoff at higher-altitude watersheds. On the other hand, rising temperatures potentially alter hydrological conditions in terms of less snowfall, higher infiltration, enhanced evapotranspiration, etc., which in turn results in decreasing streamflow trends at lower-altitude watersheds. However, these patterns are masked at mid-altitudes because the resulting positive and negative trends balance each other. To support these hypotheses, we attempted to attribute the detected trends to specific causes. For this purpose, we analysed trends of filtered daily streamflow data, as the causes for these changes might be restricted to a smaller temporal scale than the annual one. This allowed for the explicit determination of the exact days of year (DOYs) when certain streamflow trends emerge, which were then linked with the corresponding DOYs of the trends and characteristic dates of other observed variables, e.g. the average DOY when temperature crosses the freezing point in spring. Based on these analyses, an empirical statistical model was derived that was able to simulate daily streamflow trends sufficiently well. Analyses of subdaily streamflow changes provided additional insights. Finally, the present study supports many modelling approaches in the literature which found out that the main drivers of alpine streamflow changes are increased glacial melt, earlier snowmelt and lower snow accumulation in wintertime.
The results of streamflow trend studies are often characterized by mostly insignificant trends and inexplicable spatial patterns. In our study region, Western Austria, this applies especially for trends of annually averaged runoff. However, analysing the altitudinal aspect, we found that there is a trend gradient from higher-altitude to lower-altitude stations, i.e. a pattern of mostly positive annual trends at higher stations and negative ones at lower stations. At midaltitudes, the trends are mostly insignificant. Here we hypothesize that the streamflow trends are caused by the following two main processes: on the one hand, melting glaciers produce excess runoff at higher-altitude watersheds. On the other hand, rising temperatures potentially alter hydrological conditions in terms of less snowfall, higher infiltration, enhanced evapotranspiration, etc., which in turn results in decreasing streamflow trends at lower-altitude watersheds. However, these patterns are masked at mid-altitudes because the resulting positive and negative trends balance each other. To support these hypotheses, we attempted to attribute the detected trends to specific causes. For this purpose, we analysed trends of filtered daily streamflow data, as the causes for these changes might be restricted to a smaller temporal scale than the annual one. This allowed for the explicit determination of the exact days of year (DOYs) when certain streamflow trends emerge, which were then linked with the corresponding DOYs of the trends and characteristic dates of other observed variables, e.g. the average DOY when temperature crosses the freezing point in spring. Based on these analyses, an empirical statistical model was derived that was able to simulate daily streamflow trends sufficiently well. Analyses of subdaily streamflow changes provided additional insights. Finally, the present study supports many modelling approaches in the literature which found out that the main drivers of alpine streamflow changes are increased glacial melt, earlier snowmelt and lower snow accumulation in wintertime.