Transient hydrodynamic and salinity simulations in the Chesapeake Bay network

A transient network model is applied to the Chesapeake Bay and its tributary estuaries. Calibration of the model is based on only three external parameters: a friction factor that is spatially described, and two global constants required to calibrate a dynamic dispersion relationship that depends on both the local salinity gradient and hydraulic conditions. The transient hydrodynamics and the transient salinity distribution of the Bay and its tributary estuaries are simulated for the period of one month and comparisons made between calculated and observed salinities.     

Preliminary observation of a tidal intrusion front inside the mouth of the Chesapeake Bay

In situ sampling during flooding tidal flow and southerly winds shows the up-estuary translation of a surface front along the Thimble Shoal channel, which is located in the southern part of Chesapeake Bay. Currents and surface density were measured using a towed acoustic Doppler current profiler and CTD, and a ship-borne radar was used to monitor the orientation and planform of the front, which varied over time. These preliminary observations suggest that dense shelf water, upwelled to the surface along the coast near Cape Henry, Virginia, can be found well into the southern part of the bay during flood, and that the boundary between the intruding shelf water and less dense estuarine water has many of the same characteristics as tidal intrusion fronts found in smaller scale estuaries.     

Meteorological influences on sea level and water temperature in the lower Chesapeake Bay: 1992

The influence of atmospheric forcing on the flow and heat transports in the lower Chesapeake Bay and the adjacent coastal ocean were studied by comparing nontidal sea level and sea surface temperature variations in this region with meteorological data for 1992. Northeasterly and southwesterly winds caused the greatest changes in mean sea level (greater than 0.25 m) throughout the year. Northeastely winds caused a more rapid response than southwesterly winds, causing sea-level rises in less than 6 h. Barometric pressure changes typically contributed approximately 10% to extreme sea-level variations and were less influential than wind stress in most cases. Wind forcing was also responsible for summer events in which the horizontal water temperature gradient between two near-surface locations in the vicinity of the bay mouth vanished. These zero-gradient events corresponded to inflows and outflows at the bay's entrance caused by northeasterly and southwesterly winds, respectively. Wind-induced advection outside the lower Chesapeake Bay was additionally responsible for extreme heat flux variations. Heat gains and losses during the spring and fall occurred in pulsating events related to wind direction but were probably not connected to lower bay processes.     

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.     

The occurrence of ladyfish,Elops saurus,larvae in low salinity waters and another record for Chesapeake Bay

Some confusion exists concerning the early life history of ladyfish,Elops saurus, as a result of the use of confusing terms in describing the three morphologically and ecologically distinct growth phases of this species. This has resulted in conflicting data on the occurrence of stage I (leptocephali) and stage II (early metamorphic) larvae in oligohaline and mesohaline estuarine zones. The early life history ofE. saurus is reviewed relative to the collection of 17 stage II larvae at the freshwater transition zone of the James River, Virginia. There are no known spawning populations ofE. saurus north of Cape Hatteras, North Carolina. These collections represent the first recurrent record north of Cape Hatteras and for Chesapeake 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.     

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.     

Delaware Bay and Chesapeake Bay populations of the horseshoe crabLimulus polyphemus are genetically distinct

The Delaware Bay contains the world’s largest population of horseshoe crabs, which constitute an ecologically significant component of this estuarine ecosystem. The North Atlantic speciesLimulus polyphemus has an extensive geographical distribution, ranging from New England to the Gulf of Mexico. Recent assessments of the Delaware Bay population based on beach spawning and trawling data have suggested a considerable decrease in the number of adult animals since 1990. Considerable debate has centered on the accuracy of these estimates and their impact on marine fisheries management planning. Compounding this problem is the lack of information concerning the genetic structure of Atlantic horseshoe crab populations. This study assessed patterns of genetic variation within and between the horseshoe crab populations of Delaware Bay and Chesapeake Bay, using both Random Amplification of Polymorphic DNA (RAPD) and DNA sequence analysis of the mitochondrial cytochrome oxidase I gene (COI). We examined 41 animals from Delaware Bay and 14 animals from the eastern shore of Chesapeake Bay. To provide high quality, uncontaminated genomic DNA for RAPD analysis, DNA was isolated from hemocytes by direct cardiac puncture, purified by spin column chromatography, and quantified by agarose gel electrophoresis. RAPD fingerprints revealed a relative paucity of polymorphic fragments, with generally homogeneous banding patterns both within and between populations. DNA sequence analysis of 515 bases of the 5′ portion of the mitochondrial COI gene showed haplotype diversity in the Chesapeake Bay sample to be significantly higher than in the Delaware Bay sample, despite the larger size of the latter. Haplotype analysis indicates minimal contemporary gene flow between Delaware Bay and Chesapeake Bay crab populations, and further suggests that the Delaware Bay population is recovering from a recent population decline.     

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.     

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.     

A sediment chronology of the eutrophication of Chesapeake Bay

Chesapeake Bay sediments were examined for biogeochemical evidence of eutrophication trends using two mesohaline sediment cores. Measurements of ²¹⁰Pb geochronology and sediment profiles of organic carbon, nitrogen, organic phosphorus, inorganic phosphorus, and biogenis silica (BSi) were used used to develop temporal concentration trends. Recent sediments have 2–3 times as much organic carbon and nitrogen as sediments from 80 to 100 yr ago, but the increases result from both changes in organic matter deposition and time-dependent changes in organic matter decomposition rates. Despite increases in phosphorus loading, no major changes in phosphorus concentration were noted throughout most of the century; anthropogenic phosphorus deposition, though not evident in sulfidic mid-bay sediments, must occur in more oxidizing sediment environments in both the northern and southern bays. Temporal trends in BSi concentrations are much less evident and the lack of substantial increases in this century suggest that BSi inputs may be capped by late spring-summer Si limitation.     

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.     

Environmental Models and Public Stakeholders in the Chesapeake Bay Watershed

The Chesapeake Bay is not only North America's largest estuary, but it is also home to one of most comprehensive computational modeling efforts that seeks to study and integrate a very large range of the complex socio-ecological dynamics within its watershed. Known as the Chesapeake Bay Modeling System (CBMS), this suite of models are invaluable to scientists, environmental policymakers, resource managers, and local government and community leaders, all of whom are engaged in efforts to reduce the negative impacts of population growth and development in the watershed on the ecosystems and living resources of the Chesapeake Bay. Until recently, the results of the CBMS were used to guide voluntary efforts to reduce nutrient and sediment runoff into the Chesapeake Bay. However, the results are now being used to guide and evaluate the implementation of Watershed Implementation Plans at the state and county levels, which in turn are based on CBMS estimates of total maximum daily loads (TMDLs) for these sub-watershed regions. The use of the CBMS for these regulatory efforts has also increased the number of public stakeholders who are now able to directly inquire about the how the CBMS works and raise questions about its ability to accurately represent their land use decisions and practices. This "going public" of the CBMS has raised many societal and cultural issues, including the articulation of local expert and scientific knowledge, and modelers, scientists, policymakers, and resource managers are now realizing the need to understand more of the human dimensions arising from the translation and implementation of the CBMS. A multidisciplinary approach is needed to understand these rapidly emerging societal and cultural dimensions of the expanded use of the CBMS. In this paper, we draw upon our collective experience as social and natural scientists, with modeling experience, to describe a range of social, economic, political, and cultural issues that have emerged as a result of the CBMS being used to support mandatory nutrient reduction regulations (TMDL).     

Reconsidering the physics of the Chesapeake Bay estuarine turbidity maximum

A series of cruises was carried out in the estuarine turbidity maximum (ETM) region of Chesapeake Bay in 1996 to examine physical and biological variability and dynamics. A large flood event in late January shifted the salinity structure of the upper Bay towards that of a salt wedge, but most of the massive sediment load delivered by the Susquehanna River appeared to bypass the ETM zone. In contrast, suspended sediments delivered during a flood event in late October were trapped very efficiently in the ETM. The difference in sediment trapping appeared to be due to increases in particle settling speed from January to October, suggesting that the fate of sediments delivered during large events may depend on the season in which they occur. The ETM roughly tracked the limit of salt (defined as the intersection of the 1 psu isohaline with the bottom) throughout the year, but it was often separated significantly from the limit of salt with the direction of separation unrelated to the phase of the tide. This was due to a lag of ETM sediment resuspension and transport behind rapid meteorologically induced or river flow induced motion of the salt limit. Examination of detailed time series of salt, suspended sediment, and velocity collected near the limit of salt, combined with other indications, led to the conclusion that the convergence of the estuarine circulation at the limit of salt is not the primary mechanism of particle trapping in the Chesapeake Bay ETM. This convergence and its associated salinity structure contribute to strong tidal asymmetries in sediment resuspension and transport that collect and maintain a resuspendable pool of rapidly settling particles near the salt limit. Without tidal resuspension and transport, the ETM would either not exist or be greatly weakened. In spite of this repeated resuspension, sedimentation is the ultimate fate of most terrigenous material delivered to the Chesapeake Bay ETM. Sedimentation rates in the ETM channel are at least an order of magnitude greater than on the adjacent shoals, probably due to focusing mechanisms that are poorly understood.     

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.     

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.     

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.     

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.     

Influence of salinity on atrazine toxicity to a Chesapeake Bay copepod (Eurytemora affinis) and fish (Cyprinodon variegatus)

The objective of this study was to determine the influence of a range of salinities (5%., 15%., and 25%.) on the acute toxicity of atrazine to nauplii of the copepodEurytemora affinis and larvae of the sheepshead minnow,Cyprinodon variegatus. Ninety-six-hour LC50 values for the copepod were 0.5 mg 1^?1, 2.6 mg 1^?1, and 13.2 mg 1^?1 at salimities of 5%., 15%. and 25%. respectively. A comparison of LC50 values between adjacent salinities showed a statistical difference between 15%. and 25%. but not between 5%. and 15%.. Atrazine was more toxic toE. affinis at the lowest salinity. The 96-h LC50s for the sheepshead minnow were 16.2 mg 1^?1, 2.3 mg 1^?1, and 2.0 mg 1^?1 at salinities of 5%., 15%., and 25%., respectively. There was a statistical difference between LC50 values at 5%. and 15%. but not between 15%. and 25%.. In contrast toE. affinis results, atrazine was more toxic toC. variegatus at the highest salinity. The toxicity data from these species suggest that development of estuarine water quality criteria is warranted.     

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.