Evaluation of Rio Grande Management Alternatives Using a Surface‐Water/Ground‐Water Model1

Abstract: Previous investigations observed significant seepage losses from the Rio Grande to the shallow aquifer between Socorro and San Antonio, New Mexico. High‐resolution telescopic modeling was used along a 10‐km reach of the Rio Grande and associated drains and canals to evaluate several management alternatives aimed at improving river conveyance efficiency. Observed data consisted of ground‐water and surface‐water elevations, seepage rates along the Rio Grande and associated canals and drains, and borehole geology. Model calibration was achieved by adjusting hydraulic conductivity and specific storage until the output matched observed data. Sensitivity analyses indicated that the system was responsive to changes in hydrogeologic properties, especially when such alterations increased vertical connectivity between layers. The calibrated model predicted that removal of the low flow conveyance channel, a major channel draining the valley, would not only decrease river seepage by 67%, but also decrease total flow through the reach by 75%. The decreased flow through the reach would result in increased water logging and an average increase in ground‐water elevations of 1.21 meter. Simulations of the system with reduced riparian evapotranspiration rates or a relocated river channel also predicted decreased river seepage, but to a much lesser degree.     

WATER LOSSES ALONG A REACH OF THE PECOS RIVER IN NEW MEXICO1

ABSTRACT: A reach of the Pecos River, located in eastern New Mexico, was examined to evaluate losses of river flows due to evaporation, seepage, and transpiration. An accurate assessment of the water losses along this reach is critical for determining how water rights are adjudicated for water users in the Pecos basin and interstate compact accounting. Water losses significantly impact flows through critical habitat for species protected under the Endangered Species Act. Daily losses of river flows were analyzed for the study reach that extends from immediately below the Pecos River confluence with Taiban Creek to the United States Geological Survey (USGS) gage near Acme. The analysis was completed with consideration for other processes including flood wave travel times and attenuation along with stream bank storage and returns. The analysis was completed using daily stream flow data from USGS gages located along the study reach. Empirical seasonal functions were developed to relate flow loss to the flow rate in the river. The functions were ultimately developed to provide a method for comparing the effects of different river flows on the available water supply.     

WATER-SURFACE PROFILES WITHOUT ENERGY LOSS COEFFICIENTS1

ABSTRACT: This paper describes a new method of computing water-surface profiles, which does not require the assumption of hydrostatic pressure or of roughness coefficients. The method is based upon distributing entropy production values along a channel as uniformly as boundary conditions permit. The method requires a discharge capacity rating at a channel cross-section within the limits of the reach of the channel. A new theoretical relationship between the kinetic and momentum coefficients for the velocity profile is used, together with a dynamic programming technique for optimal distribution of energy losses along the channel. A computer model was developed and was used to verify the methodology for flood flow and channel data at four locations.     

MODELING PERCOLATION LOSSES FROM A PONDED FIELD UNDER VARIABLE WATER-TABLE CONDITIONS1

ABSTRACT: A numerical simulation model was developed to predict the vertical and lateral percolation losses from a ponded agricultural field. The two-dimensional steady-state unsaturated/ saturated flow equation was solved using the finite-difference technique. A constant ponding depth was maintained at the soil surface with different water table conditions in an application of the model for rice fields bordered by bunds. Field experiments were conducted for two different water table depths to collect data on the spatial distribution of volumetric soil-moisture content for model verification. The measured soil-moisture content values were found to be in close agreement with those predicted by the model. The sensitivity analysis of the model with selected hydrologic conditions shows that the model is most sensitive to the values of saturated hydraulic conductivity, but relatively less sensitive to water table depth, ponding depth, and evaporation rate from the soil surface. It implies that, in a ponded rice field condition, the lateral and vertical percolation losses are mostly governed by the hydraulic conductivity of the soil. The vertical percolation losses were almost equal to the saturated hydraulic conductivity values and, in most cases, these losses increased with deeper water table depths. The lateral percolation losses also increased with deeper water table depths; however, these losses were relatively small in comparison to the vertical percolation losses. The vertical and lateral percolation losses increased with the increase in ponding depths. The lateral percolation losses through the bund decreased when the evaporation losses increased from the soil surface. The results of this study indicate that the percolation losses from a ponded field may be predicted accurately for a wide range of soil and hydrological conditions when the values of hydraulic conductivity, evaporation rate, depth of ponding, and water table depth are accurately known.     

PARAMETER ESTIMATES FOR A GROUND WATER MODEL1

ABSTRACT: This paper presents a parameter sensitivity study of a two-dimensional flow and transport model of a contaminated site. Hydrogeological and site data from previous investigations were used for calibration. The USGS contaminant transport model (MOC) was used. After flow calibration to establish a reference model, parameters were varied to examine the effect each had on predictions of a contaminant plume. Hydrogeological parameters and a step size parameter were incrementally varied individually. Each result was compared to the reference model output to evaluate changes in concentration values and contaminant plume configuration. The study indicated that a generally predictable trend can be established for some parameters not affected by pumping or similar high stresses. Ranges were identified to relate concentration error or plume change to the amount of parameter error. Some parameter perturbations produced distorted model responses at high stress locations. Porosity and anisotropy were found to be the most influential of the model parameters studied on the plume predictions. (KEY TERMS: ground water hydrology; hydrogeology; pollution modeling; water quality; model calibration.)     

Simulating nitrogen dynamics in agricultural soils fertilized with pig slurry and urea

Within the framework of an interregional project in the Emilia Romagna region of northern Italy, the coupled MACRO-SOILN model was chosen to estimate soil protective capacity against pollutants. The aim of our study was to evaluate the model to better identify key parameters and processes that influence N losses in agricultural soils. Nitrate N content was monitored in soil under corn (Zea mays L.) fertilized with urea and/or pig slurry, in two field experiments performed on four different soils: a Fienili clay, a Barco-like silt, a Sant'Omobono silt loam, and a La Boaria silty clay soil. Measurements were compared with model predictions. For all soils, nitrate content was underestimated on average by 24 to 88% at lower N rates; it was overestimated by 1 to 104% at higher N rates. The root mean square error (RMSE) was equal to 81.1%. Simulation of crop N uptake and soil water flow, estimation of the ammonia losses at pig slurry spreading, and N transformation parameter setting were considered as possible error sources. The calibration of crop N uptake gave rise to good model efficiency index values. The RMSE for the simulation of soil water content varied between 9.8 and 20.2%. A more accurate setting of the ammonia losses and of the feces transformation parameter values could allow the RMSE for the simulation of soil nitrate content to be reduced by no more than 10 to 15%. It is possible for the model not to include the simulation of processes that could have relevant effects on the soil N dynamics.     

Nitrogen loss through lateral seepage in near-trench paddy fields

A near-trench paddy field experiment with five urea application rates (0-360 kg N ha(-1) in 90-kg increments) was conducted on a paddy soil in the Taihu Lake Region of China to elucidate N losses through lateral seepage during three rice (Oryza sativa L.) growing seasons. The total N (Nt), NH4(+) -N, and NO3(-) -N concentrations in the lateral seepage water increased with increasing N rates. The seasonal Nt fluxes by lateral seepage varied from 6.8 to 25.6 kg N ha(-1) for urea application rates of 90 to 360 kg N ha(-1). Lateral seepage accounted for 4.7 to 6.6% of the Nt applied, implying that lateral seepage was an important pathway of N loss from near-trench paddy fields. The cumulative N loss via lateral seepage was significantly related to N fertilization rate (P = 0.05). Floodwater level was also identified as a main factor affecting N losses via lateral seepage from paddy fields, as indicated by a positive linear relationship (R2 = 0.43) between floodwater level and daily lateral flow during the flooded period (P = 0.05). Under the conditions of these experiments, a shallow floodwater depth of 50 mm, urea application rates of 90 kg N ha(-1) or less, and no rainfall within 1 wk after N application reduced N losses by lateral seepage from paddy fields.     

INTERFACE BETWEEN TWO DISSIMILAR FLUIDS BENEATH A SHEETPILE COFFERDAM1

ABSTRACT: A closed form solution is presented for determining the shape and location of the interface between two dissimilar fluids (having different densities) when steady flow takes place through a homogeneous and isotropic porous medium, into a sheetpile cofferdam; the interface is assumed to be sharp and the lower fluid stationary. The solution is obtained using the inverse hodograph. Numerical results are presented in nondimensional form for various parametric conditions in the physical plane; the interface pattern, as also the seepage discharge and exit gradient distribution are shown. The critical conditions of the interface are studied.     

DESIGN AND IMPACT ANALYSIS FOR DIVERSION AT COAL CREEK MINE1

A diversion system has been designed to carry the flow from East Fork of Coal Creek around the area proposed for mining at Thunder Basin Coal Company's (TBCC) Coal Creek mine in Campbell County, Wyoming. This paper describes the field and analysis procedures necessary to prepare the diversion design and impact evaluation, and the innovative concepts developed for the diversion system design to minimize impacts on downstream channel stability. Under the proposed diversion system design, water from the East Basin of Coal Creek will be diverted at two locations. At one location, flow will be impounded by a small dam and decanted by a pump through a pipeline into East Fork at the location of the second diversion. At this location, a training dike will be placed across the stream channel to divert flows into a diversion channel. Gravity flow along the diversion channel will deliver water to a playa area which will be converted into a detention basin by placing a small dam across its southern end. Flows up to the magnitude of the 24-hour 2-year peak flow will be passed directly through the detention basin into Middle Fork with negligible attenuation of flow rates. For less frequent events, water will be stored in the detention basin in order to prevent velocities in Lower Middle Fork from exceeding the maximum permissible velocity above which scouring may occur. Evaporation and seepage losses from the diversion system were estimated to be small and should be more than offset by the addition of water from the playa drainage basin into the Coal Creek drainage. Velocities predicted for the Lower Middle Fork after-the diversion is constructed are expected to be low enough that significant erosion of the channel is not expected to occur.     

MODELING THE LONG TERM IMPACTS OF USING RIGID STRUCTURES IN STREAM CHANNEL RESTORATION1

Abstract: Natural channel designs often incorporate rigid instream structures to protect channel banks, provide grade control, promote flow deflection, or otherwise improve channel stability. The long term impact of rigid structures on natural stream processes is relatively unknown. The objective of this study was to use long term alluvial channel modeling to evaluate the effect of rigid structures on channel processes and assess current and future stream channel stability. The study was conducted on Oliver Run, a small stream in Pennsylvania relocated due to highway construction. Field data were collected for one year along the 107 m reach to characterize the stream and provide model input, calibration, and verification data. FLUVIAL-12 was used to evaluate the long term impacts of rigid structures on natural channel adjustment, overall channel stability, and changing form and processes. Based on a consideration of model limitations and results, it was concluded that the presence of rigid structures reduced channel width-to-depth ratios, minimized bed elevation changes due to long term aggradation and degradation, limited lateral channel migration, and increased the mean bed material particle size throughout the reach. Results also showed how alluvial channel modeling can be used to improve the stream restoration design effort.     

Enhancing Hyrdological Simulation Program – FORTRAN Model Channel Hydraulic Representation1

Abstract: The hydrological simulation program – FORTRAN (HSPF) is a comprehensive watershed model that employs depth‐area‐volume‐flow relationships known as the hydraulic function table (FTABLE) to represent the hydraulic characteristics of stream channel cross‐sections and reservoirs. An accurate FTABLE determination for a stream cross‐section site requires an accurate determination of mean flow depth, mean flow width, roughness coefficient, longitudinal bed slope, and length of stream reach. A method that uses regional regression equations to estimate mean flow depth, mean flow width, and roughness coefficient is presented herein. FTABLES generated by the proposed method (Alternative Method) and FTABLES generated by Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) were compared. As a result, the Alternative Method was judged to be an enhancement over the BASINS method. First, the Alternative Method employs a spatially variable roughness coefficient, whereas BASINS employs an arbitrarily selected spatially uniform roughness coefficient. Second, the Alternative Method uses mean flow width and mean flow depth estimated from regional regression equations whereas BASINS uses mean flow width and depth extracted from the National Hydrography Dataset (NHD). Third, the Alternative Method offers an option to use separate roughness coefficients for the in‐channel and floodplain sections of compound channels. Fourth, the Alternative Method has higher resolution in the sense that area, volume, and flow data are calculated at smaller depth intervals than the BASINS method. To test whether the Alternative Method enhances channel hydraulic representation over the BASINS method, comparisons of observed and simulated streamflow, flow velocity, and suspended sediment were made for four test watersheds. These comparisons revealed that the method used to estimate the FTABLE has little influence on hydrologic calibration, but greatly influences hydraulic and suspended sediment calibration. The hydrologic calibration results showed that observed versus simulated daily streamflow comparisons had Nash‐Sutcliffe efficiencies ranging from 0.50 to 0.61 and monthly comparisons had efficiencies ranging from 0.61 to 0.84. Comparisons of observed and simulated suspended sediments concentrations had model efficiencies ranging from 0.48 to 0.56 for the daily, and 0.28 to 0.70 for the monthly comparisons. The overall results of the hydrological, hydraulic, and suspended sediment concentration comparisons show that the Alternative Method yielded a relatively more accurate FTABLE than the BASINS method. This study concludes that hydraulic calibration enhances suspended sediment simulation performance, but even greater improvement in suspended sediment calibration can be achieved when hydrological simulation performance is improved. Any improvements in hydrological simulation performance are subject to improvements in the temporal and spatial representation of the precipitation data.     

SEEPAGE FLOW INTO FLORIDA LAKES1

ABSTRACT: Direct measurements indicate that subsurface seepage in the littoral zone contributed 17.5 and 2.0 percent of the total hydraulic inputs, respectively, to Lakes Conway and Apopka, Florida. Two variations of seepage measuring devices were evaluated and gave relative standard deviations of 7 and 24 percent. Measurement inaccuracies were minimized by using large diameter (0.9 cm ID) plastic tubing. For a given transect perpendicular to shore, flow patterns were reproducible over time. Seepage flows ranged from 0 to 112 1/m²-day and most were between 4 and 30 1/m²-day. The detection limit was about 0.2 1/m²-day for a one-hour collection period. Seepage occured primarily within 30 m of shore and generally decreased exponentially with distance from shore. The shape of the bottom profile influenced flow patterns; lake bottoms with steeper slopes had higher flows that were compressed within a narrower zone. After a short-term rain event at Lake Conway, seepage flows increased rapidly to 2.4 times the prerain flow for 1 h and decreased to near background within about 6 h.     

A CRITICAL APPROACH TO THE CALIBRATION OF A WATERSHED MODEL1

ABSTRACT: A complex watershed-scale water quality simulation model, the Hydrological Simulation Program-FORTRAN (HSPF) model, was calibrated for a 16 km² catchment. The simulation step size was 0.33 hours with predicted and recorded hydrologic flows compared on an annual and monthly basis during a total calibration period of four years. Unguided numerical optimization when applied alone did not yield a model parameter set with acceptable predictive capability; instead, it was necessary to apply a critical process that included sensitivity analysis, numerical optimization, and testing of derived model parameter sets to evaluate their performance for periods other than those for which they were determined. Using this critical calibration process, the model was proven to have significant predictive capability. Numerical optimization is an aid for model calibration, but it must not be used blindly.     

Simulation of Surface Water for Un‐Gauged Areas with Storage‐Attenuation Wetlands1

Abstract: A nontraditional application of the Hydrological Simulation Program – FORTRAN (HSPF) model to simulate freshwater discharge to upper Charlotte Harbor along Florida’s west coast was performed. This application was different from traditional HSPF applications in three ways. First, the domain of the model was defined based on the hydraulic characteristics of the landforms using small distributed parameter discretization. Second, broad wetland land forms, representing more than 20% of this area, were simulated as reaches with storage‐attenuation characteristics and not as pervious land segments (PERLNDs). Finally, the reach flow‐tables (F‐Tables) were configured in a unique way to be calibrated representing the uncertainty of the storage‐attenuation process. Characterizing wetlands as hydrography elements allows flow from the wetlands to be treated as a stage‐dependent flux. The study was conducted for the un‐gauged portion of the Peace and Myakka rivers in west‐central Florida. Due to low gradient tidal influences, a large portion of the basin is un‐gauged. The objective of this study was to simulate stream flow discharges and to estimate freshwater inflow from these un‐gauged areas to upper Charlotte Harbor. Two local gauging stations were located within the model domain and were used for calibration. Another gauge with a shorter period of record was used for verification. A set of global hydrologic parameters were selected and tested using the parameter optimization software (PEST) during the calibration. Model results were evaluated using PEST and well‐known statistical indices. The correlation coefficients were very high (0.899 and 0.825) for the two calibration stations. Further testing of this approach appears warranted for watersheds with significant wetlands coverage.     

Further Development of a Specific Conductivity Approach to Measure Groundwater Discharge Area within Lakes

Groundwater exchanges with most lakes are rarely quantified because there are many technical challenges to quantification. We investigated a lakebed mapping approach to infer the relative areas of groundwater exchange in 12 prairie shallow lakes and five Laurentian mixed forest shallow lakes in Minnesota, USA in 2011. We used a relatively common approach (seepage meters) to provide baseline information on the magnitude and direction of flow at four locations in each lake. To expand from point measurements to the whole‐lake scale, we explored use of specific conductivity as a cheaper and more time efficient proxy for groundwater discharge to lakes. We validated the approach at near shore stations in each lake where seepage meter measurements and specific conductivity surveys overlapped. Specific conductivity surveys provided a similar assessment of groundwater discharge compared to seepage meters for 50% of the lake‐sampling period combinations. The lakebed mapping approach, when validated for a lake with a limited number of seepage meter (or alternative methods) measurements, offers the advantages of being more time and labor efficient over the use of a similar number of seepage meter monitoring locations; seepage meters (or piezometers, for example) are costlier in terms of equipment and labor, even for single‐lake studies. We show the combined approach could provide useful baselines for understanding and mapping groundwater exchange in shallow lakes.     

EPIC tile flow and nitrate loss predictions for three Minnesota cropping systems

Subsurface tile drains are a key source of nitrate N (NO3-N) losses to streams in parts of the north central USA. In this study, the Erosion Productivity Impact Calculator (EPIC) model was evaluated by comparing measured vs. predicted tile flow, tile NO3-N loss, soil profile residual NO3-N, crop N uptake, and yield, using 4 yr of data collected at a site near Lamberton, MN, for three crop rotations: continuous corn (Zea mays L.) or CC, corn-soybean [Glycine max (L.) Merr.] or CS, and continuous alfalfa (Medicago sativa L.) or CA. Initially, EPIC was run using standard Soil Conservation Service (SCS) runoff curve numbers (CN2) for CC and CS; monthly variations were accurately tracked for tile flow (r2 = 0.86 and 0.90) and NO3-N loss (r2 = 0.69 and 0.52). However, average annual CC and CS tile flows were underpredicted by -32 and -34%, and corresponding annual NO3-N losses were underpredicted by -11 and -52%. Predicted average annual tile flows and NO3-N losses generally improved following calibration of the CN2; tile flow underpredictions were -9 and - 12%, whereas NO3-N losses were 0.6 and -54%. Adjusting a N parameter further improved predicted CS NO3-N losses. Predicted monthly tile flows and NO3-N losses for the CA simulation compared poorly with observed values (r2 values of 0.27 and 0.19); the annual drainage volumes and N losses were of similar magnitude to those measured. Overall, EPIC replicated the relative impacts of the three cropping systems on N fate.     

Streamflow,Sediment Transport,and Geomorphic Change during the 2011 Flood on the Missouri River Near Bismarck–Mandan,ND

Geomorphic change from extreme events in large managed rivers has implications for river management. A steady‐state, quasi‐three‐dimensional hydrodynamic model was applied to a 29‐km reach of the Missouri River using 2011 flood data. Model results for an extreme flow (500‐year recurrence interval [RI]) and an elevated managed flow (75‐year RI) were used to assess sediment mobility through examination of the spatial distribution of boundary or bed shear stress (τ_b) and longitudinal patterns of average τ_b, velocity, and kurtosis of τ_b. Kurtosis of τ_b was used as an indicator of planform channel complexity and can be applied to other river systems. From differences in longitudinal patterns of sediment mobility for the two flows we can infer: (1) under extreme flow, the channel behaves as a single‐thread channel controlled primarily by flow, which enhances the meander pattern; (2) under elevated managed flows, the channel behaves as multithread channel controlled by the interaction of flow with bed and channel topography, resulting in a more complex channel; and (3) for both flows, the model reach lacks a consistent pattern of deposition or erosion, which indicates migration of areas of erosion and deposition within the reach. Despite caveats and limitations, the analysis provides useful information about geomorphic change under extreme flow and potential implications for river management. Although a 500‐year RI is rare, extreme hydrologic events such as this are predicted to increase in frequency.     

A SEMI-DISTRIBUTED ADAPTIVE MODEL FOR REAL-TIME FLOOD FORECASTING1

ABSTRACT: A semi-distributed deterministic model for real-time flood forecasting in large basins is proposed. Variability of rainfall and losses in space is preserved and the effective rainfall-direct runoff model segment based on the Clark procedure is incorporated. The distribution of losses in space is assumed proportional to rainfall intensity and their evolution in time is represented by the φ-index; furthermore, an initial period without production of effective rainfall is considered. The first estimation of losses and the associated forecasts of flow are performed at the time corresponding to the first rise observed in the hydrograph. Then the forecasts of flow are corrected at each subsequent time step through the updating of the φ-index. The model was tested by using rainfall-runoff events observed on two Italian basins and the predictions of flow for lead times up to six hours agree reasonably well with the observations in each event. For example, for the coefficient of persistence, which compares the model forecasts with those generated by the no-model assumption, appreciable positive values were computed. In particular, for the larger basin with an area of 4,147 km², the mean values were 0.4, 0.4 and 0.5 for forecast lead times of two hours, four hours and six hours, respectively. Good performance of the model is also shown by a comparison of its flow predictions with those derived from a unit hydrograph based model     

COMPUTATION OF STAGE-DISCHARGE RELATIONSHIPS AFFECTED BY UNSTEADY FLOW1

ABSTRACT: The dynamic relationship between stage and discharge which is unique to a particular flood for a selected station along the river can be determined via a mathematical model based on the complete one-dimensional equations of unsteady flow, i.e., the equations for the conservation of mass and momentum of the flood wave, and the Manning equation which accounts for energy losses. By assuming the bulk of the flood wave moves as a kinematic wave, the need for spatial resolution of the flood can be eliminated, and only the time variation of either the discharge or stage at the selected station is necessary for the computation of the other. The mathematical model can be used in river forecasting to convert the forecast discharge hydrograph into a stage hydrograph which properly reflects the unique dynamic stage-discharge relationship produced by the variable energy slope of the flood discharge. The model can be used also in stream gaging to convert a recorded stage hydrograph into a discharge hydrograph which properly accounts for the effects of unsteady flow. The model is applied to several observed floods at selected stations along the Lower Mississippi, Red, and Atchafalaya Rivers. The root mean square errors between observed and computed discharges are in the range of 3 to 7 percent, values well within the accuracy of the observations. A simple, easily-applied graphical procedure is also provided for estimating the magnitude of the effect of the unsteady flow on stage-discharge ratings. As a general rule, the dynamic effect may be significant if the channel bottom slope is less than 0.001 ft/ft (about 5 ft/mi) when the rate of change of stage is greater than about 0.10 ft/hr.