SPATIAL VARIABILITY OF CONDUCTWITY AS APPLIED TO DISTRIBUTED PARAMETER INFILTRATION MODELS1

ABSTRACT: Two infiltration models, called VVSIM (variable variance stochastic infiltration model) and EVVSIM (enhanced variable variance stochastic infiltration model), are developed in this study. A distributed parameter infiltration model can estimate the amount of infiltration over a field area by computing the infiltration over zones of the field area. Hydraulic conductivity is the most important parameter determining infiltration in simulations by infiltration models. The performance of an infiltration model depends on how well the model accommodates the complicated spatial distribution of hydraulic conductivity. The two proposed models include the effects of spatial correlation of the conductivity distribution. Virtual conductivity fields are generated using the turning bands method. Monte Carlo simulations show that the proposed models give infiltration estimates more accurate than those obtained by the other models employed.     

Developing risk models of Cryptosporidium transport in soils from vegetated, tilted soilbox experiments

Transport of Cryptosporidium parvum through macroporous soils is poorly understood yet critical for assessing the risk of groundwater contamination. We developed a conceptual model of the physics of flow and transport in packed, tilted, and vegetated soilboxes during and immediately after a simulated rainfall event and applied it to 54 experiments implemented with different soils, slopes, and rainfall rates. Using a parsimonious inverse modeling procedure, we show that a significant amount of subsurface outflow from the soilboxes is due to macropore flow. The effective hydraulic properties of the macropore space were obtained by calibration of a simple two-domain flow and transport model that accounts for coupled flow in the matrix and in the macropores of the soils. Using linear mixed-effects analysis, macropore hydraulic properties and oocyst attenuation were shown to be associated with soil bulk density and rainfall rate. Macropore flow was shown to be responsible for bromide and C. parvum transport through the soil into the underlying pore space observed during the 4-h experiments. We confirmed this finding by conducting a pair of saturated soil column studies under homogeneously repacked conditions with no macropores in which no C. parvum transport was observed in the effluent. The linear mixed-effects and logistic regression models developed from the soilbox experiments provide a basis for estimating macropore hydraulic properties and the risk of C. parvum transport through shallow soils from bulk density, precipitation, and total shallow subsurface flow rate. The risk assessment is consistent with the reported occurrence of oocysts in springs or groundwater from fractured or karstic rocks protected only by shallow overlying soils.     

A STRIP MODEL APPROACH TO PARAMETERIZE A COUPLED GREEN‐AMPT KINEMATIC WAVE MODEL1

ABSTRACT: Infiltration processes at the plot scale are often described and modeled using a single effective hydraulic conductivity (Kg) value. This can lead to errors in runoff and erosion prediction. An integrated field measurement and modeling study was conducted to evaluate: (1) the relationship among rainfall intensity, spatially variable soil and vegetation characteristics, and infiltration processes; and (2) how this relationship could be modeled using Green and Ampt and a spatially distributed hydrologic model. Experiments were conducted using a newly developed variable intensity rainfall simulator on 2 m by 6 m plots in a rangeland watershed in southeastern Arizona. Rainfall application rates varied between 50 and 200 mm/hr. Results of the rainfall simulator experiments showed that the observed hydrologic response changed with changes in rainfall intensity and that the response varied with antecedent moisture condition. A distributed process based hydrologic simulation model was used to model the plots at different levels of hydrologic complexity. The measurement and simulation model results show that the rainfall runoff relationship cannot be accurately described or modeled using a single Kg value at the plot scale. Multi‐plane model configurations with infiltration parameters based on soil and plot characteristics resulted in a significant improvement over single‐plane configurations.     

Mathematical modeling of phosphorus losses from land application of hog and cattle manure

Mathematical models may provide a means to estimate phosphorus (P) losses from land application of manure. Phosphorus losses typically occur during brief episodes of runoff and erosion. Models must be able to simulate P losses during these episodes by representing the basic chemical, physical, and biological processes by which these losses occur. The mathematical model ecosys combines dynamic distributed flow of solutes and nonsolutes through runoff and erosion with convective-dispersive transport of solutes, and both biologically and thermodynamically driven transformations between solutes and nonsolutes. This model was tested against P lost in runoff, erosion, and leachate measured during 90 min of controlled rainfall at 65 mm h(-1) on soils from six sites at which different rates of manure had been applied over the previous 3 to 6 yr. Transport and transformation kinetics in the model enabled it to simulate changes of dissolved inorganic phosphorus (DIP) in runoff from >1.0 to <0.05 mg L(-1) and changes of total phosphorus (TP) in sediment from 15 to 3 mg L(-1) measured during controlled rainfall on soils with diverse P contents. Results from 60-yr model runs using these kinetics with different application rates of cattle manure indicated that (i) a positive interaction exists between annual rainfall and application rate on P losses and (ii) rates greater than 30 Mg ha(-1) yr(-1) would cause TP concentrations in water leaving the site to rise above acceptable limits. The interaction between rainfall and rate suggests that P losses from manure application at any site should be assessed under the upper range of likely rainfall intensities.     

EVALUATION OF AN INDUCED INFILTRATION MODEL AS APPLIED TO GLACIAL AQUIFER SYSTEMS1

ABSTRACT: Numerical models were used to examine the limitations of the assumptions used in an analytical induced infiltration model. The assumptions tested included negligible streambed effects, negligible areal recharge, two-dimensional ground water flow, fully penetrating rivers and wells, and constant surface water stage. For situations that deviate from the underlying assumptions, the analytical model becomes a less reliable predictor of induced infiltration. The numerical experiments show that streambed effects cannot be neglected if the streambed conductivity is more than one order of magnitude lower than the aquifer hydraulic conductivity. Areal recharge cannot be neglected if the ground water basin receives more than 5 in/yr of areal recharge. Three-dimensional flow effects due to well partial penetration cannot be neglected if the ratio of horizontal hydraulic conductivity to vertical hydraulic conductivity (K_h/K_u) is greater than 10. Surface water elevation fluctuations can significantly influence the induced infiltration rate, depending on the degree of fluctuations and the ground water hydraulic gradient.     

IMPROVEMENT OF LEACHATE PREDICTION THROUGH MUNICIPAL SOLID WASTE LAYERS1

ABSTRACT: Existing leachate prediction models over- or underestimate leachate generation by up to three orders of magnitude. Practical experiments show that channeled flow in waste leads to rapid discharge of large leachate volumes and heterogeneous moisture distributions. In order to more accurately predict leachate generation, leachate models must be improved. To predict moisture movement through waste, the one-domain water balance model HELP, and two-domain PREFLO, are tested. Experimental waste and leachate flow values are compared with model predictions. When calibrated with experimental parameters, the HELP model yields reasonable predictions of cumulative leachate flow and PREFLO provides estimates of breakthrough time. In the short term, field capacity has to be reduced to 0.12 and effective storage and hydraulic conductivity of the waste must be increased to 0.2 and 2.2 cm/s respectively. In the long term, a new modeling approach must be developed to adequately describe the moisture movement mechanisms.     

RASTER-BASED HYDROLOGIC MODELING OF SPATIALLY-VARIED SURFACE RUNOFF1

ABSTRACT: The proliferation of watershed databases in raster Geographic Information System (GIS) format and the availability of radar-estimated rainfall data foster rapid developments in raster-based surface runoff simulations. The two-dimensional physically-based rainfall-runoff model CASC2D simulates spatially-varied surface runoff while fully utilizing raster GIS and radar-rainfall data. The model uses the Green and Ampt infiltration method, and the diffusive wave formulation for overland and channel flow routing enables overbank flow storage and routing. CASC2D offers unique color capabilities to display the spatio-temporal variability of rainfall, cumulative infiltrated depth, and surface water depth as thunderstorms unfold. The model has been calibrated and independently verified to provide accurate simulations of catchment response to moving rainstorms on watersheds with spatially-varied infiltration. The model can accurately simulate surface runoff from flashfloods caused by intense thunderstorms moving across partial areas of a watershed.     

APPLICATION OF THE GREEN-AMPT INFILTRATION EQUATION TO WATERSHED MODELING1

ABSTRACT: A computer model was developed, based on the Green-Ampt infiltration equation, to computed rainfall excess for a single precipitation event. The model requires an estimate of parameters related to hydraulic conductivity, wetting front section, and fillable porosity of the soil layers. Values of parameters were estimated from soil textural averages or regression equations based on percent sand, percent clay, and porosity. Average values of effective porosity and wetting front suction were largely acceptable due to the relatively low variability and low model sensitivity to the parameters. Hydraulic conductivity was the most erratic constituent of the loss rate computation due to the high variability and the high sensitivity of the computed infiltration to the parameter. The performance of the Green-Ampt infiltration model was tested through a comparison with the SCS curve number procedure. Seven watersheds and 23 storms with precipitation of one inch or greater were used in the comparison. For storms with less than one inch of rainfall excess, the SCS curve number procedure generally gave the best results; however, for six of the seven storms with precipitation excess greater than one inch, the Green-Ampt procedure delivered better results. In this comparison, both procedures used the same initial abstractions. The separation of rainfall losses into infiltration, interception, and surface retention is, in theory, an accurate method of estimating precipitation excess. In the second phase of the study using nine watersheds and 39 storms, interception and surface retention losses were computed by the Horton equations. Green-Ampt and interception parameters were estimated from value sin the literature, while the surface retention parameter was calibrated so that the computed runoff volumes matched observed volumes. A relationship was found between the surface retention storage capacity and the 15-day antecedent precipitation index, month of year, and precipitation amount.     

A DYNAMIC SYSTEM MODEL FOR SIMULATING SEDIMENT DISCHARGE1

ABSTRACT: A dynamic sediment discharge model was developed and proposed for the simulation of watershed systems. It war developed from an expansion of splash and flow erosion relationships under steady state conditions. It was described as a general erosion model that can be reduced to forms comparable to many conceptual soil erosion and sediment yield models. The model incorporates eight parameters such as rainfall intensities, runoff rates, and previous sediment discharges. The model was tested with two small watersheds with simulation results which were very satisfactory compared to the data.     

Design and risk assessment tool for vegetative treatment areas receiving agricultural wastewater: Preliminary results

Vegetative treatment areas (VTAs) are commonly being used as an alternative method of agricultural process wastewater treatment. However, it is also apparent that to completely prevent discharge of pollutants to the surrounding environment, settling of particulates and bound constituents from overland flow through VTAs is not sufficient. For effective remediation of dissolved agricultural pollutants, VTAs must infiltrate incoming wastewater. A simple water balance model for predicting VTA soil saturation and surface discharge in landscapes characterized by sloping terrain and a shallow restrictive layer is presented and discussed. The model accounts for the cumulative effect of successive rainfall events and wastewater input on soil moisture status and depth to water table. Nash–Sutcliffe efficiencies ranged from 0.65 to 0.81 for modeled and observed water table elevations after calibration of saturated hydraulic conductivity. Precipitation data from relatively low, average, and high annual rainfall years were used with soil, site, and contributing area data from an example VTA for simulations and comparisons. Model sensitivity to VTA width and contributing area (i.e. barnyard, feedlot, silage bunker, etc.) curve number was also investigated. Results of this analysis indicate that VTAs should be located on steeper slopes with deeper, more-permeable soils, which effectively lowers the shallow water table. In sloping landscapes (>2%), this model provides practitioners an easy-to-use VTA design and/or risk assessment tool that is more hydrological process-based than current methods.     

Model Development and Calibration of a Saltwater Intrusion Model in Southern California1

Abstract: A nine‐layered confined‐unconfined flow and transport model is developed for the Alamitos saltwater intrusion barrier in Southern California. The conceptual model is based on the geological structure of the coastal aquifer system. The key parameters in the flow and transport models are calibrated using a two‐phase procedure which matches the types of data available for calibration. Because of the abundance of point measurements of hydraulic conductivity, the heterogeneous and random hydraulic conductivity field for each of the five aquifers is estimated by the geostatiscal method of natural‐neighbor‐kriging in Phase 1. In Phase 2, the longitudinal and transverse dispersivities in the transport model are estimated by a traditional inverse procedure that minimizes the least‐squares error for concentration (LSE‐CON). The minimum LSE‐CON is achieved near 15.2 and 1.52 m for the longitudinal and transverse dispersivities, respectively. Additional simulations with increasing transport parameter complexity did not yield significant improvements in LSE‐CON. Also, tracking least‐squares error for head while parametrically varying the transport parameters revealed there is a negligible interaction between predicted head and transport parameters.     

MODELING THE FOREST HYDROLOGY OF WETLAND-UPLAND ECOSYSTEMS IN FLORIDA1

ABSTRACT: Few hydrological models are applicable to pine flat-woods which are a mosaic of pine plantations and cypress swamps. Unique features of this system include ephemeral sheet flow, shallow dynamic ground water table, high rainfall and evapotranspiration, and high infiltration rates. A FLATWOODS model has been developed specifically for the cypress wetland-pine upland landscape by integrating a 2-D ground water model, a Variable-Source-Area (VAS)-based surface flow model, an evapotranspiration (ET) model, and an unsaturated water flow model. The FLATWOODS model utilizes a distributed approach by dividing the entire simulation domain into regular cells. It has the capability to continuously simulate the daily values of ground water table depth, ET, and soil moisture content distributions in a watershed. The model has been calibrated and validated with a 15-year runoff and a four-year ground water table data set from two different pine flat woods research watersheds in northern Florida. This model may be used for predicting hydrologic impacts of different forest management practices in the coastal regions.     

The effect of flow path and mixing layer on phosphorus release: physical mechanisms and temperature effects

Soil phosphorus (P) concentrations typically are greater in surface soils compared with subsurface soils. Surface soils have a greater chance to interact with runoff leading to P transport to streams. The thin surface layer where P concentrates is referred to as the mixing layer denoting where water and chemicals mix during transport. The objective of this study was to evaluate the effect of hydrologic flow paths on soluble reactive phosphorus (SRP) loss at two temperatures. Laboratory flumes were built to simulate infiltration, return flow, saturation excess, and interflow, and subsequent interaction with the mixing layer. The sandy loam soil in the flumes was kept at saturation throughout all experiments, so that biochemical effects were normalized. Flow through the flumes was maintained at 3.6 mm/h for 24 to 99 h (at 6 and 25 degrees C) with water entering and exiting the flumes at different ports (to simulate different flow paths) or as low intensity rainfall. Experiments were performed with and without an artificially created P-enriched surface layer (5 mm thick, total P increased from 1010 mg/kg in the original soil to 2310 mg/kg by addition of dissolved phosphate). Results indicated that (i) SRP release was greater in soil with a mixing layer than in soil without a mixing layer; (ii) SRP release was greater during experiments at 25 degrees C than at 6 degrees C; (iii) at 25 degrees C, SRP release was greatest when water traversed the mixing layer in the upward direction (i.e., in return flow), and by flow parallel to the mixing layer (i.e., surface runoff); and (iv) at 6 degrees C, SRP release in subsurface flow following rainfall was slightly greater than in return flow and infiltration. Our results confirmed the presence of a variable, temperature-dependent desorption process when runoff water interacted with the mixing layer. Our findings have important implications for how different water flow paths in and over the soil interact with P in the soil, and what the ultimate concentration will be in runoff and interflow.     

MODELING GROUND WATER FLOW USING FLUX BOUNDARY CONDITIONS1

ABSTRACT: Determination of the boundary conditions for modeling ground water flow is a critical point especially in regional models. Normally the regional models require model areas that are greater than the given area of interest. This work focuses on the prediction of hydraulic heads in regional models using flux boundary conditions. The model uses flux boundary conditions that were estimated using a radial flow analog and Darcy's law. The regional model that is presented uses no parameter identification (inverse estimation) procedures. In the present work, the Houston area was used. The simulation of the hydrological conditions of the Chicot and Evangeline Aquifers that underlie the Houston area were made using the available information about the geological profile in the Houston region and the current information about the existing production wells. The regional model works as a forward problem. The system parameters such as hydraulic conductivity, specific storage, and hydrological stresses were specified, and the model predicts the hydraulic head. Actual data from piezometers operated by the U.S. Geological Survey (USGS) in many places throughout Houston were used as initial conditions. Some piezometric head data were generated using the regional variable theory called kriging to supply head estimates in areas where data were unavailable. The Modular Three Dimensional Finite Difference Groundwater Flow Model developed by the USGS was used to predict the hydraulic heads. The predicted ground water heads are compared to the actual data. The results show that the model performs well for locations where data were available.     

INFILTRATION EVALUATION OF FOUR MECHANICAL RAINFALL SIMULATION TECHNIQUES IN SIERRA NEVADA WATERSHEDS1

ABSTRACT: Little quantitative site-specific infiltration, runoff and sediment transport data for Tahoe Basin soils under varying storm events or stage of development are available. Modular (Ml), F-type (M2), Impact nozzle (M3), and Impact-Fan nozzle (M4) rainfall simulators were evaluated as to their practicality and ability to characterize infiltration for the Cagwin Soil Series within the Tahoe Basin. Three slope (0–15,15–30, >30%) and four plot conditions (natural with duff [P1], natural without duff [P2], disturbed without duff [P3], and disturbed with duff [P4]) were studied. The measured data were incorporated into a modified Philip's infiltration model and multiple non-linear regression analyses were used to examine relationships between method, slope, plot condition, and infiltration characteristics.^t Simulation methods Ml and M4 produced statistically similar (P=0.01) infiltration data, as did M2 and M3 which produced lower infiltration rates. All were found suitable for use in Sierra Nevada watersheds. Ml was considered most practical. Slope had negligible effect on infiltration. The plot condition was found to significantly influence infiltration, and the effect of each plot condition was significantly different. Final infiltration rates ranged from 4.7 to 6.2 cm/hr. Thus, the Cagwin soil demonstrated moderate to high infiltration rates even when exposed to extreme storm conditions (8–10 cm/hr).     

Tillage and field scale controls on greenhouse gas emissions

There is a lack of understanding of how associations among soil properties and management-induced changes control the variability of greenhouse gas (GHG) emissions from soil. We performed a laboratory investigation to quantify relationships between GHG emissions and soil indicators in an irrigated agricultural field under standard tillage (ST) and a field recently converted (2 yr) to no-tillage (NT). Soil cores (15-cm depth) were incubated at 25 degrees C at field moisture content and 75% water holding capacity. Principal component analysis (PCA) identified that most of the variation of the measured soil properties was related to differences in soil C and N and soil water conditions under ST, but soil texture and bulk density under NT. This trend became more apparent after irrigation. However, principal component regression (PCR) suggested that soil physical properties or total C and N were less important in controlling GHG emissions across tillage systems. The CO2 flux was more strongly determined by microbial biomass under ST and inorganic N content under NT than soil physical properties. Similarly, N2O and CH4 fluxes were predominantly controlled by NO3- content and labile C and N availability in both ST and NT soils at field moisture content, and NH4+ content after irrigation. Our study indicates that the field-scale variability of GHG emissions is controlled primarily by biochemical parameters rather than physical parameters. Differences in the availability and type of C and N sources for microbial activity as affected by tillage and irrigation develop different levels and combinations of field-scale controls on GHG emissions.     

SOLUTE TRANSPORT IN THREE‐DIMENSIONAL SATURATED POROUS MEDIA HAVING SELF‐SIMILAR HYDRAULIC CONDUCTIVITY DISTRIBUTION1

ABSTRACT: Customarily, it has been assumed that hydraulic conductivity is a stationary, homogeneous stochastic process with a finite variance for stochastic analysis of solute transport in the subsurface. That the distribution of hydraulic conductivity may have a fractal behavior with long range correlations was suggested from field data analyses. This motivates us to further investigate how the fractal behavior of permeability distribution impacts solute transport in porous media. This study provides longitudinal and transverse macrodispersivity coefficients and the variance of the solute concentration. Longitudinal and transverse macrodispersivity coefficients are found to depend strongly on the fractal dimension (D) of logarithmic hydraulic conductivity (logK). The longitudinal and transverse macrodispersivity coefficients are the highest when D = 1, and the values decrease monotonically to zero at D = 2. Both coefficients correspond to the characteristic length scale of the logK distribution, thus are scale dependent parameters. The ratio of the transverse to the longitudinal macrodispersivity coefficient is on the order of 10^‐1 to 10^‐4. Concentration variance also decreases with the fractal dimension of logK. There is no spatial spreading of solute for D = 2, and the concentration variance reaches zero for this case.     

ASSESSING CHANGES IN WATERSHED FLOW REGIMES WITH SPATIALLY EXPLICIT HYDRAULIC MODELS1

ABSTRACT: Urbanization, farming, and other watershed activities can significantly alter storm hydrographs and sediment erosion rates within a watershed. These changes routinely cause severe economic and ecological problems manifested in the form of increased flooding and significant changes in channel morphology. As the activities within a watershed influence the hydrologic, hydraulic, and ecological conditions within a river, interdisciplinary approaches to predict and assess the impacts that different land uses have on streams need to be developed. An important component of this process is ascertaining how hydrologic changes induced by specific watershed activities will affect hydraulic conditions and the accompanying flood levels, sediment transport rates, and habitat conditions within a stream. A conceptual model for using spatially explicit (two‐dimensional) hydraulic models to help evaluate the impacts that changes in flow regime might have on a river is presented. This framework proposes that reproducing and quantifying flow complexity allows one to compare the hydraulic conditions within urban, urbanizing, and non‐urban streams in a more biologically and economically meaningful way. The justification, advantage, and need for such a method is argued through the results of one‐ and two‐dimensional hydraulic model studies. The implementation of this methodology in watershed urbanization studies is described.     

DEVELOPMENT OF A BACTERIAL TRANSPORT MODEL FOR COARSE SOILS1

ABSTRACT: A bacterial transport model, developed to analyze bacterial translocation in coarse-grained soils, is presented. The complex governing equation is presented first, followed by analyses of each of the major processes influencing bacterial transport. These analyses suggest simplification of the governing equation is feasible when input data on specific processes are limited or unavailable. Model parameters, including bacterial die-off, bacterial distribution, input bacterial concentration, and saturated hydraulic conductivity, were randomly generated using a procedure known to produce either a normal or log-normal distribution of random numbers. Monte Carlo simulations were completed, and the resulting output was used to generate cumulative frequency distributions showing the probability of bacterial transport beyond various soil depths. Results from these simulations indicate that bacteria have a high probability of traveling through coarse-grained soils when low clay content and soil water temperatures limit bacterial retention and die-off.