1. Introduction
The vertical mixing coefficient D controls the diapycnal component of the meridional overturning circulation (Munk and Wunsch 1998). The lack of knowledge of vertical mixing, its magnitude and spatial variation, leads to large uncertainties in our estimates of the volume, heat, freshwater, nutrient, and tracer transport and the role of the ocean in climate sensitivity and climate feedbacks (Manabe and Stouffer 1988). The along-isopycnal mixing coefficient K is used for mixing tracers such as temperature and salinity or potential vorticity along isopycnal surfaces. The value of K, like D, is thought to strongly control the overturning circulation and climate sensitivity (Gnanadesikan 1999; Sijp et al. 2006).
Munk (1966) pointed out the role of vertical mixing in the global overturning circulation using simple theoretical arguments. He discusses the balance between upwelling and vertical mixing in the thermocline. Given 25 Sv (Sv ≡ 106 m3 s−1) of diapycnal upwelling, Munk showed that the globally averaged vertical mixing coefficient D must be approximately 10−4 m2 s−1.
Zika et al. (2009), similar to the analysis of Munk (1966), explore both the along-isopycnal and diapycnal components of the Southern Ocean meridional overturning circulation. The along-isopycnal component of the overturning is “upwelled” along isopycnal layers to the surface boundary layer of the Southern Ocean. Zika et al. found that, to upwell 20–50 Sv of dense water, the along-isopycnal mixing coefficient K must be 150–450 m2 s−1 and the vertical mixing coefficient must be 0.5–1.5 × 10−4 m2 s−1 below the thermocline at latitudes of the Antarctic Circumpolar Current. What fraction of the meridional overturning circulation is upwelled through isopycnals below the thermocline or upwelled along isopycnals in the Southern Ocean is unknown.
Most observations of D have been much smaller than those predicted by Munk (1966). From the distribution of the tracer SF6, released below the thermocline in the eastern North Atlantic during the North Atlantic Tracer Release Experiment (NATRE), Ledwell et al. (1993) estimates D to be O(10−5 m2 s−1). Many studies have used microstructure to estimate vertical mixing (e.g., Gregg 1987). These have tended to reveal weak mixing O(10−5 m2 s−1) in most regions. Other in situ measurement approaches, while confirming weak mixing over much of the World Ocean, have identified stronger mixing, O(10−4 to 10−3 m2 s−1), close to rough topography and in energetic regions such as the Southern Ocean (Naveira-Garabato et al. 2004; Kunze et al. 2006; Sloyan 2005). Since direct estimates of diapycnal mixing are sparse and infrequent, it is not clear whether vertical mixing in these energetic regions is sufficient to induce global diapycnal upwelling on the order of 25 Sv.
Few direct estimates exist for K in the oceans. Ledwell et al. (1998) estimated that the SF6 tracer, released in the eastern North Atlantic at an approximate depth of 300 m, mixed horizontally at a rate on the order of 1000 m2 s−1 over the 30 months of their experiment. Horizontal mixing at the sea surface has been estimated to be equal to or greater than that determined by Ledwell (i.e., 1000–10 000 m2 s−1) (Zhurbas and Oh 2004; Marshall et al. 2006; Sallée et al. 2008).
A depth dependence of the mixing coefficient, particularly that which mixes potential vorticity or interface height [commonly referred to as κ in the Eulerian coordinate parameterization of Gent et al. (1995), with interface height being the depth of isopycnals], is reinforced by adjoint inversions (Ferraira et al. 2005) and eddy resolving models (Eden and Greatbatch 2008). It should be noted that calculations of K, in eddy resolving simulations, are not trivial, even with complete knowledge of a model’s full velocity and density fields (Eden et al. 2007).
Inverse methods are tools used to diagnose the ocean circulation from hydrographic data. Such methods seldom give insight into mixing processes. The most well known method is the box inverse method (Wunsch 1978). In the box method, oceanic sections bound regions of the oceans or entire ocean basins. In each region mass, heat, salt, and other properties are conserved in isopycnal layers. Inverse box models simultaneously solve a set of equations for various unknowns, including reference level velocities and diapycnal fluxes of properties or the vertical mixing coefficient D. Box inversions are always underdetermined, and the final solution depends largely on the choice of model and model variance. This sensitivity is particularly strong when mixing coefficients are considered as unknowns in inverse studies (Tziperman 1988). Inverse models have tended to perform better when individual diapycnal property fluxes are considered (Sloyan and Rintoul 2000) or when observed mixing coefficients are included as “knowns” (St. Laurent et al. 2001). Recent adjoint and gridded methods (Wunsch and Heimbach 2007; Herbei et al. 2008) have shown promise but have, so far, been unable to resolve the spatial structure of D and K.
Zika et al. (2010, hereafter ZMS10) present a new inverse method called the tracer-contour inverse method. This inverse method is less sensitive to error than existing methods and has been validated against the output of a numerical model. With these advances and the presence of new observations that give a better representation of the upper 2000 m of the ocean, it may now be possible to accurately infer the ocean circulation and rates of along-isopycnal and vertical mixing directly from hydrographic data. The purpose of this study is twofold: first, to demonstrate the utility of this new inverse method by applying it to ocean observations in the North Atlantic (Fig. 1) where direct estimates of vertical and along-isopycnal mixing exist and, second, to reveal the depth dependence of the mixing coefficients in that region.
This article is structured as follows: In section 2, we briefly describe the tracer-contour inverse method and the ocean climatology used. In section 3, we present the results of the tracer-contour inverse method, as applied to a region of the North Atlantic. Conclusions are given in section 4. In appendix B, we test the sensitivity of the method to the adjustment of various parameters and differing data sources.
2. The tracer-contour inverse method and hydrographic data
Here we describe the tracer-contour inverse method of ZMS10 and the data used in this study. We then explain how the inverse method is applied in the eastern North Atlantic.
a. The tracer-contour inverse method
The density variable used is neutral density (McDougall 1987), and the temperature variable used is conservative temperature (McDougall 2003). The previous two points are not necessarily critical to the analysis presented here and the variables potential temperature θ and potential density, σ1 or σ2, could be used instead.
In (2), if there is mean advection across a density surface, wγ, there must be a compensating effect of vertical mixing (1/ηγ ≈ γzz/γz and/or Dz) or nonlinear effects such as cabbeling and thermobaricity (1/ηγ). Equation (1) differs from (2) in that the vertical coordinate moves with the isopycnal surfaces, whereas the lateral coordinate is fixed in space; hence, there is no γt /γz term in (1) as there would be if depth z was the vertical coordinate. Equation (2) is similar to that used by Munk (1966), as it relates diapycnal upwelling to vertical mixing but also includes nonlinear processes (i.e., the K term on the right-hand side).
Equation (3) is the tracer-contour equation and relates mixing on an isopycnal γ to the flow on a reference pressure p0. In (3), Ψp0 is the geostrophic streamfunction on p0 (Fig. 2). The Coriolis frequency is f. The coordinate xΘ is that running perpendicular to n on the isopycnal; that is, ∫dxΘ is the integral along a temperature contour on an isopycnal. The tracer-contour specific volume anomaly is δcontour(p) = 1/ρ(S, Θ, p) − 1/ρ(Scontour, Θcontour, p), where Scontour and Θcontour are the conservative temperature and salinity values of the contour between (x2, y2) and (x1, y1); that is, δcontour is redefined for each contour.
Equation (4) is the layer conservation equation for volume and any conservative tracer C and, like (3), relates mixing on γ to flow at p0 (Fig. 2). For conservation of volume C = 1, for conservation of temperature anomaly C = Θ −
The tracer-contour inverse method combines aspects of the three major inverse modeling concepts of modern oceanography, the box, beta spiral, and Bernoulli methods, and the advective–diffusive balance concept is motivated by Munk (1966). By assuming the mixing coefficients are constant on each layer or have some spatial form, the tracer-contour inverse method is overdetermined. There is a direct relationship between lateral advection and diffusion, and the diffusion terms are leading order. In the method, there is a linear relationship between advection and vertical mixing through the Θ − S curvature (d2Θ/dS2|x,y, a quantity less susceptible to noise due to heave than Θzz and Szz are individually). Down-temperature and down-salinity gradient transports along isopycnals are explicitly considered in the tracer-contour inverse method, and it is these transports that constitute the along-isopycnal component of the thermohaline overturning circulation. In the tracer-contour inverse method, a streamfunction variable is solved for, along “sections,” allowing for direct integration with the box inverse method and the computation of section transports for comparison with, or constraint by, shipboard observations.
b. Data and setup of the inversion
The climatology of Durack and Wijffels (2010, hereafter DW10) is used in this study. Hydrographic data from shipboard observations and Argo floats, up to and including 2008, have been used in their analysis. Annual and interannual cycles as well as a linear trend have been locally fit to the data, minimizing temporal aliasing. The data was averaged by DW10 on approximate neutral density surfaces (γrf, Jackett and McDougall 2005), spaced by 0.025 kg m−3 on a 1° latitude × 2° longitude grid. The tracer-contour inverse method has also been applied to the climatology of Gouretski and Koltermann (2004, hereafter GK04), and these results are discussed in appendix B. The inversion is carried out for the NATRE region in the North Atlantic (20°–30°N, 40°–25°W) shown in Fig. 1.
Regularly spaced points on the boundary of the domain (between 20° and 30°N, 40° and 25°W), given by the grid of the climatology, are referred to as “casts,” and the eastern, northern, western, and southern boundaries of the domain are referred to as “sections.” There is a total of 34 casts. Tracer contours (contours of constant Θ and S in this study) are defined on γrf surfaces between 200 and 1800 m. As in ZMS10, contours are defined, on each surface, by the temperature and salinity at the cast locations; that is, there is about one tracer contour per cast per surface. At each point on a contour the flow across the contour, in the direction n, is related to along-isopycnal and diapycnal mixing through the scale lengths λh, 
Maps of the lateral mixing coefficient at the sea surface from Zhurbas and Oh (2004) suggest gradients of the mixing coefficient are very small in the latitude bands considered here. For this reason, we choose to find solutions where K does not vary laterally and hence ignore the ∇γK · n term in (1). It is likely that at the boundary between different mixing regimes this term will become important and the sensitivity of the method to this choice is left to future work. By default we assume a statistically steady state, hence assuming the unsteady term Θt|γ/|∇γΘ| is negligible. This choice, and sensitivity to it, are discussed in appendix B.
Volume, conservative temperature anomaly, and salinity anomaly are conserved on 28 isopycnal layers between 200 and 1800 m (4). Each “layer” in which properties are conserved may encompass more than one of the “surfaces” on which tracer contours are defined. The density range of each layer is chosen such that at least 40 tracer contours are within each layer, allowing K and D to be well resolved within each layer. The diapycnal fluxes of properties due to advection wγ are related to mixing through (2). The geostrophic advection across a section is v · m, where m is the unit vector normal to the section (Fig. 2). This advection is related to the geostrophic streamfunction along the section. The difference in geostrophic streamfunction is related, at each pair of casts, to the reference level streamfunction Ψp0 at 1000 db (C2). A streamfunction is solved for, on the reference pressure, at each of the 34 cast locations. The diffusivities K and D are solved for on each density layer.
Contours with very small along-isopycnal gradients (∇γΘ < 2 × 10−8 K m−1) are excluded. The results are generally insensitive to minor changes of the choice of reference level, depth range, minimum number of contours per layer, and the minimum along-isopycnal temperature gradient. The inversion is more sensitive to the factors considered in appendix B.
3. Results
a. Vertical mixing
The tracer-contour inverse method gives direct estimates of the mixing coefficients K and D (Figs. 3 and 4). The vertical mixing coefficient D in the upper 200–1000 m of the water column is O(10−5 m2 s−1). There is a modest increase in D toward the deepest layers at 1800 m. The deepest value of D is 2.6 ± 0.8 × 10−5 at 1800 m.
A SF6 tracer has been released in the eastern North Atlantic at an approximate depth of 300 m, and its dispersion used to infer a vertical mixing coefficient (Ledwell et al. 1993, 1998). Two observation of the effective mixing coefficient for that tracer were made: 1.2(±0.1) × 10−5 m2 s−1 in the spring of 1992 and 1.7(±0.1) × 10−5 m2 s−1 in the fall of 1993 (Fig. 3). As there are only two measurements, a suitable error range for the “long term mean” D is not possible. Promisingly, the error range of D in the upper 200–400-m depth range, as estimated in this study, almost overlaps with the spring 1992 measurements of Ledwell et al. (1993).
Microstructure measurements of a vertical mixing coefficient in the NATRE region at depths of 200–400 m are generally O(10−5 m2 s−1). For example, St. Laurent and Schmitt (1999) observed mixing of 0.1 − 1.3 × 10−5 ± 10−5 m2 s−1 considering diffusivities for density (Dρ = 0.1 × 10−5 ± 0.3 × 10−5 m2 s−1), potential temperature (Dθ = 0.8 × 10−5 ± 0.1 × 10−5 m2 s−1), and salinity (DS = 1.3 × 10−5 ± 0.1 × 10−5 m2 s−1). Unlike St. Laurent and Schmitt, we do not consider double-diffusive convection since we have taken the vertical mixing coefficient D to be the same for both Θ and S. Ferrari and Polzin (2005) find an equally low turbulent diffusivity close to 0.7 × 10−5 m2 s−1. There are large uncertainties in microstructure estimates of D owing to the assumption of a turbulent mixing efficiency. Despite this uncertainty, the estimates of D presented here agree very strongly with the microstructure estimates of D of Ferrari and Polzin (2005) and St. Laurent and Schmitt (1999). The increase with depth of D toward 1800 m, observed in this study, is also consistent with Ferrari and Polzin (2005).
Microstructure estimates of D are near instantaneous, representing turbulent motions evolving over minutes and hours. The estimates of D by Ledwell et al. (1993, 1998) represent mixing of a tracer integrated over a 6–18-month time period. The estimate of D presented here represents the long-term mean effective mixing of heat and salt across isopycnals. That the three approaches (microstructure, tracer release, and tracer-contour inverse method)—each with their respective time scales, spatial scale, and uncertainties—agree on a value of D O(10−5 m2 s−1) is encouraging. Our results support the hypothesis that, below the thermocline in the eastern North Atlantic, the D estimated by microstructure studies is equivalent to the long-term mean D and is thus an appropriate value for use in ocean circulation models. This canonical “background” value is consistent with the hypothesis that mixing is low away from boundary layers and is unable to sustain the upwelling required in the balance of Munk (1966).
b. Along-isopycnal mixing
The along-isopycnal mixing coefficient K is estimated to be O(1000 m2 s−1) within the upper 200–400 m of the water column, consistent with the North Atlantic Tracer Release Experiment studies for the same region (Ledwell et al. 1998; Sundermeyer and Price 1998). The mixing coefficient smoothly decreases with depth. Between 500- and 1000-m depth, K reduces to around 0–200 m2 s−1 (Fig. 4). These values are broadly consistent with nearby estimates from float trajectories (Joyce et al. 1998; Spall et al. 1993); 3He, 3H, and tritium tracer studies (Jenkins 1987, 1998); and other inverse studies (Armi and Stommel 1983; Zika and McDougall 2008). Ferrari and Polzin (2005) are able to derive a vertical profile of the along-isopycnal mixing coefficient using a mixing length argument and combining CTD and mooring data from NATRE and the Subduction Experiment (Joyce et al. 1998). There is uncertainty in the mixing efficiency that Ferrari and Polzin (2005) choose; thus, they are able to adjust the overall magnitude of K to fit the observations shown in Fig. 4. Our vertical profile, which cannot be adjusted and is a direct output of the inverse method, is very similar to that of Ferrari and Polzin (2005). There are no scaling parameters used in this study.
The vertical profile of K presented here is consistent with the findings of Ferraira et al. (2005), which suggested a depth dependence to the mixing coefficient K that is mostly O(1000 m2 s−1) in the upper 200–500 m of the water column, decreasing to O(100 m2 s−1) at 1000–2000-m depth. However, in the NATRE region, the estimates of Ferraira et al. (2005) are mostly negative. Recent parameterization efforts have focused on the need for a depth dependence of the along-isopycnal mixing coefficient (Eden and Greatbatch 2008; Danabasoglu and Marshall 2008), motivated largely by inferences from eddy resolving models and theoretical arguments. This study represents one of the few demonstrations of a strongly depth-dependent along-isopycnal mixing coefficient, inferred directly from observational data.
c. Geostrophic flow
Because the geostrophic streamfunction at the reference level is an output of the tracer-contour inverse method, the full depth geostrophic velocity and transport are easily inferred (Fig. 5). The mean geostrophic velocity is to the southwest throughout most of the study region. The velocity is intensified in the upper 400 m of the water column with velocities of up to 2–3 cm s−1. The flow is generally stronger in the zonal direction and intensifies in the southwest corner of the domain. This pattern is consistent with the largely westward migration of the SF6 tracer throughout the NATRE study of approximately 1.3 cm s−1 (Ledwell et al. 1998, see their Table 1). The flow is also consistent with the southeastward flowing limb of the wind-driven North Atlantic subtropical gyre in Sverdrup balance (Sverdrup 1947). As deep velocities are estimated to be small, O(0.2 cm s−1), the inferred circulation is consistent with studies that use atlas data and a deep level of no motion (e.g., Reid 1978). The typical columnar nature of the velocity field, seen in conventional inverse calculations using one-time hydrographic surveys, is not seen here partly because individual eddies are not present in the climatology.
The deep flow is in no way minimized or assumed small in this analysis. Despite this, the solution for the geostrophic velocity at 1800 m is small, O(0.2 cm s−1), relative to the upper 500 m O(2 cm s−1). However, the inferred contribution to the section transport between 1800 m and the sea floor across each section is significant, O(2–3 Sv), demonstrating the need for an appropriate estimate of the deep velocity. Some structure to deep currents is revealed, including subsurface zonal currents.
4. Conclusions
The tracer-contour inverse method has been applied to the mean hydrography of the eastern North Atlantic. Vertical profiles of the vertical mixing coefficient D and along-isopycnal mixing coefficient K have been diagnosed, as well as the mean geostrophic circulation. The mixing coefficients and circulation patterns are consistent with in situ measurements.
Vertical mixing is found to be weak, O(10−5 m2 s−1), throughout most of the water column. These low mixing values are consistent with the hypothesis that abyssal mixing is small in the subtropical oceans and is insufficient, on its own, to close the meridional overturning circulation. It is likely that the methods demonstrated in this study will be able to identify regions of intense mixing and quantify the global spatial structure of the mixing coefficients.
In the NATRE region, along-isopycnal mixing is found to be strong, O(1000 m2 s−1), near the thermocline and reduces with depth to background values O(100 m2 s−1) below 500-m depth. Our results confirm that along-isopycnal mixing varies with depth as inferred by Ferrari and Polzin (2005) and Ferraira et al. (2005), further motivating parameterizations for this effect in coarse-resolution ocean models (Eden and Greatbatch 2008; Danabasoglu and Marshall 2008). Uniquely, we reveal a depth dependence of the along-isopycnal mixing coefficient from mean hydrographic data alone.
The eastern North Atlantic is thought to be a region of low eddy kinetic energy and small diapycnal mixing and, as such, is a suitable region for determining the background K and D. The tracer-contour inverse method accurately reproduces in situ observations there, suggesting that the method can be used to identify regions of both weak and intense vertical and along-isopycnal mixing in the world’s oceans.
Acknowledgments
We thank Paul Durack and Dr. Susan Wijffels for the use of their climatology. We are grateful to Andrew Meijers, Frank Colberg, and two anonymous reviewers for their helpful comments. This work contributes to the CSIRO Climate Change Research Program and has been partially supported by the CSIRO Wealth from Oceans Flagship and the Australian Research Council Centre of Excellence for Mathematics and Statistics of Complex Systems.
REFERENCES
Armi, L., and H. Stommel, 1983: Four views of a portion of the North Atlantic subtropical gyre. J. Phys. Oceanogr., 13 , 828–857.
Danabasoglu, G., and J. Marshall, 2008: Effects of vertical variations of thickness diffusivity in an ocean general circulation model. Ocean Modell., 18 , 122–141.
Durack, P. J., and S. E. Wijffels, 2010: Fifty-year trends in global ocean salinities and their relationship to broad-scale warming. J. Climate, 23 , 4342–4362.
Eden, C., and R. J. Greatbatch, 2008: Towards a mesoscale eddy closure. Ocean Modell., 20 , 223–239.
Eden, C., R. J. Greatbatch, and J. Willebrand, 2007: A diagnosis of thickness fluxes in an eddy-resolving model. J. Phys. Oceanogr., 37 , 727–742.
Ferraira, D., J. Marshall, and P. Heimbach, 2005: Estimating eddy stresses by fitting dynamics to observations using a residual-mean ocean circulation model and its adjoint. J. Phys. Oceanogr., 35 , 1891–1910.
Ferrari, R., and K. L. Polzin, 2005: Finescale structure of the T–S relation in the eastern North Pacific. J. Phys. Oceanogr., 35 , 1437–1454.
Gent, P. R., J. Willebrand, T. J. McDougall, and J. C. McWilliams, 1995: Parameterizing eddy-induced tracer transports in ocean circulation models. J. Phys. Oceanogr., 25 , 463–474.
Gnanadesikan, A., 1999: A simple predictive model for the structure of the oceanic pycnocline. Science, 283 , 2077–2079.
Gouretski, V. V., and K. P. Koltermann, 2004: WOCE global hydrographic climatology. Berichte des Bundesamtes für Seeschifffahrt und Hydrographie Tech. Rep. 35, 49 pp.
Gregg, M. C., 1987: Diapycnal mixing in the thermocline: A review. J. Geophys. Res., 92 , 5249–5289.
Herbei, R., I. W. McKegue, and K. G. Speer, 2008: Gyres and jets: Inversion of tracer data for ocean circulation structure. J. Phys. Oceanogr., 38 , 1180–1202.
Jackett, D. R., and T. J. McDougall, 1997: A neutral density variable for the world’s oceans. J. Phys. Oceanogr., 27 , 237–263.
Jackett, D. R., and T. J. McDougall, 2005: The material derivative of neutral density. J. Mar. Res., 63 , 159–185.
Jenkins, W. J., 1987: 3H and 3He in the Beta Triangle: Observations of gyre ventilation and oxygen utilization rates. J. Phys. Oceanogr., 17 , 763–783.
Jenkins, W. J., 1998: Studying subtropical thermocline ventilation and circulation using tritium and He-3. J. Geophys. Res., 103 , 15817–15831.
Joyce, T. M., J. R. Luyten, A. Kubryakov, F. B. Bahr, and J. S. Pallant, 1998: Meso- to large-scale structure of subducting water in the subtropical gyre of the eastern North Atlantic Ocean. J. Phys. Oceanogr., 28 , 40–61.
Kunze, E., E. Firing, J. Hummon, T. Chereskin, and A. M. Thurnherr, 2006: Global abyssal mixing inferred from lowered ADCP shear and CTD strain profiles. J. Phys. Oceanogr., 36 , 1553–1576.
Ledwell, J., A. Watson, and C. Law, 1993: Evidence of slow mixing across the pycnocline from an open-ocean tracer experiment. Nature, 364 , 701–702.
Ledwell, J., A. Watson, and C. Law, 1998: Mixing of a tracer in the pycnocline. J. Geophys. Res., 103 , 21499–21529.
Lux, M., H. Mercier, and M. Arhan, 2001: Interhemispheric exchanges of mass and heat in the Atlantic Ocean in January-March 1993. Deep-Sea Res. I, 48 , 605–638.
Manabe, S., and R. J. Stouffer, 1988: Two stable equilibria of a coupled ocean–atmosphere model. J. Climate, 1 , 841–866.
Marshall, J., E. Shuckburgh, H. Jones, and C. Hill, 2006: Estimates and implications of surface eddy diffusivity in the Southern Ocean from tracer transport. J. Phys. Oceanogr., 36 , 1806–1821.
McDougall, T. J., 1984: The relative roles of diapycnal and isopycnal mixing on subsurface water mass conservation. J. Phys. Oceanogr., 14 , 1577–1589.
McDougall, T. J., 1987: Neutral surfaces. J. Phys. Oceanogr., 17 , 1950–1964.
McDougall, T. J., 2003: Potential enthalpy: A conservative oceanic variable for evaluating heat content and heat fluxes. J. Phys. Oceanogr., 33 , 945–963.
McDougall, T. J., and P. C. McIntosh, 2001: The temporal-residual-mean velocity. Part II: Isopycnal interpretation and the tracer and momentum equations. J. Phys. Oceanogr., 31 , 1222–1246.
McDougall, T. J., and A. Klocker, 2010: An approximate geostrophic streamfunction for use in density surfaces. Ocean Modell., 32 , 107–117. doi:10.1016/j.ocemod.2009.10.006.
Munk, W. H., 1966: Abyssal recipes. Deep-Sea Res., 13 , 707–730.
Munk, W. H., and K. Wunsch, 1998: Abyssal recipes II: Energetics of tidal and wind mixing. Deep-Sea Res., 45 , 1977–2010.
Reid, J. L., 1978: On the middepth circulation and salinity field in the North Atlantic Ocean. J. Geophys. Res., 83 , 5063–5067.
Sallée, J. B., K. Speer, R. Morrow, and R. Lumpkin, 2008: An estimate of Lagrangian eddy statistics and diffusion in the mixed layer of the Southern Ocean. J. Mar. Res., 66 , 441–463.
Sijp, W. P., M. Bates, and M. H. England, 2006: Can isopycnal mixing control the stability of the thermohaline circulation in ocean climate models? J. Climate, 19 , 5637–5651.
Sloyan, B. M., 2005: Spatial variability of mixing in the Southern Ocean. Geophys. Res. Lett., 32 , L18603. doi:10.1029/2005GL023568.
Sloyan, B. M., and S. R. Rintoul, 2000: Estimates of area-averaged diapycnal fluxes from basin-scale budgets. J. Phys. Oceanogr., 30 , 2320–2341.
Spall, M. A., P. L. Richardson, and J. Price, 1993: Advection and eddy-mixing in the Mediterranean Salt Tongue. J. Mar. Res., 51 , 797–818.
St. Laurent, L. C., and R. W. Schmitt, 1999: The contribution of salt fingers to vertical mixing in the North Atlantic Tracer Release Experiment. J. Phys. Oceanogr., 29 , 1404–1424.
St. Laurent, L. C., J. M. Toole, and R. W. Schmitt, 2001: Buoyancy forcing by turbulence above rough topography in the abyssal Brazil Basin. J. Phys. Oceanogr., 31 , 3476–3495.
Sundermeyer, M. A., and J. F. Price, 1998: Lateral mixing and the North Atlantic Tracer Release Experiment: Observations and numerical simulations of Lagrangian particles and a passive tracer. J. Geophys. Res., 103 , 21481–21497.
Sverdrup, H. U., 1947: Wind-driven currents in a baroclinic ocean, with application to the equatorial currents of the eastern Pacific. Proc. Natl. Acad. Sci. USA, 33 , 318–326.
Tziperman, E., 1988: Calculating the time mean oceanic general circulation and mixing coefficients from hydrographic data. J. Phys. Oceanogr., 18 , 519–525.
Wunsch, C., 1978: The North Atlantic general circulation west of 50°W determined by inverse methods. Rev. Geophys. Space Phys., 16 , 583–620.
Wunsch, C., and P. Heimbach, 2007: Practical global oceanic state estimation. Physica D, 230 , 197–208. doi:10.1016/j.physd.2006.09.040.
Zhurbas, V., and I. S. Oh, 2004: Drifter-derived maps of lateral diffusivity in the Pacific and Atlantic Oceans in relation to surface circulation patterns. J. Geophys. Res., 109 , C05015. doi:10.1029/2003JC002241.
Zika, J. D., and T. J. McDougall, 2008: Vertical and lateral mixing processes deduced from the Mediterranean Water signature in the North Atlantic. J. Phys. Oceanogr., 38 , 164–176.
Zika, J. D., B. M. Sloyan, and T. J. McDougall, 2009: Diagnosing the Southern Ocean overturning from tracer fields. J. Phys. Oceanogr., 39 , 2926–2940.
Zika, J. D., T. J. McDougall, and B. M. Sloyan, 2010: A tracer-contour inverse method for estimating ocean circulation and mixing. J. Phys. Oceanogr., 40 , 26–47.
APPENDIX A
Derivation of the Advective–Diffusive Balance Equations with a Linear Equation of State
So, for any tracer, the component of the flow across contours on isopycnals is related to along-isopycnal mixing through the isopycnal curvature of the tracer (∇γ2C/|∇γC|) and to vertical mixing through the curvature of the tracer with respect to density d2C/dγ2|x,y. Note that the curvature of Θ or S with respect to γ is proportional to the curvature of Θ with respect to S.
APPENDIX B
Uncertainties and Sensitivity Analysis
Mixing parameters and diapycnal fluxes determined using inverse methods are known to be highly sensitive to the model details and the data sources used (Tziperman 1988; Lux et al. 2001; St. Laurent et al. 2001). This sensitivity is not evident in standard statistical parameters [i.e., the singular values of the SVD solution] and is best explored through sensitivity tests. It is thus pertinent to test whether the results presented here, for the vertical and along-isopycnal mixing coefficients, are robust in terms of their sensitivity to changes to key assumptions and parameterizations, equation weighting, and data sources. To test the robustness of our results, we conduct a series of sensitivity tests, in each case, assessing the quantitative and qualitative changes made to the vertical and along-isopycnal mixing rates inferred. We show that the choice of equation weighting, parameterizations, and data source has a largely negligible effect, the last being the most influential.
In this study, the bolus velocity v* is represented as a downgradient diffusion of the thickness of isopycnal layers h. The rate of diffusion of thickness is controlled by the mixing coefficient KPV such that v*h = −KPV∇γ(h). In the analysis presented in section 3, we assume that the coefficient for thickness is the same as that for temperature and salinity (i.e., KPV = K). These choices are arbitrary, and one could equally mix potential vorticity f/h instead of h and/or allow K and KPV to vary independently. To test the influence of this parameter in the eastern North Atlantic, we remove the thickness term from the inverse method. The effect on the mixing coefficients and their standard error is negligible (dotted line, Figs. B1a,b). This suggests that the choice of parameterization is inconsequential in this region and the method would be unable to diagnose KPV independently here. However, the bolus velocity is likely to be of leading order importance in the Southern Ocean and regions of steeply sloping isopycnals (Gent et al. 1995; McDougall and McIntosh 2001).
Here a statistically steady balance is assumed such that the trend in Θ and S, on isopycnal surfaces, is negligible when compared to the mean advection, along-isopycnal mixing, and diapycnal mixing terms (i.e., we have assumed Θt|γ = 0). In the climatology of DW10, trends have been quantified and are discussed in DW10. We test the importance of the trend term by retaining the Θt|γ/∇γΘ term in Eq. (1). This term effectively represents the gradual movement of temperature contours on isopycnals. The effect on the inverse calculation in the eastern North Atlantic is negligible (dashed line, Figs. B1a,b), suggesting a steady-state balance is a valid assumption in the region considered. This assumption may break down in regions where strong changes are detected on isopycnals and/or along isopycnal gradients are particularly low.
In ZMS10 and here, the relative weighting of the contour equation (3) versus the large-scale conservation, or box equation (4), is chosen by minimizing the condition number of the matrix 𝗔 in the system 𝗔x = b. It is conventional, in box inverse modeling studies, to weight equations by a prior error estimate. As there is significant uncertainty in the estimation of prior error, the choice is still largely arbitrary. To test the sensitivity of our mixing estimates to the weighting of contour versus box equations, we simply multiply and divide by 2, the weighting given by the minimum condition number criterion. This change in the weights gives a negligible impact on the along-isopycnal mixing coefficient profile in the upper 2000 m and a small impact on the vertical mixing profile in the upper 1000 m (Figs. B1c,d). There is a nonnegligible effect on the vertical mixing coefficient D in the 1000–2000-m depth range, although they show an increase with depth for all choices of weights. There is a small amount of sensitivity to the choice of reference level, partly because this choice affects the magnitude of b, which in turn affects the relative weights of each equation.
The scale lengths, λγ and 
We have mapped the GK04 climatology onto neutral density surfaces (γn, Jackett and McDougall 1997). We interpolate onto the same grid as DW10 so that a clear comparison can be made between the two. The velocity and mixing coefficients generated by the tracer-contour inverse method are generally more noisy when using the GK04 climatology. When using GK04, the results are more sensitive to the proximity of the analysis to the mixed layer. For this reason, we avoid using data shallower than 250 m. The vertical mixing coefficient, inferred using GK04, is weak, O(10−5 m2 s−1), throughout all of the water column, increasing somewhat in the 1000–2000-m depth range (Fig. B1e). The along-isopycnal mixing coefficient K, inferred using GK04, decreases from highs O(1000 m2 s−1) in the upper 500 m of the water column to O(100 m2 s−1) below 1000 m (Fig. B1f). These results are very similar to those derived using DW10 and earlier estimates summarized in Figs. 3 and 4.
APPENDIX C
Equivalent Form of the Thermal Wind Equation
Pressure on the γrf = 26.525 kg m−3 (σ0 ≈ 26.75 kg m−3) surface from the climatology of DW10. Also marked is the region considered in this study (solid black rectangle), the site of mooring C from the Subduction Experiment (black cross, Joyce et al. 1998), and the NATRE study regions (Ledwell et al. 1998) showing the approximate SF6 tracer extent during fall 1992 (dark gray rectangle), spring 1993 (medium gray rectangle), and November 1994 (light gray rectangle).
Citation: Journal of Physical Oceanography 40, 8; 10.1175/2010JPO4360.1
Schematic showing how the tracer-contour inverse method is implemented. Mixing at points along a tracer contour between (x1, y1) and (x2, y2) on the isopycnal surface γ is related to a geostrophic streamfunction on the reference pressure p0 using (3), shown in green. The unit vector n is normal to the contour on the isopycnal. The tracer C, typically volume, conservative temperature anomaly, or salinity anomaly, is conserved in an isopycnal layer bounded by sections using (4). Part of an isopycnal layer is shown in light blue, bound by surfaces with tracer properties Cu and Cl. In practice, the layer is a complete rectangle bound laterally by all casts (from i = 1 to i = L). Properties are advected or mixed through the upper and lower bounding surfaces of the layer and across each section. The velocity across a section is related to the streamfunction on the reference surface using a form of the thermal wind equation (appendix C). The unit vector m is normal to the section. In this study, Ψp0 is solved for each cast i and the mixing coefficients K and D are solved for on each layer.
Citation: Journal of Physical Oceanography 40, 8; 10.1175/2010JPO4360.1
Vertical mixing coefficient D, as determined using the tracer-contour inverse method (solid line) and the standard error of the inversion (light gray shading). Also shown are the estimates of D from Ledwell et al. (1993, 1998) and Ferrari and Polzin (2005) and their reported uncertainties.
Citation: Journal of Physical Oceanography 40, 8; 10.1175/2010JPO4360.1
Along-isopycnal mixing coefficient K, as determined using the tracer-contour inverse method (solid line) and the standard error of the inversion (light gray shading). Also shown are the estimates of K from Joyce et al. (1998), Ledwell et al. (1998), Jenkins (1998), Armi and Stommel (1983), Spall et al. (1993), and Zika and McDougall (2008) and their reported uncertainties.
Citation: Journal of Physical Oceanography 40, 8; 10.1175/2010JPO4360.1
(top) Cross-sectional velocity v · m (cm s−1) on γrf = 26.525 kg m−3 (σ0 ≈ 26.75 kg m−3). Length of arrows represents magnitude of v · m, whereas the direction indicates whether the flow is into or out of the study region. Each section has a corresponding velocity figure starting from (left) western section moving clockwise to (right) the southern section. The gray box indicates the study region. Ledwell et al. (1998) release a SF6 tracer at σ0 = 26.75 kg m−3, which is advected to the southeast, consistent with the flow field observed here. (bottom) Cross section of velocity on neutral density surfaces v · m (positive values are directed out of the study region). As in (top), each section has a corresponding velocity figure.
Citation: Journal of Physical Oceanography 40, 8; 10.1175/2010JPO4360.1
(left) Vertical mixing coefficient D and (right) along-isopycnal mixing coefficient K, as determined using the tracer-contour inverse method. Black lines represent the estimated value, whereas the area between the gray lines represent the standard error. (a),(b) For the standard settings (solid line), no bolus velocity (dotted line) and including an unsteady term (dashed line); (c),(d) for the standard settings (solid line), increasing the weight on the contour equations (dashed line) and increasing the weight on the box equations (dotted line); and (e),(f) as in (c),(d), but using the climatology of GK04.
Citation: Journal of Physical Oceanography 40, 8; 10.1175/2010JPO4360.1
