Estimating annual generation rates of total P and total N for different land uses in Tasmania, Australia

Water quality issues have become increasingly important to Australian catchment stakeholders. As extensive nutrient sampling and modelling expertise are often absent or unattainable, simple unit-area models like Catchment Management Support System (CMSS) remain an attractive option for informing water quality management decisions. The selection of nutrient generation rates for use in CMSS is often an arbitrary assignment based on limited literature sources or expert opinion. Using a Bayesian model to estimate nutrient generation rates for the region of Tasmania, Australia, improved the rigor of CMSS modelling and in the process highlighted that dairy pastures were the most significant contributor of total phosphorus and total nitrogen loads to Tasmanian rivers.     

Evaluating stream water quality through land use analysis in two grassland catchments: impact of wetlands on stream nitrogen concentration

We evaluated the impacts of natural wetlands and various land uses on stream nitrogen concentration in two grassland-dominated catchments in eastern Hokkaido, Japan. Analyzing land use types in drainage basins, measuring denitrification potential of its soil, and water sampling in all seasons of 2003 were performed. Results showed a highly significant positive correlation between the concentration of stream NO3-N and the proportion of upland area in drainage basins in both catchments. The regression slope, which we assumed to reflect the impact on water quality, was 24% lower for the Akkeshi catchment (0.012 +/- 0.001) than for the Shibetsu catchment (0.016 +/- 0.001). In the Akkeshi catchment, there was a significant negative correlation between the proportion of wetlands in the drainage basins and stream NO3-N concentration. Stream dissolved organic nitrogen (DON) and carbon (DOC) concentrations were significantly higher in the Akkeshi catchment. Upland and urban land uses were strongly linked to increases in in-stream N concentrations in both catchments, whereas wetlands and forests tended to mitigate water quality degradation. The denitrification potential of the soils was highest in wetlands, medium in riparian forests, and lowest in grasslands; and was significant in wetlands and riparian forests in the Akkeshi catchment. The solubility of soil organic carbon (SOC) and soil moisture tended to determine the denitrification potential. These results indicate that the water environment within the catchments, which influences denitrification potential and soil organic matter content, could have caused the difference in stream water quality between the two catchments.     

Index models to evaluate the risk of phosphorus and nitrogen loss at catchment scales

This paper investigates index models as a tool to estimate the risk of N and P source strengths and loss at the catchment scale. The index models assist managers in improving the focus of remediation actions that reduce nutrient delivery to waterbodies. N and P source risk factors (e.g. soil nutrient concentrations) and transport risk factors (e.g. distance-to-streams) are used to determine the overall risk of nutrient loss for a case study in the Tuross River catchment of coastal southeast Australia. In the development of the N index model for Tuross, particulate N was considered important based on the observed event water quality data. In contrast to previous N index models, erosion and contributing distance were therefore included in the Tuross River catchment N index. Event-based water quality monitoring, and soil information, or in data-poor catchments conceptual understanding, are essential to represent catchment-scale processes. The techniques have high applicability in other catchments, and are complementary to other modelling techniques such as process-based semi-distributed modelling. Index models generally provide much more detailed spatial resolution than fully- or semi-distributed conceptual modelling approaches. Semi-distributed models can be used to quantify nutrient loads and provide overall direction to set the broad focus of management. Index models can then be used to refine on-the-ground investigations and investment priorities. In this way semi-distributed models can be combined with index models to provide a set of powerful tools to influence management decisions and outcomes.     

DAYCENT national-scale simulations of nitrous oxide emissions from cropped soils in the United States

Until recently, Intergovernmental Panel on Climate Change (IPCC) emission factor methodology, based on simple empirical relationships, has been used to estimate carbon (C) and nitrogen (N) fluxes for regional and national inventories. However, the 2005 USEPA greenhouse gas inventory includes estimates of N2O emissions from cultivated soils derived from simulations using DAYCENT, a process-based biogeochemical model. DAYCENT simulated major U.S. crops at county-level resolution and IPCC emission factor methodology was used to estimate emissions for the approximately 14% of cropped land not simulated by DAYCENT. The methodology used to combine DAYCENT simulations and IPCC methodology to estimate direct and indirect N2O emissions is described in detail. Nitrous oxide emissions from simulations of presettlement native vegetation were subtracted from cropped soil N2O to isolate anthropogenic emissions. Meteorological data required to drive DAYCENT were acquired from DAYMET, an algorithm that uses weather station data and accounts for topography to predict daily temperature and precipitation at 1-km2 resolution. Soils data were acquired from the State Soil Geographic Database (STATSGO). Weather data and dominant soil texture class that lie closest to the geographical center of the largest cluster of cropped land in each county were used to drive DAYCENT. Land management information was implemented at the agricultural-economic region level, as defined by the Agricultural Sector Model. Maps of model-simulated county-level crop yields were compared with yields estimated by the USDA for quality control. Combining results from DAYCENT simulations of major crops and IPCC methodology for remaining cropland yielded estimates of approximately 109 and approximately 70 Tg CO2 equivalents for direct and indirect, respectively, mean annual anthropogenic N2O emissions for 1990-2003.     

Characterizing land surface erosion from cesium-137 profiles in lake and reservoir sediments

Recognition of the threat to the sustainable use of the earth's resources posed by soil erosion and associated off-site sedimentation has generated an increasing need for reliable information on global rates of soil loss. Existing methods of assessing rates of soil loss across large areas possess many limitations and there is a need to explore alternative approaches to characterizing land surface erosion at the regional and global scale. The downcore profiles of 137Cs activity available for numerous lakes and reservoirs located in different areas of the world can be used to provide information on land surface erosion within the upstream catchments. The rate of decline of 137Cs activity toward the surface of the sediment deposited in a lake or reservoir can be used to estimate the rate of surface lowering associated with eroding areas within the upstream catchment, and the concentration of 137Cs in recently deposited sediment provides a basis for estimating the relative importance of surface and channel, gully, and/or subsurface erosion as a source of the deposited sediment. The approach has been tested using 137Cs data from several lakes and reservoirs in southern England and China, spanning a wide range of specific suspended sediment yield. The results obtained are consistent with other independent evidence of erosion rates and sediment sources within the lake and reservoir catchments and confirm the validity of the overall approach. The approach appears to offer valuable potential for characterizing land surface erosion, particularly in terms of its ability to provide information on the rate of surface lowering associated with the eroding areas, rather than an average rate of lowering for the entire catchment surface.     

Factors controlling sediment and phosphorus export from two Belgian agricultural catchments

Sediment and total phosphorus (TP) export vary through space and time. This study was conducted to determine the factors controlling sediment and TP export in two agricultural catchments situated in the Belgian Loess Belt. At the outlet of these catchments runoff discharge was continuously measured and suspended sediment samples were taken during rainfall events. Within the catchments vegetation type and cover, soil surface parameters, erosion features, sediment pathways, and rainfall characteristics were monitored. Total P content and sediment characteristics such as clay, organic carbon, and suspended sediment concentration were correlated. Total sediment and TP export differ significantly between the monitored catchments. Much of the difference is due to the occurrence of an extreme event in one catchment and the morphology and spatial organization of land use in the catchments. In one catchment, the direct connection between erosive areas and the catchment outlet by means of a road system contributed to a high sediment delivery ratio (SDR) at the outlet. In the other catchment, the presence of a wide valley in the center of the catchment caused sediment deposition. Vegetation also had an effect on sediment production and deposition. Thus, many factors control sediment and TP export from small agricultural catchments; some of these factors are related to the physical catchment characteristics such as morphology and landscape structure and are (semi)permanent, while others, such as vegetation cover and land use, are time dependent.     

Assessing the Impact of Areal Precipitation Input on Streamflow Simulations Using the SWAT Model1

Masih Ilyas, Shreedhar Maskey, Stefan Uhlenbrook, and Vladimir Smakhtin, 2011. Assessing the Impact of Areal Precipitation Input on Streamflow Simulations Using the SWAT Model. Journal of the American Water Resources Association (JAWRA) 47(1):179‐195. DOI: 10.1111/j.1752‐1688.2010.00502.x Abstract: Reduction of input uncertainty is a challenge in hydrological modeling. The widely used model Soil Water Assessment Tool (SWAT) uses the data of a precipitation gauge nearest to the centroid of each subcatchment as an input for that subcatchment. This may not represent overall catchment precipitation conditions well. This paper suggests an alternative – using areal precipitation obtained through interpolation. The effectiveness of this alternative is evaluated by comparing its simulations with those based on the standard SWAT precipitation input procedure. The model is applied to mountainous semiarid catchments in the Karkheh River basin, Iran. The model performance is evaluated at daily, monthly, and annual scales by using a number of performance indicators at 15 streamflow gauging stations each draining an area in the range of 590‐42,620 km². The comparison suggests that the use of areal precipitation improves model performance particularly in small subcatchments in the range of 600‐1,600 km². The modified areal precipitation input results in increased reliability of simulated streamflows in the areas of low rain gauge density. Both precipitation input methods result in reasonably good simulations for larger catchments (over 5,000 km²). The use of areal precipitation input improves the accuracy of simulated streamflows with spatial resolution and density of rain gauges having significant impact on results.     

Uncertainty Models for Estimates of Physical Characteristics of River Segments Over Large Areas

Spatially comprehensive estimates of the physical characteristics of river segments over large areas are required in many large‐scale analyses of river systems and for the management of multiple basins. Remote sensing and modeling are often used to estimate river characteristics over large areas, but the uncertainties associated with these estimates and their dependence on the physical characteristics of the segments and their catchments are seldom quantified. Using test data with varying degrees of independence, we derived analytical models of the uncertainty associated with estimates of upstream catchment area (CA), segment slope, and mean annual discharge for all river segments of a digital representation of the hydrographic network of France. Although there were strong relationships between our test data and estimates at the scale of France, there were also large relative local uncertainties, which varied with the physical characteristics of the segments and their catchments. Discharge and CA were relatively uncertain where discharge was low and catchments were small. Discharge uncertainty also increased in catchments with large rainfall events and low minimum temperature. The uncertainty of segment slope was strongly related to segment length. Our uncertainty models were consistent across large regions of France, suggesting some degree of generality. Their analytical formulation should facilitate their use in large‐scale ecological studies and simulation models.     

Influence of Catchment-Scale Military Land Use on Stream Physical and Organic Matter Variables in Small Southeastern Plains Catchments (USA)

We conducted a 3-year study designed to examine the relationship between disturbance from military land use and stream physical and organic matter variables within 12 small (<5.5 km2) Southeastern Plains catchments at the Fort Benning Military Installation, Georgia, USA. Primary land-use categories were based on percentages of bare ground and road cover and nonforested land (grasslands, sparse vegetation, shrublands, fields) in catchments and natural catchments features, including soils (% sandy soils) and catchment size (area). We quantified stream flashiness (determined by slope of recession limbs of storm hydrographs), streambed instability (measured by relative changes in bed height over time), organic matter storage [coarse wood debris (CWD) relative abundance, benthic particulate organic matter (BPOM)] and stream-water dissolved organic carbon concentration (DOC). Stream flashiness was positively correlated with average storm magnitude and percent of the catchment with sandy soil, whereas streambed instability was related to percent of the catchment containing nonforested (disturbed) land. The proportions of in-stream CWD and sediment BPOM, and stream-water DOC were negatively related to the percent of bare ground and road cover in catchments. Collectively, our results suggest that the amount of catchment disturbance causing denuded vegetation and exposed, mobile soil is (1) a key terrestrial influence on stream geomorphology and hydrology and (2) a greater determinant of in-stream organic matter conditions than is natural geomorphic or topographic variation (catchment size, soil type) in these systems.     

Integrated watershed/river catchment planning and management: a comparison of selected Canadian and United Kingdom experiences

Whilst river catchments are axiomatic units to hydrological science, their use as administrative boundaries, even within the water function, is highly dependent upon socio‐political considerations. In moves to accomplish sustainable development, with its emphasis on proactive community procedures, the catchment/basin unit is, however, receiving detailed attention—notably in England/Wales and in Canada. This paper catalogues the patterns and processes of catchment planning on both sides of the Atlantic via case studies in Ontario and southern England (Thames). Conclusions are that scale, tradition (of consultation) and information about land‐use effects on river behaviour are important controls on the rates and styles of implementation.     

Hydrologic Sensitivity to Climate and Land Use Changes in the Santiam River Basin,Oregon

Future changes in water supply are likely to vary across catchments due to a river basin's sensitivity to climate and land use changes. In the Santiam River Basin (SRB), Oregon, we examined the role elevation, intensity of water demands, and apparent intensity of groundwater interactions, as characteristics that influence sensitivity to climate and land use changes, on the future availability of water resources. In the context of water scarcity, we compared the relative impacts of changes in water supply resulting from climate and land use changes to the impacts of spatially distributed but steady water demand. Results highlight how seasonal runoff responses to climate and land use changes vary across subbasins with differences in hydrogeology, land use, and elevation. Across the entire SRB, water demand exerts the strongest influence on basin sensitivity to water scarcity, regardless of hydrogeology, with the highest demand located in the lower reaches dominated by agricultural and urban land uses. Results also indicate that our catchment with mixed rain‐snow hydrology and with mixed surface‐groundwater may be more sensitive to climate and land use changes, relative to the catchment with snowmelt‐dominated runoff and substantial groundwater interactions. Results highlight the importance of evaluating basin sensitivity to change in planning for planning water resources storage and allocation across basins in variable hydrogeologic settings.     

EFFECTS OF SCALE ON LAND USE AND WATER QUALITY RELATIONSHIPS: A LONGITUDINAL BASIN‐WIDE PERSPECTIVE1

ABSTRACT: Human land use is a major source of change in catchments in developing areas. To better anticipate the long‐term effects of growth, land use planning requires estimates of how changes in land use will affect ecosystem processes and patterns across multiple scales of space and time. The complexity of biogeochemical and hydrologic interactions within a basin makes it difficult to scale up from process‐based studies of individual reaches to watershed scales over multiple decades. Empirical models relating land use/land cover (LULC) to water quality can be useful in long‐term planning, but require an understanding of the effects of scale on apparent land use‐water quality relationships. We empirically determined how apparent relationships between water quality and LULC data change at different scales, using LIJLC data from the Willapa Bay watershed (Washington) and water quality data collected along the Willapa and North Rivers. Spatial scales examined ranged from the local riparian scale to total upstream catchment. The strength of the correlations between LTJLC data and longitudinal water quality trends varied with scale. Different water quality parameters also varied in their response to changes in scale. Intermediate scales of land use data generally were better predictors than local riparian or total catchment scales. Additional data from the stream network did not increase the strength of relationships significantly. Because of the likelihood of scale‐induced artifacts, studies quantifying land use‐water quality relationships performed at single scales should be viewed with great caution.     

Controls on catchment-scale patterns of phosphorus in soil, streambed sediment, and stream water

Many models of phosphorus (P) transfer at the catchment scale rely on input from generic databases including, amongst others, soil and land use maps. Spatially detailed geochemical data sets have the potential to improve the accuracy of the input parameters of catchment-scale nutrient transfer models. Furthermore, they enable the assessment of the utility of available, generic spatial data sets for the modeling and prediction of soil nutrient status and nutrient transfer at the catchment scale. This study aims to quantify the unique and joint contribution of soil and sediment properties, land cover, and point-source emissions to the spatial variation of P concentrations in soil, streambed sediments, and stream water at the scale of a medium-sized catchment. Soil parent material and soil chemical properties were identified as major factors controlling the catchment-scale spatial variation in soil total P and Olsen P concentrations. Soil type and land cover as derived from the generic spatial database explain 33.7% of the variation in soil total P concentrations and 17.4% of the variation in Olsen P concentrations. Streambed P concentrations are principally related to the major element concentrations in streambed sediment and P delivery from the hillslopes due to sediment erosion. During base flow conditions, the total phosphorus (<0.45 microm) concentrations in stream water are mainly controlled by the concentrations of P and the major elements in the streambed sediment.     

Prioritisation of farm scale remediation efforts for reducing losses of nutrients and faecal indicator organisms to waterways: a case study of New Zealand dairy farming

The international competitiveness of the New Zealand (NZ) dairy industry is built on low cost clover-based systems and a favourable temperate climate that enables cows to graze pastures mostly all year round. Whilst this grazed pasture farming system is very efficient at producing milk, it has also been identified as a significant source of nutrients (N and P) and faecal bacteria which have contributed to water quality degradation in some rivers and lakes. In response to these concerns, a tool-box of mitigation measures that farmers can apply on farm to reduce environmental emissions has been developed. Here we report the potential reduction in nutrient losses and costs to farm businesses arising from the implementation of individual best management practices (BMPs) within this tool-box. Modelling analysis was carried out for a range of BMPs targeting pollutant source reduction on case-study dairy farms, located in four contrasting catchments. Due to the contrasting physical resources and management systems present in the four dairy catchments evaluated, the effectiveness and costs of BMPs varied. Farm managements that optimised soil Olsen P levels or used nitrification inhibitors were observed to result in win-win outcomes whereby nutrient losses were consistently reduced and farm profitability was increased in three of the four case study farming systems. Other BMPs generally reduced nutrient and faecal bacteria losses but at a small cost to the farm business. Our analysis indicates that there are a range of technological measures that can deliver substantial reductions in nutrient losses to waterways from dairy farms, whilst not increasing or even reducing other environmental impacts (e.g. greenhouse gas emissions and energy use). Their implementation will first require clearly defined environmental goals for the catchment/water body that is to be protected. Secondly, given that the major sources of water pollutants often differed between catchments, it is important that BMPs are matched to the physical resources and management systems of the existing farm businesses.     

Hydrochemical Buffer Assessment in Agricultural Landscapes: From Local to Catchment Scale

Non-point-source pollution of surface and groundwater is a prominent environmental issue in rural catchments, with major consequences on water supply and aquatic ecosystem quality. Among surface-water protection measures, environmental or landscape management policies support the implementation and the management of buffer zones. Although a great number of studies have focused on buffer zones, quantification of the buffer effect is still a recurring question.The purpose of this article is a critical review of the assessment of buffer-zone functioning. Our objective is to provide land planners and managers with a set of variables to assess the limits and possibilities for quantifying buffer impact at the catchment scale. We first consider the scale of the local landscape feature. The most commonly used empirical method for assessing buffers is to calculate water/nutrient budgets from inflow–outflow monitoring at the level of landscape structures. We show that several other parameters apart from mean depletion of flux can be used to describe buffer functions. Such parameters include variability, with major implication for water management. We develop a theoretical framework to clarify the assessment of the buffer effect and propose a systematic analysis taking account of temporal variability. Second, we review the current assessment of buffer effects at the catchment scale according to the theoretical framework established at the local scale. Finally, we stress the limits of direct empirical assessment at the catchment scale and, in particular, we emphasize the hierarchy in hydrological processes involved at the catchment scale: The landscape feature function is constrained by other factors (climate and geology) that are of importance at a broader spatial and temporal scale.Published on line     

Application of a Three-Dimensional Water Quality Model as a Decision Support Tool for the Management of Land-Use Changes in the Catchment of an Oligotrophic Lake

While expansion of agricultural land area and intensification of agricultural practices through irrigation and fertilizer use can bring many benefits to communities, intensifying land use also causes more contaminants, such as nutrients and pesticides, to enter rivers, lakes, and groundwater. For lakes such as Benmore in the Waitaki catchment, South Island, New Zealand, an area which is currently undergoing agricultural intensification, this could potentially lead to marked degradation of water clarity as well as effects on ecological, recreational, commercial, and tourism values. We undertook a modeling study to demonstrate science-based options for consideration of agricultural intensification in the catchment of Lake Benmore. Based on model simulations of a range of potential future nutrient loadings, it is clear that different areas within Lake Benmore may respond differently to increased nutrient loadings. A western arm (Ahuriri) could be most severely affected by land-use changes and associated increases in nutrient loadings. Lake-wide annual averages of an eutrophication indicator, the trophic level index (TLI) were derived from simulated chlorophyll a, total nitrogen, and total phosphorus concentrations. Results suggest that the lake will shift from oligotrophic (TLI = 2–3) to eutrophic (TLI = 4–5) as external loadings are increased eightfold over current baseline loads, corresponding to the potential land-use intensification in the catchment. This study provides a basis for use of model results in a decision-making process by outlining the environmental consequences of a series of land-use management options, and quantifying nutrient load limits needed to achieve defined trophic state objectives.     

AN OBJECTWE TEST FOR HYDROLOGIC SCALE1

ABSTRACT: Improving the reliability of parametric hydrologic models (sometimes called cenceptual rainfall-runoff models) in the continuous simulation of runoff from ungaged catchments has been frustrated by difficulties in estimating model parameters from catchment characteristics. An underlying problem is that these models use parameters to represent catchments as a whole, whereas data on catchment characteristics are collected at multiple field locations and are difficult to transform into one measure of collective impact. Subdividing the catchment and calibrating a stochastic parametric model to estimate distributions for the parameters that covered the range of observed streamflow values was found to improve the simulations. This paper presents an optimization of the amount of subdivision to use in simulation with a version of the Stanford Watershed Model using available climatological data. The calibration process assumes that catchment heterogeneity introduces errors that can be reduced by calibrating parameters as spatial distributions rather than single values. Calibrations for three diverse small gaged catchments located in California and in Virginia found the optimal number of subdivisions to range from 4 to 25 and the optimal scale to range from 0.3 to 2.1 mi².     

A framework for assessing the impact of land use policy on community exposure to air toxics

Our research focuses on the linkage between land use planning policy and the spatial pattern of exposure to air toxics emissions. Our objective is to develop a modeling framework for assessment of the community health risk implications of land use policy. The modeling framework is not intended to be a regulatory tool for small-scale land use decisions, but a long-range planning tool to assess the community health risk implications of alternative land use scenarios at a regional or subregional scale. This paper describes the development and application of an air toxic source model for generating aggregate emission factors for industrial and commercial zoning districts as a function of permitted uses. To address the uncertainty of estimating air toxics emission rates for planned general land use or zoning districts, the source model uses an emissions probability mass function that weights each incremental permitted land use activity by the likelihood of occurrence. We thus reduce the uncertainty involved in planning for development with no prior knowledge of the specific industries that may locate within the land use district. These air toxics emission factors can then be used to estimate pollutant atmospheric mass flux from land use zoning districts, which can then be input to air dispersion and human health risk assessment models to simulate the spatial pattern of air toxics exposure risk. The model database was constructed using the California Air Toxics Inventory, 1997 US Economic Census, and land assessment records from several California counties. The database contains information on more than 200 air toxics at the 2-digit Standard Industrial Classification (SIC) level. We present a case study to illustrate application of the model. LUAIRTOX, the interactive spreadsheet model that applies our methodology to the California data, is available at http://www2.bren.ucsb.edu/~mwillis/LUAIRTOX.htm.     

Spatial variability of soil phosphorus in relation to the topographic index and critical source areas: sampling for assessing risk to water quality

A measure of soil P status in agricultural soils is generally required for assisting with prediction of potential P loss from agricultural catchments and assessing risk for water quality. The objectives of this paper are twofold: (i) investigating the soil P status, distribution, and variability, both spatially and with soil depth, of two different first-order catchments; and (ii) determining variation in soil P concentration in relation to catchment topography (quantified as the "topographic index") and critical source areas (CSAs). The soil P measurements showed large spatial variability, not only between fields and land uses, but also within individual fields and in part was thought to be strongly influenced by areas where cattle tended to congregate and areas where manure was most commonly spread. Topographic index alone was not related to the distribution of soil P, and does not seem to provide an adequate indicator for CSAs in the study catchments. However, CSAs may be used in conjunction with soil P data for help in determining a more "effective" catchment soil P status. The difficulties in defining CSAs a priori, particularly for modeling and prediction purposes, however, suggest that other more "integrated" measures of catchment soil P status, such as baseflow P concentrations or streambed sediment P concentrations, might be more useful. Since observed soil P distribution is variable and is also difficult to relate to nationally available soil P data, any assessment of soil P status for determining risk of P loss is uncertain and problematic, given other catchment physicochemical characteristics and the sampling strategy employed.     

A Heuristic Method for Determining Changes of Source Loads to Comply with Water Quality Limits in Catchments

A common land and water management task is to determine where and by how much source loadings need to change to meet water quality limits in receiving environments. This paper addresses the problem of quantifying changes in loading when limits are specified in many locations in a large and spatially heterogeneous catchment, accounting for cumulative downstream impacts. Current approaches to this problem tend to use either scenario analysis or optimization, which suffer from difficulties of generating scenarios that meet the limits, or high complexity of optimization approaches. In contrast, we present a novel method in which simple catchment models, load limits, upstream/downstream spatial relationships and spatial allocation rules are combined to arrive at source load changes. The process iteratively establishes the critical location (river segment or lake) where the limits are most constraining, and then adjusts sources upstream of the critical location to meet the limit at that location. The method is demonstrated with application to New Zealand (268,000 km2) for nutrients and the microbial indicator E. coli, which was conducted to support policy development regarding water quality limits. The model provided useful insights, such as a source load excess (the need for source load reduction) even after mitigation measures are introduced in order to comply with E. coli limits. On the other hand, there was headroom (ability to increase source loading) for nutrients. The method enables assessment of the necessary source load reductions to achieve water quality limits over broad areas such as large catchments or whole regions.