Atmospheric tropical modes are important drivers of Sahelian springtime heatwaves

Heatwaves pose a serious threat to human health worldwide but remain poorly documented over Africa. This study uses mainly the ERA5 dataset to investigate their large-scale drivers over the Sahel region during boreal spring, with a focus on the role of tropical modes of variability including the Madden–Julian Oscillation (MJO) and the equatorial Rossby and Kelvin waves. Heatwaves were defined from daily minimum and maximum temperatures using a methodology that retains only intraseasonal scale events of large spatial extent. The results show that tropical modes have a large influence on the occurrence of Sahelian heatwaves, and, to a lesser extent, on their intensity. Depending on their convective phase, they can either increase or inhibit heatwave occurrence, with the MJO being the most important of the investigated drivers. A certain sensitivity to the geographic location and the diurnal cycle is observed, with nighttime heatwaves more impacted by the modes over the eastern Sahel and daytime heatwaves more affected over the western Sahel. The examination of the physical mechanisms shows that the modulation is made possible through the perturbation of regional circulation. Tropical modes thus exert a control on moisture and the subsequent longwave radiation, as well as on the advection of hot air. A detailed case study of a major event, which took place in April 2003, further supports these findings. Given the potential predictability offered by tropical modes at the intraseasonal scale, this study has key implications for heatwave risk management in the Sahel.


Introduction
The Sahel is a tropical semi-arid region located in West Africa that experiences high temperatures during the largest part of the year (Nicholson 2018), and especially in spring (March-June, MAMJ). Recent studies (Fontaine et al. 2013;Ringard et al. 2016;Moron et al. 2016;Oueslati et al. 2017;Guichard et al. 2017; Barbier et al. 2018) have highlighted an increase of extreme temperatures, and future projections predict even worse heatwave conditions (Russo et al. 2016;Dosio 2017;Déqué et al. 2017;Xu et al. 2020;Raymond et al. 2020).
Unfortunately, this region is also one of the least economically developed in the world (e.g. Davidson et al. 2003;Tschakert 2007), implying high levels of vulnerability. An efficient way of using the limited resources available to mitigate heatwave impacts, is to focus on preventive actions, based on skilful forecasts of these hazards as already implemented in other places (e.g. WMO N°1142; Wilkinson et al. 2018). The prevention however requires knowledge and understanding of the large-scale drivers of Sahelian heatwaves, which are both currently under-documented.
Greenhouse effect of moisture, hot air advection and incoming solar radiation (at daytime) appear as the dominant processes inducing extremely hot temperatures in spring over the Sahel (e.g. Slingo et al. 2009;Guichard et al. 2009;Oueslati et al. 2017;Guigma et al. 2020). The spring season also marks the peak of monsoon convection over the Guinean region of West Africa (i.e. south of 10° N; Nguyen et al. 2011), and, although not intuitive, there are potential links between convection and Sahelian 1 3 heatwaves. Guigma et al. (2020) indeed found that largescale heatwaves over the eastern Sahel are associated with a low-level cyclonic circulation anomaly over northern tropical Africa with a decrease of precipitation over the Guinean band, versus an anticyclonic circulation anomaly and an increase of Guinean precipitation for heatwaves over the western Sahel. When convection is weakened over the Guinean band, the associated northeasterly flow indeed allows the transport of hot air from the eastern Sahel/Sahara towards the coast of the western Sahel; conversely, a convective intensification conveys moisture from the Gulf of Guinea into the Sahel, through northward penetration of monsoon flow (e.g. Lothon et al. 2008;Couvreux et al. 2010), leading to longwave warming.
At the synoptic and intraseasonal scale, the Guinean convection is strongly modulated by tropical modes of variability, namely the Madden Julian Oscillation (MJO) and equatorial waves (e.g. Gu 2009; Kamsu-Tamo et al. 2014;Berhane et al. 2015). Can these modes of convective variability by extension modulate heatwaves in the Sahel? This is a plausible scenario in view of previous studies over West Africa. Moron et al. (2018a) for instance, using a weather type approach, found a significant modulation of near surface temperature over northern tropical Africa by the MJO and equatorial Kelvin (EK) waves at the intraseasonal scale in spring. Advection of both hot air and moisture plays an important role in this modulation. Kalapureddy et al. (2010) also found that the synoptic variability over the Sahel, during the pre-(April-June) and post-monsoon (October and November) seasons, is controlled by monsoon surges with a periodicity similar to that of the African Easterly Waves (AEWs; 3-5 days). Furthermore, Couvreux et al. (2010) indicated that these surges can be stationary or westward moving, along with the AEWs. In the same logic, Mera et al. (2014) showed that synoptic and subseasonal circulation disturbances, driven by extratropical cyclones and also by EK and equatorial Rossby (ER) waves, and possibly by the MJO, lead to influx of moisture from the Gulf of Guinea into the Sahel.
Connections between heatwaves and atmospheric synoptic and intraseasonal modes of convective variability have already been evidenced in other regions. Murari et al. (2016), for example, found that the delay of the Indian monsoon onset weakens southwesterlies in the Arabian Sea, favouring clear sky days over India, and subsequently longer lasting and warmer heatwaves. In South America, the intraseasonal variability considerably modulates heatwaves in association with the Southern Atlantic Convergence Zone (Cerne and Vera 2011). In the southeastern part of Australia, Parker et al. (2014) showed that warm conditions are significantly associated with the phases 3-6 of the MJO. More recently, the summer 2018 long lasting heatwave over Northeast Asia was found by Hsu et al. (2020) to have been favoured by an unusually strong MJO over the western Pacific warm pool.
The present research aims at assessing the extent to which and in what ways Sahelian heatwaves are modulated by tropical modes during the MAMJ season. The reason for the focus on tropical modes is the potential for predictability at intraseasonal lead times. Tropical modes are indeed good sources of predictability at these scales (Moron et al. 2018b;Dias et al. 2018;Bengtsson et al. 2019;Li and Stechmann 2020;Judt 2020).
To the best of the authors' knowledge, this study is the first to use a systematic and comprehensive approach to investigate the impact of tropical modes on actual heatwaves in the Sahel. Two complementary methods are used for this purpose. The first consists of a statistical study over the 1979-2018 period while the second exemplifies the emerging overall processes with a detailed case study event in April 2003.
The rest of the manuscript is structured as follows. The data and the methods used are summarised in Sect. 2. Then, Sect. 3 presents the results obtained from the statistical study. In Sect. 4, the case study of the 2003 heatwave event is analysed in detail. Finally, conclusions and perspectives are given in Sect. 5.

Methodology
The Sahel area is hereafter defined as the continental domain limited by the coordinates 20° W, 30° E, 10° N, 20° N. Since this work aims at depicting the large-scale drivers, the analysis however extends to the entire northern half of Africa, and sometimes to an even larger domain.

Data
The main dataset used in this paper is the fifth generation of the European reanalyses ERA5 (Hersbach et al. 2020) at a resolution of 0.5° × 0.5°, covering the period 1979-2018. The reason for the choice of this dataset is that it likely offers the best representation of near surface meteorological variables (Olauson 2018;Ramon et al. 2019) as well as radiative fields over land (Martens et al. 2020), and provides a comprehensive and self-consistent set of variables for diagnostic analysis.
This dataset is thus used in this research for (i) heatwave detection from daily minimum and maximum temperature at 2 m (hereafter referred to as "Tmin" and "Tmax" respectively), (ii) tropical mode filtering from outgoing longwave radiation (OLR), and (iii) retrieving physical fields to understand the mechanisms of modulation. Previous studies (Oueslati et al. 2017;Barbier et al. 2018) have already shown the good quality of near surface thermal indices (including Tmin and Tmax) over the Sahel in ERA-Interim of which ERA5 is an improvement. OLR in ERA5 is taken as the negative of the net top-of-atmosphere (TOA) thermal radiation. Wang et al. (2017), Tall et al. (2019), Wright et al. (2020) and Hersbach et al. (2020) assessed radiative fluxes in different products and found that ERA5 shows TOA fluxes (including OLR) that are very consistent with the observed. Using observed OLR data such as the daily interpolated data from the NOAA (Liebmann and Smith 1996) does not significantly alter the results obtained with ERA5 (not shown), making it suitable for the present study.
Some of the physical variables used for analysing the mechanisms of heatwave modulation are directly available from ERA5: temperature, zonal and meridional components of wind, specific humidity, precipitable water and radiation data. As in Oueslati et al. (2017), the surface energy budget is derived from the radiation fields using the following equation: where SWR is the net shortwave radiation, LWR is the net longwave radiation, SHF is the sensible heat flux, LHF is the latent heat flux and c s is the surface heat capacity.
Radiative fluxes are hereafter counted positively when directed from the atmosphere to the surface.
Other meteorological variables necessary for understanding the physical mechanisms of heatwave modulation are not directly accessible from ERA5 and are thus derived. These include advection of heat (Eq. 2) and of specific humidity (Eq. 3) which are obtained from the horizontal components of wind speed, temperature and specific humidity: where T is temperature, q specific humidity and �� ⃗ V horizontal wind which is decomposed into its zonal (u) and meridional (v) components.
All variables of the ERA5 database described in this section are available at an hourly frequency. The radiative variables are extracted as 12-hourly accumulations (0600 UTC to 1800 UTC for daytime and 1800 UTC to 0600 UTC of day + 1 for nighttime) while the others as instantaneous entries at a 03-hourly resolution. For the latter, the daytime value is obtained by averaging over 0900, 1200, 1500 and 1800 UTC and the nighttime value over 2100, 0000, 0300 and 0600 UTC.
For the specific case of the derived variables (heat and moisture advection), they are first computed at the 03-hourly timescale before daily averaging.
In addition to the ERA5 reanalysis, various observational data have been utilised mainly for the case study (Sect. 4). Satellite estimates of radiative surface fields (incoming and outgoing shortwave and longwave radiation) were extracted from the Clouds and Earth's Radiant Energy System (CERES) Synoptic (SYN1deg) Product (Wielicki et al. 1996;Doelling et al. 2013Doelling et al. , 2016 at a spatial resolution of 1° × 1° over the time period 2000-2018. The initial 03-hourly temporal resolution was aggregated to a daily frequency. Satellite estimates of precipitable water from the Atmospheric Infrared Sounder (AIRS, Teixeira et al. 2013) aboard NASA's second Earth Observing System polar-orbiting platform was also used for the case study at 1° × 1° spatial grid. Twice-daily data provided by the ascending and descending orbits were averaged to get daily estimates. This dataset runs from 2002 to 2016. To analyse the evolution of precipitation over the case study period, daily rainfall totals from the Global Precipitation Climatology Project (GPCP, Schamm et al. 2014) were used at a resolution of 1° × 1°. Finally observed 2 m temperature and specific humidity data at the Demokeya station in Sudan located at 30.5° E-13.3° N (Ardö 2013) were also investigated. The data cover the 2002-2012 period with a temporal resolution of 30 minutes. This provides a valuable in-situ observation over the eastern Sahel, a data sparse region, where the case study is mainly investigated.

Heatwavedetection
Among the different definitions of heatwaves (Perkins 2015), a choice has been made on one allowing to capture major, large-scale events occurring over synoptic to intraseasonal scales during the spring season. Tropical modes are indeed planetary-scale disturbances and their effects are more likely to cover extensive areas than smaller isolated areas. Therefore, as in Guigma et al. (2020), Tmin and Tmax data are first 90-day highpass filtered to retain variability on the scales of interest. Then, potential heatwaves are defined at a given grid-point as spells of at least three consecutive days where the daily anomalies of Tmin or Tmax (taken separately) exceed the 90th percentile of the distribution of the corresponding calendar days. The 90th percentile for a given calendar day is derived, similarly to Russo et al. (2014) and Guigma et al. (2020), over a 31-day window centred on that day, to remove noise caused by the relative shortness of the dataset. After this step, to ensure sampling of only absolutely hot events regardless of the period of the year, potential heatwaves sampled on days where the raw data (i.e. data before deriving the anomalies) do not exceed the 75th percentile of their total distribution are discarded (Guigma et al. 2020). Heatwave intensity is defined, at each grid-cell and for each heatwave day (not event), as the 90-day highpass filtered anomaly of Tmin (for a nighttime heatwave) or Tmax (for a daytime heatwave) on that day. Next, major large-scale events are defined by adding a minimum spatial extent threshold of 600,000 km 2 using a region growing technique as in Barbier et al. (2018). The outcome is a set of daily binary masks specifying the regions where these large-scale heatwaves occur, providing the final sample of heatwaves analysed in this study.
It should be noted that heatwaves can also be defined from other thermal indices than Tmin and Tmax, taking into account additional environmental variables (moisture, wind, solar radiation) as in Guigma et al. (2020). Only Tmin and Tmax are however retained for this study because of the consistency of the results across these indices (not shown).

Tropical mode detection
In this study, the expression "tropical modes'' refers to the synoptic and intraseasonal modes of variability of the tropical atmosphere including the MJO and equatorially trapped waves. The MJO is the most important driver of intraseasonal variability of the tropical atmosphere (Madden and Julian 1971, 1972, 1994 and influences the West African Monsoon (e.g. Sossa et al. 2017). In this study, equatorial waves are explored in connection with convection (hence the expression convectively coupled equatorial waves CCEWs). Over West Africa, Kamsu-Tamo et al. (2014) and Schlueter et al. (2019a, b) demonstrated statistically and dynamically that they have an impact on convection. In the present paper, the equatorial waves under investigation are the ER and EK waves. Other waves such as AEWs and mixed Rossby-gravity waves are not included as they are not so active in the region during the season of interest (MAMJ).
Each tropical mode is characterised by a range of wavenumbers, periods and equivalent depths summarised in Table 1. A wavenumber is the number of waves (each wave comprising one ridge and one trough) that it takes to zonally circumscribe the globe. The period of a wave measures the time separating the passage of two consecutive ridges at a given location. As for the equivalent depth, a more complex notion used for ER and EK wave filtering, it refers to the depth of the shallow layer of atmospheric fluid that is required (by theory) to get the appropriate values for the time-varying and horizontal components of their motion (Wheeler and Nguyen 2015).
Tropical modes are detected from the fields of OLR, which is a good proxy for tropical convection (Arkin and Ardanuy 1989), and has been extensively used in tropical mode investigation over West Africa (e.g. Lavender and Matthews 2009;Pearson et al. 2010;Yang et al. 2018). To extract a given mode, the mean and first three harmonics of the OLR annual cycle are first removed to obtain intraseasonal anomalies. Then, the filtering in the wavenumberfrequency space consists of setting to zero any spectral coefficient outside the window corresponding to the mode of interest (Table 1), such that only the relevant coefficients are retained. The "kf_filter" function of the NCAR Command Language (NCL) was used for this operation. The output of the wavenumber-frequency decomposition is thus a three dimensional (time, longitude and latitude) filtered OLR field for each mode. The full explanation of the method can be found in Wheeler and Kiladis (1999) and Kiladis et al. (2006). No decomposition of the input data into asymmetric and symmetric components is done because the area of interest is fully located on one side of the Equator (northern hemisphere), as already done in several previous studies (e.g. Lafore et al. 2017).
In order to characterise the daily local activity of each mode, a phase and amplitude are determined as follows. Firstly, the 5° E longitude (centre of the Sahel domain) is arbitrarily set as reference longitude. Minor variations of the reference longitude do not significantly impact the results (not shown). Then, the daily mode-filtered OLR data at this longitude are averaged over the Guinean zone, i.e. for the latitudes between the Equator and 10°N (vertical line in Fig. 1), where convection is the most active in spring (Nguyen et al. 2011). Including latitudes up to 20° N does not make significant changes to the results. As in Riley et al. (2011), the resulting unidimensional timeseries and its (first order) local time derivative (both normalised by their respective standard deviations) are used to derive the trigonometric form of the corresponding mode, from which an angle and an amplitude are extracted for each day. On a given day, a mode is considered active only if its amplitude is greater than one.  Mekonnen et al. 2008) Changes to this threshold do not affect the spatial structures of the impact of tropical modes on heatwaves but the magnitude of this impact is slightly modified and is discussed in Sect. 3.3. For the active days of each mode, the angles are further binned into eight 45° wide phases labelled 1-8 such that phase 1 represents instances where the mode is the most convectively suppressed (maximum of OLR), phase 5 the most convectively enhanced (minimum of OLR) and the six others are transitory phases [see the description of the full method in Riley et al. (2011)]. Thus, a given active mode (active because passing the minimum amplitude threshold) can be either in a convectively suppressed or a convectively enhanced phase. It should be noted that, while the practice of setting a reference longitude effectively allows to detect the passage of tropical modes in the region, and to better characterise their local activity, it however comes with the caveat that, away from this longitude, the signal weakens progressively. As a consequence, towards the boundaries of the region, the modes with large wavenumbers are not as well accounted for as those with small wavenumbers.

Modulation of heatwaves by tropical modes and associated evolution of physical variables
The modulation of heatwave probability of occurrence and heatwave intensity is assessed, during the spring season, by modifying a method initially used by Xavier et al. (2014) for extreme rainfall and Hsu et al. (2017) for heatwaves. It consists of comparing heatwave probability/intensity when a given mode is active on a specific phase versus the climatological probability/intensity during the MAMJ season. The modulation of heatwave probability (MP) is thus given by the following formula: where P x is heatwave probability given the mode is in phase x, i.e. the number of days where heatwaves occur during phase x over the total number of days where the mode is in this phase; P a is the climatological probability of heatwaves, i.e. the number of heatwave days (regardless of the activity of the modes) over the total number of days. The probabilities are defined in terms of days rather than events, consistently with Hsu et al. (2017). Likewise, the modulation of heatwave intensity (MI) is defined as follows: where I x is the average heatwave intensity (defined for each heatwave day and at each grid-cell as the highpass filtered anomaly of Tmin or Tmax on that day; Sect. 2.2) given the mode is in phase x and I a the climatological intensity of heatwaves.
As will be shown in Sect. 3.1, heatwaves can be associated with distinct active tropical modes at different phases, in such a way that the superposition of two modes can lead to an amplified or reduced forcing on heatwaves. The modulation of heatwave probability by two superposed modes is investigated using the following formula: P xy : heatwave probability when the first mode is in phase x and the second mode in phase y. P a : heatwave climatological probability.
The significance of the modulation patterns is tested at each grid-point and for each phase of the modes through bootstrap resampling (Wang et al. 2008;Mazdiyasni and AghaKouchak 2015;Nissan et al. 2017;Harrington et al. 2019). Let N be the number of days where a given mode is in phase × (or one mode in phase × and another in phase y for the superposition case). Then, 1000 random samples, each of size N days are generated. Heatwave probability (or average intensity) is calculated for each sample. From these random probabilities (or average intensities), corresponding values of modulation of heatwave probability (intensity) are derived using Eqs. (4), (5) or (6) accordingly. Significance is finally tested at a 5% probability level against these 1000 random modulation values. The results of the modulation of heatwave probability are presented in Sect. 3.3, and that of heatwave intensity in Sect. 3.4.
To understand the physical mechanisms associated with the modulation of heatwaves, the composites of each of the diagnostic variables (presented in Sect. 2.1.) are derived over all active mode-phase days. To do this, the mean and the first three harmonics of the seasonal cycle for each variable are first removed. Then, they are passed to a Lanczos highpass filter (Duchon 1979) to retain variability at timescales shorter than 90 days. The filtering causes a loss of data at both edges of the input time series (also observed for heatwave detection and OLR filtering); therefore the analyses cover the 1980-2017 period. Finally, the composite mean of a given variable over a certain instance (e.g. phase × of a mode) corresponds to the average of this variable over all days of this instance. At each grid-point, the significant departure of the average values from zero is tested using the Student'st test at the 0.05 probability level. The results are shown in Sect. 3.5.
For the case study purpose, the activity of tropical modes is projected onto various physical fields, similarly to previous studies (e.g. Kiladis et al. 2006;Knippertz and Todd 2010). To do this, the daily mode-filtered OLR data are regressed onto the daily highpass filtered anomalies of the physical variables at each grid-point over the 1980-2017 period. Then, the influence of each tropical mode on the variable during the heatwave event is qualitatively estimated on a daily basis, using the predicted value of the diagnostic variable by the linear regression model. The advantage of this technique is that it allows on a given day to know the likely impact of a certain mode on a physical process (results in Sect. 4.4). Table 2 presents statistics associated with the activity of tropical modes over the Guinean region as detected at the reference longitude of 5° E (described in Sect. 2.3) focusing mainly on the most suppressed and the most enhanced phases (phases 1 and 5 respectively). During the MAMJ season, at least one mode is active 96% of the time irrespectively of the phase. Therefore the spring season is intense in terms of tropical mode activity. For each of these modes, the activity is evenly distributed across the eight phases without a preference for any of them (not shown). The number of active days is not much different from one mode to another but it is apparent that the MJO has the largest sample. The occurrences of exactly one mode (whichever it is) at a given time exclusive of any other mode are quite infrequent (19% of the time). For a given active phase of any mode, merely 10% of its passages satisfy this configuration. Therefore active tropical modes generally overlap. Consequently, an explicit analysis of the impact of these situations on Sahelian heatwaves is undertaken in Sect. 3.3. The amount of OLR variability explained by the modes is clearly different between them. Figure 1 shows the variance of the mode-filtered anomalies of OLR during the MAMJ season. The MJO, EK and ER waves by decreasing order, have an intense activity over the Guinean sector in boreal spring. For the MJO and EK waves, the variance peaks at above 200 W 2 m − 4 near the Equator (Fig. 1a, c). The MJO is active over the entire sector between 10°S and 15°N. On the other hand, the EK waves are active only within a band surrounding the Equator consistent with theory (the high variance of EK waves observed at the edges of the domain are mainly due to imperfections in the filtering technique which allows in noise from the mid-latitude eastward Rossby waves). As for the ER waves, they maintain a relatively important activity (peaks exceeding 150 W 2 m − 4 ) over the eastern and central parts of the domain ( Fig. 1.b). There is a slightly different seasonality between the modes during the spring season. The peak of convective activity related to the MJO occurs in March and April, that of EK waves in April and May whereas the ER is the most active in May (Fig. SM1). The activity of these three modes decreases significantly in June, prior to the onset of the Sahelian phase of the West African monsoon (e.g. Fitzpatrick et al. 2015). This seasonality is further shown by the monthly count of active days of each mode at the reference longitude of 5° E (Fig. 2).

Synoptic to intraseasonal heatwaves in the Sahel
Previous publications have already elaborated on the statistical characteristics of Sahelian heatwaves at different timescales. There is a general consensus on their short-lasting nature (e.g. Oueslati et al. 2017;Guigma et al. 2020) and their important intensity owing to the hot mean state of the region (red lines in Fig. 2). They are more frequent over the eastern and central Sahel, as shown by the count of springtime heatwave days at each grid-point over the 1980-2017 period in Fig. SM2 (which provides the sampling basis of heatwaves for this study). As with tropical modes, the frequency of occurrence of synoptic to intraseasonal scale heatwaves has a moderate seasonality. Figure 2 shows that daytime heatwaves are mostly observed between February and June, with a peak in April. On the other hand, nighttime heatwaves occur mainly between March and July with May standing as the most heatwave-prone month. There is thus a shift of nighttime heatwaves with respect to their daytime counterparts. The main reason for this is probably the larger dependence of Tmin (and thus nighttime heatwaves) on lowlevel moisture whose penetrations into the Sahel are more frequent closer to the onset of the monsoon Mera et al. 2014).

Modulation of heatwave probability of occurrence by tropical modes
This section first discusses the impact that each tropical mode, taken separately, has on heatwave occurrence, before analysing the outcome of the superposition of several modes.

Heatwave modulation by modes taken separately
Figure3 shows the phase-longitude diagram of the modulation of heatwave probability (MP) averaged over the Sahel band (10° N-20° N) for the MJO, ER and EK waves during the spring season. Among these three investigated modes, the MJO stands as the lead modulator of heatwave occurrence in the Sahel (Fig. 3a, b) and its eastward propagation is clearly reflected in the modulation MP pattern. Overall, phases 1 to 4 are associated with higher heatwave probability than the climatological probability of occurrence. The values of MP, averaged over the Sahel band, can be as high as 1 (i.e. heatwaves are twice more likely under these phases than according to their climatology). Reversely for phases 5-8, the probability of heatwave occurrence is lower than usual. For these phases, MP is typically below − 0.5, meaning that the probability of occurrence is less than half of its climatological value. The patterns are relatively similar for daytime and nighttime heatwaves, but it is apparent that the stripe of positive MP values is wider in Tmin heatwaves than in Tmax (red stripe in Fig. 3a vs. b), and reversely for the stripe of negative MP values (blue stripe in Fig. 3a vs. b), especially over the western Sahel. The modulation of Sahelian heatwaves by the ER wave (Fig. 3c, d) is less important than by the MJO. Heatwave probability of occurrence is mostly below normal when the ER wave is in phases 5 to 8 with values of MP generally higher than − 0.5 (i.e. of lower intensity than observed for the MJO). In phases 1 to 4, there is a higher heatwave risk than normal over much of the Sahel, with positive MP values although usually lower than 1. It should be noted that the modulation by the ER wave is more sensitive to the diurnal cycle and geographic location than the two other modes. For example, daytime heatwaves over the eastern Sahel are only marginally affected by the activity of the ER wave whereas nighttime heatwaves are much more impacted. On the other hand, over the western Sahel, the modulation is more important for daytime than nighttime heatwaves (this is valid for both positive and negative modulations).
The spatial patterns of the modulation of heatwaves by the EK wave (Fig. 3e, f) are similar to that of the MJO with marked eastward propagation in both cases. However, they present narrower stripes and lower magnitude than with the MJO. The values of MP are indeed in most cases absolutely below 0.5, making the EK wave the least important of the three modes for heatwave occurrence in the Sahel. They present only little sensitivity to the diurnal cycle.
The amplitude of tropical modes, as stated in Sect. 2.3, does not affect the spatial distribution of MP. On the other hand, with increasing amplitudes, the values of MP are slightly increased (without changing their sign). This is illustrated in Fig. SM3, using a minimum amplitude of 2. Therefore, the stronger the amplitude of tropical modes, Fig. 3 Phase-longitude diagrams of the modulation of (left panel) nighttime and (right panel) daytime large-scale heatwave probability of occurrence by the MJO, ER and EK waves during MAMJ averaged over the Sahel band (10° N-20° N). The modulation corresponds to the variable MP defined in Sect. 2.4 (Eq. (4)) as a comparison of heatwave probability of occurrence in a given active phase versus the climatological probability. Significance was tested at the probability 0.05 level, using bootstrap resampling of tropical mode active days, with 1000 repetitions the more confident the changes of heatwave probability described above are.
The main reason why the MJO has a greater influence on heatwaves than the ER and EK waves is very likely related to its spectral properties. Heatwaves are slowly varying events whose minimum duration is set to three days here. Therefore, for a given mode to influence their occurrence, it should be able to sustain conditions which are favourable to the heating of the atmosphere for a long-enough amount of time. Evidently, the longer the period of a mode the higher its ability to develop such conditions. This sensitivity of the modulation to the temporal scale corroborates previous results by Schlueter et al. (2019a) who investigated the modulation of precipitation in West Africa by tropical modes. Indeed, they found that, on the scale of three days, the MJO, ER and EK waves have an equal importance on precipitation while on the scale of a week, the effect of EK wave becomes marginal and beyond 20 days, the MJO is the only mode able to significantly contribute to precipitation variability. Besides, the additional constraint on the spatial extension of heatwaves may also favour the MJO. Consequently, the results of this study imply that the smaller the wavenumber of a mode, the more likely this mode can trigger heatwaves as it is also more able to promote warming conditions that spread over large areas.
It should be noted that the patterns of the modulation of heatwaves shown at the seasonal scale are generally consistent across the different months. However, some differences are noticeable, especially with nighttime heatwaves, for some modes. Thus, phases 1 to 4 of the ER wave lead to a decreased probability of nighttime heatwave occurrence over the central Sahel (broadly between 0° E and 20° E) in early spring (March-April) whereas the probability is increased in late spring (May-June), and reversely for phases 5 to 8 (Figs. SM4.a vs. SM4.b). For the EK wave, the decrease of heatwave probability observed in phases 5 to 8 is much more pronounced in late spring than in early spring (Figs. SM4.c vs. SM4.d).

Heatwave modulation by superposition of modes in various phases
Section 3.1 evidenced the need for special consideration of instances where multiple tropical modes are superposed. These superpositions may occur on any phase of the modes. For conciseness purpose, only the most convectively suppressed phase (phase 1) and the most convectively enhanced phase (phase 5) are investigated in the present section. However, it should be noted that, although these phases represent the extremes of the convective activity driven by the modes in the region, they do not necessarily correspond to the peak of their effects on heatwaves as shown in Fig. 3. It is therefore possible to consider other phases when assessing the effect of the superposition of multiple modes on heatwaves. Sample sizes of different combinations of modes in their phases 1 and 5 are relatively moderate, not exceeding 50 days (Table 3) but allow an insight as to how (much) these configurations impact Sahelian heatwave occurrence. The combinations can be mutually constructive (i.e. all modes on the same phase) or conflicting (e.g. one mode in an enhanced phase and another in a suppressed phase).
Figure4 presents the zonal averages over the Sahel domain (10° N-20° N) of the modulations by pairs of tropical modes. The superposition of the ER wave and the MJO gives relatively similar results for daytime and nighttime, at least over the eastern Sahel (Fig. 4a, b). When these modes co-occur in a constructive way, the resulting MP pattern is comparable to the one obtained when they occur separately and mainly consists in dipolar MP structures: for the suppressed case, an increase of heatwave probability in the western Sahel and a decrease in the eastern Sahel (red dotted lines in Fig. 4a, b), and conversely for the enhanced phase (blue solid lines in Fig. 4a, b). When the MJO and the ER wave are conflictingly superposed, there is a decrease of heatwave probability in most cases (orange dotted lines and green solid lines in Fig. 4a, b). A few exceptions to this include an increase of daytime heatwave probability over the western Sahel when an enhanced ER wave and a suppressed MJO are superposed (orange dotted line in Fig. 4a), and an increase of both daytime and nighttime heatwave probability over the eastern Sahel when the ER wave is in a suppressed phase and the MJO in an enhanced phase (green solid lines in Fig. 4a, b).
The superposition of the EK wave onto the MJO (Fig. 4c, d) leads to modulations of relatively large magnitude which generally take the sign of a singular modulation by the MJO. Similarly to the combined MJO-ER modulation, dipolar patterns are observed as a result of constructive superpositions of the MJO and the EK wave. Heatwave probability is increased over the western Sahel (especially at daytime), in opposition to a relatively small decrease over the eastern Sahel when both modes are in their suppressed phases (red dotted lines in Fig. 4c, d). When they are both in their enhanced phases, heatwave occurrence is increased over the eastern Sahel whereas over the western Sahel it is either decreased (case of daytime heatwaves; blue solid lines in Fig. 4c) or unchanged (case of nighttime heatwaves; blue solid lines in Fig. 4d). The conflicting configuration in which the MJO is in phase 1 has only minor impacts on heatwave occurrence in the Sahel (orange dotted lines in Fig 4c, d).
On the other hand, in the contingency where the MJO is in phase 5 and the EK wave in phase 1, the modulation of daytime heatwaves is very close to that observed with the MJO alone (green solid lines in Fig. 4c) whilst nighttime heatwave probability is significantly increased throughout the Sahel (MP often exceeds 2; green solid lines in Fig. 4d). Finally the superposition of the ER and EK waves generally causes only minor changes to Sahelian heatwave probability (Fig. 4e, f). However when both modes are in a suppressed phase, the probability of heatwave occurrence is significantly augmented over the western Sahel (a modulation of up to 2). In summary, the superposition of multiple modes over the African domain can substantially increase heatwave probability over the western Sahel when they are all in a convectively suppressed phase and over the eastern Sahel when they are all in a convectively enhanced phase, while the outcome for conflicting superpositions depends on the modes under consideration.

Modulation of heatwave intensity
Heatwave intensity is overall less sensible to tropical modes than the probability of occurrence. Some significant patterns are however obtained with the MJO and ER wave and are shown in Fig. 5. This figure reveals that the intensity of heatwaves is increased when the MJO and the ER waves are in phases 4 to 7, with the modulation MI often reaching 0.4. On the other hand, the intensity is decreased in phases 1 to 3 of the MJO and ER wave, and MI can also get below − 0.4. Finally, the influence of EK waves on  5)) as a comparison of heatwave average intensity in a given active phase versus the climatological intensity of heatwaves. Significance was tested at the 0.05 probability level, using bootstrap resampling of tropical mode active days, with 1000 repetitions heatwave intensity is quite minor (Fig. 5e, f). Nighttime heatwave intensity is also more affected than daytime heatwave intensity (Fig. 5a, c vs. Fig. 5b, d). From a spatial point of view, the modulation is more marked over central and the eastern Sahel than over the western Sahel.
A comparative examination of the patterns of MP (Fig. 3) versus that of MI (Fig. 5) shows that the two do not match exactly, and this can have implications for heatwave risk management, especially for nighttime heatwaves whose intensity is more impacted by tropical modes. Thus, in phases 6 and 7 of the MJO, nighttime heatwaves are less probable but are more intense (Figs. 3a vs. 5a). Conversely, they are more probable in phases 2 and 3 of the MJO but less intense. When the MJO is in phase 4, it poses a higher risk because heatwaves are both more probable and more intense. This is also the case when the ER wave is in phases 4 to 7 (Figs. 3a vs. 5a).

Potential mechanisms of modulation
Previous studies (e.g. Oueslati et al. 2017;Guigma et al. 2020) have pointed to the increase of longwave radiation at the surface due to water vapour greenhouse effect and hot air advection as processes underpinning Sahelian heatwaves. This motivates the composite analysis of thermodynamic and dynamical variables presented below for active phases of the modes. The aim here is to provide an insight as to how heatwave processes are impacted by the physical conditions created by tropical modes.
From Fig. 6a it is apparent that the patterns of the advection of heat into the Sahel, although patchy, are relatively consistent with that of the modulation of heatwave occurrence. Increased heatwave occurrence is indeed broadly associated with hotter air advection whereas cooler air advection goes with a decrease of heatwave probability (except over a narrow area around 20° E). The influence Fig. 6 Phase-longitude diagrams of the composite of daily anomalies of (left panel) heat advection in K day − 1 and (right panel) precipitable water in mm over active phases of the MJO, ER and EK waves during MAMJ averaged over the Sahel band (10° N-20° N). Significance was tested using a Student t-test at the 0.05 probability level of tropical modes on heatwave occurrence is also due to their control on water vapour and its associated longwave radiative effect. As shown in Fig. 6, over the eastern Sahel and in phases 4 and 5, the higher probability of heatwave coincides with a positive anomaly of precipitable water reaching 3 mm on zonal average (the longwave radiation and precipitable water patterns are very similar, not shown). However, the MJO-induced precipitable water anomaly has its greatest impact on heatwave intensity. Phases 4 to 7, which are associated with higher heatwave intensity, are also those for which the anomaly of precipitable water is positive; reversely, in phases 1 to 3 where heatwave intensity is lower, there is a dry anomaly especially marked over the eastern Sahel (Figs. 5a, b vs. 6b). The stronger dependence of nighttime heatwaves on water vapour (Guigma et al. 2020;Largeron et al. 2020) likely explains why their intensity is more affected by the MJO (Fig. 5a vs. b).
Associated with ER wave, an advection of hot air is present too, though mostly confined to the eastern Sahel (Fig. 6c) which could be linked to the increase of heatwave probability in phases 3 to 5 for nighttime heatwaves, and 7 to 8 for daytime heatwaves (Fig. 3c, d). Atmospheric water vapour itself is not directly linked to heatwave probability but its anomalies give an insight into other processes. Over the western Sahel, the increase of heatwave probability in phases 1 to 3 (Fig. 3c, d) can be connected to the negative anomaly of precipitable water (Fig. 6d), which is associated with reduced cloudiness and increased incoming solar radiation (not shown) which is a frequent underlying process of Sahelian heatwaves (Guigma et al. 2020). The reverse interpretation can be made for phases 5 to 8, where increased cloudiness and precipitable water (Fig. 6d) weaken incoming solar radiation (not shown) and are associated with a decrease of heatwave probability (Fig. 3c, d). As with the MJO, water vapour has a direct impact on heatwave intensity. The patterns of precipitable water (Fig. 6d) and of the modulation of heatwave intensity (Fig. 5c, d) are indeed relatively close.
The links between the EK wave-modulated patterns of heatwave occurrence and physical processes are less evident given that the anomalies of the EK wave are predominantly confined to the equatorial region as per theory (not shown). Nonetheless, it is apparent that advection of hot air (which is important in magnitude) in phases 1 to 5 plays an important role especially over the eastern Sahel (Fig. 6e), where heatwave probability is increased. On the other hand, advection of cooler air over the western Sahel in phases 4 to 7 may also explain the decrease of heatwave probability there.
While phase-longitude diagrams are efficient at summarising the physical conditions promoted by tropical modes in their different phases, they do not inform on their full horizontal distribution. To get around this issue, Fig. 7 shows for illustration purposes, the maps of these variables composited over the most convectively suppressed and enhanced phases (phases 1 and 5 respectively) of the MJO. They are superimposed with low-level (925 hPa) wind anomalies in order to assess the origins of the heat (whose climatological distribution is depicted in Lavaysse et al. 2009). The suppressed phase of the MJO is associated with a low level northeasterly wind anomaly across the Sahel (which is part of an anticyclonic cell) carrying the drier and hotter air from the Sahara desert and the eastern Sahel towards the western Sahel (Fig. 7a, c). As a result, in the western part of the Sahel, advection of hotter air (with magnitude reaching 1.5 K day − 1 ), leads to an increase of heatwave probability (Fig. 3a, b). On the other hand the arrival of air masses from higher latitudes causes drying and cooling over the eastern Sahel, reducing longwave radiation, hence the decline of heatwave probability in this region (Fig. 3a, b). Reversely, when the MJO is convectively enhanced, the low-level flow presents a cyclonic anomaly characterised by a northwesterly wind anomaly over the western Sahel/Sahara, which advects cooler air (Fig. 7d) leading to a decrease of heatwave probability (Fig. 3a, b). Over the eastern Sahel/Sahara, a southwesterly wind anomaly carries humid air from the equatorial region (Fig. 7b), causing an increase of heatwave probability (Fig. 3a, b) associated with increased longwave radiation.
These results agree with previous findings by Moron et al. (2018a) who also showed the importance of the modulation of low-level circulation by the MJO in increasing/decreasing temperature over tropical Africa. They also extend them further by documenting the associated physical processes. Guigma et al. (2020) showed that intraseasonal heatwaves over the Sahel are generally short-lived (average duration not exceeding five days), with only a small portion of them (about 10%) occurring concurrently at daytime and nighttime. Yet, in terms of impacts on health, longer lasting daytime events which persist during the nighttime are expected to be the most harmful (Schär 2016;Murage et al. 2017). Given their small sample size, the findings of the broad statistical analysis of heatwave modulation by tropical modes (Sect. 3) may not be directly applicable to them, especially since daytime and nighttime were dissociated. From this perspective, it is relevant to investigate them through a representative case study as also suggested by Guigma et al. (2020). This is the reason why the Sahelian heatwave which took place during the first half of April 2003 was selected. It is characterised by (i) a relatively long duration (more than 10 days) and (ii) a simultaneousness of daytime and nighttime heatwaves over a common wide domain. This event is also mentioned as a severe one in Oueslati et al. (2017) but, to the best of the authors' knowledge, it has never been investigated. The daytime heatwave started first, on 03 April, while the nighttime one started on 04 April, and they both ceased on 15 April. From a spatial point of view, the time-longitude diagram of temperature anomalies averaged over the Sahel band of 10° N-20° N (Fig. 8a) reveals that the heatwave kicked off over the eastern Sahel (in northern Chad/Sudan more precisely; not shown), stood there until approximately 08 April, moved towards the central Sahel before heading back to the east on its latest days. April 2003 was also marked by an important intraseasonal variability of temperature, as revealed by a wavelet decomposition of Tmin and Tmax (not shown). The variability of temperature at these scales is particularly relevant here as it suggests an important modulation by the high-frequency modes of climate variability, including the MJO and equatorial waves. Finally in April 2003, there was no incursion of extratropical perturbations into West Africa, nor any strong interannual temperature anomaly (not shown), therefore allowing to get a "pure" modulation by tropical modes.

Activity of tropical modes in April 2003
Figure8b shows a time-longitude diagram of highpass-filtered OLR anomalies averaged over the Guinean band to which the mode-filtered OLR anomalies (averaged over the same latitudinal band) are superimposed as contours. This is a common technique used in tropical mode investigations (e.g. Schreck and Molinari2011;Ventrice and Thorn-croft2012;Lafore et al. 2017).
Over the heatwave period (i.e. first half of April), two MJO packets crossed West Africa (black contours in Fig. 8b). For the first packet, it is its convectively suppressed phase which is relevant to the heatwave. It penetrated West Africa during the second half of March 2003 and reached the Indian Ocean near 10 April. As for the second packet, its convectively enhanced phase (which is the one relevant to the heatwave) entered West Africa in early April, reached

Overall physical processes during the event
An analysis of the thermodynamic budget emphasizes that two processes, namely longwave heating (Fig. 9a) and hot air advection (Fig. 9c) largely shaped the April 2003 heatwave. Indeed, during this period, the area covering Chad and Sudan is home to a strong hot air advection, reaching 1.5 K day − 1 at some places. As a result, the thermal gradient between the atmosphere and the ground is weakened, leading also to a decrease of the sensible heat flux according to ERA5 (not shown). As for the net longwave heating, it is observed over the western and central Sahel/ Sahara, and most importantly over the southeasternmost parts of the Sahel with a maximum exceeding 20 Wm − 2 . This longwave anomaly results mostly from a greenhouse effect (GHE) of water vapour, as testified by the pattern of precipitable water which presents a similar anomaly (Fig. 9b). Dust aerosol (an important constituent of the Sahelian atmosphere) was not significant in April 2003 according to MODIS data (not shown), precluding any significant contribution from this component to the observed warming.
The longwave warming anomaly at the surface is also a footprint of a reinforcement of convective activity over West Africa. Besides the positive anomaly of precipitable water, there is indeed a low-level cyclonic anomaly with monsoonal wind surges into the Sahel, advecting water vapour from the Gulf of Guinea (Fig. 9d). Furthermore, the patterns of low-level circulation, hot air advection and precipitable water observed during the heatwave period are very similar to their composite mean anomalies over the convectively enhanced phases of the MJO (Fig. 7b, d). As a reminder, heatwave probability is increased over the eastern Sahel when the MJO is in these phases (Fig. 3a,  b). As such, the April 2003 event, despite being exceptional (by its length and daytime/nighttime concomitance), adheres to the general findings presented in Sect. 3, both in terms of the physical conditions created by tropical modes and the outcome for heatwave occurrence. The next section emphasises this in more detail.

Day to day evolution of the heating processes and impact of tropical modes
A day to day analysis of the heatwave event further allows a deeper insight into the evolution of different processes shaping the heatwave and the role played by tropical modes.

An initial role of heat advection fostered by the MJO (late March − 04 April)
Figure10a shows that the onset of the heatwave over the eastern Sahel was preceded by a sustained advection of hot air in late March and during the first days of April. The low-level wind pattern in Fig. 9c shows that the hot air originated from the central and southern , where the West African Heat Low (WAHL) is typically located at this time of the year (Lavaysse et al. 2009). The projection of the activity of the different modes on this process (described in Sect. 2.4) shows that the MJO (which was in a convectively suppressed phase at that time; Fig. 8b) was its main driver. In Fig. 10a, the MJO-regressed anomalies (black contours) indeed encapsulate the heat advection. Therefore the MJO acted in suppressing convection over the WAHL area, causing a warming which is further advected to the northeastern Sahel and led to the heatwave initiation.

A progressive takeover by water vapour GHE led by the ER wave (05-08 April)
Despite a progressive decrease of heat advection after the onset, the heatwave was sustained mainly as a result of the arrival of moist air in the Sahel from the subequatorial region (Fig. 9d), causing a water vapour GHE as illustrated in Fig. 10b (using observed incoming longwave radiation at the surface, a proxy for GHE). In-situ observations at the Demokeya station (30.5° E; 13.3° N) also confirm the progressive increase of water vapour during this stage of the heatwave (Fig. SM5).
The regressed values of incoming longwave radiation (green contours in Fig. 10b) as well as of precipitable water (green contours in Fig. 10c) on the ER wave are positive and match relatively well the spatial structure of the anomalies of these variables, suggesting a preponderant role of this wave at this stage of the event. Thus, agreeingly with its expected impact on heatwaves discussed in early spring (Sect. 3), the convectively enhanced ER wave organised a low-level meridional transport of water vapour into the eastern Sahel, and, through longwave radiation, heated up the region. The statistical findings are further strengthened since consistency is shown here across different observational (satellite) datasets.

Final stage: a reinforcement of the GHE forced by the MJO (09-14 April)
The westward motion of the ER wave allowed the heating to reach the central and western parts of the Sahel. From 09  Fig. 8b). It further increased moisture through a zonal transport of water vapour from the Atlantic into the Sahel. As a result, the GHE effect reached its peak (above 30 Wm − 2 ) leading to an intensification of temperature anomalies, especially at night (Fig. 8a). The leading role of the MJO on this stage of the heatwave is discernible from the regression of the relevant processes on its timeseries (black contours in Fig. 10b, c). Its eastward propagation also dragged the region of maximum heating back towards the eastern Sahel. Again this is predicted by the statistical analysis as it is shown that, through water vapour GHE (Figs. 6b, 7b), convectively enhanced phases of the MJO heat up the eastern Sahel, augmenting heatwave probability therein (Fig. 3a, b). From approximately 15 April, the MJO-driven intensified convection caused organised rainfall events that watered large parts of the heatwave region (Fig. SM6) and led to a decrease of temperature. This eventually ended the heatwave.
It should be noted that the evolution of GHE and precipitable water depicted by ERA5 is very similar to that of CERES and AIRS respectively. Furthermore, the influence of the modes on the evolution of the physical processes is also well captured by ERA5. This is shown by the daily evolution of the regressed fields of interest spatially averaged over grid-cells where the heatwave is detected (Fig. SM7).
In summary, the April 2003 heatwave, from a large-scale point of view, was initiated by a convectively suppressed MJO that promoted hot air advection, and sustained through GHE of water vapour brought in, first by a convectively enhanced ER wave and later on by a convectively enhanced MJO. In addition to supporting the statistical characteristics of heatwave occurrence in the eastern Sahel, this case study also provided an insight on the organisation of the water vapour low-level flow by tropical modes which is an important driver of heatwaves. Thus, it appears that the MJO exerts its control mainly through zonal transport of water vapour whereas the ER wave mostly affects meridional transport.

Conclusion and perspectives
Extreme heat represents a growing threat for Sahelian populations, with the latest research stressing that by the end of the century, climate change would make the region reach worrying levels of thermal discomfort (Xu et al. 2020;Fig. 10 Time-longitude diagram of a ERA5 heat advection at the 925 hPa level. Shades represent the daily highpass filtered anomalies and contours the regression of these anomalies on the MJO (black contours) and ER waves (green contours). Units are in K day − 1 and the contour level is 0.1 (negative levels are omitted). b CERES incoming longwave radiation. Shades represent the daily highpass filtered anomalies and contours the regression of these anomalies on the MJO (black contours), the ER waves (green contours). Units are in W m − 2 and the contour levels are 2.5 and 5 (negative levels are omitted). c AIRS precipitable water. Shades represent the daily highpass filtered anomalies and contours the regression of these anomalies on the MJO (black contours) and ER waves (green contours). Units are in mm and the contour levels are 1 and 3 (negative levels are omitted). All the quantities are averaged over the 10° N-20° N latitudinal band between 15 March and 30 April 2003 Raymond et al. 2020). Therefore, undertaking actions which could mitigate this important issue is a top priority. This study contributes towards this through an analysis of the relationship between tropical modes and Sahelian heatwave events with the view of improving scientific understanding of these extremes. More precisely, it demonstrates that the MJO, ER and EK waves are important factors for the occurrence and -to a lesser extent -intensity of heatwaves in the region during the spring season. Their influence on heatwave occurrence is a function of their convective phase (enhanced or suppressed) and amplitude; it is also distinct for nighttime and daytime heatwaves and varies with the geographical location. Generally, in their convectively suppressed phases, they advect hot air from the climatological location of the WAHL towards the western coast of the Sahel, increasing heatwave probability there. When they are in convectively enhanced phases, they increase longwave warming especially over the eastern Sahel by bringing in water vapour from the Guinean region of Africa. Several tropical modes often overlap over the region. In these situations, heatwaves are more likely to occur over the western Sahel when all modes are in a convectively suppressed phase, and over the eastern Sahel when they are convectively suppressed. The statistical findings were further illustrated by the detailed analysis of a strong Sahelian heatwave event which occurred in April 2003. Hot air advection and longwave radiation, driven by the MJO and ER waves were responsible for a relatively long lasting event that had the particularity of affecting both daytime and nighttime.
From an operational point of view, the heatwaves-tropical modes relationship is important since the wave activity can be used by forecasters as an additional relevant tool to prepare warnings. Moreover, there is a potential for further work aiming at building statistical models of heatwave predictability based on the activity of tropical modes and the interaction between them. The modulation patterns at 5° E shown in the present study can already be used to forecast the likelihood of heatwave occurrence at a given location in the Sahel, taking into account the spectral properties (wavenumber and period) of the modes. For example, when a convectively suppressed ER wave is detected at 30° E, a forecaster in Senegal can expect hot conditions over his/ her region in the following 10 days. Better forecast recipes can be obtained from the monitoring of tropical mode activity by reassessing their impact when they are located at longitudes closer to the areas covered by the forecaster, but also at selected longitudes further away (because, given their wavenumbers, tropical modes can impact regions fairly remote from their central position). This therefore calls for more regional coordination between national meteorological services.
Besides, with the amelioration of weather observation networks, the increase of computational capacity and improvement of numerical models, tropical modes are now better represented and predicted (Dias et al. 2018;Janiga et al. 2018;Kim et al. 2018;Bengtsson et al. 2019). These skills can therefore be transferred to heatwave prediction at the intraseasonal scale and thus help to win time for preparedness actions. The authors are already engaged in research in that direction.
Other large-scale drivers of Sahelian heatwaves have been suggested by previous studies (Fontaine et al. 2013;Moron et al. 2018;Guigma et al. 2020). In particular, intrusions of extratropical Rossby waves often combined with tropical plumes can lead to significant increases of temperature and moisture in the Sahel as was the case in April 2010 (Largeron et al. 2020). It is therefore important that future studies also focus on the links between these large-scale features and heatwaves and on the interactions between them and tropical modes.
Finally, it may be important to assess whether the upward trend of Sahelian temperature and heatwave frequency and intensity over the recent past (Fontaine et al. 2013;Moron et al. 2016;Guichard et al. 2017;Ceccherini et al. 2017) is partly attributable to the global warming-induced increased activity of tropical modes over the same period (Song and Seo2016;Adames et al. 2017). This has implications for future projections, as the accuracy of the simulation and future evolution of springtime convection over the Guinean region by climate models might have repercussions on Sahelian heatwave predictions.
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://creat iveco mmons .org/licen ses/by/4.0/.