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.95)

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.96) is set so that $ 1/\mu \sim 1 \ \mathrm{year} $ for the surface ocean and $ \mu=0$ elsewhere. With this value, eqn. (3.96) 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.96), 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.96)

Equation(3.97) integrates the out-gassing term, $ \mu C$ , from (3.96) 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 Last update 2011-01-09