SATELLITE MICROWAVE OBSERVATIONS OF SOIL MOISTURE VARIATIONS1

ABSTRACT: Results from studies in the Illinois-Indiana and Texas-Oklahoma areas indicate that satellite microwave observations at the 1.55 cm wavelength are responsive to relative moisture variations in the near surface layer of the soil. Because significant vegetation cover absorbs the 1.55 cm microwave emission from the soil, the target area must be predominately bare soil or low density vegetation cover for meaningful measurements to result. The 25 km resolution of the satellite sensor limits application of the microwave techniques to large areas such as watersheds or agricultural districts rather than individual fields. In general, at 1.55 cm. there is an inverse relationship between microwave brightness temperature and changes in soil moisture levels (as indicated by antecedent rainfall) in agricultural regions before the planting of crops or during the early growing season when vegetation cover is sparse. Even early season observations should be of great value in deciding on the time and type of crop planting and for initial irrigation scheduling when the root zone is still in close proximity to the surface.     

SIMULATING SOIL WATER RECESSION COEFFICIENTS FOR AGRICULTURAL WATERSHEDS1

ABSTRACT: The antecedent precipitation index (API) has been a useful indicator of soil moisture conditions for watershed runoff calculations, and recent attempts to correlate this index with spaceborne microwave observations have been fairly successful. The prognostic equation for soil moisture used in some of the atmospheric general circulation models (GCM) together with Thomthwaite-Mather parameterization of actual evapotranspiration leads to API equations. The recession coefficient for API is found to depend on climatic factors as contained in potential evapotranspiration and to depend on soil texture as reflected by field capacity and permanent wilting point. A recently developed model for global insolation is used with climatological data for Wisconsin to simulate the annual trend of the recession coefficient. Good quantitative agreement is shown with the observed trends at Fennimore and Colby watersheds in Wisconsin. This study suggests that API could be a unifying concept for watershed and atmospheric general circulation modeling.     

MODIS Biophysical States and NEXRAD Precipitation in a Statistical Evaluation of Antecedent Moisture Condition and Streamflow1

Abstract: The potential of remotely sensed time series of biophysical states of landscape to characterize soil moisture condition antecedent to radar estimates of precipitation is assessed in a statistical prediction model of streamflow in a 1,420 km² watershed in south‐central Texas, Moderate Resolution Imaging Spectroradiometer (MODIS) time series biophysical products offer significant opportunities to characterize and quantify hydrologic state variables such as land surface temperature (LST) and vegetation state and status. Together with Next Generation Weather Radar (NEXRAD) precipitation estimates for the period 2002 through 2005, 16 raw and deseasoned time series of LST (day and night), vegetation indices, infrared reflectances, and water stress indices were linearly regressed against observed watershed streamflow on an eight‐day aggregated time period. Time offsets of 0 (synchronous with streamflow event), 8, and 16 days (leading streamflow event) were assessed for each of the 16 parameters to evaluate antecedent effects. The model results indicated a reasonable correlation (r² = 0.67) when precipitation, daytime LST advanced 16 days, and a deseasoned moisture stress index were regressed against log‐transformed streamflow. The estimation model was applied to a validation period from January 2006 through March 2007, a period of 12 months of regional drought and base‐flow conditions followed by three months of above normal rainfall and a flood event. The model resulted in a Nash‐Sutcliffe estimation efficiency (E) of 0.45 for flow series (in log‐space) for the full 15‐month period, ?0.03 for the 2006 drought condition period, and 0.87 for the 2007 wet condition period. The overall model had a relative volume error of ?32%. The contribution of parameter uncertainties to model discrepancy was evaluated.     

IMPROVED RAINFALL/RUNOFF ESTIMATES USING REMOTELY SENSED SOIL MOISTURE1

ABSTRACT: Remotely sensed soil moisture data measured during the Southern Great Plains 1997 (SGP97) experiment in Oklahoma were used to characterize antecedent soil moisture conditions for the Soil Conservation Service (SCS) curve number method. The precipitation‐adjusted curve number and the soil moisture were strongly related (r²= 0.70). Remotely sensed soil moisture fields were used to adjust the curve numbers and the runoff estimates for five watersheds, in the Little Washita watershed; the results ranged from 2.8 km² to 601.6 km². The soil moisture data were applied at two spatial scales, a finer one (800 m) measuring spatial resolution and a coarser one (28 km). The root mean square error (RMSE) and the mean absolute error (MAE) of the runoff estimated by the standard SCS method was reduced by nearly 50 percent when the 800 m soil moisture data were used to adjust the curve number. The coarser scale soil moisture data also significantly reduced the error in the runoff predictions with 41 percent and 28 percent reductions in MAE and RMSE, respectively. The results suggest that remote sensing of soil moisture, when combined with the SCS method, can improve rainfall runoff predictions at a range of spatial scales.     

Stochastic Transfer Function Model and Neural Networks to Estimate Soil Moisture1

Abstract: A practical methodology is proposed to estimate the three‐dimensional variability of soil moisture based on a stochastic transfer function model, which is an approximation of the Richard’s equation. Satellite, radar and in situ observations are the major sources of information to develop a model that represents the dynamic water content in the soil. The soil‐moisture observations were collected from 17 stations located in Puerto Rico (PR), and a sequential quadratic programming algorithm was used to estimate the parameters of the transfer function (TF) at each station. Soil texture information, terrain elevation, vegetation index, surface temperature, and accumulated rainfall for every grid cell were input into a self‐organized artificial neural network to identify similarities on terrain spatial variability and to determine the TF that best resembles the properties of a particular grid point. Soil moisture observed at 20 cm depth, soil texture, and cumulative rainfall were also used to train a feedforward artificial neural network to estimate soil moisture at 5, 10, 50, and 100 cm depth. A validation procedure was implemented to measure the horizontal and vertical estimation accuracy of soil moisture. Validation results from spatial and temporal variation of volumetric water content (vwc) showed that the proposed algorithm estimated soil moisture with a root mean squared error (RMSE) of 2.31% vwc, and the vertical profile shows a RMSE of 2.50% vwc. The algorithm estimates soil moisture in an hourly basis at 1 km spatial resolution, and up to 1 m depth, and was successfully applied under PR climate conditions.     

PREPLANTING SOIL MOISTURE USING PASSIVE MICROWAVE SENSORS1

ABSTRACT: Accurate assessment of preplanting soil moisture conditions is necessary for good agricultural management, and can have a significant influence on crop yield in the Texas Panhandle region. The Texas High Plains Underground Water Conservation District invests considerable time and money in developing a soil moisture deficit map each year in the hopes of achieving optimal use of irrigation water. Microwave sensors are responsive to surface soil moisture and, if used in this application, can provide timely and detailed information on root zone soil moisture. For this reason, an experiment was conducted in 1984 to evaluate the potential of aircraft-mounted passive microwave sensors. Microwave radiometer data were collected over a 2700 km² area near Lubbock, Texas, with a processed resolution of 0.32 km². These data were ground registered and converted to estimates of soil moisture using an appropriate model and land cover and soil texture information. Analyses indicate that the system provides an efficient means for mapping variations in soil moisture over large areas.     

Sensitivity of SCS Models to Curve Number Variation1

ABSTRACT: The Soil Conservation Service (SCS) models, including the TR-20 computer program and the simplified methods in TR-55, are widely used in hydrologic design. The runoff curve number (CN), which is an important input parameter to SCS models, is defined in terms of land use tretments, hydrologic, condition, antecedent soil moisture, and soil type. The objective of this study was to evaluate the sensitivity of the SCS models to errors in CN estimates. The results show that the effects of CN variation decrease as the design rainfall depth increases, such as for the larger storm events. The value and use of the sensitivity curves are demonstrated using a comparison of Landsat and conventionally derived curve numbers for three watersheds in Pennsylvania.     

USE OF SOIL MOISTURE PROBES TO ESTIMATE GROUND WATER RECHARGE AT AN OIL SPILL SITE1

Soil moisture data collected using an automated data logging system were used to estimate ground water recharge at a crude oil spill research site near Bemidji, Minnesota. Three different soil moisture probes were tested in the laboratory as well as the field conditions of limited power supply and extreme weather typical of northern Minnesota: a self‐contained reflectometer probe, and two time domain reflectometry (TDR) probes, 30 and 50 cm long. Recharge was estimated using an unsaturated zone water balance method. Recharge estimates for 1999 using the laboratory calibrations were 13 to 30 percent greater than estimates based on the factory calibrations. Recharge indicated by the self‐contained probes was 170 percent to 210 percent greater than the estimates for the TDR probes regardless of calibration method. Results indicate that the anomalously large recharge estimates for the self‐contained probes are not the result of inaccurate measurements of volumetric moisture content, but result from the presence of crude oil, or borehole leakage. Of the probes tested, the 50 cm long TDR probe yielded recharge estimates that compared most favorably to estimates based on a method utilizing water table fluctuations. Recharge rates for this probe represented 24 to 27 percent of 1999 precipitation. Recharge based on the 30 cm long horizontal TDR probes was 29 to 37 percent of 1999 precipitation. By comparison, recharge based on the water table fluctuation method represented about 29 percent of precipitation.     

DEFINING WATERSHED-SCALE EVAPORATION USING A NORMALIZED DIFFERENCE VEGETATION INDEX1

Monthly composites of the Normalized Difference Vegetation Indices (NDVI), derived from the National Oceanic and Atmospheric Administration's (NOAA) Advanced Very High Resolution Radiometer (AVILRR), were transformed linearly into monthly evaporation rates and compared with detailed hydrologic-model simulation results for five watersheds across the United States. Model-simulated monthly evaporation values showed high correlations (mean R²= .77) with NDVI-derived evaporation estimates. These latter estimates, used in a classical water balance model, resulted in equally accurate simulations of monthly runoff than when the model was run to estimate monthly evaporation via soil moisture accounting. Comparison of NDVI-derived evaporation estimates with pan data showed promise for transforming NDVI values into evaporation estimates under both wet and water-limiting conditions without resorting to the application of any kind of calibrated hydrologic models.     

WATER LOSS EQUATIONS AND COEFFICIENTS TO ESTIMATE RUNOFF FROM RAINFALL1

Infiltration models are based on physical characteristics of the soil and initial soil moisture. For a given soil it is based on the initial soil moisture distribution. A computer simulation model for flood runoff systems (FH-Model) was used to analyze 39 sets of rainfall-runoff data on four small watersheds ranging in size from 17 to 342 square kilometers located in the Yamaska River basin in Quebec. From these analyses, parameters and coefficients have been determined for a water loss (infiltration) equation. A method for determining the loss parameters, using a nonlinear least square curve fitting technique, is presented. Expressions were made to relate the loss parameters to antecedent precipitation. The equations were tested on 11 storm rainfall and runoff events on a watershed located in the same region and close agreements were found.     

Comparison of Field Performance of Multiple Soil Moisture Sensors in a Semi‐Arid Rangeland1

Abstract:  Automated electronic soil moisture sensors, such as time domain reflectometry (TDR) and capacitance probes are being used extensively to monitor and measure soil moisture in a variety of scientific and land management applications. These sensors are often used for a wide range of soil moisture applications such as drought forage prediction or validation of large‐scale remote sensing instruments. The convergence of three different research projects facilitated the evaluation and comparison of three commercially available electronic soil moisture probes under field application conditions. The sensors are all installed in shallow soil profiles in a well instrumented small semi‐arid shrub covered subwatershed in Southeastern Arizona. The sensors use either a TDR or a capacitance technique; both of which indirectly measure the soil dielectric constant to determine the soil moisture content. Sensors are evaluated over a range of conditions during three seasons comparing responses to natural wetting and drying sequences and using water balance and infiltration simulation models. Each of the sensors responded to the majority of precipitation events; however, they varied greatly in response time and magnitude from each other. Measured profile soil moisture storage compared better to water balance estimates when soil moisture in deeper layers was accounted for in the calculations. No distinct or consistent trend was detected when comparing the responses from the sensors or the infiltration model to individual precipitation events. The results underscore the need to understand how the sensors respond under field application and recognize the limitations of soil moisture sensors and the factors that can affect their accuracy in predicting soil moisture in situ.     

SOIL MOISTURE VARIATION PATTERNS OBSERVED IN HAND COUNTY,SOUTH DAKOTA1

ABSTRACT: Sail moisture data were taken during nine sampling events (1976-1978) at a test site in South Dakota as part of the ground truth used in NASA's aircraft experiments studying the microwave sensing of soil moisture. This portion of the study dealt only with the spatial variability observed with regard to the ground data. Samples were taken over three surface depths at each point, and the data reported as the mean field moisture content within each of three surface horizons. The results shed additional light on the relationship between ground sampling and remote sensing of soil moisture. First, it was found that it is best to partition data of well drained sites from poorly drained areas when attempting to characterize the surface moisture content throughout an area of varying soil and cover conditions. It was also found that the moisture coefficient of variation within a field decreased as the mean field soil moisture increased, and that the standard deviation was at a maximum in the mid-range of observed moisture conditions (15-25 percent). Within field sample variation also decreases as the sample is integrated over a greater surface depth. It was determined that a sampling intensity of 10 samples per kilometer was adequate to characterize the mean field soil moisture at all three depths along a transect in the areas of moderate to good drainage-.     

A REVISED VERSION OF PnET‐II TO SIMULATE THE HYDROLOGIC CYCLE IN SOUTHEASTERN FORESTED AREAS1

ABSTRACT: The PnET‐II model uses hydroclimatic data on maximum and minimum temperatures, precipitation, and solar radiation, together with vegetation and soil parameters, to produce estimates of net primary productivity, evapotranspiration (ET), and runoff on a monthly time step for forested areas. In this study, the PnET‐II model was employed to simulate the hydrologic cycle for 17 Southeastern eight‐digit hydrologic unit code (HUC) watersheds dominated by evergreen or deciduous tree species. Based on these control experiments, model biases were quantified and tentative revision schemes were introduced. Revisions included: (1) replacing the original single soil layer with three soil layers in the water balance routine; (2) introducing calibrating factors to rectify the phenomenon of overestimation of ET in spring and early summer months; (3) parameterizing proper values of growing degree days for trees located in different climate zones; and (4) adjusting the parameter of fast‐flow (overland flow) fraction based on antecedent moisture condition and precipitation intensity. The revised PnET‐II model, called PnET‐II3SL in this work, substantially improved runoff simulations for the 17 selected experimental sites, and therefore may offer a more powerful tool to address issues in water resources management.     

ENERGY BALANCE OF A CORN RESIDUE-COVERED FIELD DURING SNOWMELT1

ABSTRACT: Transport of agricultural chemicals in runoff and recharge waters from snowmelt and soil thawing may represent a significant event in terms of annual contaminant loadings in temperate regions. Improved understanding of the melt dynamics of shallow snowpacks is necessary to fully assess the implications for water quality. The objective of this study was to measure the energy balance components of a corn (Zea mays L.) stubble field during the melting of its snowcover. Net radiation (Rn), soil (G), sensible (H), and latent (Q) heat fluxes were measured in a field near Ames, Iowa, during the winter of 1994–1995. Energy consumed by melting including change in energy storage of the snowpack was determined as the residual of the measured energy balance. There was continuous snowcover at the field site for 71 days (maximum depth = 222 mm) followed by an open period of 11 days before additional snowfall and a second melt period. The net radiation and snow melt/energy storage change (5) terms dominated the energy balance during both measurement intervals. Peak daily sensible and latent heat fluxes were below 100 W m^?2 on all days except the last day of the second melt period. There was good agreement between predicted and measured values of H and Q during the melting of an aged snow layer but poorer agreement during the melt of fresh snow. Both snowpacks melted rapidly and coincident changes in soil moisture storage were observed. Improved estimates of Q and H, especially for partially open surfaces, will require better characterization of the surface aerodynamic properties and spatially-representative surface temperature measurements.     

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.     

Estimating Root Zone Soil Moisture at Continental Scale Using Neural Networks

This study investigates the feasibility of artificial neural networks (ANNs) to retrieve root zone soil moisture (RZSM) at the depths of 20 cm (SM20) and 50 cm (SM50) at a continental scale, using surface information. To train the ANNs to capture interactions between land surface and various climatic patterns, data of 557 stations over the continental United States were collected. A sensitivity analysis revealed that the ANNs were able to identify input variables that directly affect the water and energy balance in root zone. The data important for RZSM retrieval in a large area included soil texture, surface soil moisture, and the cumulative values of air temperature, surface soil temperature, rainfall, and snowfall. The results showed that the ANNs had high skill in retrieving SM20 with a correlation coefficient above 0.7 in most cases, but were less effective at estimating SM50. The comparison of the ANNs showed that using soil texture data improved the model performance, especially for the estimation of SM50. It was demonstrated that the ANNs had high flexibility for applications in different climatic regions. The method was used to generate RZSM in North America using Soil Moisture and Ocean Salinity (SMOS) soil moisture data, and achieved a spatial soil moisture pattern comparable to that of Global Land Data Assimilation System Noah model with comparable performance to the SMOS surface soil moisture retrievals. The models can be efficient alternatives to assimilate remote sensing soil moisture data for shallow RZSM retrieval.     

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.     

RESPONSE OF GREEN‐AMPT MEIN‐LARSEN SIMULATED RUNOFF VOLUMES TO TEMPORALLY AGGREGATED PRECIPITATION1

ABSTRACT: Data collection frequency in automated systems is user determined and can range from seconds to hours or days. Currently, there is no standard or recommended frequency interval for collecting precipitation data from automated systems for input to event‐based models such as Green‐Ampt Mein‐Larsen (GAML). Data from 47 storm events at seven locations were used to simulate the response of GAML excess rainfall to temporally aggregated precipitation data. No difference in model efficiency was recognized when comparing one‐minute interval data (R²= 1.00) to five‐minute data (R²= 1.00). Very little model efficiency was lost at a 10‐minute (R²= 0.96) interval. After 10‐minutes, decline in efficiency became more rapid with R²= 0.16 at one hour. The combined effect of time interval with respect to drainage area, hydraulic conductivity, maximum 30‐minute intensity, and total precipitation also revealed similar results.     

PRECIPITATION DIFFERENCES AMONGST GCMs USED FOR THE U.S. NATIONAL ASSESSMENT1

ABSTRACT: Two general circulation models (GCMs) used in the U.S. national assessment of the potential consequences of climate variability and change (CGCM1 and HadCM2) show a large increase in precipitation in the future over the southwestern U.S., particularly during winter. This precipitation increase is an extension of a larger region of increased precipitation in the Pacific Ocean off the west coast of North America that is associated with a deepened and southward-shifted Aleutian Low, a weaker subtropical high, and warmer sea surface temperatures (SSTs). The models differ in their simulation of precipitation anomalies over the southeastern U.S., with CGCM1 showing drier conditions and HadCM2 showing wetter conditions in the future. While both models show decreased frequency of Atlantic storms, consistent with decreased meridional and land/sea temperature gradients, the more coastal position of the storm track in CGCM1 results in less precipitation than modern along the eastern seaboard of the U.S. During summer, differences in land surface models within the two GCMs sometimes lead to differences in soil moisture that feed back to the precipitation over land due to available moisture.     

EVAPOTRANSPIRATION AND SOIL MOISTURE DYNAMICS ON A SEMIARID PONDEROSA PINE HILLSLOPE1

ABSTRACT: Detailed measurements of soil moisture and ET in semiarid forest environments have not been widely reported in the literature. In this study, soil moisture and water balance components were measured over a four‐year period on a semiarid ponderosa pine hillslope, with evapotranspiration (ET) determined as the residual of measured precipitation, runoff, and change in soil moisture storage. ET accounts for approximately 95 percent of the water budget and has a distinctly bimodal annual pattern, with peaks occurring after spring snowmelt and during the late summer monsoon season, periods that coincide with high soil moisture. Weekly growing season ET rates determined by the hillslope water balance are found to be invariably below calculated potential rates. Normalized ET rates are linearly correlated (r²= 0.62) with soil moisture; therefore, a simple linear relation is proposed. Growing season soil moisture dynamics were modeled based on this relation. Results are in fair agreement (r²= 0.63) with the observed soil moisture data over the four growing seasons; however, for two dry summers with little surface runoff, much better results (r² > 0.90) were obtained.