Development and Application of the 2010 Chesapeake Bay Watershed Total Maximum Daily Load Model

The Phase 5.3 Watershed Model simulates the Chesapeake watershed land use, river flows, and the associated transport and fate of nutrient and sediment loads to the Chesapeake Bay. The Phase 5.3 Model is the most recent of a series of increasingly refined versions of a model that have been operational for more than two decades. The Phase 5.3 Model, in conjunction with models of the Chesapeake airshed and estuary, provides estimates of management actions needed to protect water quality, achieve Chesapeake water quality standards, and restore living resources. The Phase 5.3 Watershed Model tracks nutrient and sediment load estimates of the entire 166,000 km² watershed, including loads from all six watershed states. The creation of software systems, input datasets, and calibration methods were important aspects of the model development process. A community model approach was taken with model development and application, and the model was developed by a broad coalition of model practitioners including environmental engineers, scientists, and environmental managers. Among the users of the Phase 5.3 Model are the Chesapeake watershed states and local governments, consultants, river basin commissions, and universities. Development and application of the model are described, as well as key scenarios ranging from high nutrient and sediment load conditions if no management actions were taken in the watershed, to low load estimates of an all‐forested condition.     

Estimating Contributions of Nitrate and Herbicides From Groundwater to Headwater Streams,Northern Atlantic Coastal Plain,United States1

Abstract: Groundwater transport often complicates understanding of surface‐water contamination. We estimated the regional flux of nitrate and selected herbicides from groundwater to nontidal headwater streams of the Atlantic Coastal Plain (New Jersey through North Carolina) based on late‐winter or spring base‐flow samples from 174 streams. Sampled streams were selected randomly, and flux estimates are based on resulting population estimates rather than on empirical models, which have been used previously for similar estimates. Base‐flow flux in the estimated 8,834 headwater streams of the study area are an estimated 21,200 kg/day of nitrate (as N) and 5.83, 0.565, and 20.7 kg/day of alachlor, atrazine, and metolachlor (and selected degradates), respectively. Base‐flow flux of alachlor and metolachlor is <3% of the total base‐flow flux of those compounds plus degradates. Base‐flow flux of nitrate and herbicides as a percentage of applications is typically highest in well‐drained areas and lowest in areas with abundant poor drainage and anoxic conditions. In Coastal Plain watersheds of Albemarle and Pamlico Sounds, <2% of applied nitrogen reaches headwater streams as base flow. On the Delmarva Peninsula part of the Chesapeake Bay watershed, however, more than 10% of such applications are transported through groundwater to streams, and base‐flow nitrate flux represents 70% of total nitrogen flux in headwater streams.     

Long‐Term Trends of Nutrients and Sediment from the Nontidal Chesapeake Watershed: An Assessment of Progress by River and Season

To assess historical loads of nitrogen (N), phosphorus (P), and suspended sediment (SS) from the nontidal Chesapeake Bay watershed (NTCBW), we analyzed decadal seasonal trends of flow‐normalized loads at the fall‐line of nine major rivers that account for >90% of NTCBW flow. Evaluations of loads by season revealed N, P, and SS load magnitudes have been highest in January‐March and lowest in July‐September, but the temporal trends have followed similar decadal‐scale patterns in all seasons, with notable exceptions. Generally, total N (TN) load has dropped since the late 1980s, but particulate nutrients and SS have risen since the mid‐1990s. The majority of these rises were from Susquehanna River and relate to diminished net trapping at the Conowingo Reservoir. Substantial rises in SS were also observed, however, in other rivers. Moreover, the summed rise in particulate P load from other rivers is of similar magnitude as from Susquehanna. Dissolved nutrient loads have dropped in the upland (Piedmont and above) rivers, but risen in two small rivers in the Coastal Plain affected by lagged groundwater input. In addition, analysis of fractional contributions revealed consistent N trends across the upland watersheds. Finally, total N:total P ratios have declined in most rivers, suggesting the potential for changes in nutrient limitation. Overall, this integrated study of historical data highlights the value of maintaining long‐term monitoring at multiple watershed locations.     

Effect of Watershed Subdivision and Filter Width on SWAT Simulation of a Coastal Plain Watershed1

Cho, Jaepil, Richard R. Lowrance, David D. Bosch, Timothy C. Strickland, Younggu Her, and George Vellidis, 2010. Effect of Watershed Subdivision and Filter Width on SWAT Simulation of a Coastal Plain Watershed. Journal of the American Water Resources Association (JAWRA) 46(3):586-602. DOI: 10.1111/j.1752-1688.2010.00436.x Abstract: The Soil and Water Assessment Tool (SWAT) does not fully simulate riparian buffers, but has a simple filter function that is responsive to filter strip width (FILTERW). The objectives of this study were to (1) evaluate SWAT hydrology and water quality response to changes in watershed subdivision levels and different FILTERW configurations and (2) provide guidance for selecting appropriate watershed subdivision for model runs that include the riparian buffer feature through the FILTERW parameter. Watershed subdivision level is controlled by the critical source area (CSA) which defines the minimum drainage area required to form the origin of a stream. SWAT was calibrated on a 15.7 km² subdrainage within the Little River Experimental Watershed, Georgia. The calibrated parameter set was applied to 32 watershed configurations consisting of four FILTERW representations for each of eight CSA levels. Streamflow predictions were stable regardless of watershed subdivision and FILTERW configuration. Predicted sediment and nutrient loads from upland areas decreased as CSA increased when spatial variations of riparian buffers are considered. Sediment and nutrient yield at the watershed outlet was responsive to different combinations of CSA and FILTERW depending on selected in-stream processes. CSA ranges which provide stable sediment and nutrient yields at the watershed outlet was suggested for avoiding significant modifications in selected parameter set.     

Water Quality Functions of Riparian Forest Buffers in Chesapeake Bay Watersheds

/ Maryland, Virginia, and Pennsylvania, USA, have agreed to reduce nutrient loadings to Chesapeake Bay by 40% by the year 2000. This requires control of nonpoint sources of nutrients, much of which comes from agriculture. Riparian forest buffer systems (RFBS) provide effective control of nonpoint source (NPS) pollution in some types of agricultural watersheds. Control of NPS pollution is dependent on the type of pollutant and the hydrologic connection between pollution sources, the RFBS, and the stream. Water quality improvements are most likely in areas of where most of the excess precipitation moves across, in, or near the root zone of the RFBS. In areas such as the Inner Coastal Plain and Piedmont watersheds with thin soils, RFBS should retain 50%-90% of the total loading of nitrate in shallow groundwater, sediment in surface runoff, and total N in both surface runoff and groundwater. Retention of phosphorus is generally much less. In regions with deeper soils and/or greater regional groundwater recharge (such as parts of the Piedmont and the Valley and Ridge), RFBS water quality improvements are probably much less. The expected levels of pollutant control by RFBS are identified for each of nine physiographic provinces of the Chesapeake Bay Watershed. Issues related to of establishment, sustainability, and management are also discussed.KEY WORDS: Riparian forest buffers; Chesapeake Bay; Nonpoint source pollution; Nitrogen; Phosphorus; Sediment     

Integration of SWAT and HSPF for Simulation of Sediment Sources in Legacy Sediment‐Impacted Agricultural Watersheds

A total maximum daily load for the Chesapeake Bay requires reduction in pollutant load from sources within the Bay watersheds. The Conestoga River watershed has been identified as a major source of sediment load to the Bay. Upland loads of sediment from agriculture are a concern; however, a large proportion of the sediment load in the Conestoga River has been linked to scour of legacy sediment associated with historic millpond sites. Clarifying this distinction and identifying specific segments associated with upland vs. channel sources has important implications for future management. In order to address this important question, we combined the strengths of two widely accepted watershed management models — Soil and Water Assessment Tool (SWAT) for upland agricultural processes, and Hydrologic Simulation Program FORTRAN (HSPF) for instream fate and transport — to create a novel linked modeling system to predict sediment loading from critical sources in the watershed including upland and channel sources, and to aid in targeted implementation of management practices. The model indicates approximately 66% of the total sediment load is derived from instream sources, in agreement with other studies in the region and can be used to support identification of these channel source segments vs. upland source segments, further improving targeted management. The innovated linked SWAT‐HSPF model implemented in this study is useful for other watersheds where both upland agriculture and instream processes are important sources of sediment load.     

Application of SPARROW Modeling to Understanding Contaminant Fate and Transport from Uplands to Streams

Understanding spatial variability in contaminant fate and transport is critical to efficient regional water‐quality restoration. An approach to capitalize on previously calibrated spatially referenced regression (SPARROW) models to improve the understanding of contaminant fate and transport was developed and applied to the case of nitrogen in the 166,000 km² Chesapeake Bay watershed. A continuous function of four hydrogeologic, soil, and other landscape properties significant (α = 0.10) to nitrogen transport from uplands to streams was evaluated and compared among each of the more than 80,000 individual catchments (mean area, 2.1 km²) in the watershed. Budgets (including inputs, losses or net change in storage in uplands and stream corridors, and delivery to tidal waters) were also estimated for nitrogen applied to these catchments from selected upland sources. Most (81%) of such inputs are removed, retained, or otherwise processed in uplands rather than transported to surface waters. Combining SPARROW results with previous budget estimates suggests 55% of this processing is attributable to denitrification, 23% to crop or timber harvest, and 6% to volatilization. Remaining upland inputs represent a net annual increase in landscape storage in soils or biomass exceeding 10 kg per hectare in some areas. Such insights are important for planning watershed restoration and for improving future watershed models.     

Urbanization and the Loss of Resource Lands in the Chesapeake Bay Watershed

We made use of land cover maps, and land use change associated with urbanization, to provide estimates of the loss of natural resource lands (forest, agriculture, and wetland areas) across the 168,000 km² Chesapeake Bay watershed. We conducted extensive accuracy assessments of the satellite-derived maps, most of which were produced by us using widely available multitemporal Landsat imagery. The change in urbanization was derived from impervious surface area maps (the built environment) for 1990 and 2000, from which we estimated the loss of resource lands that occurred during this decade. Within the watershed, we observed a 61% increase in developed land (from 5,177 to 8,363 km²). Most of this new development (64%) occurred on agricultural and grasslands, whereas 33% occurred on forested land. Some smaller municipalities lost as much as 17% of their forest lands and 36% of their agricultural lands to development, although in the outlying counties losses ranged from 0% to 1.4% for forests and 0% to 2.6% for agriculture. Fast-growing urban areas surrounded by forested land experienced the most loss of forest to impervious surfaces. These estimates could be used for the monitoring of the impacts of development across the Chesapeake Bay watershed, and the approach has utility for other regions nationwide. In turn, the results and the approach can help jurisdictions set goals for resource land protection and acquisition that are consistent with regional restoration goals.     

Flow,Organic, and Inorganic Sediment Yields from a Channelized Watershed in the South Carolina Lower Coastal Plain

Many small streams in coastal watersheds in the southeastern United States are modified for agricultural, residential, and commercial development. In the South Carolina Lower Coastal Plain, low‐relief topography and a shallow water table make stream channelization ubiquitous. To quantify the impacts of urbanization and stream channelization, we measured flow and sediment from an urbanizing watershed and a small forested watershed. Flow and sediment export rates were used to infer specific yields from forested and nonforested regions of the urbanizing watershed. Study objectives were to: (1) quantify the range of runoff‐to‐rainfall ratios; (2) quantify the range of specific sediment yields; (3) characterize the quantity and quality of particulate matter exported; and (4) estimate sediment yield attributable to agriculture, development, and channelization activities in the urbanizing watershed. Our results showed that the urban watershed exported over five times more sediment per unit area compared with the forested watershed. Sediment concentration was related to flow flashiness in the urban watershed and to flow magnitude in the forested watershed. Sediments from the forested watershed were dominated by organic matter, whereas mineral matter dominated sediment from the urban stream. Our results indicated that a significant shift in sediment quality and quantity are likely to occur as forested watersheds are transformed by urbanization in coastal South Carolina.     

REHABILITATION OF WATERSHEDS WITH INCISING CHANNELS1

ABSTRACT: Channel incision is a pervasive problem that threatens infrastructure, destroys arable land, and degrades environmental resources. A program initiated in 1983 is developing technology for rehabilitation of watersheds with erosion and sedimentation problems caused by incision. Demonstration projects are located in 15 watersheds in the hills of northwest Mississippi. Watershed sizes range from 0.89 to 1,590 km², and measured suspended sediment yields average about 1,100 t km^-2-yr^-1. Water quality is generally adequate to support aquatic organisms, but physical habitat conditions are poor. Rehabilitation measures, which are selected and laid out using a subjective integration of hydraulic and geotechnical stability analyses, include grade controls, bank protection, and small reservoirs. Aquatic habitat studies indicate that stone-protected stilling basins below grade-control weirs and habitats associated with drop pipes and stone spur dikes are assets to erosion-damaged streams. Additional recovery of habitat resources using modified stone stabilization designs, woody vegetation plantings, and reservoir outlets designed to provide non-zero minimum flows is under investigation.     

Improvement in HSPF's Low‐Flow Predictions by Implementation of a Power Law Groundwater Storage‐Discharge Relationship

We have enhanced the ability of a widely used watershed model, Hydrologic Simulation Program — FORTRAN (HSPF), to predict low flows by reconfiguring the algorithm that simulates groundwater discharge. During dry weather periods, flow in most streams consists primarily of base flow, that is, groundwater discharged from underlying aquifers. In this study, HSPF's groundwater storage‐discharge relationship is changed from a linear to a more general nonlinear relationship which takes the form of a power law. The nonlinear algorithm is capable of simulating streamflow recession curves that have been found in some studies to better match observed dry weather hydrographs. The altered version of HSPF is implemented in the Chesapeake Bay Program's Phase 5 Model, an HSPF‐based model that simulates nutrient and sediment loads to the Chesapeake Bay, and is tested in the upper Potomac River basin, a 29,950 km² drainage area that is part of the Bay watershed. The nonlinear relationship improved median Nash‐Sutcliffe efficiencies for log daily flows at the model's 45 calibration points. Mean absolute percent error on low‐flow days dropped in five major Potomac River tributaries by up to 12 percentage points, and in the Potomac River itself by four percentage points, where low‐flow days were defined as days when observed flows were in the lowest 5th percentile range. Percent bias on low‐flow days improved by eight percentage points in the Potomac River, from ?11 to ?3%.     

Stream Condition in Piedmont Streams with Restored Riparian Buffers in the Chesapeake Bay Watershed1

Orzetti, Leslie L., R. Christian Jones, and Robert F. Murphy, 2010. Stream Condition in Piedmont Streams with Restored Riparian Buffers in the Chesapeake Bay Watershed. Journal of the American Water Resources Association (JAWRA) 46(3):473-485. DOI: 10.1111/j.1752-1688.2009.00414.x Abstract: This study tested the efficacy of restored forest riparian buffers along streams in the Chesapeake Bay watershed by examining habitat, selected water quality variables, and benthic macroinvertebrate community metrics in 30 streams with buffers ranging from zero to greater than 50 years of age. To assess water quality we measured in situ parameters (temperature, dissolved oxygen, and conductivity) and laboratory-analyzed grab samples (soluble reactive phosphorus, total phosphorus, nitrate, ammonium, and total suspended solids). Habitat conditions were scored using the Environmental Protection Agency Rapid Bioassessment Protocols for high gradient streams. Benthic macroinvertebrates were quantified using pooled riffle/run kick samples. Results showed that habitat, water quality, and benthic macroinvertebrate metrics generally improved with age of restored buffer. Habitat scores appeared to stabilize between 10 and 15 years of age and were driven mostly by epifaunal substrate availability, sinuosity, embeddedness, and velocity depth regime. Benthic invertebrate taxa richness, percent Ephemeroptera, Plecoptera, Trichoptera minus hydropsychids (%EPT minus H), % Ephemeroptera, and the Family Biotic Index were among the metrics which improved with age of buffer zone. Results are consistent with the hypothesis that forest riparian buffers enhance instream habitat, water quality, and resulting benthic macroinvertebrate communities with noticeable improvements occurring within 5-10 years postrestoration, leading to conditions approaching those of long established buffers within 10-15 years of restoration.     

Cropland Riparian Buffers throughout Chesapeake Bay Watershed: Spatial Patterns and Effects on Nitrate Loads Delivered to Streams

We used statistical models to provide the first empirical estimates of riparian buffer effects on the cropland nitrate load to streams throughout the Chesapeake Bay watershed. For each of 1,964 subbasins, we quantified the 1990 prevalence of cropland and riparian buffers. Cropland was considered buffered if the topographic flow path connecting it to a stream traversed a streamside forest or wetland. We applied a model that predicts stream nitrate concentration based on physiographic province and the watershed proportions of unbuffered and buffered cropland. We used another model to predict annual streamflow based on precipitation and temperature, and then multiplied the predicted flows and concentrations to estimate 1990 annual nitrate loads. Across the entire Chesapeake watershed, croplands released 92.3 Gg of nitrate nitrogen, but 19.8 Gg of that was removed by riparian buffers. At most, 29.4 Gg more might have been removed if buffer gaps were restored so that all cropland was buffered. The other 43.1 Gg of cropland load cannot be addressed with riparian buffers. The Coastal Plain physiographic province provided 52% of the existing buffer reduction of Bay‐wide nitrate loads and 36% of potential additional removal from buffer restoration in cropland buffer gaps. Existing and restorable nitrate removal in buffers were lower in the other three major provinces because of less cropland, lower buffer prevalence, and lower average buffer nitrate removal efficiency.     

SPARROW Modeling of Nitrogen Sources and Transport in Rivers and Streams of California and Adjacent States,U.S.

The SPARROW (SPAtially Referenced Regressions On Watershed attributes) model was used to evaluate the spatial distribution of total nitrogen (TN) sources, loads, watershed yields, and factors affecting transport and decay in the stream network of California and portions of adjacent states for the year 2002. The two major TN sources to local catchments on a mass basis were fertilizers and manure (51.7%) and wastewater discharge (15.9%). Other sources contributed < 12%. Fertilizer use is widespread in the Central Valley region of California, and also important in several other regions because of the diversity of California agriculture. Precipitation, sand content of surficial soils, wetlands, and tile drains were important for TN movement to stream reaches. Median streamflow in the study area is about 0.04 m³/s. Aquatic losses of nitrogen were found to be most important in intermittent and small to medium sized streams (0.2‐14 m³/s), while larger streams showed less loss, and therefore are important for TN transport. Nitrogen loss in reservoirs was found to be insignificant, possibly because most of the larger ones are located upstream of nitrogen sources. The model was used to show loadings, sources, and tributary inputs to several major rivers. The information provided by the SPARROW model is useful for determining both the major sources contributing nitrogen to streams and the specific tributaries that transport the load.     

Response of Nitrogen Loading to the Chesapeake Bay to Source Reduction and Land Use Change Scenarios: A SPARROW‐Informed Analysis

In response to concerns regarding the health of streams and receiving waters, the United States Environmental Protection Agency established a total maximum daily load for nitrogen in the Chesapeake Bay watershed for which practices must be in place by 2025 resulting in an expected 25% reduction in load from 2009 levels. The response of total nitrogen (TN) loads delivered to the Bay to nine source reduction and land use change scenarios was estimated using a Spatially Referenced Regression on Watershed Attributes model. The largest predicted reduction in TN load delivered to the Bay was associated with a scenario in which the mass of TN as fertilizer applied to agricultural lands was decreased. A 25% decrease in the mass of TN applied as fertilizer resulted in a predicted reduction in TN loading to the Bay of 11.3%, which was 2.5–5 times greater than the reductions predicted by other scenarios. Eliminating fertilizer application to all agricultural land in the watershed resulted in a predicted reduction in TN load to the Bay of 45%. It was estimated that an approximate 25% reduction in TN loading to the Bay could be achieved by eliminating fertilizer applied to the 7% of subwatersheds contributing the greatest fertilizer‐sourced TN loads to the Bay. These results indicate that management strategies aimed at decreasing loading from a small number of subwatersheds may be effective for reducing TN loads to the Bay, and similar analyses are possible in other watersheds.     

Contribution of In‐Channel Processes to Sediment Yield of an Urbanizing Watershed1

Abstract: A study was conducted between September 2003 and September 2006 to obtain baseline sediment inventories and monitor sediment transport and storage along a 3.7 km length of the channel of Valley Creek within Valley Forge National Historical Park, Pennsylvania. Valley Creek is a tributary of the Schuylkill River and drains an urbanizing 60.6 km² watershed that currently has 18% impervious land cover. Numerous field methods were employed to measure the suspended sediment yield, longitudinal profile, cross‐sections, banklines, and particle size distribution of the streambed. Suspended sediment yield for the watershed was measured at a USGS gage located just upstream of the park boundary between July 2004 and July 2005, the period corresponding to field surveys of bank erosion and channel change. The estimated suspended sediment yield of 95.7 t/km²/year is representative of a year with unusually high discharge, including a storm event that produced a peak of 78 m³/s, the second highest discharge on record for the USGS gage. Based on the median annual streamflow for the 24 years of record at the USGS gage from 1983 to 2006, the median annual sediment yield is estimated to be closer to 34 t/km²/year, considerably lower than median and mean values for other sites within the region. The mass of silt, clay, and fine sand derived from bank erosion along the 3.7 km study reach during the field survey period accounts for an estimated 2,340 t, equivalent to about 43% of the suspended sediment load. The mass of fine sediment stored in the bed along the study reach was estimated at 1,500 t, with about 330 t of net erosion during the study period. Although bank erosion appears to be a potentially dominant source of sediment by comparison with annual suspended sediment load, bed sediment storage and potential for remobilization is of the same order of magnitude as the mass of sediment derived from bank erosion.     

Modeled Sources,Transport, and Accumulation of Dissolved Solids in Water Resources of the Southwestern United States1

Anning, David W., 2011. Modeled Sources, Transport, and Accumulation of Dissolved Solids in Water Resources of the Southwestern United States. Journal of the American Water Resources Association (JAWRA) 47(5):1087‐1109. DOI: 10.1111/j.1752‐1688.2011.00579.x Abstract: Information on important source areas for dissolved solids in streams of the southwestern United States, the relative share of deliveries of dissolved solids to streams from natural and human sources, and the potential for salt accumulation in soil or groundwater was developed using a SPAtially Referenced Regressions On Watershed attributes model. Predicted area‐normalized reach‐catchment delivery rates of dissolved solids to streams ranged from <10 (kg/year)/km² for catchments with little or no natural or human‐related solute sources in them to 563,000 (kg/year)/km² for catchments that were almost entirely cultivated land. For the region as a whole, geologic units contributed 44% of the dissolved‐solids deliveries to streams and the remaining 56% of the deliveries came from the release of solutes through irrigation of cultivated and pasture lands, which comprise only 2.5% of the land area. Dissolved‐solids accumulation is manifested as precipitated salts in the soil or underlying sediments, and (or) dissolved salts in soil‐pore or sediment‐pore water, or groundwater, and therefore represents a potential for aquifer contamination. Accumulation rates were <10,000 (kg/year)/km² for many hydrologic accounting units (large river basins), but were more than 40,000 (kg/year)/km² for the Middle Gila, Lower Gila‐Agua Fria, Lower Gila, Lower Bear, Great Salt Lake accounting units, and 247,000 (kg/year)/km² for the Salton Sea accounting unit.     

Monitored and Modeled Correlations of Sediment and Nutrients with Chesapeake Bay Water Clarity

This article analyzes the correlations of the observed and modeled light attenuation coefficient, Kd, with in situ total suspended solids (TSS) and chlorophyll‐a concentrations in Chesapeake Bay (CB) tidal waters, and with sediment and nutrient loads from the Chesapeake watershed. Light attenuation is closely related to in situ TSS and chlorophyll‐a concentrations, however, the strength of the correlation differs among the CB segments. There are distinct differences between saline and tidal fresh segments in the main Bay, but less distinction among saline and tidal fresh segments in the tidal tributaries. The correlation between Kd with sediment and nutrient loads is complicated by the lag times of TSS and the chlorophyll‐a responses to reductions in nutrient and sediment loads from the watershed, and also due to the diverse load sources. Three sets of model sensitivity scenarios were performed with: (1) differential sediment and nutrient loads; (2) selective sediment source types; and (3) geographically isolated inputs. The model results yield similar findings as those based on observed data and provide information regarding the effect of sediment on specific water bodies. Based on the model results a method was developed to determine sediment and nutrient load reductions needed to achieve the water clarity standards of the CB segments.     

Legacy Effects of Colonial Millponds on Floodplain Sedimentation,Bank Erosion,and Channel Morphology,Mid‐Atlantic,USA1

Abstract: Many rivers and streams of the Mid‐Atlantic Region, United States (U.S.) have been altered by postcolonial floodplain sedimentation (legacy sediment) associated with numerous milldams. Little Conestoga Creek, Pennsylvania, a tributary to the Susquehanna River and the Chesapeake Bay, is one of these streams. Floodplain sedimentation rates, bank erosion rates, and channel morphology were measured annually during 2004‐2007 at five sites along a 28‐km length of Little Conestoga Creek with nine colonial era milldams (one dam was still in place in 2007). This study was part of a larger cooperative effort to quantify floodplain sedimentation, bank erosion, and channel morphology in a high sediment yielding region of the Chesapeake Bay watershed. Data from the five sites were used to estimate the annual volume and mass of sediment stored on the floodplain and eroded from the banks for 14 segments along the 28‐km length of creek. A bank and floodplain reach based sediment budget (sediment budget) was constructed for the 28 km by summing the net volume of sediment deposited and eroded from each segment. Mean floodplain sedimentation rates for Little Conestoga Creek were variable, with erosion at one upstream site (?5 mm/year) to deposition at the other four sites (highest = 11 mm/year) despite over a meter of floodplain aggradation from postcolonial sedimentation. Mean bank erosion rates range between 29 and 163 mm/year among the five sites. Bank height increased 1 m for every 10.6 m of channel width, from upstream to downstream (R² = 0.79, p < 0.0001) resulting in progressively lowered hydraulic connectivity between the channel and the floodplain. Floodplain sedimentation and bank erosion rates also appear to be affected by the proximity of the segments to one existing milldam, which promotes deposition upstream and scouring downstream. The floodplain and bank along the 28‐km reach produced a net mean sediment loss of 5,634 Mg/year for 2004‐2007, indicating that bank erosion was exceeding floodplain sedimentation. In particular, the three segments between the existing dam and the confluence with the Conestoga River (32% of the studied reach) account for 97% of the measured net sediment budget. Future research directed at understanding channel equilibria should facilitate efforts to reduce the sediment impacts of dam removal and legacy sediment.     

EFFECTS OF LAND COVER AND GEOLOGY ON STREAM CHEMISTRY IN WATERSHEDS OF CHESAPEAKE BAY1

ABSTRACT: We measured the base‐flow stream chemistry in all the major physiographic provinces of the Chesapeake Bay drainage basin. The spatial variation of stream chemistry was closely related to differences in geology and land cover among the sampled watersheds. Some stream chemistry variables were strongly affected by geological settings in the watersheds while others were more influenced by land cover. The effects of land cover differed among chemical constituents and regions. Concentrations of Ca²⁺, Mg²⁺, pH, total alkalinity, and conductivity were mainly functions of carbonate bedrock, especially in the Great Valley. Nitrate‐N and total dissolved N were closely related to cropland and increased as the percentage of cropland increased. The rate of increase varied from region to region with the highest in the Piedmont. Na⁺ and Cl^? were mainly affected by the percentage of developed area in a watershed, especially in the Coastal Plain and Piedmont. We observed no significant effects of region or land cover on species of phosphorus because samples were collected under base flow conditions and only dissolved forms were measured. Dissolved silicate (DSi) was not related to any other water chemistry variables. DSi increased as developed area decreased and cropland increased in the Coastal Plain, but these patterns were reversed in the Piedmont. There was no consistent pattern in the spatial variation of land cover effects on the reduced forms of N, dissolved organic P, dissolved organic matter, and K⁺.