Aly I. El-Kadi

A lumped parameter model for nitrogen transformation in the unsaturated zone

Prediction of NO3 -N leaching in the unsaturated zone is difficult owing to the complex nature of processes dominating the fate of N in soils. Although many sophisticated numerical models have been developed for simulating nitrogen behavior in the soil-plant system, applying them in practice is often difficult, owing to the limited specific-site information. A lumped parameter model (LPM) for the assessment of nitrate leaching using limited data is presented. An LPM considers nitrogen transformation processes, climate data, and irrigation and fertilizer applications. The model was tested against detailed numerical models and against field data and provided realistic estimates of leachate at the bottom of the root zone. The model can be used for decision making regarding of nitrogen-N application and can provide an alternative way to estimate leachate loading to the groundwater.

Analysis of an Unconfined Aquifer Subject to Asynchronous Dual-Tide Propagation

Most published solutions for aquifer responses to ocean tides focus on the one-sided attenuation of the signal as it propagates inland. However, island aquifers experience periodic forcing from the entire coast, which can lead to integrated effects of different tidal signals, especially on narrow high-permeability islands. In general, studies disregard a potential time lag as the tidal wave sweeps around the island. We present a one-dimensional analytical solution to the ground water flow equation subject to asynchronous and asymmetric oscillating head conditions on opposite boundaries and test it on data from an unconfined volcanic aquifer in Maui. The solution considers sediment-damping effects at the coastline. The response of Maui Aquifers indicate that water table elevations near the center of the aquifer are influenced by a combination of tides from opposite coasts. A better match between the observed ground water head and the theoretical response can be obtained with the proposed dual-tide solution than with single-sided solutions. Hydraulic diffusivity was estimated to be 2.3 x 10(7) m(2)/d. This translates into a hydraulic conductivity of 500 m/d, assuming a specific yield of 0.04 and an aquifer thickness of 1.8 km. A numerical experiment confirmed the hydraulic diffusivity value and showed that the y-intercepts of the modal attenuation and phase differences estimated by regression can approximate damping factors caused by low-permeability units at the boundary.

The Courant and Peclet Number criteria for the numerical solution of the Richards Equation

Recent studies indicate that inaccuracies exist in the numerical solutions of the Richards equation. A need exists to examine such inaccuracies and to correlate them, formally, to the soil hydraulic properties and the size of the spatial and temporal mesh. In the absence of criteria for mesh design, mass conservation and the rate of convergence are usually used as bases for acceptance of the solution. It is well documented, however, that neither can guarantee an accurate solution. This paper proposes the Courant and Peclet numbers as criteria for estimating spatial and temporal increments. The study defines upper limits of these numbers for acceptable accuracy of the numerical solutions. For cases that are characterized by decreasing pressure head gradients, although the initial size of the time step may be extremely small, efficient solutions are possible through a careful choice of an expansion factor for such a step.

NETPATH-WIN: An Interactive User Version of the Mass-Balance Model, NETPATH

NETPATH-WIN is an interactive user version of NETPATH, an inverse geochemical modeling code used to find mass-balance reaction models that are consistent with the observed chemical and isotopic composition of waters from aquatic systems. NETPATH-WIN was constructed to migrate NETPATH applications into the Microsoft WINDOWS® environment. The new version facilitates model utilization by eliminating difficulties in data preparation and results analysis of the DOS version of NETPATH, while preserving all of the capabilities of the original version. Through example applications, the note describes some of the features of NETPATH-WIN as applied to adjustment of radiocarbon data for geochemical reactions in groundwater systems.

Logistic regression modeling to assess groundwater vulnerability to contamination in Hawaii, USA

Capture zone analysis combined with a subjective susceptibility index is currently used in Hawaii to assess vulnerability to contamination of drinking water sources derived from groundwater. In this study, we developed an alternative objective approach that combines well capture zones with multiple-variable logistic regression (LR) modeling and applied it to the highly-utilized Pearl Harbor and Honolulu aquifers on the island of Oahu, Hawaii. Input for the LR models utilized explanatory variables based on hydrogeology, land use, and well geometry/location. A suite of 11 target contaminants detected in the region, including elevated nitrate (>1 mg/L), four chlorinated solvents, four agricultural fumigants, and two pesticides, was used to develop the models. We then tested the ability of the new approach to accurately separate groups of wells with low and high vulnerability, and the suitability of nitrate as an indicator of other types of contamination. Our results produced contaminant-specific LR models that accurately identified groups of wells with the lowest/highest reported detections and the lowest/highest nitrate concentrations. Current and former agricultural land uses were identified as significant explanatory variables for eight of the 11 target contaminants, while elevated nitrate was a significant variable for five contaminants. The utility of the combined approach is contingent on the availability of hydrologic and chemical monitoring data for calibrating groundwater and LR models. Application of the approach using a reference site with sufficient data could help identify key variables in areas with similar hydrogeology and land use but limited data. In addition, elevated nitrate may also be a suitable indicator of groundwater contamination in areas with limited data. The objective LR modeling approach developed in this study is flexible enough to address a wide range of contaminants and represents a suitable addition to the current subjective approach.

Bioaccessible arsenic in soils of former sugar cane plantations, Island of Hawaii.

Arsenical herbicides were used extensively for emergent weed control in Hawaiian sugar cane cultivation from 1913 to about 1950. As a result, surface soil arsenic concentrations average 280 mg kg(-1) across more than 60 km(2) of former sugar plantation land in the eastern portion of the Island of Hawaii. This study was conducted to elucidate the relationship between soil properties and arsenic bioaccessibility in the iron-rich volcanic soils. Soils are predominantly Andisols, formed by weathering of basaltic lava and tephra, with pedogenic solid phases consisting of short-range order iron oxyhydroxides, allophane-like aluminosilicates, and metal-humus compounds. These reactive solid phases strongly adsorb oxyanions, such as phosphate and arsenite/arsenate. High arsenic sorption capacity limits desorption and vertical migration within the soil column and prevents contamination of the underlying groundwater aquifer, despite high arsenic loading and precipitation rates. In vitro arsenic bioaccessibility, as measured by the SBRC gastric-phase test, ranges from 2% to 35% and averages 9% of total arsenic. Bioaccessible arsenic is higher in less weathered soils (Udifolists, Typic and Lithic Hydrudands) and lower in more weathered ash-dominant soils (Acrudoxic Hydrudands). Soil weathering indicators, such as reactive iron content, are strong predictors of arsenic bioaccessibility. Based on evidence from soil mineralogy, geochemistry and arsenic speciation, as well as limited soil arsenic bioavailability/bioaccessibility comparisons, risks to human health from direct contact (soil ingestion) are significantly reduced by low arsenic bioaccessibility. Nonetheless, some soils within former sugar cane cultivation areas contain bioaccessible arsenic concentrations exceeding Hawaii Department of Health risk-based action levels, and will require mitigating actions. Even higher levels of soil arsenic contamination have been identified at former pesticide storage and mixing areas, but are generally of localized extent.

On estimating the hydraulic properties of soil, Part 1. Comparison between forms to estimate the soil-water characteristic function

Abstract The hydraulic properties of soil include the soil-water characteristic function [ h ( θ ), in which θ is water content and h is pressure head (suction)], and the hydraulic conductivity function [ K ( θ ) or K ( h )]. These functions are essential to the solution of unsaturated groundwater flow problems. A number of empirical and semiempirical forms have been proposed in the literature to estimate these functions. The present paper employs a nonlinear least-square analysis for comparison between some of the available forms, using a large number of experimental measurements of h ( θ ) for different classes of soil. Suitability of the forms for predicting the hydraulic conductivity function is examined. In the absence of accurate measurements, the paper facilitates modeling by providing estimates for the parameters of the soil-water characteristic function.

Modeling hydrocarbon biodegradation in tidal aquifers with water-saturation and heat inhibition effects

A model is developed for hydrocarbon biodegradation, which includes saturated and unsaturated flow, multi-species transport, heat transport, and bacterial growth processes. Numerical accuracy of the model was tested against analytical solutions. The model was also verified against laboratory results for a saturated-flow problem and reasonable match was obtained. Expressions are proposed for inhibition due to water content and temperature fluctuations. Bioactivities under cyclic water content variation were studied under no-flow conditions. A quantitative approach was used to reconcile some of the apparent contradictory conclusions regarding the efficiency of biodegradation of soils under wetting and drying conditions. The efficiency depends on the nature of the oxygenation process. For cases involving the presence of dissolved oxygen and the absence of O2 vapor, subjecting the soil to constant water content close to its optimal value for degradation is most efficient. However, wetting and drying can enhance degradation if O2 is only provided through aeration or direct contact between air and the medium. Also presented are the results of a typical field application of the model and a discussion of the effects of tides, saturation inhibition, and heat inhibition. Other inhibition factors, such as pH or salinity, can be easily incorporated in the formulation. The quantitative approach developed here can be used in assessing bioremediation not only in tidal aquifers but also in areas where water-table or temperature effects are of significance. The approach can be useful in the design of remediation strategies under water-flow or no-flow conditions involving water content and temperature fluctuations.

Iron amendments to reduce bioaccessible arsenic.

Former sugarcane lands on the Island of Hawaii have elevated levels of soil arsenic (As) from historical use of arsenical pesticides. The bioaccessible fraction of total As (AsTOT), a measure of the potential for human As uptake by incidental ingestion of soil, is used in the assessment of human health risk and the determination of the need for remedial action. Ferric chloride plus lime and ferrous sulfate plus lime were applied to As-contaminated soils in a field plot setting to determine the potential for reducing in vitro bioaccessible As (AsIVBA) by increasing As sequestration by the formation of additional iron (Fe) oxyhydroxides. The two Fe sources performed similarly in reducing AsIVBA over a 2-year observation period, with 30-41% reduction in AsIVBA for 0.25wt% Fe dosing (dry soil basis) and 59-63% reduction for 0.5wt% Fe dosing. Addition of phosphate (PO4) to treated and untreated soils caused a significant increase in AsIVBA. Iron-treated and control soils showed more than twice the AsIVBA after the addition of 1500mgPkg(-1). The cost of in situ treatment of As-contaminated soil with ferrous sulfate plus lime to lower AsIVBA was estimated to be an order of magnitude less than excavation and landfill disposal on the Island of Hawaii, making the technology a viable alternative when remedial action objectives were based on AsIVBA levels.

On estimating the hydraulic properties of soil, Part 2. A new empirical equation for estimating hydraulic conductivity for sands

A new empirical equation to estimate hydraulic conductivity is proposed, based on a large set of measured data for hydraulic properties of soil. The equation is simpler and more accurate than the series-parallel model. Under conditions of insufficient data, the new equation provides a good estimation of hydraulic conductivity for sands. For the same class of soils, another empirical equation is proposed to estimate the power N in the Averjanov-Irmay function.