Seasonal oxygen depletion in Chesapeake Bay

The spring freshet increases density stratification in Chesapeake Bay and minimizes oxygen transfer from the surface to the deep layer so that waters below 10 m depth experiece oxygen depletion which may lead to anoxia during June to September. Respiration in the water of the deep layer is the major factor contributing to oxygen depletion. Benthic respiration seems secondary. Organic matter from the previous year which has settled into the deep layer during winter provides most of the oxygen demand but some new production in the surface layer may sink and thus supplement the organic matter accumulated in the deep layer.     

Modeling of wind-induced destratification in Chesapeake Bay

It has been observed that storms in early fall can result in top-to-bottom mixing of Chesapeake Bay. A three-dimensional, time-dependent circulation model is used to examine this destratification process for September 1983, when extensive current and hydrographic data were available. The model bay is forced at the surface by observed hourly winds, at the ocean boundary by observed hourly surface and bottom salinities and sea level fluctuations, and at the head by observed daily discharges for a 28-d period. A second-moment, turbulence-closure submodel, with no adjustments from previous applications to its requisite coefficients, is used to calculate the vertical turbulence mixing coefficients. Comparisons with data inside the model domain indicate relative errors of 7% to 14% for sea level, 7% to 35% for current, and 11% to 21% for salinity. The tidal portion of the spectrum is modeled better than the subtidal portion. The model is used to examine both the mechanisms of wind mixing and the temporal and spatial distribution of vertical mixing within the estuary. Wind-driven internal shear is shown to be a more effective mechanism of inducing destratification than turbulence generated at the surface. The model is also used to show that the vertical temperature inversion which occurs in the fall does not affect the timing of the destratification as much as its completeness. The distribution of mid-depth vertical mixing shows highly variable values in the mid-bay region, where wind-induced mixing is dominant. This suggests that the source of oxygen to mid-bay bottom waters is similarly variable. Vertical turbulence mixing coefficients of 10^?2 cm² s^?1 (background) to 10³ cm² s^?1 were needed to simulate the September period, indicating the need for time-variable mixing in models of dissolved and suspended estuarine constituents.     

System-wide submerged aquatic vegetation model for Chesapeake Bay

A predictive model of submerged aquatic vegetation (SAV) biomass is coupled to a eutrophication model of Chesapeake Bay. Domain of the model includes the mainstem of the bay as well as tidal portions of major embayments and tributaries. Three SAV communities are modeled: ZOSTERA, RUPPIA, and FRESHWATER. The model successfully computes the spatial distribution and abundance of SAV for the period 1985–1994. Spatial distribution is primarily determined by computed light attenuation. Sensivitity analysis to reductions in nutrient and solids loads indicates nutrient controls will enhance abundance primarily in areas that presently support SAV. Restoration of SAV to areas in which it does not presently exist requires solids controls, alone or in combination with nutrient controls. For regions in which SAV populations exist at the refuge level or greater, improvements in SAV abundance are expected within 2 to 10 years of load reductions. For regions in which no refuge population exists, recovery time is unpredictable and will depend on propagule supply.     

Fish biomass size spectra in Chesapeake Bay

Biomass size spectra of pelagic fish were modeled to describe community structure, estimate potential fish production, and delineate trophic relationships in Chesapeake Bay. Spectra were constructed from midwater trawl collections each year in April, June–August, and October 1995–2000. The size spectra were bimodal: the first spectral dome corresponded to small zooplanktivorous fish, primarily bay anchovyAnchoa mitchilli; the second dome consisted of larger fish from several feeding guilds that are supported by multiple prey-predator linkages. Annual production estimates of pelagic fish, derived from a mean production to biomass ratio, varied nearly three-fold, ranging from 162 × 109 kcal (125 × 103 tons) in 1996 to 457 × 109 kcal (352 × 103 tons) in 2000. Seasonally, the biomass level and mean individual sizes of fish in the first dome increased from April to October, while the biomass level of the second dome was relatively stable. Regionally, biomass levels in the second dome were higher than biomasses in the first dome for the upper and lower Bay, but were minimal in the middle Bay where seasonal and episodic hypoxia occurs. To test a benthic-pelagic coupling hypothesis that could explain the higher biomass in the second domes for the lower and upper Bay, a cyclic size-spectrum model was fit that included only species in the zooplanktivorous-piscivorous fish guilds. The mean, normalized slope equaled ?1, indicating that zooplanktivorous fish may support piscivore production, but that a benthic-pelagic linkage is required to fully support fish production in the second dome. Interannual variability in slopes and intercepts of modeled size spectra was related to salinity, recruitment level of bay anchovy, and the primary axis of a correspondence analysis (salinity effect) on fish community structure. The spectral slope and intercept of normalized spectra were lowest in 1996, a near-record wet year. Results suggest that fish size spectra can be developed as useful indicators of ecosystem state and response to perturbations, especially if prey-predator relationships are explicitly represented.     

A pollution history of Chesapeake Bay

Present day anthropogenic fluxes of some heavy metals to central Chesapeake Bay appear to be intermediate to those of the southern California coastal region and those of Narragansett Bay. The natural fluxes, however, are in general higher. On the bases of Pb-210 and Pu-239 + 240 geochronologies and of the time changes in interstitial water compositions, there is a mixing of the upper 30 or so centimeters of the sediments in the mid-Chesapeake Bay area through bioturbation by burrowing mollusks and polychaetes. Coal, coke and charcoal levels reach one percent or more by dry weight in the deposits, primarily as a consequence of coal mining operations.     

Wind Modulation of Dissolved Oxygen in Chesapeake Bay

A numerical circulation model with a simplified dissolved oxygen module is used to examine the importance of wind-driven ventilation of hypoxic waters in Chesapeake Bay. The model demonstrates that the interaction between wind-driven lateral circulation and enhanced vertical mixing over shoal regions is the dominant mechanism for providing oxygen to hypoxic sub-pycnocline waters. The effectiveness of this mechanism is strongly influenced by the direction of the wind forcing. Winds from the south are most effective at supplying oxygen to hypoxic regions, and winds from the west are shown to be least effective. Simple numerical simulations demonstrate that the volume of hypoxia in the bay is nearly 2.5 times bigger when the mean wind is from the southwest as compared to the southeast. These results provide support for a recent analysis that suggests much of the long-term variability of hypoxia in Chesapeake Bay can be explained by variations in the summertime wind direction.     

Dissolved and particulate organic carbon in Chesapeake Bay

We measured dissolved and particulate organic carbon (DOC and POC) in samples collected along 13 transects of the salinity gradient of Chesapeake Bay. Riverine DOC and POC end-members averaged 232±19 μM and 151±53 μM, respectively, and coastal DOC and POC end-members averaged 172±19 μM and 43±6 μM, respectively. Within the chlorophyll maximum, POC accumulated to concentrations 50–150 μM above those expected from conservative mixing and it was significantly correlated with chlorophylla, indicating phytoplankton origin. POC accumulated primarily in bottom waters in spring, and primarily in surface waters in summer. Net DOC accumulation (60–120 μM) was observed within and downstream of the chlorophyll maximum, primarily during spring and summer in both surface and bottom waters, and it also appeared to be derived from phytoplankton. In the turbidity maximum, there were also net decreases in chlorophylla (?3 μg l^?1 to ?22 μg l^?1) and POC concentrations (?2 μM to ?89 μM) and transient DOC increases (9–88 μM), primarily in summer. These occurred as freshwater plankton blooms mixed with turbid, low salinity seawater, and we attribute the observed POC and DOC changes to lysis and sedimentation of freshwater plankton. DOC accumulation in both regions of Chesapeake Bay was estimated to be greater than atmospheric or terrestrial organic carbon inputs and was equivalent to ≈10% of estuarine primary production.     

Scales of nutrient-limited phytoplankton productivity in Chesapeake Bay

The scales on which phytoplankton biomass vary in response to variable nutrient inputs depend on the nutrient status of the plankton community and on the capacity of consumers to respond to increases in phytoplankton productivity. Overenrichment and associated declines in water quality occur when phytoplankton growth rate becomes nutrient-saturated, the production and consumption of phytoplankton biomass become uncoupled in time and space, and phytoplankton biomass becomes high and varies on scales longer than phytoplankton generation times. In Chesapeake Bay, phytoplankton growth rates appear to be limited by dissolved inorganic phosphorus (DIP) during spring when biomass reaches its annual maximum and by dissolved inorganic nitrogen (DIN) during summer when phytoplankton growth rates are highest. However, despite high inputs of DIN and dissolved silicate (DSi) relative to DIP (molar ratios of N∶P and Si∶P>100), seasonal accumulations of phytoplankton biomass within the salt-intruded-reach of the bay appear to be limited by riverine DIN supply while the magnitude of the spring diatom bloom is governed by DSi supply. Seasonal imbalances between biomass production and consumption lead to massive accumulations of phytoplankton biomass (often>1,000 mg Chl-a m^?2) during spring, to spring-summer oxygen depletion (summer bottom water <20% saturation), and to exceptionally high levels of annual phytoplankton production (>400 g m^?2 yr^?1). Nitrogen-dependent seasonal accumulations of phytoplankton biomass and annual production occur as a consequence of differences in the rates and pathways of nitrogen and phosphorus cycling within the bay and underscore the importance of controlling nitrogen inputs to the mesohaline and lower reaches of the bay.     

Wave-forced resuspension of upper Chesapeake Bay muds

Moored instruments were used to make observations of near bottom currents, waves, temperature, salinity, and turbidity at shallow (3.5 m and 5.5 m depth) dredged sediment disposal sites in upper Chesapeake Bay during the winters of 1990 and 1991 to investigate time-varying characteristics of resuspension processes over extended periods. Resulting time series data show the variability of two components of the suspended sediment concentration field. Background suspended sediment concentrations varied inversely with salinity and in direct relation to Susquehanna River flow. Muddy bottom sediments were also resuspended locally by both tidal currents and wind-wave forcing, resulting in short-term increases and decreases in suspended concentration, with higher peak concentrations near the bottom. In both years, episodes of wave-forced resuspension dominated tidal resuspension on an individual event basis, exceeding most tidal resuspension peaks by a factor of 3 to 5. The winds that generated the waves responsible for the observed resuspension events were not optimal for wave generation, however. Application of a simple wind-wave model showed that much greater wave-forced resuspension than that observed might be generated under the proper conditions. The consolidated sediments investigated in 1990 were less susceptible to both tidal and wave-forced resuspension than the recently deposited sediments investigated in 1991. There was also some indication that wave-forced resuspension increased erodibility of the bottom sediments on a short-term basis. Wave-forced resuspension is implicated as an important part of sediment transport processes in much of Chesapeake Bay. Its role in deeper, narrower, and more tidally energetic estuaries is not as clear, and should be investigated on a case-by-case basis.     

Observations of a Chesapeake Bay tidal front

This paper presents combined conductivity-temperature-depth (CTD), thermistor chain, current meter, and acoustic backscatter observations of a tidal front observed in the Chesapeake Bay. The data were obtained from a moored platform as the front migrated past the platform. The thermistor chain and CTD data show an interface that slopes steeply down from the surface to an asymptotic depth of 6 m, marking the bottom of the light-water pool. The thermistor chain data show much higher activity levels within the light-water pool as compared to the dense-water pool. Current meter data taken at 3 m show a pronounced shear in the currents upon crossing the frontal boundary. The acoustic backscatter from a layer of copepods positioned on the interface shows episodic occurrences of overturning at the interface. This observation is borne out by the concurrent thermistor chain data, which also show the overturning events.     

Phytoplankton reference communities for Chesapeake Bay and its tidal tributaries

Phytoplankton reference communities for Chesapeake Bay were quantified from least-impaired water quality conditions using commonly measured parameters and indicators derived from measured parameters. A binning approach was developed to classify water quality. Least-impaired conditions had relatively high water column transparency and low concentrations of dissolved inorganic nitrogen and orthophosphate. Reference communities in all seasons and salinity zones are characterized by consistently low values of chlorophylla and pheophytin coupled with relative stable proportions of the phytoplankton taxonomic groups and low biomasses of key bloom-forming species. Chlorophyll cell content was lower and less variable and average cell size and seasonal picophytoplankton biomass tended to be greater in the mesohaline and polyhaline reference communities as compared to the impaired communities. Biomass concentrations of the nano-micro phytoplankton size fractions (2–200 μm) in 12 of the 16 season-specific and salinity-specific reference communities were the same or higher than those in impaired habitat conditions, suggesting that nutrient reductions will not decrease the quantity of edible phytoplankton food available to large consumers. High (bloom) and low (bust) biomass events within the impaired phytoplankton communities showed strikingly different chlorophyll cell content and turnover rates. Freshwater flow had little effect on phytoplankton responses to water quality condition in most of the estuary. Improved water column transparency, or clarity, through the reduction of suspended sediments will be particularly important in attaining the reference communities. Significant nitrogen load reductions are also required.     

Regional geochemistry of trace elements in Chesapeake Bay sediments

The concentrations of Cr, Mn, Fe, Co, Ni, Cu, Zn, Cd, and Pb in 177 surface sediment samples from throughout Chesapeake Bay are reported. Analyses were made of both unfractionated samples and the <63 μm fractions. Analytical uncertainty, always less than ±10%, controlled reproducibility in analyses of the <63 μm fractions, but sampling variance controlled reproducibility in the unfractionated samples, especially when coarse-grained sediments were being analyzed. Sediments in the northernmost part of the bay are enriched relative to average continental crust in all elements except Cr. This reflects the composition of dissolved and suspended material being delivered to that region by the Susquehanna River. The enriched sediments appear not to be transported south of Baltimore in significant quantily. Zinc, cadmium, and lead are enriched relative to average crust throughout the bay and in most other estuaries in the eastern United States.     

Turbidity Maximum Entrapment of Phytoplankton in the Chesapeake Bay

Estuarine turbidity maxima (ETM) play an important role in zooplankton and larval fish productivity in many estuaries. Yet in many of these systems, little is known about the food web that supports this secondary production. To see if phytoplankton have the potential to be a component of the ETM food web in the Chesapeake Bay estuary a series of cruises were carried out to determine the biomass distribution and floral composition of phytoplankton in and around the ETM during the winter and spring using fluorometry, high-performance liquid chromatography (HPLC), and microscopy. Two distinct phytoplankton communities were observed along the salinity gradient. In lower salinity waters, biomass was low and the community was composed mostly of diatoms, while in more saline waters biomass was high and the community was composed mostly of mixotrophic dinoflagellates, which were often concentrated in a thin layer below the pycnocline. Phytoplankton biomass was always low in the ETM, but high concentrations of phytoplankton pigment degradation products and cellular remains were often observed suggesting that this was an area of high phytoplankton mortality and/or an area where phytoplankton derived particulate organic matter was being trapped. These results, along with a box model analysis, suggest that under certain hydrodynamic conditions phytoplankton derived organic matter can be trapped in ETM and potentially play a role in fueling secondary production.     

Hydrogen isotopes in dinosterol from the Chesapeake Bay estuary

The hydrogen isotope ratio of the dinoflagellate sterol dinosterol (4α,23,24-trimethyl-5α-cholest-22E-en-3β-ol) was measured in suspended particles and surface sediments from the Chesapeake Bay estuary in order to evaluate the influence of salinity on hydrogen isotope fractionation. D/H fractionation was found to decrease by 0.99 ± 0.23‰ per unit increase in salinity over the salinity range 10-29 PSU, a similar decrease to that observed in a variety of lipids from hypersaline ponds on Christmas Island (Kiribati). We hypothesize that the hydrogen isotopic response to salinity may result from diminished exchange of water between algal cells and their environment, lower growth rates and/or increased production of osmolytes at high salinities. Regardless of the mechanism, the consistent sign and magnitude of dinosterol δD response to changing salinity should permit qualitative to semi-quantitative reconstructions of past salinities from sedimentary dinosterol δD values.     

Long-Term Trends of Nutrients and Phytoplankton in Chesapeake Bay

Climate effects on hydrology impart high variability to water-quality properties, including nutrient loadings, concentrations, and phytoplankton biomass as chlorophyll-a (chl-a), in estuarine and coastal ecosystems. Resolving long-term trends of these properties requires that we distinguish climate effects from secular changes reflecting anthropogenic eutrophication. Here, we test the hypothesis that strong climatic contrasts leading to irregular dry and wet periods contribute significantly to interannual variability of mean annual values of water-quality properties using in situ data for Chesapeake Bay. Climate effects are quantified using annual freshwater discharge from the Susquehanna River together with a synoptic climatology for the Chesapeake Bay region based on predominant sea-level pressure patterns. Time series of water-quality properties are analyzed using historical (1945–1983) and recent (1984–2012) data for the bay adjusted for climate effects on hydrology. Contemporary monitoring by the Chesapeake Bay Program (CBP) provides data for a period since mid-1984 that is significantly impacted by anthropogenic eutrophication, while historical data back to 1945 serve as historical context for a period prior to severe impairments. The generalized additive model (GAM) and the generalized additive mixed model (GAMM) are developed for nutrient loadings and concentrations (total nitrogen—TN, nitrate?+?nitrate—NO₂?+?NO₃) at the Susquehanna River and water-quality properties in the bay proper, including dissolved nutrients (NO₂?+?NO₃, orthophosphate—PO₄), chl-a, diffuse light attenuation coefficient (K _D (PAR)), and chl-a/TN. Each statistical model consists of a sum of nonlinear functions to generate flow-adjusted time series and compute long-term trends accounting for climate effects on hydrology. We present results identifying successive periods of (1) eutrophication ca. 1945–1980 characterized by approximately doubled TN and NO₂?+?NO₃ loadings, leading to increased chl-a and associated ecosystem impairments, and (2) modest decreases of TN and NO₂?+?NO₃ loadings from 1981 to 2012, signaling a partial reversal of nutrient over-enrichment. Comparison of our findings with long-term trends of water-quality properties for a variety of estuarine and coastal ecosystems around the world reveals that trends for Chesapeake Bay are weaker than for other systems subject to strenuous management efforts, suggesting that more aggressive actions than those undertaken to date will be required to counter anthropogenic eutrophication of this valuable resource.     

Relationship between Hypoxia and Macrobenthic Production in Chesapeake Bay

Human development has degraded Chesapeake Bay's health, resulting in an increase in the extent and severity of hypoxia (≤2 mg O₂ l^-1). The Bay's hypoxic zones have an adverse effect on both community structure and secondary production of macrobenthos. From 1996 to 2004, the effect of hypoxia on macrobenthic production was assessed in Chesapeake Bay and its three main tributaries (Potomac, Rappahannock, and York Rivers). Each year, in the summer (late July???early September), 25 random samples of the benthic macrofauna were collected from each system, and macrobenthic production in the polyhaline and mesohaline regions was estimated using Edgar's allometric equation. Fluctuations in macrobenthic production were significantly correlated with dissolved oxygen. Macrobenthic production was 90 % lower during hypoxia relative to normoxia. As a result, there was a biomass loss of ~7,320–13,200 metric tons C over an area of 7,720 km², which is estimated to equate to a 20 % to 35 % displacement of the Bay's macrobenthic productivity during the summer. While higher consumers may benefit from easy access to stressed prey in some areas, the large spatial and temporal extent of seasonal hypoxia limits higher trophic level transfer, via the inhibition of macrobenthic production. Such a massive loss of macrobenthic production would be detrimental to the overall health of the Bay, as it comes at a time when epibenthic and demersal predators have high-energy demands.     

Climate Forcing and Salinity Variability in Chesapeake Bay,USA

Salinity is a critical factor in understanding and predicting physical and biogeochemical processes in the coastal ocean where it varies considerably in time and space. In this paper, we introduce a Chesapeake Bay community implementation of the Regional Ocean Modeling System (ChesROMS) and use it to investigate the interannual variability of salinity in Chesapeake Bay. The ChesROMS implementation was evaluated by quantitatively comparing the model solutions with the observed variations in the Bay for a 15-year period (1991 to 2005). Temperature fields were most consistently well predicted, with a correlation of 0.99 and a root mean square error (RMSE) of 1.5°C for the period, with modeled salinity following closely with a correlation of 0.94 and RMSE of 2.5. Variability of salinity anomalies from climatology based on modeled salinity was examined using empirical orthogonal function analysis, which indicates the salinity distribution in the Bay is principally driven by river forcing. Wind forcing and tidal mixing were also important factors in determining the salinity stratification in the water column, especially during low flow conditions. The fairly strong correlation between river discharge anomaly in this region and the Pacific Decadal Oscillation suggests that the long-term salinity variability in the Bay is affected by large-scale climate patterns. The detailed analyses of the role and importance of different forcing, including river runoff, atmospheric fluxes, and open ocean boundary conditions, are discussed in the context of the observed and modeled interannual variability.     

Discovery of the “Phantom” dinoflagellate in Chesapeake Bay

Since its discovery in natural estuarine habitat of North Carolina in 1991, the widespread impact of the toxic dinoflagellate, Pfiesteria piscicida (gen. et sp. nov.), popularly called the “phantom” dinoflagellate, on North Carolina fish stocks has been established, yet little is known about its influence outside of North Carolina estuaries. Here, we document the presence of P. piscicida in Chesapeake Bay. A fish kill was observed after inoculating an aquarium containing mummichogs with sediment samples from Jenkins Creek, a brackish creek (salinity 11‰) of the Chesapeake Bay system. P. piscicida was the cause of the kill, as supported by morphological, physiological, and histological evidence. The appearance and behavior of the algae and symptoms associated with fish mortality were consistent with those previously observed in P. piscicida-associated aquaria fish kills in North Carolina. The discovery of P. piscicida in Chesapeake Bay supports the speculation that these toxic dinoflagellates have a dramatic and far-reaching impact on fish stocks in shallow, eutrophic estuaries along the eastern United States.     

Can oyster restoration reverse cultural eutrophication in Chesapeake Bay?

We investigated the hypothesis that effects of cultural eutrophication can be reversed through natural resource restoration via addition of an oyster module to a predictive eutrophication model. We explored the potential effects of native oyster restoration on dissolved oxygen (DO), chlorophyll, light attenuation, and submerged aquatic vegetation (SAV) in eutrophic Chesapeake Bay. A tenfold increase in existing oyster biomass is projected to reduce system-wide summer surface chlorophyll by approximately 1 mg m⁻³, increase summer-average deep-water DO by 0.25 g m⁻³, add 2100 kg C (20%) to summer SAV biomass, and remove 30,000 kg d⁻¹ nitrogen through enhanced denitrification. The influence of osyter restoration on deep extensive pelagic waters is limited. Oyster restoration is recommended as a supplement to nutrient load reduction, not as a substitute.     

Hydrogen isotopes in individual alkenones from the Chesapeake Bay estuary

Hydrogen isotope ratios of individual alkenones from haptophyte algae were measured in suspended particles and surface sediment from the Chesapeake Bay (CB) estuary, eastern USA, in order to determine their relationship to water δD values and salinity. δD values of four alkenones (MeC_37:2, MeC_37:3, EtC_38:2, EtC_38:3) from particles and sediments were between −165‰ and −221‰ and increased linearly (R² = 0.7-0.9) with water δD values from the head to the mouth of the Bay. Individual alkenones were depleted in deuterium by 156-188‰ relative to water. The MeC₃₇ alkenones were consistently enriched by ∼12‰ relative to the EtC₃₈ alkenones, and the di-unsaturated alkenones of both varieties were consistently enriched by ∼20‰ relative to the tri-unsaturated alkenones. All of the increase in alkenone δD values could be accounted for by the water δD increase. Consequently, no net change in alkenone-water D/H fractionation occurred as a result of the salinity increase from 10 to 29. This observation is at odds with results from culture studies with alkenone-producing marine coccolithophorids, and from two field studies, one with a dinoflagellate sterol in the CB, and one with a wide variety of lipids in saline ponds on Christmas Island, that indicate a decline in D/H fractionation with increasing salinity. Why D/H fractionation in alkenones in the CB showed no dependence on salinity, while D/H fractionation in CB dinsoterol decreased by 1‰ per unit increase in salinity remains to be determined. Two hypotheses we consider to be valid are that (i) the assemblage of alkenone-producing haptophytes changes along the Bay and each species has a different sensitivity to salinity, such that no apparent trend in α_{alkenone-water} occurs along the salinity gradient, and (ii) greater osmoregulation capacity in coastal haptophytes may result in a diminished sensitivity of alkenone-water D/H fractionation to salinity changes.