SIMULATING THE IMPACT OF DRAINAGE DESIGN IN A COLD CLIMATE WITH ADAPT1

ABSTRACT: Surface and subsurface drainage make crop production economically viable in much of southern Minnesota because drainage allows timely field operations and protects field crops from extended periods of flooded soil conditions. However, subsurface drainage has been shown to increase nitrate/nitrogen losses to receiving waters. When engaging in drainage activities, farmers are increasingly being asked to consider, apart from the economic profit, the environmental impact of drainage. The Agricultural Drainage and Pesticide Transport model (ADAPT) was used in this study to evaluate the impact of subsurface drainage design on the soil water balance over a two‐year period during which observed drainage discharge data were available. Twelve modeling scenarios incorporated four drainage coefficients (DC), 0.64 cm/d, 0.95 cm/d, 1.27 cm/d, and 1.91 cm/d, and three drain depths, 0.84 m, 1.15 m, and 1.45 m. The baseline condition corresponded to the drainage system specifications at the field site: a drain depth and spacing of 1.45 m and 28 m, respectively (DC of 0.64 cm/d). The results of the two‐year simulation suggested that for a given drainage coefficient, soils with the shallower drains (but equal DC) generally have less subsurface drainage and can produce more runoff (but reduced total discharge) and evapotranspiration. The results also suggested that it may be possible to design for both water/nitrate/nitrogen reduction and crop water needs.     

MODELING THE IMPACT OF MULTICOMPONENT VOCs ON GROUND WATER USING THE STEFAN‐MAXWELL EQUATION1

ABSTRACT: Computer programs that model the fate and transport of organic contaminants through porous media typically use Fick's first law to calculate vapor phase diffusion. Fick's first law, however, is limited to the case of a single, dilute species diffusing into a stagnant, high concentration, bulk vapor phase. When dealing with more than one diffusing species and at higher concentrations, the multicomponent coupling effects on vapor phase diffusion and advection of the various constituents become significant. VLEACH, a one‐dimensional finite difference model developed for the U.S. Environmental Protection Agency (USEPA), is typical of the models using Fick's first law to model vapor‐phase diffusion. The VLEACH model was modified to accommodate up to 10 components and to calculate the binary diffusion coefficients for each of the components based on molecular weight, molecular volume, temperature and pressure, and to address the coupling effects on multiple component vapor phase diffusion and its impact on ground water. The resulting model was renamed MC‐CHEMSOIL. At low vapor phase concentrations, MC‐CHEMSOIL predicts identical ground water impacts (dissolved phase loading) to those from VLEACH 2.2a. At higher vapor phase concentrations, however, the relative difference between the models exceeded 20 percent.     

NITRATE CONCENTRATIONS IN ILLINOIS RIVERS 1967 TO 1974: COMPARISON AMONG REPORTING AGENCIES1

ABSTRACT: Long term data on surface water quality can sometimes be assembled by combining data collected by different agencies at different times and assuming that between agency differences in data quality are insignificant. The objective of this paper was to assess the quality of riverine nitrate (NO₃) concentrations in Illinois measured and reported by four agencies from 1967 to 1974 by comparing median values for similar sampling locations and periods. A total of 17 river reaches were identified for which two agencies reported NO₃ concentrations during similar periods. Nonparametric comparison of median values and analysis of covariance with discharge as a covariant produced similar results. Nitrate concentrations reported by the U.S. Geological Survey (USGS) from 1967 to 1971 were not statistically (P > 0.05) different from values reported by the Illinois State Water Survey (ISWS) for two of three river reaches. Additionally, NO₃ concentrations reported by USGS from 1972 to 1974 were not statistically different than concentrations reported by the Illinois Environmental Protection Agency (IEPA) for four of five river reaches. From 1969 to 1971, NO₃ concentrations reported by the Illinois Department of Public Heath and the Illinois Environmental Protection Agency (IDPH/IEPA) were less than one‐fourth the magnitude of values reported by ISWS. The median NO₃ concentrations measured by the Central Illinois Public Service (CIPS) were significantly greater than those measured by USGS and IDPH/IEPA in the three comparable sampling locations. The use of NO₃ concentrations measured by CIPS and IDPH/IEPA prior to 1972 is not recommended.     

MODELING SEDIMENT AND PHOSPHORUS LOSSES IN AN AGRICULTURAL WATERSHED TO MEET TMDLs1

ABSTRACT: This paper studies the effectiveness of alternative farm management strategies at improving water quality to meet Total Maximum Daily Loads (TMDLs) in agricultural watersheds. A spatial process model was calibrated using monthly flow, sediment, and phosphorus (P) losses (1994 to 1996) from Sand Creek watershed in south‐central Minnesota. Statistical evaluation of predicted and observed data gave r² coefficients of 0.75, 0.69, and 0.49 for flow (average 4.1 m³/s), sediment load (average 0.44 ton/ha), and phosphorus load (average 0.97 kg/ha), respectively. The calibrated model was used to evaluate the effects of conservation tillage, conversion of crop land to pasture, and changes in phosphorus fertilizer application rate on pollutant loads. TMDLs were developed for sediment and P losses based on existing water quality standards and guidelines. Observed annual sediment and P losses exceeded these TMDLs by 59 percent and 83 percent, respectively. A combination of increased conservation tillage, reduced application rates of phosphorus fertilizer, and conversion of crop land to pasture could reduce sediment and phosphorus loads by 23 percent and 20 percent of existing loads, respectively. These reductions are much less than needed to meet TMDLs, suggesting that control of sediment using buffer strips and control of point sources of phosphorus are needed for the remaining reductions.     

TEMPORAL CHANGES IN THE YADKIN RIVER FLOW VERSUS SUSPENDED SEDIMENT CONCENTRATION RELATIONSHIP1

ABSTRACT: Dynamic linear models (DLM) and seasonal trend decomposition (STL) using local regression, or LOESS, were used to analyze the 50‐year time series of suspended sediment concentrations for the Yadkin River, measured at the U.S. Geological Survey station at Yadkin College, North Carolina. A DLM with constant trend, seasonality, and a log₁₀ streamflow regressor provided the best model to predict monthly mean log₁₀ suspended sediment concentrations, based on the forecast log likelihood. Using DLM, there was evidence (odds approximately 69:1) that the log₁₀ streamflow versus log₁₀ suspended sediment concentration relationship has changed, with an approximate 20 percent increase in the log₁₀ streamflow coefficient over the period 1981 to 1996. However, sediment concentrations in the Yadkin River have decreased during the decade of the 1990s, which has been accompanied by a concomitant increase in streamflow variability. Although STL has been shown to be a versatile trend analysis technique, DLM is shown to be more suitable for discovery and inference of structural changes (trends) in the model coefficient describing the relationship between flow and sediment concentration.     

CHEMISTRY OF COALBED METHANE DISCHARGE WATER INTERACTING WITH SEMI‐ARID EPHEMERAL STREAM CHANNELS1

ABSTRACT: The objective of this study was to examine the chemistry of Coalbed Methane (CBM) discharge water reacting with semi‐arid ephemeral stream channels in the Powder River Basin, Wyoming. The study area consisted of two ephemeral streams: Burger Draw and Sue Draw. These streams are tributaries to the perennial Powder River. Samples were collected bimonthly from three CBM discharge points and seven channel locations in Burger Draw and Sue Draw. Samples were also collected bimonthly from the Powder River above and below the confluence of Burger Draw. Before sample collection, the pH and electrical conductivity (EC) were measured in the field. Samples were transported to the laboratory and analyzed for alkalinity, major cations, and anions. From the measurement of sodium (Na), calcium (Ca), and magnesium (Mg), practical sodium adsorption ratio (SARp) and true sodium adsorption ratio (SARt) were calculated. Results suggest pH and EC of CBM discharge water was 7.1 and 4.3 dS/m, respectively. The CBM discharge water consisted of higher concentrations of sodium and alkalinity compared to other components. The pH of CBM discharge water increased significantly (p = 0.000) in the downstream channel of Burger Draw from 7.1 to 8.84 before it joined with the Powder River. Dissolved calcium concentration of CBM discharge water decreased significantly (p = 0.000) in the downstream channel water. Subsequently, SARp increased approximately from 24 to 29. The SARt also increased significantly (p = 0.001) in the downstream channel water. For instance, SARt of CBM discharge water increased from 32.93 to 45.5 downstream channels after the confluence of Sue Draw with the Burger Draw. The only significant difference in water chemistry above and below the confluence of Burger Draw with the Powder River was pH, which increased from 8.36 to 8.52. The significant increase in SAR values of CBM discharge water in Burger Draw and Sue Draw tributaries suggest a careful monitoring of salinity and sodicity is needed if CBM discharge water is used for irrigation in semi‐arid environments. Results discussed in this study will be useful to downstream water users who depend on water for irrigation.     

THE INFLUENCE OF CLIMATE VARIATION ON THE ESTIMATION OF LOW FLOWS USED TO PROTECT WATER QUALITY: A NATIONWIDE ASSESSMENT1

ABSTRACT: Historical flow records are used to estimate the regulatory low flows that serve a key function in setting discharge permit limits through the National Pollutant Discharge Elimination System, which provides a nationwide mechanism for protecting water quality. Use of historical records creates an implicit connection between water quality protection and climate variability. The longer the record, the more likely the low flow estimate will be based on a broad set of climate conditions, and thus provides adequate water quality protection in the future. Unfortunately, a long record often is not available at a specific location. This analysis examines the connection between climate variability and the variability of biologically based and hydrologically based low flow estimates at 176 sites from the Hydro‐Climatic Data Network, a collection of stream gages identified by the USGS as relatively free of anthropogenic influences. Results show that a record of 10 to 20 years is necessary for satisfactory estimates of regulatory low flows. Although it is possible to estimate a biologically based low flow from a record of less than 10 years, these estimates are highly uncertain and incorporate a bias that undermines water quality protection.     

TRANSIENT RESPONSE FUNCTIONS FOR CONJUNCTIVE WATER MANAGEMENT IN THE SNAKE RIVER PLAIN,IDAHO1

ABSTRACT: Increasing demands on western water are causing a mounting need for the conjunctive management of surface water and ground water resources. Under western water law, the senior water rights holder has priority over the junior water rights holder in times of water shortage. Water managers have been reluctant to conjunctively manage surface water and ground water resources because of the difficulty of quantification of the impacts to surface water resources from ground water stresses. Impacts from ground water use can take years to propagate through an aquifer system. Prediction of the degree of impact to surface water resources over time and the spatial distribution of impacts is very difficult. Response functions mathematically describe the relationship between a unit ground water stress applied at a specific location and stream depletion or aquifer water level change elsewhere in the system. Response functions can be used to help quantify the spatial and temporal impacts to surface water resources caused by ground water pumping. This paper describes the theory of response functions and presents an application of transient response functions in the Snake River Plain, Idaho. Transient response functions can be used to facilitate the conjunctive management of surface and ground water not only in the eastern Snake River Plain basin, but also in similar basins throughout the western United States.     

ENHANCING FECAL COLIFORM TOTAL MAXIMUM DAILY LOAD MODELS THROUGH BACTERIAL SOURCE TRACKING1

ABSTRACT: Surface water impairment by fecal coliform bacteria is a water quality issue of national scope and importance. In Virginia, more than 400 stream and river segments are on the Commonwealth's 2002 303(d) list because of fecal coliform impairment. Total maximum daily loads (TMDLs) will be developed for most of these listed streams and rivers. Information regarding the major fecal coliform sources that impair surface water quality would enhance the development of effective watershed models and improve TMDLs. Bacterial source tracking (BST) is a recently developed technology for identifying the sources of fecal coliform bacteria and it may be helpful in generating improved TMDLs. Bacterial source tracking was performed, watershed models were developed, and TMDLs were prepared for three streams (Accotink Creek, Christians Creek, and Blacks Run) on Virginia's 303(d) list of impaired waters. Quality assurance of the BST work suggests that these data adequately describe the bacteria sources that are impairing these streams. Initial comparison of simulated bacterial sources with the observed BST data indicated that the fecal coliform sources were represented inaccurately in the initial model simulation. Revised model simulations (based on BST data) appeared to provide a better representation of the sources of fecal coliform bacteria in these three streams. The coupled approach of incorporating BST data into the fecal coliform transport model appears to reduce model uncertainty and should result in an improved TMDL.     

MODELING THE IMPACTS OF WATER LEVEL CHANGES ON A GREAT LAKES COMMUNITY1

ABSTRACT: Recent research that couples climate change scenarios based on general circulation models (GCM) with Great Lakes hydrologic models has indicated that average water levels are projected to decline in the future. This paper outlines a methodology to assess the potential impact of declining water levels on Great Lakes waterfront communities, using the Lake Huron shoreline at Goderich, Ontario, as an example. The methodology utilizes a geographic information system (GIS) to combine topographic and bathymetric datasets. A digital elevation surface is used to model projected shoreline change for 2050 using water level scenarios. An arbitrary scenario, based on a 1 m decline from February 2001 lake levels, is also modeled. By creating a series of shoreline scenarios, a range of impact and cost scenarios are generated for the Goderich Harbor and adjacent marinas. Additional harbor and marina dredging could cost as much as CDN $7.6 million. Lake freighters may experience a 30 percent loss in vessel capacity. The methodology is used to provide initial estimates of the potential impacts of climate change that can be readily updated as more robust climate change scenarios become available and is adaptable for use in other Great Lakes coastal communities.