Spatiotemporal Variability of Modeled Watershed Scale Surface‐Depression Storage and Runoff for the Conterminous United States

This study explores the viability of using simulated monthly runoff as a proxy for landscape‐scale surface‐depression storage processes simulated by the United States Geological Survey’s National Hydrologic Model (NHM) infrastructure across the conterminous United States (CONUS). Two different temporal resolution model codes (daily and monthly) were run in the NHM with the same spatial discretization. Simulated values of daily surface‐depression storage (treated as a decimal fraction of maximum volume) as computed by the daily Precipitation‐Runoff Modeling System (NHM‐PRMS) and normalized runoff (0 to 1) as computed by the Monthly Water Balance Model (NHM‐MWBM) were aggregated to monthly and annual values for each hydrologic response unit (HRU) in the CONUS geospatial fabric (HRU; n = 109,951) and analyzed using Spearman’s rank correlation test. Correlations between simulated runoff and surface‐depression storage aggregated to monthly and annual values were compared to identify where which time scale had relatively higher correlation values across the CONUS. Results show Spearman’s rank values >0.75 (highly correlated) for the monthly time scale in 28,279 HRUs (53.35%) compared to the annual time scale in 41,655 HRUs (78.58%). The geographic distribution of HRUs with highly correlated monthly values show areas where surface‐depression storage features are known to be common (e.g., Prairie Pothole Region, Florida).     

DEEP WELL DISPOSAL OF SECONDARY EFFLUENT IN DADE COUNTY,FLORIDA1

ABSTRACT: Data were obtained from drilling and testing of a test injection well for deep underground injection of waste water effluent from the proposed 50-million-gallon-per-day (mgd) South District Regional Wastewater Treatment Plant of the Miami-Dade Water and Sewer Authority, Dade County, Florida. The drilling operation progressed in stages, each stage coverting the strata to be sealed off by the 48-inch, 40-inch, 30-inch, and 20-inch casings, respectively. Total depth of the well is 3,200 feet. The top of the saline, cavernous, dolomitic Boulder Zone was found at 2,790 feet below the surface and is separated from the Floridan aquifer above by approximately 1,100 feet of confining limestone layers. These confining layers were determined, by packer testing, to be very effective. The transmissivity of the Boulder Zone was estimated to be 14 × 10⁶ gallons per day per foot (gpd/ft) from the data obtained from pump out tests. An 8,000-gallon-per-minute (gpm) injection test was conducted to confirm well performance under operating conditions. Based on all of the data obtained, it was concluded that underground injection into the Boulder Zone of secondary waste water effluent from the proposed treatment plant is feasible, both hydraulically and environmentally. A monitoring system was proposed to provide a record of the effects of injection on the subsurface environment.     

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.     

A stochastic approach for the operation of a wind and pumped storage plant under a deregulated environment

To improve the competiveness in the energy market, it is necessary that the wind power plants provide guaranteed power generation, although, it is not possible to forecast power availability from wind power plant accurately. This paper presents a stochastic model and solution technique for the combined operation of wind and pumped storage power plants to improve the power availability and increasing the profit considering uncertainties of wind power generation. In this model, uncertainties in wind data have been forecasted for grid connected day-ahead market using Weibull distribution model. The imbalances in the forecasted wind data and the market demand have been reduced by operating the pumped storage power plant. In this stochastic mixed integer problem, pumped storage plant can take the supply either from the grid or from the wind power plant for the pumping operation to store the energy in order to utilize this energy during peak hours for increasing the overall revenue. The reliability of the pumped storage is improved by replacing the conventional unit with the adjustable speed type pumped storage unit. In order to prove the optimality of the solution, two case studies were considered. In case study–I, scheduling is provided by operating the conventional pumped storage unit, whereas in case study–II, adjustable speed pumped storage unit has been used. It has been found that the adjustable speed pumped storage unit has further reduced the imbalance between generated power and demand. The complete approach has been formulated and implemented using AMPL software.     

VERIFICATION AND APPLICATION OF REGIONAL EQUATIONS FOR THE STORAGE FUNCTION RUNOFF MODEL1

The storage function model is a nonlinear rainfall-runoff model that has been developed for and applied to flood runoff analysis in Japan. This paper extends the model applicability by developing practical equations for estimating model parameters which are appropriate on a regional basis, i.e., so-called regional equations. Previously, the parameters were computed from historical data for a specific basin or from relationships that do not account for land use and topography. To develop the regionalized equations, model parameters were identified for 91 flood events from 22 watersheds in Japan by applying a mathematical optimization technique. Results from 39 of these events were statistically compared and regional relationships were determined as a function of land use, basin area and rainfall intensity. The utility of the estimated equations were tested by computing runoff hydrographs for lumped basins. The estimated parameters were also applied in a distributed watershed model formulation. Both applications showed acceptable results that validate the use of the regionalized relationships.     

RUNOFF SIMULATION USING RADAR RAINFALL DATA1

ABSTRACT: Rainfall data products generated with the national network of WSR-88D radars are an important new data source provided by the National Weather Service. Radar-based data include rainfall depth on an hourly basis for grid cells that are nominally 4 km square. The availability of such data enables application of improved techniques for rainfall-runoff simulation. A simple quasi-distributed approach that applies a linear runoff transform to grid-ded rainfall excess has been developed. The approach is an adaptation of the Clark conceptual runoff model, which employs translation and linear storage. Data development for, and results of, an initial application to a 4160 km² watershed in the Midwestern U.S. are illustrated.     

Parameter Uncertainty Reduction for SWAT Using Grace,Streamflow, and Groundwater Table Data for Lower Missouri River Basin1

Abstract: This study incorporates the newly available Gravity Recovery and Climate Experiment (GRACE) water storage data and water table data from well logs to reduce parameter uncertainty in Soil and Water Assessment Tool (SWAT) calibration using a SUFI2 (sequential uncertainty fitting) framework for the Lower Missouri River Basin. Model evaluations are performed in multiple stages using a multiobjective function consisting of multisite streamflow and GRACE water storage data as well as a groundwater component. Results show that (1) a model calibrated with both streamflow and GRACE data simultaneously can maintain the water balance for the whole basin, but may improperly partition surface flow and base flow. Additional inclusion of the groundwater constraint can significantly improve the model performance in groundwater hydrological processes. In our case, the estimation of specific yield of shallow aquifers has been increased to 10^?2 from previous much underestimated level (<10^?3). (2) The daily streamflow data are needed to confine the parameters related to water flow in channels such as the Manning’s coefficient, which are less sensitive to the monthly simulations. (3) Parameters are nonuniformly sensitive for different goal variables, and thus, proper specification of a prior distribution of parameters may be the key factor for global optimization algorithms to obtain stable and realistic model performance.     

TRANSIENT STORAGE ASSESSMENTS OF DYE‐TRACER INJECTIONS IN RIVERS OF THE WILLAMETTE BASIN,OREGON1

ABSTRACT: Rhodamine WT dye‐tracer injections in rivers of the Willamette Basin yield concentration‐time curves with characteristically long recession times suggestive of active transient storage processes. The scale of drainage areas contributing to the stream reaches studied in the Willamette Basin ranges from 10 to 12,000 km². A transient storage assessment of the tracer studies has been completed using the U.S. Geological Survey's One‐dimensional Transport with Inflow and Storage (OTIS) model, which incorporates storage exchange and decay functions along with the traditional dispersion and advection transport equation. The analysis estimates solute transport of the dye. It identifies first‐order decay coefficients to be on the order of 10^?5/sec for the nonconservative Rhodamine W.T. On an individual subreach basis, the first‐order decay is slower (typically by an order of magnitude) than the transient storage process, indicating that nonconservative tracers may be used to evaluate transient storage in rivers. In the transient storage analysis, a dimensionless parameter (A_s/A) expresses the spatial extent of storage zone area relative to stream cross section. In certain reaches of Willamette Basin pool‐and‐riffle, gravel‐bed rivers, this parameter was as large as 0.5. A measure of the storage exchange flux was calculated for each stream subreach in the simulation analysis. This storage exchange is shown subjectively to be higher at higher stream discharges. Hyporheic linkage between streams and subsurface flows is the probable physical mechanism contributing to a significant part of this inferred active transient storage. Hyporheic linkages are further suggested by detailed measurements of river discharge with an Acoustic Doppler Current Profiler system delineating zones in two large rivers where water alternately enters and leaves the surface channels through gravel‐and‐cobble riverbeds. Measurements show patterns of hyporheic exchange that are highly variable in time and space.     

Modelling of CO2 arrival time at Ketzin – Part I

The injection of CO₂ at the Ketzin storage site and the chemical detection of its arrival in the observation well allowed testing of different numerical simulation codes. ECLIPSE 100 (E100, black-oil simulator), ECLIPSE 300 (E300, compositional CO₂STORE) and MUFTE-UG were used for predictive modelling applying a constant injection rate of 1 kg s^?1 CO₂ and for a history match applying the actual variable injection rate which ranged from 0 to 0.7 kg s^?1 and averaged 0.23 kg s^?1. The geological model applied, is based on all available geophysical and geological information and has been the same for all programs.The results of the constant injection regime show a good agreement among the programs with a discrepancy of 21–33% for the CO₂ arrival times. However, it is determined from the comparison of the cumulative mass of CO₂ at the time of CO₂ arrival that the injection regime is an important factor for the accurate prediction of CO₂ migration within a saline aquifer. Comparing the actual variable injection regime with the simulations applying a constant injection rate the results are relatively inaccurate.Regarding the actual variable injection regime, which was evaluated using all three simulators, the computational results show a good agreement with the data actually measured at the first observation well. Here, the calculated arrival times exceeded the actual ones by 8.1% (E100), 9.2% (E300) and 17.7% (MUFTE-UG).It can be concluded that irrespective of the deviations of the simulations, due to combinations of different codes and slight differences in input parameters, all three programs are well equipped to give a reliable estimate of the arrival of CO₂. Deviations in the results mainly occur due to different input data and grid size choices done by the different modelling teams working independently of each other. Deviations of the simulations results compared to the actual CO₂ arrival time result from uncertainties in the implementation of the geological model, which was set up based on well log data and analogue studies.     

OBSERVATIONS ON KINEMATIC RESPONSE IN URBAN RUNOFF MODELS1

ABSTRACT: This paper first discusses the results of sensitivity analyses conducted on various parameters of the San Francisco Stormwater Model ta version of WREM) and the Penn State Runoff Model in terms of their impact on outflow hydrographs. The parameters considered within a idealized catchment include: basin shape, imperivous fraction, overland roughness and slope: deterntion depth; infiltration capacity; and hyetograph timing. Second, the results for the hypothetical catchment are extended to the lazzard laboratory surfaces (asphalt, grass, roofing material) as a mean of illustrating the need for changes in model structure, as opposed to continued parameter adjustment Finally the effect of altering the scale of hydraulic representation in the surface runoff and sewer transport calculations are demonstrated for two gaged watersheds in Hamburg, West Germany.     

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.     

SIMPLE MEASUREMENTS OF SOIL SURFACE MICRORELIEF CHARACTERISTICS1

ABSTRACT: Soil surface roughness and the topographic orientation of microrelief influence runoff and water storage. Distances measured following the surface of the microrelief and horizontally across the area with the microrelief can be used with the number of microridges or microhills to estimate the average cross-sectional areas of the depression, the microslope angles, and crest heights. Microridge sideslope angles which effect the rate of runoff were estimated reasonably well in field plots by the method.     

INTEGRATING GEOGRAPHIC INFORMATION SYSTEMS IN GROUND-WATER APPLICATIONS USING NUMERICAL MODELING TECHNIQUES1

ABSTRACT: The Sand and Gravel Aquifer is the sole source of potable water in Escambia County, Florida. In order to better understand the hydraulics of the aquifer, a numerical computer model of the aquifer was developed. The model applied a finite element technique which allowed for density-dependent transport and flow in three dimensions. The modeling technique was integrated with GIS to develop a system for optimal management of the resource. The GIS was the primary tool in the development of the model grid, as well as being the integral component in the modeling procedure. Multiple model grids were developed for simulating regional flow and local flow/transport phenomena. The model grids were generated by the GIS where nodal and element sequencing were recorded. The grid topology was stored in the GIS with the element numbers, node numbers, and the related hydrogeologic attributes. The Triangulated Irregular Network (TIN) module was used for transferring interpolated value between GIS coverages. TIN allowed a fit of the model grid to the physical dimensions of the aquifer and for interpolating boundary values for telescopically refined grids. Calculations between TIN surfaces provided the residuals of the dependent variable from observed TIN surfaces. Model calibration was conducted within the GIS environment through a combination of visual and relational querying. The GIS provided an integrated environment which facilitated model analyses and data storage and retrieval.     

Monitoring of the microbial community composition in saline aquifers during CO2 storage by fluorescence in situ hybridisation

This study reveals the first analyses of the composition and activity of the microbial community of a saline CO₂ storage aquifer. Microbial monitoring during CO₂ injection has been reported. By using fluorescence in situ hybridisation (FISH), we have shown that the microbial community was strongly influenced by the CO₂ injection. Before CO₂ arrival, up to 6 × 10⁶ cells ml⁻¹ were detected by DAPI staining at a depth of 647 m below the surface. The microbial community was dominated by the domain Bacteria that represented approximately 60% to 90% of the total cell number, with Proteobacteria and Firmicutes as the most abundant phyla comprising up to 47% and 45% of the entire population, respectively. Both the total cell counts as well as the counts of the specific physiological groups revealed quantitative and qualitative changes after CO₂ arrival. Our study revealed temporal outcompetition of sulphate-reducing bacteria by methanogenic archaea. In addition, an enhanced activity of the microbial population after five months CO₂ storage indicated that the bacterial community was able to adapt to the extreme conditions of the deep biosphere and to the extreme changes of these atypical conditions.     

Spatial Variability of Groundwater Depth and Quality Parameters in the National Capital Territory of Delhi

The groundwater quantity and quality scenario is of much concern in the National Capital Territory of Delhi, India, which necessitates an investigation to envisage the extent of spatial variability of groundwater depth and pollutant concentration levels in this region. Therefore, in this study, an effort was made to generate the spatial variability map of groundwater depth and quality parameters (viz. chloride, electrical conductivity, fluoride, magnesium, and nitrate). Ordinary kriging was used to analyze the spatial variability of groundwater depth and quality parameters, whereas indicator kriging was used to analyze groundwater quality parameters equal to or greater than the pollution threshold values. It was observed that the semivariogram parameters fitted well in the exponential model for water depth and in the spherical model for water quality parameters. The generated spatial variability maps indicated that in 43% of the study area, groundwater depth was within 20 m. The salinity level was higher than 2.5 dS m⁻¹ in 69% of the study area and the nitrate concentration exceeded 45 mg l⁻¹ in 36% of the area. The probability maps showed that about 24% of the area had the highest probability (0.8–1.0) of exceedence of the threshold electrical conductivity value and an area of 2% exhibited the highest probability of exceedence of the threshold value of nitrate concentration in the groundwater. The generated spatial variability and probability maps will assist water resource managers and policymakers in development of guidelines in judicious management of groundwater resources for agricultural and drinking purposes in the study area.     

COMPUTATION OF RUNOFF REDUCTION CAUSED BY FARM PONDS1

ABSTRACT: A two-parameter farm pond storage index, FPSI, was Used to adjust computed surface. runoff using the partial area runoff contribution resulting from runoff captured by farm ponds. The validity of the index method was tested by fitting a continuous accounting version of the Soil Conservation Service curve number procedure to surface runoff data from each of three watersheds, first with and then without the FPSI routine. Evapotranspiration computed with the Jensen-Haise method and rainfall were input to the model. A linear relationship was assumed between the storage index and the portion of the controlled drainage area that was contributing to runoff. Adjusting the computed runoff with the FPSI reduced the coefficient of variation of monthly measured versus computed surface runoff for each of the three watersheds. The correlation coefficients for the same comparisons were increased. The annual predicted surface runoff Was improved for 12 of the 17 station years of data tested. The farm pond storage index could be used with any surface runoff model to improve the prediction of runoff from watersheds with drainage areas greater than 1 square mile and with about 20 percent or more of the drainage area controlled by farm ponds.     

Thermal Characteristics of Stormwater Runoff from Asphalt and Sod Surfaces1

Abstract: Urban impervious surfaces absorb and store thermal energy, particularly during warm summer months. During a rainfall/runoff event, thermal energy is transferred from the impervious surface to the runoff, causing it to become warmer. As this higher temperature runoff enters receiving waters, it can be harmful to coldwater habitat. In an urban watershed, impervious asphalt surfaces (roads, parking lots, and driveways) and pervious residential lawns comprise a significant portion of the watershed area. A paired asphalt‐turfgrass sod plot was constructed to compare the thermal runoff characteristics between asphalt and turfgrass sod surfaces, to identify meteorological variables that influence these thermal characteristics, and to evaluate evaporative heat loss for runoff from asphalt surfaces. Rainfall simulations were conducted during the summers of 2004 and 2005 under a range of climatic conditions. Asphalt surface temperatures immediately prior to rainfall simulations averaged 43.6°C and decreased an average of 12.3°C over 60 min as rain cooled the surface. In contrast, presimulation sod surface temperatures averaged only 23.3°C and increased an average of 1.3°C throughout the rainfall events. Heat transferred from the asphalt to the runoff resulted in initial asphalt runoff temperatures averaging 35.0°C that decreased by an average of 4.1°C at the end of the event. Sod runoff temperatures averaged only 25.5°C and remained fairly constant throughout the simulations. Multivariable regression equations were developed to predict (1) average asphalt surface temperature (R² = 0.90) and average asphalt runoff temperature (R² = 0.92) as a function of solar radiation, rain temperature, and wind speed, and (2) average sod surface temperature (R² = 0.85) and average sod runoff temperature (R² = 0.94) as a function of solar radiation, rain temperature, rain intensity, and wind speed. Based on a heat balance analysis, existing evaporation equations developed from studies on lakes were not adequate to predict evaporation from runoff on a heated impervious surface. The combined heat from the asphalt and sod plots was an average of 38% less than the total heat had the total area consisted solely of asphalt.     

STREAM DISCHARGE PREDICTION VIA A GRID BASED SOIL WATER ROUTING WITH PADDY FIELDS1

ABSTRACT: A grid based daily hydrologic model for a watershed with paddy fields was developed to predict the stream discharge. ASCII formatted elevation, soil, and land use data supported by the GRASS Geographic Information System are used to generate distributed results such as surface runoff and subsurface flow, soil water content, and evapotranspiration. The model uses a single flow path algorithm and simulates a water balance at each grid element. A linear reservoir assumption was used to predict subsurface runoff components. The model was applied to a 75.6 km² watershed located in the middle of South Korea, and observed stream flow hydrographs from 1995 and 1996 were compared to model predictions. The stream flow predictions of 1995 and 1996 generally agreed with the observed flow, resulting in a Nash‐Sutcliffe efficiency R² of 0.60 and 0.62, respectively. The hydraulic conductivity for percolating water through the saturated layer affected baseflow generation. The levee height of the paddy influenced the time and magnitude of the surface runoff, depending on irrigation management. The model will be used for making low flow management decisions by evaluating the role of each land use to stream flow, especially in case of paddy decrease by gradual urbanization of a watershed.     

Modeling of Vertical Solute Dispersion in a Sediment Bed Enhanced by Wave‐Induced Interstitial Flow1

Abstract: Mass (solute) transport in a stream or lake sediment bed has a significant effect on chemical mass balances and microbial activities in the sediment. A “1D vertical dispersion model” is a useful tool to analyze or model solute transfer between river or lake water and a sediment bed. Under a motionless water column, solute transfer into and within the sediment bed is by molecular diffusion. However, surface waves or bed forms create periodic pressure waves along the sediment/water interface, which in turn induce flows in the pores of the sediment bed. The enhancement of solute transport by these interstitial periodic flows in the pores has been incorporated in a 1D depth‐dependent “enhanced dispersion coefficient (D_E).” Typically, D_E diminishes exponentially with depth in the sediment bed. Relationships have been developed to estimate D_E as a function of the characteristics of sediment (particle size, hydraulic conductivity, and porosity) and pressure waves (wave length and height). In this paper, we outline and illustrate the calculation of D_E as well as the penetration depth (d_p) of the flow effect. Sample applications to illustrate the computational procedure are provided for dissolved oxygen transfer into a stream gravel bed and release of phosphorus from a lake bed. The sensitivity of the results to input parameter values is illustrated, and compared with the errors obtained when interstitial flow is ignored. Maximum values of D_E near the sediment surface can be on the order of 1 cm²/s in a stream gravel bed with standing waves, and 0.001 cm²/s in a fine sand lake bed under progressive surface waves, much larger than molecular diffusion coefficients.