Variability of Oceanic Barrier Layer and Compensated Layer
Total Page:16
File Type:pdf, Size:1020Kb
1
2
3
4
5
6 Observed subseasonal variability of oceanic barrier and compensated
7 layers
8
9
10 Hailong Liu, Semyon A. Grodsky, and James A. Carton 11
12
13 Submitted to Journal of Climate 14 December 10, 2008 15 16 17 18 19 20 21 22 23 24 25 26 27Department of Atmospheric and Oceanic Science 28University of Maryland, College Park, MD 20742 29 30 31Corresponding author: Semyon Grodsky ([email protected])
1 32Abstract 33 34A monthly gridded analysis of barrier layer/compensated layer thickness covering the period
351960-2007 is explored for evidence of subseasonal variability and the causes of these
36components of mixed layer variability. The strongest variability from monthly climatology is
37found in the subpolar latitudes of the North Atlantic in winter, where it reaches 100m. A
38temperature-salinity compensated layer in the eastern North Atlantic varies interannually in
39conjunction with a North Atlantic Oscillation-like pattern of anomalous sea level pressure. Thus
40thickening/thinning of this layer occurs in conjunction with the meteorological changes
41associated with strengthening/weakening of the Azores sea level high pressure center.
42
43In the winter subpolar North Pacific a salinity stratified barrier layer exists. While further south
44along the Kuroshio extension a compensated layer exists. The thickness of both of these layers
45varies from year to year by up to 60m. A significant component of this variability is a long-term
46trend of shrinkage of the barrier layer in the subpolar gyre and growth of the compensated layer
47along the Kuroshio Extension. This trend is shown to result from trends in meteorological
48variables. In the subpolar gyre the barrier layer has shrunk in response to a strengthening of the
49Aleutian low pressure system, the resulting strengthening of dry northerly winds, and a decrease
50of precipitation. In the Kuroshio Extension the compensated layer thickness has increased in
51response to this pressure system strengthening and related amplification of the midlatitude
52westerly winds, the resulting increase of net surface heat loss, and its effect on the temperature
53and salinity of the upper ocean water masses.
54
55The tropical Pacific and Atlantic both have barrier layers with variability that is less than 20m
2 2 56but is still important in influencing mixed layer properties. In the tropical Pacific west of 160oE,
57the barrier layer thickness changes by approximately 5m in response to 10 unit change in the
58South Oscillation Index. It thickens during La Ninas as a result of the abundant rainfall and thins
59during dry El Ninos. But east of 160oE variations in zonal advection of salt becomes as important
60as local rainfall in influencing the thickness of the barrier layer.
3 3 611. Introduction
62The ocean mixed layer is a nearsurface layer of fluid with quasi-uniform properties such as
63temperature, salinity, and density. The thickness of this mixed layer and its time rate of change
64both strongly influence the ocean’s role in air-sea interaction. However, the thickness of the
65nearsurface layer of quasi-uniform temperature, MLT, may differ from the thickness of the
66nearsurface layer of quasi-uniform density, MLD. MLT may be deeper than MLD when positive
67salinity stratification forms a barrier layer (BL=MLT-MLD) isolating the shallower and deeper
68levels of the mixed layer as was originally found in the western equatorial Pacific (Lukas and
69Lindstrom, 1991). Elsewhere MLT may be shallower than MLD when negative salinity
70stratification compensates for positive temperature stratification (or the reverse situation) to form
71a Compensated Layer (CL=MLD-MLT) (Stommel and Fedorov, 1967; Weller and Plueddemann,
721996). Changes in the seasonal thicknesses of BLs and CLs from one year to the next may cause
73corresponding changes in the role of the mixed layer in air-sea interaction by altering the
74effective depth of the mixed layer or the temperature of water at the mixed layer base (e.g., Ando
75and McPhaden, 1997). Here we examine the global historical profile data set covering the period
761960-2007 for corresponding year-to-year changes in the BL/CL thickness distribution.
77
78Four studies; Sprintall and Tomczak (1992), Tomczak and Godfrey (1994), de Boyer Montegut et
79al. (2007), and Mignot et al. (2007); have provided an observational description of the seasonal
80cycle of BL/CL distribution over much of the global ocean. BLs are a persistent feature of the
81tropics as well as high latitudes during winter. In high latitudes BLs occur where freshening in
82the near-surface is produced by excess precipitation over evaporation, river discharge, or ice
83melting (de Boyer Montegut et al., 2007) and may be most evident in regions where upward
4 4 84Ekman pumping ( wEk ) acts against the effects of vertical mixing such as occurs in the subpolar
85gyre (Kara et al., 2000). BLs are also observed equatorward of the salty subtropical gyres where
86subduction leads to a shallow subsurface salinity maximum (Sato et al., 2006). The North Pacific
87provides an example of this, the result of southward subduction of salty North Pacific
88Subtropical Mode Water below fresher tropical surface water (as originally suggested by
89Sprintall and Tomczak, 1992). BLs also occur in the Southern Ocean south of the Polar Front as a
90result of near surface freshening and weak thermal stratification (e.g. de Boyer Montegut et al.,
912004).
92
93BLs also occur seasonally in the tropics in regions of high rainfall and river discharge such as the
94Arabian Sea and Bay of Bengal, where layers as thick as 20-60m have been observed (Thadathil
95et al., 2008). Similarly, BLs occur in the western Equatorial Pacific under the high precipitation
96regions of the Intertropical Convergence Zone and South Pacific Convergence Zone (Lukas and
97Lindstrom, 1991) and in the western tropical Atlantic (Pailler et al., 1999; Ffield, 2007).
98
99Much less is known about subseasonal variations in BL/CL. In their examination of mooring
100time series Ando and McPhaden (1997) show that BLs do have interannual variability in the
101central and eastern Pacific, and conclude that the major driver is precipitation variability
102associated with El Nino. At 0oN 140oW, for example, the BL thickness increased from 10m to
10340m in response to the enhanced rains of the 1982-3 El Nino. In contrast, Mignot et al. (2007)
104suggest that changes in zonal advection of high salinity water in response to El Nino winds are
105important in regulating BLs in midbasin. In the west an observational study by Cronin and
106McPhaden (2002) is able to document the response of the mixed layer to intense westerly wind
5 5 107bursts and accompanying precipitation and show how these lead to both the formation and
108erosion of BLs.
109
110CLs in contrast may result from excess evaporation over precipitation, such as occurs in the
111subtropical gyres, or by differential advection where it leads to cooler fresher surface water
112overlying warmer saltier subsurface water (Yeager and Large, 2007). de Boyer Montegut et al.
113(2004) summarize several additional possible mechanisms of CL formation, such as subduction-
114induced advection, Ekman transport, slantwise convection and density adjustment. CLs are most
115prominent in the eastern subpolar North Atlantic and in the Southern Ocean (de Boyer Montegut
116et al., 2007). In the eastern North Atlantic a CL is formed by transport of the warm and salty
117North Atlantic Current above fresher colder subpolar water. Further east the North Atlantic
118Current splits into a northern branch comprising the Norwegian and Irminger Currents, and the
119southward Canary Current, all of which also develop CLs.
120
121In this study we build on previous observational examinations of the seasonal cycle of BL/CL
122development, to explore year-to-year variability. This study is made possible by the extensive 7.9
123million hydrographic profile data set contained in the World Ocean Database 2005 (Boyer et al.,
1242006) supplemented by an additional 0.4 million profiles collected as part of the Argo observing
125program. We focus our attention primarily on the Northern Hemisphere because of its higher
126concentration of historical observations.
127
1282. Data and methods
129This study is based on the combined set of temperature and salinity vertical profiles archived in
6 6 130the World Ocean Database 2005 for the period 1960-2004 and Argo floats from 1997 to 2007.
131The first dataset was adopted from the study of variability of the oceanic mixed layer by Carton
132et al. (2008) where issues of data quality control and processing are detailed.
133
134Mixed layer depth is defined here following Carton et al. (2008) which in turn combines the
135approaches of Kara et al. (2000) and de Boyer Montegut et al. (2004). The mixed layer depth is
136the depth at which the change in temperature or density from its value at the reference depth of
13710m exceeds a specified value (for temperature: T | 0.20 C |). This reference depth is
138sufficiently deep to avoid aliasing by the diurnal signal, but shallow enough to give a reasonable
139approximation of monthly SST. The value of T | 0.20 C | is chosen following de Boyer
140Montegut et al. (2004) as a compromise between the need to account for the accuracy of mixed
141layer depth retrievals and the need to avoid sensitivity of the results to measurement error. The
142absolute temperature difference instead of the negative temperature difference is used following
143Kara et al. (2000) in order to accommodate for temperature inversions that are widespread at
144high latitudes. The specified change in density used to define the density-based mixed layer
145depth follows the variable density criterion (e.g. Monterey and Levitus, 1997) to be locally
146compatible with the specified temperature value, (i.e. ( / T ) 0.20 C ). In this study the
147thickness (or width) of either a barrier layer or compensated layer is defined as a difference of
148isothermal mixed layer depth and isopycnal mixed layer depth, MLT-MLD. As a result of these
149definitions a positive MLT-MLD difference indicates the presence of a BL while a negative
150difference indicates the presence of a CL. We compute BL/CL thicknesses for each profile which
151are then passed through a subjective quality control to eliminate outliers and binned into 2o×2o×1
152month bins without attempt to fill in empty bins.
7 7 153
154In order to quantify the relative impact of temperature and salinity stratification within BLs and
155CLs we introduce a bulk Turner Angle, defined following Yeager and Large (2007) as:
1 1 1 156Tub tan [(T S) /(T S)], where / T and / S are the
157expansion coefficients for temperature, T , and salinity, S . In this study the changes in
158temperature and salinity T and S are defined between the top, zt min(MLT, MLD) , and
159the bottom, zb max(MLT, MLD) , of either a BL or CL based on analysis of individual vertical
160profiles. The bulk Turner angle is then evaluated from spatially binned values of T and S .
161Fig. 1 shows the correspondence between values of T , S , Tub , and the presence of BLs and
o 162CLs. Angles | Tub |<45 correspond to BLs stabilized by both, temperature and salinity ( T 0 ,
o 163 S <0). A perfect BL stabilized by salinity and homogeneous in T corresponds to Tub =-45 ,
o o 164while Tub =45 corresponds to a pure thermal stratification. Angles greater than 45 correspond to
165the most frequently occurring CLs where positive temperature stratification compensates for
166negative salinity stratification (mixed layer is saltier than the thermocline). Less frequently
o 167occurring CLs are observed at high latitudes below cool and fresh mixed layers (-90 o 168-72 ). The transition point is associated with the density ratio R T / S 0.5 or 1 0 169Tub tan (3) 72 . Thus bulk Turner angle provides an alternate way of displaying BL/CL 170distribution. 171 8 8 172We explore the role that surface forcing plays in regulating mixed layer properties through 173comparison of the BL/CL distribution to fluxes from the NCEP-NCAR reanalysis of Kalnay et 174al. (1996). Satellite QuikSCAT scatterometer winds (see Liu, 2002), which begin in mid-1999, 175are used to characterize the finer scale spatial patterns of wEk . To better characterize precipitation 176in the tropics, we also examine the Climate Prediction Center Merged Analysis of Precipitation 177(CMAP) of Xie and Arkin (1997), which covers the period 1979 -present. 178 179Limitations in the historical observational coverage limit our ability to explore the mechanisms 180giving rise to subseasonal BL/CL variability. In order to address these mechanisms we repeat our 181analysis on output from a 140-year long control simulation by the NOAA/Geophysical Fluid 182Dynamics Laboratory CM2.1 global coupled climate model forced with natural and 183anthroprogenic forcings beginning in year 1860 (after a 220 year spin-up period, Delworth et al., 1842006). This coupled simulation maintains full freshwater and heat flux balances at the air-sea 185interfaces. The ocean model component is based on Modular Ocean Model V4 numerics with 18610m vertical resolution in the upper 220 m and including vertical mixing with mixing rates 187determined based on the K-profile parameterization (KPP) mixed layer scheme of Large et al. 188(1994). Comparison of the simulated and observed mixed layer shows encouraging 189correspondence in both the geographic position and seasonal timing of BL/CL variability. 190 1913. Results 192We begin with a brief review of features of the mean and seasonal distribution. 193 1943.1. Time mean and seasonal patterns Global seasonal patterns of BL/CL display many features 9 9 195revealed by previous analyses (de Boyer Montegut et al., 2007). Throughout the year there are 196persistent BLs in the tropics in areas of high precipitation (Figs. 2a, 2b). BLs are thick under the 197Intertropical Convergence Zone and the South Pacific Convergence Zone. BLs are generally 198thickest in the western side of the tropical Pacific and Atlantic Oceans reflecting higher levels of 199rain as well as (in the case of the Atlantic) Amazon river discharge. BLs are thickest on the 200eastern side of the tropical Indian Ocean due to the presence of the Java and Sumatra high 201precipitation area (Qu and Meyers, 2005). Rainfall in the southern Intertropical Convergence 202Zone (Grodsky and Carton, 2003) may contribute to freshening of the mixed layer along 10oS in 203the South Atlantic during austral winter. In midlatitude BL/CLs occur mainly in winter and early 204spring. In boreal winter BLs exceeding 60 m are observed in the North Pacific subpolar gyre 205(Fig. 2a). Similarly thick BLs occur in the Atlantic Ocean north of the Gulf Stream. In both 206locations the BLs appear coincident with a seasonal increase of precipitation. 207 208Sea surface salinity (SSS) increases drastically across the Gulf Stream front leading to a switch 209from the BL regime north of the front to a CL regime south of the front (Fig. 2a). Thick CLs 210(thicker than 30m) are also observed along the Gulf Stream due to cross-frontal transport of low 211salinity water. And even thicker CLs (thicker than 60m) are observed further northeast along the 212path of the North Atlantic Current where its warm, salty water overlies cooler, fresher water. 213Interestingly, despite the presence of warm and salt western boundary currents in both the 214Atlantic and Pacific Oceans, the winter CLs are much less pronounced in the North Pacific than 215the North Atlantic. Explanation for the basin-to-basin difference likely lies in the higher surface 216salinity of the Atlantic (Fig. 2a) and consequently larger values of S (Fig. 3a). 217 10 10 218CLs are also evident in the southern subtropical gyres of the Pacific and Atlantic Oceans as well 219as the South Indian and Southwest Pacific Oceans (Fig. 2b) south of the 30oS SSS. The presence 220of CLs in this region reflects the northward advection of cold and fresh water which subducts 221under the water of the SSS maximum (Sprintall and Tomczak, 1993; note correspondence 222between CL in Fig. 2b and S in Fig. 3b). 223 224A similar mechanism is used to explain BL formation in subtropical gyres (e.g. Sato et al., 2006). 225In the southern Indian Ocean BLs north of 30oS form as salty water from the region of the SSS 226maximum subducts northward under relatively fresh surface water (Fig. 2b). Similar meridional 227dipole-like patterns with CLs to the south and BLs to the north of local subtropical SSS maxima 228are seen during austral winter in the South Pacific and the South Atlantic in the regions of 229downward wEk (Fig. 2b). In the Northern Hemisphere the areas of CLs in the North Pacific and 230Atlantic Oceans extend well north of the downward wEk regions (Fig. 2a) suggesting that 231horizontal transport of warm salt waters by the western boundary currents plays a role. 232 233Spatial patterns of BL/CL (Fig. 2) are in close correspondence with the spatial patterns of the 234vertical changes of salinity, S , (Figs. 3a, 3b). As expected, the barrier layers are distinguished 235by a stable salinity stratification, S 0 , where salinity increases downward below the mixed 236layer. In contrast, compensated layers have unstable salinity stratification, S 0 . As discussed 237above, regions of fresh mixed layer trace major areas of precipitation (like the Intertropical 238Convergence Zone) and river runoff (the Bay of Bengal). Different type of BL is observed on the 239equatorward flanks of subtropical salinity maxima. In these areas the ocean accumulates salt due 240to an excess of evaporation over precipitation. Here the equatorward propagation of subducted 11 11 241water produces meridional dipole-like BL/CL and S structures that are most pronounced 242during local winter in the Southern Hemisphere (Figs. 2b, 3b). Here barrier layers are observed 243to the north of SSS maxima where mixed layer tops saltier water below while compensated 244layers are observed to the south where mixed layer is saltier than thermocline. 245 246The spatial patterns of the bulk Turner angle (Figs. 3e,f) indicate that majority of CL cases are o 247associated with warm, salty water overlaying colder, fresher water (thus Tub >45 ). CLs widen 248during winter leading to higher values of Tub in winter than summer. When cold, fresh water o 249overlays warmer, saltier water a CL can also be formed where Tub <-72 . This latter type of 250density compensation is observed only in limited regions of the Labrador Sea and near Antarctica 251during winter. 252 2533.2 Subseasonal variability Subseasonal variability of BL/CL thickness is similar in amplitude to 254seasonal variability (compare Fig. 4 and Fig. 2). Subseasonal variability is stronger at higher 255latitudes reflecting weaker temperature stratification. The highest variability of BL/CL thickness 256of up to 100m occurs in the North Atlantic. It is also high in the western Atlantic and Pacific due 257to the impact of subseasonal variations in precipitation and river discharge. 258 259We next consider BL/CL thickness separated by season and roughly 15-year averaging periods 260(Fig. 5). During the first period 1960-1975 thick BLs are evident during local winter in the North 261Pacific, tropical western Pacific and western Atlantic, northern Indian Ocean, and Southern 262Ocean (the latter being evident even in austral summer). CLs during this early period appear 263primarily in the eastern North Atlantic. Little can be said about the existence of BLs in the 12 12 264Southern Ocean in austral winter due to the lack of data during this period. By the last period, 2651991-2007 several changes are evident. CLs have appeared in the subtropical North Pacific 266during winter and have strengthened in the eastern North Atlantic. Strong CLs are also evident 267on the northern side of the Circumpolar Current during austral winter (in fact these may have 268existed earlier but simply not been observed). By this last period the extensive winter BLs noted 269in the North Pacific in earlier decades had declined. We next examine monthly time series 270describing four regions: (1) BLs in the subpolar North Pacific, NP/BL box (2) CLs in the 271subtropical North Pacific, NP/CL box; (3) CLs in the eastern North Atlantic; and (4) BLs in the 272western equatorial Pacific (the regions outlined in Fig. 5c). 273 274The monthly time series of the northern subpolar North Pacific BL region and the subtropical CL 275region both show a long-term trend towards thinner BLs and thicker CLs interrupted by 276occasional reversals (Figs. 6a, 6b). Indeed, the subtropical CL region actually supported a 10- 27720m thick BL prior to 1980s. One direct cause of this change from BL to CL seems to be the 278gradual deepening of the late winter-spring mixed layer in the central North Pacific noted by 279Polovina et al. (1995) and Carton et al. (2008). This observed 20 m deepening into the cooler, 280fresher sub-mixed layer water has the effect of strengthening density compensation (the ‘spice 281injection’ mechanism is discussed by Yeager and Large, 2007). Carton et al. (2008) attribute the 282cause of mixed layer deepening to changes in the atmospheric forcing associated with the 283deepening of the Aleutian sea level pressure low after 1976. These changes led to strengthening 284of the midlatitude westerlies and the ocean surface heat loss in the North Pacific, hence the 285deepening of the mixed layer. The deepening of the mixed layer has opposite impacts on the 286width of CLs and BLs. Deepening widens CLs by injecting saltier water from the mixed layer 13 13 287into fresher thermocline. In contrast, stronger atmospheric forcing normally destroys near- 288surface BLs by enhancing mixing. These mechanisms likely explain the narrowing of BLs and 289the widening of CLs in the North Pacific during recent decades (Figs. 6a, 6b). 290 291The BL and CL thicknesses in the NP/BL and NP/CL boxes both seem to be associated with 292wind changes resulting from changes in the Aleutian pressure low (Fig. 7). Widening of CLs 293(that corresponds to the inverse of regression patterns shown in Fig. 7a) is linked to 294strengthening of westerly winds and latent heat loss (hence, deeper mixed layers) in the NP/CL 295box. Winds also strengthen in the NP/BL box to the north. This wind strengthening may reduce 296the BL due to enhanced mixing. But, in addition the strengthening dry northerly winds to the east 297of the low that decreases moisture transport from the south and thus reduces precipitation that is 298vital to maintain the BL (Fig. 7b). Limitations in the historical data record limit the extent to 299which these relationships can be identified based on data analysis. To overcome this data 300limitation we examine the processes regulating the mixed layer of the CM2.1 coupled model. 301 302We begin by comparing the simulated and observed winter/spring BL/CL distributions. In the 303simulation as in the observations a BL forms in the North Pacific north of approximately 35oN 304and a CL south of the BL region (compare Fig. 2a, Fig. 8a). The southern CL region can be 305divided into two zones. The first CL zone lies in the Kuroshio-Oyashio confluence at 306approximately 40oN in the northwest. The second CL zone extends across the North Pacific 307along approximately 30oN. This latter zone corresponds to the NP/CL region previously defined 308and is collocated with the subtropical SSS maximum and is bounded to the north by the northern 309edge of the downward wEk area (Fig. 8b). Subduction to the north of the subtropical SSS 14 14 310maximum produces negative salinity stratification in the region of SSS maximum while 311subduction at the SSS maximum produces positive salinity stratification to the south (see 312examples in Fig. 9). 313 314We next explore year to year changes in the NP/CL region by contrasting March 1953 when the 315simulated CL is the thickest of the 140 year record (25m) with March 1956 when the CL has 316been replaced by a 10m thick BL (Fig. 8c). The striking difference between these two years 317seems to be the different surface forcing. In 1953 downward wEk extends throughout the 318subtropics leading to a deep mixed layer (Fig. 9a). During this year the isohaline contours along 319180oE indicate an exaggeration of the normal situation in which there is southward advection of 320high salinity water below the fresh nearsurface layer. These conditions lead to a CL layer north o o 321of 25 N and BL layer to the south. In contrast, in 1956 wEk was upward north of 28 N trapping 322the freshwater introduced through precipitation in the surface layer and thus causing a shallow 323BL to form. The lag correlations in Fig. 8d confirm that CLs in the NP/CL region thicken in 324years when the mixed layer is anomalously deep, in line with the spice injection mechanism of 325Yeager and Large (2007), and BLs in the NP/BL thin as a result of the enhanced mixing 326associated with deep mixed layers. Thus deepening of the subtropical mixed layer may explain 327the trends observed in Figs. 6a, 6b. 328 329Correspondence between temporal changes of BL/CL width and atmospheric pressure patterns 330seen in both, observations and simulations in the north Pacific suggests that similar 331correspondence with local atmospheric patterns may take place in the north Atlantic. But the 332time series of variations in CL region of the North Atlantic is striking in showing primarily year- 15 15 333to-year variations (Fig. 6c) even though the winter-spring meteorology of this region does 334exhibit decadal variations frequently referred to as the North Atlantic Oscillation (Hurrell, 1995). 335We suspect that some of this variability is the result of the smaller spatial scales of variability 336relative to the North Pacific which has been insufficiently sampled by the historical observing 337systems. During the past ten years when the data coverage has been increasingly extensive (due 338to implementation of the Argo floats) CL thickness in the North Atlantic CL region has 339undergone 15m variations from year-to-year with thick layers in 1999, 2002-3, and 2006-7 and 340thinner layers during the other years (Fig. 10a). This pattern of CL variation is related to the 341strength of the subtropical high in sea level pressure (Fig. 10a). Indeed, a regression of sea level 342pressure, latent heat flux, and winds on CL thickness in the North Atlantic CL box shows basin- 343scale relationships (Fig. 10b). Vertical expansion of CLs (more negative CL width) corresponds 344to NAO-like pattern of the atmospheric pressure with stronger Azores high and deeper Iceland 345low. This pressure pattern strengthens winds and latent heat losses in the east while weakens 346them in the central and western subtropics (Fig. 10b). Anomalously strong winds and net surface 347heat losses suggest deepening of the mixed layer that, in turn, expands width of CLs by the spice 348injection mechanism. This is in line with 2002-3 and 2006-7 data when anomalously thick CLs 349in the NA/CL box occur in phase with anomalously deep mixed layers. But in 1999 this 350relationship is opposite (Fig. 10a). This may suggest that alternative mechanisms including salt 351advection play a role. Indeed, in the north-east it is clear that advection plays a central role in 352variability of the mixed layers (Nilsen and Falck, 2006 show this process at work in the 353Norwegian Sea). 354 355Finally, we consider year-to-year variability of the BL region of the tropical Pacific. BL thickness 16 16 356in the western Pacific west of 160oE, which is generally 10-20m, varies in thickness coherently 357with the value of the Southern Oscillation Index (Fig. 11c). The eastward shift of precipitation 358during El Nino decreases precipitation west of 160oE by 20% or 1 mm/dy (in response to a 20 359unit decrease of the SOI) and as a result BLs shrink by 20% or 2-5m (Figs. 11a, 11b). Thus, in 360this region variations in BL thickness respond primarily to surface freshwater flux. East of 160oE 361BLs are thinner and occur at much shallower depths so that they are affected more strongly by 362solar radiation, advection and entrainment. 363 3644. Conclusions 365This study examines subseasonal changes of barrier and compensated layers (BLs and CLs) 366based on analysis of profiles of temperature and salinity covering the years 1960-2007. Because 367of data limitations we focus mainly on the Northern Hemisphere and tropics. The processes that 368regulate subseasonal variability of BL/CL thickness are similar to those which regulate their 369seasonal appearance: fluctuations in surface freshwater flux, Ekman pumping, and horizontal 370advection. Thus, the spatial distribution of subseasonal variability reflects aspects of the 371subseasonal variability of these forcing terms. 372 373The most striking interannual variability of BL width occurs in the rainy tropical oceans with 374deep thermoclines such as the Arabian Sea and Bay of Bengal in the Indian Ocean and the 375western tropical Pacific under the Intertropical Convergence Zone and the South Pacific 376Convergence Zone. Precipitation in these regions varies strongly interannually. During high 377precipitation years the mixed layers in these regions show capping fresh layers and thick BLs. In 378contrast, during low precipitation years mixed layer salinities increase and BL thickness 17 17 379decreases. Thus, in the western Equatorial Pacific between 120oE and 160oE, the BL thickens by 3802m to 5m during La Nina while during the El Nino the BL thins by a similar amount. In the 381central and eastern equatorial Pacific variations of BL thickness are not nearly as well correlated 382with the phase of ENSO. We think this happens because in the central basin salt advection plays 383an important role in regulating BL thickness as well, while further east the mixed layer shallows 384and mixed layer entrainment variability becomes increasingly important. 385 386In the subtropics and midlatitudes during late winter-spring we find alternating regions of CLs 387and BLs in the seasonal climatology. The northern tropics of both the Pacific and Atlantic (the 388southern edge of the subtropical gyres) show broad regions of BLs where salty subtropical 389surface water formed further north has subducted, advected equatorward, and affected the water 390properties of the winter mixed layer. Within the evaporative subtropical North Pacific and eastern 391North Atlantic we find CLs resulting from mixed layers with positive temperature stratification 392but negative salinity stratification. In the eastern subpolar North Atlantic this CL merges with a 393subpolar CL resulting from negative temperature stratification but positive salinity stratification. 394In the North Pacific, strengthening of the Aleutian pressure low during successive winters, thus 395strengthening the midlatitude westerly winds has led to deeper mixed layers, cooler SSTs, and a 396long-term increase in the thickness of this CL. Further north in the North Pacific a supolar BL 397region is evident, which does not have a counterpart in the North Atlantic. The same changes in 398meteorology which include strengthening of the Aleutian pressure low have led to an increase in 399dry northerly winds which in turn cause a thinning of this BL from ~40m before 1980s to ~ 20m 400afterwards. 401 18 18 402In the North Atlantic the broad area of CL also undergoes fluctuations in thickness associated 403with surface meteorology. In this basin the fluctuations in winter thickness also seem to vary in 404response to meteorological changes associated with changes of the Azores sea level high 405pressure. Estimation of the timescale of variability is more difficult in this basin as the analysis is 406more noisy. Here our examination focuses on the changes that are apparent in the shorter 1997- 4072007 period when the data coverage is better. During this period the thickness of the winter- 408spring CL in the northeastern Atlantic undergoes 4-5 year timescale fluctuations from 45m to 40960m in phase with the strengthening/weakening of the Azores high. 410 411Determining the basin-scale BL/CL structure tests the limits of the historical observing system. 412Further progress in understanding BL/CL variability and its role in air-sea interactions will likely 413require further exploration of models that provide reasonable simulations of observed variability. 414 415Acknowledgements We gratefully acknowledge the Ocean Climate Laboratory of the National 416Oceanographic Data Center/NOAA and the Argo Program data upon which this work is based. 417Support for this research has been provided by the National Science Foundation (OCE0351319) 418and the NASA Ocean Programs. 19 19 419References 420Ando, K., and M.J. McPhaden, 1997: Variability of surface layer hydrography in the tropical 421 Pacific Ocean. J. Geophys. Res., 102, 23063-23078. 422Boyer, T.P., J.I. Antonov, H.E. Garcia, D.R. Johnson, R.A. Locarnini, A.V. Mishonov, M.T. 423 Pitcher, O.K. Baranova, and I.V. Smolyar, 2006: World Ocean Database 2005. S. Levitus, 424 Ed., NOAA Atlas NESDIS 60, U.S. Government Printing Office, Washington, D.C., 190 pp., 425 DVDs. 426Carton, J.A., S.A. Grodsky, and H. Liu, 2008: Variability of the Oceanic Mixed Layer, 1960– 427 2004. J. Climate, 21, 1029–1047. 428Cronin, M. F., and M. J. McPhaden, 2002: Barrier layer formation during westerly wind bursts. 429 J. Geophys. Res., 107(C12), 8020, doi:10.1029/2001JC001171. 430de Boyer Montégut, C., G. Madec, A. S. Fischer, A. Lazar, and D. Iudicone, 2004: Mixed layer 431 depth over the global ocean: An examination of profile data and a profile-based climatology. 432 J. Geophys. Res., 109, C12003, doi:10.1029/2004JC002378. 433de Boyer Montégut, C., J. Mignot, A. Lazar, and S. Cravatte, 2007: Control of salinity on the 434 mixed layer depth in the world ocean: 1. General description. J. Geophys. Res., 112, C06011, 435 doi:10.1029/2006JC003953. 436Delworth, T.L., and Coauthors, 2006: GFDL's CM2 Global Coupled Climate Models. Part I: 437 Formulation and Simulation Characteristics. J. Climate, 19, 643–674. 438Ffield, A., 2007: Amazon and Orinoco River Plumes and NBC Rings: Bystanders or Participants 439 in Hurricane Events? J. Climate, 20, 316–333. 440Grodsky, S. A., and J. A. Carton, 2003: Intertropical convergence zone in the South Atlantic and 441 the equatorial cold tongue. J. Climate, 16(4), 723-733, 20 20 442Hurrell, J. W., 1995: Decadal trends in the North-Atlantic oscillation - regional temperatures and 443 precipitation. Science, 269, 676-679. 444Kalnay, E., & co-authors, 1996: The NCEP/NACR 40-year reanalysis project. Bull. Amer. 445 Meteor. Soc., 77, 437-471. 446Kara, A. B., P. A. Rochford, and H. E. Hurlburt, 2000: Mixed layer depth variability and barrier 447 layer formation over the North Pacific Ocean. J. Geophys. Res., 105(C7), 16,783–16,801. 448Large, W., J. McWilliams, and S. Doney, 1994: Oceanic Vertical Mixing: A Review and a Model 449 with a Nonlocal Boundary Layer Parameterization. Rev. Geophys., 32(4), 363-403. 450Liu, W.T., 2002: Progress in scatterometer application. J. Oceanogr., 58(1), 121-136. 451Lukas, R., and E. Lindstrom, 1991: The mixed layer of the western equatorial Pacific Ocean. J. 452 Geophys. Res., 96 (Supplement), 3343 –3357. 453Mignot, J., C. de Boyer Montégut, A. Lazar, and S. Cravatte, 2007: Control of salinity on the 454 mixed layer depth in the world ocean: 2. Tropical areas. J. Geophys. Res., 112, C10010, 455 doi:10.1029/2006JC003954 456Monterey, G., and S. Levitus, 1997: Seasonal Variability of Mixed Layer Depth for the World 457 Ocean, NOAA Atlas NESDIS 14, 100 pp., Natl. Oceanic and Atmos. Admin., Silver Spring, 458 MD. 459Nilsen, J. E. Ø., and E. Falck, 2006: Variations of mixed layer properties in the Norwegian Sea 460 for the period 1948–1999. Prog. Oceanogr., 70, 58–90. 461Pailler, K., B. Bourlès, and Y. Gouriou, 1999: The Barrier Layer in the Western Tropical Atlantic 462 Ocean. Geophys. Res. Lett., 26(14), 2069-2072. 21 21 463Polovina, J.J., G.T. Mitchum, and G.T. Evans, 1995: Decadal and basin-scale variation in MLD 464 and the impact on biological production in the central and North Pacific, 1960-88. Deep-Sea 465 Res., 42, 1701-1716. 466Qu, T., and G. Meyers, 2005: Seasonal variation of barrier layer in the southeastern tropical 467 Indian Ocean. J. Geophys. Res., 110, C11003, doi:10.1029/2004JC002816. 468Sato, K., T. Suga, and K. Hanawa, 2006: Barrier layer in the North Pacific subtropical gyre. 469 Geophys. Res Lett, 31, L05301, doi:10.1029/2003GL018590. 470Sprintall, J., and M. Tomczak, 1992: Evidence of the Barrier Layer in the Surface Layer of the 471 Tropics. J. Geophys. Res., 97(C5), 7305–7316. 472Sprintall, J., and M. Tomczak, 1993: On the formation of central water and thermocline 473 ventilation in the Southern Hemisphere. Deep Sea Res., Part I, 40, 827–848. 474Stommel, H., and K. N. Fedorov, 1967: Small-scale structure in temperature and salinity near 475 Timor and Mindanao. Tellus, 19, 306-325. 476Thadathil, P., P. Thoppil, R.R. Rao, P.M. Muraleedharan, Y.K. Somayajulu, V.V. Gopalakrishna, 477 R. Murthugudde, G.V. Reddy, and C. Revichandran, 2008: Seasonal Variability of the 478 Observed Barrier Layer in the Arabian Sea. J. Phys. Oceanogr., 38, 624–638. 479Tomczak, M., and J. S. Godfrey, 1994: Regional Oceanography: An Introduction, Pergamon 480 Press, New York, 422 pp. (pdf edition is available at 481 http://gyre.umeoce.maine.edu/physicalocean/Tomczak/regoc/pdfversion.html) 482Weller, R. A., and A. J. Plueddemann, 1996: Observations of the vertical structure of the oceanic 483 boundary layer. J. Geophys. Res., 101, 8789–8806. 22 22 484Xie, P., and P.A. Arkin, 1997: Global Precipitation: A 17-Year Monthly Analysis Based on Gauge 485 Observations, Satellite Estimates, and Numerical Model Outputs. Bull. Amer. Meteor. Soc., 486 78, 2539–2558. 487Yeager, S.G., and W.G. Large, 2007: Observational Evidence of Winter Spice Injection. J. Phys. 488 Oceanogr., 37, 2895–2919. 23 23 489Figure captions. 490 Figure 1. Climatological cold season barrier layer/compensated layer thickness versus the 491bulk Turner angle evaluated using temperature ( T ) and salinity ( S ) difference between the 492top and bottom of a barrier or compensated layer. Vertical bars show the mean and the standard 493deviation for consecutive 22.50 interval of the Turner angle. Grey dots are the cold season data 494that aggregate January-March data from the Northern Hemisphere and July-September data from 495the Southern Hemisphere. The Turner angle range -720 to 450 corresponds to BL. CL occurs 496outside this interval. 497 Figure 2. Climatological (a) January-March and (b) July-September barrier layer width 498(positive) and compensated layer width (negative). Climatological sea surface salinity ( 35 psu 499solid contours, below 35 psu dashed contours) and downward wind-driven vertical velocity 500(hatching) are overlain. Salinity data are from Boyer et al. (2006). 501 Figure 3. (a,b) Salinity ( S ) and (c,d) temperature ( T ) difference between the top zt = 502min[MLT, MLD] and the bottom zb =max[MLT, MLD] of barrier/compensated layer, (e,f) bulk 503Turner angle calculated from S and T between the same two depths. (left) January - March 504(JFM) values, (right) July - September (JAS) values. Turner angles in the range from -72o to 45o 505corresponds to barrier layers, while compensated layers occur outside this range 506 507 Figure 4. Standard deviation (STD) of monthly anomalous BL/CL thickness. 508 Figure 5. Quasi-decadal mean barrier layer (positive) and compensated layer (negative) 509thickness in (left) northern winter and (right) austral winter. Rectangles in (c) show locations of 510the North Pacific barrier layer box (NP/BL 140oE-140oW, 50o-60oN), North Pacific compensated 511layer box (NP/CL 140oE-160oW, 25o-35oN), equatorial Pacific barrier layer box (EP/BL 120oE- 24 24 512160oE, 5oS-5oN), and North Atlantic compensated layer box (NA/CL 30oW-10oW, 40o-60oN) 513 Figure 6. Box averaged BL/CL width in (a) North Pacific barrier layer region, (b) North 514Pacific compensated layer region, and (c) North Atlantic compensated layer region. Thin lines 515are monthly average values, bold lines are January-March average values. JFM values are shown 516if at least 10 data points are available. Box locations are shown in Fig. 5c. 517 Figure 7. Time regression of JFM barrier layer thickness in (a) North Pacific compensated 518layer box on latent heat flux (Wm-2/m, shading), winds (ms-1/m, arrows), and sea level pressure 519(mbar/m); (b) North Pacific barrier layer box on precipitation (mm h-1/m). See also Fig.5c for 520box locations. Atmospheric parameters are from the NCEP/NCAR Reanalysis. 521 Figure 8. CM2.1 model simulations in the north Pacific. (a) March climatological BL/CL 522width with climatological sea surface height (SSH, CI=0.2m) and sea surface salinity (SSS, bold 523contours) overlain. CL and BL boxes are drawn in red and blue, respectively. (b) Zonally 524averaged March climatological Ekman pumping (positive – downward is shaded gray). (c) Time 525series of the box-averaged CL width. Two contrasting model years with thick (1953) and thin 526(1956) CL are marked by ‘o’. (d) Lagged correlation of anomalous (blue) box-averaged CL 527width and MLD, (red) box-averaged BL width and MLD. (e) Time regression of March box- 528averaged CL width on wind stress and Ekman pumping (mdy-1/m). 529 Figure 9. Latitude-depth section of salinity along 180oE for March of the two contrasting 530model years with thick (1953) and thin (1956) CL. Ekman pumping (positive – downward is 531shaded gray) is shown against the right-hand vertical axes. 532 Figure 10. (a) JFM compensated layer width (solid) and monthly mixed layer depth 533(compensated layer depth range between MLT and MLD is shaded) averaged in the North 534Atlantic compensated layer box, JFM anomalous mean sea level pressure (dashed) averaged 25 25 53540oW-30oW, 30oN-40oN. (b) Time regression of the 1997-2007 JFM compensated layer width on 536latent heat flux (Wm-2/m, shaded), anomalous mean sea level pressure (mbar/m, contours), and 537winds (ms-1/m, arrows). 538 Figure 11. Lag regression of SOI index on 5oS-5oN averaged (a) anomalous barrier layer 539width, (b) precipitation. Time mean BL width (solid), salinity (dashed), and precipitation are 540shown against right hand vertical axes. Precipitation is Xie and Arkin (1997). (c) Time series of 541annual running mean SOI (shaded) and anomalous BL width averaged in the western equatorial 542Pacific box (120oE-160oE, 5oS-5oN). See Fig. 5c for the box location. 26 26