<figure><img src="image_1.png" /></figure>

HAL
open science

# Northeast Pacific warm blobs sustained via extratropical atmospheric teleconnections

Jian Shi, Hao Huang, Alexey Fedorov, Neil J. Holbrook, Yu Zhang, Ruiqiang

Ding, Yongyue Luo, Shengpeng Wang, Jiajie Chen, Xi Hu, et al.

► To cite this version:

Jian Shi, Hao Huang, Alexey Fedorov, Neil J. Holbrook, Yu Zhang, et al. Northeast Pacific warm blobs sustained via extratropical atmospheric teleconnections. Nature Communications, 2024, 15, 10.1038/s41467-024-47032-x. insu-04725893

HAL Id: insu-04725893

https://insu.hal.science/insu-04725893v1

Submitted on 8 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Northeast Pacific warm blobs sustained via extratropical atmospheric teleconnections

Received: 6 April 2023

Accepted: 18 March 2024

Published online: 02 April 2024

<figure><img src="image_4.png" /></figure>

Check for updates

Jian Shi Ⓓ1,2,3, Hao Huang<fnref n="1" />,², Alexey V. Fedorov Ⓓ4,5, Neil J. Holbrook Ⓓ6,7, Yu Zhang Ⓓ1,3✉, Ruiqiang Ding Ⓓ8, Yongyue Luo Ⓓ1,2, Shengpeng Wang³, Jiajie Chen Ⓓ1,2, Xi Hu<fnref n="1" />,², Qinyu Liu<fnref n="1" />, Fei Huang Ⓓ1,2,3✉ & Xiaopei Lin Ⓓ1,2,3

Large-scale marine heatwaves in the Northeast Pacific (NEP), identified here and previously as ‘warm blobs’, have devastating impacts on regional ecosystems. An anomalous atmospheric ridge over the NEP is known to be crucial for maintaining these warm blobs, also causing abnormally cold temperatures over North America during the cold season. Previous studies linked this ridge to teleconnections from tropical sea surface temperature anomalies. However, it was unclear whether teleconnections from the extratropics could also contribute to the ridge. Here we show that planetary wave trains, triggered by increased rainfall and latent heat release over the Mediterranean Sea accompanied by decreased rainfall over the North Atlantic, can transport wave energy to the NEP, guided by the westerly jet, and induce a quasi-barotropic ridge there. Our findings provide insights into extratropical teleconnections sustaining the NEP ridge, offering a source of potential predictability for the warm blobs and temperature fluctuations over North America.

Large-scale marine heatwaves in the NEP region, identified here and previously as ‘warm blobs’, have attracted wide attention since the long-lived, so-called ‘Blob’ event of 2013–2016¹⁻³ and another event in 2019–2020⁴⁻⁶. Since the probability of such marine heatwaves is expected to increase with global warming⁷,⁸, it is critical to understand their driving mechanisms to better understand warm blob potential predictability⁹. Both atmospheric (e.g., anomalous atmospheric pressure ridge and clouds) and oceanic (e.g., salinity-related mixed layer depth variation) processes have been suggested to play important roles in the warm blob formation<fnref n="1" />,²,¹⁰⁻¹².

Among these key processes, a strong and persistent higher-than-normal sea level pressure (SLP)<fnref n="1" /> over the NEP is thought to be critical by modulating surface latent heat flux as well as wind-driven ocean currents and mixing in winter<fnref n="1" />,¹⁰. A number of studies suggest that atmospheric teleconnections originating from the tropics can trigger

this anomalous pressure ridge over the NEP²,¹³⁻¹⁸. For example, Wang et al.¹⁶ proposed a cross-Pacific pathway of Rossby wave energy at 850 hPa during winter 2013/14 in response to anomalous convection over the Philippine Sea, which contributed to the anomalous ridge over the Gulf of Alaska. Shi et al.¹⁸ discussed potential teleconnections to the NEP triggered by different types of El Niño events. However, unlike these important processes originating in the tropics, possible contributions from midlatitudes to the anomalous ridge have received much less attention.

Rossby wave trains in the mid-latitudes, accompanied by jet stream variations, play essential roles in generating surface air temperature variations in boreal autumn and winter¹⁹⁻²⁴. Moreover, the importance of Rossby wave-induced effects has increased during the past two decades²⁵. Here we will demonstrate the importance of extratropical wave trains, specifically those originating from the North

<figure><img src="image_5.png" /></figure>

<footnote n="1">¹Frontier Science Center for Deep Ocean Multispheres and Earth System (FDOMES) and Physical Oceanography Laboratory, Ocean University of China, Qingdao, China. ²College of Oceanic and Atmospheric Sciences, Ocean University of China, Qingdao, China. ³Laoshan Laboratory, Qingdao, China.</footnote>

<footnote n="4">⁴Department of Earth and Planetary Sciences, Yale University, New Haven, USA. ⁵LOCEAN/IPSL, Sorbonne University, Paris, France. ⁶Institute for Marine and Antarctic Studies, University of Tasmania, Hobart, Tasmania, Australia. ⁷Australian Research Council Centre of Excellence for Climate Extremes, University of Tasmania, Hobart, Tasmania, Australia. ⁸State Key Laboratory of Earth Surface Processes and Resource Ecology, Beijing Normal University, Beijing, China.</footnote>

<footnote>e-mail: zhangyu@ouc.edu.cn; huangf@ouc.edu.cn</footnote>
Atlantic (NATL) and Mediterranean Sea regions, on the long-lasting ridge over the NEP, thus providing insights into the role of remote teleconnections in maintaining the NEP ridge.

## Results

Using data for 13 warm blob events that peak during the boreal cold season (from November to March; Supplementary Table 1 and Fig. 1)⁵,²⁶, we firstly consider the composite evolution of geopotential height anomalies at 500 hPa from September to the following February (Fig. 1). An anomalous ridge over the NEP emerges in September (Fig. 1a), somewhat weakens and shifts westward in October (Fig. 1b), reaches its maximum intensity in November (Fig. 1c), and weakens and moves northward in the following months (Fig. 1d–f). This ridge has a quasi-barotropic structure and extends from the surface to the upper troposphere (Supplementary Fig. 2). Aside from the ridge, a trough centered north of the Great Lakes deepens (Fig. 1 and Supplementary Fig. 2), entraining cold air over North America¹⁴,¹⁶. In November, prominent geopotential height anomalies are also found over the eastern NATL, featured by a meridional dipole (Supplementary Fig. 2). This atmospheric pattern may be related to the North Atlantic Oscillation (NAO), with 9 out of 13 warm blob events (-70%) occurring during the positive NAO phase. On the planetary scale, such ridge and trough structures resemble the tropical Northern Hemisphere pattern²⁷ in winter (Fig. 1d, e), which is potentially important for the warm blobs¹⁰. To elucidate the intensification of the ridge during the lifetime of the warm blobs, we mainly focus on the results in November hereafter due to the strongest intensity of the anomalous ridge (Fig. 1c).

To explore the potential relationship between the anomalous ridge and midlatitude wave trains, we analyze meridional wind

anomalies (shading in Fig. 2a) at 300 hPa in November for warm blobs peaking during the cold season. Two prominent wave trains appear to originate over the NATL, possibly associated with the NAO. One propagates along the subtropical jet stream with a larger wave number and another is along the subpolar jet stream with a smaller wave number. The jet streams act as waveguides for the transportation of energy downstream. The two wave trains converge near Japan, intensifying the wave amplitude and further propagating towards the NEP.

Rossby wave energy, indicated by the wave activity flux (WAF; “Methods”), greatly intensifies over the Mediterranean Sea (vectors in Fig. 2a), implying a potential Rossby wave source (RWS) in this region. Although the two wave trains somewhat weaken over the Eurasian continent, they converge and strengthen again near Japan and travel northeastward to the NEP, which appears to be important for the development and persistence of the anomalous ridge above the warm blobs (Fig. 1 and Supplementary Fig. 2).

For comparison, meridional wind anomalies and WAF are also shown for October (Supplementary Fig. 3a). Although the wave train structure is clear, it is much weaker compared to that in November (Fig. 2a). Consistently, the WAF is much weaker, especially over the Mediterranean Sea and Europe (Supplementary Fig. 3a). Note that the NATL and Mediterranean regions originated circumnavigating wave trains were not reported in Wang et al.¹⁶, but rather they emphasized the wave origin near Japan. Although similar wave trains have been documented in previous studies²³,²⁸⁻³¹, none of these focus on the wave propagation towards the NEP. Nevertheless, a recent study implied that the subtropical wave train is closely connected to intense Rossby wave activity over the Mediterranean region and anomalous warming over the NEP²⁸.

<figure><img src="image_6.png" /><figcaption>**Fig. 1** | Geopotential height anomalies at 500 hPa associated with warm blobs. a–f Composite geopotential height anomalies (shading, in gpm) at 500 hPa from September through February. Stippling indicates exceeding a 0.1 significance level</figcaption></figure>

based on the two-tailed Student's t-test. Warm blob events used for computations are listed in Supplementary Table 1. Purple boxes represent the study area for the warm blobs.
<figure><img src="image_7.png" /><figcaption>**Fig. 2 |** Rossby wave trains and their source for warm blobs peaking during the cold season. **a** Composite meridional wind anomalies (shading, in m/s), mean zonal winds (green contour, in m/s), and WAF (vectors, in m²/s²) at 300 hPa. **b** RWS (shading, in 10⁻¹¹/s²) at 300 hPa. Stippling in panels (a) and (b) indicates meridional</figcaption></figure>

wind anomalies and RWS exceeding a 0.1 significance level based on the two-tailed Student's t-test, respectively. The green and black boxes in (b) outline the Mediterranean region (35°N–45°N, 0°–30°E) and the NATL (48°N–65°N, 26°W–0°), respectively. WAF wave activity flux, RWS Rossby wave source.

Rossby waves are usually triggered when strong convection (or rainfall) generates anomalous condensation (latent) heating, inducing anomalous flow divergence and vorticity in the upper troposphere³²,³³, which often occurs over tropical oceans³⁴⁻³⁸ and the NATL³⁹⁻⁴¹. East of the NATL, Rossby wave trains often bifurcate into two branches over the Mediterranean²⁸, acting as a pivotal region for Rossby wave propagation and generation²⁹,⁴². The reasons for this wave train bifurcation are: (1) the wave train near the Mediterranean is located outside of the area with strong absolute vorticity gradient so that it is partially reflected²⁸, and (2) the wave train is sensitive to the origin of the Rossby wave³⁰.

To uncover the origin of the Rossby wave train that contributes to the NEP warm blobs, we calculate the RWS ("Methods"; shading in Fig. 2b) at 300 hPa in November for warm blobs peaking in the cold season. A large area with a positive RWS is identified over the mid-latitude eastern NATL while negative RWS values are detected over the Mediterranean region²³,⁴³. The two RWS areas correspond to anomalous convergence and divergence in the upper troposphere, respectively (not shown). We mainly focus here on these two RWS regions due to their stronger WAF. In addition, a band of positive RWS is found from southern China to south of Japan, which is consistent with the WAF trajectory along the subtropical jet (Fig. 2a). Further north, another area with a negative RWS is also significant over northern China (Fig. 2b).

In contrast, the RWSs over the NATL and Mediterranean regions are much weaker in October (Supplementary Fig. 3b) compared to November (Fig. 2b), in accordance with the weaker WAF (Supplementary Fig. 3a). Two other RWS areas over China also change greatly in location (Supplementary Fig. 3b). Although there are other centers of RWS, they are variable in time or cover a relatively small area (Fig. 2b and Supplementary Fig. 3b). Further, since the RWS accounts for the location of Rossby wave emanation induced by anomalous vertical motion, upper-tropospheric divergence, and vorticity³², we conjecture that the above-mentioned major RWS regions are potentially related to local rainfall anomalies.

To assess the rainfall variations associated with the NEP warm blobs, we estimated the rainfall anomalies for warm blobs peaking during the cold season (Fig. 3a). In November, significant negative rainfall anomalies are identified over the eastern NATL, in association with the NAO-related positive geopotential height anomalies (Supplementary Fig. 2). Moreover, we find significantly enhanced rainfall over the Mediterranean region (Fig. 3a). In fact, 10 out of the 13 warm blob events (-77%) correspond with higher-than-normal rainfall over the Mediterranean region (for example, the 2013/14, 2015, and 2019 events). To better understand the rainfall anomalies over the Mediterranean region, we show both the monthly rainfall climatology and the variability (Supplementary Fig. 4). Overall, this region has abundant rainfall from October to February, with a maximum of around 90 mm in November. However, a relatively limited amount of precipitation falls from June to August. These distinct seasonal differences are typical for the Mediterranean climate. In terms of variability, year-to-year rainfall variations are large between November and February. Rainfall anomalies in October are much weaker, especially over the Mediterranean region (Supplementary Fig. 5a), again suggesting that the Mediterranean region may play an important role in sustaining the wave train and anomalous ridge over the NEP in November. In addition, a meridional rainfall dipole (Fig. 3a) associated with the two RWS centers (Fig. 2b) is found over China and surrounding seas.

According to previous studies³²,³³, anomalous diabatic heating sustained by rainfall/convection is essential for exciting atmospheric Rossby waves. Hence, we show a composite apparent heat source ("Methods") at 300 hPa for warm blobs peaking during the cold season (Fig. 3b). Prominent heating and cooling anomalies are found over the Mediterranean and NATL regions, respectively, consistent with the rainfall anomalies (Fig. 3a). However, the two rainfall and RWS centers over China and surrounding seas are not robust in terms of anomalous heating (Fig. 3b).

Further, to verify the roles of the Mediterranean and NATL regions in triggering the wave trains that contribute to the long-lasting anomalous ridge over the NEP, we impose anomalous heating/cooling in November within an atmospheric linear baroclinic model (LBM;
"Methods")⁴⁴. Background atmospheric conditions in November given in the LBM are presented in Supplementary Fig. 6. The vertical structure of the apparent heat source over the Mediterranean region (black box in Fig. 3) exhibits anomalous heating induced by the positive rainfall anomalies throughout the whole troposphere, with a maximum magnitude at 300 hPa (Fig. 4a). After the simulation reaches a quasi-equilibrium state, a wave train originating from the Mediterranean region is identified in the extratropics with the WAF propagating eastward to the NEP (Fig. 4b). Over the NEP (purple box in Fig. 4b), an anomalous ridge is evident at 500 hPa (Fig. 4b). At the surface, the anomalous high-pressure system shifts a bit northwestward (Fig. 4c), with easterly anomalies blowing against climatological midlatitude southwesterlies (Supplementary Fig. 6c), which are favorable for inducing the NEP warm blobs by reducing latent heat loss from the ocean to the atmosphere¹,⁵,²⁶.

In terms of the detailed evolution, after imposing heating over the Mediterranean region, positive geopotential height anomalies appear over the NEP approximately 16 days after the heating and then develop and become stable in the following days (Supplementary Fig. 7). Moreover, Rossby wave energy can disperse further downstream from the anomalous ridge, amplifying the trough over North America and leading to a displacement of the polar vortex¹⁶. In addition, when the anomalous heating is centered at 700 hPa³⁴, the positive geopotential height anomalies are weaker and shift westward (Supplementary Fig. 8). Hence, the LBM experiment reproduces the anomalous ridge reasonably well over the NEP forced by anomalous heating over the Mediterranean region.

Further examining the role of the NATL forcing, we impose cold anomalies over the NATL (Fig. 4d). This results in positive geopotential height anomalies at 500 hPa over the NEP (Fig. 4e), which resemble the observations very well (Fig. 1c). However, surface westerly wind anomalies associated with a cyclonic circulation anomaly in the study

area are not favorable for the warm blobs (Fig. 4f). The difference between circulation anomalies at 500 hPa and those at the surface may arise from the baroclinic response of the LBM.

On the other hand, we also notice that the WAF converges near East Asia (Fig. 2a), which might further boost the wave train. To test this possibility, we impose anomalous cooling/heating according to the structure of the apparent heat source over this region (Supplementary Figs. 10 and 11). Although a weak anomalous ridge emerges in the heating experiment (Supplementary Fig. 11), the heating (Fig. 3b) in these two regions is scattered and not robust. In this study, we have mainly explored the roles of the Mediterranean and NATL regions, which appear to be the origin of extratropical teleconnections, rather than a midway booster, that sustain the anomalous ridge.

Next, to confirm the driving role of the Mediterranean and NATL regions in exciting the wave trains and anomalous ridge over the NEP, we conduct atmospheric model experiments using the sophisticated CAM5 (i.e., version 5.0 of the Community Atmosphere Model; "Methods"). Results show positive geopotential height anomalies over the NEP and subpolar North America in October (Fig. 5b) and November (Fig. 5c) when the anomalous warming pattern is imposed in the Mediterranean region (Fig. 5a). The location of the modeled anomalous ridge is shifted eastward compared to the observations (Fig. 1). When the sea surface temperature (SST) warming pattern associated with the largest enhanced rainfall in the Mediterranean region is prescribed, the anomalous ridge over the NEP is stronger, with a similar location to the observations (left column of Supplementary Fig. 12), suggesting that the Mediterranean region can play an important role in driving this anomalous ridge. When composite SST anomalies during the warm blobs are superimposed onto the climatology in the NATL (Fig. 5d), the anomalous ridge is robust in November (Fig. 5f). In the Exp_Natl_NAO experiment (right column of Supplementary Fig. 12), tripole-like SST anomalies with a broad cooling in the eastern NATL are

<figure><img src="image_8.png" /><figcaption>**Fig. 3** | Rainfall and diabatic heating anomalies in November for warm blobs peaking during the cold season. Composite (a) rainfall (shading, in mm) and b 300-hPa diabatic heating ($q_1$; shading, in K/day) anomalies. Stippling indicates</figcaption></figure>

exceeding a 0.1 significance level based on the two-tailed Student's t-test. The green and black boxes mark the Mediterranean and North Atlantic, respectively.
<figure><img src="image_9.png" /></figure>

**Fig. 4** | Steady atmospheric circulation response to anomalous heating/cooling over the Mediterranean and NATL regions in the LBM. **a, d** Vertical profiles of $q_1$ anomaly (blue line) and the imposed heating (red line) in the LBM over the Mediterranean and NATL regions, respectively. Note that the heating was multiplied by a factor of two to obtain a more robust response. **b** Geopotential height anomalies (contours, in gpm) at 500 hPa and WAF (vectors, in m²/s²) at 300 hPa and **c** SLP

anomalies (contours, in hPa) and surface winds anomalies (vectors, in m/s) in the Mediterranean experiment. Shading indicates the imposed heating. Purple boxes represent the study area for the warm blobs. **e, f** Same as (**b, c**), but for the NATL experiment. HGT geopolitical height, NATL North Atlantic, LBM linear baroclinic model, WAF wave activity flux, SLP sea level pressure.

prescribed, associated with the positive phase of the NAO. A prominent anomalous ridge over the NEP and polar regions is found in October and November, resembling the observations (Fig. 1). Thus, the effects of both the Mediterranean and NATL regions on the NEP anomalous ridge are captured by the CAM5 simulations.

## Discussion

In summary, this study provides insights into the development of the anomalous NEP ridge from the perspective of extratropical wave trains, highlighting downstream climate impacts reaching the NEP from the Mediterranean and NATL regions. The wave train dynamics could provide a potentially important source of predictability for the anomalous ridge and the resultant cold-season NEP warm blobs as well as North American temperature anomalies. We emphasize that this mechanism prominent in November may not be applicable for other winter months because the background state plays a crucial role in the generation and guidance of Rossby waves, influencing the establishment of teleconnection patterns. Moreover, the contribution from the tropics should not be neglected, as it could be amplified when combined with extratropical forcing (e.g., from the Mediterranean; Supplementary Fig. 13).

## Methods

### Observational and reanalysis datasets

For SST, we use monthly data from the National Oceanic and Atmospheric Administration (NOAA) Extended Reconstructed SST version 5 (ERSST v5), gridded at 2° × 2°⁴⁵. For precipitation, wind, SLP, temperature, and geopotential height, we use the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis⁴⁶. The horizontal resolution of ERA5 data is 0.25° × 0.25°. The data prior to 1959 are from an ERA5 preliminary version⁴⁷. We also use velocity potential, divergence, and relative vorticity data from the Japanese 55-year Reanalysis (JRA-55) at 2.5° × 2.5° resolution⁴⁸. We analyze the period from 1951 to 2021 in this study. Anomalies of variables are calculated relative to a baseline climatology from 1979 to 2021.

### Composite analysis of the NEP warm blobs

The study area of the warm blobs is 40°N–50°N, 160°W–135°W. Following the definition from previous studies⁵,²⁶, a warm blob event is identified when normalized monthly SST anomalies averaged over the study area (i.e., the blob index) are greater than 1.0 for no fewer than 5 months with at most 1-month interruption. Accordingly, 13 warm blob events (Supplementary Table 1) peaking during the cold season
<figure><img src="image_10.png" /><figcaption>Fig. 5 | Response of geopotential height to the prescribed SST forcing in CAM5. a, d Prescribed SST anomalies (shading, in °C) relative to climatology in Exp_Med and Exp_Natl, respectively. b, c Geopotential height anomalies (shading, in gpm) at</figcaption></figure>

300 hPa in Exp_Med relative to CTRL in October and November. e, f As in (b, c), but for Exp_Natl. SST sea surface temperature, CAM5 Community Atmosphere Model 5.0.

(from November to March) are selected to perform our composite analysis. The blob index is calculated after removing linear trends$^{5,26}$; such an approach is widely used, although the SST variations do have nonlinear features$^{49–52}$. Peak intensity is given by the blob index when reaching the maximum. The return period is computed according to refs. 7,53. Note that due to the period of the JRA-55 dataset, we use 11 warm blob events after 1958 when calculating the RWS.

## Calculation of WAF

To analyze the propagation of Rossby wave energy in the Northern Hemisphere, the horizontal WAF is computed following Takaya and Nakamura$^{54}$:

$$W = \frac{P \cos \phi}{2|U|} \cdot \begin{pmatrix} \frac{U}{a^2 \cos^2 \phi} \left[ \left( \frac{\partial \psi}{\partial \lambda} \right)^2 - \psi' \frac{\partial^2 \psi}{\partial \lambda^2} \right] \\ + \frac{V}{a^2 \cos \phi} \left[ \frac{\partial \psi'}{\partial \lambda} \frac{\partial \psi}{\partial \phi} - \psi' \frac{\partial^2 \psi}{\partial \lambda \partial \phi} \right], \\ + \frac{U}{a^2 \cos \phi} \left[ \frac{\partial \psi'}{\partial \lambda} \frac{\partial \psi}{\partial \phi} - \psi' \frac{\partial^2 \psi}{\partial \lambda \partial \phi} \right] \\ + \frac{V}{a^2} \left[ \left( \frac{\partial \psi}{\partial \phi} \right)^2 - \psi' \frac{\partial^2 \psi}{\partial \phi^2} \right] \end{pmatrix} \quad (1)$$

where W (unit: m$^2$ s$^{-2}$) denotes the horizontal WAF; P is pressure/ (1000 hPa); U (unit: m s$^{-1}$) represents the basic flow with the zonal component U and meridional component V; $\psi(=\Phi/f)$ is the stream function; a is the radius of the Earth; $\phi$ and $\lambda$ denote latitude and longitude, respectively; $\Phi$ is geopotential height, and $f=2\Omega \sin \phi$ is the Coriolis parameter; where $\Omega$ is the rotation rate of the Earth. Primes indicate anomalies relative to climatology.

## Calculation of RWS

To locate the regions that can trigger Rossby waves in the upper troposphere, the Rossby wave source (RWS) is calculated at 300 hPa following Sardeshmukh and Hoskins$^{32}$:

$$RWS = -\nabla_H \cdot \{u_X(f + \zeta)\}' \quad (2)$$

where $u=(u, v)$ denotes the horizontal wind velocity vector and subscript $\chi$ indicates its divergent component. $\zeta$ and $f$ are relative vorticity and planetary vorticity, respectively; $\nabla_H$ is the horizontal gradient operator. Primes indicate anomalies.

## Calculation of the apparent heat source

To isolate the anomalous heating and cooling pattern associated with precipitation anomalies, we calculate the atmospheric diabatic heating ($q_1$) as below$^{55}$ based on data from the ERA5 dataset:

$$q_1 = \frac{\partial T}{\partial t} - (\omega \sigma - \mathbf{V} \cdot \nabla_H T) \quad (3)$$

where T is air temperature; and $\omega$ is the vertical velocity in pressure coordinates. Static stability is expressed as $\sigma = (RT/c_p p) - (\partial T / \partial p)$, where $c_p$ is the specific heat capacity of air at constant pressure; R denotes the specific gas constant, and p is pressure. V = (u, v) is the horizontal wind vector, and $\nabla$ is the horizontal gradient operator. Here, $q_1$ represents total diabatic heating, including radiation, latent and sensible heat fluxes, as well as subgrid-scale heat flux convergence$^{55}$.

## Experimental setup in the LBM

The atmospheric linear baroclinic model (LBM) is widely used to explore linear atmospheric dynamics$^{44}$. The LBM is based on a linearized version of the primitive equations. Given a basic state $\bar{X}$, a steady response $X$ follows an equation written in a matrix form as

$$L(\bar{X})X = \mathbf{F} \quad (4)$$

where $\mathbf{F}$ indicates a forcing vector, and $L$ is a linear dynamical operator related to the perturbed primitive equations$^{44}$. In this study, the LBM is utilized to examine whether prescribing heating/cooling over the Mediterranean and NATL can drive a Rossby wave train that sustains an atmospheric ridge over the NEP in accordance with the observations. The model is run at T42L20, with a nominal horizontal resolution of 2.8° × 2.8° and 20 vertical sigma levels. The model consists of basic
equations linearized with respect to the atmospheric mean state in November derived from the NCEP reanalysis²⁶. Dissipation in the model includes (1) a bi-harmonic horizontal diffusion with a damping timescale of 1 day for the shortest waves, (2) a weak vertical diffusion (a damping timescale of 1000 days) to remove noise arising from finite differencing, (3) Newtonian damping, and (4) Rayleigh friction represented by a linear drag, with a timescale of 1 day applied only to the lower boundary layers and the uppermost two levels⁵⁷. The simulations are run for 20 days to reach a quasi-equilibrium atmospheric state⁵⁷⁻⁵⁹. Output variables, including geopotential height, winds, and sea level pressure, are averaged over days 16–18 of the simulations, which yields a steady response to the prescribed heating/cooling. Note that the imposed heating/cooling in the Mediterranean and NATL experiments is multiplied by a factor of two to obtain a robust response as we calculate mean values over these regions.

## Experimental setup in CAM5

CAM5 model was developed by NCAR. Five experiments were conducted: a control experiment (CTRL), two Mediterranean experiments (Exp_Med and Exp_Med_rain), and two NATL experiments (Exp_Natl and Exp_Natl_NAO; Supplementary Table 2). We use the f19_g16 model configuration (i.e., 2.5° × 1.9° atmospheric resolution) for the experiments. The prescribed SST⁶⁰ anomalies are computed from the merged products of HadISST version 1⁶¹ and NOAA OISST version 2⁶². The CTRL run is forced by monthly climatological SST globally over 1981–2010. For the Exp_Med run, composite SST anomalies during the warm blobs are superimposed onto climatological SST in the Mediterranean region (Fig. 5a). For the Exp_Med_rain run, SST anomalies averaged in 2013 and 2019 (when the Mediterranean has the largest enhanced rainfall) are superimposed onto climatological SST in the Mediterranean region (Supplementary Fig. 10a). For the Exp_Natl run, composite SST anomalies during the warm blobs are superimposed onto climatological SST in the NATL (Fig. 5d). For the Exp_Natl_NAO run, SST anomalies in 1993 with the largest NAO are superimposed onto climatological SST in the NATL (Supplementary Fig. 10d). For these sensitivity experiments, SST climatology is prescribed in the regions outside of the above forcing area. Each experiment is initialized from January 1, 1979 and integrated for 25 years, with the last 20 years averaged for our analyses.

## Data availability

The ERSS v5 data are available at https://psl.noaa.gov/data/gridded/data.noaa.ersst.v5.html. The ERA5 data can be accessed at https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels-monthly-means?tab form, https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels-monthly-means?tab=overview, https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels-monthly-means-preliminary-back-extension?tab=overview, and https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels-monthly-means-preliminary-back-extension?tab=overview. The JRA-55 data are available at https://rda.ucar.edu/datasets/ds628.1/index.html#cgi-bin/datasets/getWebList?dsnum=628.1&action=customize&disp=. The model output and source data used in this study are available on Zenodo at https://zenodo.org/doi/10.5281/zenodo.10473150.

## Code availability

This article uses Python (https://anaconda.org/) and NCL to perform the reported analyses. The codes that support the findings of this study are available at https://doi.org/10.5281/zenodo.10494745.

## References

1. Bond, N. A., Cronin, M. F., Freeland, H. & Mantua, N. Causes and impacts of the 2014 warm anomaly in the NE Pacific. Geophys. Res. Lett. 42, 3414–3420 (2015).

2. Di Lorenzo, E. & Mantua, N. Multi-year persistence of the 2014/15 North Pacific marine heatwave. Nat. Clim. Change 6, 1042–1047 (2016).

3. Hu, Z.-Z., Kumar, A., Jha, B., Zhu, J. & Huang, B. Persistence and predictions of the remarkable warm anomaly in the northeastern Pacific Ocean during 2014–16. J. Clim. 30, 689–702 (2017).

4. Amaya, D. J., Miller, A. J., Xie, S. P. & Kosaka, Y. Physical drivers of the summer 2019 North Pacific marine heatwave. Nat. Commun. 11, 1903 (2020).

5. Chen, Z., Shi, J., Liu, Q., Chen, H. & Li, C. A persistent and intense marine heatwave in the Northeast Pacific during 2019–2020. Geophys. Res. Lett. 48, e2021GL093239 (2021).

6. Scannell, H. A., Johnson, G. C., Thompson, L., Lyman, J. M. & Riser, S. C. Subsurface evolution and persistence of marine heatwaves in the Northeast Pacific. Geophys. Res. Lett. 47, e2020GL090548 (2020).

7. Laufkötter, C., Zscheischler, J. & Frölicher, T. L. High-impact marine heatwaves attributable to human-induced global warming. Science 369, 1621–1625 (2020).

8. Wang, S. et al. Changing ocean seasonal cycle escalates destructive marine heatwaves in a warming climate. Environ. Res. Lett. 17, 054024 (2022).

9. Holbrook, N. J. et al. Keeping pace with marine heatwaves. Nat. Rev. Earth Environ. 1, 482–493 (2020).

10. Liang, Y.-C., Yu, J.-Y., Saltzman, E. S. & Wang, F. Linking the tropical Northern Hemisphere pattern to the Pacific warm blob and Atlantic cold blob. J. Clim. 30, 9041–9057 (2017).

11. Schmeisser, L., Bond, N. A., Siedlecki, S. A. & Ackerman, T. P. The role of clouds and surface heat fluxes in the maintenance of the 2013–2016 Northeast Pacific marine heatwave. J. Geophys. Res. Atmos. 124, 10772–10783 (2019).

12. Zhi, H., Lin, P., Zhang, R.-H., Chai, F. & Liu, H. Salinity effects on the 2014 warm "Blob" in the Northeast Pacific. Acta Oceanol. Sin. 38, 24–34 (2019).

13. Anderson, B. T., Gianotti, D. J. S., Furtado, J. C. & Di Lorenzo, E. A. decadal precession of atmospheric pressures over the North Pacific. Geophys. Res. Lett. 43, 3921–3927 (2016).

14. Hartmann, D. L. Pacific sea surface temperature and the winter of 2014. Geophys. Res. Lett. 42, 1894–1902 (2015).

15. Seager, R. et al. Causes of the 2011–14 California drought. J. Clim. 28, 6997–7024 (2015).

16. Wang, S. Y., Hipps, L., Gillies, R. R. & Yoon, J. H. Probable causes of the abnormal ridge accompanying the 2013–2014 California drought: ENSO precursor and anthropogenic warming footprint. Geophys. Res. Lett. 41, 3220–3226 (2014).

17. Holbrook, N. J. et al. A global assessment of marine heatwaves and their drivers. Nat. Commun. 10, 2624 (2019).

18. Shi, J., Fedorov, A. V. & Hu, S. North Pacific temperature and precipitation response to El Niño-like equatorial heating: sensitivity to forcing location. Clim. Dyn. 53, 2731–2741 (2019).

19. Rodwell, M. J. & Hoskins, B. J. Monsoons and the dynamics of deserts. Q. J. R. Meteorol. Soc. 122, 1385–1404 (1996).

20. Rodwell, M. J. & Hoskins, B. J. Subtropical anticyclones and summer monsoons. J. Clim. 14, 3192–3211 (2001).

21. Branstator, G. Circumglobal teleconnections, the jet stream waveguide, and the North Atlantic Oscillation. J. Clim. 15, 1893–1910 (2002).

22. Wakabayashi, S. & Kawamura, R. Extraction of major teleconnection patterns possibly associated with the anomalous summer climate in Japan. J. Meteorol. Soc. Jpn. 82, 1577–1588 (2004).

23. Watanabe, M. Asian jet waveguide and a downstream extension of the North Atlantic Oscillation. J. Clim. 17, 4674–4691 (2004).

24. Li, C. & Sun, J. Role of the subtropical westerly jet waveguide in a southern China heavy rainstorm in December 2013. Adv. Atmos. Sci. 32, 601–612 (2015).
25. Kornhuber, K., Petoukhov, V., Petri, S., Rahmstorf, S. & Coumou, D. Evidence for wave resonance as a key mechanism for generating high-amplitude quasi-stationary waves in boreal summer. *Clim. Dyn.* **49**, 1961–1979 (2017).

26. Chen, Z., Shi, J. & Li, C. Two types of warm blobs in the Northeast Pacific and their potential effect on the El Niño. *Int. J. Climatol.* **41**, 2810–2827 (2021).

27. Mo, K. C. & Livezey, R. E. Tropical-extratropical geopotential height teleconnections during the Northern Hemisphere winter. *Mon. Weather Rev.* **114**, 2488–2515 (1986).

28. An, X. et al. Possible dynamic mechanisms of high- and low-latitude wave trains over Eurasia and their impacts on air pollution over the North China Plain in early winter. *J. Geophys. Res. Atmos.* **127**, e2022JD036732 (2022).

29. Hu, K., Huang, G., Wu, R. & Wang, L. Structure and dynamics of a wave train along the wintertime Asian jet and its impact on East Asian climate. *Clim. Dyn.* **51**, 4123–4137 (2018).

30. Jiao, Y., Wu, R. & Song, L. Individual and combined impacts of two Eurasian wave trains on intraseasonal East Asian winter monsoon variability. *J. Geophys. Res. Atmos.* **124**, 4530–4548 (2019).

31. Song, J., Li, C. & Zhou, W. High and low latitude types of the downstream influences of the North Atlantic Oscillation. *Clim. Dyn.* **42**, 1097–1111 (2014).

32. Sardeshmukh, P. D. & Hoskins, B. J. The generation of global rotational flow by steady idealized tropical divergence. *J. Atmos. Sci.* **45**, 1228–1251 (1988).

33. Qin, J. & Robinson, W. A. On the Rossby wave source and the steady linear response to tropical forcing. *J. Atmos. Sci.* **50**, 1819–1823 (1993).

34. Hoskins, B. J. & Karoly, D. J. The steady linear response of a spherical atmosphere to thermal and orographic forcing. *J. Atmos. Sci.* **38**, 1179–1196 (1981).

35. Enomoto, T., Hoskins, B. J. & Matsuda, Y. The formation mechanism of the Bonin high in August. *Q. J. R. Meteorol. Soc.* **129**, 157–178 (2003).

36. Hu, S. & Fedorov, A. V. Indian Ocean warming can strengthen the Atlantic meridional overturning circulation. *Nat. Clim. Change* **9**, 747–751 (2019).

37. Hu, S. & Fedorov, A. V. Indian Ocean warming as a driver of the North Atlantic warming hole. *Nat. Commun.* **11**, 4785 (2020).

38. Franzke, C., Feldstein, S. B. & Lee, S. Synoptic analysis of the Pacific-North American teleconnection pattern. *Q. J. R. Meteorol. Soc.* **137**, 329–346 (2011).

39. Wallace, J. M. & Gutzler, D. S. Teleconnections in the geopotential height field during the Northern Hemisphere winter. *Mon. Weather Rev.* **109**, 784–812 (1981).

40. Czaja, A. & Frankignoul, C. Observed impact of Atlantic SST anomalies on the North Atlantic Oscillation. *J. Clim.* **15**, 606–623 (2002).

41. Li, S. & Bates, G. T. Influence of the Atlantic Multidecadal Oscillation on the winter climate of East China. *Adv. Atmos. Sci.* **24**, 126–135 (2007).

42. An, X. et al. The combined effect of two westerly jet waveguides on heavy haze in the North China Plain in November and December 2015. *Atmos. Chem. Phys.* **20**, 4667–4680 (2020).

43. Nie, Y., Zhang, Y., Yang, X. Q. & Ren, H. L. Winter and summer Rossby wave sources in the CMIP5 models. *Earth Space Sci.* **6**, 1831–1846 (2019).

44. Watanabe, M. & Kimoto, M. Atmosphere-ocean thermal coupling in the North Atlantic: a positive feedback. *Q. J. R. Meteorol. Soc.* **126**, 3343–3369 (2000).

45. Huang, B. et al. Extended reconstructed sea surface temperature version 5 (ERSSTv5): upgrades, validations, and intercomparisons. *J. Clim.* **30**, 8179–8205 (2017).

46. Hersbach, H. et al. The ERA5 global reanalysis. *Q. J. R. Meteorol. Soc.* **146**, 1999–2049 (2020).

47. Bell, B. et al. The ERA5 global reanalysis: preliminary extension to 1950. *Q. J. R. Meteorol. Soc.* **147**, 4186–4227 (2021).

48. Kobayashi, S. et al. The JRA-55 Reanalysis: general specifications and basic characteristics. *J. Meteorol. Soc. Jpn.* **93**, 5–48 (2015).

49. Mantua, N. J., Hare, S. R., Zhang, Y., Wallace, J. M. & Francis, R. C. Pacific interdecadal climate oscillation with impacts on salmon production. *Bull. Am. Meteorol. Soc.* **78**, 1069–1079 (1997).

50. Trenberth, K. & Shea, D. Atlantic Hurricanes and natural variability in 2005. *Geophys. Res. Lett.* **33**, L12704 (2006).

51. Zhang, Y., Wallace, J.-M. & Battisti, D.-S. ENSO-like interdecadal variability: 1900–93. *J. Clim.* **10**, 1004–1020 (1996).

52. Silva, E. N. S. & Anderson, B. T. Northeast Pacific marine heatwaves linked to Kuroshio Extension variability. *Commun. Earth Environ.* **4**, 367 (2023).

53. Volpi, E., Fiori, A., Grimaldi, S., Lombardo, F. & Koutsoyiannis, D. One hundred years of return period: strengths and limitations. *Water Resour. Res.* **51**, 8570–8585 (2015).

54. Takaya, K. & Nakamura, H. A formulation of a phase-independent wave-activity flux for stationary and migratory quasigeostrophic eddies on a zonally varying basic flow. *J. Atmos. Sci.* **58**, 608–627 (2001).

55. Yanai, M., Esbensen, S. & Chu, J. H. Determination of bulk properties of tropical cloud clusters from large-scale heat and moisture budgets. *J. Atmos. Sci.* **30**, 611–627 (1973).

56. Kanamitsu, M. et al. NCEP-DOE AMIP-II reanalysis (R-2). *Bull. Am. Meteorol. Soc.* **83**, 1631–1644 (2002).

57. Watanabe, M. & Jin, F.-F. A moist linear baroclinic model: coupled dynamical-convective response to El Niño. *J. Clim.* **16**, 1121–1139 (2003).

58. An, X. et al. Effect of rainfall-induced diabatic heating over southern China on the formation of wintertime haze on the North China Plain. *Atmos. Chem. Phys.* **22**, 725–738 (2022).

59. Lu, R. & Lin, Z. Role of subtropical precipitation anomalies in maintaining the summertime meridional teleconnection over the Western North Pacific and East Asia. *J. Clim.* **22**, 2058–2072 (2009).

60. Hurrell, J. W., Hack, J. J., Shea, D., Caron, J. M. & Rosinski, J. A new sea surface temperature and sea ice boundary dataset for the Community Atmosphere Model. *J. Clim.* **21**, 5145–5153 (2008).

61. Rayner, N. A. et al. Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century. *J. Geophys. Res. Atmos.* **108**, 4407 (2003).

62. Reynolds, R. W. & Smith, T. M. Improved global sea surface temperature analyses using optimum interpolation. *J. Clim.* **7**, 929–948 (1994).

## Acknowledgements

This work is supported by the National Natural Science Foundation of China (42006013, 42230405, and 42276006 to J.S.; 42075024 to F.H.), Laoshan Laboratory (No. LSKJ202202503 to J.S.; No. LSKJ202202602 to Y.Z.), the National Key Research and Development Program of China (2019YFA0607004 to F.H.), China Postdoctoral Science Foundation (2021M703034 to Y.Z.), the ARC Centre of Excellence for Climate Extremes (CE170100023 to N.J.H.), NASA (80NSSC21K0558 to A.V.F.), NOAA (NA20OAR4310377 to A.V.F.), and the ARCHANGE project (ANR-18-MPGA-0001, France to A.V.F.). We greatly thank Marine Big Data Center of the Institute for Advanced Ocean Study of Ocean University of China for providing computational resources. We thank Prof. Shang-Ping Xie for his insightful comments and edits for this paper. We also thank Drs. Hong Wang, Boqi Liu, Zedong Liu, Peng Liu, Xiadong An, and Tengfei Yu for their helpful suggestions and comments on model set-up, experimental design, and result interpretation.
**Author contributions**

J.S. conceived the study and wrote the manuscript. H.H. plotted the figures and conducted the analyses. F.H., A.V.F., N.J.H., and Y.Z. provided additional thoughts and helped revise the manuscript. Y.L., J.C., and X.H. helped improve the figures. R.D., S.W., Q.L., and X.L. were involved in interpreting the results.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-024-47032-x.

**Correspondence** and requests for materials should be addressed to Yu Zhang or Fei Huang.

**Peer review information** *Nature Communications* thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024