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.     

Suspended sediment, nitrogen and phosphorus concentrations and exports during storm-events to the Tuross estuary, Australia

This paper presents a process for estimating pollutant loads from water quality data, to improve catchment-scale modelling in the region for resource management purposes. It describes a program to estimate suspended sediment, total and dissolved nitrogen and phosphorus loads to the Tuross estuary from the Tuross River catchment (1810 km(2)) of coastal southeast Australia. Event-based water quality sampling results obtained during storm events in 2005 are presented. Event 1, during July 2005 was the largest storm event in terms of peak flow for 3.5 years. Other events monitored were also in July, November and December 2005. The early July 2005 event had a flow-weighted mean suspended sediment (SS) concentration during the first 4 days of 63 mg L(-1). Of the events monitored, this was unusual as it was preceded by drought and had the largest SS concentrations (peaking at 180 mg L(-1)) during the rising-stage. In contrast, the November event had a much lower flow-weighted SS mean (28 mg L(-1)), even though peak flow magnitudes were similar. The July and November 2005 events had peak flows of 12,360 and 11,330 ML day(-1). Low-cost rising-stage siphon samplers were used to collect samples during the rapidly rising phase of these events. The use of such samplers and consideration of time-lead/lag flow adjustments, quantified using cross-correlation analysis to account for hysteresis effects, were incorporated into the load estimation techniques. The technique is a potentially useful approach for understanding relationships between water quality concentrations and flow for modelling catchment source strengths and transport processes.     

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.     

MODERNIZATION OF NATIONAL WEATHER SERVICE RIVER FORECASTING TECHNIQUES1

The National Weather Service is nearing the conclusion of a five year period of transition from index type catchment modelling to the use of conceptual hydrologic models. The decision to make this technological change was based on an extensive research project in which various catchment models were tested in a wide variety of basins and their strong and weak points ascertained. This project is described. Some of the problems involved in the changeover, which are discussed, are practical parameter optimization methods, computer requirements for the more complex technology, data requirements, fitting of the catchment model to major river systems, training of personnel and staffing problems.     

Distributed Hydrologic Modeling in the Suburban Landscape: Assessing Parameter Transferability from Gauged Reference Catchments1

Mittman, Tamara, Lawrence E. Band, Taehee Hwang, and Monica Lipscomb Smith, 2012. Distributed Hydrologic Modeling in the Suburban Landscape: Assessing Parameter Transferability from Gauged Reference Catchments. Journal of the American Water Resources Association (JAWRA) 48(3): 546-557. DOI: 10.1111/j.1752-1688.2011.00636.x Abstract: Distributed, process-based models of catchment hydrologic response are potentially useful tools for the assessment of Low Impact Development (LID) techniques in urbanized catchments. Their application is often limited, however, by the lack of continuous streamflow records to calibrate poorly constrained parameters. This article examines the transferability of soil and groundwater parameters from a forested reference catchment to a nearby suburban catchment. We use the Regional Hydro-Ecologic Simulation System (RHESSys) to develop hydrologic models of one gauged forested and one ungauged suburban catchment within the Baltimore Ecosystem Study (BES) study area. We use a parameter uncertainty framework to calibrate soil and groundwater parameters for the forested catchment, and discrete measurements of streamflow from the suburban catchment to assess parameter transferability. Results indicate that the transfer of soil and groundwater parameters from forested reference to nearby suburban catchments is viable, with performance measures for the suburban catchment often exceeding those for the forested catchment. We propose that the simplification of hydrologic processes in urbanized catchments may account for the increase in model performance in the suburban catchment.     

Calibration Catchment Selection for Flood Regionalization Modeling1

Wan Jaafar, Wan Zurina, and Dawei Han, 2012. Calibration Catchment Selection for Flood Regionalization Modeling. Journal of the American Water Resources Association (JAWRA) 48(4): 698‐706. DOI: 10.1111/j.1752‐1688.2012.00648.x Abstract: There are two unsolved problems in flood regionalization model development related to the quantity and quality of calibration catchments: (1) how many calibration catchments should be used? and (2) how to select the calibration catchments? This study explores these two questions through a case study on the median annual maximum flood (QMED) model in the United Kingdom. It has been found that the chance of developing a good QMED model decreases significantly when the number of calibration catchments drops below a critical number (e.g., 60 in the case study). However, no significant improvement is achieved if the number of calibration catchments is above it. This number could be used as a benchmark for choosing randomly selected calibration catchments. Across a broad range of calibration catchment numbers, there are good and poor calibrated models regardless of calibration catchment numbers. High quality models could be developed from a small number of calibration catchments and also poor models from a large number of calibration catchments. This indicates that the number of calibration catchments may not be the dominating factor for developing a high quality regionalization model. Instead, the information content could be more important. The study has demonstrated that the standard deviation values between the best and poorest groups are distinctive and could be used in choosing appropriate calibration catchments.     

Spatial Variability of Nitrate Concentrations Under Diverse Conditions in Tributaries to a Lake Watershed1

Abstract: Nitrate‐nitrogen (NO₃‐N) concentrations in stream water often respond uniquely to changes in inter‐annual conditions (e.g., biological N uptake and precipitation) in individual catchments. In this paper, we assess (1) how the spatial distribution of NO₃‐N concentrations varies across a dense network of nonnested catchments and (2) how relationships between multiple landscape factors [within whole catchments and hydrologically sensitive areas (HSAs) of the catchments] and stream NO₃‐N are expressed under a variety of annual conditions. Stream NO₃‐N data were collected during two synoptic sampling events across >55 tributaries and two synoptic sampling periods with >11 tributaries during summer low flow periods. Sample tributaries drain mixed land cover watersheds ranging in size from 0.150 to 312 km² and outlet directly to Cayuga Lake, New York. Changes in NO₃‐N concentration ratios between each sampling event suggest a high degree of spatial heterogeneity in catchment response across the Cayuga Lake Watershed, ranging from 0.230 to 61.4. Variations in NO₃‐N concentrations within each of the large synoptic sampling events were also high, ranging from 0.040 to 8.7 mg NO₃‐N/l (March) and 0.090 to 15.5 mg NO₃‐N/l (October). Although Pearson correlation coefficients suggest that this variability is related to multiple landscape factors during all four sampling events, partial correlations suggest percentage of row crops in the catchments as the only similar factor in March and October and catchment area as the only factor during summer low flows. Further, the strength of the relationships is typically lower in the HSAs of catchment. Advancing current understanding of such variations and relationships to landscape factors across multiple catchments – and under a variety of biogeochemical and hydrological conditions – is important, as (1) nitrate continues to be employed as an indicator of regional aquatic ecosystem health and services and (2) a unified framework approach for understanding individual catchment processes is a rapidly evolving focus for catchment‐based science and management.     

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.     

Approaches to the implementation of the Water Framework Directive: targeting mitigation measures at critical source areas of diffuse phosphorus in Irish catchments

The Water Framework Directive (WFD) has initiated a shift towards a targeted approach to implementation through its focus on river basin districts as management units and the natural ecological characteristics of waterbodies. Due to its role in eutrophication, phosphorus (P) has received considerable attention, resulting in a significant body of research, which now forms the evidence base for the programme of measures (POMs) adopted in WFD River Basin Management Plans (RBMP). Targeting POMs at critical sources areas (CSAs) of P could significantly improve environmental efficiency and cost effectiveness of proposed mitigation strategies. This paper summarises the progress made towards targeting mitigation measures at CSAs in Irish catchments. A review of current research highlights that knowledge related to P export at field scale is relatively comprehensive however; the availability of site-specific data and tools limits widespread identification of CSA at this scale. Increasing complexity of hydrological processes at larger scales limits accurate identification of CSA at catchment scale. Implementation of a tiered approach, using catchment scale tools in conjunction with field-by-field surveys could decrease uncertainty and provide a more practical and cost effective method of delineating CSA in a range of catchments. Despite scientific and practical uncertainties, development of a tiered CSA-based approach to assist in the development of supplementary measures would provide a means of developing catchment-specific and cost-effective programmes of measures for diffuse P. The paper presents a conceptual framework for such an approach, which would have particular relevance for the development of supplementary measures in High Status Waterbodies (HSW). The cost and resources necessary for implementation are justified based on HSWs' value as undisturbed reference condition ecosystems.     

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.     

Investigation of Biogeochemical Functional Proxies in Headwater Streams Across a Range of Channel and Catchment Alterations

Historically, headwater streams received limited protection and were subjected to extensive alteration from logging, farming, mining, and development activities. Despite these alterations, headwater streams provide essential ecological functions. This study examines proxy measures of biogeochemical function across a range of catchment alterations by tracking nutrient cycling (i.e., inputs, processing, and stream loading) with leaf litter fall, leaf litter decomposition, and water quality parameters. Nutrient input and processing remained highest in second growth forests (the least altered areas within the region), while recently altered locations transported higher loads of nutrients, sediments, and conductivity. Biogeochemical functional proxies of C and N input and processing significantly, positively correlated with rapid assessment results (Pearson coefficient = 0.67–0.81; P = 0.002–0.016). Additionally, stream loading equations demonstrate that N and P transport, sediment, and specific conductivity negatively correlated with rapid assessment scores (Pearson coefficient = 0.56–0.81; P = 0.002–0.048). The observed increase in stream loading with lower rapid assessment scores indicates that catchment alterations impact stream chemistry and that rapid assessments provide useful proxy measures of function in headwater ecosystems. Significant differences in nutrient processing, stream loading, water quality, and rapid assessment results were also observed between recently altered (e.g., mined) headwater streams and older forested catchments (Mann–Whitney U = 24; P = 0.01–0.024). Findings demonstrate that biogeochemical function is reduced in altered catchments, and rapid assessment scores respond to a combination of alteration type and recovery time. An analysis examining time and economic requirements of proxy measurements highlights the benefits of rapid assessment methods in evaluating biogeochemical functions.     

Dissolved organic nitrogen regulation in freshwaters

Dissolved organic nitrogen (DON) has been hypothesized to play a major role in N cycling in a variety of ecosystems. Our aim was to assess the seasonal and concentration relationships between dissolved organic carbon (DOC), DON, and NO3- within 102 streams and 16 lakes within catchments of differing complexity situated in Wales. Further, we aimed to assess whether patterns of land use, soil type, and vegetation gave consistent trends in DON and dissolved inorganic nitrogen (DIN) relationships over a diverse range of catchments. Our results reinforce that DON constitutes a significant component of the total dissolved N pool typically representing 40 to 50% of the total N in streams and lakes but sometimes representing greater than 85% of the total dissolved N. Generally, the levels of DON were inversely correlated with the concentration of DIN. In contrast to DIN concentrations, which showed distinct seasonality, DON showed no consistent seasonal trend. We hypothesize that this reflects differences in the bioavailability of these two N types. The amount of DON, DOC, and DIN was significantly related to soil type with higher DON export from Histosol-dominated catchments in comparison with Spodosol-dominated watersheds. Vegetation cover also had a significant effect on DON concentrations independent of soil type with a nearly twofold decrease in DON export from forested catchments in comparison with nonforested watersheds. Due to the diversity in catchment DON behavior, we speculate that this will limit the adoption of DON as a broad-scale indicator of catchment condition for use in monitoring and assessment programs.     

Non-point source mercury emission from the Idrija Hg-mine region: GIS mercury emission model

A mercury emission model was developed to estimate non-point source mercury (Hg) emissions occurring over the year from the Idrijca River catchment, draining the area of the world's second largest Hg mine in Idrija, Slovenia. Site-specific empirical correlations between the measured Hg emission fluxes and the parameters controlling the emission (comprising substrate Hg content, soil temperature, solar radiation and soil moisture) were incorporated into the mercury emission model developed using Geographic Information System technology. In this way, the spatial distribution and significance of the most polluted sites that need to be properly managed was assessed. The modelling results revealed that annually approximately 51 kg of mercury are emitted from contaminated surfaces in the catchment (640 km(2)), highlighting that emission from contaminated surfaces contributes significantly to the elevated Hg concentrations in the ambient air of the region. Very variable meteorological conditions in the modelling domain throughout the year resulted in the high seasonal and spatial variations of mercury emission fluxes observed. Moreover, it was found that mercury emission fluxes from surfaces in the Idrija region are 3-4 fold higher than the values commonly used in models representing emissions from global mercuriferous belts. Sensitivity and model uncertainty analysis indicated the importance of knowing not only the amount but also the type of mercury species and their binding in soils in future model development.     

Roles of land use resolution and unit-area load rates in assessment of diffuse nutrient emissions

In contrast to its counterparts in Europe and North America, the Australian National Pollutant Inventory (NPI) includes estimates of aggregated emissions of nutrients (total nitrogen and total phosphorus) from catchments and facilities. Sparse or inadequate data limit the extent to which nutrient exports may be estimated from direct observations. The paucity of data for calibration and simulation limits the use of sophisticated models in most Australian catchments. Therefore, a simple unit-area load model-Catchment Management Support System (CMSS)-was selected to estimate aggregated catchment emissions for the NPI. Estimates from models like CMSS are sensitive to spatial and categorical resolution of land uses identified within the catchment and to nutrient generation rates selected for each land use category. Using three Hawkesbury-Nepean subcatchments, we show that while high spatial resolution of land use mapping is useful, only four or five major land use categories with carefully selected generation rates were required to estimate potential nutrient exports sufficiently well and to determine subcatchments contributing most. Nutrient emission estimates proved to be highly dependent on selection of generation rates so a bootstrap technique was adopted to reduce subjectivity and to improve estimates of confidence limits. This led to a specification of new generation rates for Natural, Unimproved pasture, Rural and Urban land uses and to establishment of uncertainty limits.     

MODELING ACID MINE DRAINAGE ON A WATERSHED SCALE FOR TMDL CALCULATIONS1

ABSTRACT: The Cheat River of West Virginia is impaired by acid mine drainage (AMD). Fifty‐five of its river segments were placed on the 303(d) list, which required calculations of total maximum daily load (TMDL) to meet the water quality criteria for pH, total iron, aluminum, manganese, and zinc. An existing watershed model was enhanced to simulate AMD as nonpoint source load. The model divided a watershed into a network of catchments and river segments. Each catchment was divided into soil layers, which could contain pyrite, calcite and other minerals. A kinetic expression was used to simulate pyrite oxidation as a function of oxygen in the soil voids. Oxygen in the soil voids was consumed by pyrite oxidation and replenished by earth breathing. The by‐products of pyrite oxidation were calculated according to its mass action equations. Chemical equilibrium was used to account for the speciation of ferrous and ferric irons and precipitation of metal hydroxides. Simulated hydrology and water quality were compared to available data. The USEPA used the calibrated model to calculate the TMDLs in the Cheat River Watershed.     

Comparison of atrazine losses in three small headwater catchments

Understanding the processes causing herbicide transport to surface waters is crucial to determine mitigation options to reduce these losses. To this end, we investigated the atrazine (2-chloro-4-ethylamino-6-isopropylamino-1,3,5-triazine) transport in three agricultural catchments (1.1-2.1 km2) in the watershed of Lake "Greifensee" (Switzerland). In 1999, atrazine application data were recorded for all three catchments. Time proportional samples were taken at a high temporal resolution at the catchment outlets. Extremely wet conditions caused large relative losses from the catchments, ranging between 0.6 and 3.5% of the amount applied. Most of the atrazine load was due to event-driven diffuse losses from the fields. Farmyard runoff contributed less but caused the highest concentrations (up to 31 microg L(-1)) in the brooks. The maximum concentrations due to diffuse losses varied between 1.2 and 8.2 microg L(-1) among the catchments. Despite different absolute concentration levels, the concentration time-series were very similar. It seems that the travel-times within the catchments were mainly controlled by the rainfall pattern with little influence of the catchment properties. These properties, however, caused the relative losses to vary by a factor of 6 between the catchments. This variability could be partly explained by differences in the connectivity of the fields to the brooks and by their hydrological soil properties. A comparison of the losses from the three catchments with those from the entire watershed of Lake Greifensee demonstrated that they were representative for the larger area. Hence, the study results provide a good data set to evaluate distributed models predicting herbicide losses.     

Nitrogen removal in valley bottom wetlands: assessment in headwater catchments distributed throughout a large basin

Although the reduction of nutrient loading between uplands and streams is sometimes considered evidence of the effect of wetlands acting as buffer zones, the influence of valley bottom wetlands (VBWs) on NO(3)(-) loading has seldom been assessed at the catchment scale. The objective of this study was to quantify the impact of VBWs on NO(3)(-) concentrations in streams in the Brittany region of France. We analyzed the spatial variation in NO(3)-N concentrations in 18 headwater catchments located in a 400-km(2) basin, with varying topographic, climatic, and agricultural intensity conditions. Approximately every 10 d, water was sampled during the high flow season. We investigated the relationships between the mean NO(3)(-) concentration and different characteristics of the catchments: (i) the amount of effective rainfall, i.e., the combined effect of precipitation and actual evapotranspiration on discharge and chemical dilution, (ii) the intensity of farming, i.e., the area used for farming in the catchments and the surplus of the agricultural N budget, and (iii) the relative area of VBWs. Although the first two characteristics were the main factors controlling N concentration variability, a step-by-step regression allowed us to attribute a significant part of the NO(3)(-) concentration decrease to the increase of VBW area in each catchment. For an increase of VBW area from 11 to 16%, the NO(3)-N concentration decreased from 5.3 to 4.2 mg L(-1). Therefore in this basin, VBWs reduced the NO(3)(-) concentrations in streams with sources in agricultural fields by 30%. This work demonstrates the contribution of natural VBWs to NO(3)(-) removal at the catchment scale compared to other sources of variation, which is a current need for integrating water quality criteria into wetland management.     

SELECTED CHALLENGES IN RUNOFF GENERATION RESEARCH IN FORESTS FROM THE HILLSLOPE TO HEADWATER DRAINAGE BASIN SCALE1

ABSTRACT: When faced with practical forest land management issues such as the impacts of logging or forest conversion to other land uses, planners ideally require a comprehensive understanding of within drainage basin hydrological processes to determine the most vulnerable areas to increased storm runoff and erosion. Land managers in particular need to know the source areas and magnitude of inputs to the storm hydrograph, in terms of water quantity, sediment and solute transport; and the routing of such hydrographs from headwater to larger drainage basins. The latter includes an overall assessment at various scales of the impacts of forest disturbance and conversion on the water balance. This paper will focus on runoff generation in terms of identifying the various pathways and source areas. Such aspects will be linked with the need for a more comprehensive effort towards the field testing of so-called ‘physically based’ models of runoff generation. Some of the controversial issues arising from the difficulties in reconciling results from hydrochemical investigations with complementary hydrometric studies will be highlighted. Subsequently, attention will be given to topographic-wetness models, which have promising applications in forestland management. In addition, alternative simple models for application at the catchment scale will be assessed. The latter is in recognition that at smaller scales, heterogeneity both in time and space of soil hydraulic properties demand a greater number of parameters in modelling. Such considerations can even prove an obstacle in terms of the confident application of ‘physically based’ models.     

EFFECT OF CATCHMENT SIZE ON PHOSPHORUS EXPORT1

ABSTRACT: We tested the common assumption, made when expressing phosphorus export on an areal basis, that this export is a linear function of catchment area and found it wanting. The data show that in agricultural catchments, TP (total phosphorus) export varies as the 0.77 power of drainage basin area, resulting in a reduction in phosphorus delivery per unit area with increasing catchment size. Following further division of catchments according to agricultural practice, we found that this spatial scale effect is restricted to row crops and pastures. We present simple statistical models to allow a comparison of TP export from catchments of different size. Such models are not needed for nonrow crops, mixed agricultural and forested catchments, where TP export is a linear function of catchment size.