COMPARISON OF METHODS TO ESTIMATE DEEP PERCOLATION RATES1

ABSTRACT: Deep percolation rates are normally estimated from a water balance. Results are presented of a study undertaken to evaluate three alternative methods of estimating percolation below the root zone when knowledge about the history of applied water and evapotranspiration are not available. The alternative methods are: 1) use of Darcy's equation to calculate deep percolation rate; 2) measurement of the soil temperature prof and calculation of the deep percolation rate from the shape of the temperature depth curve; and 3) measurement of the tritium concentration in the soil water and its relationship to the history of the tritium concentration in rainfall. At the principal study site, the Darcy velocity of flow ranged from 9 cm per year determined by the temperature method, to 40 cm per year determined by the tritium method. Darcy's equation to calculate seepage rates resulted in an estimation of deep seepage of 18 cm per year. An average deep percolation rate at the principal study site of 22 cm per year was determined using the average of all three methods. Results for other sites based on the temperature method indicated a lower seepage rate.     

Characterizing scale- and location-dependent correlation of water retention parameters with soil physical properties using wavelet techniques

Understanding the correlation between soil hydraulic parameters and soil physical properties is a prerequisite for the prediction of soil hydraulic properties from soil physical properties. The objective of this study was to examine the scale- and location-dependent correlation between two water retention parameters (alpha and n) in the van Genuchten (1980) function and soil physical properties (sand content, bulk density [Bd], and organic carbon content) using wavelet techniques. Soil samples were collected from a transect from Fuxin, China. Soil water retention curves were measured, and the van Genuchten parameters were obtained through curve fitting. Wavelet coherency analysis was used to elucidate the location- and scale-dependent relationships between these parameters and soil physical properties. Results showed that the wavelet coherence between alpha and sand content was significantly different from red noise at small scales (8-20 m) and from a distance of 30 to 470 m. Their wavelet phase spectrum was predominantly out of phase, indicating negative correlation between these two variables. The strong negative correlation between alpha and Bd existed mainly at medium scales (30-80 m). However, parameter n had a strong positive correlation only with Bd at scales between 20 and 80 m. Neither of the two retention parameters had significant wavelet coherency with organic carbon content. These results suggested that location-dependent scale analyses are necessary to improve the performance for soil water retention characteristic predictions.     

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.     

THE FLEXIBLE TOLERANCE METHOD FOR ESTIMATING HYDROLOGIC PARAMETERS IN THE ROOT ZONE1

In this study, a constrained minimization method, the flexible tolerance method, was used to solve the optimization problems for determining hydrologic parameters in the root zone: water uptake rate, spatial root distribution, infiltration rate, and evaporation. Synthetic soil moisture data were first generated using the Richards' equation and its associated initial and boundary conditions, and these data were then used for the inverse analyses. The results of inverse simulation indicate the following. If the soil moisture data contain no noise, the rate of estimated water uptake and spatial root distribution parameters are equal to the true values without using constraints. If there is noise in the observed data, constraints must be used to improve the quality of the estimate results. In the estimation of rainfall infiltration and surface evaporation, interpolation methods should be used to reduce the number of unknowns. A fewer number of variables can improve the quality of inversely estimated parameters. Simultaneous estimation of spatial root distribution and water uptake rate or estimation of evaporation and water uptake rate is possible. The method was used to estimate the water uptake rate, spatial root distribution, infiltration rate, and evaporation using long‐term soil moisture data collected from Nebraska's Sand Hills.     

PHOSPHORUS RELEASE AND ASSIMILATORY CAPACITY OF TWO LOWER MISSISSIPPI VALLEY FRESHWATER WETLAND SOILS1

ABSTRACT: Phosphorus fluxes and water quality functions of a bottomland hardwood and freshwater marsh wetland soil were compared. The effect of soil physicochemical conditions, phosphorus loading rate, and diffusive exchange between soils and the overlying food water column on phosphorus release and retention were studied. The predominantly mineral swamp forest soil displayed greater phosphorus sorption potential than the organic freshwater marsh soil. Moreover, due to its low bulk density (0.11 g cm^?3), the freshwater marsh soil surface area required for phosphorus retention is very large compared to the bottomland hardwood wetland soil. For both wetlands, soil redox status affected P release and assimilatory capacity. The more reducing the soils, the smaller their phosphorus retention capacity (greater their release). Phosphorus removal from the overlying water column into the wetland soils followed a first-order kinetic model. Under similar hydrological conditions, phosphorus was found to diffuse 1.2 times faster to the bottom. land hardwood soil than in the freshwater marsh soil. Results indicate that while the bottomland hardwood wetland soil will serve as a sink for phosphorus entering such wetland, phosphorus will be released and exported from the freshwater marsh soil into adjacent ecosystems.     

VARIABILITY OF SOIL WATER PROPERTIES AND CROP YIELD IN A SLOPED WATERSHED1

ABSTRACT: Spatial distribution of soil and water properties and the correlations between them and crop yield were determined for a natural rainfall environment. Hydraulic conductivity, soil texture, water retention, and soil-water flux were variables used to investigate their relationship to crop yield using multiple regression techniques. Variations in crop yields on a watershed with a 3 to 4 percent slope and moderately erosive soils were related to soil-water characteristics and soil properties along slope and with depth. Climatic conditions to sustain crop growth and yield ranged from inadequate soil water in 1983 to adequate soil water in 1984. Crop yield was predicted with models using both available and measured soil-water content. Available water content provided a better model for the prediction of water yield and does not require field measurements of actual soil-water content. Soil water holding capacity was more significant for predicting crop yield in soils with moderate to high silt content than infiltrability of water into the soil.     

ESTIMATING PARAMETERS OF EV1 DISTRIBUTION FOR FLOOD FREQUENCY ANALYSIS1

ABSTRACT: The parameters of the extreme value type 1 distribution were estimated for 55 annual flood data sets by seven methods. These are the methods of (1) moments, (2) probability weighted moments, (3) mixed moments, (4) maximum likelihood estimation, (5) incomplete means, (6) principle of maximum entropy, and (7) least squares. The method of maximum likelihood estimation was found to be the best and the method of incomplete means the worst. The differences between the methods of principle of maximum entropy, probability weighted moments, moments, and least squares were only minor. The difference between these methods and the method of maximum likelihood was not pronounced.     

SOIL MOISTURE SENSORS FOR CONTINUOUS MONITORING1

ABSTRACT: Certain physical and chemical properties of soil vary with soil water content. The relationship between these properties and water content is complex and involves both the pore structure and constituents of the soil solution. One of the most economical techniques to quantify soil water content involves the measurement of electrical resistance of a dielectric medium that is in equilibrium with the soil water content. The objective of this research was to test the reliability and accuracy of fiberglass soil-moisture electrical resistance sensors (ERS) as compared to gravimetric sampling and Time Domain Reflectometry (TDR). The response of the ERS was compared to gravimetric measurements at eight locations on the USDA-ABS Walnut Gulch Experimental Watershed. The comparisons with TDR sensors were made at three additional locations on the same watershed. The high soil rock content (≥45 percent) at seven locations resulted in consistent overestimation of soil water content by the ERS method. Where rock content was less than 10 percent, estimation of soil water was within 5 percent of the gravimetric soil water content. New methodology to calibrate the ERS sensors for rocky soils will need to be developed before soil water content values can be determined with these sensors.     

Investigation of the Curve Number Method For Surface Runoff Estimation In Tropical Regions

This study tests the applicability of the curve number (CN) method within the Soil and Water Assessment Tool (SWAT) to estimate surface runoff at the watershed scale in tropical regions. To do this, surface runoff simulated using the CN method was compared with observed runoff in numerous rainfall‐runoff events in three small tropical watersheds located in the Upper Blue Nile basin, Ethiopia. The CN method generally performed well in simulating surface runoff in the studied watersheds (Nash‐Sutcliff efficiency [NSE] > 0.7; percent bias [PBIAS] < 32%). Moreover, there was no difference in the performance of the CN method in simulating surface runoff under low and high antecedent rainfall (PBIAS for both antecedent conditions: ~30%; modified NSE: ~0.4). It was also found that the method accurately estimated surface runoff at high rainfall intensity (e.g., PBIAS < 15%); however, at low rainfall intensity, the CN method repeatedly underestimated surface runoff (e.g., PBIAS > 60%). This was possibly due to low infiltrability and valley bottom saturated areas typical of many tropical soils, indicating that there is scope for further improvements in the parameterization/representation of tropical soils in the CN method for runoff estimation, to capture low rainfall‐intensity events. In this study the retention parameter was linked to the soil moisture content, which seems to be an appropriate approach to account for antecedent wetness conditions in the tropics.     

Measuring sorption of hydrophilic organic compounds in soils by an unsaturated transient flow method

Determination of sorption of hydrophilic, weakly sorbing organic compounds in soil by conventional batch methods using a slurried suspension is often prone to considerable errors because small changes in the solution concentration on equilibration must be accurately determined. This difficulty is exacerbated for compounds susceptible to degradation, which also decreases the solution concentration. The objective of this study was to determine sorption of hydrophilic pesticides by applying an unsaturated transient flow method, which enables determination of sorption at sufficiently small solution to soil ratios. The method makes use of piston-like displacement of the antecedent solution in equilibrium with sorbed phase when pesticide-free water is infiltrated into a soil column spiked with a pesticide. Pesticide sorption and the solution concentration are inferred from a plot of total pesticide content per unit mass of soil vs. water content in a region where the antecedent solution is accumulated. Thus, extraction of solution from relative dry soil is unnecessary. We tested this method for two hydrophilic pesticides, monocrotophos [dimethyl (E)-1-methyl-2-(methyl-carbamoyl) vinyl phosphate] and dichlorvos (2,2-dichlorovinyl dimethyl phosphate). The sorption coefficient, K(d), obtained for monocrotophos was slightly lower than that by batch method (K(d) = 0.10 vs. 0.19 L kg(-1)), whereas for dichlorvos, a compound highly susceptible to degradation, the unsaturated flow method yielded a much smaller K(d) (0.19 vs. 3.22 L kg(-1)). The K(d) values for both compounds were consistent with the observed retardation in the pesticide displacement in the columns. The proposed method is more representative of field conditions and particularly suitable for weakly sorbing organic compounds in soils.     

Influence of flooding on phosphorus mobility in manure-impacted soil

Agricultural lands are often used for constructing stormwater treatment areas (STAs) to abate nutrient loading to adjacent aquatic systems. Flooding agricultural lands to create STAs could stimulate a significant release of phosphorus (P) from soil to the water column. To assess the suitability of agricultural lands, specifically those impacted by animal operations, for the construction of STAs, soils from different components of the New Palm-Newcomer dairies (Nubbin Slough Basin, Okeechobee, Florida, USA) were collected by horizon and their P retention and release capacities estimated. In general, P released from A-horizon soil under flooded (anaerobic) conditions was greater than under drained (aerobic) conditions due to redox effect on iron (Fe) and consequent P releases. However, the P released from Bh-horizon soil was greater under aerobic conditions than under anaerobic conditions, possibly due to excessive aluminum (Al) content in the horizon. Double acid-extractable calcium (Ca), magnesium (Mg), Al, and P explained 87% of the variability in P release under aerobic conditions, and 80% of that under anaerobic conditions. The P release maxima indicated a high solubility of P in A-horizon soil from both active and abandoned dairies (13 and 8% of the total P, respectively), suggesting that these soils could function as potential sources of P to the overlying water column when used in STA construction. Preestablishment of vegetative communities or chemical amendment, however, could ameliorate high P flux from soil to the water column.     

ESTIMATION OF REGIONAL EVAPOTRANSPIRATION FROM ARID AND SEMI-ARID SURFACES1

ABSTRACT: Evapotranspiration (ET) from oasis and Gobi surfaces in the Heihe River region of China was estimated by Kotoda (1986) and Advection-Aridity (Brutsaert and Stricker, 1979) models. The ET estimates were compared with eddy-correlation ET estimates. The ET estimated by both models are in good agreement with ET obtained by means of eddy-correlation method for the oasis surface while underestimating ET as compared to the eddy-correlation method. For the Gobi surface, the models yielded obviously overestimates of evapotranspiration. In order to estimate evapotranspiration from arid surfaces, the Kotoda model was modified by introducing the surface moisture availability parameter a from Barton (1979). The modified Kotoda model yielded ET estimates that were very close to that from the eddy-correlation method for the Gobi surface. The modified Kotoda model was used to estimate evapotranspiration from the Heihe River watershed, an area with complicated topography and land use, and the results compared with those from a water balance method. A sensitivity analysis of the modified model was performed. The results show that the modified Kotoda model can reflect the relationship between the actual evapotranspiration and the main controlling factors on it for both wet and arid surfaces reasonably. From this study, it can be concluded that the modified Kotoda model is applicable fro the estimation of regional evapotranspiration from areas with complicated topography and land use.     

ARE SOILS LIKE SPONGES?1

ABSTRACT: The purpose of this paper is to explore the validity of the old analogy that “soil is like a sponge.” Laboratory experiments were conducted to measure two hydrologic properties of porous media: drainage under gravity, and water potential curves. A tipping bucket rain gage connected to a data logger was used to measure the rate at which water drained under the force of gravity from a trough filled with four saturated porous media ‐ cellulose sponges, topsoil, peat, and a medium sand. Pressure plate techniques were used to determine water potential curves for soil materials and sponges. In terms of relative cumulative discharge from the trough, sponges were intermediate between peat and topsoil. Because of their tremendous water‐holding capacity, sponges discharged more than 2.5 times as much water as did peat. The water potential curve for sponges was fairly flat, like that of topsoil, but the high water content across all pressures (0.30–15.0 bars) indicated some similarity to peat. The results of these experiments suggest that the general patterns of water retention and release in soil materials and sponges are similar and vary only in degree.     

EVALUATION OF CURVE NUMBER PROCEDURES TO PREDICT RUNOFF IN GLEAMS1

ABSTRACT: Proper selection of curve number values will improve the capability of the SCS-Curve Number procedure in predicting runoff. Both CREAMS and GLEAMS models use the Smith and Williams (1980) approach of converting CNII (curve number value for average antecedent moisture conditions) into CNI (curve number value for dry antecedent moisture conditions) in calculating the soil retention parameter (S). CREAMS and GLEAMS have been found to under predict runoff because of the internal conversion of CNII to CNI. This study shows modifications of the GLEAMS model using CMI without converting it to CM and it also shows the seasonal curve number approaches with and without converting CNII to CNI. Results indicate that using CNII without internal conversion to CNI provides better runoff and erosion predictions than the original version of GLEAMS and versions with seasonal curve numbers when tested with four years of field data in the Coastal Plain physiographic region of Maryland.     

POTENTIAL IMPACT OF IRRIGATION ON WATER RESOURCES IN THE LITTLE WABASH BASIN OF ILLINOIS1

The level of water demand for supplemental irrigation and the impact of such demand on water supplies were estimated, as a function of the price of corn (Zea Mays L.). The method of estimation was based on an economic analysis of irrigation practice which assumed constant irrigation costs, profit maximizing behavior on the part of irrigators, and which was deliberately structured to underestimate the level of irrigation water use. The analysis was applied to and used data from the Little Wabash basin in Illinois. No irrigation was predicted at a corn price below $3.50 per bushel. Between $3.50 and about $6.50 per bushel, irrigation was estimated to be profitable for a small region of the basin where acceptable groundwater was available. Above about $6.50 to $7.50, irrigation was found to be profitable in the remainder of the basin, where impoundment storage was required. The potential impact on the water resources of the basin is significant. For a corn price between $3.50 and about $6.50, the probability of a shortfall, defined as the event where the potential demand exceeds the supply, was estimated to be between 2 percent and 20 percent during the growing season. Above about $7.50, this probability was found to be about one-third. The development of policies to control withdrawals for irrigation and other uses is endorsed.     

SUSCEPTIBILITY OF INDIANA WATERSHEDS TO HERBICIDE CONTAMINATION1

ABSTRACT: An index of watershed susceptibility to surface water contamination by herbicides could be used to improve source water assessments for public drinking water supplies, prioritize watershed restoration projects, and direct funding and educational efforts to areas where the greatest environmental benefit can be realized. The goal of this study is to use streamflow and herbicide concentration data to develop and evaluate a method for estimating comparative watershed susceptibility to herbicide loss. United States Geological Survey (USGS) concentration data for five relatively water soluble herbicides (alachlor, atrazine, cyanazine, metolachlor, and simazine) were analyzed for 16 Indiana watersheds. Correlation was assessed between observed herbicide losses and: (1) a herbicide runoff index using GIS‐based land use, soil type, SCS runoff curve number, tillage practice, herbicide use estimates, and combinations of these factors; and (2) predicted herbicide losses from a non‐point source pollution model (NAPRA‐Web, an Internet‐based interface for GLEAMS). The highest adjusted R²value was found between herbicide concentration and the runoff curve number alone, ranging from 0.25 to 0.56. Predictions from the simulation model showed a poorer correlation with observed herbicide loss. This indicates potential for using the runoff curve number as a simple herbicide contamination susceptibility index.     

ACCUMULATION OF AS,NI, CU,AND PB IN RETENTION AND RECHARGE BASINS SOILS FROM URBAN RUNOFF1

ABSTRACT: The accumulation of arsenic, nickel, copper, and lead in the soil profile was determined beneath five urban storm-water retention/recharge basins used by the Fresno Metropolitan Flood Control District, California. Soils were sampled from the surface to the first zone of saturation and compared with soils from an adjacent un-contaminated control site. These elements were found to be accumulating in the first few centimeters of basin soil and are important to the effectiveness of a specific best management practice, i.e., the retention and recharge of urban storm water. Study basins in use since 1962, 1965, and 1969 had lead contents in the 0–2 cm soil depth interval‘of 570, 670, and 1400 mg Pb/kg soil, respectively. The median indigenous soil lead concentration was 4.6 mg/kg soil. The practice of removing excess flood runoff water from two basins by pumping apparently is a factor in reducing the accumulation rate of these elements in the surface soils of the basins.     

Polyacrylamide distribution in columns of organic matter-removed soils following surface application

Knowledge of how polyacrylamide (PAM) penetrates and distributes in a soil profile after application in irrigation water is important for understanding PAM conditioning depth and evaluating its environmental effects. Little is known, however, about PAM distribution in soil because of the difficulty in quantifying PAM content in natural soils. By using a recently modified substrate-borne PAM quantification method, PAM distribution in columns of organic matter-removed soils was determined. Results showed that penetration of PAM into the soil was affected by salt level of irrigation water, soil texture, initial soil water content, water application method, and other factors. Polyacrylamide penetration depth was about one-eighth to one-half of the water penetration depth, with a particularly high PAM retention in the top few centimeters of the soil. Under different experimental conditions, the PAM retained in the top 0 to 2 cm of soil ranged from 16 to 95% of the total applied amount. More favorable solution-soil contact conditions, longer solution-soil contact time, and lower initial soil moisture caused much more PAM retention in the top few centimeters of the soil. High sorptive affinity of PAM on soil is the main reason for its low penetration into the soil. Although these results were not obtained from natural soils, they are still helpful in improving our understanding of PAM transport behavior in soils.     

Effects of calibration on L-THIA GIS runoff and pollutant estimation

Urbanization can result in alteration of a watershed's hydrologic response and water quality. To simulate hydrologic and water quality impacts of land use changes, the Long-Term Hydrologic Impact Assessment (L-THIA) system has been used. The L-THIA system estimates pollutant loading based on direct runoff quantity and land use based pollutant coefficients. The accurate estimation of direct runoff is important in assessing water quality impacts of land use changes. An automated program was developed to calibrate the L-THIA model using the millions of curve number (CN) combinations associated with land uses and hydrologic soil groups. L-THIA calibration for the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana was performed using land use data for 1991 and daily rainfall data for six months of 1991 (January 1-June 30) to minimize errors associated with use of different temporal land use data and rainfall data. For the calibration period, the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The calibrated CN values were used for validation of the model for the same year (July 1-December 31), and the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The Nash-Sutcliffe coefficient was 0.52 for January 1, 1991 to December 31, 1991 using uncalibrated CN values. As shown in this study, the use of better input parameters for the L-THIA model can improve accuracy. The effects on direct runoff and pollutant estimation of the calibrated CN values in the L-THIA model were investigated for the LEC. Following calibration, the estimated average annual direct runoff for the LEC watershed increased by 34%, total nitrogen by 24%, total phosphorus by 22%, and total lead by 43%. This study demonstrates that the L-THIA model should be calibrated and validated prior to application in a particular watershed to more accurately assess the effects of land use changes on hydrology and water quality.     

EFFECT OF SERIAL CORRECTION ON GROUND WATER QUALITY SAMPLING FREQUENCY1

ABSTRACT: Methods of calculating uncertainty in estimates of serial correlation coefficients, and correcting for bias in short and medium length (less than 50 data point) records, are presented. Uncertainty and bias in the estimation of serial correlation coefficients for ground water quality data is shown to be considerable and to result in inaccurate calculation of the sampling frequencies for monitoring purposes. The methods are applied to a ground water data set consisting of 87 monthly measurements of nitrate concentrations. The variation in serial correlation coefficients with variation of record length is examined. The optimum sampling frequencies for detection of changes in ground water nitrate concentrations are estimated.