Impact of hydrographic data assimilation on the Atlantic meridional overturning circulation

Impact of hydrographic data assimilation on the Atlantic meridional overturning circulation G. C. Smith, K. Haines, T. Kanzow, and S. Cunningham Env. Sys. Sci. Centre, Univ. Reading, Reading, UK National Oceanographic Centre, Southampton, UK now at: Environment Canada, Montreal, Canada now at: IFM-GEOMAR, Univ. Kiel, Kiel, Germany Received: 13 October 2009 – Accepted: 31 October 2009 – Published: 16 November 2009 Correspondence to: G.C. Smith (gregory.smith@ec.gc.ca) Published by Copernicus Publications on behalf of the European Geosciences Union.


Abstract
The poleward ocean heat transports in the North Atlantic controlled by the Atlantic Meridional Overturning Circulation (AMOC), play a key role in regional climate.If the AMOC can be initialized in numerical models through ocean assimilation this may help improve the predictability of North Atlantic climate variability on timescales out to a few years.Here we make an initial step toward the development of an ocean assimilation system that can determine the AMOC to support climate predictions.A detailed comparison is presented of 1 • and 1/4 • resolution global model simulations with and without sequential data assimilation to the observations and transport estimates from the RAPID/MOCHA mooring array across 26.5 • N in the Atlantic.Comparisons of modelled water properties with the observations from the merged RAPID boundary arrays demonstrate the ability of in situ data assimilation to accurately constrain the east-west density gradient between these mooring arrays.However, the presence of an unconstrained "western boundary wedge" between Abaco Island and the RAPID mooring site WB2 (16 km offshore) leads to the intensification of an erroneous southwards flow in this region when in situ data are assimilated.The result is an overly intense southward upper mid-ocean transport (0-1100 m) as compared to the estimates derived from the RAPID array.
Correction of upper layer zonal density gradients is found to compensate mostly for a weak subtropical gyre circulation in the free model run (i.e. with no assimilation).
Despite the important changes to the density structure and transports in the upper layer imposed by the assimilation, very little change is found in the amplitude and subseasonal variability of the AMOC.This shows that assimilation of upper layer density information projects mainly on the gyre circulation with little effect on the AMOC at 26 • N due to the absence of corrections to density gradients below 2000 m (the maximum depth of Argo).
The sensitivity to initial conditions was explored through two additional experiments using a climatological initial condition.These experiments showed that the weak bias Introduction

Conclusions References
Tables Figures

Introduction
The Atlantic Meridional Overturning Circulation (AMOC) involves a northward movement of warm upper waters accompanied by a southward movement of cold waters at depth.This circulation transports about 1 PW of heat northwards at 26 • N, which accounts for about half the northward energy transport by the global ocean at this latitude.This heat is released to the atmosphere in the North Atlantic and Arctic where it leads to higher temperatures over continents, notably in northern Europe.There is evidence that natural variations in the AMOC can affect climate over decadal timescales (e.g.Knight et al., 2005).Furthermore, forecasts presented by the Intergovernmental Panel on Climate Change (Meehl et al., 2007) suggest that the AMOC may slow down in response to greenhouse gas forcing, reducing the northward transport of heat by the Atlantic Ocean, and leading to a cooling of northern Europe that could offset anthropogenic warming.Evidence from palaeoclimate records suggests that the AMOC can undergo very rapid transitions, such as a total shutdown in little more than a decade (Jansen et al., 2007).It is possible that increasing levels of greenhouse gases could trigger such a rapid change, with potentially serious consequences for societies in Europe and other regions surrounding the Atlantic basin.
In the face of such risks, there is an obvious need for better, more quantitative, forecasts of the future behaviour of the AMOC: an AMOC prediction system could Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version Interactive Discussion potentially provide a valuable early warning of imminent climate change.There is evidence that changes in the AMOC, and related climate impacts, may be predictable on decadal timescales (Griffies and Bryan, 1997;Collins et al., 2003;Collins et al., 2006).This evidence is based on studies that assume near perfect knowledge of the initial ocean state and a perfect model.An important research question is whether potential predictability translates into actual prediction skill, when forecasts are initialised using real ocean observations.This is a much harder challenge because of limited observational data and problems arising from model bias.Nevertheless, recent studies (Smith et al., 2007;Keenlyside et al., 2008) have shown evidence that ocean initialization does lead to enhanced skill in decadal forecasts of global mean surface temperature, relative to reference forecasts that use changing radiative forcing but are not initialized.There is also evidence of particular skill in the North Atlantic (Pohlmann et al., 2009).These results are encouraging for the prospects of building a decadal forecasting capability.However uncertainty exists regarding the state of the ocean, which leads to uncertainty in the ocean initial conditions required for climate forecasts, resulting in larger forecast errors.The synthesis of ocean models and observations through data assimilation has been shown to constrain the ocean state estimate, which in turn leads to improved forecasting skill (Vidard et al., 2007).Corrections to upper layer temperature and salinity through assimilation should lead to improved density gradients and hence geostrophic currents and transports.Indeed, assimilation has been shown to improve the representation of the AMOC mean intensity and temporal variability (Balmaseda et al., 2007), as well as the agreement with published estimates from five transatlantic sections (Bryden et al., 2005).However the CLIVAR-GSOP (Global Synthesis and Observations Panel) intercomparison of the AMOC in various reanalysis products revealed large disagreements (K öhl, 2006).These differences were attributed to a combination of factors including variations in: ocean model configurations and resolutions, assimilation methods, and observational datasets used.An additional major source of uncertainty is related to the paucity of historical ocean observations that may only impose a weak constraint on the system.Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version

Interactive Discussion
Further advances in predicting the AMOC and its impacts on climate require better climate models, better observations of the ocean state, and improvements to ocean syntheses.In this context, the first direct continuous measurements of the AMOC, as provided by the RAPID array (Cunningham et al., 2007), are of particular importance.In the first place, these new measurements can be used to evaluate existing analyses and forecasts (Baehr et al., 2009), offering an unprecedented opportunity to assess the representation of the AMOC and related quantities.In the second place, the data can be assimilated to generate improved estimates of initial conditions for coupled models and perhaps improvements in forecast skill.In this paper we focus on the first of these steps, using a global ocean GCM at two different resolutions (1 • and 1/4 • ) to quantify the AMOC in these models and assess the impact of hydrographic data assimilation.The syntheses presented here assimilate all available in situ hydrography (e.g.Argo, XBT) apart from data from the RAPID array which is used for an independent evaluation.The assimilation of in situ temperature and salinity observations allow for corrections to the hydrography in mid-basin down to a maximum depth of 2000 m.The AMOC array is used to test the model circulations against these results and allows us to assess what additional information could be gained from assimilating RAPID array data.
Section 2 introduces the model, the data available for assimilation and the assimilation method used.Section 3 describes the RAPID array data and discusses how the RAPID transports are calculated, and also presents different ways of comparing these data with the model results.Section 4 presents the results of a number of different assimilation experiments designed to understand what aspects of the AMOC can be recovered from the assimilation experiments.Sections 5 and 6 look at the effect of initial conditions and resolution respectively, on the model AMOC.Section 7 provides a discussion and conclusions on how the RAPID array data can best be used with models.Introduction

Conclusions References
Tables Figures

Back Close
Full

Model description
The numerical model used here is the NEMO coupled ice-ocean model (Madec, 2008) version 2.3, based on the OPA9 ocean model (Madec et al., 1998) and the LIM2.0 sea ice model (Louvain sea Ice Model: Fichefet and Maqueda, 1997;Goosse and Fichefet, 1999).The ocean model is a primitive equation z-level model making use of the hydrostatic and Boussinesq approximations.The model employs a free surface (Roullet and Madec, 2000) with partial cell topography (Adcroft et al., 1997).The version used here has a tri-polar "ORCA" grid and 46 levels in the vertical, with thicknesses ranging from 6 m at the surface to 250 m at the ocean floor.
Two different model configurations have been used.The first has a global 1 • resolution with a tropical refinement to 1/3 • (ORCA1).The second is a global eddy-permitting model at 1/3 • resolution (ORCA025).Both configurations have been developed through the DRAKKAR Consortium (Barnier et al., 2007) and use model parameter settings as defined in Barnier et al. (2008) and Penduff et al. (2009).Both configurations employ an energy-enstrophy conserving momentum advection scheme (Barnier et al., 2008) and a Laplacian diffusion.Horizontal viscosity is parameterized with a biharmonic operator in ORCA025, and a Laplacian in ORCA1.Additionally, ORCA1 makes use of the Gent and McWilliams (1990) mixing parameterization.Vertical mixing is parameterized using a one-equation turbulent kinetic energy scheme (Blanke and Delecluse, 1991).
A detailed evaluation and comparison of these model configurations with satellite sea surface height observations is given by Penduff et al. (2009).Surface atmospheric forcing for the model is obtained from the DRAKKAR Forcing Set 3 (DFS3; Brodeau et al., 2009).DFS3 is a hybrid dataset making use of temperature, humidity and winds from the ERA40 atmospheric reanalysis  and the ECWMF operational analyses thereafter.Short and long wave radiative fluxes are obtained from the Common Ocean Reference Experiment (CORE) dataset (Large and Yeager, 2004) Project.Precipitation and snow fields have also been taken from CORE and have been subsequently modified to provide a more balanced global freshwater budget.

Assimilation scheme
Here, we use the S(T ) assimilation scheme (Smith and Haines, 2009), which is a sequential scheme for hydrographic data based on an Optimal Interpolation approach.
Temperature profiles (T ) are assimilated along with a salinity balancing increment as advocated by Troccoli and Haines (1999).Salinity profiles (S) are then assimilated along isotherms (i.e.S(T )) as advocated by Haines et al. (2006).By evaluating model errors as S(T ) rather than on depth levels (S(z)), the large dynamical variability present in S(z) is removed (i.e.heaving of water up and down), allowing a more accurate assessment of water mass errors.Second, S(T ) increments can be spread along isotherms so that corrections to a particular water mass do not influence adjacent water masses, which may have uncorrelated errors.Model-data differences (innovations) are evaluated at the closest model time step, and assimilation increments are calculated every 5 days (73 cycles per year).Increments are then introduced into the model evenly over the following day (known as Incremental Analysis Updating, see Bloom et al., 1996).Smith and Haines (2009) give more details of the assimilation and report on experiments all performed with the 1 • resolution ORCA1 model looking at the impact of Argo data in a 3 year reanalysis of the period 2002-2004.The only change to the assimilation system here from that presented in Smith and Haines (2009) is the addition of balancing increments for sea surface height and geostrophic velocity to reduce variability induced by the model adjustments.
In situ temperature and salinity observations are obtained from the UK MetOffice quality controlled ENACT/ENSEMBLES dataset (Ingleby and Huddleston, 2007).This dataset is largely composed of observations from the World Ocean Database 2005 (Boyer et al., 2006) and supplemented by data from the GTSPP (Global Temperature and Salinity Profile Program) and Argo autonomous profiling floats (Gould, 2005) such, the dataset includes all available hydrographic observations, including those from shipboard Expendable Bathythermographs (XBTs) and Conductivity-Temperature-Depth (CTD) measurements, as well as observations from mooring arrays (such as TAO and PIRATA).The operational quality control system used for the UK MetOffice FOAM assimilation products (Martin et al., 2007) is used to perform a number of consistency checks on the data.These include buddy checks, track checking, testing for density inversions and thinning of the data.A detailed description of the data included and the quality control processing is presented in Ingleby and Huddleston (2007).The version of the dataset used here (EN3 v1c) differs from previous versions of the dataset (e.g.EN2, EN3 v1a and EN3 v1b) in several ways (see www.hadobs.comfor details); in particular the Argo data from SOLO/FSI floats have been omitted due to errors in the reported pressure levels.

Experiment design
The simulations here focus on the period from January 2002 to April 2005, allowing comparison with the first complete year of the RAPID AMOC array results (deployed in April 2004).Initial conditions for most experiments with both the 1 • and 1/4 • resolution models were taken from multi-decadal control simulations.Control simulations begin in 1958 following a 10-year spin up from the World Ocean Atlas 2001 (WOA; Conkright et al., 2002) water mass properties, and employ surface forcing from DFS3 (Brodeau et al., 2009).The 1/4 • control run is designated CTL025 and was performed by the DRAKKAR consortium (Barnier et al., 2007) as run G70, while the 1 • control run we designate CTL, which has been referred to previously as ORCA1-R70 (Barnier et al., 2007;Penduff et al., 2009).
A description of experiments and their motivation is set out below (and summarized in Table 1), with all experiments referring to the 1 • resolution model except the last pair; the overall effect of assimilation in correcting errors was tested by starting assimilation in January 2002 (SYN) from the control run (CTL) conditions.Section 4 Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc

Printer-friendly Version
Interactive Discussion shows the improvements to the modelled temperature and salinity fields in the Atlantic provided by assimilation and the effect on transports at 26 • N as compared to estimates from the first year of the RAPID array (4 March-5 April).The initial response and adjustment of the AMOC is discussed in Sect. 5.
-CTL contains significant errors associated with drifts and biases developed over the 44 years of simulation prior to 2002.To investigate the influence of starting from better water masses two experiments (with and without assimilation) are made starting in January 2002 initialized from the WOA climatology and are presented in Sect. 5 (SYN-IC and CTL-IC respectively).
-To illustrate the sensitivity of the modelled transports near coastal boundaries, three 1-year experiments starting from CTL in January 2004 are presented, where the assimilation increments are tapered near coastal boundaries using a nearcoast tapering factor, c, calculated as Where the distance, d , of a model grid point to the nearest land point is less than the tapering length scale, L, both temperature and salinity increments are multiplied by this tapering factor prior to being applied to the model.Results are presented in Sect. 4 using length scales of 1000 km (SYN-T1000), 500 km (SYN-T500) and 0 km (SYN-T0, no tapering).
-Finally, the impact of resolution is studied using the 1/4 • resolution control There are a number of assumptions required to which the resulting transport estimates may be sensitive.In particular, the calculation of the barotropic circulation components, which cannot be monitored directly by density measurements, and the calculation of transports very close to the western boundary, make comparisons with modelled results, even for the higher resolution 1/4 • model, difficult.We therefore need to review the key measurements and the calculations made from them to put the model results in context.

AMOC calculation from the model
Here we describe the methods used to calculate the model equivalent of the RAPID array observations and transport estimates shown in C07.In both model resolutions the continental slope along the western boundary at 26.5 wall adjacent to several shallow grid points representing the Blake Plateau and Bahamas, with a well-defined Florida Strait to the west.
The "Mid-ocean transport" in C07 is calculated from zonal density gradients obtained by merging the temperature and salinity data from separate moorings into single profiles for the eastern and western boundaries.The mid-Atlantic ridge data is not used.
At the western boundary, the merged density profile contains data mainly from a tall mooring array, referred to as WB2, located roughly 16 km offshore (Fig. 2).Geostrophic transports obtained from these merged density profiles are combined with current meter transports covering a "western boundary wedge" (WBW) along the continental slope to give the mid-ocean transport.
To simulate this configuration, model currents from the first shallow grid point on the shelf directly adjacent to the continental slope and those westward to Abaco Island (Bahamas) are used to represent the WBW observed by current meters (Fig. 2).The hydrography in the first deep model column adjacent to the western continental slope is then compared to the RAPID array's merged western boundary profile.The eastern-most model grid point at each depth is used to compare to the RAPID array's merged eastern profile.These boundary profile comparisons are examined in Fig. 4 (temperature and salinity) and Fig. 5 (east and west boundary densities,ρ E and ρ W respectively).The separate contributions of the WBW and offshore transports are examined in Fig. 6.
Here we are interested in how changes in the cross-basin density gradient due to assimilation in the different model runs is manifest in meridional transports.Ignoring the mid-Atlantic ridge, the geostrophic meridional transport, V T W , can be obtained by integrating the thermal wind equation upwards from the deepest observed depth, D, as where L x is the distance between the western profile and the easternmost grid point, g is the gravitational acceleration, f the Coriolis parameter and ρ 0 is a reference density.Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version

Interactive Discussion
Figure 5 shows the transports obtained using (Eq.2) along with an assumption of no motion at the deepest observed level, i.e.V T W (D) = 0. Due to the omission of the deep flow of Antarctic bottom water, and therefore the level of no motion applied in this calculation, these transports do not represent the total meridional flow and differ from the transports of C07 (their Fig. 1).However a comparison of the water properties at the east and west boundaries in the model to those observed, will help us understand the extent to which assimilation constrains the zonal density gradient contribution to the mid-ocean transports.
In addition we examine the three AMOC components defined by C07: the windinduced surface Ekman flow, the Florida Current transport and the mid-ocean transport.For this purpose, we interpret the timeseries from C07 as "best estimates" of the various components (and AMOC), and thus compare to exact model values of these transports rather than estimates using the geostrophic methodology employed by C07 1 The Ekman volume transport, T E k , is obtained by zonal integral of the wind stress at 26.5 • N, τ(x), as where L B is the full basin width from the Bahamas to Morocco.As both model configurations have a well-defined Florida Strait bounded by the Bahama Islands and Florida, the volume transport in the Florida Current can be obtained by simply integrating the meridional velocities vertically and zonally across the strait.Finally, the mid-ocean transport is calculated here as the vertical and zonal integral of the meridional volume flux from the Bahamas to the eastern boundary, minus the Ekman volume transport.Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version

Interactive Discussion
The mid-ocean transport can be further decomposed in the vertical into transports in the thermocline recirculation (0-800 m), the intermediate layer (800-1100 m), the upper North Atlantic Deep Water (NADW; 1100-3000 m) and the lower NADW (3000-4730 m).Note that in this vertical decomposition only the thermocline layer transport (0-800 m) has the Ekman transport subtracted.These depth ranges were selected for comparison with the estimates of C07.The lower depth limit for the lower NADW corresponds to the deepest observed zonal density gradient from the RAPID array and omits a small northward transport of Antarctic Bottom Water.
Note that all calculations presented here have been repeated using model grid points slightly to the north and south of 26.5 • N, and no appreciable differences were found.In addition, several experiments were performed with the 1 • resolution model configuration to adjust the depth of the sills and islands around the Bahamas to more accurately reflect the actual bathymetry.While this changes some details of the local circulation, the meridional transports around 26 • N were not substantially affected.

Evaluation of model experiments at 26.5 • N
We first present the general hydrography of the control and assimilation experiments (Sect.4.1).Section 4.2 compares the model experiments to RAPID boundary array observations to assess the impact of hydrographic assimilation on mean geostrophicallyderived volume transports.Section 4.3 makes comparisons using the full meridional velocities from the model.Sections 4.4 and 4.5 look at the time evolution of the "midoceanic transports" and the AMOC components respectively, showing the main result found here that the gyre intensity can be greatly improved through assimilation and yet with little impact on the AMOC itself.The timescales of these adjustments and the impact of initial conditions are explored in Sect. 5. Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version Interactive Discussion

Basic evaluation of model experiments
A detailed evaluation of the S(T ) assimilation method and the impact of assimilation on the modelled temperature and salinity across the various global basins is provided by Smith and Haines (2009) and thus only a brief description is provided here.Assimilation is found to eliminate many biases present in the free-model run and leads to a significant improvement in the agreement with independent in situ hydrographic observations.Gemmell et al. (2008Gemmell et al. ( , 2009) ) show that use of the S(T ) assimilation method leads to a better representation of certain mode waters and water mass properties.
In addition, they show that this assimilation method provides an ocean state estimate which is in better agreement with in situ temperature and salinity observations than many other synthesis products.As such, the results presented here should be considered as being strongly constrained to data.
To illustrate the spatial structure of changes imposed by the assimilation, we show the 2004 annual-mean Atlantic temperature cross-section at 26.5 • N between Abaco Island (Bahamas) and Morocco for CTL and SYN as anomalies from the World Ocean Atlas 2001 climatology (Conkright et al., 2002) in Fig. 3.The CTL simulation exhibits a large-scale cold bias between 400-1000 m, associated with a weak subtropical gyre circulation.Smith and Haines (2009) demonstrate that this bias is mainly of a dynamical origin, and may be corrected through a vertical heave of water masses without requiring a modification to water mass properties (i.e S(T )). Figure 3  in the subpolar gyre, which results in an overly large penetration of Antarctic Bottom Water at depth (AABW).As there are no observations assimilated below 2000 m, the assimilation is unable to directly correct this model error locally.While some upstream correction may be possible (e.g.North of Denmark Strait or in the Labrador Sea), the SYN experiment is not long enough to see this effect.We show in Sect.4.4 that the presence of this error at depth has an important impact on the modelled volume transports at 26.5 • N and results in a weak AMOC.

Boundary profiles
Figures 4 and 5 shown the mean model properties at the western and eastern boundaries compared against the first-year RAPID array merged boundary profiles (averaged over April 2004-March 2005 corresponding to C07).Leaving aside the upper 300 m, the cold, fresh anomaly between 300-1000 m associated with the weak gyre circulation in CTL is seen at the western boundary (Fig. 4a and b) to represent errors of as much as 5 • C and 1 psu.These errors result in a positive density bias of nearly 0.4 kg/m 3 (Fig. 5a).The assimilation provides a correction to both temperature and salinity errors bringing experiment SYN into excellent agreement with both climatology and the RAPID array measurements (Fig. 4a and b).Indeed, the discrepancy between the model and the RAPID array at the western boundary is reduced to less than 0.1 kg/m 3 .
The correction of the western boundary densities also leads to a significant increase in the zonal density gradient in the upper layer (Fig. 5c), intensifying the subtropical gyre.
Between 1500-3500 m (Fig. 4, lower panels), there is a warm, saline error associated with weak NADW formation that can be seen on both sides of the basin in CTL, and which is only partially corrected in SYN.This bias in the NADW properties results in a negative density anomaly of up to 0.05 kg/m 3 (Fig. 5a).As there are no observations assimilated below 2000 m, there is little improvement of this error in SYN over CTL below 2000 m.
Below 4000 m the cold, fresh anomaly in CTL can be seen clearly in Fig. 4 and 5 at the western boundary.Although there is no direct assimilation below 2000 m, there is Introduction

Conclusions References
Tables Figures

Back Close
Full This results in a slight improvement in zonal density gradient below 3500 m.However, changes in density at the eastern boundary contribute to a poorer agreement with RAPID between 2500-3500 m (Fig. 5c).
The implication of these changes in the zonal density gradient for offshore geostrophic transports (i.e neglecting Florida Current and western boundary wedge contributions), can be inferred by application of the thermal wind equation (Eq.2) with a level of no motion applied at the deepest observed level by the RAPID Array.The resulting transport estimates are shown in Fig. 5d.One can see that the improvements at mid-depth result in a stronger southward flow for SYN over the upper 1000 m, correcting the weak flow in CTL and providing an excellent agreement with the equivalent RAPID array estimates.Interestingly, there is also some slight change in the deeper waters.

Vertical profile of meridional transports
The actual meridional transport between Abaco Island (the Bahamas) and Morocco depends on the full zonal pressure gradient evaluated on the two boundaries, which can differ from that determined using the difference between the merged boundary arrays, because the western array (mostly mooring WB2) is roughly 16 km offshore (Fig. 2).Thus, the analysis in the previous section omits a small "western boundary wedge" (WBW) between RAPID mooring WB2 and Abaco Island.The transport through this WBW is measured by the RAPID array using a series of moorings equipped with current meters.To more carefully assess the impact of assimilation on the mid-ocean transport, we separate the contributions from the WBW and the rest of the basin (i.e. from WB2 to the eastern boundary, or in the model from the continental shelf break to the eastern boundary) and compare transport estimates from the RAPID array with the full model velocities (Fig. 6).The mean offshore transports (i.e. the mid-ocean transport excluding the WBW) obtained using the full model velocities (Fig. 6a) show the main features noted above Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc

Printer-friendly Version
Interactive Discussion from the geostrophic estimates (Fig. 5d).In the thermocline recirculation (0-800 m), CTL has a notably weak mean southwards flow (−9.8 Sv), whereas the assimilation has intensified the gyre circulation in SYN (−19.0Sv), bringing it in close agreement with the RAPID estimates of C07 (−19.1 Sv).However, the mean transport in the WBW in SYN (−6.7 Sv) is poorly constrained2 , and leads to the intensification of an erroneous southwards flow3 in the upper 300 m present in CTL (−4.8 Sv), increasing the discrepancy with the current meter measurements from RAPID (+2.2 Sv) (Table 2).The increased transport in the WBW in SYN is related to an equivalent increase in the Florida Current transport (discussed later in Sect.4.5).As a result of this overly strong WBW transport, the mid-ocean transport in the thermocline recirculation in SYN is too large (−25.7 Sv).In the CTL the errors in the offshore and WBW transports tend to cancel providing a total transport (−14.6 Sv) in better agreement with those obtained from RAPID (−16.9Sv).These results are summarized in Table 2.
In contrast to the large changes seen in the upper layer, both CTL and SYN show very similar transports below 1500 m (Fig. 6, Table 3).While CTL and SYN have roughly the correct amplitude of upper NADW (1100-3000 m) transport (−12.7 Sv and −13.4 Sv respectively compared to −10.7 Sv for RAPID), both simulations have a significantly weaker southward flow of lower NADW (below 3000 m; −1.6 Sv and −1.8 Sv respectively) than found by C07 (−7.8 Sv).Indeed, a weak northward flow is present below 4000 m which is not present in the estimates of C07.Baehr et al. (2009)  a similar feature in their model simulation and attribute the difference to the omission of densities along the flank of the mid-Atlantic ridge in the zonal density gradient calculation from C07, leading to errors in their geostrophic flow estimates.An updated transport estimate from the RAPID array (Kanzow et al., 2009) that uses densities along the mid-Atlantic ridge confirms this result.However, even compared to the improved estimates of Kanzow et al. (2009;−5.1 Sv), both CTL and SYN show significantly weak transports of the lower NADW.Mean transports for CTL, SYN and RAPID for the various layers are summarized in Table 3.In summary, the assimilation provides a significant correction to the bias in the subtropical gyre resulting in an intensified offshore transport in better agreement with the RAPID array estimates of C07.However, an erroneous surface flow along the east coast of the Bahamas is not corrected through assimilation.These results suggest that to accurately recover the upper layer transports using assimilation, it is important that the assimilation constrains the densities close to boundaries as well, and not just in the deeper waters.This is of particular importance along the western boundary, where model errors may be uncorrelated with those offshore.Indeed, a recent study by Kanzow et al. (2009) has suggested that there is a significant anticorrelation of flows in the western boundary which compensate for offshore transports leading to a significant reduction in the variability in total meridional transport in the upper layer.

Time series of layer transports
We will now look at the time variability of transports.Following C07, the mid-ocean transports shown here include the western boundary wedge, but exclude the Florida Current.Also, as in Sect.4.3 (Fig. 6) we use the full model velocity fields to assess transports.
As noted in Sect.4.3, the assimilation induces a strengthening and overestimation of the southward flow in the thermocline recirculation (Fig. 7a, 0-800 m) due to errors in the WBW (Fig. 6b), despite an improvement in offshore transports (Fig. 6a).In addition to the change in amplitude, the variability in SYN is larger on sub-monthly to Introduction

Conclusions References
Tables Figures

Back Close
Full The impact of assimilation on the deeper waters is relatively minor.For the upper NADW (Fig. 7c, 1100-3000 m), both simulations have a weak sub-annual variability, which is not improved by the assimilation.The lower NADW transport (Fig. 7d, below 3000 m) also shows a weak sub-annual variability with little influence by assimilation.In addition, the RAPID observations show a strong intensification of southward flow in the spring and fall of 2004, which is not captured by the model simulations.However, both model simulations do show a southward pulse in early 2005 in agreement with the RAPID timeseries of C07.

Time series of AMOC components
We now turn our attention to the impact of data assimilation on the transport timeseries of AMOC components (Florida Current, Ekman and mid-ocean transports) as determined by the RAPID monitoring system (and described in C07). Figure 8 shows the three AMOC components along with the total AMOC (mean values are summarized in Table 4).The "mid-ocean" transport includes the full transport from Abaco Island to Morocco integrated down to 1100 m and is equivalent to the sum of transports shown in Fig. 7a and b.
As expected, the Ekman transports are nearly identical between C07 and the model experiments, with small instantaneous differences (1-2 Sv) occurring due to Introduction

Conclusions References
Tables Figures

Back Close
Full the different wind stress products used (QuikScat and ECMWF analyses respectively).Some additional minor differences are present between experiments due to the dependence of the wind stress formulation on the surface currents.Due to the erroneous southward WBW transport (Fig. 6b), the southward mid-ocean transport in SYN is too strong as compared to RAPID (Fig. 8d) by roughly 8 Sv.The intensification of the mid-ocean transport in SYN over that in CTL is accompanied by an increase in both the mean and variability of the northward flow of the Florida Current (Fig. 8a).In particular, significant sub-annual variability of the Florida Current transport is excited by the assimilation.The mean amplitude of the Florida Current transport in SYN (33.3 Sv) is in better agreement with observed cable measurements (31.7 Sv) as compared to CTL (23.3 Sv), although the variations are not correlated.
The large changes in the amplitude and variability of the mid-ocean and Florida Current transports in SYN show a nearly perfect compensation, such that there is remarkably little change to the overturning circulation (Fig. 8b).Indeed, the mean amplitude of the overturning in SYN (11.2 Sv) is only 0.3 Sv larger than that of CTL (10.9 Sv), despite a 10.0 Sv increase in the mean amplitude of the Florida Current transport.Therefore the changes imposed by the assimilation to the upper layer density appear to project almost entirely onto the subtropical gyre circulation, with little influence on the overturning.Moreover, because the assimilation does not adequately constrain the transports in the WBW, sub-monthly variations in the mid-ocean and Florida Current transports are not captured either.
To help understand the compensation between the mid-ocean and Florida Current transports and why the assimilation appears to have little impact on the AMOC, three additional experiments were performed.As the most significant changes from the assimilation occur along the western boundary, in these experiments the assimilation increments are reduced near coastal boundaries to reduce the impact of assimilation on the zonal density gradients.The experiments SYN-T1000, SYN-T500 and SYN-T0 are made using Eq. ( 1) with length scales of 1000 km, 500 km and 0 km respectively.In addition, these simulations were initalized from CTL in January 2004 to allow a direct

Conclusions References
Tables Figures

Back Close
Full Figure 9 shows the resulting AMOC and components from experiments CTL, SYN-T1000, SYN-T500 and SYN-T0.Reducing the impact of assimilation near the coastal boundaries has little effect on the overturning circulation strength.However it has a much larger effect on the gyre variability as seen in the Florida Current and midocean transports.Thus, reduced assimilation near the boundary reduces the aliasing of offshore variability (e.g.eddies and Rossby waves) onto the total gyre transports.

Sensitivity to initial conditions
The assimilation experiments discussed above have all been initialized from the spunup control run and thus begin with large water mass errors, which have developed in CTL since the initialization in 1958.Here we seek to examine the timescales on which the AMOC responds in experiments that begin with more correct density gradients that we then seek to maintain through assimilation.A better knowledge of the timescales of AMOC adjustment will highlight the physical processes involved in the adjustments and shed light on how best to maintain a balanced model circulation.
Two additional simulations are considered: CTL-IC uses no assimilation and is initialized from the WOA climatology in January 2002 and SYN-IC is the same as CTL-IC except with data assimilation.Timeseries of the Florida Current transport, the upper mid-ocean transport and the AMOC for these four simulations are shown in Fig. 10 for the period January 2002-January 2003.After a fast initial spin-up of the circulation in CTL-IC taking only a few days, the AMOC variability (Fig. 10e) is nearly identical to that of CTL but differs by being persistently higher by about 2 Sv bringing the results closer to those of RAPID.Also CTL-IC initially has a much stronger subtropical gyre circulation, which is reflected in the stronger southward upper mid-ocean transport (Fig. 10c) and northward Florida Current transport (Fig. 10a).Over a period of about 6 months,

Conclusions References
Tables Figures

Back Close
Full  3) than CTL (1.6 Sv).Moreover, these differences are still detectable after 3 years (not shown).
A comparison of the transports for SYN and SYN-IC is shown in the right panels of Figs. 10 and 11.Over the first month the temperature and salinity fields for both SYN and SYN-IC undergo significant modification and converge rapidly, resulting in very similar upper mid-ocean transports (Fig. 10d).The resulting transports have a similar strength to the initial transports in CTL-IC, with values considerably high compared to observations (see Table 3).The Florida Current transport in SYN-IC reaches a high value over 30 Sv immediately following initialization as in CTL-IC, but remains high instead of dropping back on a 6 month timescale as in CTL-IC.In constrast, it takes about 6 months for the Florida Current transport in SYN to get above 30 Sv (Fig. 10b).During this period of Florida Current adjustment in SYN, there is a strong southward component to the gyre circulation, and a weak Florida Current transport which together result in a very weak AMOC of about 5 Sv (Fig. 10f).During this adjustment the stronger southward upper mid-ocean flow is compensated by a northward barotropic flow (mostly along the western boundary) which affects the deeper transports seen in Fig. 11d and f.Note that the transports for SYN-IC and CTL-IC for the deep layers below 1100 m have very similar amplitudes over most of the first year because there is little correction to the zonal density gradients in SYN-IC at these depths.

Conclusions References
Tables Figures

Back Close
Full To assess the robustness of the AMOC response to model resolution we have two model simulations of the 1/4 • resolution global NEMO configuration: CTL025 (no assimilation) and SYN025 (with assimilation).The AMOC and its three components (Florida Current transport, Ekman transport, mid-ocean transport) for CTL025 and SYN025 are shown in Fig. 12. CTL025 shows similar areas of agreement and disagreement with the estimates of C07 as the 1 • resolution model (which is why the previous experiments have been presented at the lower resolution).For instance, CTL025 has a weak Florida Current transport (Fig. 12a) and a mid-ocean transport broadly consistent with C07 (Fig. 12d), which together result in a weak AMOC (Fig. 12b).
The assimilation again reduces the errors in temperature and salinity (not shown) resulting in a stronger gyre circulation, a more intense Florida Current transport (Fig. 12a), and a stronger southward mid-ocean circulation (Fig. 12d).As with the 1 • resolution model, in SYN025 the more intense gyre circulation is accompanied by a stronger southwards circulation in the WBW, creating a disagreement with the estimates of C07.Despite this significant modification to the gyre caused by assimilation, only a relatively small change can be seen in the AMOC, present from April to August.This difference in AMOC between CTL025 and SYN025 can be mostly attributed to changes between 1100-3000 m associated with a pulse of southward flowing upper NADW over the same period (Fig. 13c).As such, the improvements in AMOC in in deep water formation at higher latitude will have had time to influence the AMOC at 26.5 • N.

Discussion
Clearly, there are still many challenges in constraining the meridional overturning circulation in ocean or climate models.We have demonstrated that assimilation of upper layer hydrographic data, even with the Argo dataset, is not sufficient to constrain the Atlantic AMOC at 26.5 • N. A major reason is because it fails to constrain the deep density gradients below 2000 m.A consequence of this is that circulation changes resulting from assimilation are reflected in gyre transports rather than AMOC transports.To assess the impact of modifying the deeper water masses, we did perform some preliminary experiments assimilating the RAPID array data itself into the model (not shown).However, when the deep RAPID array observations are assimilated using standard Gaussian covariance scales, large changes in dynamic height can develop around the RAPID moorings resulting in a deterioration of the local circulation.A more careful approach may thus be needed, for example, involving boundary specific covariances associated with lower NADW variability and bias.One of the key aims for developing an ocean assimilation system that can determine the AMOC is to support climate predictions.The poleward ocean heat transports controlled by the AMOC play a key role in controlling climate in and downstream of the North Atlantic.If the AMOC can be initialised through ocean assimilation this may help predict some aspects of North Atlantic climate variability on timescales out to a few years.The results presented here (in Sect.5), show that an improved representation of deep water masses in the initial conditions (i.e.CTL vs. CTL-IC) can persist for periods longer than 3 years and result in an improved AMOC at 26.5 • N.However, simply correcting the upper layer zonal density gradient does not result in a better AMOC.
Moreover, when the assimilation disturbs the sensitive balance in the gyre circulation between the mid-ocean and the Florida Current a significant deterioration of the AMOC Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version Interactive Discussion can occur.Even if a particular gyre intensity in the initial condition is not supported by the model and a bias develops (as found here), as long as the balance between the mid-ocean and the Florida Current transport is maintained, then the AMOC is unaffected.The results here, particularly on the persistence of the better AMOC with improved initial lower NADW distributions, do suggest that AMOC predictions of several years may be possible provided the deeper waters are correctly initialized and a balanced gyre is maintained.This is promising for the development of decadal forecasting capabilities.
These results also suggest that while the RAPID/MOCHA monitoring array at 26.5 • N provides valuable information about the state of the overturning and heat transport in the North Atlantic, it may not be the ideal observing system by itself to aid in initializing coupled model AMOC predictions.An observing system providing additional information about the state of dense waters upstream in the subpolar gyre could be of great benefit.For example, this could be achieved through the occasional profiling to 4-5 km depth by a small subset of Argo floats.However, the results presented here do show some evidence of the benefits of assimilation of Argo at higher latitudes on the AMOC at 26.5 • N, consistent with assertions made by previous studies (e.g.Balmaseda et al., 2007).Further work using long reanalyses are required to evaluate the extent to which assimilation can improve dense water formation, and thus improve the AMOC.Finally, it has been suggested that sequential data assimilation may induce gravity wave noise which could induce spurious mixing that would negatively impact on the model AMOC.Here, we obtain only a small difference in AMOC with assimilation (i.e. between CTL and SYN), despite a radical impact on the gyre circulation.This is evidence that the assimilation is not adversely affecting the water masses through spurious mixing.

Conclusions
Here, we have presented a detailed comparison of model simulations with and without sequential data assimilation to the observations and transport estimates from the RAPID/MOCHA mooring array across 26.5 • N in the Atlantic.Comparisons of modelled water properties with the observations from the merged RAPID boundary arrays demonstrate the ability of in situ data assimilation to tightly constrain the east-west density gradient between these mooring arrays.However, the presence of an unconstrained "western boundary wedge" between Abaco Island (the Bahamas) and the RAPID mooring site WB2 (16 km offshore) leads to the intensification of an erroneous southwards flow in this region when in situ data are assimilated.The result is an overly intense southward upper mid-ocean transport (0-1100 m) as compared to the estimates derived from the RAPID array (C07).In addition, the sub-seasonal variability of both the mid-ocean and Florida Current transports is not well recovered by the assimilation.This highlights the importance of processes near the western boundary in setting the zonal pressure gradient (and the AMOC), rather than the zonal density gradient as observed close to, but not at, the western boundary.
Correction of upper layer zonal density gradients is found to compensate mostly for a weak subtropical gyre circulation in the free model run (i.e. with no assimilation).The rapid intensification of the zonal density gradient (and southwards thermal wind), leads to a fast barotropic response with a northward flow focused near the western boundary.The northward barotropic flow results in a reduction in southwards transport of the upper and lower NADW and hence a slowing of the AMOC.Gradually, over a period of 6 months the Florida Current transport intensifies in response to the stronger zonal density gradient allowing the northward barotropic circulation to slacken and the AMOC to return to its previous amplitude.
Despite the important changes to the density structure and transports in the upper layer imposed by the assimilation, very little change occurs in the AMOC after the initial 6 month adjustment is made to the gyre amplitude.This shows that assimilation Introduction

Conclusions References
Tables Figures

Back Close
Full of upper layer density information is absorbed mainly by the gyre circulation with little effect on the AMOC due to the absence of corrections to density gradients below 2000 m (the maximum depth of Argo).An assimilation experiment in a higher resolution eddy-permitting 1/4 • resolution model confirmed that the primary response to upper layer hydrographic assimilation is in the strength of the gyre and not the AMOC.
This response was also found to be robust to changes in model bathymetry near the western boundary (not shown), and details of sea surface height and geostrophic velocity balancing increments applied by the assimilation (not shown).
The sensitivity to initial conditions was explored through use of two additional experiments using a climatological initial condition.These experiments showed that the weak bias in gyre intensity in the control simulation (without data assimilation) develops over the period of about 6 months, but does so independently from the overturning, with no change to AMOC.However, differences in the properties and volume of NADW persisted throughout the 3 year simulations resulting in a difference of several Sverdrups in AMOC intensity.The persistence of these dense water anomalies and their influence on the AMOC is promising for the development of decadal forecasting capabilities.However, the results found here suggest that the deeper waters must be accurately initialized in order to constrain the AMOC.This highlights the need for deep ocean profiling floats which can be used to constrain the circulation and properties of the dense waters present at depth.Further study is required using longer simulations to fully understand the extent to which the deeper waters at 26.5 • N are constrained by assimilation at higher latitudes.Introduction

Conclusions References
Tables Figures

(3
CTL025) with a 1/4 • assimilation experiment initialized from the control in January 1987 (SYN025).A comparison with the first year of RAPID array transport estimates(March 2004-April 2005)  is presented in Sect.6. Description of RAPID array and transport calculationsThe RAPID/MOCHA monitoring array was deployed across 26.5 • N in the Atlantic ocean starting in April 2004 and is presently funded until at least 2013.The array consists of moorings on the western boundary of the Atlantic, on either side of the mid-Atlantic ridge, and on the eastern boundary slope (Fig.1).The moorings provide continuous measurements of the hydrography in these areas and are supplemented by current meter measurements at the western boundary (seeKanzow et al., 2008 for a detailed description of the mooring array and instrumentation).The transports calculated from this array can be combined with the cable data monitoring of the Florida Current transport between Florida and the Bahamas(Baringer and Larson, 2001), and estimates of Ekman transport from wind data, to give a total AMOC at 26.5• N in the Atlantic.The mooring data are of course only available after moorings have been recovered, normally on an annual timescale.When the current study was performed, observations from the first year of deployment only were available.Cunningham et al. (2007) (C07 hereafter) andKanzow et al. (2007) describe how meridional transports and the AMOC can be calculated from the array and other data.
shows clearly that this bias has been corrected by assimilation and is no longer present in SYN.Since the deepest observations in this region are from Argo, which only provide observations over the upper 2000 m, the direct influence of temperature and salinity assimilation does not penetrate below this depth.As such, model biases below 2000 m present in CTL can only be affected indirectly; for example, by improved transports or the advection of corrected waters below 2000 m.In particular, a cold anomaly found in CTL below 3000 m on the western side of the mid-Atlantic ridge can be seen to persist in SYN.This anomaly is likely associated with insufficient formation of dense NADW Introduction of the isopycnals along the western boundary below 3000 m in SYN.
and probably reflects aliasing of mesoscale and Rossby wave variability from offshore onto the western boundary densities.The sub-monthly fluctuations in both model runs show little correlation with those from the RAPID array where variability is mostly associated with waves close to the western boundary.On longer timescales, both model runs display a minimum in southward thermocline flow in October in agreement with the RAPID array observations, and longer model simulations show that this minimum is a robust feature of the model seasonal cycle (not shown).In the intermediate waters (Fig.7b, 800-1000 m), both CTL and SYN have roughly the same amplitude of flow as found by C07, with SYN again showing a slightly larger variability.
. Initialization was done in 2004 instead of 2002 to reduce computational costs.As such, SYN-T0 and SYN differ only in that they have different start dates (January 2004 and January 2002 respectively).
gradients relax until the Florida Current transport and the upper midocean transport both reach values within 1 Sv of the CTL simulation, with the Florida Current transport being 1 Sv stronger and the upper mid-ocean return flow being 1 Sv weaker, providing for the extra 2 Sv of overturning.This adjustment of the gyre circulation is independent of the AMOC and is a further example of the decoupling of the gyre and AMOC components discussed in Sect.4.5.Figure11(left panels) show the mid-ocean transport breakdown by depth layers.The persistently stronger southward flow of lower NADW below 3000 m and weaker southward flow of upper NADW (1100-3000 m) for CTL-IC are clear.In particular, the stronger lower NADW transport at the end of 2002 in CTL-IC of about 3 Sv is in better agreement with the RAPID mean value (7.8 Sv; Table SYN025 over CTL025 seen in the summer of 2004 are most likely associated with upstream improvements in the model state.Hence, as with the 1 • resolution model, the 1/4• resolution model shows a distinct separation of the gyre and overturning circulations in the Atlantic, with the direct effect of upper ocean hydrography assimilation impacting mostly the gyre circulation only, with only an indirect effect present in the deeper layers over longer timescales.Note that some differences between CTL025 and SYN025 (such as the pulse of upper NADW) may be expected because SYN025 was initialized from CTL025 in 1987 (rather than in 2002 as in SYN) and thus changes

Fig. 1 .
Fig. 1.Distribution of RAPID/MOCHA moorings along 26.5 • N in the subtropical North Atlantic.Mooring locations are shown in the upper panel along with a cross-section illustrating the mooring structure in the lower panel.The current meter moorings west of WB2 are not shown here for clarity.

Fig. 2 .Fig. 5 .Fig. 7 .Fig. 13 .
Fig. 2. RAPID/MOCHA moorings near the western boundary (off Abaco, the Bahamas).Density sensors, bottom pressure recorders and current meters are denoted as crosses, squares and circles, respectively.The dots at WBA and WB0 indicate the part of the water column covered by Acoustic Doppler Current Profiler measurements.
, which are derived from the International Satellite Cloud Climatology 2672 Introduction

Table 3 .
Annual-mean meridional layer transports at 26.5• N (April 2004-March 2005).Transports are split into four depth ranges, corresponding to the thermocline recirculation (0-800 m), the intermediate water (800-1100 m), the upper NADW (1100-3000 m) and the lower NADW (below 3000 m).Model values are calculated by integrating the meridional velocities across the full basin (from Abaco to Morocco).The RAPID values are quoted from C07.