On the Role of Spatial,Temporal, and Climatic Forces on Stream Sediment Loading from Rural and Urban Ecosystems

In this study, we characterize the greatest sediment loading events by their sediment delivery behavior; dominant climate, watershed, and antecedent conditions; and their seasonal distribution for rural and urban land uses. The study area is Paradise Creek Watershed, a mixed land use watershed in northern Idaho dominated by saturation excess processes in the upstream rural area and infiltration excess in the downstream urban area. We analyzed 12 years of continuous streamflow, precipitation, and watershed data at two monitoring stations. We identified 137 sediment loading events in the upstream rural section of the watershed and 191 events in the downstream urban section. During the majority of these events conditions were transport limited and the sediment flush occurred early in the event, generally in the first 20% of elapsed event time. Statistical analysis including two dozen explanatory variables showed peak discharge, event duration, and antecedent baseflow explained most of the variation in event sediment load at both stations and for the watershed as a whole (R² = 0.73‐0.78). In the rural area, saturated soils combined with spring snowmelt in March led to the greatest loading events. The urban area load contribution peaked in January, which could be a re‐suspension of streambed sediments from the previous water year. Throughout the study period, one event contributed, on average, 33% of the annual sediment load but only accounted for 2% of the time in a year.     

Weighting Nitrogen and Phosphorus Pixel Pollutant Loads to Represent Runoff and Buffering Likelihoods

Watershed models often estimate annual nitrogen (N) or phosphorus (P) pollutant loads in rural areas with export coefficient (EC) (kg/ha/yr) values based on land cover, and in urban areas as the product of spatially uniform event mean concentration (EMC) (mg/L) values and runoff volume. Actual N and P nonpoint source (NPS) pollutant loading has more spatial complexity due to watershed variation in runoff likelihood and buffering likelihood along surface and subsurface pathways, which can be represented in a contributing area dispersal area (CADA) NPS model. This research develops a CADA NPS model to simulate how watershed properties of elevation, land cover, and soils upslope and downslope of each watershed pixel influence nutrient loading. The model uses both surface and subsurface runoff indices (RI), and surface and subsurface buffer indices (BI), to quantify the runoff and buffering likelihood for each watershed pixel, and generate maps of weighted EC and EMC values that identify NPS pollutant loading hotspots. The research illustrates how CADA NPS model maps and pixel loading values are sensitive to the spatial resolution and accuracy of elevation and land cover data, and model predictions can represent the lower and upper bounds of NPS loading. The model provides managers with a tool to rapidly visualize, rank, and investigate likely areas of high nutrient export.     

From agricultural intensification to conservation: sediment transport in the Raccoon River, Iowa, 1916-2009

Fluvial sediment is a ubiquitous pollutant that negatively affects surface water quality and municipal water supply treatment. As part of its routine water supply monitoring, the Des Moines Water Works (DMWW) has been measuring turbidity daily in the Raccoon River since 1916. For this study, we calibrated daily turbidity readings to modern total suspended solid (TSS) concentrations to develop an estimation of daily sediment concentrations in the river from 1916 to 2009. Our objectives were to evaluate long-term TSS patterns and trends, and relate these to changes in climate, land use, and agricultural practices that occurred during the 93-yr monitoring period. Results showed that while TSS concentrations and estimated sediment loads varied greatly from year to year, TSS concentrations were much greater in the early 20th century despite drier conditions and less discharge, and declined throughout the century. Against a backdrop of increasing discharge in the Raccoon River and widespread agricultural adaptations by farmers, sediment loads increased and peaked in the early 1970s, and then have slowly declined or remained steady throughout the 1980s to present. With annual sediment load concentrated during extreme events in the spring and early summer, continued sediment reductions in the Raccoon River watershed should be focused on conservation practices to reduce rainfall impacts and sediment mobilization. Overall, results from this study suggest that efforts to reduce sediment load from the watershed appear to be working.     

Modeling Streamflow and Water Quality Sensitivity to Climate Change and Urban Development in 20 U.S. Watersheds

Watershed modeling in 20 large, United States (U.S.) watersheds addresses gaps in our knowledge of streamflow, nutrient (nitrogen and phosphorus), and sediment loading sensitivity to mid‐21st Century climate change and urban/residential development scenarios. Use of a consistent methodology facilitates regional scale comparisons across the study watersheds. Simulations use the Soil and Water Assessment Tool. Climate change scenarios are from the North American Regional Climate Change Assessment Program dynamically downscaled climate model output. Urban and residential development scenarios are from U.S. Environmental Protection Agency's Integrated Climate and Land Use Scenarios project. Simulations provide a plausible set of streamflow and water quality responses to mid‐21st Century climate change across the U.S. Simulated changes show a general pattern of decreasing streamflow volume in the central Rockies and Southwest, and increases on the East Coast and Northern Plains. Changes in pollutant loads follow a similar pattern but with increased variability. Ensemble mean results suggest that by the mid‐21st Century, statistically significant changes in streamflow and total suspended solids loads (relative to baseline conditions) are possible in roughly 30‐40% of study watersheds. These proportions increase to around 60% for total phosphorus and total nitrogen loads. Projected urban/residential development, and watershed responses to development, are small at the large spatial scale of modeling in this study.     

Water Quality Signals from Rural Land Use and Exurbanization in a Mountain Landscape: What's Clear and What's Confounded?

In mountainous landscapes with high climatic and geomorphic variability, how do rural land uses and exurbanization alter hydrology and water quality? We evaluated effects of rural land use and exurbanization on streamflows, suspended sediment concentrations and loads, specific conductance, and summer water temperatures in 12 streams and rivers within the Upper Little Tennessee River basin in the southern Appalachian Mountains. Eleven streams featured low levels of development (>61% forest cover) but differed in land use patterning, basin size, annual precipitation, and watershed morphology. One urban stream, located within the largest town in the basin, provided the high development comparative endpoint. Even low levels of rural development and exurbanization were associated with substantial increases in suspended sediment concentrations, sediment loads, and summer stream temperature daily maxima and diurnal variation. Observed summer temperature increases were much larger than would be expected due to global climate change over the next century. Specific conductance was idiosyncratic among the smaller streams. These water quality changes were not accompanied by streamflow changes that were discernible amid the high natural variation in precipitation and geomorphology. The water quality findings suggest the need for applying the best management practices, including riparian buffers, to even low levels of rural development.     

NUTRIENT CONTRIBUTION OF NONPOINT SOURCE RUNOFF IN THE LAS VEGAS VALLEY1

ABSTRACT: A Geographic Information System (GIS) based non‐point source runoff model is developed for the Las Vegas Valley, Nevada, to estimate the nutrient loads during the years 2000 and 2001. The estimated nonpoint source loads are compared with current wastewater treatment facilities loads to determine the non‐point source contribution of total phosphorus (TP), total nitrogen (TN), and total suspended solids (TSS) on a monthly and annual time scale. An innovative calibration procedure is used to estimate the pollutant concentrations for different land uses based on available water quality data at the outlet. Results indicate that the pollutant concentrations are higher for the Las Vegas Valley than previous published values for semi‐arid and arid regions. The total TP and TN loads from nonpoint sources are approximately 15 percent and 4 percent, respectively, of the total load to the receiving water body, Lake Mead. The TP loads during wet periods approach the permitted loads from the wastewater treatment plants that discharge into Las Vegas Wash. In addition, the GIS model is used to track pollutant loads in the stream channels for one of the subwatersheds. This is useful for planning the location of Best Management Practices to control nonpoint pollutant loads.     

BUILDUP,WASHOFF, AND EVENT MEAN CONCENTRATIONS1

ABSTRACT: Most watershed water quality simulation models require the user to specify pollutant buildup and washoff rate parameters for pollutants, by land use. Buildup and washoff rates are difficult to measure directly, and only limited guidance and few observed data are available from the literature. Many studies, however, report storm event mean concentrations (EMCs). These EMCs must arise as a result of the buildup and washoff processes, but typically represent the net contribution from a variety of pervious and impervious surfaces. This paper explores the relationship between EMCs and buildup/washoff parameters. An assumption of the mathematical form of the buildup/washoff relationship gives an algebraic expression for the EMC consistent with model assumptions. This yields techniques to separate observed EMCs into contributions from different land uses and from pervious and impervious surfaces. Given this relationship, numerical optimization may be used to estimate site specific values of buildup and washoff parameters from observed storm EMCs for use in modeling. Use of this approach helps ensure that model parameters are consistent with observed data, providing a rational starting point for final model calibration. Several site examples demonstrate use of the method.     

WATERSHED OPTIMIZATION OF AGRICULTURAL BEST MANAGEMENT PRACTICES: CONTINUOUS SIMULATION VERSUS DESIGN STORMS1

ABSTRACT: Nonpoint source (NPS) models and expert opinions are often used to prescribe best management practices (BMPs) for controlling NPS pollution. An optimization algorithm (e.g., a genetic algorithm, or GA) linked with a NPS model (e.g., Annualized AGricultural Nonpoint Source pollution model, or AnnAGNPS), can be used to more objectively prescribe BMPs and to optimize NPS pollution control measures by maximizing pollutant reduction and net monetary return from a watershed. Pollutant loads from design storms and annual loads from a continuous simulation can both be used for optimizing BMP schemes. However, which strategy results in a better solution (in terms of providing water quality protection) for a watershed is not clear. The specific objective of the study was to determine the differences in watershed pollutant loads, in an experimental watershed in Pennsylvania, resulting from optimization analyses performed using pollutant loads from a series of five 2‐yr 24‐hr storm events, a series of five 5‐yr 24‐hr storm events, and cumulative pollutant loads from a continuous simulation of five years of weather data. For each of these three different event alternatives, 100 near optimal solutions (BMP schemes) were generated. Sediment (Sed), sediment nitrogen (SedN), dissolved N (SolN), sediment organic carbon (SedOC), and sediment phosphorus (SedP) loads from a different five‐year period (an evaluation period) suggest that the optimal BMP schemes resulting from the use of annual cumulative pollutant loads from a continuous simulation of five years of weather data provide smaller cumulative NPS pollutant loads at the watershed outlet.     

DRY AND WET WEATHER FLOW NUTRIENT LOADS FROM A LOS ANGELES WATERSHED1

Effective watershed management requires an accurate assessment of the pollutant loads from the associated point and nonpoint sources. The importance of wet weather flow (WWF) pollutant loads is well known, but in semi‐arid regions where urbanization is significant the pollutant load in dry weather flow (DWF) may also be important. This research compares the relative contributions of potential contaminants discharged in DWF and WWF from the Ballona Creek Watershed in Los Angeles, California. Models to predict DWF and WWF loads of total suspended solids, biochemical oxygen demand, nitrate‐nitrogen, nitrite‐nitrogen, ammonia‐nitrogen, total Kjeldahl nitrogen, and total phosphorus from the Ballona Creek Watershed for six water years dating from 1991 to 1996 were developed. The contaminants studied were selected based on data availability and their potential importance in the degradation of Ballona Creek and Santa Monica Bay beneficial uses. Wet weather flow was found to contribute approximately 75 percent to 90 percent of the total annual flow volume discharged by the Ballona Creek Watershed. Pollutant loads are also predominantly due to WWF, but during the dry season, DWF is a more significant contributor. Wet weather flow accounts for 67 to 98 percent of the annual load of the constituents studied. During the dry season, however, the portion attributable to DWF increases to greater than 40 percent for all constituents except biochemical oxygen demand and total suspended solids. When individual catchments within the watershed are considered, the DWF pollutant load from the largest catchment is similar to the WWF pollutant load in two other major catchments. This research indicates WWF is the most significant source of nonpoint source pollution load on an annual basis, but management of the effects of the nonpoint source pollutant load should consider the seasonal importance of DWF.     

Influence of Sampling Frequency on Estimation of Annual Total Phosphorus and Total Suspended Solids Loads1

Abstract: The determination of sediment and nutrient loads is typically based on the collection and analysis of grab samples. The frequency and regularity of traditional sampling may not provide representation of constituent loading, particularly in systems with flashy hydrology. At two sites in the Little Bear River, Utah, continuous, high‐frequency turbidity was used with surrogate relationships to generate estimates of total phosphorus and total suspended solids concentrations, which were paired with discharge to estimate annual loads. The high frequency records were randomly subsampled to represent hourly, daily, weekly, and monthly sampling frequencies and to examine the effects of timing, and resulting annual load estimates were compared to the reference loads. Higher frequency sampling resulted in load estimates that better approximated the reference loads. The degree of bias was greater at the more hydrologically responsive site in the upper watershed, which required a higher sampling frequency than the lower watershed site to achieve the same level of accuracy in estimating the reference load. The hour of day and day of week of sampling impacted load estimation, depending on site and hydrologic conditions. The effects of sampling frequency on the determination of compliance with a water quality criterion were also examined. These techniques can be helpful in determining necessary sampling frequency to meet the objectives of a water quality monitoring program.     

EFFECT OF WATERSHED SUBDIVISION ON SWAT FLOW,SEDIMENT, AND NUTRIENT PREDICTIONS1

ABSTRACT: The size, scale, and number of subwatersheds can affect a watershed modeling process and subsequent results. The objective of this study was to determine the appropriate level of subwatershed division for simulating flow, sediment, and nutrients over 30 years for four Iowa watersheds ranging in size from 2,000 to 18,000 km² with the Soil and Water Assessment Tool (SWAT) model. The results of the analysis indicated that variation in the total number of subwatersheds had very little effect on streamflow. However, the opposite result was found for sediment, nitrate, and inorganic P; the optimal threshold subwatershed sizes, relative to the total drainage area for each watershed, required to adequately predict these three indicators were found to be around 3, 2, and 5 percent, respectively. Decreasing the size of the subwatersheds below these threshold levels does not significantly affect the predicted levels of these environmental indicators. These threshold subwatershed sizes can be used to optimize input data preparation requirements for SWAT analyses of other watersheds, especially those within a similar size range. The fact that different thresholds emerged for the different indicators also indicates the need for SWAT users to assess which indicators should have the highest priority in their analyses.     

SPATIAL AND TEMPORAL VARIABILITY IN SUSPENDED SEDIMENT YIELDS OF NORTH ISLAND BASINS,NEW ZEALAND1

ABSTRACT: Specific annual suspended sediment yields and their standard deviations are presented for 47 basins of North Island, New Zealand. Most of the variance in yields is explained by catchment mean rainfall. Rivers with similar flow range have similar suspended sediment concentration ratings, independent of differing watershed lithology and regolith, except for six basins having an abundance of soft fine sediments. Prediction equations for yield and its standard deviation are derived for four essentially arbitrary regions. AU feature rainfall as the independent variable. Differences between regions may owe to variations in intensity, frequency, and duration patterns of storms and, in one area, to bed material size as well. The temporal distribution of annual yields from a basin m be modeled by a two-parameter lognormal function: the prediction equations above may be used to evaluate this function at a site for which suspended sediment data are unavailable.     

Development of a SWAT Patch for Better Estimation of Sediment Yield in Steep Sloping Watersheds1

Abstract: The watershed scale Soil and Water Assessment Tool (SWAT) model divides watersheds into smaller subwatersheds for simulation of rainfall‐runoff and sediment loading at the field level and routing through stream networks. Typically, the SWAT model first needs to be calibrated and validated for accurate estimation through adjustment of sensitive input parameters (i.e., Curve Number values, USLE P, slope and slope‐length, and so on). However, in some instances, SWAT‐simulated results are greatly affected by the watershed delineation and Digital Elevation Models (DEM) cell size. In this study, the SWAT ArcView GIS Patch II was developed for steep sloping watersheds, and its performance was evaluated for various threshold values and DEM cell size scenarios when delineating subwatersheds using the SWAT model. The SWAT ArcView GIS Patch II was developed using the ArcView GIS Avenue program and Spatial Analyst libraries. The SWAT ArcView GIS Patch II improves upon the SWAT ArcView GIS Patch I because it reflects the topographic factor in calculating the field slope‐length of Hydrologic Response Units in the SWAT model. The simulated sediment value for 321 subwatersheds (watershed delineation threshold value of 25 ha) is greater than that for 43 subwatersheds (watershed delineation threshold value of 200 ha) by 201% without applying the SWAT ArcView GIS Patch II. However, when the SWAT ArcView GIS Patch II was applied, the difference in simulated sediment yield decreases for the same scenario (i.e., difference in simulated sediment with 321 subwatersheds and 43 subwatersheds) was 12%. The simulated sediment value for DEM cell size of 50 m is greater than that for DEM cell size of 10 m by 19.8% without the SWAT ArcView GIS Patch II. However, the difference becomes smaller (3.4% difference) between 50 and 10 m with the SWAT ArcView GIS Patch II for the DEM scenarios. As shown in this study, the SWAT ArcView GIS Patch II can reduce differences in simulated sediment values for various watershed delineation and DEM cell size scenarios. Without the SWAT ArcView GIS Patch II, variations in the SWAT‐simulated results using various watershed delineation and DEM cell size scenarios could be greater than those from input parameter calibration. Thus, the results obtained in this study show that the SWAT ArcView GIS Patch II should be used when simulating hydrology and sediment yield for steep sloping watersheds (especially if average slope of the subwatershed is >25%) for more accurate simulation of hydrology and sediment using the SWAT model. The SWAT ArcView GIS Patch II is available at http://www.EnvSys.co.kr/~swat for free download.     

Coupling PCSWMM and WASP to Evaluate Green Stormwater Infrastructure Impacts to Storm Sediment Loads in an Urban Watershed

We coupled rainfall–runoff and instream water quality models to evaluate total suspended solids (TSS) in Wissahickon Creek, a mid‐sized urban stream near Philadelphia, Pennsylvania. Using stormwater runoff and instream field data, we calibrated the model at a subdaily scale and focused on storm responses. We demonstrate that treating event mean concentrations as a calibration parameter rather than a fixed input can substantially improve model performance. Urban stormwater TSS concentrations vary widely in time and space and are difficult to represent simply. Suspended and deposited sediment pose independent stressors to stream biota and model results suggest that both currently impair stream health in Wissahickon Creek. Retrofitting existing detention basins to prioritize infiltration reduced instream TSS loads by 20%, suggesting that infiltration mitigates sediment more effectively than detention. Infiltrating stormwater from 30% of the watershed reduced instream TSS loads by 47% and cut the frequency of TSS exceeding 100 mg/L by half. Settled loads and the frequency of high TSS values were reduced by a smaller fraction than suspended loads and duration at high TSS values. A widely distributed network of infiltration‐focused projects is an effective stormwater management strategy to mitigate sediment stress. Coupling rainfall–runoff and water quality models is an important way to integrate watershed‐wide impacts and evaluate how management directly affects urban stream health.     

Integrated Modular Modeling of Water and Nutrients From Point and Nonpoint Sources in the Patuxent River Watershed1

Abstract: We present a simple modular landscape simulation model that is based on a watershed modeling framework in which different sets of processes occurring in a watershed can be simulated separately with different models. The model consists of three loosely coupled submodels: a rainfall‐runoff model (TOPMODEL) for runoff generation in a subwatershed, a nutrient model for estimation of nutrients from nonpoint sources in a subwatershed, and a stream network model for integration of point and nonpoint sources in the routing process. The model performance was evaluated using monitoring data in the watershed of the Patuxent River, a tributary to the Chesapeake Bay in Maryland, from July 1997 through August 1999. Despite its simplicity, the landscape model predictions of streamflow, and sediment and nutrient loads were as good as or better than those of the Hydrological Simulation Program‐Fortran model, one of the most widely used comprehensive watershed models. The landscape model was applied to predict discharges of water, sediment, silicate, organic carbon, nitrate, ammonium, organic nitrogen, total nitrogen, organic phosphorus, phosphate, and total phosphorus from the Patuxent watershed to its estuary. The predicted annual water discharge to the estuary was very close to the measured annual total in terms of percent errors for both years of the study period (≤2%). The model predictions for loads of nutrients were also good (20‐30%) or very good (<20%) with exceptions of sediment (40%), phosphate (36%), and organic carbon (53%) for Year 1.     

WATER QUALITY IN AGRICULTURAL,URBAN, AND MIXED LAND USE WATERSHEDS1

ABSTRACT: Water quality and nonpoint source (NPS) pollution are important issues in many areas of the world, including the Inner Bluegrass Region of Kentucky where urban development is changing formerly rural watersheds into urban and mixed use watersheds. In watersheds where land use is mixed, the relative contributions of NPS pollution from rural and urban land uses can be difficult to separate. To better understand NPS pollution sources in mixed use watersheds, surface water samples were taken at three sites that varied in land use to examine the effect of land use on water quality. Within the group of three watersheds, one was predominately agriculture (Agricultural), one was predominately urban (Urban), and a third had relatively equal representation of both types of land uses (Mixed). Nitrogen (N), phosphorus (P), total suspended solids (TSS), turbidity, pH, temperature, and streamflow were measured for one year. Comparisons are made among watersheds for concentration and fluxes of water quality parameters. Nitrate and orthophosphate concentrations were found to be significantly higher in the Agricultural watershed. Total suspended solids, turbidity, temperature, and pH, were found to be generally higher in the Urban and Mixed watersheds. No differences were found for streamflow (per unit area), total phosphorus, and ammonium concentrations among watersheds. Fluxes of orthophosphate were greater in the Agricultural watershed that in the Urban watershed while fluxes of TSS were greater in the Mixed watershed when compared to the Agricultural watershed. Fluxes of nitrate, ammonium, and total phosphorus did not vary among watersheds. It is apparent from the data that Agricultural land uses are generally a greater source of nutrients than the Urban land uses while Urban land uses are generally a greater source of suspended sediment.     

流域污染物通量测算方法研究

流域水系内污染物通量不仅能够用于评价各类污染源的水污染物入河负荷,也是对流域污染特征,水污染物在河流水体中复杂迁移、转化过程的最直观反应。准确测算流域水系内污染物跨界通量及其时空分布是进行流域水环境风险预警和风险管理的重要前提之一。针对目前多种污染物通量测算方法在进行污染物年通量估算时,结果不确定性大这一突出问题,以流域水质监测站年内逐日流量、悬浮颗粒物监测数据作为悬浮颗粒物年通量参考值,基于以月、半月、周为周期的监测策略,将逐日流量、悬浮颗粒物监测数据重新筛选抽样构造,由此,系统分析了不同流域集水面积、污染通量监测频次和目前常用通量估算方法对污染物年通量估算不确定性的影响。所得方法和结论可为进一步制定流域污染物通量的测算规范提供方法指引和技术支持。     

Physically Based,Hydrologic Model Results Based on Three Precipitation Products1

Abstract: The main objective of the study is to examine the accuracy of and differences among simulated streamflows driven by rainfall estimates from a network of 22 rain gauges spread over a 2,170 km² watershed, NEXRAD Stage III radar data, and Tropical Rainfall Measuring Mission (TRMM) 3B42 satellite data. The Gridded Surface Subsurface Hydrologic Analysis (GSSHA), a physically based, distributed parameter, grid‐structured, hydrologic model, was used to simulate the June‐2002 flooding event in the Upper Guadalupe River watershed in south central Texas. There were significant differences between the rainfall fields estimated by the three types of measurement technologies. These differences resulted in even larger differences in the simulated hydrologic response of the watershed. In general, simulations driven by radar rainfall yielded better results than those driven by satellite or rain‐gauge estimates. This study also presents an overview of effects of land cover changes on runoff and stream discharge. The results demonstrate that, for major rainfall events similar to the 2002 event, the effect of urbanization on the watershed in the past two decades would not have made any significant effect on the hydrologic response. The effect of urbanization on the hydrologic response increases as the size of the rainfall event decreases.     

Fecal Coliform Export From Four Coastal North Carolina Areas1

Abstract: Fecal coliform (FC) bacteria in coastal waters impair the use of these waters for shellfish harvesting and recreation. This study was designed to quantify and compare FC levels and export in two coastal watersheds with different land uses. Continuous monitoring of rainfall and discharge at three sites in the Jumping Run Creek watershed and one site in the Pettiford Creek watershed were conducted during a 4.5‐year period. Primary land use in the drainage area of one of the three Jumping Run Creek sites is low density industrial, while the other two are residential. Land use in the Pettiford Creek watershed is managed national forest. Nonstorm or base‐flow grab and flow‐proportional storm‐event samples were collected and analyzed for turbidity, conductivity, suspended sediment, nitrogen, phosphorus, and FC. Geometric mean FC levels for the Jumping Run Creek monitoring sites ranged from 593 to 2,096 mpn/100 ml, while the mean level at the Pettiford Creek site was 191 mpn/100 ml. Levels of most other parameters were greater in storm discharge from the Jumping Run Creek sites as compared to Pettiford Creek indicating that pollutant export from a watershed increases with development. Statistical analysis of the monitoring data suggested that FC levels in stormwater samples consistently increased with storm rainfall, but were not consistently correlated with any other parameter, including total suspended solids. Multivariate analysis indicated that the weekly FC export for each of the four sites was lowest during the December‐February quarter. Export was highest during the spring and summer at the Jumping Run Creek sites, while for the Pettiford Creek site, FC export was highest during September‐November. The cause of the seasonal variability was unknown but was thought to be associated with human activity in the watersheds.