Modeling Effects of Natural Flow Restoration on Metals Fate and Transport in a Mountain Stream Impacted by Mine Waste1

Abstract: The effects of natural flow restoration on metals fate and transport in the Upper Tenmile Creek Watershed, Montana, were modeled using the Water Quality Analysis Simulation Program developed by the U.S. Environmental Protection Agency (USEPA). This 50‐km² watershed has over 150 historic abandoned mines, including mine waste rock and tailings, as well as adits discharging acid mine drainage, and is the primary drinking water supply for the City of Helena. Water supply diversions almost completely dewater some stream reaches during summer low flows, but the city is considering a new drinking water source and restoration of natural flows in Tenmile Creek as part of acid mine drainage remediation and broader aquatic habitat restoration. One dimensional steady‐state simulation of total recoverable cadmium, copper, lead, and zinc in the mainstem was performed, and the model was calibrated to June 2000 base‐flow data. Representative low‐flows in August and high‐flow snowmelt conditions in June were modeled using mean monthly natural flow estimates from the U.S. Geological Survey and representative USEPA metals concentrations data. The modeling showed that total recoverable metals concentrations, and especially loads, can vary significantly among input locations and over time in the watershed. Some data gaps limit evaluation of variability and increase uncertainty in several locations. Model results indicated, however, that natural low‐ and high‐flow restoration by itself can reduce some metals concentrations in the mainstem compared to June 2000 values, which were influenced by significant water diversion. Some values (such as Zn) may still exceed standards during natural August low flow due to the remaining high concentrations and loads in the primary inputs to the mainstem. Others (such as Cu) can increase during high flow due to remaining mine waste sources and loading of particulate Cu associated with erosion and transport of solids. Greater than 50% reduction in concentrations and loads from some of the main tributaries may be necessary to meet all standards, especially for potential particulate loads with higher flows in June.     

Applications of Turbidity Monitoring to Forest Management in California

Many California streams have been adversely affected by sedimentation caused by historic and current land uses, including timber harvesting. The impacts of timber harvesting and logging transportation systems on erosion and sediment delivery can be directly measured, modeled, or inferred from water quality measurements. California regulatory agencies, researchers, and land owners have adopted turbidity monitoring to determine effects of forest management practices on suspended sediment loads and water quality at watershed, project, and site scales. Watershed-scale trends in sediment discharge and responses to current forest practices may be estimated from data collected at automated sampling stations that measure turbidity, stream flow, suspended sediment concentrations, and other water quality parameters. Future results from these studies will provide a basis for assessing the effectiveness of modern forest practice regulations in protecting water quality. At the project scale, manual sampling of water column turbidity during high stream flow events within and downstream from active timber harvest plans can identify emerging sediment sources. Remedial actions can then be taken by managers to prevent or mitigate water quality impacts. At the site scale, manual turbidity sampling during storms or high stream flow events at sites located upstream and downstream from new, upgraded, or decommissioned stream crossings has proven to be a valuable way to determine whether measures taken to prevent post-construction erosion and sediment production are effective. Turbidity monitoring at the project and site scales is therefore an important tool for adaptive management. Uncertainty regarding the effects of current forest practices must be resolved through watershed-scale experiments. In the short term, this uncertainty will stimulate increased use of project and site-scale monitoring.     

Test of APEX for nine forested watersheds in East Texas

Hydrologic/water quality models are increasingly used to explore management and policy alternatives for managing water quality and quantity from intensive silvicultural practices with best management practices (BMPs) in forested watersheds due to the limited number of and cost of conducting watershed monitoring. The Agricultural Policy/Environmental eXtender (APEX) model was field-tested using 6 yr of data for flow, sediment, nutrient, and herbicide losses collected from nine small (2.58 to 2.74 ha) forested watersheds located in southwest Cherokee County in East Texas. Simulated annual average stream flow for each of the nine watersheds was within +/- 7% of the corresponding observed values; simulated annual average sediment losses were within +/- 8% of measured values for eight out of nine watersheds. Nash-Sutcliffe efficiency (EF) values ranged from 0.68 to 0.94 based on annual stream flow comparison and from 0.60 to 0.99 based on annual sediment comparison. Similar to what was observed, simulated flow, sediment, organic N, and P were significantly increased on clear-cut watersheds compared with the control watersheds. APEX reasonably simulated herbicide losses, with an EF of 0.73 and R(2) of 0.74 for imazapyr, and EF of 0.65 and R(2) of 0.68 for hexazinone based on annual values. Overall, the results show that APEX was able to predict the effects of silvicultural practices with BMPs on water quantity and quality and that the model is a useful tool for simulating a variety of responses to forest conditions.     

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.     

A FRAMEWORK FOR PHOSPHORUS TRANSPORT MODELING IN THE LAKE OKEECHOBEE WATERSHED1

ABSTRACT: A modeling framework was developed to determine phosphorus loadings to Lake Okeechobee from watersheds located north of the lake. This framework consists of the land-based model CREAMS-WT, the in-stream transport model QUAL2E, and an interface procedure to format the land-based model output for use by the in-stream model. QUAL2E hydraulics and water quality routines were modified to account for flow routing and phosphorus retention in both wetlands and stream channels. Phosphorus loadings obtained from previous applications of CREAMS-WT were used by QUAL2E, and calibration and verification showed that QUAL2E accurately simulated seasonal and annual phosphorus loadings from a watershed. Sensitivity and uncertainty analyses indicated that the accuracy of monthly loadings can be improved by using better estimates of in-stream phosphorus decay rates, ground water phosphorus concentrations, and runoff phosphorus concentrations as input to QUAL2E.     

WATERSHED AND LAKE WATER QUALITY ASSESSMENT: AN INTEGRATED MODELING APPROACH1

Watershed planning groups and action agencies seek to understand how lake water quality responds to changes in watershed management. This study developed and demonstrated the applicability of an integrated modeling approach for providing this information. An integrated model linking watershed conditions to water-quality of the receiving lake incorporated the following components: (1) an event-based AGNPS model to estimate watershed pollutant losses; (2) annualization of AGNPS results to produce annual lake pollutant loadings; (3) a base flow separation package, SAM, to estimate base flow; (4) estimates of nutrients in base flow and point sources; and (5) linkage of watershed loadings directly to EUTROMOD lake water quality algorithms. Results are presented for Melvern Lake, a 28-km² multipurpose reservoir with a 900-km² agricultural watershed in east central Kansas. Reasonable estimates of current lake quality were attained using an average phosphorus availability factor of 31 percent to calibrate model results to measured in-lake phosphorus. Comparison of a range of possible scenarios, including all cropland changed to no-till (best case) and all CRP and good-condition grasslands changed to cropland (worst case), indicated only a (4 percent change for in-lake phosphorus and a (2 percent change for chlorophyll a. These results indicated that this watershed is not sensitive to projected changes in land use and management.     

Nitrate and sediment fluxes from a California rangeland watershed

Long-term water quality records for assessing natural variability, impact of management, and that guide regulatory processes to safeguard water resources are rare for California oak woodland rangelands. This study presents a 20-yr record (1981-2000) of nitrate-nitrogen (NO(3)-N) and suspended sediment export from a typical, grazed oak woodland watershed (103 ha) in the northern Sierra Nevada foothills of California. Mean annual precipitation over the 20-yr period was 734 mm yr(-1) (range 366-1205 mm yr(-1)). Mean annual stream flow was 353 mm y(-1) (range 87-848 mm yr(-1)). Average annual stream flow was 48.1 +/- 16% of precipitation. Mean annual NO(3)-N export was 1.6 kg ha(-1) yr(-1) (range 0.18-3.6 kg ha(-1) yr(-1)). Annual NO(3)-N export significantly (P < 0.05) increased with increasing annual stream flow and precipitation. Mean daily NO(3)-N export was 0.004 kg ha(-1) d(-1) (range 10(-5) to 0.55 kg ha(-1) d(-1)). Mean annual suspended sediment export was 198 kg ha(-1) yr(-1) (range 23-479 kg ha(-1) yr(-1)). There was a positive relationship (P < 0.05) between annual suspended sediment export, annual stream flow and precipitation. Mean daily suspended sediment export was 0.54 kg ha(-1) d(-1) (range 10(-4) to 155 kg ha(-1) d(-1)). Virtually no sediment was exported during the dry season. The large variation in daily and annual fluxes highlights the necessity of using long-term records to establish quantitative water quality targets for rangelands and demonstrates the difficulty of designing a water quality monitoring program for these ecosystems.     

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.     

Direct and Indirect Effects of Forest Harvesting on Sediment Yield in Forested Watersheds of the United States

Managed forests generally produce high water quality, but degradation is possible via sedimentation if proper management is not implemented during forest harvesting. To mitigate harvesting effects on total watershed sediment yield, it is necessary to understand all processes that contribute to these effects. Forest harvesting best management practices (BMPs) focus almost exclusively on overland sediment sources, whereas in‐and‐near stream sources go unaddressed although they can contribute substantially to sediment yield. Thus, we propose a new framework to classify forest harvesting effects on stream sediment yield according to their direct and indirect processes. Direct effects are those caused by erosion and sediment delivery to surface water from overland sources (e.g., forest roads). Indirect effects are those caused by a shift in hydrologic processes due to tree removal that accounts for increases in subsurface and surface flows to the stream such that alterations in water quality are not predicated upon overland sediment delivery to the stream, but rather in‐stream processes. Although the direct/indirect distinction is often implicit in forest hydrology studies, we have formalized it as a conceptual model to help identify primary drivers of sediment yield after forest harvesting in different landscapes. Based on a literature review, we identify drivers of these effects in five regions of the United States, discuss current forest management BMPs, and identify research needs.     

Dry Weather Water Quality Loadings in Arid,Urban Watersheds of the Los Angeles Basin,California, USA1

Abstract: Dry weather runoff in arid, urban watersheds may consist entirely of treated wastewater effluent and/or urban nonpoint source runoff, which can be a source of bacteria, nutrients, and metals to receiving waters. Most studies of urban runoff focus on stormwater, and few have evaluated the relative contribution and sources of dry weather pollutant loading for a range of constituents across multiple watersheds. This study assessed dry weather loading of nutrients, metals, and bacteria in six urban watersheds in the Los Angeles region of southern California to estimate relative sources of each constituent class and the proportion of total annual load that can be attributed to dry weather discharge. In each watershed, flow and water quality were sampled from storm drain and treated wastewater inputs, as well as from in‐stream locations during at least two time periods. Data were used to calculate mean concentrations and loads for various sources. Dry weather loads were compared with modeled wet weather loads under a range of annual rainfall volumes to estimate the relative contribution of dry weather load. Mean storm drain flows were comparable between all watersheds, and in all cases, approximately 20% of the flowing storm drains accounted for 80% of the daily volume. Wastewater reclamation plants (WRP) were the main source of nutrients, storm drains accounted for almost all the bacteria, and metals sources varied by constituent. In‐stream concentrations reflected major sources, for example nutrient concentrations were highest downstream of WRP discharges, while in‐stream metals concentrations were highest downstream of the storm drains with high metals loads. Comparison of wet vs. dry weather loading indicates that dry weather loading can be a significant source of metals, ranging from less than 20% during wet years to greater than 50% during dry years.     

An Ecologically-Sustainable Surface Water Withdrawal Framework for Cropland Irrigation: A Case Study in Alabama

Agricultural production in the state of Alabama, USA, is mostly rain-fed, because of which it is vulnerable to drought during growing season. Since Alabama receives a significant portion of its annual precipitation during winter months, the goal of this study was to evaluate the feasibility of water withdrawal from streams during winter months for irrigation in the growing season. The Soil and Water Assessment Tool (SWAT) was used to estimate the quantity of water that can be sustainably withdrawn from streams during winter high flow periods. The model was successfully calibrated and validated for surface runoff, base flow, and total stream flow. The stream flows generated by the model at several locations within the watershed were then used to examine how much water can be sustainably withdrawn from streams of various orders (first, second and third). Although there was a considerable year-to-year variability in the amount of water that can be withdrawn, a 16-year average showed that first, second, and third order streams can irrigate about 11.6, 10.3, and 10.6% of their drainage areas, respectively. The percentage of drainage area that can be irrigated was not a function of stream order.     

DETERMINING NUTRIENT EXPORT COEFFICIENTS AND SOURCE LOADING UNCERTAINTY USING IN‐STREAM MONITORING DATA1

ABSTRACT: Over a three‐year period, flow and nutrients were monitored at 13 sites in the upper North Bosque River watershed in Texas. Drainage areas above sampling sites differed in percent of dairy waste application fields, forage fields, wood/range, and urban land area. A multiple regression approach was used to develop total phosphorus (TP) and total nitrogen (TN) export coefficients for the major land uses in these heterogeneous drainage areas. The largest export coefficients were associated with dairy waste application fields followed by urban, forage fields, and wood/range. An empirical model was then established to assess nutrient contribution by major sources using developed export coefficients and point source loadings from municipal wastewater treatment. This model was verified by comparison of estimated loadings to measured in‐stream data. Monte Carlo simulation techniques were applied to provide an uncertainty analysis for nutrient loads by source, based on the variance associated with each export coefficient. The largest sources of nutrients contributing to the upper North Bosque River were associated with dairy waste application fields and forage fields, while the greatest relative uncertainty in source contribution was associated with loadings from urban and wood/range land uses.     

USE OF A HYDROLOGIC MODEL IN A BASIN-WIDE WATER ALLOCATION PROCEEDING1

ABSTRACT: The Montana Department of Natural Resources and Conservation developed a hydrologic model to help analyze the effects of allocating water for consumptive and instream uses in the upper Missouri River basin of Montana. The model, a PC-based FORTRAN program, uses a mass-balance approach to compute monthly streamflows, reservoir operations, hydropower production, and irrigation and municipal water uses throughout the 54,000 square mile basin for a 59-year base period. Simulation results are presented as monthly mean and percentile-exceedence values. The model was run for baseline conditions and six hypothetical water-allocation alternatives. Results were used by staff resource area specialists to assess potential impacts to water quantity and distribution, water rights, water quality, stream channel form, fisheries, wildlife, recreation, hydropower production, and economics. These analyses were presented to the public and the decision-making board in an environmental impact statement (EIS). Though, in many instances, the model did not allow for detailed, site-specific analyses, the model was an important tool and its simulation results formed the hydrologic basis for the EIS.     

Sources and Transport of Phosphorus to Rivers in California and Adjacent States,U.S., as Determined by SPARROW Modeling

The SPARROW (SPAtially Referenced Regression on Watershed attributes) model was used to simulate annual phosphorus loads and concentrations in unmonitored stream reaches in California, U.S., and portions of Nevada and Oregon. The model was calibrated using de‐trended streamflow and phosphorus concentration data at 80 locations. The model explained 91% of the variability in loads and 51% of the variability in yields for a base year of 2002. Point sources, geological background, and cultivated land were significant sources. Variables used to explain delivery of phosphorus from land to water were precipitation and soil clay content. Aquatic loss of phosphorus was significant in streams of all sizes, with the greatest decay predicted in small‐ and intermediate‐sized streams. Geological sources, including volcanic rocks and shales, were the principal control on concentrations and loads in many regions. Some localized formations such as the Monterey shale of southern California are important sources of phosphorus and may contribute to elevated stream concentrations. Many of the larger point source facilities were located in downstream areas, near the ocean, and do not affect inland streams except for a few locations. Large areas of cultivated land result in phosphorus load increases, but do not necessarily increase the loads above those of geological background in some cases because of local hydrology, which limits the potential of phosphorus transport from land to streams.     

STREAM FLOW IN RELATION TO OHIA FOREST DECLINE ON THE ISLAND OF HAWAII1

ABSTRACT: Dieback of the ohia forest over a large part of the Hilo watershed and adjacent areas has resulted in a severe loss of the overstory crown component of the vegetation. The decline could cause serious damage to the watershed. To evaluate possible changes in conditions in the Hilo area since the decline began, stream flow, water quality, and precipitation data from 1929 to 1980 were collected and analyzed. The limited data available do not indicate that the ohia decline has resulted in a significant change in either annual stream flow or peak stream flow of a stream discharging from the ohia forest in the Hilo area. Suspended sediment production of two streams remains well within accepted levels as does the chemical makeup of components dissolved in the streams.     

Baseline requirements can hinder trades in water quality trading programs: Evidence from the Conestoga watershed

The U.S. Environmental Protection Agency (USEPA) and the U.S. Department of Agriculture (USDA) are promoting point/nonpoint trading as a way of reducing the costs of meeting water quality goals. Farms can create offsets by implementing management practices such as conservation tillage, nutrient management and buffer strips. To be eligible to sell offsets or credits, farmers must first comply with baseline requirements. USEPA guidance recommends that the baseline for nonpoint sources be management practices that are consistent with the water quality goal. A farmer would not be able to create offsets until the minimum practice standards are met. An alternative baseline is those practices being implemented at the time the trading program starts, or when the farmer enters the program. The selection of the baseline affects the efficiency and equity of the trading program. It has major implications for which farmers benefit from trading, the cost of nonpoint source offsets, and ultimately the number of offsets that nonpoint sources can sell to regulated point sources. We use a simple model of the average profit-maximizing dairy farmer operating in the Conestoga watershed in Pennsylvania to evaluate the implications of baseline requirements on the cost and quantity of offsets that can be produced for sale in a water quality trading market, and which farmers benefit most from trading.     

EFFECTIVENESS OF TIMBER HARVEST PRACTICES FOR CONTROLLING SEDIMENT RELATED WATER QUALITY IMPACTS1

ABSTRACT: Timber harvest best management practices (BMPs) in Washington State were evaluated to determine their effectiveness at achieving water quality standards pertaining to sediment related effects. A weight‐of‐evidence approach was used to determine BMP effectiveness based on assessment of erosion with sediment delivery to streams, physical disturbance of stream channels, and aquatic habitat conditions during the first two years following harvest. Stream buffers were effective at preventing chronic sediment delivery to streams and physical disturbance of stream channels. Practices for ground‐based harvest and cable yarding in the vicinity of small streams without buffers were ineffective or only partially effective at preventing water quality impacts. The primary operational factors influencing BMP effectiveness were: the proximity of ground disturbing activities to streams; presence or absence of designated stream buffers; the use of special timber falling and yarding practices intended to minimize physical disturbance of stream channels; and timing of harvest to occur during snow cover or frozen ground conditions. Important site factors included the density of small streams at harvest sites and the steepness of inner stream valley slopes. Recommendations are given for practices that provide a high confidence of achieving water quality standards by preventing chronic sediment delivery and avoiding direct stream channel disturbance.     

BIOECONOMIC ANALYSIS OF SELECTED CONSERVATION PRACTICES ON SOIL EROSION AND FRESHWATER FISHERIES1

ABSTRACT: Farmers can generate environmental benefits (improved water quality and fisheries and wildlife habitat), but they may not be able to quantify them. Furthermore, farmers may reduce their incomes from managing lands to produce these positive externalities but receive little monetary compensation in return. This study simulated the relationship between agricultural practices, water quality, fish responses to suspended sediment and farm income within two small watersheds, one of a cool water stream and one of a warm water stream. Using the Agricultural Drainage and Pesticide Transport (ADAPT) model, this study related best management practices (BMPs) to calculated instream suspended sediment concentrations by estimating sediment delivery, runoff, base flow, and streambank erosion to quantify the effects of suspended sediment exposure on fish communities. By implementing selected BMPs in each watershed, annual net farm income declined $18,000 to $28,000 (1 to 3 percent) from previous levels. “Lethal” fish events from suspended sediments in the cool water watershed decreased by 60 percent as conservation tillage and riparian buffers increased. Despite reducing suspended sediments by 25 percent, BMPs in the warm water watershed did not reduce the negative response of the fisheries. Differences in responses (physical and biological) between watersheds highlight potential gains in economic efficiency by targeting BMPs or by offering performance based “green payments.”     

EUTROPHICATION MEASURES FOR SMALL LAKE WATER QUALITY MANAGEMENT1

ABSTRACT: Continuous measures of water quality responsiveness to changes in nutrient loadings are developed for use by environmental regulators attempting to achieve the greatest degree of water quality improvement consistent with competing uses of scarce resources. Recent contributions to the literature, based upon the nutrient loading concept, provide statistical relationships applicable to broad categories of lakes. In order to meet the requirements of water quality regulators, equations based upon USEPA survey data are developed in terms of the functional forms formulated and tested in the literature. Since recent authors have observed that much existing survey data are inappropriate for the formulation and testing of nutrient loading relationships, we have utilized results from the literature that work best with the existing data base, on grounds that such data are most inexpensively available to water quality managers. The nature of existing data limitations is discussed, and approaches to the revision of current results in light of improvements in the data base are indicated.     

Particulate phosphorus transport within stream flow of an agricultural catchment

There is interest in quantifying phosphorus (P) loss from intensively grazed dairy landscapes to identify key pathways and target remediation methods. The Bog Burn drains a dairying catchment in Southland, New Zealand, and has been monitored at fortnightly intervals over a 12-mo period at four sites for suspended sediment (SS), dissolved reactive phosphorus (DRP), and total phosphorus (TP). Time-integrated samplers, deployed at 0.6 median water depth at each site (calculated from previous year's flow data), collected sediment samples, which were analyzed for SS, bioavailable phosphorus (BAP), and TP. Mean concentrations of DRP and TP in stream flow and BAP and TP in sediment were generally highest in summer or autumn (0.043 mg DRP L(-1), 0.160 mg TP L(-1), 173 mg BAP kg(-1), 2228 mg TP kg(-1)) and lowest in winter or spring (0.012 mg DRP L(-1), 0.034 mg TP L(-1), 6 mg BAP kg(-1), 711 mg TP kg(-1)), while loads were highest in winter. Analysis of (137)Cs concentrations in trapped sediment, topsoil, subsoil, and stream bed and bank sediment indicated that trapped sediment was derived from topsoil and entered the stream either through tile drainage or, to a lesser extent, overland flow. Because concentrations of DRP and TP in stream flow are in excess of recommended limits for good water quality (>0.01 mg DRP L(-1), 0.033 mg TP L(-1)), management should focus on the topsoil and specifically on decreasing P loss via tile drainage. This is best achieved by decreasing soil Olsen P concentrations, especially because, on average, Olsen P concentrations in the catchment were above the agronomic optimum.