Skip to content

Latest commit

 

History

History
245 lines (175 loc) · 10.6 KB

global_oce_biogeo.rst

File metadata and controls

245 lines (175 loc) · 10.6 KB
Raw

Biogeochemistry Simulation

(in directory: :filelink:`verification/tutorial_global_oce_biogeo/`)

Overview

This model overlays the dissolved inorganic carbon biogeochemistry model (:filelink:`pkg/dic`) over a 2.8o global physical model. The physical model has 15 levels, and is forced with a climatological annual cycle of surface wind stresses (Trenberth et al. 1989 :cite:`trenberth:89`, surface heat and freshwater fluxes (Jiang et al. 1999 :cite:`jiang:99`) with additional relaxation toward climatological sea surface temperature and salinity (Levitus and Boyer (1994a,b) :cite:`levitus:94a,levitus:94b`). It uses the Gent and McWilliams (1990) :cite:`gen-mcw:90` eddy parameterization scheme, has an implicit free-surface, implicit vertical diffusion and uses the convective adjustment scheme.

The biogeochemical model considers the coupled cycles of carbon, oxygen, phosphorus and alkalinity. A simplified parameterization of biological production is used, limited by the availability of light and phosphate. A fraction of this productivity enters the dissolved organic pool pool, which has an e-folding timescale for remineralization of 6 months (following Yamanaka and Tajika 1997 :cite:`yamanaka:97`). The remaining fraction of this productivity is instantaneously exported as particulate to depth (Yamanaka and Tajika 1997 :cite:`yamanaka:97`) where it is remineralized according to the empirical power law relationship determined by Martin et al. (1987]) :cite:`martin:87`. The fate of carbon is linked to that of phosphorus by the Redfield ratio. Carbonate chemistry is explicitly solved (see Follow et al. 2006) :cite:`follows:06`) and the air-sea exchange of CO2 is parameterized with a uniform gas transfer coefficient following Wanninkhof (1992) :cite:`wannink:92`. Oxygen is also linked to phosphorus by the Redfield ratio, and oxygen air-sea exchange also follows Wanninkhof (1992) :cite:`wannink:92`. For more details see Dutkiewicz et al. (2005) :cite:`dutkiewicz:05`.

The example setup described here shows the physical model after 5900 years of spin-up and the biogeochemistry after 2900 years of spin-up. The biogeochemistry is at a pre-industrial steady-state (atmospheric ppmv is kept at 278). Five tracers are resolved: dissolved inorganic carbon (DIC), alkalinity (ALK), phosphate (PO4), dissolved organic phosphorus (DOP) and dissolved oxygen (O2).

Modeled annual mean air-sea CO2

Modeled annual mean air-sea CO2 flux (mol C m-2 y-1) for pre-industrial steady-state. Positive indicates flux of CO2 from ocean to the atmosphere (out-gassing), contour interval is 1 mol C m-2 y-1.

Equations Solved

The physical ocean model velocity and diffusivities are used to redistribute the 5 tracers within the ocean. Additional redistribution comes from chemical and biological sources and sinks. For any tracer A:

\frac{\partial A}{\partial t}=- \nabla  \cdot (\vec{u^{*}} A)+ \nabla  \cdot
  (\mathbf{K}\nabla A)+S_A \nonumber

where \vec{u^{*}} is the transformed Eulerian mean circulation (which includes Eulerian and eddy-induced advection), \mathbf{K} is the mixing tensor, and S_A are the sources and sinks due to biological and chemical processes.

The sources and sinks are:

\begin{aligned}
S_{DIC} & =  F_{CO_2} + V_{CO_2} + r_{C:P} S_{PO_4}  + J_{Ca} \\
S_{ALK} & =  V_{ALK}-r_{N:P} S_{PO_4}  + 2 J_{Ca}  \\
S_{PO_4}& =  -f_{DOP} J_{prod} - \frac{\partial F_P}{\partial z} + \kappa_{remin} [DOP]\\
S_{DOP} & =  f_{DOP} J_{prod} -\kappa_{remin} [DOP] \\
S_{O_2} & = \left\{ \begin{array}{ll}
               -r_{O:P} S_{PO_4} & \mbox{if $O_2>O_{2crit}$} \\
                0  & \mbox{if $O_2<O_{2crit}$}
                      \end{array}
              \right. \end{aligned}

where:

  • F_{CO_2} is the flux of CO2 from the ocean to the atmosphere
  • V_{CO_2} is “virtual flux” due to changes in DIC due to the surface freshwater fluxes
  • r_{C:P} is Redfield ratio of carbon to phosphorus
  • J_{Ca} includes carbon removed from surface due to calcium carbonate formation and subsequent cumulation of the downward flux of CaCO_3
  • V_{ALK} is “virtual flux” due to changes in alkalinity due to the surface freshwater fluxes
  • r_{N:P} Redfield ratio is nitrogen to phosphorus
  • f_{DOP} is fraction of productivity that remains suspended in the water column as dissolved organic phosphorus
  • J_{prod} is the net community productivity
  • \frac{\partial F_P}{\partial z} is the accumulation of remineralized phosphorus with depth
  • \kappa_{remin} is rate with which DOP remineralizes back to PO_4
  • F_{O_2} is air-sea flux of oxygen
  • r_{O:P} is Redfield ratio of oxygen to phosphorus
  • O_{2crit} is a critical level below which oxygen consumption if halted

These terms (for the first four tracers) are described more in Dutkiewicz et al. (2005) :cite:`dutkiewicz:05` and by McKinley et al. (2004) :cite:`mckinley:04` for the terms relating to oxygen.

Code configuration

The modifications to the code (in :filelink:`verification/tutorial_global_oce_biogeo/code`) are:

The input fields needed for this run (in :filelink:`verification/tutorial_global_oce_biogeo/input`) are:

Running the example

As the model is set up to run in the verification experiment, it only runs for 4 timesteps (2 days) and outputs data at the end of this short run. For a more informative run, you will need to run longer. As set up, this model starts from a pre-spun up state and initializes physical fields and the biogeochemical tracers from the pickup files.

Physical data (e.g., S,T, velocities etc) will be output as for any regular ocean run. The biogeochemical output are: