Catchment-scale Richards equation-based modeling of evapotranspiration via boundary condition switching and root water uptake schemes

In arid and semiarid climate catchments, where annual evapotranspiration (ET) and rainfall are typically comparable, modeling ET is important for proper assessment of water availability and sustainable land use management. The aim of the present study is to assess different parsimonious schemes for representing ET in a process-based model of coupled surface and subsurface flow. A simplified method for computing ET based on a switching procedure for the boundary conditions of the Richards equation at the soil surface is compared to a sink term approach that includes root water uptake, root distribution, root water compensation, and water and oxygen stress. The study site for the analysis is a small pasture catchment in southeastern Australia. A comprehensive sensitivity analysis carried out on the parameters of the sink term shows that the maximum root depth is the dominant control on catchment-scale ET and streamflow. Comparison with the boundary condition switching method demonstrates that this simpler scheme (only one parameter) can successfully reproduce ET when the vegetation root depth is shallow (not exceeding approximately 50 cm). For deeper rooting systems, the switching scheme fails to match the ET fluxes and is affected by numerical artifacts, generating physically unrealistic soil moisture dynamics. It is further shown that when transpiration is the dominant contribution to ET, the inclusion of oxygen stress and root water compensation in the model can have a considerable effect on the estimation of both ET and streamflow; this is mostly due to the water fluxes associated with the riparian zone. Key Points: Simple, parsimonious ET schemes for integrated hydrological models are assessed Boundary condition switching is suitable only for shallow root depths Oxygen stress and root water compensation influence riparian zone ET dynamics

Newtonian nudging for a Richards equation-based distributed model

This report describes a series of simulations conducted with a hydrological model, CATHY, to test a recently implemented data assimilation technique, Newtonian nudging

Assessment and formulation of data assimilation techniques for a 3D Richards equation-based hydrological model

The main objectives of de DAUFIN project are: to develop a unifying modeling framework applicable at the catchment scale and based on rigorous conservation equations for the study of hydrological processes in the stream channel, land surface, soil, and groundwater components of a river basin; to implement data assimilation methodologies within this modeling framework and for other control models to enable the optimal use of remote sensing, ground-based, and simulation data; to test and apply the models and the data assimilation methods at various catchment scales, including hillslopes and subcatchment of the Ourthe water shed in Belgium and the entire Meuse river basin, one of the major basins in Europe with a drainage area of 33000 km² that comprises the Ourthe

Literature review on NAPL contamination and remediation

Remediation of polluted soils and groundwater is of major concern due to the increasing number of contaminated aquifers. Subsurface aquifers constitute one of the most important sources of drinkable water. In recent years, water needs have been increasing due to increases in development and human population. Several sorts of contaminants can be found in groundwater: metal ions, pesticides, aliphatic and aromatic hydrocarbons, polycyclic hydrocarbons, chlorinated hydrocarbons, etc. The toxicity of these compounds varies and so do guidelines that establish allowable concentration levels in drinking water. Among the aforementioned types of compounds, a particular importance is assumed by those which exist as a separate phase when their concentrations in water exceed a certain limit. The transport behavior and dynamics of multiphase contaminants are very different from their dissolved counterparts, and are very difficult both to describe and to model. Several phenomena can take place, such as organic phase trapping, formation of ganglia and pools of contaminant, sorption, hysteresis in both soil imbibition and drainage, capillarity, fingering, and mass-transfer. In such cases, our ability to describe and predict the fate of a contaminant plume in which a separate organic phase occurs is limited, and research within this field is quite open. Much effort has been devoted in trying to describe the characteristics of the phenomena occuring in multiphase systems, and several models and formulations have been proposed for predicting the fate of contaminants when present in such systems (see Miller et al. 1997) for a review on multiphase modeling in porous media). Work has also been done for modeling human intervention techniques for containing and/or reducing soil contaminantion (NRC, 1994), such as pumping, clean water-air-steam injection, soil heating, surfactants, biological methods, etc. Finally, much work has also been done on the numerical solution of mathematical models whose complexity does not allow for an analytical solution. Among the dozens of remediation methods which have been proposed and which are strongly dependent on site environmental conditions, biological methods are achieving increasing importance, due to their “naturalness" and their low costs (NRC, 1993) . It has been noticed that soil microorganisms are able to degrade several classes of compounds, in particular those which partition between an aqueous and an organic phase, or sometimes also gaseous phase, for e.g., hydrocarbons, chlorinated compounds, pesticides. These compounds, or better said, their fractions dissolved in water, are liable to be metabolized by subsurface microrganisms which have the capability to degrade the compounds and to transform them into carbon dioxide and/or other compounds, which are less toxic or unnoxious. Several laboratory and field studies have been conducted for assessing and evaluating the capability and the limits of soil microorganisms to degrade several classes of contaminants (Mayer et al., 1994, 1995, 1996, 1997) . Much work has also been devoted to modeling biodegration of groundwater contaminants. The outline of this report is as follows: section 2 gives a brief description of the characteristics and properties of NAPLs, including a review of the literature with regards to formulations and modeling; section 3 discusses biodegradation of contaminants and past efforts at modeling biodegradation; section 4 surveys specific remediation technologies and experiences; and section 5 discusses open issues for further research. In the final section possible lines of research for the second phase of the PhD program are indicated

Physically based modeling in catchment hydrology at 50: Survey and outlook

Integrated, process-based numerical models in hydrology are rapidly evolving, spurred by novel theories in mathematical physics, advances in computational methods, insights from laboratory and field experiments, and the need to better understand and predict the potential impacts of population, land use, and climate change on our water resources. At the catchment scale, these simulation models are commonly based on conservation principles for surface and subsurface water flow and solute transport (e.g., the Richards, shallow water, and advection-dispersion equations), and they require robust numerical techniques for their resolution. Traditional (and still open) challenges in developing reliable and efficient models are associated with heterogeneity and variability in parameters and state variables; nonlinearities and scale effects in process dynamics; and complex or poorly known boundary conditions and initial system states. As catchment modeling enters a highly interdisciplinary era, new challenges arise from the need to maintain physical and numerical consistency in the description of multiple processes that interact over a range of scales and across different compartments of an overall system. This paper first gives an historical overview (past 50 years) of some of the key developments in physically based hydrological modeling, emphasizing how the interplay between theory, experiments, and modeling has contributed to advancing the state of the art. The second part of the paper examines some outstanding problems in integrated catchment modeling from the perspective of recent developments in mathematical and computational science

Examination of the seepage face boundary condition in subsurface and coupled surface/subsurface hydrological models

A seepage face is a nonlinear dynamic boundary that strongly affects pressure head distributions, water table fluctuations, and flow patterns. Its handling in hydrological models, especially under complex conditions such as heterogeneity and coupled surface/subsurface flow, has not been extensively studied. In this paper, we compare the treatment of the seepage face as a static (Dirichlet) versus dynamic boundary condition, we assess its resolution under conditions of layered heterogeneity, we examine its interaction with a catchment outlet boundary, and we investigate the effects of surface/subsurface exchanges on seepage faces forming at the land surface. The analyses are carried out with an integrated catchment hydrological model. Numerical simulations are performed for a synthetic rectangular sloping aquifer and for an experimental hillslope from the Landscape Evolution Observatory. The results show that the static boundary condition is not always an adequate stand-in for a dynamic seepage face boundary condition, especially under conditions of high rainfall, steep slope, or heterogeneity; that hillslopes with layered heterogeneity give rise to multiple seepage faces that can be highly dynamic; that seepage face and outlet boundaries can coexist in an integrated hydrological model and both play an important role; and that seepage faces at the land surface are not always controlled by subsurface flow. The paper also presents a generalized algorithm for resolving seepage face outflow that handles heterogeneity in a simple way, is applicable to unstructured grids, and is shown experimentally to be equivalent to the treatment of atmospheric boundary conditions in subsurface flow models

Implementation of a catchment hydrologic model for the Brisy subcatchment of the Ourthe watershed, and generation of a dataset for a 240-day storm-interstorm sequence

This report describes the generation of a synthetic dataset needed for testing and verification of more simplified modeling approaches which aim to develop models applicable at large catchment and river basin scales. The work is carried out within the framework of a European project (DAUFIN) on developing data assimilation methodologies and a unified framework for hydrological modeling of catchment and river basin flow processes

Rapporto di ricerca bibliografica (Stato dell'arte dei modelli di flusso e trasporto in mezzi porosi)

In questo primo rapporto vengono presentate le equazioni fondamentali che reggono i fenomeni di flusso e trasporto in mezzi porosi insieme ai metodi numerici che verranno utilizzati per la soluzione di tali equazioni. I modelli matematici in questione sono basati su equazioni differenziali a derivate parziali che impongono il bilancio di massa sia per il fluido che per il soluto (inquinante disciolto in acqua). Queste equazioni vengono scritte in forma generale per un mezzo poroso tridimensionale; in dipendenza dal tipo di applicazione è possibile adottare modelli mono o bidimensionali che portano a semplificazioni notevoli. L'equazione di flusso è sviluppata per il caso di mezzi porosi a saturazione variabile e può essere quindi utilizzata contemporaneamente nella zona insatura (suoli superficiali) e satura (falde freatiche e artesiane). Nell'equazione di trasporto si considerano i processi di dispersione, diffusione e avvezione, insieme ad alcune fenomenologie di interazione chimico-fisica tra il soluto e la matrice porosa. Accanto a queste equazione, si descrive anche un modello, a scala di bacino, di afflussi-deflussi superficiali accoppiato con un modello di infiltrazione. Questo approccio viene tiene conto di fenomeni importanti qualora vi sia una stretta correlazione tra il moto dell'acqua in superficie e il moto dell'acqua nella zona insatura

Description of a hydrologic dataset for the Brisy subcatchment

This report describes the dataset for the Brisy subcatchment in south eastern Belgium, which is a subcatchment of the Ourthe catchment, itself a subcatchment of the Meuse river basin. The data preparation, organization, and processing steps undertaken for both the Meuse basin and the Brisy subcatchment will be detailed

Comparison of a 3-D Richards equation-based model with the hillslope-storage Boussinesq model: A test case for nine characteristic hillslopes

In this study, a numerical simulation for these nine characteristic hillslopes is carried out using the three-dimensional, distributed model CATHY: a coupled model describing both subsurface flow and surface runoff. The model is capable of accurately describing flow characteristic, and thus provides us with the opportunity of studying spatially distributed phenomena like water table or storage dynamics in detail. The objective of the study is to compare the CATHY simulations with a recently developed extensions to the Boussinesq model that uses the width function to account for hillslopes shape