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.     

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-.     

Spatial Dynamics of Soil Moisture in a Complex Terrain in the Semi‐Arid Loess Plateau Region,China1

Abstract: Soil moisture is an important hydrological variable in reforestation practices in a water‐limited region of the Loess Plateau of northwestern China. The objective of this study was to quantify the spatial dynamics of soil moisture on a complex terrain. During 2004‐2006, a total of 313 sample points in two kinds of grid (2 × 2 m and 20 × 20 m) were arranged for soil moisture measurements (two soil layers: 0‐30 and 30‐60 cm) with Time Domain Reflectometry. The geostatistical properties of soil moisture patterns, the variance and correlation structure of the soil moisture, and the effects of terrain factors on soil moisture were analyzed. The results suggested that our sampling grid captured the spatial variability of soil moisture distributions for this complex terrain. Principal Component Analysis and Cluster Analysis statistics showed that soil moisture decreased as slope gradient increased; that sunny aspects (112.5°‐292.5°) had relatively lower soil moisture than did shady aspects (292.5°‐112.5°); that soil moisture was lowest in the SWW direction and highest in the NWN direction; and that hillslope aspect was the main factor affecting soil moisture in the 0‐ to 30‐cm soil layer, whereas the main factor for the 30‐ to 60‐cm layer was slope gradient. It was found that the relative values of soil moisture for steep slopes (>36%) with shady aspect (292.5°‐112.5°), gentle slopes (<36%) with sunny aspect (112.5°‐292.5°), and steep slopes with sunny aspect were 99, 82, and 80, respectively – assuming a soil moisture value of 100 for gentle slopes with shady aspect. The results of this study are expected to be relevant to and useful for reforestation planning and design, parameterization of distributed hydrology models, and land productivity assessment in the study region.     

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.     

SPATIAL DISTRIBUTION OF POINT SOIL MOISTURE ESTIMATES USING LANDSAT TM DATA AND FUZZY-C CLASSIFICATION1

ABSTRACT: Many hydrologic models have input data requirements that are difficult to satisfy for all but a few well-instrumented, experimental watersheds. In this study, point soil moisture in a mountain watershed with various types of vegetative cover was modeled using a generalized regression model. Information on sur-ficial characteristics of the watershed was obtained by applying fuzzy set theory to a database consisting of only satellite and a digital elevation model (DEM). The fuzzy-c algorithm separated the watershed into distinguishable classes and provided regression coefficients for each ground pixel. The regression model used the coefficients to estimate distributed soil moisture over the entire watershed. A soil moisture accounting model was used to resolve temporal differences between measurements at prototypical measurement sites and validation sites. The results were reasonably accurate for all classes in the watershed. The spatial distribution of soil moisture estimates corresponded accurately with soil moisture measurements at validation sites on the watershed. It was concluded that use of the regression model to distribute soil moisture from a specified number of points can be combined with satellite and DEM information to provide a reasonable estimation of the spatial distribution of soil moisture for a watershed.     

Determination of Frequency for Remeasuring Ground and Vegetation Cover Factor Needed for Soil Erosion Modeling

Determining a remeasurement frequency of variables over time is required in monitoring environmental systems. This article demonstrates methods based on regression modeling and spatio-temporal variability to determine the time interval to remeasure the ground and vegetation cover factor on permanent plots for monitoring a soil erosion system. The spatio-temporal variability methods include use of historical data to predict semivariograms, modeling average temporal variability, and temporal interpolation by two-step kriging. The results show that for the cover factor, the relative errors of the prediction increase with an increased length of time interval between remeasurements when using the regression and semivariogram models. Given precision or accuracy requirements, appropriate time intervals can be determined. However, the remeasurement frequency also varies depending on the prediction interval time. As an alternative method, the range parameter of a semivariogram model can be used to quantify average temporal variability that approximates the maximum time interval between remeasurements. This method is simpler than regression and semivariogram modeling, but it requires a long-term dataset based on permanent plots. In addition, the temporal interpolation by two-step kriging is also used to determine the time interval. This method is applicable when remeasurements in time are not sufficient. If spatial and temporal remeasurements are sufficient, it can be expanded and applied to design spatial and temporal sampling simultaneously.     

Soil Moisture Data Assimilation Using Support Vector Machines and Ensemble Kalman Filter1

Abstract: A hybrid data assimilation (DA) methodology that combines two state‐of‐the‐art techniques, support vector machines (SVMs) and ensemble Kalman filter (EnKF), is applied for soil moisture DA in this work. The SVM methodology provides a statistically sound and robust approach to solving the inverse problem, and thus to building statistical models. EnKF is an extension of the Kalman Filter (KF), a well‐known tool in prediction updating. In the present research, ground measurements were used to build a SVM‐type soil moisture predictor. Subsequent observations and their statistics were assimilated to update predictions from the SVM model by coupling it with EnKF. In this way, both model predictions and ground data, as well as their statistics, are fused thus minimizing the prediction error and making the predictions and observations statistically consistent. The results are shown for two approaches; one in which update is done at every time step and the other which assumes that data is only available at alternate time steps (in window of 10 time steps) and hence update is performed at those occasions. The SVM‐EnKF coupling is shown to improve soil moisture forecasts in an example using data from the Soil Climate Analysis Network site at Ames, Iowa.     

KALMAN FILTER CALCULATION OF SAMPLING FREQUENCY WHEN DETERMINING ANNUAL MEAN SOLUTE CONCENTRATIONS

ABSTRACT: A linear filter (Kalman filter) technique was used with a Streamflow-concentration model the minimize surface water quality sampling frequencies when determining annual mean solute concentrations with a predetermined allowable error. The Kalman filter technique used the stream discharge interval as a replacement for the more commonly used time interval. Using filter computations, the measurement error variance was minimized within the sample size constraints. The Kalman filter application proposed here is applicable only under several conditions including: monitoring is solely to estimate annual mean concentration; discharge measurement errors are negligible; the Streamflow-concentration model is valid; and monthly samples reflect the total variance of the solute in question.     

INFLUENCE OF DIFFERENT TEMPORAL SAMPLING STRATEGIES ON ESTIMATING TOTAL PHOSPHORUS AND SUSPENDED SEDIMENT CONCENTRATION AND TRANSPORT IN SMALL STREAMS1

ABSTRACT: Various temporal sampling strategies are used to monitor water quality in small streams. To determine how various strategies influence the estimated water quality, frequently collected water quality data from eight small streams (14 to 110 km²) in Wisconsin were systematically subsampled to simulate typically used strategies. These subsets of data were then used to estimate mean, median, and maximum concentrations, and with continuous daily flows used to estimate annual loads (using the regression method) and volumetrically weighted mean concentrations. For each strategy, accuracy and precision in each summary statistic were evaluated by comparison with concentrations and loads of total phosphorus and suspended sediment estimated from all available data. The most effective sampling strategy depends on the statistic of interest and study duration. For mean and median concentrations, the most frequent fixed period sampling economically feasible is best. For maximum concentrations, any strategy with samples at or prior to peak flow is best. The best sampling strategy to estimate loads depends on the study duration. For one‐year studies, fixed period monthly sampling supplemented with storm chasing was best, even though loads were overestimated by 25 to 50 percent. For two to three‐year load studies and estimating volumetrically weighted mean concentrations, fixed period semimonthly sampling was best.     

Calibrating a Basin‐Scale Groundwater Model to Remotely Sensed Estimates of Groundwater Evapotranspiration

Remotely sensed vegetation indices correspond to canopy vigor and cover and have been successfully used to estimate groundwater evapotranspiration (ET_g) over large spatial and temporal scales. However, these data do not provide information on depth to groundwater (dtgw) necessary for groundwater models (GWM) to calculate ET_g. An iterative approach is provided that calibrates GWM to ET_g derived from Landsat estimates of the Enhanced Vegetation Index (EVI). The approach is applied to different vegetation groups in Mason Valley, Nevada over an 11‐year time span. An uncertainty analysis is done to estimate the resulting mean and 90% confidence intervals in ET_g to dtgw relationships to quantify errors associated with plant physiologic complexity, species variability, and parameter smoothing to the 100 m GWM‐grid, temporal variability in soil moisture and nonuniqueness in the solution. Additionally, a first‐order second moment analysis shows ET_g to dtgw relationships are almost exclusively sensitive to estimated land surface, or maximum, ET_g despite relatively large uncertainty in extinction depths and hydraulic conductivity. The EVI method of estimating ET_g appears to bias ET_g during years with exceptionally wet spring/summer conditions. Excluding these years improves model performance significantly but highlights the need to develop a methodology that accounts not only on quantity but timing of annual precipitation on phreatophyte greenness.     

SAMPLING FRAME FOR IMPROVING PEBBLE COUNT ACCURACY IN COARSE GRAVEL‐BED STREAMS1

ABSTRACT: Improved sampling techniques are needed to increase the accuracy of pebble‐count particle‐size distributions used for stream studies in gravel‐bed streams. However, pebble counts are prone to operator errors introduced through subjective particle selection, serial correlation, and inaccurate particle‐size measurements. Errors in particle‐size measurements can be minimized by using a gravel template. Operator influence on particle selection can be minimized by using a sampling frame, 60 by 60 cm, in which sampling points are identified by the cross points of thin elastic bands. Serial correlation can be minimized by adjusting the spacing between the cross points and setting it equal to the dominant large particle size (=D₉₅). In a field test in a cobble‐bed stream, the sampling frame developed in this study produced slightly coarser size distributions, particularly in the cobble range, than the traditional heel‐to‐toe walk that selects particles with a blind touch at the tip of the boot. The sampling frame produced more similar sampling results between two operators than heel‐to‐toe walks. The difference between the two sampling methods is attributed to an unbiased selection of fine and coarse particles when using the sampling frame.     

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.     

Investigating the Fate and Transport of Escherichia coli in the Charles River,Boston, Using High‐Resolution Observation and Modeling1

Abstract: The processes affecting the fate and transport of Escherichia coli in surface waters were investigated using high‐resolution observation and modeling. The concentration patterns in Boston’s Charles River were observed during four sampling events with a total of 757 samples, including two spatial surveys with two along‐river (1,500 m length) and three across‐river (600 m length) transects at approximately 25‐m intervals, and two temporal surveys at a fixed location (Community Boating) over seven days at hourly intervals. The data reveal significant spatial and temporal structure at scales not resolved by typical monitoring programs. A mechanistic, time‐variable, three‐dimensional coupled hydrodynamic and water quality model was developed using the ECOMSED and RCA modeling frameworks. The computational grid consists of 3,066 grid cells with average length dimension of 25 m. Forcing functions include upstream and downstream boundary conditions, Stony Brook, and Muddy River (major tributaries) combined sewer overflow (CSO) and non‐CSO discharge and wind. The model generally reproduces the observed spatial and temporal patterns. This includes the presence and absence of a plume in the study area under similar loading, but different hydrodynamic conditions caused by operation of the New Charles River Dam (downstream) and wind. The model also correctly predicts an episode of high concentrations at the time‐series station following seven days of no rainfall. The model has an overall root mean square error (RMSE) of 250 CFU/100 ml and an error rate (above or below the USEPA‐recommended single sample criteria value of 235 CFU/100 ml) of 9.4%. At the time series station, the model has an RMSE of 370 CFU/100 ml and an error rate of 15%.     

WATER QUALITY SAMPLING SCHEMES FOR VARIABLE FLOW CANALS AT REMOTE SITES1

ABSTRACT: Growing interest in water quality has resulted in the development of monitoring networks and intensive sampling for various constituents. Common purposes are regulatory, source and sink understanding, and trend observations. Water quality monitoring involves monitoring system design; sampling site instrumentation; and sampling, analysis, quality control, and assurance. Sampling is a process to gather information with the least cost and least error. Various water quality sampling schemes have been applied for different sampling objectives and time frames. In this study, a flow proportional composite sampling scheme is applied to variable flow remote canals where the flow rate is not known a priori. In this scheme, historical weekly flow data are analyzed to develop high flow and low flow sampling trigger volumes for auto‐samplers. The median flow is used to estimate low flow sampling trigger volume and the five percent exceedence probability flow is used for high flow sampling trigger volume. A computer simulation of high resolution sampling is used to demonstrate the comparative bias in load estimation and operational cost among four sampling schemes. Weekly flow proportional composite auto‐sampling resulted in the least bias in load estimation with competitive operational cost compared to daily grab, weekly grab sampling and time proportional auto‐sampling.     

Simulating Spatial Pattern and Dynamics of Military Training Impacts for Allocation of Land Repair Using Images

The land management of US Army installations requires information on land conditions and their history for planning future military training activities and allocation of land repair. There is thus a strong need for methodology development to estimate the land conditions and cumulative military training impacts for the purpose of repair and restoration. In this study, we simulated at Fort Riley, USA, spatial patterns and temporal dynamics of military training impacts on land conditions quantified as percent ground cover using an image-aided spatial conditional co-simulation algorithm. Moreover, we estimated the historical percent ground cover as a measure of the cumulative impacts, and then calculated the allocation of land repair and restoration based on both current and historical land conditions. In addition, we developed a loss function method for allocation of land repair and restoration. The results showed: (1) this co-simulation algorithm reproduced spatial and temporal variability of percent ground cover and provided estimates of uncertainties with the correlation coefficients and root mean square errors between the simulated and observed values varying from 0.63 to 0.88 and from 23% to 78%, respectively; (2) with and without the cumulative impacts, the obtained spatial patterns of the land repair categories were similar, but their land areas differed by 5% to 40% in some years; (3) the combination of the loss function with the co-simulation made it possible to estimate and computationally propagate the uncertainties of land conditions into the uncertainties of expected cost loss for misallocation of land repair and restoration; and (4) the loss function, physical threshold, and probability threshold methods led to similar spatial patterns and temporal dynamics of the land repair categories, however, the loss function increased the land area by 5% to 30% for intense and moderate repairs and decreased the area by 5% to 30% for no repairs and light repairs for most of the years. This approach provided the potential to improve and automate the existing land rehabilitation and maintenance (LRAM) system used for the land management of the U.S. Army installations, and it can be applied to the management of other civil lands and environments. In conclusion, this study overcame the important gaps that exist in the methodological development and application for simulating land conditions and cumulative impacts due to human activities, and also in the methods for the allocation of land for repair and restoration.     

PREDICTING FECAL COLIFORM BACTERIA LEVELS IN THE CHARLES RIVER,MASSACHUSETTS, USA1

In Massachusetts, the Charles River Watershed Association conducts a regular water quality monitoring and public notification program in the Charles River Basin during the recreational season to inform users of the river's health. This program has relied on laboratory analyses of river samples for fecal coliform bacteria levels, however, results are not available until at least 24 hours after sampling. To avoid the need for laboratory analyses, ordinary least squares (OLS) and logistic regression models were developed to predict fecal coliform bacteria concentrations and the probabilities of exceeding the Massachusetts secondary contact recreation standard for bacteria based on meteorological conditions and streamflow. The OLS models resulted in adjusted R²s ranging from 50 to 60 percent. An uncertainty analysis reveals that of the total variability of fecal coliform bacteria concentrations, 45 percent is explained by the OLS regression model, 15 percent is explained by both measurement and space sampling error, and 40 percent is explained by time sampling error. Higher accuracy in future bacteria forecasting models would likely result from reductions in laboratory measurement errors and improved sampling designs.     

ESTIMATES OF MEAN CHLOROPHYLL-a CONCENTRATION: PRECISION,ACCURACY, AND SAMPLING DESIGN1

ABSTRACT: Data from 56 north-temperate lakes and reservoirs are used to develop models predicting temporal variance as a function of the mean chlorophyll-a concentration. Trophy, as estimated by mean chlorophyll-a concentration, is shown to have little effect on the sampling effort required to achieve a pre-determined level of precision for lakes sampled year-round. Collecting ten observations results in a coefficient of variation that averaged 20 percent; collecting more than ten observations yields increasingly marginal improvements in precision. The same guidelines apply to mesotrophic or eutrophic lakes sampled in the summer, whereas oligotrophic lakes sampled in the summer require fewer observations to achieve the same level of precision. The bias resulting from collecting too few observations is minimized if five or more observations are collected.     

Evaluation of the European Space Agency Climate Change Initiative Soil Moisture Product over China Using Variance Reduction Factor

In this study, we evaluated the European Space Agency Climate Change Initiative soil moisture product v02.1 (ESA CCI SM v02.1) using in situ observations collected at 547 stations in China from 1991 to 2013. A conventional validation was first conducted, and the triple collocation errors of ESA CCI SM and the European Centre for Medium Range Weather Forecasting reanalysis data were approximately 0.053 and 0.050 m³/m³, respectively. To obtain more reliable validation results, the average soil moisture of the in situ observations per ESA CCI SM pixel was also used as the validation sites. Variance reduction factor (VRF) was adopted to quantify the accuracy of the soil moisture validation sites, and the average VRF was estimated at 4.88%. The validation results were enhanced by excluding validation sites with VRF errors greater than 5% from the statistical analysis. Although the ESA CCI SM underestimated the in situ observations with a Bias of 0.05 m³/m³, a moderately high correlation coefficient of 0.44 and a relatively small unbiased root‐mean‐square difference of 0.05 m³/m³ were observed. This study provides information on the utilization of ESA CCI SM for ecological protection, climate change, and hydrological forecasting. It also suggests the adoption of VRF for future error corrections of satellite‐based products.     

时域反射仪结合土钻法测定土壤容重

土壤水分是土地持续利用、水资源规划与管理、环境化学、节水农业技术研究的基础。介绍了时域反射仪（TDR）及其测定土壤水分的方法及应用，并利用时域反射仪结合土钻法测定土壤容重，与经典的环刀法测定土壤容重进行了比较，结果表明：两种测定结果存在一定的差异，TDR结合土钻法测定土壤容重能连续测定而且稳定性好。     

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.