RETRASO, a code for modeling reactive transport in saturated and unsaturated porous media

The code RETRASO (REactive TRAnsport of SOlutes) simulates reactive transport of dissolved and gaseous species in non-isothermal saturated or unsaturated problems. Possible chemical reactions include aqueous complexation (including redox reactions), sorption, precipitation-dissolution of minerals and gas dissolution. Various models for sorption of solutes on solids are available, from experimental relationships (linear KD, Freundlich and Langmuir isotherms) to cation exchange and surface complexation models (constant capacitance, diffuse layer and triple layer models). Precipitation-dissolution and aqueous complexation can be modelled in equilibrium or according to kinetic laws. For the numerical solution of the reactive transport equations it uses the Direct Substitution Approach. The use of the code is demonstrated by three examples. The first example models various sorption processes in a smectite barrier. The second example models a complex chemical system in a two dimensional cross-section. The last example models pyrite weathering in an unsaturated medium

Selection of a numerical unsaturated flow code for tilted capillary barrier performance evaluation

Capillary barriers consisting of tilted fine-over-coarse layers have been suggested as landfill covers as a means to divert water infiltration away from sensitive underground regions under unsaturated flow conditions, especially for arid and semi-arid regions. Typically, the HELP code is used to evaluate landfill cover performance and design. Unfortunately, due to its simplified treatment of unsaturated flow and its essentially one-dimensional nature, HELP is not adequate to treat the complex multidimensional unsaturated flow processes occurring in a tilted capillary barrier. In order to develop the necessary mechanistic code for the performance evaluation of tilted capillary barriers, an efficient and comprehensive unsaturated flow code needs to be selected for further use and modification. The present study evaluates a number of candidate mechanistic unsaturated flow codes for application to tilted capillary barriers. Factors considered included unsaturated flow modeling, inclusion of evapotranspiration, nodalization flexibility, ease of modification, and numerical efficiency. A number of unsaturated flow codes are available for use with different features and assumptions. The codes chosen for this evaluation are TOUGH2, FEHM, and SWMS{_}2D. All three codes chosen for this evaluation successfully simulated the capillary barrier problem chosen for the code comparison, although FEHM used a reduced grid. The numerical results are a strong function of the numerical weighting scheme. For the same weighting scheme, similar results were obtained from the various codes. Based on the CPU time of the various codes and the code capabilities, the TOUGH2 code has been selected as the appropriate code for tilted capillary barrier performance evaluation, possibly in conjunction with the infiltration, runoff, and evapotranspiration models of HELP. 44 refs

Multi-source data assimilation for physically based hydrological modeling of an experimental hillslope

Data assimilation has recently been the focus of much attention for integrated surface–subsurface hydrological models, whereby joint assimilation of water table, soil moisture, and river discharge measurements with the ensemble Kalman filter (EnKF) has been extensively applied. Although the EnKF has been specifically developed to deal with nonlinear models, integrated hydrological models based on the Richards equation still represent a challenge, due to strong nonlinearities that may significantly affect the filter performance. Thus, more studies are needed to investigate the capabilities of the EnKF to correct the system state and identify parameters in cases where the unsaturated zone dynamics are dominant, as well as to quantify possible tradeoffs associated with assimilation of multi-source data. Here, the CATHY (CATchment HYdrology) model is applied to reproduce the hydrological dynamics observed in an experimental two-layered hillslope, equipped with tensiometers, water content reflectometer probes, and tipping bucket flow gages to monitor the hillslope response to a series of artificial rainfall events. Pressure head, soil moisture, and subsurface outflow are assimilated with the EnKF in a number of scenarios and the challenges and issues arising from the assimilation of multi-source data in this real-world test case are discussed. Our results demonstrate that the EnKF is able to effectively correct states and parameters even in a real application characterized by strong nonlinearities. However, multi-source data assimilation may lead to significant tradeoffs: the assimilation of additional variables can lead to degradation of model predictions for other variables that are otherwise well reproduced. Furthermore, we show that integrated observations such as outflow discharge cannot compensate for the lack of well-distributed data in heterogeneous hillslopes.</p

Prediction of tilted capillary barrier performance

Capillary barriers, consisting of tilted fine-over-coarse layers under unsaturated conditions, have been suggested as landfill covers to divert water infiltration away from sensitive underground regions, especially for arid and semi-arid regions. The Hydrological Evaluation of Landfill Performance (HELP) computer code is an evaluation tool for landfill covers used by designers and regulators. HELP is a quasi-two-dimensional model that predicts moisture movement into and through the underground soil and waste layers. Processes modeled within HELP include precipitation, runoff, evapotranspiration, unsaturated vertical drainage, saturated lateral drainage, and leakage through liners. Unfortunately, multidimensional unsaturated flow phenomena that are necessary for evaluating tilted capillary barriers are not included in HELP. Differences between the predictions of the HELP and those from a multidimensional unsaturated flow code are presented to assess the two different approaches. Comparisons are presented for the landfill covers including capillary barrier configurations at the Alternative Landfill Cover Demonstration (ALCD) being conducted at Sandia

Comparison of Ross` capillary barrier diversion formula with detailed numerical simulations

Ross developed an analytical relationship to calculate the diversion length of a tilted fine-over-coarse capillary barrier. Oldenburg and Pruess compared simulation results using upstream and harmonic weighting to the diversion length predicted by Ross formula with mixed results; the qualitative agreement is reasonable but the quantitative comparison is poor, especially for upstream weighting. The proximity of the water table to the fine-coarse interface at breakthrough is a possible reason for the poor agreement. In the present study, the Oldenburg and Pruess problem is extended to address the water table issue. When the water table is sufficiently far away from the interface at breakthrough, good qualitative and quantitative agreement is obtained using upstream weighting

Prediction of water retention of soils from the humid tropics by the nonparametric k-nearest neighbor approach

Nonparametric approaches such as the k-nearest neighbor (k-NN) approach are considered attractive for pedotransfer modeling in hydrology; however, they have not been applied to predict water retention of highly weathered soils in the humid tropics. Therefore, the objectives of this study were: to apply the k-NN approach to predict soil water retention in a humid tropical region; to test its ability to predict soil water content at eight different matric potentials; to test the benefit of using more input attributes than most previous studies and their combinations; to discuss the importance of particular input attributes in the prediction of soil water retention at low, intermediate, and high matric potentials; and to compare this approach with two published tropical pedotransfer functions (PTFs) based on multiple linear regression (MLR). The overall estimation error ranges generated by the k-NN approach were statistically different but comparable to the two examined MLR PTFs. When the best combination of input variables (sand + silt + clay + bulk density + cation exchange capacity) was used, the overall error was remarkably low: 0.0360 to 0.0390 m(3) m(-3) in the dry and very wet ranges and 0.0490 to 0.0510 m(3) m(-3) in the intermediate range (i.e., -3 to -50 kPa) of the soil water retention curve. This k-NN variant can be considered as a competitive alternative to more classical, equation-based PTFs due to the accuracy of the water retention estimation and, as an added benefit, its flexibility to incorporate new data without the need to redevelop new equations. This is highly beneficial in developing countries where soil databases for agricultural planning are at present sparse, though slowly developing

Spatiotemporal Analysis of Groundwater Recharge Trends and Variability in Northern Taiwan

In this study, the base flow estimation method was used to assess long‐term changes of groundwater recharge in Northern Taiwan. The Mann‐Kendall test was used to examine the characteristics of the trends. This was followed by trend slope calculation and change‐point analysis. The annual groundwater recharge was found to exhibit a significant upward trend for the Fushan and Hengxi stations (Tamsui river basin). On the other hand, the Ximen Bridge station (Lanyang river basin) recorded a significant downward trend. Calculations showed that the rate of change for the Fengshan and Touqian river basins was small (less than 10%). However, that for the following stations was greater than 30%: Fushan, Hengxi, Ximen Bridge, and Niudou (also in the Lanyang river basin). The results of the change‐point analysis further indicated a significant change‐point for the annual recharge at Fushan, Hengxi, and Ximen Bridge stations in 1999, 1983, and 2001, respectively. The findings can be used for regional hydrological studies and as reference for water resource planning