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.     

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.     

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.     

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.     

Toward Explaining Nitrogen and Phosphorus Trends in Chesapeake Bay Tributaries, 1992–2012

Understanding trends in stream chemistry is critical to watershed management, and often complicated by multiple contaminant sources and landscape conditions changing over varying time scales. We adapted spatially referenced regression (SPARROW) to infer causes of recent nutrient trends in Chesapeake Bay tributaries by relating observed fluxes during 1992, 2002, and 2012 to contemporary inputs and watershed conditions. The annual flow‐normalized nitrogen flux to the bay from its watershed declined by 14% to 127,000 Mg (metric tons) between 1992 and 2012, due primarily (more than 80% of the decline) to reduced point sources. The remainder of the decline was due to reduced atmospheric deposition (13%) and urban nonpoint sources. Agricultural inputs, which contribute most nitrogen to the bay, changed little, although trends in the average nitrogen yield (flux per unit area) from cropland and pasture to streams in some settings suggest possible effects of evolving nutrient applications or other land management practices. Point sources of phosphorus to local streams declined by half between 1992 and 2012, while nonpoint inputs were relatively unchanged. Annual phosphorus delivery to the bay increased by 9% to 9,570 Mg between 1992 and 2012, however, due mainly to reduced retention in the Susquehanna River at Conowingo Reservoir.     

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.     

A Geospatial Methodology to Identify Locations of Concentrated Runoff from Agricultural Fields

A geospatial methodology has been developed that utilizes high resolution lidar‐derived DEMs to help track runoff from agricultural fields and identify areas of potential concentrated flow through vegetated riparian areas in the Coastal Plain of Virginia. Points of concentrated flow are identified across 74 agricultural fields within the Virginia portion of the Chesapeake Bay watershed. On average, 70% of the surface area of the agricultural fields analyzed drains through less than 20 m of the field margin, and on average 81% of the field surface area drains through 1% or less of the field margin. Within the riparian buffer, locations that were predicted by the geospatial model to have high levels of concentrated flow were found to exhibit evidence of channelization. Results indicate that flow concentration and channelized flow through vegetated riparian areas may be common along the margin of agricultural fields, resulting in vegetated riparian areas that are less effective at sediment trapping than assumed. Additional results suggest that the regulations governing the location and width of vegetated riparian may not be sufficient to achieve goals for reducing sediment and nutrient runoff from nonpoint agricultural sources. Combined with the increasing availability of lidar‐derived DEMs, the geospatial model presented has the potential to advance management practices aimed at reducing nonpoint source pollution leaving agricultural fields.     

Total Maximum Daily Load Criteria Assessment Using Monitoring and Modeling Data

Applications of Total Maximum Daily Load (TMDL) criteria for complex estuarine systems like Chesapeake Bay have been limited by difficulties in estimating precisely how changes in input loads will impact ambient water quality. A method to deal with this limitation combines the strengths of the Chesapeake Bay's Water Quality Sediment Transport Model (WQSTM), which simulates load response, and the Chesapeake Bay Program's robust historical monitoring dataset. The method uses linear regression to apply simulated relative load responses to historical observations of water quality at a given location and time. Steps to optimize the application of regression analysis were to: (1) determine the best temporal and spatial scale for applying the WQSTM scenarios, (2) determine whether the WQSTM method remained valid with significant perturbation from calibration conditions, and (3) evaluate the need for log transformation of both dissolved oxygen (DO) and chlorophyll a (CHL) datasets. The final method used simple linear regression at the single month, single WQSTM grid cell scale to quantify changes in DO and CHL resulting from simulated load reduction scenarios. The resulting linear equations were applied to historical monitoring data to produce a set of “scenario‐modified” DO or CHL concentration estimates. The utility of the regression method was validated by its ability to estimate progressively increasing attainment in support of the 2010 Chesapeake Bay TMDL.     

Sources of Suspended-Sediment Flux in Streams of the Chesapeake Bay Watershed: A Regional Application of the SPARROW Model1

Brakebill, John W., Scott W. Ator, and Gregory E. Schwarz, 2010. Sources of Suspended-Sediment Flux in Streams of the Chesapeake Bay Watershed: A Regional Application of the SPARROW Model. Journal of the American Water Resources Association (JAWRA) 46(4): 757-776. DOI: 10.1111/j.1752-1688.2010.00450.x Abstract: We describe the sources and transport of fluvial suspended sediment in nontidal streams of the Chesapeake Bay watershed and vicinity. We applied SPAtially Referenced Regressions on Watershed attributes, which spatially correlates estimated mean annual flux of suspended sediment in nontidal streams with sources of suspended sediment and transport factors. According to our model, urban development generates on average the greatest amount of suspended sediment per unit area (3,928 Mg/km²/year), although agriculture is much more widespread and is the greatest overall source of suspended sediment (57 Mg/km²/year). Factors affecting sediment transport from uplands to streams include mean basin slope, reservoirs, physiography, and soil permeability. On average, 59% of upland suspended sediment generated is temporarily stored along large rivers draining the Coastal Plain or in reservoirs throughout the watershed. Applying erosion and sediment controls from agriculture and urban development in areas of the northern Piedmont close to the upper Bay, where the combined effects of watershed characteristics on sediment transport have the greatest influence may be most helpful in mitigating sedimentation in the bay and its tributaries. Stream restoration efforts addressing floodplain and bank stabilization and incision may be more effective in smaller, headwater streams outside of the Coastal Plain.     

Computing Atmospheric Nutrient Loads to the Chesapeake Bay Watershed and Tidal Waters

Application of integrated Chesapeake Bay models of the airshed, watershed, and estuary support air and water nitrogen controls in the Chesapeake. The models include an airshed model of the Mid‐Atlantic region which tracks the estimated atmospheric deposition loads of nitrogen to the watershed, tidal Bay, and adjacent coastal ocean. The three integrated models allow tracking of the transport and fate of nitrogen air emissions, including deposition in the Chesapeake watershed, the subsequent uptake, transformation, and transport to Bay tidal waters, and their ultimate influence on Chesapeake water quality. This article describes the development of the airshed model, its application to scenarios supporting the Chesapeake Total Maximum Daily Load (TMDL), and key findings from the scenarios. Key findings are that the atmospheric deposition loads are among the largest input loads of nitrogen in the watershed, and that the indirect nitrogen deposition loads to the watershed, which are subsequently delivered to the Bay are larger than the direct loads of atmospheric nitrogen deposition to Chesapeake tidal waters. Atmospheric deposition loads of nitrogen deposited in coastal waters, which are exchanged with the Chesapeake, are also estimated. About half the atmospheric deposition loads of nitrogen originate from outside the Chesapeake watershed. For the first time in a TMDL, the loads of atmospheric nitrogen deposition are an explicit part of the TMDL load reductions.     

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.     

Monitoring Regional Riparian Forest Cover Change Using Stratified Sampling and Multiresolution Imagery1

Claggett, Peter R., Judy A. Okay, and Stephen V. Stehman, 2010. Monitoring Regional Riparian Forest Cover Change Using Stratified Sampling and Multiresolution Imagery. Journal of the American Water Resources Association (JAWRA) 46(2):334-343. DOI: 10.1111/j.1752-1688.2010.00424.x Abstract: The Chesapeake Bay watershed encompasses 165,760 km² of land area with 464,098 km of rivers and streams. As part of the Chesapeake Bay restoration effort, state and federal partners have committed to restoring 26,000 miles (41,843 km) of riparian forest buffers. Monitoring trends in riparian forest buffers over large areas is necessary to evaluate the efficacy of these restoration efforts. A sampling approach for estimating change in riparian forest cover from 1993/1994 to 2005 was developed and implemented in Anne Arundel County, Maryland, to exemplify a method that could be applied throughout the Bay watershed. All stream reaches in the county were stratified using forest cover change derived from Landsat imagery. A stratified random sample of 219 reaches was selected and forest cover change within the riparian buffer of each sampled reach was interpreted from high-resolution aerial photography. The estimated footprint of gross change in riparian forest cover (i.e., the sum of gross gain and gross loss) for the county was 1.83% (SE = 0.22%). Stratified sampling taking advantage of a priori knowledge of locations of change proved to be a practical and efficient protocol for estimating riparian forest buffer change at the county scale and the protocol would readily extend to much broader scale monitoring.     

Spatial and Temporal Patterns of Low Streamflow and Precipitation Changes in the Chesapeake Bay Watershed

Spatial and temporal patterns in low streamflows were investigated for 183 streamgages located in the Chesapeake Bay Watershed for the period 1939–2013. Metrics that represent different aspects of the frequency and magnitude of low streamflows were examined for trends: (1) the annual time series of seven‐day average minimum streamflow, (2) the scaled average deficit at or below the 2% mean daily streamflow value relative to a base period between 1939 and 1970, and (3) the annual number of days below the 2% threshold. Trends in these statistics showed spatial cohesion, with increasing low streamflow volume at streamgages located in the northern uplands of the Chesapeake Bay Watershed and decreasing low streamflow volume at streamgages in the southern part of the watershed. For a small subset of streamgages (12%), conflicting trend patterns were observed between the seven‐day average minimum streamflow and the below‐threshold time series and these appear to be related to upstream diversions or the influence of reservoir‐influenced streamflows in their contributing watersheds. Using multivariate classification techniques, mean annual precipitation and fraction of precipitation falling as snow appear to be broad controls of increasing and decreasing low‐flow trends. Further investigation of seasonal precipitation patterns shows summer rainfall patterns, driven by the Atlantic Multidecadal Oscillation, as the main driver of low streamflows in the Chesapeake Bay Watershed.     

Twenty‐One‐Year Simulation of Chesapeake Bay Water Quality Using the CE‐QUAL‐ICM Eutrophication Model

The CE‐QUAL‐ICM (Corps of Engineers Integrated Compartment Water Quality Model) eutrophication model was applied in a 21‐year simulation of Chesapeake Bay water quality, 1985‐2005. The eutrophication model is part of a larger model package and is forced, in part, by models of atmospheric deposition, watershed flows and loads, and hydrodynamics. Results from the model are compared with observations in multiple formats including time series plots, cumulative distribution plots, and statistical summaries. The model indicates only one long‐term trend in computed water quality: light attenuation deteriorates circa 1993 through the end of the simulation. The most significant result is the influence of physical processes, notably stratification and associated effects (e.g., anoxic volume), on computed water quality. Within the application period, physical effects are more important determinants of year‐to‐year variability in computed water quality than external loads.     

Interactions of warming and altered nutrient load timing on the phenology of oxygen dynamics in Chesapeake Bay

The effects of nutrient loading on estuaries are well studied, given the multitude of negative water quality and ecosystem effects that have been attributed to excess nitrogen and phosphorus. A current gap in this knowledge involves the sensitivity of seasonal cycles of estuarine biogeochemical processes to direct (warming) and indirect influences (nutrient load timing) of climate change. We used a coupled hydrologic–biogeochemical model to investigate changes in the phenology of hypoxia and related biogeochemical processes in Chesapeake Bay under three different hydrologic regimes. Shifts to earlier nutrient load timing during idealized simulations reduced the overall annual hypoxic volume, resulting from discernable, but relatively small reductions in phytoplankton biomass and both sediment and water-column respiration. Simulated increases in water temperature caused an increase in spring/early summer hypoxic volume associated with elevated respiration rates, but an associated exhaustion of organic matter in the early summer caused a decrease in late summer/fall hypoxic volume due to lowered respiration. Warming effects on hypoxia were larger than nutrient timing effects in scenarios where warming was restricted to spring and when it was applied to all months of the year. These idealized simulations begin the process of understanding the potential impacts of future climatic changes in the seasonal timing of key biogeochemical processes associated with eutrophication.     

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.     

EVALUATION OF COASTAL PLAIN CONSERVATION BUFFERS USING THE RIPARIAN ECOSYSTEM MANAGEMENT MODEL1

ABSTRACT: Riparian buffers are increasingly important as watershed management tools and are cost‐shared by programs such as Conservation Reserve that are part of the USDA Conservation Buffer Initiative. Riparian buffers as narrow as 4.6m (15ft) are eligible for cost‐share by USDA. The Riparian Ecosystem Management Model (REMM) provides a tool to judge water quality improvement by buffers and to set design criteria for nutrient and sediment load reduction. REMM was used for a Coastal Plain site to simulate 14 different buffers ranging from 4.6 m to 51.8 m (15 to 170 ft) with three different types of vegetation (hardwood trees, pine trees, and perennial grass) with two water and nutrient loads. The load cases were low sediment/low nutrient‐typical of a well managed agricultural field and low sediment/high nutrient‐typical of liquid manure application to perennial forage crops. Simulations showed that the minimum width buffer (4.6 m) was inadequate for control of nutrients under either load case. The minimum width buffer that is eligible for cost share assistance on a field with known water quality problems (10.7 m, 35 ft) was projected to achieve at least 50 percent reduction of N, P, and sediment in the load cases simulated.     

Deriving Chesapeake Bay Water Quality Standards

Achieving and maintaining the water quality conditions necessary to protect the aquatic living resources of the Chesapeake Bay and its tidal tributaries has required a foundation of quantifiable water quality criteria. Quantitative criteria serve as a critical basis for assessing the attainment of designated uses and measuring progress toward meeting water quality goals of the Chesapeake Bay Program partnership. In 1987, the Chesapeake Bay Program partnership committed to defining the water quality conditions necessary to protect aquatic living resources. Under section 303(c) of the Clean Water Act, States and authorized tribes have the primary responsibility for adopting water quality standards into law or regulation. The Chesapeake Bay Program partnership worked with U.S. Environmental Protection Agency to develop and publish a guidance framework of ambient water quality criteria with designated uses and assessment procedures for dissolved oxygen, water clarity, and chlorophyll a for Chesapeake Bay and its tidal tributaries in 2003. This article reviews the derivation of the water quality criteria, criteria assessment protocols, designated use boundaries, and their refinements published in six addendum documents since 2003 and successfully adopted into each jurisdiction's water quality standards used in developing the Chesapeake Bay Total Maximum Daily Load.     

Estimating the Extent of Impervious Surfaces and Turf Grass across Large Regions

The ability of researchers to accurately assess the extent of impervious and pervious developed surfaces, e.g., turf grass, using land‐cover data derived from Landsat satellite imagery in the Chesapeake Bay watershed is limited due to the resolution of the data and systematic discrepancies between developed land‐cover classes, surface mines, forests, and farmlands. Estimates of impervious surface and turf grass area in the Mid‐Atlantic, United States that were based on 2006 Landsat‐derived land‐cover data were substantially lower than estimates based on more authoritative and independent sources. New estimates of impervious surfaces and turf grass area derived using land‐cover data combined with ancillary information on roads, housing units, surface mines, and sampled estimates of road width and residential impervious area were up to 57 and 45% higher than estimates based strictly on land‐cover data. These new estimates closely approximate estimates derived from authoritative and independent sources in developed counties.     

Resource protection for waterbirds in Chesapeake bay

Many living resources in the Chesapeake Bay estuary have deteriorated over the past 50 years. As a result, many governmental committees, task forces, and management plans have been established. Most of the recommendations for implementing a bay cleanup focus on reducing sediments and nutrient flow into the watershed. We emphasize that habitat requirements other than water quality are necessary for the recovery of much of the bay's avian wildlife, and we use a waterbird example as illustration. Some of these needs are: (1) protection of fast-eroding islands, or creation of new ones by dredge deposition to improve nesting habitat for American black ducks(Anas rubripes), great blue herons(Ardea herodias), and other associated wading birds; (2) conservation of remaining brackish marshes, especially near riparian areas, for feeding black ducks, wading birds, and wood ducks(Aix sponsa); (3) establishment of sanctuaries in open-water, littoral zones to protect feeding and/or roosting areas for diving ducks such as canvasbacks(Aythya valisineria) and redheads(Aythya americana), and for bald eagles(Haliaeetus leucocephalus); and (4) limitation of disturbance by boaters around nesting islands and open-water feeding areas. Land (or water) protection measures for waterbirds need to include units at several different spatial scales, ranging from “points” (e.g., a colony site) to large-area resources (e.g., a marsh or tributary for feeding). Planning to conserve large areas of both land and water can be achieved following a biosphere reserve model. Existing interagency committees in the Chesapeake Bay Program could be more effective in developing such a model for wildlife and fisheries resources.