ABSTRACT
Subsurface contamination due to excessive nutrient surpluses is a persistent and widespread problem in agricultural areas across Europe. The vulnerability of a particular location to pollution from reactive solutes, such as nitrate, is determined by the interplay between hydrologic transport and biogeochemical transformations. Current studies on the controls of subsurface vulnerability do not consider the transient behaviour of transport dynamics in the root zone. Here, using state-of-the-art hydrologic simulations driven by observed hydroclimatic forcing, we demonstrate the strong spatiotemporal heterogeneity of hydrologic transport dynamics and reveal that these dynamics are primarily controlled by the hydroclimatic gradient of the aridity index across Europe. Contrasting the space-time dynamics of transport times with reactive timescales of denitrification in soil indicate that ~75% of the cultivated areas across Europe are potentially vulnerable to nitrate leaching for at least one-third of the year. We find that neglecting the transient nature of transport and reaction timescale results in a great underestimation of the extent of vulnerable regions by almost 50%. Therefore, future vulnerability and risk assessment studies must account for the transient behaviour of transport and biogeochemical transformation processes.
ABSTRACT
We study the transport behavior of a passive scalar in a two-dimensional (2D) time-independent Gaussian random velocity field by efficient and highly accurate numerical simulations. The model under consideration has been used in order to gain basic understanding of transport processes in incompressible flow through heterogeneous porous media. The velocity field is derived from the linearized solution of the Darcy equation with a Gauss-distributed log-hydraulic conductivity. The transport of a passive scalar is studied by a high precision random-walk method, which allows for a systematic nonperturbative study of the ensemble and effective dispersion coefficients. The conclusive numerical results validate the range of applicability of the perturbation theory and the consistency of nonperturbative approaches to the transport problem in a random medium. Furthermore, we observe closed streamlines in incompressible 2D Gaussian random fields, which restricts the direct applicability of the simulation method for transport in heterogeneous porous media, and questions the results of similar studies that do not observe this phenomenon.
ABSTRACT
With the aid of integral transforms, analytical solutions for the transport of a decay chain in homogenous porous media are derived. Unidirectional steady-state flow and radial steady-state flow in single and multiple porosity media are considered. At least in Laplace domain, all solutions can be written in closed analytical formulae. Partly, the solutions can also be inverted analytically. If not, analytical calculation of the steady-state concentration distributions, evaluation of temporal moments and numerical inversion are still possible. Formulae for several simple boundary conditions are given and visualized in this paper. The derived novel solutions are widely applicable and are very useful for the validation of numerical transport codes.