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.     

Distribution and movement of shortnose sturgeon (Acipenser brevirostrum) in the Chesapeake Bay

During a reward program for Atlantic sturgeon (Acipenser oxyrinchus), 40 federally endangered shortnose sturgeon (Acipenser brevirostrum) were captured and reported by commercial fishers between January 1996 and January 2000 from the Chesapeake Bay. Since this is more than double the number of published records of shortnose sturgeon in the Chesapeake Bay between 1876 and 1995, little information has been available on distributions and movement. We used fishery dependent data collected during the reward program to determine the distribution of shortnose sturgeon in the Chesapeake Bay. Sonically-tagged shortnose sturgeon in the Chesapeake Bay and Delaware River were tracked to determine if individuals swim through the Chesapeake and Delaware Canal. Shortnose sturgeon were primarily distributed within the upper Chesapeake Bay. The movements of one individual, tagged within the Chesapeake Bay and later relocated in the canal and Delaware River, indicated that individuals traverse the Chesapeake and Delaware Canal.     

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.     

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.     

Spatial and temporal differences of Atlantic menhaden (Brevoortia tyrannus) recruitment across major drainages (1966–2004) of the Chesapeake Bay watershed

Atlantic menhaden (Brevoortia tyrannus) is well known for its commercial and ecological importance and has been historically declining in the Chesapeake Bay (Maryland), one of its principal nursery habitats along the eastern coast. Using data from the Striped Bass Seine Survey of the Maryland Department of Natural Resources (2003), we evaluated how the distribution of Atlantic menhaden has changed from 1966 to 2004 for 12 river drainages. We observed significant or marginally significant declines in 42% of the drainages, with drainages of the northern Bay showing the majority of those declines. Continued recruitment to several drainages of the Bay may partly explain why the adult spawning population is not declining. We determined if temporal changes in abundance were related to changes in salinity or water quality for five major drainages of the watershed. For one of these drainages, the Patuxent River, differences in productivity across sites largely explained differences in abundance. For the four remaining drainages, differences in recruitment could not be explained by productivity or salinity gradients. While reducing nitrogen loading and enhancing water clarity may improve Atlantic menhaden production, we suggest that the role of offshore processes on large-scale declines has been largely neglected and studies on larval ingression are necessary for further elucidation of spatial and temporal patterns of juvenile distribution in the Chesapeake Bay.     

Potential Salinity and Temperature Futures for the Chesapeake Bay Using a Statistical Downscaling Spatial Disaggregation Framework

Estuaries are productive and ecologically important ecosystems, incorporating environmental drivers from watersheds, rivers, and the coastal ocean. Climate change has potential to modify the physical properties of estuaries, with impacts on resident organisms. However, projections from general circulation models (GCMs) are generally too coarse to resolve important estuarine processes. Here, we statistically downscaled near-surface air temperature and precipitation projections to the scale of the Chesapeake Bay watershed and estuary. These variables were linked to Susquehanna River streamflow using a water balance model and finally to spatially resolved Chesapeake Bay surface temperature and salinity using statistical model trees. The low computational cost of this approach allowed rapid assessment of projected changes from four GCMs spanning a range of potential futures under a high CO₂ emission scenario, for four different downscaling methods. Choice of GCM contributed strongly to the spread in projections, but choice of downscaling method was also influential in the warmest models. Models projected a ~2–5.5 °C increase in surface water temperatures in the Chesapeake Bay by the end of the century. Projections of salinity were more uncertain and spatially complex. Models showing increases in winter-spring streamflow generated freshening in the Upper Bay and tributaries, while models with decreased streamflow produced salinity increases. Changes to the Chesapeake Bay environment have implications for fish and invertebrate habitats, as well as migration, spawning phenology, recruitment, and occurrence of pathogens. Our results underline a potentially expanded role of statistical downscaling to complement dynamical approaches in assessing climate change impacts in dynamically challenging estuaries.     

On the hydrographic distribution in the source region of the Delaware coastal current

The buoyant discharge from Delaware Bay forms two separate branches of residual outflow near the bay mouth, one along each shore. Upon exiting the bay, the branch along the Delaware shore turns right to form the southward flowing Delaware coastal current along the inner continental shelf off the Delaware, Maryland, and Virginia coasts. CTD and thermosalinograph, data collected at the mouth of Delaware Bay over two semidiurnal tidal cycles are used to examine the hydrographic distribution at the source region of the Delaware coastal current. In this region the buoyant source water of the coastal current, is largely detached from the shoreline and confined to the top 15 m of the water column over much of the tidal cycles. The core of the coastal current's source water, as defined by the point of salinity minimum, is located over the deep channel well offshore of the Delaware coast. The separation between this buoyant water and the more saline waters right along the Delaware coast and that in the central part of the bay mouth are marked by regions of high horizontal salinity gradients. The horizontal salinity gradients around the inshore and offshore boundaries of the source water of the coastal current are intensified during the flood tide, and clearly defined fronts (with a change of 3‰ over a distance of 150 m) are present at the offshore boundary near the end of the flood tide. The structure of the mean flow and the distribution of the brackish coastal current on the inner continental shelf contribute to the persistence of stratification in the source region off the Delaware shore throughout the ebb and flood tides. In contrast, the ebb-induced stratification in the region off the New Jersey shore is quickly destroyed with the onset of the flood current.     

The impact of anoxia on Mn fluxes in the Chesapeake Bay

Release of Mn into bottom waters of the Chesapeake Bay occurs in massive quantities during summertime periods of anoxia. The measured flux is much greater than can be accounted for by diffusion from pore waters under normal concentration gradients and is enhanced by dissolution of Mn-oxides under low Eh conditions. Extrapolating flux rates to the area of the Bay which goes anoxic during the summer indicates that this internal cycling accounts for four times as much Mn entering Bay waters as the total annual input (soluble and suspended) from the Susquehanna River.     

Dispersion and recruitment of crab larvae in the Chesapeake Bay plume: Physical and biological controls

As part of the Microbial Exchanges and Coupling in Coastal Atlantic Systems (MECCAS) Project, crab larvae were collected in the shelf waters off Chesapeake Bay in June and August 1985 and April 1986. We conducted hydrographic (temperature, salinity, nutrients) and biological (chlorophyll, copepods) mapping in conjunction with Eulerian and Lagrangian time studies of the vertical distribution of crab larvae in the Chesapeake Bay plume. These abundance estimates are used with current meter records and drifter trajectories to infer mechanisms of larval crab dispersion to the shelf waters and recruitment back into Chesapeake Bay. The highest numbers of crab larvae were usually associated with the Chesapeake Bay plume, suggesting that it was the dominant source of crab larvae to shelf waters. Patches of crab larvae also were found in the higher salinity shelf waters, and possibly were remnants of previous plume discharge events. The distribution of crab larvae in the shelf waters changed on 1–2 d time scales as a consequence of both variations in the discharge rate of the Chesapeake Bay plume and local wind-driven currents. Downwelling-favorable winds (NW) intensified the coastal jet and confined the plume and crab larvae along the coast. In April during a downwelling event (when northwesterly winds predominated), crab zoeae were transported southward along the coast at speeds that at times exceeded 168 km d⁻¹. During June and August the upwelling-favorable winds (S, SW) opposed the anticyclonic turn of the plume and, via Ekman circulation, forced the plume and crab larvae to spread seaward. Plume velocities during these conditions generally were less than 48 km d⁻¹. The recruitment of crab larvae to Chesapeake Bay is facilitated in late summer by the dominance of southerly winds, which can reverse the southward flow of shelf waters. Periodic downwelling-favorable winds can result in surface waters and crab larvae moving toward the entrance of Chesapeake Bay. Approximately 27% of the larval crabs spend at least part of the day in bottom waters, which have a residual drift toward the bay mouth. There appears to be a variety of physical transport mechanisms that can enhance the recruitment of crab larvae into Chesapeake Bay.     

Distribution ofVibrio parahaemolyticus in Chesapeake Bay during the summer season

The distribution ofVibrio parahaemolyticus in Chesapeake Bay during the warmer weather of the summer months was examined. This species was found throughout the Chesapeake Bay and its tributaries, even in areas of very low salinity. Counts of this species ranged from 0.04 per 100 ml to 46 per 100 ml in the water column and 2.03 to ≥2.4×10³ per 100 cc of sediment. A variety of physical, chemical and bacteriological properties associated with the incidence and distribution ofV. parahaemolyticus were examined and salinity was found to be the major influence among the factors examined. Correlation and regression analysis showed that the population size of this species increased with increasing salinity in the estuary.     

Expression of the Estuarine species minimum in littoral fish assemblages of the Lower Chesapeake Bay Tributaries

Species richness declines to a minimum (artenminimum) in the oligohaline reach of estuaries and other large bodies of brackish water. To date, observations of this feature in temperate estuaries have been largely restricted to benthic macroinvertebrates. Five years of seine data collected during the summers of 1990–1995 in the major tidal tributaries to the lower Chesapeake Bay were examined to see if this feature arose in estuarine fish assemblages. Estimates of numerical species richness (alpha diversity) and rates of species turnover between sites (beta diversity) were generated via rarefaction and detrended correspondence analysis. Two spatial attributes of the distribution of littoral fish species along salinity gradients in the tributaries of the lower Chesapeake Bay were revealed: (1) a species richness depression in salinities of 8–10% and (2) a peak in the rate of species turnover associated with the tidal freshwater interface (salinities of 0–2%). Expression of the minimum is influenced by the physical length of the salinity gradient and the interaction between a species’ salinity preferences and tendency to make long excursions from favorable habitats.     

Long-Term Trends in Submersed Aquatic Vegetation (SAV) in Chesapeake Bay, USA, Related to Water Quality

Chesapeake Bay supports a diverse assemblage of marine and freshwater species of submersed aquatic vegetation (SAV) whose broad distributions are generally constrained by salinity. An annual aerial SAV monitoring program and a bi-monthly to monthly water quality monitoring program have been conducted throughout Chesapeake Bay since 1984. We performed an analysis of SAV abundance and up to 22 environmental variables potentially influencing SAV growth and abundance (1984–2006). Historically, SAV abundance has changed dramatically in Chesapeake Bay, and since 1984, when SAV abundance was at historic low levels, SAV has exhibited complex changes including long-term (decadal) increases and decreases, as well as some large, single-year changes. Chesapeake Bay SAV was grouped into three broad-scale community-types based on salinity regime, each with their own distinct group of species, and detailed analyses were conducted on these three community-types as well as on seven distinct case-study areas spanning the three salinity regimes. Different trends in SAV abundance were evident in the different salinity regimes. SAV abundance has (a) continually increased in the low-salinity region; (b) increased initially in the medium-salinity region, followed by fluctuating abundances; and (c) increased initially in the high-salinity region, followed by a subsequent decline. In all areas, consistent negative correlations between measures of SAV abundance and nitrogen loads or concentrations suggest that meadows are responsive to changes in inputs of nitrogen. For smaller case-study areas, different trends in SAV abundance were also noted including correlations to water clarity in high-salinity case-study areas, but nitrogen was highly correlated in all areas. Current maximum SAV coverage for almost all areas remain below restoration targets, indicating that SAV abundance and associated ecosystem services are currently limited by continued poor water quality, and specifically high nutrient concentrations, within Chesapeake Bay. The nutrient reductions noted in some tributaries, which were highly correlated to increases in SAV abundance, suggest management activities have already contributed to SAV increases in some areas, but the strong negative correlation throughout the Chesapeake Bay between nitrogen and SAV abundance also suggests that further nutrient reductions will be necessary for SAV to attain or exceed restoration targets throughout the bay.     

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.     

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.     

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.     

210Pb balance in Long Island Sound

A material balance is constructed for excess ²¹⁰Pb (relative to ²²⁶Ra) as a test of the retentivity of Long Island Sound for a reactive heavy metal. Excess ²¹⁰Pb is supplied to Long Island Sound chiefly by direct atmospheric deposition [1 ± 0.2(dis·min^?1)cm^?2·yr^?1]. Rivers supply less than 20% of the atmospheric flux, and other inputs, from open ocean waters, ²²⁶Ra decay, groundwater seepage, and sewage discharge, appear to be negligible. The total input of excess ²¹⁰Pb represents approximately the flux required to maintain the inventory of excess ²¹⁰Pb measured in sediment cores from central Long Island Sound; that is, excess ²¹⁰Pb is lost from Long Island Sound chiefly by radioactive decay. The retention of excess ²¹⁰Pb within Long Island Sound is achieved in two steps: a rapid removal of soluble ²¹⁰Pb onto suspended particles and the ongoing entrapment of particles in the basin by the residual bottom-water influx from the east.     

On the spatial structure of currents across the Chesapeake and Delaware Canal

The Chesapeake and Delaware (C&;D) Canal is a man-made waterway connecting two of the largest estuaries on the east coast of the United States: Chesapeake Bay and Delaware Bay. A set of current meter data collected during April–May 1975 along two cross-sections of the C&;D Canal was used to examine the spatial distributions of the currents at tidal and subtidal time scales. The different responses of the Chesapeake and Delaware Bays to tidal and wind forcing produce significant differences in sea level fluctuations between the two ends of the canal. These alongcanal surface slopes produce significant barotropic current fluctuations at both tidal (semidiurnal and diurnal) and subtidal (2-d to 3-d) time scales. Under the influence of bottom friction, the barotropic currents near the surface are stronger than those at depth, but these currents do not exhibit significant lateral variations across the canal. On the other hand, the long-term flow in the canal exhibits strong lateral variability with eastward flow off the south shore of the canal and westward flow off the north shore of the canal. The lateral structure of the long-term flow may carry significant implications for the long-term exchange of material between the two bays.     

Estimation of net physical transport and hydraulic residence times for a coastal plain estuary using box models

A box model based on salinity distributions and freshwater inflow measurements was developed and used to estimate net non-tidal physical circulation and hydraulic residence times for Patuxent River estuary, Maryland, a tributary estuary of Chesapeake Bay. The box model relaxes the usual assumption that salinity is at steady-state, an important improvement over previous box model studies, yet it remains simple enough to have broad appeal. Average monthly 2-dimensional net non-tidal circulation and residence times for 1986–1995 are estimated and related to river flow and salt water inflow as estimated by the box model. An important result is that advective exchange at the estuary mouth was not correlated with Patuxent River flow, most likely due to effects of offshore salinity changes in Chesapeake Bay. The median residence time for freshwater entering at the head of the estuary was 68 d and decreased hyperbolically with increasing river flow to 30 d during high flow. Estimates of residence times for down-estuary points of origin showed that, from the head of the estuary to its mouth, control of flushing changed from primarily river flow to other factors regulating the intensity of gravitational circulation.     

Estimating the mean temperature and salinity of the Chesapeake Bay mouth

We used an extensive temperature and salinity data set to develop a statistically meaningful way of estimating mean temperature and salinity from discrete measurements in the mouth of Chesapeake Bay. From April 1992 to December 1998, the Center for Coastal Physical Oceanography completed 73 monthly hydrographic sections at high spring tide across the mouth of Chesapeake Bay. Time series of area weighted mean bay mouth temperature (MBMT) and salinity (MBMS) were calculated. We found that at any time the temperature at any location in the section correlated with the MBMT with a r² of 0.95 or better. A similar analysis for salinity showed that the best correlation was about 0.9 with many locations below 0.8. A correlation between MBMT and temperature at a nearby tide station indicated it was possible to estimate MBMT from the temperature at the tide station to ±0.74°C (90% confidence interval). Salinity was not measured at the tide station, but the correlation at a location in the section similar to the tide station indicates that MBMS can be estimated with an error of ±1.5 (90% confidence interval).