Tropospheric pathways of the late‑winter ENSO teleconnection to Europe

The late-winter signal associated with the El Niño-Southern Oscillation (ENSO) over the European continent is unsettled. Two main anomalous patterns of sea-level pressure (SLP) can be identified: a “wave-like” pattern with two opposite-signed anomalies over Europe, and a pattern showing a single anomaly (“semi-isolated”). In this work, potential paths of the tropospheric ENSO teleconnection to Europe and their role in favoring a more wave-like or semi-isolated pattern are explored. Outputs from historical runs of two versions of the MPI-ESM coupled model, which simulate these two types of patterns, are examined. A novel ray-tracing approach that accounts for zonal asymmetries in the background flow is used to test potential propagation paths in these simulations and in observations; three source regions are considered: the tropical Pacific, the North America/North Atlantic, and the tropical Atlantic. The semi-isolated pattern is suggested to be related to the well-known Rossby wave train emanating from the tropical Pacific, either via a split over northern North America or via reflection due to inhomogeneities in the background flow. The wave-like pattern, in turn, appears to be related to a secondary wave train emerging from the tropical Atlantic. The competition between these two pathways contributes to determining the actual surface response.


Introduction
El Niño-Southern Oscillation (ENSO) has been associated with a late-winter (January-March) signal in the North Atlantic-European (NAE) region, a "canonical" sea-level pressure (SLP) dipole between middle and high latitudes (e.g.Brönnimann 2007).The western part of this dipole, which is located over the North Atlantic, is robust and has been shown to be mostly driven by tropospheric processes (e.g.Mezzina et al. 2020Mezzina et al. , 2021a)).In contrast, controversy still exists concerning the eastern part of the signal, located over the European continent, since disagreement is present in both observations (Fig. 1, top) and models (Fig. 1, bottom), also depending on the methodology and season used.In some cases, a "wave-like" pattern with two anomalies of alternating sign over Europe is present (Fig. 1a, e.g.Toniazzo and Scaife 2006;Hardiman et al. 2019).In others, a pattern with a single anomaly, negative for El Niño (e.g.Fraedrich and Müller 1992), is visible, sometimes appearing as an extension of the mid-latitude lobe of the canonical dipole (Fig. 1d, e), sometimes as a detached center (Fig. 1b,  c, f); hereafter, we will refer to this pattern as "semi-isolated".The semi-isolated pattern is usually accompanied by an upper-level anomaly of the same sign over northern Europe (e.g.Blackmon et al. 1983;Brönnimann 2007;Brönnimann et al. 2007;García-Serrano et al. 2011;Mezzina et al. 2020), but the relationship between the lower-level and upper-level signatures has not been settled, nor their overall nature, robustness and dynamics.
It is well-known that the dominant feature of the ENSO teleconnection to the Northern Hemisphere (NH) extra-tropics is a large-scale tropospheric Rossby wave train emanating from the tropical Pacific and propagating at upper levels, with a first center of action over the Aleutian Low (cyclonic for El Niño), a second one of opposite sign over Canada, and finally, as it bends south-eastward, a tail over eastern North America (again cyclonic for El Niño; see Trenberth et al. 1998 for a review).The surface projection of this tail has been related to the western portion of the canonical dipole with a vertically-tilted structure (e.g.Mezzina et al. 2020Mezzina et al. , 2021a)).Similarly, the signal over Europe may be related to this "main" wave train through a split occurring over northern North America, as suggested by García-Serrano et al. (2011) following theoretical considerations by Hoskins and Karoly (1981) and Karoly et al. (1989), or via reflection over eastern North America by zonal inhomogeneities in the background flow (e.g.Branstator (1983Branstator ( , 1985)); Trenberth et al. (1998) and references therein).A secondary wave train is also a plausible hypothesis: Toniazzo and Scaife (2006) identified, for strong El Niños, a source region in the tropical Atlantic possibly emerging from the response of the Walker circulation (see also García-Serrano et al. 2017), and linked it to the wave-like pattern over Europe.Additional ENSOrelated wave sources have also been detected over the Gulf of Mexico-Caribbean Sea region (e.g.Rodríguez-Fonseca et al. 2016;Fereday et al. 2018;Ayarzagüena et al. 2018;Hardiman et al. 2019), suggesting the presence of a third wave train.On the other hand, the semi-isolated pattern has been proposed to be linked to the downward propagation of ENSO anomalies from the polar stratosphere (e.g.Cagnazzo and Manzini 2009).While we acknowledge a possible role of the stratosphere (see Mezzina et al. 2021b), in this work we pursue a comprehensive understanding of the ENSO-NAE teleconnection in terms of tropospheric dynamics, which appear to dominate in the western North Atlantic.
We will explore potential paths for the tropospheric ENSO teleconnection to Europe and try to reconcile the different observed/simulated responses with theoretical basis of linear Rossby wave propagation using a novel ray tracing approach.Ray tracing is a concept borrowed from geometrical optics to describe the propagation of wave energy along "rays" aligned with the local group velocity.It was first consistently applied to atmospheric Rossby waves by Hoskins and Karoly (1981), who developed it in the framework of a zonally-symmetric basic state, a strong assumption to considerably simplify the equations.It was soon pointed out, however, that zonal asymmetries in the background flow, particularly those related to the local jets, can affect propagation by reflecting the wave trains (see Trenberth et al. 1998 for a review), but the extension of the theory to a "fully wavy" (i.e. in both zonal and meridional directions) mean flow is extremely complex ( Karoly 1983).On the other hand, as noted by Branstator (1983), a longitudinallyvarying flow can be viewed as a series of sub-regions in which the zonally-symmetric theory is locally valid, an argument used by Hoskins and Ambrizzi (1993) to qualitatively interpret results from a barotropic model with an observed mean state in terms of ray propagation.In this study, we have further extended this approach to develop a ray tracing algorithm that locally calculates the group velocity from a realistic, horizontally inhomogeneous flow, but in which reflection is treated similarly to the zonally-symmetric case.This approach, that is novel to the best of our knowledge, allows us to visually represent ray paths in a realistic flow, overcoming the limitations of using a strictly zonally-symmetric flow, albeit in a simplified manner.
We will consider outputs from two versions of the same state-of-the-art coupled model, both with a realistic mean flow, in which the ENSO response over Europe is different, and examine potential tropospheric pathways by using our ray-tracing approach and launching rays from several regions to test the various hypotheses related to the main and secondary wave trains.Our results may help understand the processes relevant for the ENSO teleconnection to the European continent and highlight which models' aspects need more attention in order for this teleconnection to be properly simulated and predicted (e.g.Dawson et al. 2011;Li et al. 2020).

Models, data and methods
We examine outputs from the CMIP5 historical runs (1850-2005, 3 members) of two versions of the MPI-ESM coupled model, with same horizontal resolution in the atmosphere (T63/1.9• ), same top (0.01 hPa), but different vertical resolution: 47 levels in the low-resolution (LR) version and 95 in the mixed-resolution (MR) one (Giorgetta et al. 2013).Despite the different vertical resolution, both versions have been shown to capture accurately the mean state of the stratospheric polar vortex (Butler et al. 2016).Further details on the model, which has been long used as a seasonal forecast system, can be found in Baehr et al. (2015) and Domeisen et al. (2015) for LR, and Dobrynin et al. (2018) for MR.
The forced ENSO response is estimated by computing member-concatenated linearly detrended anomalies and applying linear regression onto the Nião3.4(N3.4) index obtained from sea surface temperature (SST).Results from the models are compared in Sect. 4 with observational ones using the NOAA-20CR reanalysis (Compo et al. 2011).January-to-March (JFM) is the target season.Statistical significance is assessed with a two-tailed t-test at the 95% confidence level.

"Hybrid" ray tracing
We consider large-scale, stationary Rossby waves and assume a constant zonal wavenumber k, as in the zonallysymmetric case (e.g.Hoskins and Ambrizzi 1993).The horizontal components of the group velocity, which determine the direction of the ray, are computed at each step from the climatological zonal wind U as: where * = − 2 U y 2 is the meridional gradient of absolute vorticity, and K 2 s = k 2 + l 2 = * ∕U is the total stationary wavenumber.As observed by Hoskins and Karoly (1981), the slope of the ray is proportional to the meridional wavenumber l and inversely proportional to the zonal wavenumber k, so that shorter waves (larger k) are expected to have a more zonal trajectory.
In the framework of a zonally-averaged U , while the waves are generally refracted towards higher values of K s , several special cases can be identified (Hoskins and Ambrizzi 1993): all rays shall turn before reaching latitudes at which * (and hence K s ) is zero and the theory cannot be applied close to latitudes with U = 0 (hence K s → ∞ ); in fact, note also that for K s to be real, positive * and U are required.Furthermore, lines with K s = k and l = 0 act as turning latitudes beyond which the rays would decay.Therefore, in the zonally-symmetric case, maps of l 2 (for a fixed k) feature latitudes with l 2 = 0 that act as barriers and reflect all the rays (Fig. 5, bottom, for k = 3 ; see also Fig. S2 in the Supplementary Material for k = 2, 4 ).However, when l 2 is computed from the actual mean flow, 2D irregular regions of negative l 2 appear, instead of simple turning latitudes (Fig. 5, top; see also Fig. S1), and it is not straightforward to determine the ray trajectories in the proximity of these "forbidden" areas.In particular, the theory does not predict the behaviour of a ray approaching a forbidden zone from an west-east direction, since after simple meridional reflection the ray may still be inside the region of negative l 2 , unlike in the zonally-symmetric case. (1) To overcome this problem, our hybrid approach consists in treating the forbidden regions as if they were local turning latitudes from the zonally-symmetric case.Namely, we first build "adjusted" maps of positive l 2 by replacing point-by- point the negative values with the average from the nearest neighbors with l 2 > 0 , thus obtaining maps of positive l 2 in which the rays would be allowed to propagate freely, with no forbidden regions.Then, we insert latitudinal-but zonallyasymmetric-barriers to roughly represent the original negative regions (for k = 3 , compare the hatched regions and the thick horizontal lines in Fig. 6; see Fig. S1 for k = 2, 4 ).That is, a ray is free to propagate everywhere, with group velocity computed from the positively-adjusted l 2 , until it eventu- ally hits a barrier and undergoes total reflection (same angle of incidence).Despite its simplicity, this approach allows to bypass the unclear aspects of a ray entering an irregular region of negative l 2 , while maintaining some important aspects related to the zonally-asymmetric mean flow (see Sect. 3.3).
The basic state used in the analysis is the JFM climatology of the 200-hPa zonal wind smoothed with longitudinal running windows of 15 • .Computations are carried out in spherical coordinates, since the ray trajectory and curvature do not depend on the projection used ( Hoskins and Karoly 1981;Hoskins and Ambrizzi 1993).

Extra-tropical response
The deepening of the Aleutian Low related to the main wave train in the North Pacific is captured by both models (Fig. 2, bottom), in agreement with observations (e.g.Zhang et al. 2016;Mezzina et al. 2020).In contrast, the SLP dipole over the North Atlantic is fully present only in MR (Fig. 2c), whereas in LR it is confined to North America, west of 60°W (Fig. 2d).Over the European continent, the extended negative lobe of the North Atlantic dipole dominates in MR (Fig. 2c), yielding an semi-isolated pattern, while in LR a more wave-like pattern is present, since a positive anomaly centered over the Mediterranean basin is accompanied by a weaker and non-significant anomaly of opposite sign over north-eastern Europe (Fig. 2d).This difference is partly reflected at upper levels, in the 200-hPa geopotential height  2b), which is also present in MR but shifted to the west and confined to lower latitudes (Fig. 2a).
Moreover, a small but significant negative anomaly is evident over northern Europe in MR, roughly at the same location of the SLP anomaly (cf.Fig. 2a and c), as observed in previous studies (e.g.Blackmon et al. 1983;Brönnimann 2007;García-Serrano et al. 2011), suggesting a relationship between this upper-level center of action and the SLP semi-isolated pattern.In contrast, the upperlevel response in the North Pacific and North America/ western North Atlantic regions differs very little between the two models (Fig. 2, top).In LR, it is slightly shifted westward, particularly the negative center over North America, which may explain why the SLP signal in the Atlantic is limited to the western part of the basin, assuming that the mid-latitude lobe of the canonical SLP dipole is related to the tail of the main wave train (e.g.Mezzina et al. 2020Mezzina et al. , 2021a)).
These two versions of MPI-ESM constitute a good framework to investigate the reasons for the different ENSO signals reported over Europe, since they can be used to represent a more semi-isolated response, in the case of MR, or a wave-like pattern, in the case of LR.

Forcing, tropical response and mean flow
A distinct extra-tropical ENSO response in coupled simulations such as the ones examined could arise from variations in the oceanic forcing, but it does not appear to be the case here, since the SST patterns associated with the N3.4 index are very similar (Fig. 3  Concerning a possible stratospheric influence, the ENSO response in the lower stratosphere is comparable in the two models (Fig. 4, bottom), which both show realistic (cf.Domeisen et al. 2019;Mezzina et al. 2020) anomalies in the 50-hPa geopotential height-although less significant in LR-that can be interpreted as a weakening and displacement of the polar vortex towards the North Atlantic sector (e.g.Mezzina et al. 2021b).Note also that the response in the middle stratosphere (10 hPa) is less significant in MR (Fig. 4a ), and not significant over the northern North Atlantic in either of the two models (Fig. 4, top).Regarding a possible delayed impact, both models show no significant response in November-December (not shown).We thus also discard stratospheric processes as a main source of the MR/ LR differences.
Alternatively, a different path for the tropospheric teleconnection is a reasonable hypothesis to explain the distinct surface response, since the propagation of large-scale stationary perturbations such as the main ENSO wave train is modulated by the atmospheric mean flow.Indeed, while both models show a realistic upper-level climatological zonal wind ( U200; Fig. 3a, b) compared to observations (e.g.Mezzina et al. 2020), they do exhibit some significant differences: in LR, the North Atlantic extra-tropical jet is more zonal and the subtropical jet is weaker, while the North Pacific jet is comparable, particularly over the western part of the basin (Fig. 3c).These differences in the zonal wind may alter the propagation of the tropospheric ENSO anomalies and lead to a different European response in MR and LR.In the next section, we will address this hypothesis with the ray tracing approach described in Sect.2.2, considering various source regions.

Ray tracing results
For a fixed zonal wavenumber k, the meridional wavenumber l 2 computed from the actual U displays some longitudinal differences between the two models, in particular concerning the forbidden regions ( l 2 < 0 ) over the Mediterranean (Fig. 5, top, for k = 3 ; see Fig. S1 for k = 2, 4 ).In contrast, maps of l 2 obtained from  the zonally-averaged U are rather similar (Fig. 5, bot- tom, and Fig. S2); hence, an approach that considers the longitudinally-varying flow is essential in this case.Using the "hybrid" ray tracing method, we now examine the pathways of large scale waves, with k = 2, 3, 4 .These wavenumbers are chosen as they have been associated with tropical-extratropical teleconnections (e.g.Hoskins and Ambrizzi 1993), and correspond to the approximate longitudinal scale of the upper-level geopotential height anomalies.The rays are launched from three key regions selected according to the hypotheses from previous studies mentioned in the Introduction and in agreement with anomalous Rossby wave sources diagonsed in the models (e.g.Sardeshmukh and Hoskins 1988;Qin and Robinson 1993;Mezzina et al. 2021a, not shown).

Tropical Pacific (TP)
The main ENSO wave train is known to propagate towards higher latitudes with an eastward-arching route, reaching the western North Atlantic.Here, we consider the hypothesis of a split over North America, with part of the wave energy diverted to Europe rather than bending towards the North Atlantic (García-Serrano et al. 2011).Several source points in the central tropical Pacific have been selected (Fig. 6, top), none of which located west of the Date Line, as a region of negative l 2 in the western North Pacific would inhibit propagation after a few integration steps (see also Li et al. 2020).All rays are reflected at high latitudes, according to the regional waveguide ( Hoskins andAmbrizzi 1993, Shaman andTziperman 2005) but at different longitudes and with different angles depending on the wavenumber, so that longer waves ( k = 2 , light green) roughly follow the arching structure of the main wave train, crossing Canada and reaching the western North Atlantic, while shorter waves ( k = 4 , dark green) bend at more eastern longitudes, around Greenland, and seem to reach Europe, with k = 3 (green) display- ing an intermediate route.In principle, the main wave train should contain contributions from all these wavelengths, and the fact that the trajectories are rather similar over the North Pacific/western North America but are well separated in the North Atlantic supports the possibility that a split such as that described by García-Serrano et al. (2011), who detected a bifurcation over eastern North America, could occur.This mechanism could be at play for both ENSO-related patterns (Fig. 6), since there are no major differences between the paths in the two models.However, the rays that are actually reaching Europe are those with small wavenumber ( k = 4 ), which are close in spatial scale to the nega- tive center of action found in MR at upper levels (Fig. 6a, shading).In contrast, the positive anomaly in LR appears to be too large for such wavenumber (Fig. 6b, shading), suggesting that the split mechanism may be more relevant for the semi-isolated pattern in MR rather than the wave-like pattern in LR.

North America/North Atlantic (NANA)
A series of launching points over eastern North America has been chosen to investigate the hypothesis of a secondary wave source over the Gulf of Mexico-Caribbean Sea region emerging from inter-basin effects, as proposed by Ayarzagüena et al. (2018) and Hardiman et al. (2019) (Fig. 6, middle).Rays launched from south of 30 • N would propagate almost zonally due to the small values of l 2 (see Fig. S1) and end up trapped along the subtropical jet, hence the points have been distributed between 30 • N and 50 • N. Longer waves ( k = 2 , light green) propagate rather meridi- onally and turn equatorward at high latitudes, where they merge with shorter ones ( k = 4 , dark green) after reflection around 60 • N.Although the end of the trajectories partially matches the location of the European anticyclone in LR, its amplitude is comparable to that of the high-latitude anomaly over Fig. 6d,shading).Such a relatively strong anomaly seems more likely to correspond to a wave train emitted by a closer source.
In contrast, it cannot be excluded that waves of intermediate size ( k = 3 , green) may be related to the negative anomaly appearing over north-eastern Europe in MR (Fig. 6c).An alternative conjecture to that of a secondary wave train is that inhomogeneities in the basic state cause the main wave train to propagate through a wavier trajectory and experience a reflection (Branstator 1983(Branstator , 1985;;Trenberth et al. 1998) that is not captured by our simple approach.

Tropical Atlantic (TA)
The secondary wave source suggested by Toniazzo and Scaife (2006) would be located in the tropical Atlantic.While they considered a source close to the Equator, we have chosen starting points at slightly higher latitudes, above the zonal barrier related to the equatorial easterlies (Fig. 6, bottom).The difference in the initial slope of the ray trajectories is less striking than in the previous cases, as all rays start with a quite meridional orientation.In both models, shorter waves ( k = 4 , dark green) are trapped in the subtropical jet, while longer ones ( k = 2 , light green) keep propagating towards high latitudes.Interestingly, intermediate waves ( k = 3 , green) are also trapped in the subtrop- ics in MR (Fig. 6e), but not in LR (Fig. 6f), where, thanks to a discontinuity in the forbidden region, they propagate northwards.This difference in the propagation of k = 3 is consistent with the center of action of positive sign over the Mediterranean that dominates in LR, whose relatively strong amplitude also agrees well with a source in the TA.This positive anomaly may cancel out a possible negative signal over Europe, thus explaining the difference in the upperlevel and surface patterns of the models.

Discussion and conclusions
We have examined outputs from two models, MPI-ESM MR and LR, which represent examples where the surface ENSO signal over Europe is a semi-isolated and a wave-like pattern, respectively.
Based on our ray-tracing diagnostic, we suggest that the wave-like pattern is related to an upper-level anticyclone (for El Niño) over southern Europe belonging to a secondary wave train originating in the tropical Atlantic, which only reaches the region if the mean flow allows for the propagation of intermediate-scale ( k = 3 ) waves.This is consistent with the results of Toniazzo and Scaife (2006), who identified a wave source in the tropical Atlantic linked to a wavelike response but did not specify the scale of the waves, nor did they provide further details on their ray tracing.
The semi-isolated pattern, in turn, seems to be linked to an upper-level negative anomaly (for El Niño) over northern Europe emerging from the split and/or reflection of the main wave train.This upper-level circulation associated with the semi-isolated pattern was already observed in very early studies, such as Blackmon et al. (1983), who forced an atmospheric model with perpetual El Niño January conditions.Branstator (1985) compared Blackmon et al.'s pattern to forced steady solutions of a linear barotropic model with different mean states and found that it was better reproduced when the actual, zonally-asymmetric background flow was used: examining the distribution of normalized wavenumber vectors, he noticed a more meandering path in the energy propagation compared to the zonally-symmetric case, with a series of troughs and ridges crossing mid-latitudes and ultimately reaching Europe.That inhomogeneities in the background flow could be relevant for the propagation of the ENSO wave train to Europe is a hypothesis that was later reported in Trenberth et al. (1998)'s review (see their Fig. 8).We suggest that other regions previously suggested to play a role in the North Atlantic, such as the Gulf of Mexico-Caribbean Sea region (e.g.Ayarzagüena et al. 2018;Hardiman et al. 2019), may be interpreted in terms of this Branstator-type reflection, rather than being considered as sources of an additional wave train.Supporting this interpretation, note that the center of action over the western North Atlantic is usually weaker than that at subpolar latitudes, and its extension into North America at upper levels is upstream of these suggested wave sources.Our results are not conclusive regarding the Branstator-type reflection, since linear ray tracing does not identify a reflecting region, and further investigation is needed, possibly with the inclusion of non-linear effects.
Concerning the split hypothesis, it had been already noticed by Hoskins and Karoly (1981) that rays emanating from a source would follow different paths depending on the wavenumber and, as reported by Branstator (1983), "this dependence of propagation characteristics on wavenumber is manifested by a split in the energy as it disperses from the source".García-Serrano et al. (2011) returned to this idea in the context of the ENSO teleconnection to Europe in an effort to understand the semi-isolated pattern, since they had noticed, in their reanalysis dataset (ERA40), the presence of the corresponding (negative) upper-level anomaly.They tested the hypothesis by computing the wave activity flux, as in Karoly et al. (1989), and observed two branches separating over eastern North America: one bending equatorward as the wave train tail, the other proceeding towards the North Atlantic and turning equatorwards later, around eastern Europe.Our results support this hypothesis by indicating that the first part of the main wave train path, from the tropical Pacific to Canada, is similar for all the wavenumbers considered, but a split may occur later, with shorter waves (i.e.larger wavenumbers) indeed reaching Europe.The split is potentially present in both models discussed here, but the approximate k = 4 scale of the semi-isolated pat- tern over Europe in MR suggests that this pathway is more likely related to this type of response, as García-Serrano et al. (2011) speculated, rather than to the wave-like pattern.
Our modeling results are consistent with those obtained using the NOAA-20CR reanalysis over 1901-2014.The ENSO surface response in this dataset resembles the semiisolated pattern (Fig. 1b for DJF; see Mezzina et al. 2020 for JFM), and the upper-level response is the expected samesigned anomaly over northern Europe (Fig. 7, shading).The ray trajectories are similar to those found in the models (see also Fig. S3 for l 2 ) and again support the split hypoth- esis (Fig. 7a) and endorse the possibility of Branstatortype reflection of the main wave train over North America (Fig. 7b).However, mixed results are found when considering sources in the tropical Atlantic, although the overall trajectories are similar to the ones in the models.All rays with k = 2 (light green) are propagating towards Europe, and all those with k = 4 (dark green) are reflected equator- ward, but for k = 3 (green) both options occur, with the rays equally distributed (3 out of 6 are reflected).These results suggest that in the real world the wave source in the tropical Atlantic is effective in triggering a secondary wave train, but the amount of energy reaching Europe and its competition with the main wave train determine the observed signal.
We finally stress that the "hybrid" ray tracing approach used here is not strictly valid and caution is advised in the interpretation of the results.Despite that, we believe that it may constitute a valid tool for further applications beyond ENSO teleconnections, alternative or complementary to other diagnostics (e.g.Branstator 1985;García-Serrano et al. 2011;Benassi et al. 2021).

Fig. 2
Fig. 2 Linear regression onto the N3.4 index of 200-hPa geopotential height (top) and sea-level pressure (bottom) anomalies in the two model's versions: MR (left) and LR (right).JFM.Contours indicate statistically significant areas at the 95% confidence level bottom, shading), though with slightly different amplitude.Likewise, the anomalous tropical upper-level divergence, represented by the 200-hPa velocity potential, which is realistic compared to observations (e.g.García-Serrano et al. 2017), is only slightly stronger in LR (Fig. 3 bottom, contours) and thus unlikely the cause for the different signal over Europe.

Fig. 3
Fig. 3 Top: 200-hPa zonal wind climatology in a MR b LR c MR-LR.Contours indicate statistically significant areas at the 95% confidence level.Bottom: linear regression onto N3.4 of SST (shad-

Fig. 6
Fig. 6 Ray paths for k = 2 (light green), 3 (green), 4 (dark green) in MR (left) and LR (right) from different source regions: TP (top), NANA (middle), TA (bottom).Central trajectories are plotted with thick arrows to indicate the mean propagating path of the envelope.Filled contours show the linear regression of 200-hPa geopotential

Fig. 7
Fig. 7As in Fig.6, but for NOAA-20CR.In panel (c), the bottom line of sources in Fig.6e, f is discarded due to the higher latitudes of the equatorial zonal barrier, and the thick arrow is omitted for k = 3