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

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

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

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

Assessment of alternative land management practices using hydrological simulation and a decision support tool: Arborea agricultural region, Sardinia

Quantifying the impact of land use on water supply and quality is a primary focus of environmental management. In this work we apply a semidistributed hydrological model (SWAT) to predict the impact of different land management practices on water and agricultural chemical yield over a long period of time for a study site situated in the Arborea region of central Sardinia, Italy. The physical processes associated with water movement, crop growth, and nutrient cycling are directly modeled by SWAT. The model simulations are used to identify indicators that reflect critical processes related to the integrity and sustainability of the ecosystem. Specifically we focus on stream quality and quantity indicators associated with anthropogenic and natural sources of pollution. A multicriteria decision support system is then used to develop the analysis matrix where water quality and quantity indicators for the rivers, lagoons, and soil are combined with socio-economic variables. The DSS is used to assess four options involving alternative watersheds designated for intensive agriculture and dairy farming and the use or not of treated wastewater for irrigation. Our analysis suggests that of the four options, the most widely acceptable consists in the transfer of intensive agricultural practices to the larger watershed, which is less vulnerable, in tandem with wastewater reuse, which rates highly due to water scarcity in this region of the Mediterranean. More generally, the work demonstrates how both qualitative and quantitative methods and information can assist decision making in complex settings

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

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