Evaluation of a general model for multimodal unsaturated soil hydraulic properties

Many soils and other porous media exhibit dual- or multi-porosity type features. In a previous study (Seki et al., 2022) we presented multimodal water retention and closed-form hydraulic conductivity equations for such media. The objective of this study is to show that the proposed equations are practically useful. Specifically, dual-BC (Brooks and Corey)-CH (common head) (DBC), dual-VG (van Genuchten)-CH (DVC), and KO (Kosugi)$_1$BC$_2$-CH (KBC) models were evaluated for a broad range of soil types. The three models showed good agreement with measured water retention and hydraulic conductivity data over a wide range of pressure heads. Results were obtained by first optimizing water retention parameters and then optimizing the saturated hydraulic conductivity (K_s) and two parameters (p, q) or (p, r) in the general hydraulic conductivity equation. Although conventionally the tortuosity factor p is optimized and (q, r) fixed, sensitivity analyses showed that optimization of two parameters (p+r, qr) is required for the multimodal models. For 20 soils from the UNSODA database, the average $R^2$ for log (hydraulic conductivity) was highest (0.985) for the KBC model with r=1 and optimization of (Ks, p, q). This result was almost equivalent (0.973) to the DVC model with q=1 and optimization of (Ks, p, r); both were higher than $R^2$ for the widely used Peters model (0.956) when optimizing (Ks, p, a, $\omega$). The proposed equations are useful for practical applications while mathematically being simple and consistent.Comment: To be published in Journal of Hydrology and Hydromechanic

Potential of Mid-tropospheric Water Vapor Isotopes to Improve Large-Scale Circulation and Weather Predictability

Recent satellite techniques have uncovered detailed tropospheric water vapor isotope patterns on a daily basis, yet the significance of water isotopes on weather forecasting has remained largely unknown. Here, we perform a proof‐of‐concept observing system simulation experiment to show that mid‐tropospheric water isotopes observed by the Infrared Atmospheric Sounding Interferometer (IASI) can substantially improve weather forecasts through non‐local impacts on the convective heating structure and large‐scale circulation. Assimilating IASI isotopes can improve wind, humidity, and temperature fields by more than 10% at mid‐troposphere compared to only assimilating conventional non‐isotopic observations. These improvements are about two‐thirds of assimilating simultaneous IASI water vapor observations. The improvements can be attributed more to thermodynamic (phase change) effects than dynamic (transport) effects of water isotopes. Furthermore, isotopic observations produce additional 3%–4% improvements to the fields constrained by the conventional observations and simultaneous IASI water vapor observations, demonstrating the unique characteristics of water isotopes. Plain Language Summary Accurate weather forecasting has tremendous socio‐economic benefits by saving lives from natural hazards and affecting numerous sectors including water resources, energy, and agriculture. Although recent satellite techniques have enabled observing detailed water isotope patterns (e.g., HDΟ and Η$_{2}$$^{18}$O) in the atmosphere, it has not been incorporated in operational weather forecasting. Here, we show that water isotopes can substantially improve weather forecasts by improving the heating structure and large‐scale circulation. Satellite‐observed isotopes can improve wind, humidity, and temperature fields up to 3%–4% compared to utilizing conventional non‐isotopic observations and concurrent water vapor observations by the same satellite. We anticipate that our results will facilitate further modeling developments in isotopic processes and benefit the societal sectors by improving operational weather forecasting

Closed-form hydraulic conductivity equations for multimodal unsaturated soil hydraulic properties

Closed-form expressions of the hydraulic conductivity function for linearly superposed subretention (multimodal) functions were derived for arbitrary sets of the Brooks and Corey (BC), van Genuchten (VG), and Kosugi (KO) water retention models. The generalized Mualem hydraulic conductivity model was evaluated using the mathematical approach of Priesack and Durner. Three types of modification to the multimodel were also proposed. Firstly, the derived conductivity equations can be simplified when the submodel parameters, hbi for the BC model, alpha i-1 for the VG model, and hmi for the KO model have the same (common) value (denoted as CH). Secondly, as in the case of the modified single VG and KO models, a hypothetical air-entry head near saturation can be introduced for the multimodal VG and KO models to prevent unrealistic reductions in the hydraulic conductivity near saturation when the VG n parameter approaches its lower limit of n = 1. Furthermore, the multimodal hydraulic conductivity functions become a simple sum of conductivity subfunctions when the exponent r is unity (such as for Burdine's model), which leads to independent tortuosity effects for each submodel. The models are illustrated for two soils: a highly aggregated Kumamoto Andisol and a relatively unimodal dune sand. The dual-(BC, VG, KO) and the VG(1)BC(2) models equally represented the water retention data of the Andisol, with similar hydraulic conductivity curves. The dual-BC-CH, dual-VG-CH, and VG(1)BC(2)-CH models fitted the water retention data of the dune sand well, with the hydraulic conductivity curves of the dual-porosity model being similar to those of the Fayer and Simmons (FS) model

Mobility of Cr, Pb, Cd, Cu and Zn in a loamy sand soil : a comparative study

Interest in soil contamination has been growing in recent years due to the ongoing degradation of soil environments. Therefore, the development of remediation techniques and the study of contaminant sorption and migration are areas of intense research. In this study, the authors sought to evaluate the scenario of co-contamination of a loamy sand soil by multiple heavy metals. To that end, the sorption and transport of five metals—Cr, Pb, Cd, Cu and Zn—was evaluated using representative samples of a soil from the north of Portugal. The tests were conducted in batch and continuous systems using single- and multiple-metal acid solutions to evaluate the effect of metal competition. In accordance with the type of assay—batch or continuous—Langmuir or Convection Dispersion Two-Site Nonequilibrium models were adjusted to explain the sorption/transport data. FTIR analyses were performed on the final samples of the continuous systems. Generally, the results revealed good fitting of the tested models for the metals in competitive and noncompetitive scenarios, with the exception of Zn that was originally present in soil samples at higher concentrations. As expected, the influence of competition was observed in both batch and continuous systems, but with different tendencies. The FTIR spectra also revealed a strong influence of clay minerals and organic matter on the sorption of the metals.The PhD grants of Bruna Fonseca and Hugo Figueiredo and the research grant of Joana Rodrigues were financially supported by Fundacao para a Ciencia e Tecnologia, Ministerio da Ciencia e Tecnologia, Portugal and Fundo Social Europeu (FSE)

Removal of hexavalent chromium of contaminated soil by coupling electrokinetic remediation and permeable reactive biobarriers

PURPOSE: In this study, a novel and ecological alternative have been developed to treat soils contaminated with hexavalent chromium coupling two well-known systems: electrokinetic remediation and permeable reactive biobarriers. The electric field promotes the electromigration of the hexavalent chromium oxyanions towards the anode. The biobarriers were placed before the anode electrode, in order to promote the reduction and retention of the chromium migrating in its direction. Thus, this technology provided a global treatment to soil removal without subsequent treatments of the contaminated effluents. METHODS: The electrokinetic system was coupled with two different permeable reactive biobarriers composed by Arthrobacter viscosus bacteria, supported either in activated carbon or zeolite. An electric field of 10 V was applied and two different treatment times of 9 and 18 days were tested. RESULTS: Removal values of 60% and 79% were obtained when electrokinetic treatment was coupled with zeolite and activated carbon biobarriers, respectively, for a test period of 18 day. The reduction of hexavalent chromium to trivalent chromium was around 45% for both systems. CONCLUSIONS: In this work, two types of biobarriers were efficiently coupled to electrokinetic treatment to decontaminate soil with Cr(VI). Furthermore, the viability of the new coupling technology developed (electrokinetic + biobarriers) to treat low-permeability polluted soils was demonstrated.This work was supported by the Spanish Ministry of Science and Innovation (CTQ2008-03059/PPQ), Xunta de Galicia (08MDS034314PR). The authors are grateful to the Spanish Ministry of Science and Innovation for providing financial support for Marta Pazos under the Ramon y Cajal program and the Fundacao para a Ciencia e Tecnologia, Ministerio da Ciencia e Tecnologia, Portugal through the PhD grant of Bruna Fonseca (SFRH/BD/27780/2006)

Measured and Predicted Solute Transport in a Tile Drained Field

Most solute transport measurement techniques are tedious and require extensive soil excavation. A field experiment was conducted to evaluate whether surface transport properties determined by a nondestructive time domain reflectometry (TDR) technique could be used to accurately predict tile flux concentrations. A 14 by 14 m field plot selected above a 1.1-m deep tile drain was studied. Low electrical conductivity (EC) water was sprinkled on the plot surface, and after reaching a steady-state condition, a pulse of calcium chloride solution (16.3 cm) with an EC of 23 dS m−1 was applied through the same sprinklers. Time domain reflectometry equipment was used to record the change in EC of surface (∼ top 2 cm) soil at 45 locations. The EC of the tile drainage flow was measured continuously with an EC probe. The surface convective lognormal transfer (CLT) function parameters, log mean irrigation depth, μI, and its standard deviation, σI, were found to be 3.44 and 0.94 [ln(cm)], respectively, for a reference depth of 110 cm. These surface parameters were used in a one-dimensional (1-D) CLT model and in a two-dimensional (2-D) model (CLT vertical function combined with exponential horizontal transfer function) to predict the tile flux concentrations. The 1-D CLT model predicted an earlier arrival time of chemicals to the tile drain than observed values. The root mean square error, RMSE, of the 1-D CLT predictions was 0.123, and the coefficient of efficiency, E, was −0.47. The 2-D model predictions of tile flux concentrations were similar to the observed values. The root mean squared errors (RMSE) and E were 0.023 and 0.94, respectively. The findings suggest that in this field soil, the surface solute transport properties determined by TDR could be combined with a 2-D transport model to make reasonable predictions of tile flux concentrations