Home Contact Us Site Map  
 
       
    next up previous contents
Next: 3.19.2 Code configuration Up: 3.19 Sensitivity of Air-Sea Previous: 3.19 Sensitivity of Air-Sea   Contents

Subsections


3.19.1 Overview of the experiment

We describe an adjoint sensitivity analysis of out-gassing from the ocean into the atmosphere of a carbon-like tracer injected into the ocean interior (see Hill et al. [2002]).


3.19.1.1 Passive tracer equation

For this work MITgcm was augmented with a thermodynamically inactive tracer, $ C$. Tracer residing in the ocean model surface layer is out-gassed according to a relaxation time scale, $ \mu$. Within the ocean interior, the tracer is passively advected by the ocean model currents. The full equation for the time evolution

$\displaystyle \frac{\partial C}{\partial t} \, = \, -U\cdot \nabla C \, - \, \mu C \, + \, \Gamma(C) \,+ \, S$ (3.110)

also includes a source term $ S$. This term represents interior sources of $ C$ such as would arise due to direct injection. The velocity term, $ U$, is the sum of the model Eulerian circulation and an eddy-induced velocity, the latter parameterized according to Gent/McWilliams (Gent and McWilliams [1990]; Gent et al. [1995]). The convection function, $ \Gamma$, mixes $ C$ vertically wherever the fluid is locally statically unstable.

The out-gassing time scale, $ \mu$, in eqn. (3.112) is set so that $ 1/\mu \sim 1 \ \mathrm{year} $ for the surface ocean and $ \mu=0$ elsewhere. With this value, eqn. (3.112) is valid as a prognostic equation for small perturbations in oceanic carbon concentrations. This configuration provides a powerful tool for examining the impact of large-scale ocean circulation on $ CO_2 $ out-gassing due to interior injections. As source we choose a constant in time injection of $ S = 1 \,\, {\rm mol / s}$.


3.19.1.2 Model configuration

The model configuration employed has a constant $ 4^\circ \times 4^\circ$ resolution horizontal grid and realistic geography and bathymetry. Twenty vertical layers are used with vertical spacing ranging from 50 m near the surface to 815 m at depth. Driven to steady-state by climatological wind-stress, heat and fresh-water forcing the model reproduces well known large-scale features of the ocean general circulation.


3.19.1.3 Out-gassing cost function

To quantify and understand out-gassing due to injections of $ C$ in eqn. (3.112), we define a cost function $ {\cal J} $ that measures the total amount of tracer out-gassed at each timestep:

$\displaystyle {\cal J}(t=T)=\int_{t=0}^{t=T}\int_{A} \mu C \, dA \, dt$ (3.111)

Equation(3.113) integrates the out-gassing term, $ \mu C$, from (3.112) over the entire ocean surface area, $ A$, and accumulates it up to time $ T$. Physically, $ {\cal J} $ can be thought of as representing the amount of $ CO_2 $ that our model predicts would be out-gassed following an injection at rate $ S$. The sensitivity of $ {\cal J} $ to the spatial location of $ S$, $ \frac{\partial {\cal J}}{\partial S}$, can be used to identify regions from which circulation would cause $ CO_2 $ to rapidly out-gas following injection and regions in which $ CO_2 $ injections would remain effectively sequestered within the ocean.

next up previous contents
Next: 3.19.2 Code configuration Up: 3.19 Sensitivity of Air-Sea Previous: 3.19 Sensitivity of Air-Sea   Contents
mitgcm-support@mitgcm.org
Copyright © 2006 Massachusetts Institute of Technology