Climatological Annual Cycle of the Salinity Budgets of the Subtropical Maxima

Benjamin K. Johnson Department of Atmospheric and Oceanic Science, University of Maryland, College Park, College Park, Maryland

Search for other papers by Benjamin K. Johnson in
Current site
Google Scholar
PubMed
Close
,
Frank O. Bryan National Center for Atmospheric Research, Boulder, Colorado

Search for other papers by Frank O. Bryan in
Current site
Google Scholar
PubMed
Close
,
Semyon A. Grodsky Department of Atmospheric and Oceanic Science, University of Maryland, College Park, College Park, Maryland

Search for other papers by Semyon A. Grodsky in
Current site
Google Scholar
PubMed
Close
, and
James A. Carton Department of Atmospheric and Oceanic Science, University of Maryland, College Park, College Park, Maryland

Search for other papers by James A. Carton in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

Six subtropical salinity maxima (Smax) exist: two each in the Pacific, Atlantic, and Indian Ocean basins. The north Indian (NI) Smax lies in the Arabian Sea while the remaining five lie in the open ocean. The annual cycle of evaporation minus precipitation (EP) flux over the Smax is asymmetric about the equator. Over the Northern Hemisphere Smax, the semiannual harmonic is dominant (peaking in local summer and winter), while over the Southern Hemisphere Smax, the annual harmonic is dominant (peaking in local winter). Regardless, the surface layer salinity for all six Smax reaches a maximum in local fall and minimum in local spring. This study uses a multidecade integration of an eddy-resolving ocean circulation model to compute salinity budgets for each of the six Smax. The NI Smax budget is dominated by eddy advection related to the evolution of the seasonal monsoon. The five open-ocean Smax budgets reveal a common annual cycle of vertical diffusive fluxes that peak in winter. These Smax have regions on their eastward and poleward edges in which the vertical salinity gradient is destabilizing. These destabilizing gradients, in conjunction with wintertime surface cooling, generate a gradually deepening wintertime mixed layer. The vertical salinity gradient sharpens at the base of the mixed layer, making the water column susceptible to salt finger convection and enhancing vertical diffusive salinity fluxes out of the Smax into the ocean interior. This process is also observed in Argo float profiles and is related to the formation regions of subtropical mode waters.

Denotes Open Access content.

The National Center for Atmospheric Research is sponsored by the National Science Foundation.

Corresponding author address: Benjamin Johnson, Department of Atmospheric and Oceanic Science, University of Maryland, College Park, 3424 Computer and Space Science Bldg., College Park, MD 20742. E-mail: bjohnson@atmos.umd.edu

Abstract

Six subtropical salinity maxima (Smax) exist: two each in the Pacific, Atlantic, and Indian Ocean basins. The north Indian (NI) Smax lies in the Arabian Sea while the remaining five lie in the open ocean. The annual cycle of evaporation minus precipitation (EP) flux over the Smax is asymmetric about the equator. Over the Northern Hemisphere Smax, the semiannual harmonic is dominant (peaking in local summer and winter), while over the Southern Hemisphere Smax, the annual harmonic is dominant (peaking in local winter). Regardless, the surface layer salinity for all six Smax reaches a maximum in local fall and minimum in local spring. This study uses a multidecade integration of an eddy-resolving ocean circulation model to compute salinity budgets for each of the six Smax. The NI Smax budget is dominated by eddy advection related to the evolution of the seasonal monsoon. The five open-ocean Smax budgets reveal a common annual cycle of vertical diffusive fluxes that peak in winter. These Smax have regions on their eastward and poleward edges in which the vertical salinity gradient is destabilizing. These destabilizing gradients, in conjunction with wintertime surface cooling, generate a gradually deepening wintertime mixed layer. The vertical salinity gradient sharpens at the base of the mixed layer, making the water column susceptible to salt finger convection and enhancing vertical diffusive salinity fluxes out of the Smax into the ocean interior. This process is also observed in Argo float profiles and is related to the formation regions of subtropical mode waters.

Denotes Open Access content.

The National Center for Atmospheric Research is sponsored by the National Science Foundation.

Corresponding author address: Benjamin Johnson, Department of Atmospheric and Oceanic Science, University of Maryland, College Park, 3424 Computer and Space Science Bldg., College Park, MD 20742. E-mail: bjohnson@atmos.umd.edu

1. Introduction

The broad patterns of surface layer salinity (SLS) are set by the interaction of surface waters with the atmosphere (Wüst 1935). The descending branch of the Hadley cell drives anticyclonic flow over the ocean basins, causing evaporation to exceed precipitation, elevating SLS and creating six subtropical salinity maxima (hereinafter Smax). Five of the Smax lie primarily in the open ocean: the Pacific and Atlantic Smax are located around 25°N and 20°S; the south Indian Smax is around 30°S (Gordon et al. 2015). The north Indian (NI) Smax lies in the Arabian Sea and Gulf of Oman around 18°N. It is bounded by coasts to the north and west. For brevity, we will refer to each Smax by the first initial of its hemisphere and ocean basin, for example, NA Smax denotes the North Atlantic subtropical salinity maxima. The Smax are bounded by different isohaline values given in Table 1. We selected these values as those that define the maximum volume that is confined to the subtropical region of each basins’ corresponding hemisphere.

Table 1.

The isohaline contour bounding each Smax.

Table 1.

The annual cycle of evaporation minus precipitation (EP) flux over the subtropical regions is asymmetric about the equator (Large and Yeager 2009). Over the Northern Hemisphere (NH) Smax, the semiannual harmonic dominates (peaking in local summer and winter), while over the Southern Hemisphere (SH) Smax, the annual harmonic dominates (peaking in local winter). In spite of this asymmetry in EP forcing, the SLS of all six Smax reaches a maximum in local fall and minimum in local spring [open-ocean Smax (Gordon et al. 2015), NI Smax (Rao and Sivakumar 2003)]. On the eastward and poleward edges of each open-ocean Smax, the vertical gradient of salinity is destabilizing and the net surface heat flux (SHF) changes sign during the transition from summer to winter. Bingham et al. (2012) depict seasonal SLS variability in these areas (see their Fig. 4), which suggests that the position of the Smax within each basin influences SLS variability in the surrounding regions.

This study employs an eddy-resolving ocean general circulation model (OGCM) to compute salinity budgets for all six Smax to address two objectives: determine the importance of individual ocean processes in balancing climatological mean EP forcing and to discover why, in spite of their different annual cycles of EP, the annual cycle of SLS for all six Smax is dominated by the annual harmonic. Following Bryan and Bachman (2015), who studied NA Smax, these budgets are computed within a control volume that is bounded by an isohaline surface in order to eliminate the mean advection term from the budget. NA Smax is perhaps the most well studied Smax, as it was the site of the Salinity Processes in the Upper Ocean Regional Study (SPURS; Lindstrom et al. 2015) and is a region of high-coverage by the Ship of Opportunity Program (SOOP; Reverdin et al. 2007). Previous observation-based studies have emphasized the importance of different processes in the subtropical salinity balance: horizontal advection by Ekman and geostrophic flow (Yu 2011), vertical fluxes [Qu et al. (2011), who only examined NA Smax; Bingham et al. (2012)], and eddy fluxes [Gordon and Giulivi (2014), who only examined NA Smax].

While OGCMs have the ability to simulate regularly gridded fields of oceanographic variables, including diffusion terms that are difficult to measure in situ, they experience a particular challenge that is relevant to the salinity budget. Of the two components that determine the density and stability of seawater, potential temperature θ and salinity S, the latter has been more difficult for OGCMs to model with fidelity (Durack et al. 2013). Because modeled S distributions can be biased—even in this precise model configuration, as discussed by Bryan and Bachman (2015)—this study reinforces its key results with comparisons to Argo float data.

2. Background

The subsurface structure of each Smax creates gradients that constrain large-scale patterns of vertical mixing. The open-ocean Smax each have a subsurface salinity maximum known as subtropical underwater (STUW; O’Connor et al. 2002). The STUW are located westward and equatorward of the sea surface salinity maxima (shown schematically in Fig. 1). This displacement creates contrasting regions of along the edges of the open-ocean Smax. On their eastward and poleward edges, is destabilizing—S decreases with depth, creating a top-heavy salinity distribution. Conversely, on their westward and equatorward edges, is stabilizing—S increases with depth, creating a bottom-heavy salinity distribution. Unlike the open-ocean Smax, the subsurface salinity maximum beneath NI Smax is located approximately beneath the surface maximum, and thus is weak.

Fig. 1.
Fig. 1.

A schematic of an open-ocean Smax. Solid lines represent isohaline contours with increasing contour thickness representing increased salinity. The dark (light) shaded region depicts that in the mixed layer eastward and poleward (westward and equatorward) of the sea surface salinity (SSS) max, ∂zS is stabilizing (destabilizing) because of the displacement of STUW with respect to the SSS maximum.

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

In addition to these regions of contrasting , the stratification of subtropical oceans is affected by an annual cycle of SHF that changes sign. All of the Smax experience net heating during the summer and net cooling during the winter (Fig. 2); however, the annual cycle of SHF over NI Smax is distinct. It has the shortest period of surface cooling (typically from mid-November to early February) and its period of surface heating is interrupted by a relatively weak period during local summer. Consequently, its mixed layer exhibits maximum depth in January and has a secondary maximum in July (Weller et al. 2002). The open-ocean Smax have mixed layers that gradually deepen throughout local winter and rapidly shoal in local spring. For the open-ocean Smax, the coldest wintertime SSTs exist on their eastward and poleward edges.

Fig. 2.
Fig. 2.

The climatological annual cycle of surface heat flux (W m−2) over each Smax in CESM. The x axis is offset by six months for the SH Smax (top x-axis labels) with respect to the NH Smax (bottom x-axis labels) to account for the hemispheric phase shift of the seasons. The color legend is: NA Smax (light purple), SA Smax (light green), NI Smax (dark green), SI Smax (light blue), NP Smax (dark blue), and SP Smax (dark purple).

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

Since and the vertical gradient of potential temperature can compensate for each other (i.e., a column with destabilizing can be stabilized by ), a measure of their composite influence on stability is the density ratio
e1
where and are the thermal expansion and haline contraction coefficients, respectively. The density ratio in the surface layer changes throughout the year as EP fluxes alter and surface cooling in the autumn and winter weakens (warming in the spring and summer strengthens) .

A notable manifestation of this air–sea interaction is the generation of various types of subtropical mode waters (MWs). In a review of subtropical MWs, Hanawa and Talley (2001) develop a classification system organizing them into types, two of which are relevant to this study: type I, the archetype of which is Eighteen Degree MW (Worthington 1959), which forms on the equatorward side of western boundary current extensions because of convection induced by intense winter cooling and is relatively uniform in potential density , and type II, the archetype of which is Madiera MW (Käse et al. 1985; Siedler et al. 1987), which forms on the eastward and poleward edges of the Smax and has variable and density-compensating and (Toyama and Suga 2010).

Type II MWs are of particular interest in regard to the subtropical cells (STCs; McCreary and Lu 1994), which are shallow overturning circulations that connect the subtropical and tropical oceans. STCs provide a mechanism by which the properties of subtropical oceans can propagate subsurface through the Central Waters of the thermocline (e.g., Schmitt 1981; Sprintall and Tomczak 1993) to the tropical oceans, where they can potentially influence the frequency or intensity of El Niño–Southern Oscillation events (e.g., Gu and Philander 1997; Cane et al. 1997). The formation region of type II MWs coincides with the source region of the descending branch of the STCs (e.g., Joyce et al. 1998; Johnson 2006), and these MWs have been identified within the vicinity of NA Smax (Käse et al. 1985), NP Smax (Hautala and Roemmich 1998), and SP Smax (Tsuchiya and Talley 1996).

Type II MW volumes vary significantly throughout the year. After decaying to minimum thickness in the summer, they build throughout the fall and reach peak thickness in the late winter as surface cooling weakens , which, in conjunction with the destabilizing component, decreases the stability of the column and enhances mixing. This wintertime buildup deepens the mixed layer and increases S at depth. The sharpens below the MW layer near the base of the mixed layer, which, in conjunction with the weakened , decreases . This enhances S fluxes from the mixed layer into the thermocline (Yeager and Large 2004) and makes water columns susceptible to salt finger convection (Johnson 2006). Salt finger convection is a type of double-diffusive convection arising in water columns in which θ and S both decrease with depth and is sufficiently low (typically when 1 < < 2). Since the diffusivity of θ is roughly two orders of magnitude greater than the diffusivity of S, water parcels lose heat much faster than they lose salt and become negatively buoyant (Stern 1960), creating finger-like convective plumes. Salt finger convection acts to reduce the potential energy stored in the top-heavy S distribution by transporting salt downward in the convective plumes at a faster rate than they transport heat (which diffuses out of the walls of the convective plumes). This allows salt finger convection to transport density upgradient in a fluid; increase its stability; and, in the context of this discussion, hasten the decay of the MW. Salt finger convection has been inferred from in situ data beneath the type II MW formation regions of the NA (St. Laurent and Schmitt 1999; Kolodziejczyk et al. 2015), SP (Wong and Johnson 2003; Kolodziejczyk and Gaillard 2013), and NP (Sugimoto and Hanawa 2007; Toyama and Suga 2010).

3. Numerical experiment and results

a. Model description

We used the ocean (POP2; Smith et al. 2010) and ice (CICE4; Bailey et al. 2009) components of the Community Earth System Model (CESM) in an eddy-resolving oceanic tripole grid (Murray 1996) with zonal resolution of 0.1°, meridional resolution of 0.1° cos(lat), and 62 vertical levels (10-m spacing above 160 m and gradually increasing spacing toward the sea floor).

The equation of state is computed using the method of McDougall et al. (2003). Freshwater forcing is imposed as a virtual salt flux. Vertical mixing is represented using the K-profile parameterization (KPP) of Large et al. (1994). Turbulent shear mixing and double diffusive convection are represented by increasing vertical diffusivity above the background diffusivity κ. Shear mixing is a function of the gradient Richardson number. Salt finger convection is a function of , as defined in (1), according to the empirical parameterization of St. Laurent and Schmitt (1999):
e2
where and is 1 × 10−2 m2 s−1 and 7 × 10−3 m2 s−1 for salt and heat, respectively. The mixed layer depth is computed as the shallowest local maximum of the buoyancy gradient.

b. Model experiment

We forced the ocean and ice components of CESM in accordance with the Coordinated Ocean–Ice Reference Experiments (CORE; Griffies et al. 2009) guidelines. The COREv2 dataset is composed of a blend of NCEP reanalysis (Kalnay et al. 1996) and a variety of satellite-based reconstructions of surface atmospheric state and flux fields (Large and Yeager 2009). COREv2 includes interannually varying forcing data, spanning the years 1948–2009 and a single “normal year,” which is an exactly repeating synthetic year derived from the interannually varying data. The interannual dataset incorporates observational and reanalysis data from a variety of sources with irregular start and stop times. In 1979, COREv2 begins incorporating data from the Microwave Sounding Unit (MSU; Spencer 1993), the Climate Prediction Center Merged Analysis of Precipitation (CMAP) product (Xie and Arkin 1996), and Global Precipitation Climatology Project (GPCP) data (Huffman et al. 1997). Coastal runoff is imposed using the Dai and Trenberth (2002) climatology. The model was initialized with temperature and salinity distributions from the World Ocean Circulation Experiment (WOCE) climatology (Gouretski and Koltermann 2004) and spun up for 15 model years with COREv2 normal-year forcing. This spinup was shared by the experiment described in Bryan and Bachman (2015), but the integrations branched after spinup: Bryan and Bachman continued by extending their integration for another 5 years under normal-year forcing, and the integration in this study continued for another 33 years under interannually varying forcing corresponding to the years 1977–2009. Since the transition from normal-year forcing to interannual forcing makes the forcing fields discontinuous, we started the interannual forcing in 1977 to mitigate the effects of these perturbations and restricted our analysis to the output from 1979 to 2009. Output was saved as 5-day averages and converted to monthly averages using a weighted mean.

c. Model analysis

To derive our salinity budget, we start with the model’s salinity conservation equation:
e3
where S is the salinity, κ is the vertical diffusivity, κΓ is the KPP counter gradient flux, and is the horizontal hyper-diffusivity. The surface boundary conditions for (3) are
e4
e5
where E is evaporation, P is precipitation, and is the model’s constant reference salinity used to convert freshwater fluxes to a virtual salt flux. We compute the integral of (3) upon climatological monthly mean output (denoted by the overbar) within a volume V that is bounded in the ocean interior by an isohaline surface ∂V, upon which , and at the surface by a time-dependent outcrop area A, within which :
e6
Since the first term of the RHS of (6) is a covariance, we can partition it into time-mean and eddy terms:
e7
The term describing time-mean velocity advection vanishes, resulting in
e8
We denote the terms of (8) from left to right as tendency, eddy advection, vertical diffusion, EP, and horizontal diffusion. Since CESM imposes no-flux boundary conditions at solid surfaces, the volume described here can intersect coastlines and the seafloor. We use (8) to compute the budget on the climatological monthly mean output of the years 1979–2009 of the integration. Bryan and Bachman (2015) demonstrate that since eddy mixing is a two-stage process where eddies first stir the fluid, enhancing local θ/S gradients, and small-scale mixing then acts upon these gradients, computing the budget on climatological monthly mean data allows for the inference of the total mixing induced by the eddies. Finally, to evaluate model bias, we compare the outcrop area for each Smax to the Monthly Isopycnal and Mixed-Layer Ocean Climatology (MIMOC; Schmidtko et al. 2013). Since MIMOC is largely based on Argo data, we use only the model years that overlap with the Argo era (2000–09) during this comparison.

d. Model results and discussion

1) Climatological description of the Smax

Consistent with observations of STUW, subsurface isohaline contours are displaced westward and equatorward of the surface outcrops (Fig. 3) of each open-ocean Smax (resembling Fig. 1). The subsurface isohaline contour is not displaced for NI Smax. The STUW of SA Smax is displaced more equatorward and less westward in comparison with the other Smax because it abuts the coast of Brazil. In the vicinity of the open-ocean Smax, the mean mixed layer depths (MLDs) range from around 50 m in the summer months to around 100–130 m in the winter months. The simulation’s annual SLS cycles, as indicated by outcrop extent (Fig. 4), are roughly in phase with MIMOC observations although the amplitude of the MIMOC annual cycles is larger for all Smax. The phase of these annual SLS cycles is also consistent with in situ observations [NI Smax (Rao and Sivakumar 2003), open-ocean Smax (Gordon et al. 2015)] with seasonal maxima (minima) in local fall (spring). The model tends to have a fresh bias (Fig. 5), with smaller outcrops for each Smax, with the exception of NP Smax. The open-ocean Smax outcrops tend to have a spatial bias, with the MIMOC outcrops extending farther eastward and poleward. This bias suggests the model may overestimate vertical diffusive fluxes in regions that have destabilizing .

Fig. 3.
Fig. 3.

The model climatological mean SLS. Black contours are the mean position of isohaline surface outcrops for each Smax. Gray contours are isohalines at 100 m depth. Each of the open-ocean Smax contours at 100 m depth is positioned westward and equatorward of the surface outcrop. The NI Smax surface and 100-m contours are at roughly the same position. White lines depict transects along which depth profiles are shown in Figs. 8 and 9.

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

Fig. 4.
Fig. 4.

The annual cycle of outcrop extent anomaly (m2) from annual mean value for CESM (black) and MIMOC (light blue).

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

Fig. 5.
Fig. 5.

Model climatological mean SLS (psu, shaded) and comparison of climatological mean Smax outcrops of the integration (black contours) and MIMOC (light blue). Black contours are CESM; white contours are MIMOC. In all cases, with the exception of NP Smax, the MIMOC outcrops are larger, with the largest discrepancy in Smax, where MIMOC outcrop area is nearly double that of the model. The white areas along the eastward and poleward edges of each MIMOC outcrop cover the paths of the Argo floats over the 2-yr periods described in section 4.

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

2) Climatological Smax salinity budgets

We consider the percentage of the climatological annual mean EP that is balanced by each budget term (Table 2). Eddy advection balances the largest percent of EP in NI Smax (−77%), SI Smax (−70%), NP Smax (−60%), and NA Smax (−55%). Vertical diffusion balances the largest percent of EP in SP Smax (−60%) and SA Smax (−55%). The tendency term is small but of the same sign as EP in all six Smax, indicating that the mean salinity of each Smax increased during the integration. While this is consistent with the amplification of SLS variance (Durack and Wijffels 2010) in response to a strengthening hydrological cycle (Held and Soden 2006), this result could also be a consequence of model drift (Durack et al. 2013).

Table 2.

Climatological annual mean magnitude of each budget term and the percent of EP balanced by each of the remaining budget terms (Sv psu). Since the mean advection term should vanish from (7), it provides an estimate of the error in (8). We include its magnitude as an error estimate in the rightmost column.

Table 2.

Turning our attention to the annual cycles of each budget (Fig. 6), the SH Smax exhibit EP minima in the late spring (November for SA and SI; December for SP) and maxima in the winter (June for SP; July for SA and SI). The NH Smax exhibit an annual cycle of EP dominated by the semiannual harmonic with wintertime maxima in January and summertime maxima in June (NP) or July (NA and NI). In spite of this difference in the annual cycle of EP, the annual cycle of net S flux out of the Smax into the ocean interior (which is the sum of the horizontal diffusion, vertical diffusion, and eddy advection terms) is roughly in phase for all of the Smax. In each of the open-ocean Smax, the magnitude of net flux is greatest in the winter while the net flux for NI Smax peaks in the late fall. The magnitude of the vertical diffusion term peaks in the winter for all of the open-ocean Smax. The magnitude of the eddy advection term nearly doubles in the winter for SI Smax and is slightly enhanced for NA and SA Smax. It does not have a noticeable seasonality in NP and SP Smax. In NI Smax, the magnitude of eddy advection reaches a minimum in April and maximum in October. This is likely related to the onset of the southwest monsoon, which develops during the summer and dissipates during the fall and is associated with freshening in the Arabian Sea (Rao and Sivakumar 2003).

Fig. 6.
Fig. 6.

The climatological annual cycle of the salinity budget (Sv psu) for each Smax. The net flux term is the sum of horizontal diffusion, eddy advection, and vertical diffusion.

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

3) Vertical diffusion in the subtropics

Since the annual cycle of S flux from NI Smax into the ocean interior is governed by a distinct process, eddy advection, we continue by examining why the magnitude of the vertical diffusion term peaks in the late winter for the open-ocean Smax. However, since the vertical diffusion term peaks in winter in the open-ocean Smax and is approximately constant throughout the year in the NI Smax, we analyze NI Smax in parallel with the open-ocean Smax because this contrast demonstrates how the annual cycle of vertical diffusion is influenced by the annual cycle of SHF.

Figure 7 of Bingham et al. (2012) shows that in the regions eastward and poleward of the open-ocean Smax, a vertical flux term , where is the seasonal variation of the mixed layer depth, dominates their SLS budget in winter and spring. While this process is captured in the entrainment term, Bingham et al. describe it as a “seasonal heaving of the mixed layer.” We begin by interpolating the vertical diffusion term in the model output to the depth of the mixed layer (Fig. 7), which demonstrates that on the eastward and poleward regions of the Smax, there is indeed a vertical flux of S through the base of the mixed layer into the deeper ocean. Since this process coincides with wintertime surface cooling and a decrease of , we also compute and use (2) to determine the regions where salt finger convection enhances diffusivity at the base of the mixed layer.

Fig. 7.
Fig. 7.

The simulation’s climatological mean vertical diffusive flux of salt though the base of the mixed layer (psu m s−1, shaded). Black and gray contours are isohaline contours as in Fig. 3. Light blue contours are regions where salt finger convection penetrates the base of the mixed layer. In the vicinity of the Smax, these regions of salt finger convection are on the eastward and poleward edges of Smax and correspond with regions of negative vertical diffusive salt fluxes through the mixed layer base (magenta-shaded regions).

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

In the NA, NP, and SP, these areas of S flux from the mixed layer into the ocean interior coincide with the regions of type II MW formation. MWs are typically characterized using the planetary component of potential vorticity,
e9
where f is the Coriolis parameter. The variable Q is used in the study of MW since it quantifies vertical stability, is approximately conserved in geostrophic flow, and is generally much larger in magnitude than the relative potential vorticity component (Hanawa and Talley 2001). To account for the cross-equatorial change in the sign of f, we use |Q| in computations but refer to it as Q within the text. Type II MWs have a “core” (depth of minimum Q) found around 75–100 m during the late winter. While studies use different values of Q to demarcate the boundaries of mode waters, we use the value of Wong and Johnson (2003), where |Q| < 3 × 10−10 m−1 s−1, to examine the extent of low Q water in all basins.

Since the contrasting regions of stability are on diagonally opposite sides of the open-ocean Smax, we create depth transects through each Smax (Figs. 8, 9) along the paths in Fig. 3. Low Q water forms in the surface layer along each transect however, the waters beneath these low Q regions generate salt finger convection only on the eastward and poleward sides of the open-ocean Smax. When low Q water forms on the westward and equatorward side of the Smax, as visible in NP and SP, the waters underlying these low Q volumes do not generate salt finger convection because is stabilizing above the STUW. Similarly, the waters within NI Smax are not susceptible to salt finger convection because is weak. We also depict the seasonal evolution of the mixed layer depth and the salt finger convection diffusivity at intervals during the winter and spring. The NH (SH) open-ocean Smax show mixed layer deepening from January through March (July–September) and rapid shoaling from March through May (September–November). Salt finger convection forms above the mixed layer base in January (July), penetrates the mixed layer base in March (September), and dissipates below the shoaled mixed layer in May (November). This process is depicted in a simplified schematic in Fig. 10.

Fig. 8.
Fig. 8.

Transect plots through each NH Smax along the paths shown in Fig. 3. Climatological January planetary potential vorticity (m−1 s−1, shaded), isohaline bounding each Smax (white contours), and volumes where |Q| < 3.0 × 10−10 m−1 s−1 (magenta). Climatological mixed layer depth for January, March, and May is shown by mustard yellow contours ranging from dark–thick to light–thin. Climatological areas of salt finger convection for January, March, and May are shown by blue contours ranging from dark–thick to light–thin.

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

Fig. 9.
Fig. 9.

As in Fig. 8, but for the SH with July, September, and November replacing January, March, and May.

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

Fig. 10.
Fig. 10.

A simplified schematic of the seasonal deepening of the mixed layer. The solid lines depict isohaline contours. Dashed–dotted lines surrounding the white volumes depict low Q layers. Dotted lines surrounding dark gray volumes depict areas of salt finger convection.

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

Though we have not computed a buoyancy budget in this study, it is important to note that the locations of these regions where the waters become susceptible to salt finger convection are significant because of the connection between type II mode waters and the descending branch of the STCs. This seasonal propagation of salt finger convection provides a diapycnal mixing pathway that transforms near-surface water into deeper, subsurface Central Water. Furthermore, since it mixes density upgradient, the newly transformed Central Waters can escape re-entrainment when the mixed layer deepens in subsequent winters.

4. Argo float selection and analysis

The mixed layer deepening mechanism described in section 3 is significant because it demonstrates how the annual cycle of the salinity budget, type II MW formation, diapycnal mixing induced by salt finger convection, and the descending branch of the STC may be connected. However, these results come with the caveat that the model has a low SLS bias in this region. We reinforce our results using Argo (Roemmich et al. 2009) profiles to demonstrate that this process can be observed in situ. This analysis follows Johnson (2006) for SP and Toyama and Suga (2010) for NP and shows that the process occurs on the eastward and poleward edges of all five open-ocean Smax.

a. Float selection and methodology

We select Argo floats along the eastward and poleward edges of each MIMOC Smax outcrop. To provide assurance that the observed processes are not merely single-year anomalies (since this study attempts to explain processes in the climatological annual cycle), we select Argo floats that remain in the region long enough to provide approximately 2 years of data (73 depth profiles). As a consequence of this requirement, the Argo data are reported with varying start and stop dates. Table 3 displays the ID numbers and duration of data analyzed from each float. The white areas in Fig. 5 depict the path each float traversed during this duration. We follow the analysis of section 3, using (2) to parameterize the salt finger convection diffusivity and (9) to detect the presence of low Q water masses using the same value as in section 3. There are many methods for determining the depth of the mixed layer from in situ data (de Boyer Montégut et al. 2004). We employ the criterion of Monterey and Levitus (1997), who use a threshold method specifically for the study of subtropical MWs: the base of the mixed layer is the pressure at which a profile’s increases from its linearly extrapolated surface value by a set interval, kg m−3.

Table 3.

The float IDs and data duration for the Argo data included in section 4. Each float provides 73 profiles (approximately 2 years) of data within the eastward and poleward region of each Smax.

Table 3.

b. Argo results

The Argo data (Fig. 11) from the open-ocean Smax depict a buildup of a low Q layer within mixed layer as it gradually deepens in response to surface cooling throughout the winter months. Salt finger convection diffusivity peaks beneath the low Q layer and above the base of the mixed layer as the mixed layer restratifies in the late spring. Salt finger convection diffusivity persists year-round but weakens throughout the summer months. The open-ocean Smax have relatively deep maximum mixed layer depths: NP and SP ~ 160 dbar, SA and NA ~ 225 dbar, and SI ~ 240 dbar. In the NI Smax profiles, the maximum mixed layer depth is much shallower (reaching ~90 dbar); the salt finger convection diffusivity is both weak and inconstant.

Fig. 11.
Fig. 11.

Seventy-three depth profiles from each of the Argo floats with data durations listed in Table 3. Depths where |Q| < 3.0 × 10−10 m−1 s−1 are depicted in magenta. Salt finger convection diffusivity as estimated by (2) (blue shading). Depth of the mixed layer depicted with a mustard yellow thick line.

Citation: Journal of Physical Oceanography 46, 10; 10.1175/JPO-D-15-0202.1

5. Summary and conclusions

This study aims to determine the relative importance of internal ocean processes for balancing climatological annual mean EP fluxes and to explain why, in spite of their different annual cycles of EP, the SLS of all six Smax reaches a maximum in local fall and minimum in local spring. Eddy advection plays a primary role in the climatological balance of EP fluxes, from a high of balancing 77% of climatological annual mean EP for NI Smax to a low of 41% for SP Smax. Vertical diffusion also plays a significant role in balancing EP fluxes, from a high of 60% for SP Smax to a low of 28% for NI Smax. The climatological annual mean tendency term is of the same sign as EP fluxes for all six Smax, indicating that they all increased in mean S during the integration. This is consistent with the observed pattern amplification of SLS (Durack and Wijffels 2010) but may be a consequence of model drift (Durack et al. 2013). The horizontal diffusion term is negligible in the budgets of all six Smax.

While the SLS of all six Smax reaches a maximum (minimum) in local fall (spring), the salinity budget of NI Smax is distinct: its annual cycle of net flux into the ocean interior is dominated by eddy advection. The annual cycle of SHF for NI Smax is governed by the seasonal evolution of the monsoon circulation, and thus it experiences a short period of surface cooling (typically from mid-November to early February) and relatively shallow mixed layers. The annual cycle of net flux in the open-ocean Smax budgets is governed by a combination of eddy advection and vertical diffusion for SI Smax and vertical diffusion for the remaining Smax.

The westward and equatorward displacement of STUW with respect to the surface outcrop of the open-ocean Smax creates contrasting regions of stability in their vicinity. On the equatorward and westward regions is stabilizing, while on the eastward and poleward regions is destabilizing. These regions of destabilizing experience longer periods of wintertime surface cooling when compared to NI Smax (typically 4–6 months) and develop thicker layers of near-surface low Q water. This process gradually deepens the mixed layer and sharpens below the “core” of the low Q layer and above the base of the mixed layer. This makes waters susceptible to salt finger convection and enhances vertical diffusive fluxes of S from the Smax. The mixing induced by salt finger convection provides a diapycnal mixing pathway for near-surface low Q waters to transform into subsurface Central Waters. These regions have been identified as type II MW formation regions in NA (Joyce et al. 1998), NP (Ladd and Thompson 2000), and SP (Wong and Johnson 2003) and coincide with the source water regions of the descending branch of the STCs. Though our model results come with the important caveat that CESM has a low SLS bias in these regions compared to observations, this seasonal process is also observed in Argo floats.

This finding, that low Q water is produced in response to wintertime surface cooling in regions of destabilizing adjacent to the open-ocean Smax, has broader implications for our understanding of MW formation. A comparison of our Fig. 7 to Plate 5.4.3 of Hanawa and Talley (2001) shows that tends to be stabilizing or neutral in Hanawa and Talley’s type I MW regions, and the columns are not susceptible to salt finger convection except in the case of SA subtropical MW (SA STMW). Unlike the other open-ocean Smax, SA Smax is located at the westward extent of the basin, and thus large areas within the central SA have destabilizing and are susceptible to salt finger convection. This suggests that the type II MW formation process may generate low Q water in the central SA even though the classification implies that intense cooling associated with the Brazil Current is the source of low Q for SA STMW.

Two other irregularities support this notion:

  1. Type I mode waters are typically uniform in . SA STMW is an outlier in this case because its is variable (Hanawa and Talley 2001), a property generally associated with type II MWs.

  2. Type I mode waters are typically associated with western boundary currents that experience intense wintertime SHF. Kallberg et al. (2005) estimate that the Kuroshio (adjacent to NP STMW) has a maximum wintertime SHF of −401 W m−2 and the Gulf Stream Current (adjacent to NA STMW) has a maximum wintertime SHF of −396 W m−2. In contrast, the Brazil Current (adjacent to SA Smax) has a maximum wintertime SHF of only −170 W m−2.

The original study of MW in SA (Provost et al. 1999) was based on a sparse dataset that predates the Argo era. This suggests that further study may reveal a more detailed distribution of MWs in SA, analogous to recent findings in the NP (Oka and Qiu 2012). The existence of type II MW in the central SA would provide further evidence of the connection between the STCs and type II MWs. The source water region for the descending branch of the SA STC (Fig. 1 of Snowden and Molinari 2003) corresponds to the type II MW formation region in Fig. 7.

Our results come with two important caveats: the S fluxes due to salt finger convection are not differentiated from other contributions to vertical mixing in the diagnosed budgets, and the model has a low S bias in the regions eastward and poleward of the Smax. It should be possible to use the KPP mixing scheme to explicitly quantify the magnitude of S fluxes due to salt finger convection by exporting contributions to κ from each physical process in addition to the total vertical diffusion term. Additionally, St. Laurent and Schmitt (1999) developed (2) by considering mixing schemes in parameter space and found that in low Richardson number conditions, turbulent shear mixing best fit observations, while in high Richardson number and low conditions, salt finger mixing best fit observations. Our configuration of CESM (version 1.0.4) does not screen for low Richardson numbers before imposing (2). It would be prudent to run a sensitivity study in which low Richardson numbers are screened to see if this helps to mitigate the SLS bias in the vicinity of the Smax.

Acknowledgments

Two awards supported this research: S.A.G. was supported by NNX12AF68G and F.O.B. was supported by NNX10AC16G. The model integration described in this paper was carried out on Yellowstone (Computational and Information Systems Laboratory 2012; ark:/85065/d7wd3xhc) provided by NCAR’s Computational and Information Systems Laboratory, sponsored by the National Science Foundation. Furthermore, we thank David Bailey for his assistance in keeping the model integration progressing. We found the comments from the two anonymous reviewers and the guidance of the editor to be invaluable.

REFERENCES

  • Bailey, D., M. Holland, E. Hunke, B. Lipscomb, B. Briegleb, C. Bitz, and J. Schramm, 2009: Community Ice CodE (CICE) user’s guide version 4.0. NCAR Doc., 22 pp. [Available online at http://www.cesm.ucar.edu/models/cesm1.0/cice/ice_usrdoc.pdf.]

  • Bingham, F. M., G. R. Foltz, and M. J. McPhaden, 2012: Characteristics of the seasonal cycle of surface layer salinity in the global ocean. Ocean Sci., 8, 915929, doi:10.5194/os-8-915-2012.

    • Search Google Scholar
    • Export Citation
  • Bryan, F. O., and S. Bachman, 2015: Isohaline salinity budget of the North Atlantic salinity maximum. J. Phys. Oceanogr., 45, 724736, doi:10.1175/JPO-D-14-0172.1.

    • Search Google Scholar
    • Export Citation
  • Cane, M. A., and Coauthors, 1997: Twentieth-century sea surface temperature trends. Science, 275, 957960, doi:10.1126/science.275.5302.957.

    • Search Google Scholar
    • Export Citation
  • Computational and Information Systems Laboratory, 2012: Yellowstone: IBM iDataPlex System (Climate Simulation Laboratory). National Center for Atmospheric Research, accessed 17 August 2016. [Available online at http://n2t.net/ark:/85065/d7wd3xhc.]

  • Dai, A., and K. E. Trenberth, 2002: Estimates of freshwater discharge from continents: Latitudinal and seasonal variations. J. Hydrometeor., 3, 660687, doi:10.1175/1525-7541(2002)003<0660:EOFDFC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • de Boyer Montégut, C., G. Madec, A. S. Fischer, A. Lazar, and D. Iudicone, 2004: Mixed layer depth over the global ocean: An examination of profile data and a profile-based climatology. J. Geophys. Res., 109, C12003, doi:10.1029/2004JC002378.

    • Search Google Scholar
    • Export Citation
  • 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, 43424362, doi:10.1175/2010JCLI3377.1.

    • Search Google Scholar
    • Export Citation
  • Durack, P. J., S. E. Wijffels, and T. P. Boyer, 2013: Long-term salinity challenges and implications for the global water cycle. Ocean Circulation and Climate, G. Siedler et al., Eds., Academic Press, 727–758.

  • Gordon, A. L., and C. F. Giulivi, 2014: Ocean eddy freshwater flux convergence into the North Atlantic subtropics. J. Geophys. Res. Oceans, 119, 33273335, doi:10.1002/2013JC009596.

    • Search Google Scholar
    • Export Citation
  • Gordon, A. L., C. F. Giulivi, J. Busecke, and F. M. Bingham, 2015: Differences among subtropical surface salinity patterns. Oceanography, 28, 3239, doi:10.5670/oceanog.2015.02.

    • Search Google Scholar
    • Export Citation
  • Gouretski, V., and K. Koltermann, 2004: WOCE Global Hydrographic Climatology. BSH Tech. Rep. 35, Bundesamtes für Seeschifffahrt und Hydrographie, 52 pp. [Available online at www.bsh.de/de/Produkte/Buecher/Berichte_/Bericht35/Bericht1.pdf.]

  • Griffies, S. M., and Coauthors, 2009: Coordinated Ocean-ice Reference Experiments (COREs). Ocean Modell., 26, 146, doi:10.1016/j.ocemod.2008.08.007.

    • Search Google Scholar
    • Export Citation
  • Gu, D. F., and S. G. H. Philander, 1997: Interdecadal climate fluctuations that depend on exchanges between the tropics and extratropics. Science, 275, 805807, doi:10.1126/science.275.5301.805.

    • Search Google Scholar
    • Export Citation
  • Hanawa, K., and L. D. Talley, 2001: Mode waters. Ocean Circulation and Climate, G. Siedler et al., Eds., Academic Press, 373–386.

  • Hautala, S. L., and D. H. Roemmich, 1998: Subtropical mode water in the Northeast Pacific Basin. J. Geophys. Res., 103, 13 05513 066, doi:10.1029/98JC01015.

    • Search Google Scholar
    • Export Citation
  • Held, I. M., and B. J. Soden, 2006: Robust responses of the hydrological cycle to global warming. J. Climate, 19, 56865699, doi:10.1175/JCLI3990.1.

    • Search Google Scholar
    • Export Citation
  • Huffman, G. J., and Coauthors, 1997: The Global Precipitation Climatology Project (GPCP) Combined Precipitation Dataset. Bull. Amer. Meteor. Soc., 78, 520, doi:10.1175/1520-0477(1997)078<0005:TGPCPG>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Johnson, G. C., 2006: Generation and initial evolution of a mode water θS anomaly. J. Phys. Oceanogr., 36, 739751, doi:10.1175/JPO2895.1.

    • Search Google Scholar
    • Export Citation
  • 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, 4061, doi:10.1175/1520-0485(1998)028<0040:MTLSSO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Kallberg, P., P. Berrisford, B. Hoskins, A. Simmons, S. Uppala, S. Lamy-Thépaut, and R. Hine, 2005: ERA-40 Atlas. ECMWF Tech. Rep. 19, 191 pp.

  • Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437471, doi:10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Käse, R. H., W. Zenk, T. B. Sanford, and W. Hiller, 1985: Currents, fronts and eddy fluxes in the Canary Basin. Prog. Oceanogr., 14, 231257, doi:10.1016/0079-6611(85)90013-8.

    • Search Google Scholar
    • Export Citation
  • Kolodziejczyk, N., and F. Gaillard, 2013: Variability of the heat and salt budget in the subtropical southeastern Pacific mixed layer between 2004 and 2010: Spice injection mechanism. J. Phys. Oceanogr., 43, 18801898, doi:10.1175/JPO-D-13-04.1.

    • Search Google Scholar
    • Export Citation
  • Kolodziejczyk, N., G. Reverdin, and A. Lazar, 2015: Interannual variability of the mixed layer winter convection and spice injection in the eastern subtropical North Atlantic. J. Phys. Oceanogr., 45, 504525, doi:10.1175/JPO-D-14-0042.1.

    • Search Google Scholar
    • Export Citation
  • Ladd, C., and L. Thompson, 2000: Formation mechanisms for North Pacific central and eastern subtropical mode waters. J. Phys. Oceanogr., 30, 868887, doi:10.1175/1520-0485(2000)030<0868:FMFNPC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Large, W. G., and S. G. Yeager, 2009: The global climatology of an interannually varying air-sea flux data set. Climate Dyn., 33, 341364, doi:10.1007/s00382-008-0441-3.

    • Search Google Scholar
    • Export Citation
  • Large, W. G., J. C. McWilliams, and S. C. Doney, 1994: Oceanic vertical mixing: A review and a model with a nonlocal boundary layer parametrization. Rev. Geophys., 32, 363403, doi:10.1029/94RG01872.

    • Search Google Scholar
    • Export Citation
  • Lindstrom, E., F. Bryan, and R. Schmitt, 2015: SPURS: Salinity Processes in the Upper-ocean Regional Study—The North Atlantic Experiment Introduction. Oceanography, 28, 1419, doi:10.5670/oceanog.2015.01.

    • Search Google Scholar
    • Export Citation
  • McCreary, J. P., and P. Lu, 1994: Interaction between the subtropical and equatorial ocean circulations: The subtropical cell. J. Phys. Oceanogr., 24, 466497, doi:10.1175/1520-0485(1994)024<0466:IBTSAE>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • McDougall, T. J., D. R. Jackett, D. G. Wright, and R. Feistel, 2003: Accurate and computationally efficient algorithms for potential temperature and density of seawater. J. Atmos. Oceanic Technol., 20, 730741, doi:10.1175/1520-0426(2003)20<730:AACEAF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Monterey, G., and S. Levitus, 1997: Seasonal Variability of Mixed Layer Depth for the World Ocean. NOAA Atlas NESDIS 14, 102 pp. [Available online at ftp://ftp.nodc.noaa.gov/pub/data.nodc/woa/PUBLICATIONS/Atlas14.pdf.]

  • Murray, R. J., 1996: Explicit generation of orthogonal grids for ocean models. J. Comput. Phys., 126, 251273, doi:10.1006/jcph.1996.0136.

    • Search Google Scholar
    • Export Citation
  • O’Connor, B. M., R. A. Fine, K. A. Maillet, and D. B. Olson, 2002: Formation rates of subtropical underwater in the Pacific Ocean. Deep-Sea Res. I, 49, 15711590, doi:10.1016/S0967-0637(02)00087-0.

    • Search Google Scholar
    • Export Citation
  • Oka, E., and B. Qiu, 2012: Progress of North Pacific mode water research in the past decade. J. Oceanogr., 68, 520, doi:10.1007/s10872-011-0032-5.

    • Search Google Scholar
    • Export Citation
  • Provost, C., C. Escoffier, K. Maamaatuaiahutapu, A. Kartavtseff, and V. Garcon, 1999: Subtropical mode waters in the South Atlantic Ocean. J. Geophys. Res., 104, 21 03321 049, doi:10.1029/1999JC900049.

    • Search Google Scholar
    • Export Citation
  • Qu, T. D., S. Gao, and I. Fukumori, 2011: What governs the North Atlantic salinity maximum in a global GCM? Geophys. Res. Lett., 38, L07602, doi:10.1029/2011GL046757.

    • Search Google Scholar
    • Export Citation
  • Rao, R. R., and R. Sivakumar, 2003: Seasonal variability of sea surface salinity and salt budget of the mixed layer of the north Indian Ocean. J. Geophys. Res., 108, 3009, doi:10.1029/2001JC000907.

    • Search Google Scholar
    • Export Citation
  • Reverdin, G., E. Kestenare, C. Frankignoul, and T. Delcroix, 2007: Surface salinity in the Atlantic Ocean (30°S–50°N). Prog. Oceanogr., 73, 311340, doi:10.1016/j.pocean.2006.11.004.

    • Search Google Scholar
    • Export Citation
  • Roemmich, D., and Coauthors, 2009: The Argo Program: Observing the global ocean with profiling floats. Oceanography, 22, 3443, doi:10.5670/oceanog.2009.36.

    • Search Google Scholar
    • Export Citation
  • Schmidtko, S., G. C. Johnson, and J. M. Lyman, 2013: MIMOC: A global monthly isopycnal upper-ocean climatology with mixed layers. J. Geophys. Res. Oceans, 118, 16581672, doi:10.1002/jgrc.20122.

    • Search Google Scholar
    • Export Citation
  • Schmitt, R. W., 1981: Form of the temperature-salinity relationship in the central water: Evidence for double-diffusive mixing. J. Phys. Oceanogr., 11, 10151026, doi:10.1175/1520-0485(1981)011<1015:FOTTSR>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Siedler, G., A. Kuhl, and W. Zenk, 1987: The Madiera Mode Water. J. Phys. Oceanogr., 17, 15611570, doi:10.1175/1520-0485(1987)017<1561:TMMW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Smith, R., and Coauthors, 2010: The Parallel Ocean Program (POP) reference manual. Rep. LAUR-10-01853, 141 pp. [Available online at http://www.cesm.ucar.edu/models/cesm1.0/pop2/doc/sci/POPRefManual.pdf.]

  • Snowden, D. P., and R. L. Molinari, 2003: Subtropical cells in the Atlantic ocean: An observational summary. Interhemispheric Water Exchange in the Atlantic Ocean, G. J. Goni and P. Malanotte-Rizzoli, Eds., Elsevier Oceanography Series, Vol. 68, Elsevier, 287312, doi:10.1016/S0422-9894(03)80151-4.

  • Spencer, R. W., 1993: Global oceanic precipitation from the MSU during 1979–91 and comparisons to other climatologies. J. Climate, 6, 13011326, doi:10.1175/1520-0442(1993)006<1301:GOPFTM>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Sprintall, J., and M. Tomczak, 1993: On the formation of central water and thermocline ventilation in the Southern Hemisphere. Deep-Sea Res. I, 40, 827848, doi:10.1016/0967-0637(93)90074-D.

    • Search Google Scholar
    • Export Citation
  • St. Laurent, L., and R. W. Schmitt, 1999: The contribution of salt fingers to vertical mixing in the North Atlantic Tracer Release Experiment. J. Phys. Oceanogr., 29, 14041424, doi:10.1175/1520-0485(1999)029<1404:TCOSFT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Stern, M. E., 1960: The salt-fountain and thermohaline convection. Tellus, 12A, 172175, doi:10.1111/j.2153-3490.1960.tb01295.x.

  • Sugimoto, S., and K. Hanawa, 2007: Further evidence for non-reemergence of winter SST anomalies in the North Pacific eastern subtropical mode water area. J. Oceanogr., 63, 625635, doi:10.1007/s10872-007-0055-0.

    • Search Google Scholar
    • Export Citation
  • Toyama, K., and T. Suga, 2010: Vertical structure of North Pacific mode waters. Deep-Sea Res. II, 57, 11521160, doi:10.1016/j.dsr2.2009.12.004.

    • Search Google Scholar
    • Export Citation
  • Tsuchiya, M., and L. D. Talley, 1996: Water-property distributions along an eastern Pacific hydrographic section at 135W. J. Mar. Res., 54, 541564, doi:10.1357/0022240963213583.

    • Search Google Scholar
    • Export Citation
  • Weller, R. A., A. S. Fischer, D. L. Rudnick, C. C. Eriksen, T. D. Dickey, J. Marra, C. Fox, and R. Leben, 2002: Moored observations of upper-ocean response to the monsoons in the Arabian Sea during 1994–1995. Deep-Sea Res. II, 49, 21952230, doi:10.1016/S0967-0645(02)00035-8.

    • Search Google Scholar
    • Export Citation
  • Wong, A. P. S., and G. C. Johnson, 2003: South Pacific Eastern Subtropical Mode Water. J. Phys. Oceanogr., 33, 14931509, doi:10.1175/1520-0485(2003)033<1493:SPESMW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Worthington, L. V., 1959: The 18° water in the Sargasso Sea. Deep-Sea Res., 5, 297305, doi:10.1016/0146-6313(58)90026-1.

  • Wüst, G., 1935: The Stratosphere of the Atlantic Ocean: Scientific Results of the German Atlantic Expedition of the Research Vessel “Meteor” 1925–1927. Amervid, 112 pp.

  • Xie, P. P., and P. A. Arkin, 1996: Analyses of global monthly precipitation using gauge observations, satellite estimates, and numerical model predictions. J. Climate, 9, 840858, doi:10.1175/1520-0442(1996)009<0840:AOGMPU>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Yeager, S. G., and W. G. Large, 2004: Late-winter generation of spiciness on subducted isopycnals. J. Phys. Oceanogr., 34, 15281547, doi:10.1175/1520-0485(2004)034<1528:LGOSOS>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Yu, L. S., 2011: A global relationship between the ocean water cycle and near-surface salinity. J. Geophys. Res., 116, C10025, doi:10.1029/2010JC006937.

    • Search Google Scholar
    • Export Citation
Save
  • Bailey, D., M. Holland, E. Hunke, B. Lipscomb, B. Briegleb, C. Bitz, and J. Schramm, 2009: Community Ice CodE (CICE) user’s guide version 4.0. NCAR Doc., 22 pp. [Available online at http://www.cesm.ucar.edu/models/cesm1.0/cice/ice_usrdoc.pdf.]

  • Bingham, F. M., G. R. Foltz, and M. J. McPhaden, 2012: Characteristics of the seasonal cycle of surface layer salinity in the global ocean. Ocean Sci., 8, 915929, doi:10.5194/os-8-915-2012.

    • Search Google Scholar
    • Export Citation
  • Bryan, F. O., and S. Bachman, 2015: Isohaline salinity budget of the North Atlantic salinity maximum. J. Phys. Oceanogr., 45, 724736, doi:10.1175/JPO-D-14-0172.1.

    • Search Google Scholar
    • Export Citation
  • Cane, M. A., and Coauthors, 1997: Twentieth-century sea surface temperature trends. Science, 275, 957960, doi:10.1126/science.275.5302.957.

    • Search Google Scholar
    • Export Citation
  • Computational and Information Systems Laboratory, 2012: Yellowstone: IBM iDataPlex System (Climate Simulation Laboratory). National Center for Atmospheric Research, accessed 17 August 2016. [Available online at http://n2t.net/ark:/85065/d7wd3xhc.]

  • Dai, A., and K. E. Trenberth, 2002: Estimates of freshwater discharge from continents: Latitudinal and seasonal variations. J. Hydrometeor., 3, 660687, doi:10.1175/1525-7541(2002)003<0660:EOFDFC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • de Boyer Montégut, C., G. Madec, A. S. Fischer, A. Lazar, and D. Iudicone, 2004: Mixed layer depth over the global ocean: An examination of profile data and a profile-based climatology. J. Geophys. Res., 109, C12003, doi:10.1029/2004JC002378.

    • Search Google Scholar
    • Export Citation
  • 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, 43424362, doi:10.1175/2010JCLI3377.1.

    • Search Google Scholar
    • Export Citation
  • Durack, P. J., S. E. Wijffels, and T. P. Boyer, 2013: Long-term salinity challenges and implications for the global water cycle. Ocean Circulation and Climate, G. Siedler et al., Eds., Academic Press, 727–758.

  • Gordon, A. L., and C. F. Giulivi, 2014: Ocean eddy freshwater flux convergence into the North Atlantic subtropics. J. Geophys. Res. Oceans, 119, 33273335, doi:10.1002/2013JC009596.

    • Search Google Scholar
    • Export Citation
  • Gordon, A. L., C. F. Giulivi, J. Busecke, and F. M. Bingham, 2015: Differences among subtropical surface salinity patterns. Oceanography, 28, 3239, doi:10.5670/oceanog.2015.02.

    • Search Google Scholar
    • Export Citation
  • Gouretski, V., and K. Koltermann, 2004: WOCE Global Hydrographic Climatology. BSH Tech. Rep. 35, Bundesamtes für Seeschifffahrt und Hydrographie, 52 pp. [Available online at www.bsh.de/de/Produkte/Buecher/Berichte_/Bericht35/Bericht1.pdf.]

  • Griffies, S. M., and Coauthors, 2009: Coordinated Ocean-ice Reference Experiments (COREs). Ocean Modell., 26, 146, doi:10.1016/j.ocemod.2008.08.007.

    • Search Google Scholar
    • Export Citation
  • Gu, D. F., and S. G. H. Philander, 1997: Interdecadal climate fluctuations that depend on exchanges between the tropics and extratropics. Science, 275, 805807, doi:10.1126/science.275.5301.805.

    • Search Google Scholar
    • Export Citation
  • Hanawa, K., and L. D. Talley, 2001: Mode waters. Ocean Circulation and Climate, G. Siedler et al., Eds., Academic Press, 373–386.

  • Hautala, S. L., and D. H. Roemmich, 1998: Subtropical mode water in the Northeast Pacific Basin. J. Geophys. Res., 103, 13 05513 066, doi:10.1029/98JC01015.

    • Search Google Scholar
    • Export Citation
  • Held, I. M., and B. J. Soden, 2006: Robust responses of the hydrological cycle to global warming. J. Climate, 19, 56865699, doi:10.1175/JCLI3990.1.

    • Search Google Scholar
    • Export Citation
  • Huffman, G. J., and Coauthors, 1997: The Global Precipitation Climatology Project (GPCP) Combined Precipitation Dataset. Bull. Amer. Meteor. Soc., 78, 520, doi:10.1175/1520-0477(1997)078<0005:TGPCPG>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Johnson, G. C., 2006: Generation and initial evolution of a mode water θS anomaly. J. Phys. Oceanogr., 36, 739751, doi:10.1175/JPO2895.1.

    • Search Google Scholar
    • Export Citation
  • 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, 4061, doi:10.1175/1520-0485(1998)028<0040:MTLSSO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Kallberg, P., P. Berrisford, B. Hoskins, A. Simmons, S. Uppala, S. Lamy-Thépaut, and R. Hine, 2005: ERA-40 Atlas. ECMWF Tech. Rep. 19, 191 pp.

  • Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437471, doi:10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Käse, R. H., W. Zenk, T. B. Sanford, and W. Hiller, 1985: Currents, fronts and eddy fluxes in the Canary Basin. Prog. Oceanogr., 14, 231257, doi:10.1016/0079-6611(85)90013-8.

    • Search Google Scholar
    • Export Citation
  • Kolodziejczyk, N., and F. Gaillard, 2013: Variability of the heat and salt budget in the subtropical southeastern Pacific mixed layer between 2004 and 2010: Spice injection mechanism. J. Phys. Oceanogr., 43, 18801898, doi:10.1175/JPO-D-13-04.1.

    • Search Google Scholar
    • Export Citation
  • Kolodziejczyk, N., G. Reverdin, and A. Lazar, 2015: Interannual variability of the mixed layer winter convection and spice injection in the eastern subtropical North Atlantic. J. Phys. Oceanogr., 45, 504525, doi:10.1175/JPO-D-14-0042.1.

    • Search Google Scholar
    • Export Citation
  • Ladd, C., and L. Thompson, 2000: Formation mechanisms for North Pacific central and eastern subtropical mode waters. J. Phys. Oceanogr., 30, 868887, doi:10.1175/1520-0485(2000)030<0868:FMFNPC>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Large, W. G., and S. G. Yeager, 2009: The global climatology of an interannually varying air-sea flux data set. Climate Dyn., 33, 341364, doi:10.1007/s00382-008-0441-3.

    • Search Google Scholar
    • Export Citation
  • Large, W. G., J. C. McWilliams, and S. C. Doney, 1994: Oceanic vertical mixing: A review and a model with a nonlocal boundary layer parametrization. Rev. Geophys., 32, 363403, doi:10.1029/94RG01872.

    • Search Google Scholar
    • Export Citation
  • Lindstrom, E., F. Bryan, and R. Schmitt, 2015: SPURS: Salinity Processes in the Upper-ocean Regional Study—The North Atlantic Experiment Introduction. Oceanography, 28, 1419, doi:10.5670/oceanog.2015.01.

    • Search Google Scholar
    • Export Citation
  • McCreary, J. P., and P. Lu, 1994: Interaction between the subtropical and equatorial ocean circulations: The subtropical cell. J. Phys. Oceanogr., 24, 466497, doi:10.1175/1520-0485(1994)024<0466:IBTSAE>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • McDougall, T. J., D. R. Jackett, D. G. Wright, and R. Feistel, 2003: Accurate and computationally efficient algorithms for potential temperature and density of seawater. J. Atmos. Oceanic Technol., 20, 730741, doi:10.1175/1520-0426(2003)20<730:AACEAF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Monterey, G., and S. Levitus, 1997: Seasonal Variability of Mixed Layer Depth for the World Ocean. NOAA Atlas NESDIS 14, 102 pp. [Available online at ftp://ftp.nodc.noaa.gov/pub/data.nodc/woa/PUBLICATIONS/Atlas14.pdf.]

  • Murray, R. J., 1996: Explicit generation of orthogonal grids for ocean models. J. Comput. Phys., 126, 251273, doi:10.1006/jcph.1996.0136.

    • Search Google Scholar
    • Export Citation
  • O’Connor, B. M., R. A. Fine, K. A. Maillet, and D. B. Olson, 2002: Formation rates of subtropical underwater in the Pacific Ocean. Deep-Sea Res. I, 49, 15711590, doi:10.1016/S0967-0637(02)00087-0.

    • Search Google Scholar
    • Export Citation
  • Oka, E., and B. Qiu, 2012: Progress of North Pacific mode water research in the past decade. J. Oceanogr., 68, 520, doi:10.1007/s10872-011-0032-5.

    • Search Google Scholar
    • Export Citation
  • Provost, C., C. Escoffier, K. Maamaatuaiahutapu, A. Kartavtseff, and V. Garcon, 1999: Subtropical mode waters in the South Atlantic Ocean. J. Geophys. Res., 104, 21 03321 049, doi:10.1029/1999JC900049.

    • Search Google Scholar
    • Export Citation
  • Qu, T. D., S. Gao, and I. Fukumori, 2011: What governs the North Atlantic salinity maximum in a global GCM? Geophys. Res. Lett., 38, L07602, doi:10.1029/2011GL046757.

    • Search Google Scholar
    • Export Citation
  • Rao, R. R., and R. Sivakumar, 2003: Seasonal variability of sea surface salinity and salt budget of the mixed layer of the north Indian Ocean. J. Geophys. Res., 108, 3009, doi:10.1029/2001JC000907.

    • Search Google Scholar
    • Export Citation
  • Reverdin, G., E. Kestenare, C. Frankignoul, and T. Delcroix, 2007: Surface salinity in the Atlantic Ocean (30°S–50°N). Prog. Oceanogr., 73, 311340, doi:10.1016/j.pocean.2006.11.004.

    • Search Google Scholar
    • Export Citation
  • Roemmich, D., and Coauthors, 2009: The Argo Program: Observing the global ocean with profiling floats. Oceanography, 22, 3443, doi:10.5670/oceanog.2009.36.

    • Search Google Scholar
    • Export Citation
  • Schmidtko, S., G. C. Johnson, and J. M. Lyman, 2013: MIMOC: A global monthly isopycnal upper-ocean climatology with mixed layers. J. Geophys. Res. Oceans, 118, 16581672, doi:10.1002/jgrc.20122.

    • Search Google Scholar
    • Export Citation
  • Schmitt, R. W., 1981: Form of the temperature-salinity relationship in the central water: Evidence for double-diffusive mixing. J. Phys. Oceanogr., 11, 10151026, doi:10.1175/1520-0485(1981)011<1015:FOTTSR>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Siedler, G., A. Kuhl, and W. Zenk, 1987: The Madiera Mode Water. J. Phys. Oceanogr., 17, 15611570, doi:10.1175/1520-0485(1987)017<1561:TMMW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Smith, R., and Coauthors, 2010: The Parallel Ocean Program (POP) reference manual. Rep. LAUR-10-01853, 141 pp. [Available online at http://www.cesm.ucar.edu/models/cesm1.0/pop2/doc/sci/POPRefManual.pdf.]

  • Snowden, D. P., and R. L. Molinari, 2003: Subtropical cells in the Atlantic ocean: An observational summary. Interhemispheric Water Exchange in the Atlantic Ocean, G. J. Goni and P. Malanotte-Rizzoli, Eds., Elsevier Oceanography Series, Vol. 68, Elsevier, 287312, doi:10.1016/S0422-9894(03)80151-4.

  • Spencer, R. W., 1993: Global oceanic precipitation from the MSU during 1979–91 and comparisons to other climatologies. J. Climate, 6, 13011326, doi:10.1175/1520-0442(1993)006<1301:GOPFTM>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Sprintall, J., and M. Tomczak, 1993: On the formation of central water and thermocline ventilation in the Southern Hemisphere. Deep-Sea Res. I, 40, 827848, doi:10.1016/0967-0637(93)90074-D.

    • Search Google Scholar
    • Export Citation
  • St. Laurent, L., and R. W. Schmitt, 1999: The contribution of salt fingers to vertical mixing in the North Atlantic Tracer Release Experiment. J. Phys. Oceanogr., 29, 14041424, doi:10.1175/1520-0485(1999)029<1404:TCOSFT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Stern, M. E., 1960: The salt-fountain and thermohaline convection. Tellus, 12A, 172175, doi:10.1111/j.2153-3490.1960.tb01295.x.

  • Sugimoto, S., and K. Hanawa, 2007: Further evidence for non-reemergence of winter SST anomalies in the North Pacific eastern subtropical mode water area. J. Oceanogr., 63, 625635, doi:10.1007/s10872-007-0055-0.

    • Search Google Scholar
    • Export Citation
  • Toyama, K., and T. Suga, 2010: Vertical structure of North Pacific mode waters. Deep-Sea Res. II, 57, 11521160, doi:10.1016/j.dsr2.2009.12.004.

    • Search Google Scholar
    • Export Citation
  • Tsuchiya, M., and L. D. Talley, 1996: Water-property distributions along an eastern Pacific hydrographic section at 135W. J. Mar. Res., 54, 541564, doi:10.1357/0022240963213583.

    • Search Google Scholar
    • Export Citation
  • Weller, R. A., A. S. Fischer, D. L. Rudnick, C. C. Eriksen, T. D. Dickey, J. Marra, C. Fox, and R. Leben, 2002: Moored observations of upper-ocean response to the monsoons in the Arabian Sea during 1994–1995. Deep-Sea Res. II, 49, 21952230, doi:10.1016/S0967-0645(02)00035-8.

    • Search Google Scholar
    • Export Citation
  • Wong, A. P. S., and G. C. Johnson, 2003: South Pacific Eastern Subtropical Mode Water. J. Phys. Oceanogr., 33, 14931509, doi:10.1175/1520-0485(2003)033<1493:SPESMW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Worthington, L. V., 1959: The 18° water in the Sargasso Sea. Deep-Sea Res., 5, 297305, doi:10.1016/0146-6313(58)90026-1.

  • Wüst, G., 1935: The Stratosphere of the Atlantic Ocean: Scientific Results of the German Atlantic Expedition of the Research Vessel “Meteor” 1925–1927. Amervid, 112 pp.

  • Xie, P. P., and P. A. Arkin, 1996: Analyses of global monthly precipitation using gauge observations, satellite estimates, and numerical model predictions. J. Climate, 9, 840858, doi:10.1175/1520-0442(1996)009<0840:AOGMPU>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Yeager, S. G., and W. G. Large, 2004: Late-winter generation of spiciness on subducted isopycnals. J. Phys. Oceanogr., 34, 15281547, doi:10.1175/1520-0485(2004)034<1528:LGOSOS>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Yu, L. S., 2011: A global relationship between the ocean water cycle and near-surface salinity. J. Geophys. Res., 116, C10025, doi:10.1029/2010JC006937.

    • Search Google Scholar
    • Export Citation
  • Fig. 1.

    A schematic of an open-ocean Smax. Solid lines represent isohaline contours with increasing contour thickness representing increased salinity. The dark (light) shaded region depicts that in the mixed layer eastward and poleward (westward and equatorward) of the sea surface salinity (SSS) max, ∂zS is stabilizing (destabilizing) because of the displacement of STUW with respect to the SSS maximum.

  • Fig. 2.

    The climatological annual cycle of surface heat flux (W m−2) over each Smax in CESM. The x axis is offset by six months for the SH Smax (top x-axis labels) with respect to the NH Smax (bottom x-axis labels) to account for the hemispheric phase shift of the seasons. The color legend is: NA Smax (light purple), SA Smax (light green), NI Smax (dark green), SI Smax (light blue), NP Smax (dark blue), and SP Smax (dark purple).

  • Fig. 3.

    The model climatological mean SLS. Black contours are the mean position of isohaline surface outcrops for each Smax. Gray contours are isohalines at 100 m depth. Each of the open-ocean Smax contours at 100 m depth is positioned westward and equatorward of the surface outcrop. The NI Smax surface and 100-m contours are at roughly the same position. White lines depict transects along which depth profiles are shown in Figs. 8 and 9.

  • Fig. 4.

    The annual cycle of outcrop extent anomaly (m2) from annual mean value for CESM (black) and MIMOC (light blue).

  • Fig. 5.

    Model climatological mean SLS (psu, shaded) and comparison of climatological mean Smax outcrops of the integration (black contours) and MIMOC (light blue). Black contours are CESM; white contours are MIMOC. In all cases, with the exception of NP Smax, the MIMOC outcrops are larger, with the largest discrepancy in Smax, where MIMOC outcrop area is nearly double that of the model. The white areas along the eastward and poleward edges of each MIMOC outcrop cover the paths of the Argo floats over the 2-yr periods described in section 4.

  • Fig. 6.

    The climatological annual cycle of the salinity budget (Sv psu) for each Smax. The net flux term is the sum of horizontal diffusion, eddy advection, and vertical diffusion.

  • Fig. 7.

    The simulation’s climatological mean vertical diffusive flux of salt though the base of the mixed layer (psu m s−1, shaded). Black and gray contours are isohaline contours as in Fig. 3. Light blue contours are regions where salt finger convection penetrates the base of the mixed layer. In the vicinity of the Smax, these regions of salt finger convection are on the eastward and poleward edges of Smax and correspond with regions of negative vertical diffusive salt fluxes through the mixed layer base (magenta-shaded regions).

  • Fig. 8.

    Transect plots through each NH Smax along the paths shown in Fig. 3. Climatological January planetary potential vorticity (m−1 s−1, shaded), isohaline bounding each Smax (white contours), and volumes where |Q| < 3.0 × 10−10 m−1 s−1 (magenta). Climatological mixed layer depth for January, March, and May is shown by mustard yellow contours ranging from dark–thick to light–thin. Climatological areas of salt finger convection for January, March, and May are shown by blue contours ranging from dark–thick to light–thin.

  • Fig. 9.

    As in Fig. 8, but for the SH with July, September, and November replacing January, March, and May.

  • Fig. 10.

    A simplified schematic of the seasonal deepening of the mixed layer. The solid lines depict isohaline contours. Dashed–dotted lines surrounding the white volumes depict low Q layers. Dotted lines surrounding dark gray volumes depict areas of salt finger convection.

  • Fig. 11.

    Seventy-three depth profiles from each of the Argo floats with data durations listed in Table 3. Depths where |Q| < 3.0 × 10−10 m−1 s−1 are depicted in magenta. Salt finger convection diffusivity as estimated by (2) (blue shading). Depth of the mixed layer depicted with a mustard yellow thick line.