The Chester River, a tributary of Chesapeake Bay, provides critical habitats for numerous living species and oyster aquaculture, but faces increasing anthropogenic stresses due to excessive nutrient ...loading and hypoxia occurrence. An application of the Integrated Compartment Water Quality Model (ICM), coupled with the Finite-Volume Community Ocean Model (FVCOM), was carried out to study the controlling mechanisms and interannual variability in hypoxia occurrence from 2002 to 2011. Our study shows that hypoxia occurs mostly in the main stem in July, followed by August and June. On an interannual scale, 2005 had the highest hypoxia occurrence with an accumulative hypoxia volume of about 10 km3-days, whereas 2008 had the lowest occurrence with an accumulative hypoxia volume of about 1 km3-days. Nutrient loading is the predominant factor in determining the intensity and interannual variability in hypoxia in the Chester River estuary, followed by stratification and saltwater intrusion. Phosphorus has been found to be more efficient in controlling hypoxia occurrence than nitrogen due to their different limiting extent. On a local scale, the Chester River estuary is characterized by several meanders, and at certain curvatures helical circulation is formed due to centrifugal forces, leading to better reaeration and dissolved oxygen (DO) supply to the deeper layers. Our study provides valuable information for nutrient management and restoration efforts in the Chester River.
To protect the aquatic living resources of Chesapeake Bay, the Chesapeake Bay Program partnership has developed guidance for state water quality standards, which include ambient water quality ...criteria to protect designated uses (DUs), and associated assessment procedures for dissolved oxygen (DO), water clarity/underwater bay grasses, and chlorophyll-a. For measuring progress toward meeting the respective states' water quality standards, a multimetric attainment indicator approach was developed to estimate combined standards attainment. We applied this approach to three decades of monitoring data of DO, water clarity/underwater bay grasses, and chlorophyll-a data on annually updated moving 3-year periods to track the progress in all 92 management segments of tidal waters in Chesapeake Bay. In 2014–2016, 40% of tidal water segment-DU-criterion combinations in the Bay (n = 291) are estimated to meet thresholds for attainment of their water quality criteria. This index score marks the best 3-year status in the entire record. Since 1985–1987, the indicator has followed a nonlinear trajectory, consistent with impacts from extreme weather events and subsequent recoveries. Over the period of record (1985–2016), the indicator exhibited a positive and statistically significant trend (p < 0.05), indicating that the Bay has been recovering since 1985. Patterns of attainment of individual DUs are variable, but improvements in open water DO, deep channel DO, and water clarity/submerged aquatic vegetation have combined to drive the improvement in the Baywide indicator in 2014–2016 relative to its long-term median. Finally, the improvement in estimated Baywide attainment was statistically linked to the decline of total nitrogen, indicating responsiveness of attainment status to the reduction of nutrient load through various management actions since at least the 1980s.
•Chesapeake Bay's water quality history was assessed by using an indicator framework.•The indicator has a positive long-term trend (p < 0.05) and reached its peak in 2014–2016.•The indicator was responsive to extreme weather events but can recover afterwards.•Improvement of indicator score in 2014–2016 over its long-term average was driven by open water and deep channel dissolved oxygen.•The improvement in Baywide attainment was statistically linked to the decline of total nitrogen input.
Display omitted
The relative dominance of controlling factors in saltwater intrusion depends on the time scale in which the variability in characterized. Long-term time series data of salinity and temperature has ...been collected in the Chester River, Chesapeake Bay, showing considerable variabilities of saltwater intrusion with time scales from days to years. In this study, the Finite-Volume Community Ocean Model (FVCOM) was applied to the data over a decadal time scale from 2002 to 2011 to identify the drivers and their relative importance in controlling the variability of saltwater intrusion with different time scale ranging from events of a few days, through seasonal to interannual variations. FVCOM successfully reproduced the observed variations in saltwater intrusion over the entire simulation period, and thus completed time series data from monthly samples to daily and even at a higher temporal resolution for statistical analysis. The results shows that river discharge is the primary factor controlling the saltwater intrusion variability on interannual time scales, but sea surface level is more dominating on seasonal variations, while wind has more important influence on short-time scales during weather events, mostly through modulation of sea surface level on regional scales rather than local wind strain and Ekman transport. Tide has limited influence due to short tidal excursion and amplitude in the Chester Estuary. The estuarine geometry has a significant influence on the two-layer estuarine circulation and saltwater intrusion in the Chester Estuary. The upper estuary is the shallow narrow channel where two-layer estuarine circulation is limited, the meanders in the mid-estuary also restrict the development of the two-layer circulation due to helical circulation, while the lower estuary has stronger stratification and two-layer circulation as compared to the upper and mid-estuary.
Display omitted
•The relative dominance of controlling factors in saltwater intrusion changes over time scale.•The extension of the saltwater intrusion is a seasonal dispersive phenomenon.•Wind remotely controls saltwater intrusion through regional sea surface level.•Meanders restrict saltwater intrusion and gravitational circulation.
Eutrophication and hypoxia represent an ever-growing stressor to estuaries and coastal ecosystems due to population growth and climate change. Understanding water quality dynamics in shallow water ...systems is particularly challenging due to the complex physical and biogeochemical dynamics and interactions among them. Within shallow waters, benthic microalgae can significantly contribute to autotrophic primary production, generate organic matter, increase dissolved oxygen consumption, and alter nutrient fluxes at the sediment–water interface, yet they have received little attention in modeling applications. A state-of-the-art modeling system, the Semi-Implicit Cross-Scale Hydroscience Integrated System Model (SCHISM), coupled with the Integrated Compartment Model (ICM) of water quality and benthic microalgae, has been implemented in the Corsica River estuary, a tributary to Chesapeake Bay, to study benthic microalgal impact on water quality in shallow water systems. The model simulation has revealed a broad impact of benthic microalgae, ranging from sediment–water interface fluxes to water column dynamics, and the effects are observed from near-field to far-field monitoring stations. High-frequency variability and non-linearity dominate benthic microalgal dynamics, sediment oxygen demand, and nutrient fluxes at the sediment–water interface. Resource competition and supply determine the spatial scope of benthic microalgal impacts on far-field stations and the whole estuary system. Our study shows that benthic microalgae are a significant factor in shallow water dynamics that needs adequate attention in future observation and modeling applications.
Understanding shallow water biogeochemical dynamics is a challenge in coastal regions, due to the presence of highly variable land-water interface fluxes, tight coupling with sediment processes, ...tidal dynamics, and diurnal variability in biogeochemical processes. While the deployment of continuous monitoring devices has improved our understanding of high-frequency (12 - 24 hours) variability and spatial heterogeneity in shallow regions, mechanistic modeling of these dynamics has lagged behind conceptual and empirical models. The inherent complexity of shallow water systems is represented in the Corsica River estuary, a small basin within the Chesapeake Bay ecosystem, where abundant monitoring data have been collected from long-term monitoring stations, continuous monitoring sensors, synoptic sensor surveys, and measurements of sediment-water fluxes. A state-of-the-art modeling system, the Semi-implicit Cross-scale Hydroscience Integrated System Model (SCHISM), was applied to the Corsica domain with a high-resolution grid and nutrient loads from the most recent version of the Chesapeake Bay watershed model. The Corsica SCHISM model reproduced observed high-frequency variability in dissolved oxygen, as well as seasonal variability in chlorophyll-a and sediment-water fluxes. Time-series signal analyses using Empirical Model Decomposition and spectral analysis revealed that the diurnal and M2 tide frequencies are the dominant high-frequency modes and physical transport contributes a larger share to dissolved oxygen budgets than biogeochemical processes on an hourly time scale. Heterogeneity and patchiness in dissolved oxygen resulting from phytoplankton distributions and geometry-driven eddies amplify the physical transport effect, and on longer time scales oxygen is controlled more by photosynthesis and respiration. Our simulation demonstrates that interactions among physical and biological dynamics generate complex high-frequency variability in water quality and non-linear reposes to nutrient loading and environmental forcing in shallow water systems.
Low dissolved oxygen (DO) conditions are a recurring issue in waters of Chesapeake Bay, with detrimental effects on aquatic living resources. The Chesapeake Bay Program partnership has developed ...criteria guidance supporting the definition of state water quality standards and associated assessment procedures for DO and other parameters, which provides a binary classification of attainment or impairment. Evaluating time series of these two outcomes alone, however, provides limited information on water quality change over time or space. Here we introduce an extension of the existing Chesapeake Bay water quality criterion assessment framework to quantify the amount of impairment shown by space-time exceedance of DO criterion (“attainment deficit”) for a specific tidal management unit (i.e., segment). We demonstrate the usefulness of this extended framework by applying it to Bay segments for each three-year assessment period between 1985 and 2016. In general, the attainment deficit for the most recent period assessed (i.e., 2014-2016) is considerably worse for deep channel (DC; n = 10) segments than open water (OW; n = 92) and deep water (DW; n = 18) segments. Most subgroups -- classified by designated uses, salinity zones, or tributary systems -- show better (or similar) attainment status in 2014-2016 than their initial status (1985-1987). Some significant temporal trends (p < 0.1) were detected, presenting evidence on the recovery for portions of Chesapeake Bay with respect to DO criterion attainment. Significant, improving trends were observed in seven OW segments, four DW segments, and one DC segment over the 30 three-year assessment periods (1985-2016). Likewise, significant, improving trends were observed in 15 OW, five DW, and four DC segments over the recent 15 assessment periods (2000-2016). Subgroups showed mixed trends, with the Patuxent, Nanticoke, and Choptank Rivers experiencing significant, improving short-term (2000-2016) trends while Elizabeth experiencing a significant, degrading short-term trend. The general lack of significantly improving trends across the Bay suggests that further actions will be necessary to achieve full attainment of water quality standards. Insights revealed in this work are critical for understanding the dynamics of the Bay ecosystem and for further assessing the effectiveness of management initiatives aimed toward Bay restoration.
The 2010 Chesapeake Bay Total Maximum Daily Load was established for the water quality and ecological restoration of the Chesapeake Bay. In 2017, the latest science, data, and modeling tools were ...used to develop revised Watershed Implementation Plans (WIPs). In this article, we examine the vulnerability of the Chesapeake Bay watershed to the combined pressures of climate change and growth in population, agricultural intensity, and economic activity for the 60‐year period 1995–2055. The results will be used to revise WIPs, as needed, to account for expected increases in loads. Assessing changes relative to 1995 for the years 2025, 2035, 2045, and 2055, mean annual precipitation increases of 3.11%, 4.21%, 5.34%, and 6.91%, respectively, air temperature increases of 1.12, 1.45, 1.84, and 2.12°C, respectively, and potential evapotranspiration increases of 3.36%, 4.43%, 5.54%, and 6.35%, respectively, are projected. Population in the watershed is expected to grow by 3.5 million between 2025 and 2055. Watershed model results show incremental increases in streamflow (2.3%–6.2%), nitrogen (2.6%–10.8%), phosphorus (4.5%–26.7%), and sediment (3.8%–18.8%) loads to the tidal Bay due to climate change. Growth in population, agricultural intensity, development, and economic activity resulted in relatively smaller increases in loads compared to climate change.
•CART models were developed to predict nutrient limitation in Chesapeake Bay.•The approach satisfactorily reproduced bioassay-based nutrient limitation patterns.•The Bay showed more space of ...nitrogen-limitation in 2007–2017 than 1992–2002.•Nutrient limitation patterns in the Bay appear to vary with hydrologic conditions.•Further nutrient reductions are needed to achieve a less nutrient-saturated system.
Understanding the temporal and spatial roles of nutrient limitation on phytoplankton growth is necessary for developing successful management strategies. Chesapeake Bay has well-documented seasonal and spatial variations in nutrient limitation, but it remains unknown whether these patterns of nutrient limitation have changed in response to nutrient management efforts. We analyzed historical data from nutrient bioassay experiments (1992–2002) and data from long-term, fixed-site water-quality monitoring program (1990–2017) to develop empirical approaches for predicting nutrient limitation in the surface waters of the mainstem Bay. Results from classification and regression trees (CART) matched the seasonal and spatial patterns of bioassay-based nutrient limitation in the 1992–2002 period much better than two simpler, non-statistical approaches. An ensemble approach of three selected CART models satisfactorily reproduced the bioassay-based results (classification rate = 99%). This empirical approach can be used to characterize nutrient limitation from long-term water-quality monitoring data on much broader geographic and temporal scales than would be feasible using bioassays, providing a new tool for informing water-quality management. Results from our application of the approach to 21 tidal monitoring stations for the period of 2007–2017 showed modest changes in nutrient limitation patterns, with expanded areas of nitrogen-limitation and contracted areas of nutrient saturation (i.e., not limited by nitrogen or phosphorus). These changes imply that long-term reductions in nitrogen load have led to expanded areas with nutrient-limited phytoplankton growth in the Bay, reflecting long-term water-quality improvements in the context of nutrient enrichment. However, nutrient limitation patterns remain unchanged in the majority of the mainstem, suggesting that nutrient loads should be further reduced to achieve a less nutrient-saturated ecosystem.
Display omitted
In 2020, the Chesapeake Bay Program moved to offset impacts from climate change for the 30‐year period from 1995 through 2025 by having its seven watershed jurisdictions (Delaware, Maryland, New ...York, Pennsylvania, Virginia, West Virginia, and the District of Columbia) apply additional nutrient pollutant reduction practices. The climate change assessment was performed with integrated models of the Chesapeake watershed, airshed, and estuary. Scenarios run for the years 2025, 2035, 2045, and 2055 estimated effects from the different future climatic conditions. This article presents the results of that assessment and is intended to provide a guide to assist other modeling practitioners in assessing climate change impacts in coastal watersheds. Major influences of climate change that were quantified include increases in precipitation volume, potential evapotranspiration, watershed nutrient loads, tidal water temperature, and sea level. Minor influences quantified in the climate change analysis include changes in nutrient speciation and increases in wet deposition of nitrogen, CO2, rainfall intensity, tidal wetland loss, up‐estuary salt intrusion, and phytoplankton biomass. To offset climate change impacts from 1995 to 2025 on water quality, the scenarios indicate an additional 2.3 million and 0.3 million kg of nitrogen and phosphorus per annum, respectively, will need to be reduced beyond what is called for in the Chesapeake Total Maximum Daily Load.
A coupled estuarine hydrodynamic model and water quality model were used to analyze differences in destratification and anoxia/hypoxia reduction by wind directions in the north-south oriented ...Chesapeake estuary, USA. The predominant cross-channel bathymetry in the Bay's anoxic center is asymmetric with a steeper and narrower shoal on the eastern shore than on the western shore, which modifies wind-induced circulation differently for two opposite wind directions. Model experiments of winds for 2-day at 8 m/s indicated that, for a stratified water over the aforementioned asymmetric bottom topography, the easterly wind caused greater destratification and hypoxia reduction than the westerly wind. This is a result of differential modulations on the two wind-induced cross-channel circulations by the asymmetric cross channel bathymetry. The downwelling along the gentle slope in the easterly wind was characterized with stronger baroclinicity than the downwelling along the steep slope (nearly perpendicular to surfaces of constant density) in the westerly wind. On the broad slope, there undergo greater contrasting density readjustments to the vorticity changes around the bottom boundary layer (BBL) during upslope and downslope motions. During the upslope condition, the flow in BBL tends to decelerate under adverse pressure gradient which leads to a stable condition in the outer layer; whereas, during the downslope condition, the BBL tends to accelerate under favourable pressure gradient, which leads to unstable condition in the outer layer of the large scale flow. Overall, the easterly wind caused greater anoxia reduction than the westerly wind during the entire wind period. A similar case was found for northerly versus southerly winds in the early stages of the wind period; modulated by the aforementioned bathymetry on the wind-induced cross-channel circulation, the northerly wind caused greater anoxia reduction than the southerly wind. However, as wind continues, the wind-induced along-channel circulation influences a larger area of greater hypoxia in the mid-Bay, by which the southerly wind causes a greater destratification and anoxia reduction than the northerly wind. This can supersede the greater destratification and anoxia reduction by the northerly wind under the bathymetry-affected cross-channel circulation.
Display omitted
•Asymmetric cross-channel bathymetry modifies wind-induced lateral circulation.•Bed skewed to west favors E (or N) wind over W (or S) wind to reduce stratification.•The wind causing stronger destratification generally educes more hypoxia.•Wind-induced cross-channel circulation is important in early period of wind event.•Along-channel straining by S wind has greater anoxia removal in late wind period.