Multiobjective Calibration of a Global Biogeochemical Ocean Model Against Nutrients, Oxygen, and Oxygen Minimum Zones

Global biogeochemical ocean models rely on many parameters, which govern the interaction between individual components, and their response to the physical environment. They are often assessed/calibrated against quasi-synoptic data sets of dissolved inorganic tracers. However, a good fit to one observation might not necessarily imply a good match to another. We investigate whether two different metrics—the root-mean-square error to nutrients and oxygen and a metric measuring the overlap between simulated and observed oxygen minimum zones (OMZs)—help to constrain a global biogeochemical model in different aspects of performance. Three global model optimizations are carried out. Two single-objective optimizations target the root-mean-square metric and a sum of both metrics, respectively. We then present and explore multiobjective optimization, which results in a set of compromise solutions. Our results suggest that optimal parameters for denitrification and nitrogen fixation differ when applying different metrics. Optimization against observed OMZs leads to parameters that enhance fixed nitrogen cycling; this causes too low nitrate concentrations and a too high global pelagic denitrification rate. Optimization against nutrient and oxygen concentrations leads to different parameters and a lower global fixed nitrogen turnover; this results in a worse fit to OMZs. Multiobjective optimization resolves this antagonistic effect and provides an ensemble of parameter sets, which help to address different research questions. We finally discuss how systematic model calibration can help to improve models used for projecting climate change and its effect on fisheries and climate gas emissions

An adjoint method for the assimilation of statistical characteristics into eddy-resolving ocean models

The study investigates perspectives of the parameter estimation problem with the adjoint method in eddy-resolving models. Sensitivity to initial conditions resulting from the chaotic nature of this type of model limits the direct application of the adjoint method by predictability. Prolonging the period of assimilation is accompanied by the appearance of an increasing number of secondary minima of the cost function that prevents the convergence of this method. In the framework of the Lorenz model it is shown that averaged quantities are suitable for describing invariant properties, and that secondary minima are for this type of data transformed into stochastic deviations. An adjoint method suitable for the assimilation of statistical characteristics of data and applicable on time scales beyond the predictability limit is presented. The approach assumes a greater predictability for averaged quantities. The adjoint to a prognostic model for statistical moments is employed for calculating cost function gradients that ignore the fine structure resulting from secondary minima. Coarse resolution versions of eddy-resolving models are used for this purpose. Identical twin experiments are performed with a quasigeostrophic model to evaluate the performance and limitations of this approach in improving models by estimating parameters. The wind stress curl is estimated from a simulated mean stream function. A very simple parameterization scheme for the assimilation of second-order moments is shown to permit the estimation of gradients that perform efficiently in minimizing cost functions

Metabolic suppression in thecosomatous pteropods as an effect of low temperature and hypoxia in the eastern tropical North Pacific

Author Posting. © The Author(s), 2011. This is the author's version of the work. It is posted here by permission of Springer for personal use, not for redistribution. The definitive version was published in Marine Biology 159 (2012): 1955-1967, doi:10.1007/s00227-012-1982-x.Many pteropod species in the eastern tropical north Pacific Ocean migrate vertically each day, transporting organic matter and respiratory carbon below the thermocline. These migrations take species into cold (15-10ºC) hypoxic water (< 20 µmol O2 kg-1) at depth. We measured the vertical distribution, oxygen consumption and ammonia excretion for seven species of pteropod, some of which migrate and some which remain in oxygenated surface waters throughout the day. Within the upper 200 meters of the water column, changes in water temperature result in a ~60-75% reduction in respiration for most species. All three species tested under hypoxic conditions responded to low O2 with an additional ~35-50% reduction in respiratory rate. Combined, low temperature and hypoxia suppress the metabolic rate of pteropods by ~80-90%. These results shed light on the ways in which expanding regions of hypoxia and surface ocean warming may impact pelagic ecology.This work was funded by National Science Foundation grants to K. Wishner and B. Seibel (OCE – 0526502 and OCE – 0851043) and to K. Daly (OCE – 0526545), the University of Rhode Island, and the Rhode Island Experimental Program to Stimulate Competitive Research Fellowship program.2013-06-3

Biogeochemical protocols and diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP)

The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948–2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF6) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation

Processes and patterns of oceanic nutrient limitation

Microbial activity is a fundamental component of oceanic nutrient cycles. Photosynthetic microbes, collectively termed phytoplankton, are responsible for the vast majority of primary production in marine waters. The availability of nutrients in the upper ocean frequently limits the activity and abundance of these organisms. Experimental data have revealed two broad regimes of phytoplankton nutrient limitation in the modern upper ocean. Nitrogen availability tends to limit productivity throughout much of the surface low-latitude ocean, where the supply of nutrients from the subsurface is relatively slow. In contrast, iron often limits productivity where subsurface nutrient supply is enhanced, including within the main oceanic upwelling regions of the Southern Ocean and the eastern equatorial Pacific. Phosphorus, vitamins and micronutrients other than iron may also (co-)limit marine phytoplankton. The spatial patterns and importance of co-limitation, however, remain unclear. Variability in the stoichiometries of nutrient supply and biological demand are key determinants of oceanic nutrient limitation. Deciphering the mechanisms that underpin this variability, and the consequences for marine microbes, will be a challenge. But such knowledge will be crucial for accurately predicting the consequences of ongoing anthropogenic perturbations to oceanic nutrient biogeochemistry. © 2013 Macmillan Publishers Limited. All rights reserved

Distribution and Habitat Associations of Billfish and Swordfish Larvae across Mesoscale Features in the Gulf of Mexico

Ichthyoplankton surveys were conducted in surface waters of the northern Gulf of Mexico (NGoM) over a three-year period (2006–2008) to determine the relative value of this region as early life habitat of sailfish (Istiophorus platypterus), blue marlin (Makaira nigricans), white marlin (Kajikia albida), and swordfish (Xiphias gladius). Sailfish were the dominant billfish collected in summer surveys, and larvae were present at 37.5% of the stations sampled. Blue marlin and white marlin larvae were present at 25.0% and 4.6% of the stations sampled, respectively, while swordfish occurred at 17.2% of the stations. Areas of peak production were detected and maximum density estimates for sailfish (22.09 larvae 1000 m−2) were significantly higher than the three other species: blue marlin (9.62 larvae 1000 m−2), white marlin (5.44 larvae 1000 m−2), and swordfish (4.67 larvae 1000 m−2). The distribution and abundance of billfish and swordfish larvae varied spatially and temporally, and several environmental variables (sea surface temperature, salinity, sea surface height, distance to the Loop Current, current velocity, water depth, and Sargassum biomass) were deemed to be influential variables in generalized additive models (GAMs). Mesoscale features in the NGoM affected the distribution and abundance of billfish and swordfish larvae, with densities typically higher in frontal zones or areas proximal to the Loop Current. Habitat suitability of all four species was strongly linked to physicochemical attributes of the water masses they inhabited, and observed abundance was higher in slope waters with lower sea surface temperature and higher salinity. Our results highlight the value of the NGoM as early life habitat of billfishes and swordfish, and represent valuable baseline data for evaluating anthropogenic effects (i.e., Deepwater Horizon oil spill) on the Atlantic billfish and swordfish populations

Man and the Last Great Wilderness: Human Impact on the Deep Sea

The deep sea, the largest ecosystem on Earth and one of the least studied, harbours high biodiversity and provides a wealth of resources. Although humans have used the oceans for millennia, technological developments now allow exploitation of fisheries resources, hydrocarbons and minerals below 2000 m depth. The remoteness of the deep seafloor has promoted the disposal of residues and litter. Ocean acidification and climate change now bring a new dimension of global effects. Thus the challenges facing the deep sea are large and accelerating, providing a new imperative for the science community, industry and national and international organizations to work together to develop successful exploitation management and conservation of the deep-sea ecosystem. This paper provides scientific expert judgement and a semi-quantitative analysis of past, present and future impacts of human-related activities on global deep-sea habitats within three categories: disposal, exploitation and climate change. The analysis is the result of a Census of Marine Life – SYNDEEP workshop (September 2008). A detailed review of known impacts and their effects is provided. The analysis shows how, in recent decades, the most significant anthropogenic activities that affect the deep sea have evolved from mainly disposal (past) to exploitation (present). We predict that from now and into the future, increases in atmospheric CO2 and facets and consequences of climate change will have the most impact on deep-sea habitats and their fauna. Synergies between different anthropogenic pressures and associated effects are discussed, indicating that most synergies are related to increased atmospheric CO2 and climate change effects. We identify deep-sea ecosystems we believe are at higher risk from human impacts in the near future: benthic communities on sedimentary upper slopes, cold-water corals, canyon benthic communities and seamount pelagic and benthic communities. We finalise this review with a short discussion on protection and management methods

Correction for Lebrato et al., Global variability in seawater Mg:Ca and Sr:Ca ratios in the modern ocean

4 pages, 5 figures.-- Correction Global variability in seawater Mg:Ca and Sr:Ca ratios in the modern ocean; Proceedings of the National Academy of Sciences of the USA 117(36): 22281-22292 (2020); doi: 10.1073/pnas.1918943117; http://hdl.handle.net/10261/221953The authors wish to note the following: “This study’s seawater Sr:Ca values were systematically low as a consequence of normalization to another published low value for the International Association for the Physical Sciences of the Oceans (IAPSO) (1). IAPSO has been used at the Ocean Drilling Program, Texas A&M University (ODP-TAMU) (http://www-odp.tamu.edu/), and is still being used as the primary standard for elemental composition of seawater/interstitial water. Consequently, our seawater value of Sr:Ca = 8.28 mmol:mol was systematically low by approx. 3.70%, if we accept seawater Sr:Ca 8.60 mmol:mol as the recommended value for IAPSO North Atlantic surface water salinity standard. The uncertainty budget should be expanded including the uncertainty of IAPSO composition. The largest contribution to expanded uncertainty of our data comes from the uncertainty of the IAPSO reference composition, which is 3.29% using all published values. This will result in 3.30% (1 SD) expanded uncertainty for seawater Sr:Ca (and 0.5%, for seawater Mg:Ca) of the entire data set with respect to accuracy. We have corrected all seawater Sr:Ca values with a factor of 1.0243 in all our tables (e.g., SI Appendix, Table S1 averages) and in the figures (Fig. 4, Fig. 5), where a ratio was used. Note that the seawater Sr:Ca % changes are small, thus changes are hardly noticeable on large displays (e.g., Figures), but they can be seen in the tables and averages/SD calculations. Seawater Sr:Ca ratios are also corrected in the main text where relevantPeer reviewe

A vertically resolved model for phytoplankton aggregation

This work presents models of the vertical distribution and flux of phytoplankton aggregates, including changes with time in the distribution of aggregate sizes and sinking speeds. The distribution of sizes is described by two parameters, the mass and number of aggregates, which greatly reduces the computational cost of the models. Simple experiments demonstrate the effects of aggregation on the timing and depth distribution of primary production and export. A more detailed ecological model is applied to sites in the Arabian Sea; it demonstrates that aggregation can be important for deep sedimentation even when its effect on surface concentrations is small, and it presents the difference in timing between settlement of aggregates and fecal pellets

The 5th edition of the World Health Organization classification of haematolymphoid tumours: lymphoid neoplasms

We herein present an overview of the upcoming 5th edition of the World Health Organization Classification of Haematolymphoid Tumours focussing on lymphoid neoplasms. Myeloid and histiocytic neoplasms will be presented in a separate accompanying article. Besides listing the entities of the classification, we highlight and explain changes from the revised 4th edition. These include reorganization of entities by a hierarchical system as is adopted throughout the 5th edition of the WHO classification of tumours of all organ systems, modification of nomenclature for some entities, revision of diagnostic criteria or subtypes, deletion of certain entities, and introduction of new entities, as well as inclusion of tumour-like lesions, mesenchymal lesions specific to lymph node and spleen, and germline predisposition syndromes associated with the lymphoid neoplasms