Abstract

Physical processes that control ENSO are relatively fast. For instance, it takes only several months for a Kelvin wave to cross the Pacific basin (T_k ≈ 2 months), while Rossby waves travel the same distance in about half a year. Compared to such short time scales, the typical periodicity of El Niño is much longer (T ≈ 2–7 yr). Thus, ENSO is fundamentally a low-frequency phenomenon in the context of these faster processes. Here, the author takes advantage of this fact and uses the smallness of the ratio ε _k = T_k /T to expand solutions of the ocean shallow-water equations into power series (the actual parameter of expansion also includes the oceanic damping rate). Using such an expansion, referred to here as the low-frequency approximation, the author relates thermocline depth anomalies to temperature variations in the eastern equatorial Pacific via an explicit integral operator. This allows a simplified formulation of ENSO dynamics based on an integro-differential equation. Within this formulation, the author shows how the interplay between wind stress curl and oceanic damping rates affects 1) the amplitude and periodicity of El Niño and 2) the phase lag between variations in the equatorial warm water volume and SST in the eastern Pacific. A simple analytical expression is derived for the phase lag. Further, applying the low-frequency approximation to the observed variations in SST, the author computes thermocline depth anomalies in the western and eastern equatorial Pacific to show a good agreement with the observed variations in warm water volume. Ultimately, this approach provides a rigorous framework for deriving other simple models of ENSO (the delayed and recharge oscillators), highlights the limitations of such models, and can be easily used for decadal climate variability in the Pacific.

Abstract

Changes in background zonal wind in the tropical Pacific are often invoked to explain changes in ENSO properties. However, the sensitivity of ENSO to mean zonal winds has been thoroughly explored only in intermediate coupled models (following Zebiak and Cane), not in coupled GCMs. The role of mean meridional winds has received even less attention. Accordingly, the goal of this study is to examine systematically the effects of both zonal (equatorial) and meridional (cross-equatorial) background winds on ENSO using targeted experiments with a comprehensive climate model (CESM). Changes in the mean winds are generated by imposing heat flux forcing in two confined regions at a sufficient distance north and south of the equator. We find that the strengthening of either wind component reduces ENSO amplitude, especially eastern Pacific SST variability, and inhibits meridional swings of the intertropical convergence zone (ITCZ). The effect of zonal winds is generally stronger than that of meridional winds. A stability analysis reveals that the strengthening of zonal and meridional winds weakens the ENSO key positive feedbacks, specifically the zonal advection and thermocline feedbacks, which explains these changes. Zonal wind enhancement also intensifies mean upwelling and hence dynamical damping, leading to a further weakening of El Niño events. Ultimately, this study argues that the zonal and, to a lesser extent, meridional wind strengthening of the past decades may have contributed to the observed shift of El Niño characteristics after the year 2000.

Abstract

Adiabatic theories of ocean circulation and density structure have a long tradition, from the concept of the ventilated thermocline to the notion that deep ocean ventilation is controlled by westerly winds over the Southern Ocean. This study explores these ideas using a recently developed Lagrangian ocean model (LOM), which simulates ocean motions by computing trajectories of water parcels. A unique feature of the LOM is its capacity to model ocean circulations in the adiabatic limit, in which water parcels exactly conserve their densities when they are not in contact with the ocean surface. The authors take advantage of this property of the LOM and consider the circulation and stratification that develop in an ocean with a fully adiabatic interior (with both isopycnal and diapycnal diffusivities set to zero). The ocean basin in the study mimics that of the Atlantic Ocean and includes a circumpolar channel. The model is forced by zonal wind stress and a density restoring at the surface.

Despite the idealized geometry, the relatively coarse model resolution, and the lack of atmospheric coupling, the nondiffusive ocean maintains a density structure and meridional overturning that are broadly in line with those observed in the Atlantic Ocean. These are generated by just a handful of key water pathways, including shallow tropical cells described by ventilated thermocline theory; a deep overturning cell in which sinking North Atlantic Deep Water eventually upwells in the Southern Ocean before returning northward as Antarctic Intermediate Water; a Deacon cell that results from a topographically steered and corkscrewing circumpolar current; and weakly overturning Antarctic Bottom Water, which is effectively ventilated only in the Southern Hemisphere.

The main conclusion of this study is that the adiabatic limit for the ocean interior provides the leading-order solution for ocean overturning and density structure, with tracer diffusion contributing first-order perturbations. Comparing nondiffusive and diffusive experiments helps to quantify the changes in stratification and circulation that result from adding a moderate amount of tracer diffusion in the ocean model, and these include an increase in the amplitude of the deep meridional overturning cell of several Sverdrups, a 10%–20% increase in Northern Hemispheric northward heat transport, a stronger stratification just below the main thermocline, and a more realistic bottom overturning cell.

Abstract

The dynamics of El Niño–Southern Oscillation (ENSO) are studied in terms of the balance between energy input from the winds (via wind power) and changes in the storage of available potential energy in the tropical ocean. Presently, there are broad differences in the way global general circulation models simulate the dynamics, magnitude, and phase of ENSO events; hence, there is a need for simple, physically based metrics to allow for model evaluation. This energy description is a basinwide, integral, quantitative approach, ideal for intermodel comparison, that assesses model behavior in the subsurface ocean. Here it is applied to a range of ocean models and data assimilations within ENSO spatial and temporal scales. The onset of an El Niño is characterized by a decrease in wind power that leads to a decrease in available potential energy, and hence a flatter thermocline. In contrast, La Niña events are preceded by an increase in wind power that leads to an increase in the available potential energy and a steeper thermocline. The wind power alters the available potential energy via buoyancy power, associated with vertical mass fluxes that modify the slope of the isopycnals. Only a fraction of wind power is converted to buoyancy power. The efficiency of this conversion γ is estimated in this study at 50%–60%. Once the energy is delivered to the thermocline it is subject to small, but important, diffusive dissipation. It is estimated that this dissipation sets the e-folding damping rate α for the available potential energy on the order of 1 yr⁻¹. The authors propose to use the efficiency γ and the damping rate α as two energy-based metrics for evaluating dissipative properties of the ocean component of general circulation models, providing a simple method for understanding subsurface ENSO dynamics and a diagnostic tool for exploring differences between the models.

Abstract

Variations in the warm water volume (WWV) of the equatorial Pacific Ocean are considered a key element of the dynamics of the El Niño–Southern Oscillation (ENSO) phenomenon. WWV, a proxy for the upper-ocean heat content, is usually defined as the volume of water with temperatures greater than 20°C. It has been suggested that the observed variations in WWV are controlled by interplay among meridional, zonal, and vertical transports (with vertical transports typically calculated as the residual of temporal changes in WWV and the horizontal transport divergence). Here, the output from a high-resolution ocean model is used to calculate the zonal and meridional transports and conduct a comprehensive analysis of the mass balance above the 25 kg m⁻³ σ_θ surface (approximating the 20°C isotherm). In contrast to some earlier studies, the authors found that on ENSO time scales variations in the diapycnal transport across the 25 kg m⁻³ isopycnal are small in the eastern Pacific and negligible in the western and central Pacific. In previous observational studies, the horizontal transports were estimated using Ekman and geostrophic dynamics; errors in these approximations were unavoidably folded into the estimates of the diapycnal transport. Here, the accuracy of such estimates is assessed by recalculating mass budgets using the model output at a spatial resolution similar to that of the observations. The authors show that errors in lateral transports can be of the same order of magnitude as the diapycnal transport itself. Further, the rate of change of WWV correlates well with wind stress curl (a driver of meridional transport). This relationship is explored using an extended version of the Sverdrup balance, and it is shown that the two are correlated because they both have the ENSO signal and not because changes in WWV are solely attributable to the wind stress curl.

Abstract

Interactions between the tropical oceans and atmosphere permit a spectrum of natural modes of oscillation whose properties—period, intensity, spatial structure, and direction of propagation—depend on the background climatic state (i.e., the mean state). This mean state can be described by parameters that include the following: the time-averaged intensity τ of the Pacific trade winds, the mean depth (H) of the thermocline, and the temperature difference across the thermocline (ΔT). A stability analysis by means of a simple coupled ocean–atmosphere model indicates two distinct families of unstable modes. One has long periods of several years, involves sea surface temperature variations determined by vertical movements of the thermocline that are part of the adjustment of the ocean basin to the fluctuating winds, requires a relatively deep thermocline, and corresponds to the delayed oscillator. The other family requires a shallow thermocline, has short periods of a year or two, has sea surface temperature variations determined by advection and by entrainment across the thermocline, and is associated with westward phase propagation. For the modes to be unstable, both families require that the background zonal wind exceed a certain intensity. An increase in ΔT, and in H beyond a certain value, are stabilizing. For intermediate values of H, between large values that favor the one mode and small values that favor the other, the modes are of a hybrid type with some properties of each family. The observed Southern Oscillation has been of this type for the past few decades, but some paleorecords suggest that, in the distant past, the oscillation was strictly of the delayed oscillator type and had a very long period on the order of a decade.

Abstract

Westerly wind bursts (WWBs)—brief but strong westerly wind anomalies in the equatorial Pacific—are believed to play an important role in El Niño–Southern Oscillation (ENSO) dynamics, but quantifying their effects is challenging. Here, we investigate the cumulative effects of WWBs on ENSO characteristics, including the occurrence of extreme El Niño events, via modified coupled model experiments within Community Earth System Model (CESM1) in which we progressively reduce the impacts of wind stress anomalies associated with model-generated WWBs. In these “wind stress shaving” experiments we limit momentum transfer from the atmosphere to the ocean above a preset threshold, thus “shaving off” wind bursts. To reduce the tropical Pacific mean state drift, both westerly and easterly wind bursts are removed, although the changes are dominated by WWB reduction. As we impose progressively stronger thresholds, both ENSO amplitude and the frequency of extreme El Niño decrease, and ENSO becomes less asymmetric. The warming center of El Niño shifts westward, indicating less frequent and weaker eastern Pacific (EP) El Niño events. Removing most wind burst–related wind stress anomalies reduces ENSO mean amplitude by 22%. The essential role of WWBs in the development of extreme El Niño events is highlighted by the suppressed eastward migration of the western Pacific warm pool and hence a weaker Bjerknes feedback under wind shaving. Overall, our results reaffirm the importance of WWBs in shaping the characteristics of ENSO and its extreme events and imply that WWB changes with global warming could influence future ENSO.

Abstract

This study compares the impacts of Arctic sea ice decline on the Atlantic meridional overturning circulation (AMOC) in two configurations of the Community Earth System Model with different horizontal resolution. In a suite of model experiments, we impose radiative imbalance at the ice surface, replicating a loss of sea ice cover comparable to that observed during 1979–2014, and we find dramatic differences in the AMOC response between the two models. In the lower-resolution configuration, the AMOC weakens by about one-third over the first 100 years, approaching a new quasi-equilibrium. By contrast, in the higher-resolution configuration, the AMOC weakens by ~10% during the first 20–30 years followed by a full recovery driven by invigorated deep water formation in the Labrador Sea and adjacent regions. We investigate these differences using a diagnostic AMOC stability indicator, which reflects the AMOC freshwater transport in and out of the basin and hence the strength of the basin-scale salt-advection feedback. This indicator suggests that the AMOC in the lower-resolution model is less stable and more sensitive to surface perturbations, as confirmed by hosing experiments mimicking Arctic freshening due to sea ice decline. Differences between the models’ mean states, including the Atlantic Ocean mean surface freshwater fluxes, control the differences in AMOC stability. Our results demonstrate that the AMOC stability indicator is indeed useful for evaluating AMOC sensitivity to perturbations. We emphasize that, despite the differences in the long-term adjustment, both models simulate a multidecadal AMOC weakening caused by Arctic sea ice decline, relevant to climate change.

Abstract

Global climate models frequently exhibit cold biases in tropical sea surface temperature (SST) in the central and eastern equatorial Pacific. Here, Lagrangian particle back trajectories are used to investigate the source regions of the water that upwells along the equator in the IPSL climate model to test and confirm the hypothesis that the SST biases are caused by remote biases advected in from the extratropics and to identify the dominant source regions. Water in the model is found to be sourced primarily from localized regions along the western and eastern flanks of the subtropical gyres. However, while the model SST bias is especially large in the northwestern subtropical Pacific (about −5°C), it is found that the eastern subtropics contribute to the equatorial bias the most. This is due to two distinct subsurface pathways connecting these regions to the equator. The first pathway, originating in the northwestern subtropical Pacific, has relatively long advection time scales close to or exceeding 60 yr, wherein particles recirculate around the subtropical gyres while descending to approximately 500 m before then shoaling toward the equatorial undercurrent. The second pathway, from the eastern subtropics, has time scales close to 10 yr, with particles following a shallow and more direct route to the equator within the upper 200 m. The deeper and longer pathway taken by the western subtropical water ensures that vertical mixing can erode the bias. Ultimately, it is estimated that relatively confined regions in the eastern subtropics of both hemispheres control approximately half of the equatorial bias.

Abstract

This study investigates the excitation of decadal variability and predictability of the ocean climate state in the North Atlantic. Specifically, initial linear optimal perturbations (LOPs) in temperature and salinity that vary with depth, longitude, and latitude are computed, and the maximum impact on the ocean of these perturbations is evaluated in a realistic ocean general circulation model. The computations of the LOPs involve a maximization procedure based on Lagrange multipliers in a nonautonomous context. To assess the impact of these perturbations four different measures of the North Atlantic Ocean state are used: meridional volume and heat transports (MVT and MHT) and spatially averaged sea surface temperature (SST) and ocean heat content (OHC). It is shown that these metrics are dramatically different with regard to predictability. Whereas OHC and SST can be efficiently modified only by basin-scale anomalies, MVT and MHT are also strongly affected by smaller-scale perturbations. This suggests that instantaneous or even annual-mean values of MVT and MHT are less predictable than SST and OHC. Only when averaged over several decades do the former two metrics have predictability comparable to the latter two, which highlights the need for long-term observations of the Atlantic meridional overturning circulation in order to accumulate climatically relevant data. This study also suggests that initial errors in ocean temperature of a few millikelvins, encompassing both the upper and deep ocean, can lead to ~0.1-K errors in the predictions of North Atlantic sea surface temperature on interannual time scales. This transient error growth peaks for SST and OHC after about 6 and 10 years, respectively, implying a potential predictability barrier.