Refine
Year of publication
Document Type
- Scientific journal article peer-reviewed (25) (remove)
Way of publication
- Open Access (2)
Keywords
- runoff generation (3)
- Diffusive-wave equation (2)
- Runoff generation (2)
- Shallow-water equations (2)
- Surface runoff (2)
- Zero-inertia equation (2)
- deep drainage (2)
- microtopography (2)
- rainfall variability (2)
- rainfall/runoff simulation (2)
Institute
Reactive transport controls on sandy acid sulfate soils and impacts on shallow groundwater quality
(2014)
Disturbance or drainage of potential acid sulfate soils (PASS) can result in the release of acidity and degradation of infrastructure, water resources, and the environment. Soil processes affecting shallow groundwater quality have been investigated using a numerical code that integrates (bio)geochemical processes with water, solute, and gas transport. The patterns of severe and persistent acidification (pH < 4) in the sandy, carbonate-depleted podzols of a coastal plain could be reproduced without calibration, based on oxidation of microcrystalline pyrite after groundwater level decrease and/or residual groundwater acidity, due to slow vertical solute transport rates. The rate of acidification was limited by gas phase diffusion of oxygen and hence was sensitive to soil water retention properties and in some cases also to oxygen con-sumption by organic matter mineralization. Despite diffusion limitation, the rate of oxidation in sandy soils was rapid once pyrite-bearing horizons were exposed, even to a depth of 7.5 m. Groundwater level movement was thus identified as an important control on acidification, as well as the initial pyrite content. Increase in the rate of Fe(II) oxidation lead to slightly lower pH and greater accumulation of Fe(III) phases, but had little effect on the overall amount of pyrite oxidized. Aluminosilicate (kaolinite) dissolution had a small pH-buffering effect but lead to the release of Al and associated acidity. Simulated dewatering scenarios highlighted the potential of the model for risk assessment of (bio)geochemical impacts on soil and groundwater over a range of temporal and spatial scales.
Flash droughts are caused by a rapid depletion of soil moisture, and they severely affect vegetation growth and agricultural production. Notwithstanding the growing importance of flash droughts under the warming climate, drivers of flash droughts across the Europe are not well understood. Here we estimate the changes in flash droughts characteristics across Europe using the ERA5 reanalysis dataset for 1950–2019 period. We find a substantial increase in the frequency and spatial extent of flash droughts across Europe (with 79% of the total area) during the growing season with at-least one fourth of domain showing two-fold increase in the recent decades. Increased occurrence of flash drought is largely attributed to frequent occurrence of warmer and drier compound extremes, with a sharp gradient of changes being noticed in Mediterranean and Central European regions. Compound meteorological extremes causing the flash drought events across Europe are pre-dominantly driven by the recent climate warming. With unabated greenhouse gas emissions and current pace of climate warming, Europe is likely to face an increased occurrence of flash droughts, requiring prompt response for effective drought adaptation and management strategies.
Rainfall is a spatiotemporally variated process and one of the key elements to accurately capture both catchment runoff response and floodplain extents. Flash floods are the result of intense rainfall, typically associated to highly variable rain in both space and time, such as convective storms. In this work, the extent within-storm variability affects runoff and flooding is explored. The Kan catchment (Tehran, Iran) is used as base topography for the simulations. The allochthonous nature of floods in the catchment and how they interact with the effects of storm variability are further investigated. For this, 300 synthetic rainfall signals with different hyetograph variabilities are generated and imposed on a 1D/2D hydrodynamic model. Additionally, a set of simulations with different levels of spatial variability are performed. The results suggest that temporal and spatial variability affect the runoff response in different degrees. Peak discharge and hydrograph shapes, as well as flooded areas, are affected. The effect of storm temporal variability is shown to be significantly higher than storm spatial variability and storm properties such as return period, duration, and volume. Further on the influence of storm spatiotemporal variability on stream discharge and flood response is seen to be strongly dependent on the location within the drainage network at which it is assessed.
Rainfall is a spatiotemporally varied process and key to accurately capturing catchment runoff and determining flood response. Flash flood response of a catchment can be strongly governed by a rainfall’s spatiotemporal variability and is influenced by storm movement which drives a continuous spatiotemporal change throughout a rainfall event. In this work, the sensitivity of runoff and flooded areas to rainfall movement are assessed in the Kan catchment (Iran). The allochthonous nature of floods in the catchment and how they interact with the effects of rainfall movement are investigated. Fifty synthetic rain hyetographs are generated and traversed over the catchment under different velocities and directions and used to force a 1D/2D hydrodynamic model. The results suggest rainfall movement affects the runoff response in different degrees. Peak discharge, hydrograph shapes and flooded areas are affected. Storms with higher velocities result in higher peaks and faster onsets of runoff and consequently higher flooded areas in comparison to slower storms. The direction of the movement also plays a role. Storms moving along the average direction of the stream result in higher peaks and flooded areas. The relevance of storm direction is greater for slow moving storms. Additionally, the influence of rainfall movement is modulated by hyetograph structure, and the allochthonous behavior is greatly dependent on the location within the drainage network at which it is assessed.
Runoff generation from rainfall events is a complex, spatial and temporally dependent process strongly governed, among other factors, by catchment surface topography. Although it is widely known that many catchments experience morphological evolution, it is often ignored in analysis for different reasons ranging from simplification to lack of data. However, young catchments and early landscapes (such as those which are affected by natural or anthropogenic disturbances) do exhibit topography changes which in turn affect catchment hydrodynamics, hydrology and in particular runoff. In this work, we study the runoff generation and hydrodynamics of the Hühnerwasser artificial catchment (Brandenburg, Germany) during a period of erosion-based topographical changes (2006–2010). Nine Digital Elevation Models from such period were used as topography over which physically-based simulations were performed. The results suggest that topographic evolution in this catchment mostly affects the onset of runoff, whereas peak discharges and receding hydrograph limbs are less affected. These differences in hydrological signatures can be explained through the changes in the spatial distribution of runoff hydrodynamics and their impact on surface runoff connectivity. Relatively small topographical differences produce changing ponding conditions and modify flowpaths which becomes evident only through inspection of the spatial distribution of hydrodynamic variables. Moreover, the study shows that in order for simulations to be able to capture such responses, appropriate computational mesh and topographical data resolution are critical, since connectivity itself can be greatly affected by low resolution data or representation.
Surface flows of hydrological interest, including overland flow, runoff, river and channel flow and flooding have received significant attention from modellers in the past 30 years. A growing effort to address these complex environmental problems is in place in the scientific community. Researchers have stud-ied and favoured a plethora of techniques to approach this issue, ranging from very simple empirically-based mathematical models, to physically-based, deductive and very formal numerical integration of systems of partial-differential equations. In this work, we review two families of methods: cell-based simulators – later called Cellular Automata – and Finite Volume solvers for the Zero-Inertia equation, which we show to converge into a single methodology given appropriate choices. Furthermore, this convergence, mathematically shown in this work, can also be identified by critically reviewing the exist-ing literature, which leads to the conclusion that two methods originating from different reasoning and fundamental philosophy, fundamentally converge into the same method. Moreover, acknowledging such convergence allows for some generalisation of properties of numerical schemes such as error behaviour and stability, which, importantly, is the same for the converging methodology, a fact with practical implications. Both the review of existing literature and reasoning in this work attempts to aid in the effort of synchronising and cross-fertilizing efforts to improve the understanding and the outlook of Zero-Inertia solvers for surface flows, as well as to help in clarifying the possible confusion and parallel develop-ments that may arise from the use of different terminology originating from historical reasons. Moreover, synchronising and unifying this knowledge-base can help clarify model capabilities, applicability and modelling issues for hydrological modellers, specially for those not deeply familiar with the mathematical and numerical details.
To simulate the impacts of within-storm rainfall variabilities on fast hydrological processes, long precipi-tation time series with high temporal resolution are required. Due to limited availability of observed data such time series are typically obtained from stochastic models. However, most existing rainfall models are limited in their ability to conserve rainfall event statistics which are relevant for hydrological pro-cesses. Poisson rectangular pulse models are widely applied to generate long time series of alternating precipitation events durations and mean intensities as well as interstorm period durations. Multiplicative microcanonical random cascade (MRC) models are used to disaggregate precipitation time series from coarse to fine temporal resolution.
To overcome the inconsistencies between the temporal structure of the Poisson rectangular pulse model and the MRC model, we developed a new coupling approach by introducing two modifications to the MRC model. These modifications comprise (a) a modified cascade model (“constrained cascade”) which preserves the event durations generated by the Poisson rectangular model by constraining the first and last interval of a precipitation event to contain precipitation and (b) continuous sigmoid func-tions of the multiplicative weights to consider the scale-dependency in the disaggregation of precipita-tion events of different durations. The constrained cascade model was evaluated in its ability to dis-aggregate observed precipitation events in comparison to existing MRC models. For that, we used a 20-year record of hourly precipitation at six stations across Germany. The constrained cascade model showed a pronounced better agreement with the observed data in terms of both the temporal pattern of the precipitation time series (e.g. the dry and wet spell durations and autocorrelations) and event char-acteristics (e.g. intra-event intermittency and intensity fluctuation within events). The constrained cas-cade model also slightly outperformed the other MRC models with respect to the intensity-frequency relationship.
To assess the performance of the coupled Poisson rectangular pulse and constrained cascade model, precipitation events were stochastically generated by the Poisson rectangular pulse model and then disaggregated by the constrained cascade model. We found that the coupled model performs satisfac-torily in terms of the temporal pattern of the precipitation time series, event characteristics and the in-tensity-frequency relationship.
Governed by a series of non-linear feedback mechanisms among water, vegetation and decomposers, carbon storage within wetlands is important on a global scale. However, the effect that climatic fluctuations have on those mechanisms is not well documented. In this study, we introduce a mechanistic model connecting hydrology, vegetation and microbial biomass to investigate how changes in the climate signal propagate through wetland ecosystems, via vegetation and microbial dynamics, and attempt to quantify how net rates of wetland carbon metabolism change in response to a changing climate. Our particular focus is the dryland–wetland systems found in south-west Western Australia (SWWA), as they are expected to be sensitive to projected climatic changes due to their close linkage to the seasonal water delivery pattern. The model simulations investigate wetland carbon retention under different hydro-climatological conditions ranging across a regional gradient in the dryness index. The results indicate that short term and long term vegetation responses may be counter-intuitive due to adaptability in the water uptake strategy of the vegetation community partially decoupling biomass from water availability. Furthermore, changes in water delivery are not a good indicator for overall changes in wetland metabolism, defined as the net rate of carbon assimilation, due to the dominance of the soil carbon storages and their sensitivity to heightened bacterial metabolism rates with increasing temperatures. The results highlight that an optimum combination of water supply and vegetation leads to a higher percentage of carbon being stored in soils, therefore increasing the resistance of the carbon storage to changes in precipitation. The model presented here provides a first step to explain how changing patterns of rainfall, temperature and evapotranspiration can change carbon cycling characteristics and the carbon retention efficiency of dryland–wetlands.
Transport processes are often proposed to explain the sorting of rock particles on arid hillslopes, where mean rock particle size often decreases in the downslope direction. Here we show that in situ fragmentation of rock particles can also produce similar patterns. A total of 93,414 rock particles were digitized from 880 photographs of the surface of three mesa hills in the Great Sandy Desert, Australia. Rock particles were characterized by the projected Feret's diameter and circularity. Distance from the duricrust cap was found to be a more robust explanatory variable for diameter than the local hillslope gradient. Mean diameter decreased exponentially downslope, while the fractional area covered by rock particles decreased linearly. Rock particle diameters were distributed lognormally, with both the location and scale parameters decreasing approximately linearly downslope. Rock particle circularity distributions showed little change; only a slight shift in the mode to more circular particles was noted to occur downslope. A dynamic fragmentation model was used to assess whether in situ weathering alone could reproduce the observed downslope fining of diameters. Modeled and observed size distributions agreed well and both displayed a preferential loss of relatively large rock particles and an apparent approach to a terminal size distribution of the rocks downslope. We show this is consistent with a size effect in material strength, where large rocks are more susceptible to fatigue failure under stress than smaller rocks. In situ fragmentation therefore produces qualitatively similar patterns to those that would be expected to arise from selective transport.
Microtopography (MT) can govern runoff dynamics as a net result of local heterogeneities in the flow paths and ponding. This in turn controls the development of the surface water layer that connects and flows downslope. It is therefore important to understand which microtopographic features affect runoff generation dynamics and its macroscopic—hillslope scale—hydrological signatures (e.g., hydrographs, runoff and infiltration volumes). In this study, we numerically solve 2D overland flow from a single rain pulse on 1,460 idealized hillslopes with different slopes and sinusoidal microtopographies and different infiltration capacities. We assess hydrodynamic distributions, hydrographs and hydrological indices to assess the effects of MT and infiltration on the (local) hydrodynamic and (larger scale) hydrologic responses in terms of surface runoff regimes. The results show that MT enhances infiltration and that infiltration and runoff depend in a strong non-linear way on slope and the properties of MT. Three regimes of influence of MT were identified: one in which MT plays a negligible role but there is a high sensitivity to the infiltration capacity curve, a second regime in which hydrological partitioning is highly sensitive to MT and the infiltration capacity curve, and a third regime in which MT increases infiltration, but the response is insensitive to particular features, and more affected by the average slopes. The regimes are the product of the interplay between small (MT) and large scale (slope) properties. Furthermore, the results suggest that hydrological signatures can be interpreted and explained by the spatiotemporal variation of surface connectivity.