1. Introduction
The eastern tropical Pacific plays a major role in global climate variability. Sea surface temperature (SST) variability of the cold tongue couples to the tropical atmosphere via the Bjerknes feedback that can amplify initially small anomalies; it is a fundamental part of the El Niño phenomenon. Cold tongue SST variability thus has global impacts on weather and climate (Alexander et al. 2002; Jong et al. 2016; Yeh et al. 2018; Claar et al. 2018; Cai et al. 2019).
Vertical exchange, from either advection along sloping isopycnals or diapycnal (cross-isopycnal) motion and mixing, can have a large impact on overlying SST by transporting cold waters to the surface. While nearly continuous measurements of horizontal velocities are available at four moorings along the equator in the Pacific from the Tropical Atmosphere Ocean (TAO) buoy array/the Tropical Pacific Observing System (TPOS) [McPhaden et al. (1998, 2010); also see recent review of TPOS by Smith et al. (2019)], vertical velocity w is challenging to measure directly. Vertical velocity w can be derived from divergence estimates over large control boxes (Wyrtki 1981; Bryden and Brady 1985; Johnson et al. 2001; Meinen et al. 2001) or from moored arrays (Weisberg and Qiao 2000), through a calculation involving small differences between large numbers with resulting large uncertainty (Johnson et al. 2001). Bryden and Brady (1985) use a diagnostic model based on hydrographic sections to derive geostrophic velocity to quantify the three-dimensional circulation in the upper equatorial Pacific. Meinen et al. (2001) use geostrophy and Ekman balance to calculate vertical velocities over a large box spanning the cold tongue. The latter two studies conclude that the diabatic (cross-isothermal) part of the vertical velocity (wci) is a modest fraction of the total vertical velocity, though the error bars in Meinen et al. (2001) demonstrate the large uncertainty of these observational estimates.
Diabatic processes in the eastern Pacific have been shown to play an important role in driving warm water volume (WWV) variability through diabatic upwelling across the 20°C isotherm, which serves as a proxy for the depth of the thermocline (Meinen and McPhaden 2001; Lengaigne et al. 2012; Huguenin et al. 2020). WWV variability in turn is correlated to the Niño-3.4 index at lead times of 3–6 months, making it important in ENSO prediction (Meinen and McPhaden 2000; Clarke et al. 2007). Not all studies agree on the contribution of wci to WWV variability, however. Brown and Federov (2010) and Bosc and Delcroix (2008) argue that wci varies little interannually. With the current observing array, this important part of the circulation is difficult to constrain. Here, we make use of a high-resolution global ocean model to characterize the diabatic processes and their seasonal and subseasonal variability.
In our previous work we examined the modulation of wci with El Niño–Southern Oscillation (ENSO) using the saved heat budget from the high-resolution global ocean model Parallel Ocean Program version 2 (POP2) at 0.1° horizontal resolution (Smith et al. 2010). We found that the diabatic upwelling across the thermocline is dominated by the vertical heat flux divergence produced by turbulent mixing. Interannual variations in vertical shear between the Equatorial Undercurrent (EUC) and the South Equatorial Current (SEC) results in diabatic upwelling being almost entirely shut down during El Niño and strengthened during La Niña (Deppenmeier et al. 2021).
The tropical Pacific also displays strong seasonal and subseasonal variability that interact with the longer time scales of ENSO (Legeckis 1977; Chelton et al. 2000; Fiedler and Talley 2006; Kessler 2006; Willett et al. 2006). The seasonal cycle of SST in the eastern equatorial Pacific is largely driven by the competing influences of surface fluxes and ocean vertical mixing (Wang and McPhaden 1999; Moum et al. 2013). Using moored platforms for high-frequency temperature measurements (χpods) that can be processed to infer the dissipation rate of temperature variance, Moum et al. (2013) show that, on the equator at 140°W, surface heating dominates between February and May, leading to rising SST [Fig. 1a, reproduced from Fig. 3b of Moum et al. (2013)]. In contrast, between July and October, turbulent cooling from below exceeds surface heating, leading to a negative SST tendency. Moum et al. (2013) attribute the observed prolonged cooling in September–November to horizontal advection processes not captured by their simplified heat budget depicted in Fig. 1a.

The 2005–11 mean seasonal cycle of SST (black), turbulent heat flux between 20 and 60 m (Jq, blue), and surface heat flux (SHF, red) at 0°, 140°W (a) from χPOD data (contains gaps) (see Moum et al. 2013, their Fig. 1f) and (b) from the POP2 model. The error bars on Jq in (a) correspond to two different ways of estimating the turbulent heat flux, whereas the bars in (b) show the standard deviation.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

The 2005–11 mean seasonal cycle of SST (black), turbulent heat flux between 20 and 60 m (Jq, blue), and surface heat flux (SHF, red) at 0°, 140°W (a) from χPOD data (contains gaps) (see Moum et al. 2013, their Fig. 1f) and (b) from the POP2 model. The error bars on Jq in (a) correspond to two different ways of estimating the turbulent heat flux, whereas the bars in (b) show the standard deviation.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
The 2005–11 mean seasonal cycle of SST (black), turbulent heat flux between 20 and 60 m (Jq, blue), and surface heat flux (SHF, red) at 0°, 140°W (a) from χPOD data (contains gaps) (see Moum et al. 2013, their Fig. 1f) and (b) from the POP2 model. The error bars on Jq in (a) correspond to two different ways of estimating the turbulent heat flux, whereas the bars in (b) show the standard deviation.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
Mixing on the equator is also heavily modulated on subseasonal time scales connected to the passage of tropical instability waves (TIWs) (Lien et al. 2008; Moum et al. 2009; Inoue et al. 2012; Holmes and Thomas 2015; Inoue et al. 2019). TIWs result from baroclinic and barotropic instabilities in the eastern equatorial current system and can be detected (e.g., in SST or meridional velocity) as long westward traveling waves with periods between 13 and 40 days (Legeckis 1977; Miller et al. 1985; Halpern et al. 1988; Qiao and Weisberg 1995). TIW presence is strongly modulated by the seasonal cycle (strongest TIW occur during boreal fall and winter, and weaker TIWs in boreal spring) and ENSO state (with stronger TIWs presence during La Niña than during El Niño conditions). TIWs can be detected by examining variability in the appropriate frequency band for a variable that would be influenced by the passage of the TIWs, such as variations in meridional velocity. Lien et al. (2008) use a numerical model and Lagrangian floats to show that the leading edge of TIWs enhances turbulent heat fluxes at the base of surface mixed layer by orders of magnitude. Using χpod measurements, Moum et al. (2009) show that TIWs increase vertical shear in the ocean, resulting in instability and turbulence that cools the surface and warms the EUC waters. Inoue et al. (2012) demonstrate strong modulation of turbulent mixing depending on the phase of a passing TIWs which Holmes and Thomas (2015) attribute to increased shear from horizontal vortex stretching. Cherian et al. (2021) demonstrate that a similar mechanism is important off the equator.
Here, we investigate the seasonal cycle and subseasonal variability of wci and its physical drivers in and across the eastern Pacific cold tongue in a high-resolution ocean model.
2. Datasets and methodology
a. Model simulation and observational data
We use output from a global 0.1° horizontal resolution POP2 simulation forced with interannually varying JRA55-do surface fluxes with 3-hourly temporal resolution between 1958 and 2018 (Kobayashi et al. 2015; Tsujino et al. 2018) described in Bryan and Bachman (2015) and Deppenmeier et al. (2021). The vertical resolution is 10 m in the upper 200 m and then increases toward the bottom. The full heat budget averaged over five days during runtime is available for the years 1983–2018. We also use TAO mooring data (Hayes et al. 1991; McPhaden et al. 1998, 2010) provided by the Global Tropical Moored Buoy Array Project Office of NOAA/Pacific Marine Environmental Laboratory for validation of the model simulation along three longitudes that span the cold tongue (170°, 140°, 110°W).
Comparing the observational estimate of heat fluxes (Fig. 1a) to POP2 (Fig. 1b) reveals differences in both the surface heating flux (red) as well as the (cooling) turbulent heat flux (blue). The surface heat flux differences stem from differing estimates of the solar penetration flux at the base of the mixed layer and the difference in the net longwave flux (Figs. S1 and S2 in the online supplemental material). This includes a bias in the diagnosis of the solar penetration term, which arises from the POP2 estimates being based on 5-day average output data while the Moum et al. (2013) estimate is based on hourly data. This is not a bias in the model itself, which includes a resolved diurnal cycle of solar heating, but is a limitation of the sampling available for a posteriori analysis.
A comparison of the model time mean and interannual time scale (1983–2018) vertical structure of temperature and zonal velocity to observations was made in Deppenmeier et al. (2021). The model mean thermocline and the EUC are biased deep, with larger biases at 110°W than at 140° or 170°W. We found larger biases in the temperature structure (deep thermocline) during La Niña than during El Niño. The model represented the ENSO-cycle variations (anomalies in thermocline and EUC depth and strength) well, again with better representation at 140° or 170°W than at 110°W.
For the seasonal cycle a comparison of the vertical structure of the model zonal current and temperature to observations shows that the model performs better in the central Pacific than farther east (Fig. 2). For the analysis here, we find that the vertical shear above the EUC (above the eastward peak in the zonal velocity) is reasonably well simulated throughout the year, even at 110°W.

Seasonal cycle of (top) zonal velocity and (bottom) temperature from POP2 (solid lines, 1988–2018 continuous) and the TAO array (dashed lines, 1988–2018, different gaps for different mooring locations) for three different locations.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Seasonal cycle of (top) zonal velocity and (bottom) temperature from POP2 (solid lines, 1988–2018 continuous) and the TAO array (dashed lines, 1988–2018, different gaps for different mooring locations) for three different locations.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
Seasonal cycle of (top) zonal velocity and (bottom) temperature from POP2 (solid lines, 1988–2018 continuous) and the TAO array (dashed lines, 1988–2018, different gaps for different mooring locations) for three different locations.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
POP2 simulates seasonally varying TIW activity, but the TIWs’ kinetic energy is reduced by a factor of 2 compared to observations (Figs. S3 and S4). This is partly due to weak meridional shear north of the equator due to a weak North Equatorial Countercurrent and a weak South Equatorial Current. These biases are attributed to shortcomings in the JRA-55do forcing product (Sun et al. 2019). Additionally, in a forced ocean model the occurrence of TIW does not temporally and spatially align with the imprint the observed TIW left on the forcing field. This misalignment of air–sea fluxes damps the SST anomalies and thus the TIW (Renault et al. 2020; Rai et al. 2021).
b. Water mass transformation framework
The density and specific heat of seawater is given by ρ0 and cp, respectively, and I is the downward radiation flux into the water column as a function of depth z. In the model, the vertical structure of I is based on climatological chlorophyll levels and calculated according to Ohlmann (2003) (see Fig. S2). The parameterized vertical and lateral diffusive heat fluxes are abbreviated with J and Hdiff, respectively.
Evaluating Eq. (1) using 5-day-averaged terms of the model full heat budget allows calculation of the diabatic component of the vertical circulation (wci), and importantly also allows us to attribute the cross-isothermal velocities to physical processes such as vertical mixing (
The attribution of wci to physical processes such as solar penetration and vertical mixing depends on the divergence (vertical derivative) of those heat fluxes. Thus, peak wci does not necessarily coincide with peak heat flux, but instead with its maximum vertical gradient. For solar penetration, the downward flux reduces monotonically with depth: its divergence maintains one sign, and

Schematic of representative cross-equatorial structure in the cold tongue. Isotherms are represented by the blue lines, with darker shades of blue indicating colder water. Solar heat flux (I; yellow arrows and profile) is maximum at the surface and decreases with depth. Heat flux due to turbulent mixing (J; blue swirls and profile) is large in the high shear, high temperature gradient region above the EUC, then decreases in the surface mixed layer. The wci is proportional to the vertical gradient of these fluxes. Red arrows show total wci (in the direction of the temperature gradient). Turbulent mixing, solar penetration, and wci are active on and across the equator.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Schematic of representative cross-equatorial structure in the cold tongue. Isotherms are represented by the blue lines, with darker shades of blue indicating colder water. Solar heat flux (I; yellow arrows and profile) is maximum at the surface and decreases with depth. Heat flux due to turbulent mixing (J; blue swirls and profile) is large in the high shear, high temperature gradient region above the EUC, then decreases in the surface mixed layer. The wci is proportional to the vertical gradient of these fluxes. Red arrows show total wci (in the direction of the temperature gradient). Turbulent mixing, solar penetration, and wci are active on and across the equator.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
Schematic of representative cross-equatorial structure in the cold tongue. Isotherms are represented by the blue lines, with darker shades of blue indicating colder water. Solar heat flux (I; yellow arrows and profile) is maximum at the surface and decreases with depth. Heat flux due to turbulent mixing (J; blue swirls and profile) is large in the high shear, high temperature gradient region above the EUC, then decreases in the surface mixed layer. The wci is proportional to the vertical gradient of these fluxes. Red arrows show total wci (in the direction of the temperature gradient). Turbulent mixing, solar penetration, and wci are active on and across the equator.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
In the eastern Pacific cold tongue, vertical mixing is closely connected to the destabilizing shear between the eastward EUC and overlying westward SEC (Fig. 2) (Sun et al. 1998; Smyth and Moum 2013; Smyth et al. 2013). The resulting heat flux from mixing J also depends on the temperature gradient, which is large above the thermocline but much smaller in the roughly 50-m-thick surface mixed layer (Fig. 2, noting that instantaneous mixed layer thickness can be larger than these seasonal averages). Thus J, owing to vertical mixing, peaks near the base of the mixed layer so its divergence takes both signs (Figs. 3b,c; see also Fig. 7 in Deppenmeier et al. 2021). The
The combined impact of
Since wci is by definition perpendicular to isotherms, in the mixed layer, where the temperature gradient is more horizontal (meridional) rather than vertical; wci can be dominantly in the poleward direction (Fig. 3a).
3. Results
The seasonal evolution of SST and its drivers in the model is similar in character to that inferred from the Moum et al. (2013) measured temperature variance at 140°W, with timing of the transition from a surface flux dominated heat budget to a vertical turbulent flux dominated heat budget occurring at the same time of year as in the observations (cf. Figs. 1a and 1b). During September–November, however, in the model strong turbulent heat fluxes occurs (cooling) while the SST is relatively constant (Fig. 1b) while observational results suggest weak turbulent heat fluxes during this time (Fig. 1a). Both surface and turbulent terms remain large in the model representation at this time (Fig. 1b). Moum et al. (2013) suggested that horizontal advection might be significant at this time, but noted the difficulty in estimating horizontal SST gradients consistent with the measured currents that are reliable below 25 m in depth.
a. The seasonal cycle of wci
The diabatic upwelling wci in the tropical Pacific cold tongue thermocline exhibits a distinct seasonal cycle (Fig. 4). The seasonally varying thermocline wci, defined as the monthly average wci between 20° and 24°C, indicates that the strongest water mass transformation from below-thermocline to near-surface waters takes place between boreal summer and winter, and is reduced in boreal spring.

Seasonally stratified mean thermocline diabatic upwelling (approximated by averaging over 20°–24°C). Stippled vertical lines indicate the locations of the TAO moorings at 155°, 140°, 125°, and 110°W.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Seasonally stratified mean thermocline diabatic upwelling (approximated by averaging over 20°–24°C). Stippled vertical lines indicate the locations of the TAO moorings at 155°, 140°, 125°, and 110°W.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
Seasonally stratified mean thermocline diabatic upwelling (approximated by averaging over 20°–24°C). Stippled vertical lines indicate the locations of the TAO moorings at 155°, 140°, 125°, and 110°W.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
In all seasons, the largest values of wci in the thermocline are found east of 155°W (Fig. 4) and are concentrated near the equator. The latitudinal extent of these large values of wci has large temporal variability. Notably, wci is confined closely to the equator in April, whereas it extends 2° away from the equator in both hemispheres in the other seasons.
While the region of large wci in the thermocline is much reduced in April (Fig. 4b), there is equatorial diabatic upwelling in the water column throughout the year. However, it occurs above the thermocline (Fig. 5d) where wci is predominantly horizontal. Temperatures above the 24°C isotherm are more likely to be in the mixed layer rather than in the thermocline (Figs. 6c,f,i).

Seasonally stratified cross-equatorial cross sections of diabatic upwelling according to Eq. (1) along 140°W: (left) total, (center) induced by vertical mixing, and (right) solar penetration. Text insets are the average total diabatic upwelling on the equator at 140°W. Hatching indicates that the isotherm is present less than 50% of the time. Contours in the central column indicate absolute vertical shear. Horizontal dashed lines indicate that 20° and 24°C isotherms.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Seasonally stratified cross-equatorial cross sections of diabatic upwelling according to Eq. (1) along 140°W: (left) total, (center) induced by vertical mixing, and (right) solar penetration. Text insets are the average total diabatic upwelling on the equator at 140°W. Hatching indicates that the isotherm is present less than 50% of the time. Contours in the central column indicate absolute vertical shear. Horizontal dashed lines indicate that 20° and 24°C isotherms.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
Seasonally stratified cross-equatorial cross sections of diabatic upwelling according to Eq. (1) along 140°W: (left) total, (center) induced by vertical mixing, and (right) solar penetration. Text insets are the average total diabatic upwelling on the equator at 140°W. Hatching indicates that the isotherm is present less than 50% of the time. Contours in the central column indicate absolute vertical shear. Horizontal dashed lines indicate that 20° and 24°C isotherms.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Seasonal cycle of (a),(d),(g) total diabatic upwelling (wci); (b),(e),(h) vertical shear; and (c),(f),(i) temperature for three different longitudes in the cold tongue along the equator. The strength of the contribution of solar and vertical mixing induced diabatic wci is indicated by dots (solar higher than 0.3 m day−1) and slashes (vertical mixing higher than 0.3 m day−1). Note that wci and absolute vertical shear are shown in temperature coordinates, while the temperature is shown in depth coordinates.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Seasonal cycle of (a),(d),(g) total diabatic upwelling (wci); (b),(e),(h) vertical shear; and (c),(f),(i) temperature for three different longitudes in the cold tongue along the equator. The strength of the contribution of solar and vertical mixing induced diabatic wci is indicated by dots (solar higher than 0.3 m day−1) and slashes (vertical mixing higher than 0.3 m day−1). Note that wci and absolute vertical shear are shown in temperature coordinates, while the temperature is shown in depth coordinates.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
Seasonal cycle of (a),(d),(g) total diabatic upwelling (wci); (b),(e),(h) vertical shear; and (c),(f),(i) temperature for three different longitudes in the cold tongue along the equator. The strength of the contribution of solar and vertical mixing induced diabatic wci is indicated by dots (solar higher than 0.3 m day−1) and slashes (vertical mixing higher than 0.3 m day−1). Note that wci and absolute vertical shear are shown in temperature coordinates, while the temperature is shown in depth coordinates.
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
To understand the variability of thermocline diabatic upwelling we decompose wci into its physical processes according to Eq. (1) (bottom row, see discussion in section 2b) along 140°W. The signal of wci throughout the seasonal cycle, much like on interannual time scales (Deppenmeier et al. 2021), is dominated by two physical drivers: vertical mixing (Fig. 5, center column) and solar penetration (Fig. 5, right column).
In the mixed layer near the surface (also indicated through its proximity to the hatching in Fig. 5), the warming effect of solar penetration (Figs. 5c,f,i,l) dominates wci. Close to the surface
Below the mixed layer
We find that interannual wci variability is driven by
b. Zonal dependence of seasonal diabatic upwelling
The evolution of diabatic upwelling wci in the eastern Pacific along the equator is similar at 110°, 125°, and 140°W (Fig. 6). The total diabatic upwelling wci (colors, Figs. 6a,d,g) is stronger in the east (110°W, Fig. 6a) and weakens toward the west (140°W, Fig. 6g). Consistent with the increase of SST to the west along the equator (Figs. 6c,f,i), the maximum wci occurs at warmer temperatures in the west than in the east.
At 110°, 125°, and 140°W, two wci maxima exist in the year: one between June and September, and one between February and April (Figs. 6a,d,g). The bulk of the diabatic upwelling is driven by vertical mixing which is modulated by the seasonal cycle of vertical shear (Figs. 6b,e,h). In boreal summer (JJA) the maxima in wci occur at the same time as maxima in vertical shear for all longitudes. Note that it is the vertical derivative of the turbulent heat flux J, which impacts
Between February and April, however, the second maxima of wci at 125° and 140°W does not correspond to a maxima in shear (Figs. 6d,e,g,h). Instead, during these months large wci maxima occur in warmer temperature classes during the time of the large contribution from solar penetration (dots in Fig. 6).
We conclude that there are two regimes of diabatic upwelling over the seasonal cycle on the equator: in February–April, solar penetration contributes to water mass transformation in warm temperature classes close to the surface (solar regime). In all other months, vertical mixing driven by the mean vertical shear of the EUC enables diabatic upwelling. The temperatures in which wci is active in the solar driven regime are warm and above the thermocline. This wci is hence better understood as an intensification of poleward water mass transformation, rather than strong diabatic upwelling in the vertical.
c. Tropical instability waves impact on wci
On subseasonal scales, there is evidence that TIWs assert a rectified effect on the diabatic limb of the upwelling on the equator and to the north. We demonstrate this for an example year (2010), for which we show Hovmöller diagrams of time against longitude for SST together with TIW kinetic energy, wci,

Hovmöller diagrams of (a) sea surface temperature, (b) total diabatic upwelling wci, (c) vertical-mixing-driven diabatic upwelling
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Hovmöller diagrams of (a) sea surface temperature, (b) total diabatic upwelling wci, (c) vertical-mixing-driven diabatic upwelling
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
Hovmöller diagrams of (a) sea surface temperature, (b) total diabatic upwelling wci, (c) vertical-mixing-driven diabatic upwelling
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Hovmöller diagrams of (a) sea surface temperature, (b) total diabatic upwelling wci, (c) vertical-mixing-driven diabatic upwelling
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Hovmöller diagrams of (a) sea surface temperature, (b) total diabatic upwelling wci, (c) vertical-mixing-driven diabatic upwelling
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
Hovmöller diagrams of (a) sea surface temperature, (b) total diabatic upwelling wci, (c) vertical-mixing-driven diabatic upwelling
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
The diabatic upwelling signal stems from both the mixing
d. Diabatic upwelling variability across time scales
This study and the results in Deppenmeier et al. (2021) demonstrate that diabatic upwelling is strong throughout the eastern tropical Pacific and varies on subseasonal to interannual time scales. We see comparable amplitude of wci variability for subseasonal (specifically TIW modulated), seasonal, and interannual (ENSO) time scales (Figs. 9a–c). For the subseasonal (TIW) time scale we separately examine high and low TIW kinetic energy periods using a 200 cm2 s−2 threshold (the results are insensitive within the cutoffs of 150–250 cm2 s−2). For the interannual time scale we define La Niña and El Niño phases via Niño-3.4 SST anomalies that exceed ±0.4°C for at least 6 months (Deppenmeier et al. 2021).
Large diabatic upwelling in the seasonal cycle occurs in two regimes. In March, the shallow maximum near 25°C reflects primarily solar penetration (Figs. 6d and 9a) during this season of weak cloudiness (Klein and Hartmann 1993; Yu and Mechoso 1999; Xu et al. 2005). In all other seasons maximum diabatic upwelling occurs near 22°C below most of the solar influence and just above the thermocline, driven by vertical mixing (Figs. 6e and 9a). Both these seasonal maxima are comparable to or larger than the ENSO or TIW peaks (cf. Figs. 9a–c).

Profiles of total diabatic upwelling wci at 0°, 125°W averaged over different time scales: (a) the climatological seasonal cycle; (b) El Niño, La Niña, and neutral ENSO conditions; and (c) high and low TIW kinetic energy (TIW KE), both for the whole year (dashed) and for SON months only (solid). (d) Maximum wci for the seasonal cycle (letters), high and low TIW KE (triangles), and ENSO phases (dots); (e) maximum shear in the column split by time scales; and (b) mean SST for the different time scales. Symbols in (d)–(f) are color coded as in the legends for (a)–(c).
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1

Profiles of total diabatic upwelling wci at 0°, 125°W averaged over different time scales: (a) the climatological seasonal cycle; (b) El Niño, La Niña, and neutral ENSO conditions; and (c) high and low TIW kinetic energy (TIW KE), both for the whole year (dashed) and for SON months only (solid). (d) Maximum wci for the seasonal cycle (letters), high and low TIW KE (triangles), and ENSO phases (dots); (e) maximum shear in the column split by time scales; and (b) mean SST for the different time scales. Symbols in (d)–(f) are color coded as in the legends for (a)–(c).
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
Profiles of total diabatic upwelling wci at 0°, 125°W averaged over different time scales: (a) the climatological seasonal cycle; (b) El Niño, La Niña, and neutral ENSO conditions; and (c) high and low TIW kinetic energy (TIW KE), both for the whole year (dashed) and for SON months only (solid). (d) Maximum wci for the seasonal cycle (letters), high and low TIW KE (triangles), and ENSO phases (dots); (e) maximum shear in the column split by time scales; and (b) mean SST for the different time scales. Symbols in (d)–(f) are color coded as in the legends for (a)–(c).
Citation: Journal of Physical Oceanography 52, 11; 10.1175/JPO-D-21-0316.1
As the thermocline and EUC begin to deepen in early boreal summer, strong shear squared S2 becomes large at shallow depths and vertical mixing dominates in the otherwise solar-influenced upper region of the water column. The vertical mixing
We show that increased S2 also contributes to increased vertical-mixing-driven
During La Niña conditions wci has values as large as 0.6 m day−1, while during El Niño conditions wci maxima are only about half as big and occur at warmer temperatures (Fig. 9b). It is noteworthy that the time scales we investigate are connected. TIW occurrence varies both seasonally and interannually, being practically absent in March–May and during El Niño, and enhanced during September–February and during La Niña. Similarly, ENSO evolution is also connected to the seasonal cycle. Whether TIWs dominantly influence the variability we see on other time scales is a matter of active research.
To place the wci variability in context, we compare it to SST variability across the time scales (Fig. 9f). SST variability is largest (by construction) between ENSO phases (the Niño-3.4 index exceeds ±0.4 for several months). This is followed by the strength of the variability of the seasonal cycle, which displays the largest variability in S2 and wci. TIW phases also strongly influence SST, with low TIW kinetic energy usually being accompanied by warmer temperatures and high TIW kinetic energy by colder temperatures. We do not find a one to one relationship between the strength of wci variability and SST variability for these difference processes.
4. Conclusions
Diabatic upwelling wci varies across time scales. Most noticeably, large variations in wci are seen in the seasonal cycle. In all seasons but March–May (MAM) the maximum in wci is driven predominantly by vertical mixing induced by the strong vertical shear above the equatorial undercurrent. In MAM, wci is confined to warm temperature classes, and solar penetration dominates. In these warm temperature classes positive wci consists of poleward water mass transformation in the mixed layer—equivalent to solar heating shifting warmer isotherms toward the equator—rather than vertical motion of
While most of the seasonal cycle variability can be explained by changes in the seasonal mean quantities, such as differences in S2, the results here suggest that TIWs can rectify to induce seasonal variations in wci. Part of this signal contained in the vertical mixing component is induced by eddy stirring, and part of it appears as the covariance term in Eq. (1) from covariability of velocity and temperature during the passage of a TIW. Strengthened turbulent heat fluxes during the passage of TIWs have also been described by Lien et al. (2008), Moum et al. (2009), Inoue et al. (2012), and Holmes and Thomas (2015) on the equator and Cherian et al. (2021) off the equator. Here, we find evidence that passing TIW rectify to create diabatic upwelling that results in increased water mass transformation on the equator.
Observations of vertical mixing are clearly important to understand all components of the tropical Pacific circulation and heat budget. Here we demonstrate that the impact of vertical mixing variability on diabatic upwelling is not dominated by a specific time scale, rather there is sizeable variability at each of the subseasonal, seasonal, and interannual time scales. Our estimate of diabatic upwelling is based on a model simulation which captures the main aspects of the circulation reasonably well, but also displays some biases such as a too deep thermocline and weak TIW kinetic energy. This highlights the importance of observationally constraining the diabatic component of the tropical Pacific circulation and underscores the need for targeted observational efforts.
Acknowledgments.
This study was supported by the National Center for Atmospheric Research, which is a major facility sponsored by the National Science Foundation under Cooperative Agreement 1852977 and by the NOAA Climate Variability Program under Grant Agreement NA18OAR4310399. LT was supported by NOAA Climate Variability Program under Grant Agreement NA18OAR4310400. This work is part of the CVP TPOS 2020 pre-field modeling effort led by Dr. Sandy Lucas and is PMEL Contribution Number 5336. Computing resources were provided by the Climate Simulation Laboratory at NCAR’s Computational and Information Systems Laboratory (CISL). TAO observational data were made available by the GTMBA Project Office of NOAA/PMEL and accessed through the website https://www.pmel.noaa.gov/tao/drupal/disdel/. The analysis benefited from the open source packages dask (Rocklin 2015), xarray (Hoyer and Hamman 2017), and xgcm (Abernathy et al. 2020).
Data availability statement.
The data underlying the figures in this manuscript are available via https://zenodo.org/record/6762930. TAO observational data can be accessed through https://www.pmel.noaa.gov/tao/drupal/disdel/.
REFERENCES
Abernathy, R., and Coauthors, 2020: xgcm: v0.6.0. Zenodo, https://doi.org/10.5281/zenodo.5752560.
Alexander, M. A., I. Bladé, M. Newman, J. R. Lanzante, N.-C. Lau, and J. D. Scott, 2002: The atmospheric bridge: The influence of ENSO teleconnections on air–sea interaction over the global oceans. J. Climate, 15, 2205–2231, https://doi.org/10.1175/1520-0442(2002)015<2205:TABTIO>2.0.CO;2.
Bosc, C., and T. Delcroix, 2008: Observed equatorial Rossby waves and ENSO-related warm water volume changes in the equatorial Pacific Ocean. J. Geophys. Res., 113, C06003, https://doi.org/10.1029/2007JC004613.
Brown, J., and A. Federov, 2010: Estimating the diapycnal transport contribution to warm water volume variations in the tropical Pacific Ocean. J. Climate, 23, 221–237, https://doi.org/10.1175/2009JCLI2347.1.
Bryan, F., and S. Bachman, 2015: Isohaline salinity budget of the North Atlantic salinity maximum. J. Phys. Oceanogr., 45, 724–736, https://doi.org/10.1175/JPO-D-14-0172.1.
Bryden, H. L., and E. C. Brady, 1985: Diagnostic model of the three-dimensional circulation in the upper equatorial Pacific Ocean. J. Phys. Oceanogr., 15, 1255–1273, https://doi.org/10.1175/1520-0485(1985)015<1255:DMOTTD>2.0.CO;2.
Cai, W., and Coauthors, 2019: Pantropical climate interactions. Science, 363, eaav4236, https://doi.org/10.1126/science.aav4236.
Chelton, D. B., F. J. Wentz, C. L. Gentemann, R. A. de Szoeke, and M. G. Schlax, 2000: Satellite microwave SST observations of transequatorial tropical instability waves. Geophys. Res. Lett., 27, 1239–1242, https://doi.org/10.1029/1999GL011047.
Cherian, D., D. Whitt, R. Holmes, R.-C. Lien, S. Bachman, and W. Large, 2021: Off-equatorial deep-cycle turbulence forced by tropical instability waves in the equatorial Pacific. J. Phys. Oceanogr., 51, 1575–1593, https://doi.org/10.1175/JPO-D-20-0229.1.
Claar, D. C., L. Szostek, J. M. McDevitt-Irwin, J. J. Schanze, and J. K. Baum, 2018: Global patterns and impacts of El Niño events on coral reefs: A meta-analysis. PLOS ONE, 13, e0190957, https://doi.org/10.1371/journal.pone.0190957.
Clarke, A. J., S. Van Gorder, and G. Colantuono, 2007: Wind stress curl and ENSO discharge/recharge in the equatorial Pacific. J. Phys. Oceanogr., 37, 1077–1091, https://doi.org/10.1175/JPO3035.1.
Deppenmeier, A.-L., F. O. Bryan, W. S. Kessler, and L. Thompson, 2021: Modulation of cross-isothermal velocities with ENSO in the tropical Pacific cold tongue. J. Phys. Oceanogr., 51, 1559–1574, https://doi.org/10.1175/JPO-D-20-0217.1.
Fiedler, P. C., and L. D. Talley, 2006: Hydrography of the eastern tropical Pacific: A review. Prog. Oceanogr., 69, 143–180, https://doi.org/10.1016/j.pocean.2006.03.008.
Groeskamp, S., S. M. Griffies, D. Iudicone, R. Marsh, A. G. Nurser, and J. D. Zika, 2019: The water mass transformation framework for ocean physics and biogeochemistry. Annu. Rev. Mar. Sci., 11, 271–305, https://doi.org/10.1146/annurev-marine-010318-095421.
Halpern, D., R. A. Knox, and D. S. Luther, 1988: Observations of 20-day period meridional current oscillations in the upper ocean along the Pacific equator. J. Phys. Oceanogr., 18, 1514–1534, https://doi.org/10.1175/1520-0485(1988)018<1514:OODPMC>2.0.CO;2.
Hayes, S., L. Mangum, J. Picaut, A. Sumi, and K. Takeuchi, 1991: TOGA-TAO: A moored array for real-time measurements in the tropical Pacific Ocean. Bull. Amer. Meteor. Soc., 72, 339–347, https://doi.org/10.1175/1520-0477(1991)072<0339:TTAMAF>2.0.CO;2.
Hieronymus, M., J. Nilsson, and J. Nycander, 2014: Water mass transformation in salinity–temperature space. J. Phys. Oceanogr., 44, 2547–2568, https://doi.org/10.1175/JPO-D-13-0257.1.
Holmes, R. M., and L. N. Thomas, 2015: The modulation of equatorial turbulence by tropical instability waves in a regional ocean model. J. Phys. Oceanogr., 45, 1155–1173, https://doi.org/10.1175/JPO-D-14-0209.1.
Hoyer, S., and J. Hamman, 2017: xarray: N-D labeled arrays and datasets in Python. J. Open Res. Software, 5, 10, https://doi.org/10.5334/jors.148.
Huguenin, M. F., R. M. Holmes, and M. H. England, 2020: Key role of diabatic processes in regulating warm water volume variability over ENSO events. J. Climate, 33, 9945–9964, https://doi.org/10.1175/JCLI-D-20-0198.1.
Inoue, R., R.-C. Lien, and J. Moum, 2012: Modulation of equatorial turbulence by a tropical instability wave. J. Geophys. Res., 117, C10009, https://doi.org/10.1029/2011JC007767.
Inoue, R., R.-C. Lien, J. N. Moum, R. C. Perez, and M. C. Gregg, 2019: Variations of equatorial shear, stratification, and turbulence within a tropical instability wave cycle. J. Geophys. Res. Oceans, 124, 1858–1875, https://doi.org/10.1029/2018JC014480.
Johnson, G. C., M. J. McPhaden, and E. Firing, 2001: Equatorial Pacific Ocean horizontal velocity, divergence, and upwelling. J. Phys. Oceanogr., 31, 839–849, https://doi.org/10.1175/1520-0485(2001)031<0839:EPOHVD>2.0.CO;2.
Jong, B.-T., M. Ting, and R. Seager, 2016: El Niño’s impact on California precipitation: Seasonality, regionality, and El Niño intensity. Environ. Res. Lett., 11, 054021, https://doi.org/10.1088/1748-9326/11/5/054021.
Kessler, W. S., 2006: The circulation of the eastern tropical Pacific: A review. Prog. Oceanogr., 69, 181–217, https://doi.org/10.1016/j.pocean.2006.03.009.
Klein, S. A., and D. L. Hartmann, 1993: The seasonal cycle of low stratiform clouds. J. Climate, 6, 1587–1606, https://doi.org/10.1175/1520-0442(1993)006<1587:TSCOLS>2.0.CO;2.
Kobayashi, S., and Coauthors, 2015: The JRA-55 reanalysis: General specifications and basic characteristics. J. Meteor. Soc. Japan, 93, 5–48, https://doi.org/10.2151/jmsj.2015-001.
Legeckis, R., 1977: Long waves in the eastern equatorial Pacific Ocean: A view from a geostationary satellite. Science, 197, 1179–1181, https://doi.org/10.1126/science.197.4309.1179.
Lengaigne, M., U. Hausmann, G. Madec, C. Menkes, J. Vialard, and J. Molines, 2012: Mechanisms controlling warm water volume interannual variations in the equatorial Pacific: Diabatic versus adiabatic processes. Climate Dyn., 38, 1031–1046, https://doi.org/10.1007/s00382-011-1051-z.
Lien, R.-C., E. A. d’Asaro, and C. E. Menkes, 2008: Modulation of equatorial turbulence by tropical instability waves. Geophys. Res. Lett., 35, L24607, https://doi.org/10.1029/2008GL035860.
McPhaden, M. J., and Coauthors, 1998: The tropical ocean-global atmosphere observing system: A decade of progress. J. Geophys. Res., 103, 14 169–14 240, https://doi.org/10.1029/97JC02906.
McPhaden, M. J., A. J. Busalacchi, and D. L. Anderson, 2010: A toga retrospective. Oceanography, 23, 86–103, https://doi.org/10.5670/oceanog.2010.26.
Meinen, C. S., and M. J. McPhaden, 2000: Observations of warm water volume changes in the equatorial Pacific and their relationship to El Niño and La Niña. J. Climate, 13, 3551–3559, https://doi.org/10.1175/1520-0442(2000)013<3551:OOWWVC>2.0.CO;2.
Meinen, C. S., and M. J. McPhaden, 2001: Interannual variability in warm water volume transports in the equatorial Pacific during 1993–99. J. Phys. Oceanogr., 31, 1324–1345, https://doi.org/10.1175/1520-0485(2001)031<1324:IVIWWV>2.0.CO;2.
Meinen, C., M. J. McPhaden, and G. C. Johnson, 2001: Vertical velocities and transports in the equatorial Pacific during 1993–99. J. Phys. Oceanogr., 31, 3230–3248, https://doi.org/10.1175/1520-0485(2001)031<3230:VVATIT>2.0.CO;2.
Miller, L., D. R. Watts, and M. Wimbush, 1985: Oscillations of dynamic topography in the eastern equatorial Pacific. J. Phys. Oceanogr., 15, 1759–1770, https://doi.org/10.1175/1520-0485(1985)015<1759:OODTIT>2.0.CO;2.
Moum, J. N., R.-C. Lien, A. Perlin, J. Nash, M. Gregg, and P. Wiles, 2009: Sea surface cooling at the equator by subsurface mixing in tropical instability waves. Nat. Geosci., 2, 761–765, https://doi.org/10.1038/ngeo657.
Moum, J. N., A. Perlin, J. D. Nash, and M. J. McPhaden, 2013: Seasonal sea surface cooling in the equatorial Pacific cold tongue controlled by ocean mixing. Nature, 500, 64–67, https://doi.org/10.1038/nature12363.
Niiler, P., and J. Stevenson, 1982: The heat budget of tropical ocean warm-water pools. J. Mar. Res., 40, 465–480.
Nurser, A., R. Marsh, and R. G. Williams, 1999: Diagnosing water mass formation from air–sea fluxes and surface mixing. J. Phys. Oceanogr., 29, 1468–1487, https://doi.org/10.1175/1520-0485(1999)029<1468:DWMFFA>2.0.CO;2.
Ohlmann, J. C., 2003: Ocean radiant heating in climate models. J. Climate, 16, 1337–1351, https://doi.org/10.1175/1520-0442-16.9.1337.
Perez, R. C., M. F. Cronin, and W. S. Kessler, 2010: Tropical cells and a secondary circulation near the northern front of the equatorial Pacific cold tongue. J. Phys. Oceanogr., 40, 2091–2106, https://doi.org/10.1175/2010JPO4366.1.
Qiao, L., and R. H. Weisberg, 1995: Tropical instability wave kinematics: Observations from the tropical instability wave experiment. J. Geophys. Res., 100, 8677–8693, https://doi.org/10.1029/95JC00305.
Rai, S., M. Hecht, M. Maltrud, and H. Aluie, 2021: Scale of oceanic eddy killing by wind from global satellite observations. Sci. Adv., 7, eabf4920, https://doi.org/10.1126/sciadv.abf4920.
Renault, L., S. Masson, T. Arsouze, G. Madec, and J. C. McWilliams, 2020: Recipes for how to force oceanic model dynamics. J. Adv. Model. Earth Syst., 12, e2019MS001715, https://doi.org/10.1029/2019MS001715.
Rocklin, M., 2015: Dask: Parallel computation with blocked algorithms and task scheduling. Proc. 14th Python in Science Conf., Austin, TX, SciPy, 126–132, https://doi.org/10.25080/Majora-7b98e3ed-013.
Smith, N., and Coauthors, 2019: Tropical Pacific observing system. Front. Mar. Sci., 6, 31, https://doi.org/10.3389/fmars.2019.00031.
Smith, R., and Coauthors, 2010: The Parallel Ocean Program (POP) reference manual: Ocean component of the Community Climate System Model (CCSM). Los Alamos Tech. Rep. LAUR-10-01853, 140 pp.
Smyth, W., and J. Moum, 2013: Marginal instability and deep cycle turbulence in the eastern equatorial Pacific Ocean. Geophys. Res. Lett., 40, 6181–6185, https://doi.org/10.1002/2013GL058403.
Smyth, W., J. Moum, L. Li, and S. Thorpe, 2013: Diurnal shear instability, the descent of the surface shear layer, and the deep cycle of equatorial turbulence. J. Phys. Oceanogr., 43, 2432–2455, https://doi.org/10.1175/JPO-D-13-089.1.
Sun, C., W. D. Smyth, and J. N. Moum, 1998: Dynamic instability of stratified shear flow in the upper equatorial Pacific. J. Geophys. Res., 103, 10 323–10 337, https://doi.org/10.1029/98JC00191.
Sun, Z., H. Liu, P. Lin, Y. Tseng, J. Small, and F. Bryan, 2019: The modeling of the north equatorial countercurrent in the community Earth system model and its oceanic component. J. Adv. Model. Earth Syst., 11, 531–544, https://doi.org/10.1029/2018MS001521.
Toole, J. M., H.-M. Zhang, and M. J. Caruso, 2004: Time-dependent internal energy budgets of the tropical warm water pools. J. Climate, 17, 1398–1410, https://doi.org/10.1175/1520-0442(2004)017<1398:TIEBOT>2.0.CO;2.
Tsujino, H., and Coauthors, 2018: JRA-55 based surface dataset for driving ocean–sea-ice models (JRA55-do). Ocean Modell., 130, 79–139, https://doi.org/10.1016/j.ocemod.2018.07.002.
Walin, G., 1982: On the relation between sea-surface heat flow and thermal circulation in the ocean. Tellus, 34, 187–195, https://doi.org/10.3402/tellusa.v34i2.10801.
Wang, W., and M. J. McPhaden, 1999: The surface-layer heat balance in the equatorial Pacific Ocean. Part I: Mean seasonal cycle. J. Phys. Oceanogr., 29, 1812–1831, https://doi.org/10.1175/1520-0485(1999)029<1812:TSLHBI>2.0.CO;2.
Warner, S. J., and J. N. Moum, 2019: Feedback of mixing to ENSO phase change. Geophys. Res. Lett., 46, 13 920–13 927, https://doi.org/10.1029/2019GL085415.
Weisberg, R. H., and L. Qiao, 2000: Equatorial upwelling in the central Pacific estimated from moored velocity profilers. J. Phys. Oceanogr., 30, 105–124, https://doi.org/10.1175/1520-0485(2000)030<0105:EUITCP>2.0.CO;2.
Willett, C. S., R. R. Leben, and M. F. Lavín, 2006: Eddies and tropical instability waves in the eastern tropical Pacific: A review. Prog. Oceanogr., 69, 218–238, https://doi.org/10.1016/j.pocean.2006.03.010.
Wyrtki, K., 1981: An estimate of equatorial upwelling in the Pacific. J. Phys. Oceanogr., 11, 1205–1214, https://doi.org/10.1175/1520-0485(1981)011<1205:AEOEUI>2.0.CO;2.
Xu, H., S.-P. Xie, and Y. Wang, 2005: Subseasonal variability of the southeast Pacific stratus cloud deck. J. Climate, 18, 131–142, https://doi.org/10.1175/JCLI3250.1.
Yeh, S.-W., and Coauthors, 2018: ENSO atmospheric teleconnections and their response to greenhouse gas forcing. Rev. Geophys., 56, 185–206, https://doi.org/10.1002/2017RG000568.
Yu, J.-Y., and C. R. Mechoso, 1999: Links between annual variations of Peruvian stratocumulus clouds and of SST in the eastern equatorial Pacific. J. Climate, 12, 3305–3318, https://doi.org/10.1175/1520-0442(1999)012<3305:LBAVOP>2.0.CO;2.