Assimilating high-resolution salinity data into a model of a partially mixed estuary

[1] A three-dimensional circulation model of the Chesapeake Bay is used to validate a simple data assimilation scheme, using high-resolution salinity data acquired from a ship-towed undulating vehicle (a Scanfish). The simulation period spans the entire year of 1995 during which the high-resolution Scanfish data were available in July and October, lasting a few days each. Since Scanfish data were irregularly distributed in time and space, only salinity fields are nudged in the model for simplicity. Model improvements through data assimilation are evaluated from a pair of experiments: one with data assimilation and one without. Data from scattered Chesapeake Bay Program monitoring stations and a few stations maintained by the National Ocean Service inside the bay are used independently to check the model performance. In general, the simple assimilation scheme leads to visibly improved density structures in the upper and middle reaches of the bay. The improvement in the lower bay is equally pronounced after data assimilation but diminishes in a shorter timescale because of faster flushing from the adjacent coastal ocean. Moderate to weak nudging normally enhances the gravitational circulation. Strong nudging may produce transient overshooting, during which the gravitational circulation is renewed vigorously

Chesapeake Bay Nitrogen Fluxes Derived From a Land-Estuarine Ocean Biogeochemical Modeling System: Model Description, Evaluation, and Nitrogen Budgets

The Chesapeake Bay plays an important role in transforming riverine nutrients before they are exported to the adjacent continental shelf. Although the mean nitrogen budget of the Chesapeake Bay has been previously estimated from observations, uncertainties associated with interannually varying hydrological conditions remain. In this study, a land-estuarine-ocean biogeochemical modeling system is developed to quantify Chesapeake riverine nitrogen inputs, within-estuary nitrogen transformation processes and the ultimate export of nitrogen to the coastal ocean. Model skill was evaluated using extensive in situ and satellite-derived data, and a simulation using environmental conditions for 2001–2005 was conducted to quantify the Chesapeake Bay nitrogen budget. The 5 year simulation was characterized by large riverine inputs of nitrogen (154 × 109 g N yr−1) split roughly 60:40 between inorganic:organic components. Much of this was denitrified (34 × 109 g N yr−1) and buried (46 × 109 g N yr−1) within the estuarine system. A positive net annual ecosystem production for the bay further contributed to a large advective export of organic nitrogen to the shelf (91 × 109 g N yr−1) and negligible inorganic nitrogen export. Interannual variability was strong, particularly for the riverine nitrogen fluxes. In years with higher than average riverine nitrogen inputs, most of this excess nitrogen (50–60%) was exported from the bay as organic nitrogen, with the remaining split between burial, denitrification, and inorganic export to the coastal ocean. In comparison to previous simulations using generic shelf biogeochemical model formulations inside the estuary, the estuarine biogeochemical model described here produced more realistic and significantly greater exports of organic nitrogen and lower exports of inorganic nitrogen to the shelf

Predicting the distribution of Vibrio vulnificus in Chesapeake Bay

Vibrio vulnificus is a gram-negative pathogenic bacterium endemic to coastal waters worldwide, and a leading cause of seafood related mortality. Because of human health concerns, understanding the ecology of the species and potentially predicting its distribution is of great importance. We evaluated and applied a previously published qPCR assay to water samples (n = 235) collected from the main-stem of the Chesapeake Bay (2007 – 2008) by Maryland and Virginia State water quality monitoring programs. Results confirmed strong relationships between the likelihood of Vibrio vulnificus presence and both temperature and salinity that were used to develop a logistic regression model. The habitat model demonstrated a high degree of concordance (93%), and robustness as subsequent bootstrapping (n=1000) did not change model output (P > 0.05). We forced this empirical habitat model with temperature and salinity predictions generated by a regional hydrodynamic modeling system to demonstrate its utility in future pathogen forecasting efforts in the Chesapeake Bay

Forecasting Prorocentrum minimum blooms in the Chesapeake Bay using empirical habitat models

Aquaculturists, local beach managers, and other stakeholders require forecasts of harmful biotic events, so they can assess and respond to health threats when harmful algal blooms (HABs) are present. Based on this need, we are developing empirical habitat suitability models for a variety of Chesapeake Bay HABs to forecast their occurrence based on a set of physical-biogeochemical environmental conditions, and start with the dinoflagellate Prorocentrum minimum (also known as P. cordatum).To identify an optimal set of environmental variables to forecast P. minimum blooms, we first assumed a linear relationship between the environmental variables and the inverse of the logistic function used to forecast the likelihood of bloom presence, and repeated the method using more than 16,000 combinations of variables. By comparing goodness-of-fit, we found water temperature, salinity, pH, solar irradiance, and total organic nitrogen represented the most suitable set of variables. The resulting algorithm forecasted P. minimum blooms with an overall accuracy of 78%, though with a significant variability ~ 30-90% depending on region and season. To understand this variability and improve model performance, we incorporated nonlinear effects into the model by implementing a generalized additive model. Even without considering interactions between the five variables used to train the model, this yielded an increase in overall model accuracy (~ 81%) due to the model’s ability to refine the regions in which P. minimum blooms occurred. Including nonlinear interactions increased the overall model accuracy even further (~ 85%) by accounting for seasonality in the interaction between solar irradiance and water temperature. Our findings suggest that the influence of predictors of these blooms change in time and space, and that model complexity impacts the model performance and our interpretation of the driving factors causing P. minimum blooms. Apart from their forecasting potential, our results may be particularly useful when constructing explicit relationships between environmental conditions and P. minimum presence in mechanistic models

Pelagic Functional Group Modeling: Progress, Challenges and Prospects

In this paper, we review the state of the art and major challenges in current efforts to incorporate biogeochemical functional groups into models that can be applied on basin-wide and global scales, with an emphasis on models that might ultimately be used to predict how biogeochernical cycles in the ocean will respond to global warming. We define the term biogeochemical functional group to refer to groups of organisms that mediate specific chemical reactions in the ocean. Thus, according to this definition, functional groups have no phylogenetic meaning-these are composed of many different species with common biogeochemical functions. Substantial progress has been made in the last decade toward quantifying the rates of these various functions and understanding the factors that control them. For some of these groups, we have developed fairly sophisticated models that incorporate this understanding, e.g. for diazotrophs (e.g. Trichodesmium), silica producers (diatoms) and calcifiers (e.g. coccolithophorids and specifically Emiliania huxleyi). However, current representations of nitrogen fixation and calcification are incomplete, i.e., based primarily upon models of Trichodesmium and E huxleyi, respectively, and many important functional groups have not yet been considered in open-ocean biogeochemical models. Progress has been made over the last decade in efforts to simulate dimethylsulfide (DMS) production and cycling (i.e., by dinoflagellates and prymnesiophytes) and denitrification, but these efforts are still in their infancy, and many significant problems remain. One obvious gap is that virtually all functional group modeling efforts have focused on autotrophic microbes, while higher trophic levels have been completely ignored. It appears that in some cases (e.g., calcification), incorporating higher trophic levels may be essential not only for representing a particular biogeochemical reaction, but also for modeling export. Another serious problem is our tendency to model the organisms for which we have the most validation data (e.g., E huxleyi and Trichodesmium) even when they may represent only a fraction of the biogeochemical functional group we are trying to represent. When we step back and look at the paleo-oceanographic record, it suggests that oxygen concentrations have played a central role in the evolution and emergence of many of the key functional groups that influence biogeochemical cycles in the present-day ocean. However, more subtle effects are likely to be important over the next century like changes in silicate supply or turbulence that can influence the relative success of diatoms versus dinoflagellates, coccolithophorids and diazotrophs. In general, inferences drawn from the paleo-oceanographic record and theoretical work suggest that global warming will tend to favor the latter because it will give rise to increased stratification. However, decreases in pH and Fe supply could adversely impact coccolithophorids and diazotrophs in the future. It may be necessary to include explicit dynamic representations of nitrogen fixation, denitrification, silicification and calcification in our models if our goal is predicting the oceanic carbon cycle in the future, because these processes appear to play a very significant role in the carbon cycle of the present-day ocean and they are sensitive to climate change. Observations and models suggest that it may also be necessary to include the DMS cycle to predict future climate, though the effects are still highly uncertain. We have learned a tremendous amount about the distributions and biogeochemical impact of bacteria in the ocean in recent years, yet this improved understanding has not yet been incorporated into many of our models. All of these considerations lead us toward the development of increasingly complex models. However, recent quantitative model intercomparison studies suggest that continuing to add complexity and more functional groups to our ecosystem models may lead to decreases in predictive ability if the models are not properly constrained with available data. We also caution that capturing the present-day variability tells us little about how well a particular model can predict the future. If our goal is to develop models that can be used to predict how the oceans will respond to global warming, then we need to make more rigorous assessments of predictive skill using the available data

Challenges associated with modeling low-oxygen waters in Chesapeake Bay: a multiple model comparison

As three-dimensional (3-D) aquatic ecosystem models are used more frequently for operational water quality forecasts and ecological management decisions, it is important to understand the relative strengths and limitations of existing 3-D models of varying spatial resolution and biogeochemical complexity. To this end, 2-year simulations of the Chesapeake Bay from eight hydrodynamic-oxygen models have been statistically compared to each other and to historical monitoring data. Results show that although models have difficulty resolving the variables typically thought to be the main drivers of dissolved oxygen variability (stratification, nutrients, and chlorophyll), all eight models have significant skill in reproducing the mean and seasonal variability of dissolved oxygen. In addition, models with constant net respiration rates independent of nutrient supply and temperature reproduced observed dissolved oxygen concentrations about as well as much more complex, nutrient-dependent biogeochemical models. This finding has significant ramifications for short-term hypoxia forecasts in the Chesapeake Bay, which may be possible with very simple oxygen parameterizations, in contrast to the more complex full biogeochemical models required for scenario-based forecasting. However, models have difficulty simulating correct density and oxygen mixed layer depths, which are important ecologically in terms of habitat compression. Observations indicate a much stronger correlation between the depths of the top of the pycnocline and oxycline than between their maximum vertical gradients, highlighting the importance of the mixing depth in defining the region of aerobic habitat in the Chesapeake Bay when low-oxygen bottom waters are present. Improvement in hypoxia simulations will thus depend more on the ability of models to reproduce the correct mean and variability of the depth of the physically driven surface mixed layer than the precise magnitude of the vertical density gradient

United States contributions to the Second International Indian Ocean Expedition (US IIOE-2)

From the Preface: The purpose of this document is to motivate and coordinate U.S. participation in the Second International Indian Ocean Expedition (IIOE-2) by outlining a core set of research priorities that will accelerate our understanding of geologic, oceanic, and atmospheric processes and their interactions in the Indian Ocean. These research priorities have been developed by the U.S. IIOE-2 Steering Committee based on the outcomes of an interdisciplinary Indian Ocean science workshop held at the Scripps Institution of Oceanography on September 11-13, 2017. The workshop was attended by 70 scientists with expertise spanning climate, atmospheric sciences, and multiple sub-disciplines of oceanography. Workshop participants were largely drawn from U.S. academic institutions and government agencies, with a few experts invited from India, China, and France to provide a broader perspective on international programs and activities and opportunities for collaboration. These research priorities also build upon the previously developed International IIOE-2 Science Plan and Implementation Strategy. Outcomes from the workshop are condensed into five scientific themes: Upwelling, inter-ocean exchanges, monsoon dynamics, inter-basin contrasts, marine geology and the deep ocean. Each theme is identified with priority questions that the U.S. research community would like to address and the measurements that need to be made in the Indian Ocean to address them.We thank the following organizations and programs for financial contributions, support and endorsement: the U.S. National Oceanic and Atmospheric Administration; the U.S. Ocean Carbon and Biogeochemistry program funded by the National Science Foundation and the National Aeronautics and Space Administration; the NASA Physical Oceanography Program; Scripps Institution of Oceanography; and the Indo-US Science and Technology Forum

Assessment of skill and portability in regional marine biogeochemical models: Role of multiple planktonic groups

Application of biogeochemical models to the study of marine ecosystems is pervasive, yet objective quantification of these models' performance is rare. Here, 12 lower trophic level models of varying complexity are objectively assessed in two distinct regions (equatorial Pacific and Arabian Sea). Each model was run within an identical one-dimensional physical framework. A consistent variational adjoint implementation assimilating chlorophyll-a, nitrate, export, and primary productivity was applied and the same metrics were used to assess model skill. Experiments were performed in which data were assimilated from each site individually and from both sites simultaneously. A cross-validation experiment was also conducted whereby data were assimilated from one site and the resulting optimal parameters were used to generate a simulation for the second site. When a single pelagic regime is considered, the simplest models fit the data as well as those with multiple phytoplankton functional groups. However, those with multiple phytoplankton functional groups produced lower misfits when the models are required to simulate both regimes using identical parameter values. The cross-validation experiments revealed that as long as only a few key biogeochemical parameters were optimized, the models with greater phytoplankton complexity were generally more portable. Furthermore, models with multiple zooplankton compartments did not necessarily outperform models with single zooplankton compartments, even when zooplankton biomass data are assimilated. Finally, even when different models produced similar least squares model-data misfits, they often did so via very different element flow pathways, highlighting the need for more comprehensive data sets that uniquely constrain these pathways

A road map to IndOOS-2 better observations of the rapidly warming Indian Ocean

Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Bulletin of the American Meteorological Society 101(11), (2020): E1891-E1913, https://doi.org/10.1175/BAMS-D-19-0209.1The Indian Ocean Observing System (IndOOS), established in 2006, is a multinational network of sustained oceanic measurements that underpin understanding and forecasting of weather and climate for the Indian Ocean region and beyond. Almost one-third of humanity lives around the Indian Ocean, many in countries dependent on fisheries and rain-fed agriculture that are vulnerable to climate variability and extremes. The Indian Ocean alone has absorbed a quarter of the global oceanic heat uptake over the last two decades and the fate of this heat and its impact on future change is unknown. Climate models project accelerating sea level rise, more frequent extremes in monsoon rainfall, and decreasing oceanic productivity. In view of these new scientific challenges, a 3-yr international review of the IndOOS by more than 60 scientific experts now highlights the need for an enhanced observing network that can better meet societal challenges, and provide more reliable forecasts. Here we present core findings from this review, including the need for 1) chemical, biological, and ecosystem measurements alongside physical parameters; 2) expansion into the western tropics to improve understanding of the monsoon circulation; 3) better-resolved upper ocean processes to improve understanding of air–sea coupling and yield better subseasonal to seasonal predictions; and 4) expansion into key coastal regions and the deep ocean to better constrain the basinwide energy budget. These goals will require new agreements and partnerships with and among Indian Ocean rim countries, creating opportunities for them to enhance their monitoring and forecasting capacity as part of IndOOS-2.We thank the World Climate Research Programme (WCRP) and its core project on Climate and Ocean: Variability, Predictability and Change (CLIVAR), the Indian Ocean Global Ocean Observing System (IOGOOS), the Intergovernmental Oceanographic Commission of UNESCO (IOC-UNESCO), the Integrated Marine Biosphere Research (IMBeR) project, the U.S. National Oceanic and Atmospheric Administration (NOAA), and the International Union of Geodesy and Geophysics (IUGG) for providing the financial support to bring international scientists together to conduct this review. We thank the members of the independent review board that provided detailed feedbacks on the review report that is summarized in this article: P. E. Dexter, M. Krug, J. McCreary, R. Matear, C. Moloney, and S. Wijffels. PMEL Contribution 5041. C. Ummenhofer acknowledges support from The Andrew W. Mellon Foundation Award for Innovative Research.2021-05-0

Modeling the seasonal autochthonous sources of dissolved organic carbon and nitrogen in the upper Chesapeake Bay

In this paper we investigate the seasonal autochthonous sources of dissolved organic carbon (DOC) and nitrogen (DON) in the euphotic zone at a station in the upper Chesapeake Bay using a new mass-based ecosystem model. Important features of the model are: (1) carbon and nitrogen are incorporated by means of a set of fixed and varying C:N ratios; (2) dissolved organic matter (DOM) is separated into labile, semi-labile, and refractory pools for both C and N; (3) the production and consumption of DOM is treated in detail; and (4) seasonal observations of light, temperature, nutrients, and surface layer circulation are used to physically force the model. The model reasonably reproduces the mean observed seasonal concentrations of nutrients, DOM, plankton biomass, and chlorophyll a. The results suggest that estuarine DOM production is intricately tied to the biomass concentration, ratio, and productivity of phytoplankton, zooplankton, viruses, and bacteria. During peak spring productivity phytoplankton exudation and zooplankton sloppy feeding are the most important autochthonous sources of DOM. In the summer when productivity peaks again, autochthonous sources of DOM are more diverse and, in addition to phytoplankton exudation, important ones include viral lysis and the decay of detritus. The potential importance of viral decay as a source of bioavailable DOM from within the bulk DOM pool is also discussed. The results also highlight the importance of some poorly constrained processes and parameters. Some potential improvements and remedies are suggested. Sensitivity studies on selected parameters are also reported and discussed