Stiftung Alfred-Wegener-Institut f ¨ur Polar- und Meeresforschung
State estimation in support of an iron fertilization experiment in the Antarctic Polar Frontal Zone.
Martin Losch, Volker Strass, and Boris Cisewski
Alfred-Wegener-Institut f ¨ur Polar- und Meeresforschung, Bremerhaven, Germany, email: Martin.Losch@awi.de
1 Overview
The mutual availability of a high-resoltion data set and a suitable numerical ocean model makes a multivariate data-model analysis possible. In general, combining models and data with data assimilation or state estimation techniques is promising when both data and models separately ex- hibit skill. However, in oceanography more often than not, data and in particular sub- surface data are sparse and the prediction skill of ocean models tends to be poor on long time scales. We present a state esti- mation experiment, in which we use high- resolution hydrographic, tracer and veloc- ity data from the European Iron Fertiliza- tion EXperiment (EIFEX) to constrain a high-resolution coupled ecosystem-ocean circulation model of the experimental site in Atlantic sector of the Antarctic Polar Frontal Zone.
2 EIFEX: European Iron Fertilization EXperi- ment
EIFEX was aimed at testing the hypothe- sis that iron limits phytoplankton blooms in the Southern Ocean. For the open ocean experiment, a patch with a diame- ter of 15 km inside of a cyclonic, mesoscale eddy in the polar frontal zone was fertilized on February 12–13 and February 26–27, 2004 with dissolved iron. Subsequently the oceanic response was monitored. The eddy was identified with the help of in-situ mea- surements (CTD sensor and ship mounted ADCP) and satellite altimetry. It extended over 85 km by 120 km. with the center near 49◦24’S 02◦15’E in the South Atlantic
During the course of the experiment both hydrographic and dynamic parameters and bio-geochemical quantities were measured at CTD stations inside and outside the fer- tilized patch and along the ship track. Air- borne LIDAR-data were use to track the lo- cation of the bloom.
3 MITgcm and state estimation
The M.I.T. general circulation model (MITgcm Marshall et al., 1997) has been adapted for use with the Tangent linear and Adjoint Model Compiler (TAMC), and its successor TAF (Transformation of Algorithms in Fortran, Giering and Kaminski, 1998). Efficient (w.r.t. CPU/memory), exact (w.r.t. the model’s tran- sient state) derivative code can be generated for up-to-date versions of the MIT- gcm and its newly developed packages in a wide range of configurations (He- imbach et al., 2002, 2005). Here, the MITgcm is configured to cover the exper- imental region of the EIFEX cruise of appoximately 150 km by 194 km with a mean horizontal grid spacing of approximately 3.6 km; vertical grid spacing in- creases from 10 m near the surface to 25 m at 500 m depth. The integration time spans the length of the experiment (39 days).
0 20 40 60 80 100 120 140 160
10−10 10−9 10−8 10−7 10−6 10−5 10−4 10−3 10−2 10−1 100
iteration
total T/S−data U/V−data
T/S initial conditions buoyancy flux wind stress τ open boundaries
Figure 1: Evolution of different contribu- tions to the objective function J with iteration number of the BFGS-descend algorithm. Hy- drographic data from CTD-stations and ship- mounted ADCP-current measurements are used to assimilate the model using the variational data assimilation technique (state estimation, “4D- VAR”). During the minimization of the objec- tive function J that describes the quadratic de- viation of the model from data and also penal- izes deviations from the first guess, initial con- ditions for temperature and salinity, open bound- ary conditions for temperature, salinity, and ve- locity, and surface fluxes of heat, fresh water, and momentum are adjusted to give the best fit to observations.
data
40’ 2oE 20’ 40’ 3oE 20’
20’
50oS
40’
20’
49oS
first guess solution
40’ 2oE 20’ 40’ 3oE 20’
20’
50oS
40’
20’
49oS
salinity
33.75 33.8 33.85 33.9 33.95
optimized solution
40’ 2oE 20’ 40’ 3oE 20’
20’
50oS
40’
20’
49oS
Figure 2: Comparison with in-situ data shows that data assimilation improves the position of the eddy by adjusting open boundary conditions, initial hydrography, and (to a lesser extend) surface flux boundary conditions. Left:
surface salinity from hydrographic measurements, average over first 10 days of the experiment. Center: surface salinity of modeled eddy on day 5 without data assimilation; initial conditions and boundary conditions for this first guess are obtained by interpolating and extrapolating all available data; the amount of data is not sufficient to allow for time dependent boundary conditions, which appears to be the major problem for this solution. Right: surface salinity of modeled eddy on day 5 after full time-dependent state estimation; the eddy has moved southward and away from the boundary and its position is now in much better agreement with observations (Left).Hind-Cast Mesoscale Altimetry - Mar 20, 2004
1˚E 1˚E
2˚E 2˚E
3˚E 3˚E
4˚E 4˚E
51˚S 51˚S
50˚S 50˚S
49˚S 49˚S
48˚S 48˚S
-30 -25 -20 -15 -10 -5 0 5 10 15 20 25 30 Sea Surface Height Anomaly (cm)
Sea Surface Height [cm]
−30−25−20−15−10 −5 0 5 10 15 20 25 30 first guess SSH
1oE 2oE 3oE 4oE 51oS
50oS 49oS 48oS
Sea Surface Height [cm]
−30−25−20−15−10 −5 0 5 10 15 20 25 30 optimized SSH
1oE 2oE 3oE 4oE 51oS
50oS 49oS 48oS
Figure 3: Left: sea surface height (SSH) anomaly from satellite altimetry on day 37 of the experiment (not used in the assimila- tion; source: http://www-ccar.colorado.edu/˜realtime/gsfc_
global-real-time_ssh). Center: SSH of modeled eddy on day 37 with- out data assimilation. Right: SSH of modeled eddy on day 37 after full time- dependent state estimation; the eddy extends further south than before, but the comparison of absolute SSH with height anomalies is ambiguous.
measurements
5 10 15 20 25 30 35 40
−150
−100
−50 0
optimized solution
5 10 15 20 25 30 35 40
−150
−100
−50 0
temperature [degC]
3.8 4 4.2 4.4
temperature trend (warming) in the mixed layer
first guess solution
day of experiment
5 10 15 20 25 30 35 40
−150
−100
−50 0
measurements
z [m]
5 10 15 20 25 30 35 40
−150
−100
−50 0
salinity
33.87 33.875 33.88 33.885 33.89 33.895 33.9
first guess solution
z [m]
day of experiment
5 10 15 20 25 30 35 40
−150
−100
−50 0
salinity trend (freshening) in the mixed layer
optimized solution
z [m]
5 10 15 20 25 30 35 40
−150
−100
−50 0
Figure 4: Horizontally averaged temperature (left) and salinity (right) in the fertilized patch; data (top panels), optimized model (middle panels), first guess model (bottom panels). The optimization improves the description of the mixed layer depth dramatically, so that the optimized solution describes the warming and freshening of the mixed layer more accurately than the first guess solution.
4 Ecosystem Model REcoM: Comparison to Observations
In our study we use a newly developed regulated ecosystem model (REcoM, Schartau et al., 2007) that is based on an approach of Geider et al. (1998). It explicitly decouples carbon and nitrogen fluxes and does not rely on a fixed Redfield ratio. For Southern Ocean applications, REcoM has been extended to account for diatom blooms, opal export, and iron explicitly (Hohn et al., 2007).
Four additional state variables have been added: silicate, iron, and biogenic silicate compounds in phytoplankton and detritus.
This model is coupled to the physical circulation. Initial conditions for the 16 state variables of REcoM are estimated from observations and 1D-sensitivity and tuning experiments (not shown).
40’ 2oE 20’ 40’ 3oE 20’
20’
50oS
40’
20’
49oS
2
2 4
4 6 4 26−Feb−2004
0 0.5 1 1.5 2 2.5 3
40’ 2oE 20’ 40’ 3oE 20’
20’
50oS
40’
20’
49oS
2 2
2 22 4
6
01−Mar−2004
0 0.5 1 1.5 2 2.5 3
40’
2oE 20’ 40’ 3oE 20’
20’
50oS
40’
20’
49oS
2
2
2 2
4 4
6 10−Mar−2004
0 0.5 1 1.5 2 2.5 3
40’
2oE 20’ 40’ 3oE 20’
20’
50oS
40’
20’
49oS
2
2 2
2 2 2
2
2
4 4
4 4
6 6 6 18−Mar−2004
0 0.5 1 1.5 2 2.5 3
Figure 5: Model surface chlorophyll concentration (in mg m−3) on selected days (15, 19, 28, and 36 days after fertil- ization). Overlaid are contours of normalized LIDAR-derived fluourescence giving an impression of the observed bloom lo- cation.
0 10 20 30
0.05 0.1 0.15 0.2 0.25
∫(Chl a) dz [g/m2], top 100 m
model inpatch model outpatch obs inpatch obs outpatch
0 10 20 30
1 1.5 2 2.5 3
∫PON dz [g/m2], top 100 m
0 10 20 30
4 6 8 10 12 14 16
∫POC dz [g/m2], top 100 m
0 10 20 30
2550 2555 2560 2565 2570
∫DIC dz [g/m2], top 100 m
Figure 6: Integral of top 100 m of observed and modeled chlorophyll a (gChl m−2) and POC (gC m−2) concentrations in the fertilized patch.
Both particulate organic matter (POM) and the uptake of dis- solved inorganic carbon (DIC) appear to be slightly underesti- mated during the peak of the observed bloom, consistent with the too early maximum of chlorophyll concentration in the model. Note that POM outside the patch drops immediately after the beginning of the intergration contrary to observations.
This drop of approximately 4 gC m−2 is a consequence of the constant sinking velocity imposed in REcoM. DIC outside the patch increases in contrast to observations.
5 Results: Nutrient Budgets and Flux Estimates from the Coupled Ecosytem Circulation Model
vertical integral of iron induced nutrient difference [g/m2], day 39
40’
2oE 20’ 40’
3oE 20’
20’
50oS
40’
20’
49oS
carbon
−20 −10 0
40’
2oE 20’ 40’
3oE 20’
20’
50oS
40’
20’
49oS
nitrogen
−4 −2 0
40’
2oE 20’ 40’
3oE 20’
20’
50oS
40’
20’
49oS
silicate
−8 −6 −4 −2 0
Figure 7: Estimated carbon, nitrogen, and silica consumption through bio- logical activity induced by iron fer- tilization: difference (deficit) of nu- trients at the end of the integration for experiments with and without iron fertilization. From carbon invento- ries in the run with iron fertilizta- tion, we esimate a total change of DIC of 10–15 gC m−2 over 38 days (Fig. 6), which falls within the range of estimates by Smetacek et al. (2008) when surface gas exchange is ex- cluded from the calculation. The DIC difference between runs with and without iron fertilization, integrated to 100 m depth, peaks at 16.8 gC m−2. The value increases to 23.6 gC m−2 if integrated to 500 m (bottom of the do- main, see figure).
−0.8 −0.6 −0.4 −0.2 0
−500
−450
−400
−350
−300
−250
−200
−150
−100
−50 0
inside fertilized patch
gC m−2 d−1
−0.8 −0.6 −0.4 −0.2 0−500
−450
−400
−350
−300
−250
−200
−150
−100
−50 0 outside fertilized patch
gC m−2 d−1
z [m]
no sinking vDet = 5m/d vDet = 10m/d vDet = 20m/d vDet = 50m/d vDet = 100m/d vDet = 500m/d vDet = 1000m/d
Figure 8: Vertical flux of particulate organic carbon (POC, in mmol C m−2 s−1) averaged over the fertilized patch and the last 10 days of the experiment. Here we estimate the effect of the sinking velocity. Changing this constant parameter has a large effect for small values (below 100 m/d) but becomes less important for large values.
Conclusions: This study provided the rare opportunity to investigate the oceanic response to an iron fertilisation experiment over a period of about 38 days within a mesoscale eddy at the Southern Polar Front. In a first step we showed that the 3D global state estimate from hydrography and velocity data are in good agreement to the observa- tions. The result form the ecosystem model REcoM are promising, because the modeled chlorophyll concentrations coincide with the LIDAR-derived fluorescence data and the derived nutrient budgets and flux estimates are on the same order as direct estimates from observations.
References
Geider, R. J., MacIntyre, H. L., and Kana, T. M.
(1998). A dynamic regulatory model of phyto- planktonic acclimation to light, nutrients, and temperature. Limnol. Oceanogr., 43(4):679–
694.
Giering, R. and Kaminski, T. (1998). Recipes for adjoint code construction. ACM Trans. Math.
Softw., 24(4):437–474.
Heimbach, P., Hill, C., and Giering, R. (2002).
Automatic generation of efficient adjoint code for a parallel Navier-Stokes solver. In J.J.
Dongarra, P.M.A. Sloot and C.J.K. Tan, edi-
tor, Computational Science – ICCS 2002, vol- ume 2331, part 3 of Lecture Notes in Computer Science, pages 1019–1028. Springer-Verlag, Berlin (Germany).
Heimbach, P., Hill, C., and Giering, R. (2005). An efficient exact adjoint of the parallel MIT gen- eral circulation model, generated via automatic
differentiation. Future Generation Computer Systems (FGCS), 21(8):1356–1371.
Hohn, S., V¨olker, C., and Gladrow, W.-D. (2007).
Is strong silicification of diatoms in the South- ern Ocean a physiological response to iron lim- itation? submitted.
Marshall, J., Hill, C., Perelman, L., and Adcroft, A. (1997). Hydrostatic, quasi-hydrostatic, and nonhydrostatic ocean modeling. J. Geophys.
Res., 102(C3):5733–5752.
Schartau, M., Engel, A., Schr¨oter, J., Thoms, S., V¨olker, C., and Wolf-Gladrow, D. (2007).
Modelling carbon overconsumption and the
formation of extracellular particulate organic carbon. Biogeosciences, 4:433–454.
Smetacek, V., Strass, V. H., Klaas, C., Assmy, P., et al. (2008). Massive carbon flux to the deep sea from an iron-fertilized phytoplankton bloom in the Southern Ocean. submitted.