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.     

Long-term pCO2 dynamics in rivers in the Chesapeake Bay watershed

Streams and rivers are major exporters of C and other dissolved materials from watersheds to coastal waters. In streams and rivers, substantial amounts of terrigenous organic C is metabolized and degassed as CO₂ to the atmosphere. A long-term evaluation of CO₂ dynamics in streams is essential for understanding factors controlling CO₂ dynamics in streams in response to changes in climate and land-use. Long-term changes in the partial pressure of CO₂ (pCO₂) were computed in the Anacostia River and the lower Potomac River in the Chesapeake Bay watershed. Long-term estimates were made using routine monitoring data of pH, total alkalinity, and dissolved nutrients from 1985 to 2006 at 14 stations. Longitudinal variability in pCO₂ dynamics was also investigated along these rivers downstream of the urban Washington D.C. metropolitan area. Both rivers were supersaturated with CO₂ with respect to atmospheric CO₂ levels (392 μatm) and the highly urbanized Anacostia waters (202–9694 μatm) were more supersaturated than the Potomac waters (557–3800 μatm). Long-term variability in pCO₂ values may be due to changes in river metabolism and organic matter and nutrient loadings. Both rivers exchange significant amounts of CO₂ with the atmosphere (i.e., Anacostia at 0.2–72 mmol m⁻² d⁻¹ and Potomac at 0.12–24 mmol m⁻² d⁻¹), implying that waterways receiving organic matter and nutrient subsidies from urbanized landscapes have the potential to increase river metabolism and atmospheric CO₂ fluxes along the freshwater–estuarine continuum.     

Matrix-based modeling of blue crab population dynamics with applications to the Chesapeake Bay

Blue crab (Callinectes sapidus) plays an important ecological and economic role in estuaries from South America to New England. It supports a large commercial fishery in the United States with approximately one third of the landings taken from Chesapeake Bay. I developed a stage-based matrix model of the blue crab population to address three key questions: What is the ability of blue crab populations to support sustainable exploitations? What stages of the life cycle are most important in regulating the dynamics of the populations? And specific of the Chesapeake Bay, what is the importance of a winter dredge fishery in determining long-term sustainability of the crab population? The model indicated that with the current pattern of exploitation blue crab populations are able to sustain a total instantaneous mortality rate (Z)～0.7. If the natural mortality rate is estimated for a maximum life expectancy of 8 yr, this translates to moderate levels of exploitation (F<0.32). This value is less than the current estimate of exploitation in Chesapeake Bay (0.9–1.1) indicating that the level of exploitation in this system needs to be reduced to avoid overfishing. Transitions to and from small age-1 crabs were shown to be important in regulating the overall dynamics of the population. The egg production realized by large adults was also shown to be an important regulatory process. The model indicated that reductions in the winter dredge fishery would have a substantial role in ensuring the long-term sustainability of the population. Reductions in other sectors of the fishery are also required to ensure sustainability.     

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.     

Temporal and spatial dynamics of urea uptake and regeneration rates and concentrations in Chesapeake Bay

We examined the temporal and spatial variability of urea concentrations and urea uptake and regeneration rates collected on cruises along the longitudinal axis of the Chesapeake Bay between 1972 and 1998. Interannually, mean Bay-wide surface urea concentrations ranged between 0.49 and 0.91 μg-at N l^?1 with a nearly 50% decrease in surface concentrations observed between 1988 and 1998. Concentrations of urea from samples collected within ～1 m of the bottom were generally higher and much more varable than surface samples. Seasonally, two different patterns were observed in mean Bay-wide surface urea concentrations. Urea concentrations from near surface waters exhibited a clear summer peak for 1988 through 1994, while for 1973 and 1996 to 1998 a distinct winter-spring peak in concentration was observed. Urea concentrations from deeper waters showed a similar seasonal trend each year with peak concentrations measured in spring. Spatially, urea concentrations in the surface waters decreased in a conservative-type pattern from 0.91 μg-at N I^?1 at the freshwater end member to 0.46 μg-at N I^?1 at the ocean end member. Mean Bay-wide surface urea uptake rates displayed a seasonal pattern throughout the data set with maximum uptake rates (up to 0.33 μg-at N I^?1 h^?1) consistently observed during summer. Mean Bay-wide surface regeneration rates were highest but most variable during fall (1.63±0.82 μg-at N I^?1 h^?1). Mean urea uptake and regeneration rates displayed opposing spatial trends along the axis of the Bay with uptake rates being lowest in the North Bay where regeneration rates were highest. The average temporal and spatial patterns of urea concentration in Chesapeake Bay appear to reflect a balance between external inputs and internal biological recycling.     

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.     

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.     

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.     

Chesapeake Bay nutrient and plankton dynamics: III. The annual cycle of dissolved silicon

Silicic acid (H₄SiO₄) flux from the sediment, H₄SiO₄ concentration and river flow were used to obtain an annual dissolved silicon budget for Chesapeake Bay. H₄SiO₄ concentrations vary seasonally in the estuary: for a 12-year period, mean H₄SiO₄ concentrations in the mesohaline region were high both in spring and in late summer to early fall, and were low in late spring—occasionally approaching levels potentially limiting to diatom growth. Most of the annual allochthonous H₄SiO₄ supply to the estuary derives from the three major rivers, but regenerative H₄SiO₄ flux from the sediment to the water column exceeds the total riverine input by a factor of at least five. Sediment H₄SiO₄ efflux exhibits seasonality and averages approximately 2–3 mol Si m^?2 yr^?1. The high rates of sediment dissolution and efflux appear to maintain high levels of H₄SiO₄ in the mesohaline region, and Si-limitation of diatom growth there seems unlikely. The relative rates of biogenic silica formation and dissolution do not vary synchronously: seasonal variations in diatom productivity, sedimentary release of H₄SiO₄ and river flow all contribute to the observed late winter and late summer seasonal maxima and late spring minimum in water column H₄SiO₄ concentrations. If the only source of Si to support sedimentary H₄SiO₄ efflux is biogenic particulate silica recently deposited from the water column and this silica in turn was produced by diatoms in a ratio of 8C:1 Si, the minimum annual primary production by diatoms is at least 260 g C m^?2, approximately half of annual total plankton primary production. This estimate would be revised upwards according to the amount of particulate biogenic silica dissolving in the water column. Burial of biogenic silica amounts to from 2 to 84% of the sediment efflux of H₄SiO₄, depending on location in the bay. On an annual basis, burial represents from 60 to 100% of fluvial H₄SiO₄-Si inputs.     

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.     

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.     

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.     

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.     

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.     

Scientific Bases for Numerical Chlorophyll Criteria in Chesapeake Bay

In coastal ecosystems with long flushing times (weeks to months) relative to phytoplankton growth rates (hours to days), chlorophyll a (chl-a) integrates nutrient loading, making it a pivotal indicator with broad implications for ecosystem function and water-quality management. However, numerical chl-a criteria that capture the linkage between chl-a and ecosystem impairments associated with eutrophication (e.g., hypoxia, water clarity and loss of submerged aquatic vegetation, toxic algal blooms) have seldom been developed despite the vulnerability of these ecosystems to anthropogenic nutrient loading. Increases in fertilizer use, animal wastes, and population growth in the Chesapeake Bay watershed since World War II have led to increases in nutrient loading and chl-a. We describe the development of numerical chl-a criteria based on long-term research and monitoring of the bay. Baseline chl-a concentrations were derived using statistical models for historical data from the 1960s and 1970s, including terms to account for the effects of climate variability. This approach produced numerical chl-a criteria presented as geometric means and 90th percentile thresholds to be used as goals and compliance limits, respectively. We present scientific bases for these criteria that consider specific ecosystem impairments linked to increased chl-a, including low dissolved oxygen (DO), reduced water clarity, and toxic algal blooms. These multiple lines of evidence support numerical chl-a criteria consisting of seasonal mean chl-a across salinity zones ranging from 1.4 to 15 mg m^?3 as restoration goals and corresponding thresholds ranging from 4.3 to 45 mg m^?3 as compliance limits. Attainment of these goals and limits for chl-a is a precondition for attaining desired levels of DO, water clarity, and toxic phytoplankton prior to rapid human expansion in the watershed and associated increases of nutrient loading.     

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.     

The trophic consequences of oyster stock rehabilitation in Chesapeake Bay

There is mounting speculation that overharvesting of oyster stocks (Crassostrea virginica) in Chesapeake Bay may be a factor contributing to the decline in water quality and shifts in the dominance of species inhabiting the estuary. The trophic consequences of increasing the oyster population may be addressed using a simple quasi-equilibrium, mass action model of the exchanges transpiring in the Chesapeake mesohaline ecosystem. According to output from the model, increasing oyster abundance would decrease phytoplankton productivity as well as stocks of pelagic microbes, ctenophores, medusae, and particulate organic carbon. Recently acquired field data on phytoplankton productivity, bacterioplankton, and labile organic carbon in the vicinity of rafted oyster aquaculture support model predictions. The model also indicates that more oysters should increase benthic primary production, fish stocks, and mesozooplankton densities. Hence, augmenting the oyster community by restoring beds or introducing raft culture represents a potentially significant adjunct to the goal of mitigating eutrophication through curtailment of nutrient inputs. *** DIRECT SUPPORT *** A01BY059 00005     

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.