scholarly journals Diel light cycles affect phytoplankton competition in the global ocean

2021 ◽  
Author(s):  
Ioannis Tsakalakis ◽  
Michael J. Follows ◽  
Stephanie Dutkiewicz ◽  
Christopher L. Follett ◽  
Joseph J. Vallino

Light, essential for photosynthesis, is present in two periodic cycles in nature: seasonal and diel. Although seasonality of light is typically resolved in ocean ecosystem and biogeochemistry models because of its significance for seasonal succession and biogeography of phytoplankton, the diel light cycle is generally not resolved. Here we use a three-dimensional global ocean model and compare high temporal resolution simulations with and without diel light cycles. The model simulates 15 phytoplankton types of different cell size, encompassing two broad ecological strategies: small cells with high nutrient affinity (gleaners) and larger cells with high maximal growth rate (opportunists). Both are grazed by zooplankton and limited by nitrogen, phosphorus and iron. Simulations show that diel cycles of light induce diel cycles in phytoplankton populations and limiting nutrients in the global ocean. Diel nutrient cycles are associated with higher concentration of limiting nutrients by up to 200% at low latitudes (-40 to 40), a process that increases opportunists biomass by up to 50%. Size classes with the highest maximal growth rates from both gleaner and opportunist groups are favored the most by diel light cycles. This mechanism weakens as latitude increases because the effects of the seasonal cycle dominate over those of the diel cycle. The present work shows that resource competition under diel light cycles has a significant impact on phytoplankton biogeography, indicating the necessity of resolving diel processes in global ocean models.

2017 ◽  
Author(s):  
Ralph Timmermann ◽  
Sebastian Goeller

Abstract. A Regional Antarctic and Global Ocean (RAnGO) model has been developed to study the interaction between the world ocean and the Antarctic ice sheet. The coupled model is based on a global implementation of the Finite Element Sea-ice Ocean Model (FESOM) with a mesh refinement in the Southern Ocean, particularly in its marginal seas and in the sub-ice shelf cavities. The cryosphere is represented by a regional setup of the ice flow model RIMBAY comprising the Filchner-Ronne Ice Shelf and the grounded ice in its catchment area up to the ice divides. At the base of the RIMBAY ice shelf, melt rates from FESOM's ice-shelf component are supplied. RIMBAY returns ice thickness and the position of the grounding line. The ocean model uses a pre-computed mesh to allow for an easy adjustment of the model domain to a varying cavity geometry. RAnGO simulations with a 20th-century climate forcing yield realistic basal melt rates and a quasi-stable grounding line position close to the presently observed state. In a centennial-scale warm-water-inflow scenario, the model suggests a substantial thinning of the ice shelf and a local retreat of the grounding line. The potentially negative feedback from ice-shelf thinning through a rising in-situ freezing temperature is more than outweighed by the increasing water column thickness in the deepest parts of the cavity. Compared to a control simulation with fixed ice-shelf geometry, the coupled model thus yields a slightly stronger increase of ice-shelf basal melt rates.


2014 ◽  
Vol 11 (10) ◽  
pp. 2635-2643 ◽  
Author(s):  
R. Death ◽  
J. L. Wadham ◽  
F. Monteiro ◽  
A. M. Le Brocq ◽  
M. Tranter ◽  
...  

Abstract. Southern Ocean (SO) marine primary productivity (PP) is strongly influenced by the availability of iron in surface waters, which is thought to exert a significant control upon atmospheric CO2 concentrations on glacial/interglacial timescales. The zone bordering the Antarctic Ice Sheet exhibits high PP and seasonal plankton blooms in response to light and variations in iron availability. The sources of iron stimulating elevated SO PP are in debate. Established contributors include dust, coastal sediments/upwelling, icebergs and sea ice. Subglacial meltwater exported at the ice margin is a more recent suggestion, arising from intense iron cycling beneath the ice sheet. Icebergs and subglacial meltwater may supply a large amount of bioavailable iron to the SO, estimated in this study at 0.07–0.2 Tg yr−1. Here we apply the MIT global ocean model (Follows et al., 2007) to determine the potential impact of this level of iron export from the ice sheet upon SO PP. The export of iron from the ice sheet raises modelled SO PP by up to 40%, and provides one plausible explanation for seasonally very high in situ measurements of PP in the near-coastal zone. The impact on SO PP is greatest in coastal regions, which are also areas of high measured marine PP. These results suggest that the export of Antarctic runoff and icebergs may have an important impact on SO PP and should be included in future biogeochemical modelling.


2003 ◽  
Vol 15 (1) ◽  
pp. 41-46 ◽  
Author(s):  
ROBIN ROBERTSON ◽  
AIKE BECKMANN ◽  
HARTMUT HELLMER

In certain regions of the Southern Ocean, tidal energy is believed to foster the mixing of different water masses, which eventually contribute to the formation of deep and bottom waters. The Ross Sea is one of the major ventilation sites of the global ocean abyss and a region of sparse tidal observations. We investigated M2 tidal dynamics in the Ross Sea using a three-dimensional sigma coordinate model, the Regional Ocean Model System (ROMS). Realistic topography and hydrography from existing observational data were used with a single tidal constituent, the semi-diurnal M2. The model fields faithfully reproduced the major features of the tidal circulation and had reasonable agreement with ten existing tidal elevation observations and forty-two existing tidal current measurements. The differences were attributed primarily to topographic errors. Internal tides were generated at the continental shelf/slope break and other areas of steep topography. Strong vertical shears in the horizontal velocities occurred under and at the edges of the Ross Ice Shelf and along the continental shelf/slope break. Estimates of lead formation based on divergence of baroclinic velocities were significantly higher than those based on barotrophic velocities, reaching over 10% at the continental shelf/slope break.


2016 ◽  
Author(s):  
Stephen M. Griffies ◽  
Gokhan Danabasoglu ◽  
Paul J. Durack ◽  
Alistair J. Adcroft ◽  
V. Balaji ◽  
...  

Abstract. The Ocean Model Intercomparison Project (OMIP) aims to provide a framework for evaluating, understanding, and improving the ocean and sea-ice components of global climate and earth system models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses these aims in two complementary manners: (A) by providing an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing, (B) by providing a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) offering details for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows that of the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II have become the standard method to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP (Scenario MIP), as well as the ocean-sea ice OMIP simulations. The bulk of this paper offers scientific rationale for saving these diagnostics.


2013 ◽  
Vol 10 (7) ◽  
pp. 12551-12570 ◽  
Author(s):  
R. Death ◽  
J. L. Wadham ◽  
F. Monteiro ◽  
A. M. Le Brocq ◽  
M. Tranter ◽  
...  

Abstract. Southern Ocean (SO) marine primary productivity (PP) is strongly influenced by the availability of iron in surface waters, which is thought to exert a significant control upon atmospheric CO2 concentrations on glacial/interglacial timescales. The zone bordering the Antarctic Ice Sheet exhibits high PP and seasonal plankton blooms in response to light and variations in iron availability. The sources of iron stimulating elevated SO PP are in debate. Established contributors include dust, coastal sediments/upwelling, icebergs and sea ice. Subglacial meltwater exported at the ice margin is a more recent suggestion, arising from intense iron cycling beneath the ice sheet. Icebergs and subglacial meltwater may supply a large amount of bioavailable iron to the SO, estimated in this study at 0.07–1.0 Tg yr−1. Here we apply the MIT global ocean model (Follows et al., 2007) to determine the potential impact of this level of iron export from the ice sheet upon SO PP. The export of iron from the ice sheet raises modelled SO PP by up to 40%, and provides one plausible explanation for very high seasonally observed PP in the near-coastal zone. The impact on SO PP is greatest in coastal regions, which are also areas of high observed marine PP. These results suggest that the export of Antarctic runoff and icebergs may have an important impact on SO PP and should be included in future biogeochemical modelling.


2020 ◽  
Vol 13 (11) ◽  
pp. 5465-5483
Author(s):  
Clément Bricaud ◽  
Julien Le Sommer ◽  
Gurvan Madec ◽  
Christophe Calone ◽  
Julie Deshayes ◽  
...  

Abstract. Ocean biogeochemical models are key tools for both scientific and operational applications. Nevertheless the cost of these models is often expensive because of the large number of biogeochemical tracers. This has motivated the development of multi-grid approaches where ocean dynamics and tracer transport are computed on grids of different spatial resolution. However, existing multi-grid approaches to tracer transport in ocean modelling do not allow the computation of ocean dynamics and tracer transport simultaneously. This paper describes a new multi-grid approach developed for accelerating the computation of passive tracer transport in the Nucleus for European Modelling of the Ocean (NEMO) ocean circulation model. In practice, passive tracer transport is computed at runtime on a grid with coarser spatial resolution than the hydrodynamics, which reduces the CPU cost of computing the evolution of tracers. We describe the multi-grid algorithm, its practical implementation in the NEMO ocean model, and discuss its performance on the basis of a series of sensitivity experiments with global ocean model configurations. Our experiments confirm that the spatial resolution of hydrodynamical fields can be coarsened by a factor of 3 in both horizontal directions without significantly affecting the resolved passive tracer fields. Overall, the proposed algorithm yields a reduction by a factor of 7 of the overhead associated with running a full biogeochemical model like PISCES (with 24 passive tracers). Propositions for further reducing this cost without affecting the resolved solution are discussed.


2020 ◽  
Author(s):  
Casimir de Lavergne ◽  
Clément Vic ◽  
Gurvan Madec ◽  
Fabien Roquet ◽  
Amy Waterhouse ◽  
...  

<p>Vertical mixing is often regarded as the Achilles’ heel of ocean models. In particular, few models include a comprehensive and energy-constrained parameterization of mixing by internal ocean tides. Here, we present an energy-conserving mixing scheme which accounts for the local breaking of high-mode internal tides and the distant dissipation of low-mode internal tides. The scheme relies on four static two-dimensional maps of internal tide dissipation, constructed using mode-by-mode Lagrangian tracking of energy beams from sources to sinks. Each map is associated with a distinct dissipative process and a corresponding vertical structure. Applied to an observational climatology of stratification, the scheme produces a global three-dimensional map of dissipation which compares well with available microstructure observations and with upper-ocean finestructure mixing estimates. Implemented in the NEMO global ocean model, the scheme improves the representation of deep water-mass transformation and obviates the need for a constant background diffusivity.</p>


Ocean Science ◽  
2009 ◽  
Vol 5 (4) ◽  
pp. 403-419 ◽  
Author(s):  
C. Skandrani ◽  
J.-M. Brankart ◽  
N. Ferry ◽  
J. Verron ◽  
P. Brasseur ◽  
...  

Abstract. In the context of stand alone ocean models, the atmospheric forcing is generally computed using atmospheric parameters that are derived from atmospheric reanalysis data and/or satellite products. With such a forcing, the sea surface temperature that is simulated by the ocean model is usually significantly less accurate than the synoptic maps that can be obtained from the satellite observations. This not only penalizes the realism of the ocean long-term simulations, but also the accuracy of the reanalyses or the usefulness of the short-term operational forecasts (which are key GODAE and MERSEA objectives). In order to improve the situation, partly resulting from inaccuracies in the atmospheric forcing parameters, the purpose of this paper is to investigate a way of further adjusting the state of the atmosphere (within appropriate error bars), so that an explicit ocean model can produce a sea surface temperature that better fits the available observations. This is done by performing idealized assimilation experiments in which Mercator-Ocean reanalysis data are considered as a reference simulation describing the true state of the ocean. Synthetic observation datasets for sea surface temperature and salinity are extracted from the reanalysis to be assimilated in a low resolution global ocean model. The results of these experiments show that it is possible to compute piecewise constant parameter corrections, with predefined amplitude limitations, so that long-term free model simulations become much closer to the reanalysis data, with misfit variance typically divided by a factor 3. These results are obtained by applying a Monte Carlo method to simulate the joint parameter/state prior probability distribution. A truncated Gaussian assumption is used to avoid the most extreme and non-physical parameter corrections. The general lesson of our experiments is indeed that a careful specification of the prior information on the parameters and on their associated uncertainties is a key element in the computation of realistic parameter estimates, especially if the system is affected by other potential sources of model errors.


Author(s):  
A. Paul ◽  
S. Mulitza ◽  
J. Pätzold ◽  
T. Wolff

Sign in / Sign up

Export Citation Format

Share Document