OPTIMAL SPATIAL ALLOCATION OF WASTE MANAGEMENT PRACTICES TO REDUCE PHOSPHORUS POLLUTION IN A WATERSHED

Phosphorus pollution from excessive litter application and municipal discharges causes eutorphication of lakes in the Eucha-Spavinaw watershed in eastern Oklahoma and western Arkansas. Consequent algae blooms impair the taste of drinking water supply drawn from the watershed and reduce the recreational values of the lakes. The paper shows how GIS data based biophysical modeling can be used to derive spatially optimal, least-cost allocation of agricultural management practices to be combined with optimal wastewater treatment activity from the point source in order to achieve socially optimal phosphorus load in the watershed. The optimal level of phosphorus load is determined by equating marginal abatement with marginal damage cost. Transportation activities in the model allow for transportation of litter within and out of the watershed. Results show uniform regulation of litter application is excessively costly relative to measures that encourage adoption of management practices that equate marginal abatement costs across pollution sources. The results also show that change in the land use patterns in a long-run and using alum based litter additives in short-run are economically efficient management options.Environmental Economics and Policy,

The Economics of Efficient Phosphorus Abatement in a Watershed

This study presents a method to determine efficient environmental targets at a watershed level. Efficient targets are devised by estimating abatement cost and cost of environmental damages and minimizing their sum. The method was applied to a case study of phosphorus pollution in a watershed in Oklahoma. Several cumulative scenarios with alternative abatement options were simulated and efficient targets were determined. As the number of abatement options at disposal to agricultural sources increased, their optimal abatement expanded relative to the abatement at the point source. Efficient targets were found to be dependent on the choice of policy that stimulates abatement.efficiency, environmental targets, phosphorus pollution, watershed, Resource /Energy Economics and Policy,

Economic Analysis of Management Practices to Reduce Phosphorus Load to Lake Eucha and Spavinaw

Changes in management practices are often proposed to reduce phosphorus loading from a watershed due to over application of poultry litter. This study determines the choice, location, and level of each best management practice in the watershed to meet a Total Maximum Daily Load and margins of safety at least cost.best management practice, phosphorus runoff, poultry litter, Target MOTAD, Environmental Economics and Policy,

Hydrology and Sedimentology of Dynamic Rill Networks Volume II: Hydrologic Model for Dynamic Rill Networks

A comprehensive model has been developed for use in modeling the hydrologic response of rill network systems. The model, which is called HYMODRIN, is composed of both a hydrologic runoff component and a hydraulic channel routing component. The hydrologic component of the model uses a Green Ampt infiltration approach linked with a nonlinear reservoir runoff model. The channel routing component of the model is baaed on a finite element solution of the diffusion wave equations. In order to account for backwater effects the model employs a dual level iteration scheme. The model may be used in either a stand alone mode or as part of a comprehensive integrated rill erosion model. In the latter case, the hydrologic data for the rill network and the associated interrill flow areas is provided by a geographic-hydrologic interface model called GHIM. This model accepts data from a digital elevation model and translates it into a form compatible with the hydrologic model. This report contains the theoretical development and operating instructions for both GHIM and HYMODRIN. Computer listings for both programs are provided

Flow and transport experiments for a streambank seep originating from a preferential flow pathway

Streambank seeps commonly originate from localized heterogeneity or preferential flow pathways (PFPs) in riparian floodplains. However, limited field data have been reported on ground water seep flows and solute transport to seeps from PFPs. The objective of this research was to build upon previous floodplain-scale investigations of PFPs by analyzing seep discharge and transport characteristics through a single PFP. An important research question was whether this PFP could be conceptualized as a homogeneous, one-dimensional flow path. Streambank seep discharge measurements were obtained by inducing a hydraulic head in a trench injection system. Also, co-injection of Rhodamine WT (RhWT) and a potassium chloride (KCl) tracer over a 60-min period was used to investigate transport dynamics. Seep discharge and breakthrough curves for electrical conductivity (EC) and RhWT were measured at the streambank using a lateral flow collection device. The breakthrough curves were fit to one-dimensional convective-dispersion equations (CDEs) to inversely estimate solute transport parameters. The PFP from which the seep originated was clean, coarse gravel (6% by mass less than 2.0 mm) surrounded by gravel with finer particles (20% by mass less than 2.0 mm). Located approximately 2 m from the trench, the seep (50 cm by 10 cm area) required at least 40 cm of hydraulic head for flow to emerge at the streambank. At a higher hydraulic head of 125 cm, seep discharge peaked at 3.5 L/min. This research verified that localized PFPs can result in the rapid transport of water (hydraulic conductivity on the order of 400 m/d) and solutes once reaching a sufficient near-bank hydraulic head. A one-dimensional equilibrium CDE was capable of simulating the EC (R2 = 0.94) and RhWT (R2 = 0.91) breakthrough curves with minimal RhWT sorption (distribution coefficient, Kd, equal to 0.1 cm3/g). Therefore, the PFP could be conceptualized as a one-dimensional, homogenous flow and transport pathway. These results are consistent with previous research observing larger-scale phosphorus transport

PREFERENTIAL FLOW EFFECTS ON SUBSURFACE CONTAMINANT TRANSPORT IN ALLUVIAL FLOODPLAINS

For sorbing contaminants, transport from upland areas to surface water systems is typically considered to be due to surface runoff, with negligible input from subsurface transport assumed. However, certain conditions can lead to an environment where subsurface transport to streams may be significant. The Ozark region, including parts of Oklahoma, Arkansas, and Missouri, is one such environment, characterized by cherty, gravelly soils and gravel bed streams. Previous research identified a preferential flow path (PFP) at an Ozark floodplain along the Barren Fork Creek in northeastern Oklahoma and demonstrated that even a sorbing contaminant, i.e., phosphorus, can be transported in significant quantities through the subsurface. The objective of this research was to investigate the connectivity and floodplain-scale impact of subsurface physical heterogeneity (i.e., PFPs) on contaminant transport in alluvial floodplains in the Ozarks. This research also evaluated a hypothesis that alluvial groundwater acts as a transient storage zone, providing a contaminant sink during high stream flow and a contaminant source during stream baseflow. The floodplain and PFP were mapped with two electrical resistivity imaging techniques. Low-resistivity features (i.e., less than 200 Ω-m) corresponded to topographical depressions on the floodplain surface, which were hypothesized to be relict stream channels with fine sediment (i.e., sand, silt, and clay) and gravel deposits. The mapped PFP, approximately 2 m in depth and 5 to 10 m wide, was a buried gravel bar with electrical resistivity in the range of 1000 to 5000 Ω-m. To investigate the PFP, stream, and groundwater dynamics, a constant-head trench test was installed with a conservative tracer (Rhodamine WT) injected into the PFP at approximately 85 mg/L for 1.5 h. Observation wells were installed along the PFP and throughout the floodplain. Water table elevations were recorded real-time using water level loggers, and water samples were collected throughout the experiment. Results of the experiment demonstrated that stream/aquifer interaction was spatially non-uniform due to floodplain-scale heterogeneity. Transport mechanisms included preferential movement of Rhodamine WT along the PFP, infiltration of Rhodamine WT into the alluvial groundwater system, and then transport in the alluvial system as influenced by the floodplain-scale stream/aquifer dynamics. The electrical resistivity data assisted in predicting the movement of the tracer in the direction of the mapped preferential flow pathway. Spatially variable PFPs, even in the coarse gravel subsoils, affected water level gradients and the distribution of tracer into the shallow groundwater system

The hydraulic conductivity structure of gravel-dominated vadose zones within alluvial floodplains

The floodplains of many gravel-bed streams have a general stratigraphy that consists of a layer of topsoil covering gravel-dominated subsoil. Previous research has demonstrated that this stratigraphy can facilitate preferential groundwater flow through focused linear features, such as paleochannels, or gravelly regions within the vadose zone. These areas within the floodplain vadose zone may provide a route for interactions between the floodplain surface and alluvial groundwater, effectively extending the hyporheic zone across the floodplain during high stream stage. The objective of this research was to assess the structure and scale of texture heterogeneity within the vadose zone within the gravel subsoils of alluvial floodplains using resistivity data combined with hydraulic testing and sediment sampling of the vadose zone. Point-scale and broad-scale methodologies in combination can help us understand spatial heterogeneity in hydraulic conductivity without the need for a large number of invasive hydraulic tests. The evaluated sites in the Ozark region of the United States were selected due to previous investigations indicating that significant high conductivity flow zones existed in a matrix which include almost no clay content. Data indicated that resistivity corresponded with the fine content in the vadose zone and subsequently corresponds to the saturated hydraulic conductivity. Statistical analysis of resistivity data, and supported by data from the soil sampling and permeameter hydraulic testing, identified isolated high flow regions and zones that can be characterized as broad-scale high hydraulic conductivity features with potentially significant consequences for the migration of water and solutes and therefore are of biogeochemical and ecological significance

Effect of Parameter Distributions on Uncertainty Analysis of Hydrologic Models

Increasing concern about the accuracy of hydrologic and water quality models has prompted interest in procedures for evaluating the uncertainty associated with these models. If a Monte Carlo simulation is used in an uncertainty analysis, assumptions must be made relative to the probability distributions to assign to the model input parameters. Some have indicated that since these parameters can not be readily determined, uncertainty analysis is of limited value. In this article we have evaluated the impact of parameter distribution assumptions on estimates of model output uncertainty. We conclude that good estimates of the means and variances of the input parameters are of greater importance than the actual form of the distribution. This conclusion is based on an analysis using the AGNPS model

A REPEATED MEASURES ANALYSIS OF THE EFFECT OF VEGETATIVE BUFFERS ON CONTAMINANT RUNOFF FROM BERMUDAGRASS TURF

A repeated measures analysis was conducted on a set of data from a multi-year study to assess the effect of vegetative buffers on the surface runoff of selected herbicides and nutrients. Multiplicative models describing the observed behavior of runoff concentration over time for buffered and non-buffered plots were fitted on a log-transformed scale using linear mixed models with PROC MIXED in PC SAS version 6.11. A spatial power covariance structure was used. Additional models for contaminant mass flow rates were fitted to evaluate the effect of buffers on total runoff mass