Differential Imprints of Distinct ENSO Flavors in Global Patterns of Very Low and High Seasonal Precipitation

The effects of El Niño's two distinct flavors, East Pacific (EP) and Central Pacific (CP)/Modoki El Niño, on global climate variability have been studied intensively in recent years. Most of these studies have made use of linear multivariate statistics or composite analysis. Especially the former assumes the same type of linear statistical dependency to apply across different phases of the El Niño–Southern Oscillation (ENSO), which appears not necessarily a justified assumption. Here, we statistically evaluate the likelihood of co-occurrences between very high or very low seasonal precipitation sums over vast parts of the global land surface and the presence of the respective EP and CP types of both, El Niño and La Niña. By employing event coincidence analysis, we uncover differential imprints of both flavors on very low and very high seasonal precipitation patterns over distinct regions across the globe, which may severely affect, among others, agricultural and biomass production or public health. We particularly find that EP periods exhibit statistically significant event coincidence rates with hydrometeorological anomalies at larger spatial scales, whereas sparser patterns emerge along with CP periods. Our statistical analysis confirms previously reported interrelations for EP periods and uncovers additional distinct regional patterns of very high/low seasonal precipitation, such as increased rainfall over Central Asia alongside CP periods that have to our knowledge not been reported so far. Our results demonstrate that a thorough distinction of El Niño and La Niña into their two respective flavors could be crucial for understanding the emergence of strong regional hydrometeorological anomalies and anticipating their associated ecological and socioeconomic impacts.

The effects of El Niño's two distinct flavors, East Pacific (EP) and Central Pacific (CP)/Modoki El Niño, on global climate variability have been studied intensively in recent years. Most of these studies have made use of linear multivariate statistics or composite analysis. Especially the former assumes the same type of linear statistical dependency to apply across different phases of the El Niño-Southern Oscillation (ENSO), which appears not necessarily a justified assumption. Here, we statistically evaluate the likelihood of co-occurrences between very high or very low seasonal precipitation sums over vast parts of the global land surface and the presence of the respective EP and CP types of both, El Niño and La Niña. By employing event coincidence analysis, we uncover differential imprints of both flavors on very low and very high seasonal precipitation patterns over distinct regions across the globe, which may severely affect, among others, agricultural and biomass production or public health. We particularly find that EP periods exhibit statistically significant event coincidence rates with hydrometeorological anomalies at larger spatial scales, whereas sparser patterns emerge along with CP periods. Our statistical analysis confirms previously reported interrelations for EP periods and uncovers additional distinct regional patterns of very high/low seasonal precipitation, such as increased rainfall over Central Asia alongside CP periods that have to our knowledge not been reported so far. Our results demonstrate that a thorough distinction of El Niño and La Niña into their two respective flavors could be crucial for understanding the emergence of strong regional hydrometeorological anomalies and anticipating their associated ecological and socioeconomic impacts.

INTRODUCTION
The positive (El Niño) and negative (La Niña) phases of the El Niño-Southern Oscillation (ENSO) are known to be associated with wide-spread anomalies in the mean hydrometeorological conditions at various distant parts of the Earth. These long-ranged interactions are often referred to as teleconnections (Trenberth, 1997;Neelin et al., 2003;Domeisen et al., 2019). In this context, recent findings indicate that there exist two distinct types or flavors of El Niño phases, usually referred to as the East Pacific (EP) or canonical El Niño and the Central Pacific (CP) El Niño or El Niño Modoki (Ashok et al., 2007;Kao and Yu, 2009;Marathe et al., 2015), respectively. It has been shown that these two flavors are possibly associated with distinct hydrometeorological responses in certain regions (Taschetto and England, 2009). Examples include reduced rainfall over eastern Australia (Chiew et al., 1998) or Southern Africa (Ratnam et al., 2014) only during EP periods as well as increased precipitation over the tropical regions of Africa (Preethi et al., 2015) or the western Indo-Pacific Oceans (Weng et al., 2011;Feng and Chen, 2014) during CP periods. The main reason for such differential responses are different longitudinal positions of the strongest ENSO related sea surface temperature (SST) anomalies in the tropical Pacific Ocean alongside the different flavors, which results in different modifications of large-scale atmospheric circulation patterns (Ashok and Yamagata, 2009;Domeisen et al., 2019).
A similar discrimination into two types has been suggested for La Niña phases as well (Kug and Ham, 2011), even though their respective imprints on SST patterns seem less distinct than for El Niño (Kao and Yu, 2009;Ren and Jin, 2011). It therefore remains an open problem to provide further statistical and/or dynamical evidence either in favor or against such a distinction of different La Niña flavors in analogy to El Niño (Chen et al., 2015;Zhang et al., 2015;Levine et al., 2016). Still, differential hydrometeorological responses alongside specific La Niña phases have recently been identified especially across the Pacific (Shinoda et al., 2011;Magee et al., 2017;Song et al., 2017;Hidayat et al., 2018), which provides good reason for discriminating La Niña into two types as well.
Most previous studies (including those mentioned above) on the teleconnective impacts of different ENSO phases and flavors have either applied linear statistical tools, such as (partial) correlation analysis (Diaz et al., 2001;Weng et al., 2011;Preethi et al., 2015;Magee et al., 2017;Hidayat et al., 2018) or empirical orthogonal function (EOF) analysis (Taschetto and England, 2009), or investigated composites (i.e., mean spatial patterns for a specific type of ENSO period) of the corresponding climate observable of interest (Feng and Chen, 2014;Hoell et al., 2014;Song et al., 2017). Specifically, the former methods share the common limitation of focusing on linear or average interdependencies between ENSO and possible response variables. At the same time, global climate change has been projected to lead to an increase in the strength and frequency of both, climate extremes (Easterling et al., 2000;Karl and Trenberth, 2003) as well as extreme ENSO phases (Cai et al., 2014(Cai et al., , 2015. This calls for more systematically assessing possible statistical as well as dynamical/mechanistic linkages between these two findings (Allan and Soden, 2008). Accordingly, the present work aims to identify and quantitatively characterize spatial patterns of markedly wet/dry seasons that have an elevated probability to co-occur with certain types (EP or CP) of ENSO phases, following upon previous findings that ENSO can have large-scale effects on rainfall patterns at both, global and regional scales (Ropelewski and Halpert, 1987;Dai and Wigley, 2000).
Several strategies to distinguish East Pacific (EP) from Central Pacific (CP) ENSO events have been proposed in the recent past (Hendon et al., 2009). One prominent example is the ENSO Modoki Index that is computed as weighted average SST anomalies over three specific regions (165 • E to 140 • W and 10 • S to 10 • N, 110 • W to 70 • W, and 15 • S to 5 • N, as well as 125 • E to 145 • E and 10 • S to 20 • N) in the equatorial Pacific (Ashok et al., 2007). Other approaches have used empirical orthogonal functions (Kao and Yu, 2009;Graf and Zanchettin, 2012) or combinations of the Nino3 and Nino4 indices (Hu et al., 2011;Kim et al., 2011) to provide the desired categorization. Wiedermann et al. (2016) recently compiled a synthesis of eight previous studies that used either of the aforementioned frameworks and identified several ENSO periods during which these techniques proposed either ambiguous, mutually inconsistent, or incomplete classifications, with certain ENSO periods remaining entirely unassigned. In order to fill those gaps, Wiedermann et al. (2016) proposed a new index based on spatial correlation structures among global surface air temperature anomalies, which can be conveniently studied in terms of the transitivity of so-called functional climate networks (Tsonis et al., 2006;Donner et al., 2017). This transitivity index confirmed the flavors of all El Niños between 1953 and 2013 that had been classified in mutual agreement across the existing literature (Wiedermann et al., 2016). It also allowed to assign types to those cases where former work yielded incomplete or ambiguous categorizations and thereby provided a consistent, comprehensive, and complete classification of the respective flavors. Moreover, Wiedermann et al. (2016) showed that the transitivity index also naturally distinguishes La Niña episodes into two corresponding types, thereby providing a unique advantage over other classification schemes that mainly focus on El Niño phases alone. Even though the existing literature on a discrimination of different La Niña periods is comparatively scarce, the transitivity index confirmed the results of at least two recent studies (Tedeschi et al., 2013;Yuan and Yan, 2013) and again provides a consistent classification for years that previously had types ambiguously or not all assigned (Wiedermann et al., 2016). For both reasons, i.e., the comprehensiveness of the classification and the ability to distinguish La Niña into two types as well, we will directly use the classified periods from Wiedermann et al. (2016) (see Table 1 for an overview) for the purpose of our present study.
Based on this categorization of El Niño and La Niña phases, we quantify the likelihood of simultaneous or timedelayed co-occurrences of strong seasonal wet/dry anomalies at a local scale across vast parts of the global land surface with a certain type of ENSO phase. Specifically, we consider seasonal precipitation sums for boreal fall (September-November, SON),  Wiedermann et al. (2016) on the basis of a comprehensive literature synthesis and a consistent classification using the network-based transitivity index.

2014*
The given years correspond to the onset-year of each ENSO event such that, e.g., 1951 indicates the 1951/1952 El Niño event. Years marked with an asterisk are classified using the same methodology as in Wiedermann et al. (2016) to extend the data to the period of study considered in this work.
winter (December-February, DJF), and spring (March-May, MAM) as key seasons of the developmental cycle of El Niño and/or La Niña conditions. By contrast, we omit the boreal summer season, even though large scale hydrometeorological conditions during this season could still be distinctively affected by different (approaching or withdrawing) ENSO phases. The reason for this choice is that there are various cases where El Niño and La Niña conditions occur in subsequent years so that a unique attribution of anomalies to any of these two phases would be hardly possible. In our present analysis, strong wet (dry) anomalies will be defined as seasonal precipitation sums exceeding the empirical 80th (falling below the empirical 20th) percentile of the distribution of all values on record for a given season and location.
To statistically quantify co-occurrences between different types of ENSO phases and very wet/dry seasons, we employ event coincidence analysis (ECA) (Donges et al., 2011. Put in simple terms, ECA counts the fraction of events of one type (here, some particularly wet/dry season at a certain location) that coincide with those of another type (here, some type of ENSO period) while, in contrast to other conceptually related approaches like event synchronization (Quian Quiroga et al., 2002;Malik et al., 2010;Boers et al., 2014), allowing for a precise control of the relative timing between them (Wolf et al., 2020). This statistical framework has already been successfully applied to quantify the likelihood of climatic (extreme) events possibly triggering certain ecological or socioeconomic responses, such as extreme annual (Rammig et al., 2015) and daily (Siegmund et al., 2016a) tree growth or flowering dates (Siegmund et al., 2016b), the outbreak of epidemics  or armed conflicts Ide et al., 2020).
In the course of this work, we first evaluate significant event coincidence rates for the canonical (EP) El Niño. This allows us to demonstrate the consistency of our approach by comparing the obtained spatial patterns with a variety of previously reported results. We continue by studying co-occurrences between strong seasonal precipitation anomalies and the so far less intensively studied CP El Niño periods to highlight differences in their teleconnectivity patterns as compared to their EP counterparts. Ultimately, we also study La Niña periods and demonstrate that the most remarkable large-scale spatial patterns of strong seasonal rainfall anomalies non-randomly co-occurring with this negative ENSO phase are associated with the corresponding EP flavor. This finding suggests that in the light of recent discussions on the existence of two distinct types of La Niña periods (Kao and Yu, 2009;Kug and Ham, 2011;Ren and Jin, 2011;Zhang et al., 2015), it is indeed meaningful to consider a global impactbased distinction between one type that significantly affects seasonal wet/dry patterns globally and another that exhibits less spatially coherent impacts. Wherever appropriate, we also briefly discuss possible ecological and socioeconomic consequences of the identified seasonal precipitation anomalies.

GPCC Rainfall Data
We utilize gridded monthly precipitation data provided by the Global Precipitation Climatology Centre (GPCC) at a spatial resolution of 2.5 • × 2.5 • (Schneider et al., 2018). Since reliable discriminations of El Niño and La Niña periods into their respective EP and CP flavors are so far mainly available for the second half of the twentieth century (Graf and Zanchettin, 2012;Yuan and Yan, 2013;Wiedermann et al., 2016;Freund et al., 2019), we restrict our analysis to the period from 1951 to 2016 (the most recent year in the GPCC data set). We derive separate records for three seasons s by aggregating the precipitation sums of the corresponding 3-month periods from September to November (SON, 1951(SON, -2015, December to February (DJF, 1951(DJF, /52-2015, and March to May (MAM, 1952(MAM, -2016 to cover all seasons associated with the 1951/52 El Niño and all following ENSO phases while ensuring the same length of all series. This results in three time series P i,s (t) per grid cell i with M = 65 annual values each. Note that the density of stations from which the GPCC data has been derived varies between 0 and more than 100 per grid cell and year (Lorenz and Kunstmann, 2012), which generally results in a lower accuracy and reliability of the data for those areas with only few stations (Rudolf et al., 1994). We therefore consider only grid cells with at least one station present for 95% of the study period in a specific season s (SON, DJF, or MAM). In addition, we exclude those grid cells where the average precipitation sum in a specific season is below 3 cm (i.e., 1 cm per month) since this particular choice of threshold ensures the exclusion of deserts from our analysis (Chatterton et al., 1971;Houston, 2006;Thomas and Nigam, 2018). Both preprocessing steps yield a total number of N SON = 1763, N DJF = 1, 610, and N MAM = 1, 736 valid grid cells in SON, DJF, and MAM (see Supplementary Figure 5 for details on the spatial distribution of valid grid cells).
We note that in the context of the present work, using 3month precipitation sums provides information fully equivalent to that obtained when using the 3-month standardized precipitation index SPI-3 (Guttman, 1999;Svoboda et al., 2012), which is a solely precipitation-based characteristic commonly employed in drought-related studies. Specifically, SPI-3 values can be derived by a sophisticated monotonic transformation of the 3-month precipitation sums, which does not change the rank order of the values and, hence, the timing of the events considered in this study.

Classification of ENSO Periods
We use the classification of East Pacific (EP) and Central Pacific (CP) flavors of recent El Niño and La Niña episodes from Wiedermann et al. (2016) (Table 1). This classification has been based on a comprehensive literature review of eight (two) studies that distinguished certain El Niños (La Niñas) between 1953 and 2010 into their two respective types. Based on this compilation, Wiedermann et al. (2016) identified 11 ENSO periods for which previous works yielded mutually consistent results and another 10 periods for which (a subset) of previous studies yielded incomplete, ambiguous, or mutually inconsistent classifications.
In order to achieve a consistent and comprehensive classification of both, El Niño and La Niña flavors, including also the ENSO periods without previous consensus, Wiedermann et al. (2016) introduced the so-called transitivity index, which reflects distinct characteristics (in terms of abundance and localization) of ENSO's teleconnections during EP and CP periods (Radebach et al., 2013;Wiedermann et al., 2016). In brief, this index is obtained from 1-year sliding-window lagzero absolute correlation matrices between time series at all pairs of grid cells in the global daily surface air temperature anomaly field (Kistler et al., 2001). The most relevant information of these matrices is contained in the 0.5% strongest absolute correlations (Donges et al., 2009;Radebach et al., 2013;Wiedermann et al., 2016), yielding thresholds for each considered time window below which the matrix coefficients are put to zero. The thus obtained sparse matrices are then considered as weighted adjacency matrices of so-called functional climate networks (Donges et al., 2009;Radebach et al., 2013;Wiedermann et al., 2016) with positive entries indicating a strong statistical relationship between climate variability in two grid cells. The transitivity index (Newman, 2003;Antoniou and Tsompa, 2008) then describes the degree to which such strong relationships among triples of grid cells are transitive, i.e., the fraction of cases in which connections between two pairs (i, j) and (i, k) of grid cells sharing a common member are accompanied by a third connection between the remaining pair (j, k). In order to provide a representative property, this fraction is further weighted by the respective cell-sizes and strengths of the involved links, i.e., the corresponding values of the absolute cross-correlation (Saramäki et al., 2007). The transitivity index quantifies the (dispersed or localized) spatiotemporal distribution of links in the climate network and directly reflects distinct characteristics in the temporal evolution of spatial autocorrelations and global teleconnections that are unique to specific ENSO flavors (Radebach et al., 2013;Wiedermann et al., 2016).
Given the presence of either El Niño or La Niña conditions as indicated by the Oceanic Nino Index (ONI), the transitivity index then indicates EP phases by a strong peak co-occurring with the respective ENSO period. By contrast, CP phases can be identified by the absence of such a transitivity peak with values of the index close to its baseline (see Wiedermann et al., 2016 for a comprehensive description and interpretation of the framework as well as all necessary mathematical details).
The transitivity index confirmed the classification of all 11 consistently reported El Niño periods from earlier works and additionally provided a comprehensive classification for those periods that were previously classified inconsistently or incompletely. While most other recent approaches for classifying ENSO flavors were mostly tailored to El Niño events, the transitivity index also provides a consistent classification of La Niña periods (Wiedermann et al., 2016), thereby making it particularly useful for our present study.
In total, Wiedermann et al. (2016) identified six EP El Niño and sixteen CP El Niños between 1951 and 2014, which is largely consistent with a recent study by Freund et al. (2019) that used objective supervised machine learning to obtain a similar discrimination. Similarly, Wiedermann et al. (2016) found seven EP La Niñas and 11 CP La Niñas. Since Wiedermann et al. (2016) only provide classifications for all ENSO periods prior to the year 2014, we apply their methodology to extend the event categorization to our present study period. We specifically identify the 2014/2015 and 2015/2016 El Niño periods as CP and EP types, respectively, which is again consistent with the recently proposed classification by Freund et al. (2019). An overview of all ENSO periods and types that are ultimately used in this study is given in Table 1.
Based on the classification in Table 1, we create binary indicator time series for the four different ENSO phases (see Figure 1), e.g., an EP El Niño series X EPEN (t) with X EPEN (t) = 1 if t marks the onset-year of an EP El Niño and X EPEN (t) = 0 otherwise (solid lines in Figure 1A). Correspondingly, we obtain the event series X CPEN (t) for CP El Niños (dashed lines in Figure 1A). The same procedure is applied to La Niña periods, resulting in two event series X EPLN (t) and X CPLN (t), respectively ( Figure 1B).

Data Preprocessing
We identify years with seasons s (DJF, SON, or MAM, see above) exhibiting extraordinary high or low precipitation amounts from the corresponding time series P s,i (t) for each grid cell i, individually. Specifically, we consider values above (below) the 80th (20th) percentile p + s,i (p − s,i ) in each of the time series P s,i (t) as extraordinary high (low) seasonal precipitation sums (Figure 2A). We choose these particular thresholds to ensure the presence of a sufficient number of particularly dry and wet seasons that is comparable with the number of different types of ENSO periods in the considered study period. It has been FIGURE 2 | (A) Exemplary seasonal precipitation time series P s,i (t) (colored bars) as well as corresponding thresholds p + si and p − si computed as the 80th and 20th percentile of P s,i (t), respectively. Precipitation signals above p + si (blue bars) are considered very wet while those below p − s,i (red bars) are considered very dry for the subsequent analysis. (B) The two corresponding event series P + s,i (t) and P − s,i (t) with P ± s,i (t) = 1 if a corresponding event occurs at time t and P ± s,i (t) = 0 otherwise.
checked that the results obtained in this work do not change qualitatively if more restrictive or loose thresholds are applied (see Supplementary Material for details). According to these considerations, we obtain six binary (indicator) time series P ± s,i (t) for each GPCC grid cell i, where P + s,i (t) = 1 (P − s,i (t) = 1) indicates the presence of a very high (very low) seasonal precipitation sum at grid cell i during season s in year t (Figure 2). By following the above procedure, no further deseasonalization of the precipitation data is necessary, since the grid cell specific seasonality of precipitation is already taken into account. Furthermore, the events are defined for each grid cell independent of all the others, which means that the specific characteristics of rainfall variability and strength at one location do not influence the definition of events at other locations as it would be the case when considering, e.g., a global (location-independent) threshold for seasonal precipitation.

Event Coincidence Analysis
Event coincidence analysis (ECA) is a statistical tool that quantifies the empirical likelihood of co-occurrences of events in two series (Donges et al., 2011Rammig et al., 2015). To complement other conceptually related approaches like event synchronization (Quian Quiroga et al., 2002;Malik et al., 2010;Boers et al., 2014), ECA has been developed based on analytical considerations on paired point processes and combines a precise control of the relative timing of (instantaneous or mutually lagged) events that are considered synchronous with analytical confidence bounds on the obtained event coincidence rates (Donges et al., 2011. Both aforementioned features provide some considerable benefit of ECA in the specific context of the present work not shared by other similar methods. Particularly the precise control of the relative timing between events (along with the specific ability to chose a global coincidence interval; Wolf et al., 2020) is crucial for our analysis. It ensures that the we only consider precipitation anomalies that occur simultaneously (for SON or DJF) or exactly a year subsequent (for MAM) to the onset of a respective ENSO phase.
In the context of our work, ECA provides for each grid cell the fraction (called event coincidence rate) of EP/CP ENSO periods that co-occur with very high or low precipitation sums in SON or DJF of the same year or MAM of the year following the onset of an El Niño or La Niña phase. Hence, the event coincidence rate ECR ± s,i,• for one pair of ENSO and precipitation event series is given by . (2) Here, X • (t) represents one of the four time series indicating the presence of EP and CP flavors of El Niño and La Niña (see above). While interpreting t as the calendar year, the offset τ reads τ = 0 for SON and DJF and τ = 1 for MAM. Note, that our present analysis studies so-called trigger coincidence rates  that quantify the likelihood of a given ENSO period and phase to be followed by a specific strong/wet precipitation signal. In contrast, a complementary definition of ECR (denoted precursor coincidence rate; Donges et al., 2016) would address the inverse problem of quantifying likelihoods that a given precipitation event (that could arise through a variety of conditions and drivers) is preceded by a specific ENSO period and phase. However, since our present work only considers instantaneous coincidences (see Donges et al., 2016 for details), i.e., coincidences between events in the same ENSO period, the two types of coincidence rates only differ in their normalization [i.e., the denominator of Equation (2)]. They therefore provide essentially similar information, except for possible differences in the associated statistical significance resulting from different numbers of events (i.e., different sample sizes).
To assess the statistical significance of the empirical event coincidence rates, we assume both involved event sequences to be distributed randomly, independently and uniformly (Donges et al., 2011. A corresponding p-value is derived analytically from the probability distribution of event coincidence rates that would occur by chance only. We consider an empirical event coincidence rate as statistically significant if its associated p-value is smaller than a confidence level of α = 0.05 .

Seasonal Wet/Dry Patterns and EP El Niño
We first investigate co-occurrences of EP El Niño periods and very wet/dry seasons. Figures 3A,C,E highlight areas with significant event coincidence rates between EP El Niños and very dry (red squares) and wet seasons (blue squares) in SON, DJF, and MAM, respectively.
During those SON seasons that correspond to the developing stages of EP El Niños, we find an elevated probability of very dry conditions over Indonesia, the Philippines, and the southwestern Pacific Islands as well as over northern South America and the northern Amazon Basin (Figure 3A). Droughts in the latter region have been previously linked to an increased risk of biomass loss in the Amazon which normally serves as a long-term carbon sink (Phillips et al., 2009;Lewis et al., 2011).
For the same season (SON), we also observe an increased likelihood of very wet conditions along the west coast of North America (Figure 3A). Similarly, unusually wet conditions frequently emerge over Ecuador and southeastern South America in SON and DJF (Figures 3A,C). We further observe more spatially confined regions with wet conditions over parts of the Chilean Andes in SON (Figure 3A), which may result in an increased risk for the occurrence of floods in this area (Bookhagen and Strecker, 2012;Boers et al., 2014). All these observations agree well with previous studies (Diaz et al., 2001).
Coinciding with EP El Niños, we also observe more frequent wet conditions over the Mediterranean region and East Africa during SON (Figure 3A), both of which have been previously reported in local case studies (Camberlin et al., 2001;Shaman and Tziperman, 2010). The observed tendency toward very dry conditions in southwestern Africa during DJF ( Figure 3C) has also been reported recently (Hoell et al., 2014).
For MAM seasons, we observe pronounced large-scale patterns of significant event coincidence rates between EP El Niños and low precipitation sums in Northeast Brazil (Figure 3E), which is consistent with previous studies (Kane, 1997). Furthermore, the dry conditions over the Philippines that are observed during SON (see above) and DJF also persist into the MAM season ( Figure 3E). In addition, strong MAM rainfall occurring alongside EP El Niños is most prominently observed in the southeastern United States, which is again consistent with previous works that used composite analysis to determine North American weather patterns associated with El Niño conditions (Ropelewski and Halpert, 1986).
Taken together, the results of ECA are overall in good agreement with previously reported interrelations between El Niño and global precipitation patterns, which have mostly been identified using linear statistical tools such as correlation analysis (Camberlin et al., 2001;Phillips et al., 2009), principal component analysis (Diaz et al., 2001), or composites based on seasonal averages (Shaman and Tziperman, 2010;Hoell et al., 2014). Thus, we conclude that the application of ECA to unveil strong ENSO related hydrometeorological anomalies provides consistent results when compared to previous studies. Our results also imply that strong responses of seasonal precipitation to canonical (EP) El Niño conditions show similar spatial patterns as the average statistical interdependency between ENSO-related indices and the corresponding hydrometeorological observables. However, we note that most previous studies have not discriminated between the two El Niño flavors. Thus, the agreement between our results for EP El Niños and the existing literature suggests that the FIGURE 3 | Statistically significant event coincidence rates (ECR) between EP (A,C,E) and CP (B,D,F) El Niños and very dry (red squares) or very wet (blue squares) conditions for the three seasons SON, DJF, and MAM. Dry/wet periods are defined by seasonal precipitation sums below/above the 20th/80th percentile of all years from 1951 to 2016. Only significant clusters of at least two adjacent grid-cells are shown. Different shades of red/blue indicate increasing magnitudes of significant ECR between the respective El Niño events and strong/weak seasonal precipitation. Yellow areas indicate grid cells with non-significant event coincidence rates. White areas over land indicate insufficient quality of the GPCC dataset (i.e., grid cells excluded from our analysis).
previously observed (linear) effects might be dominated by the (on average stronger; Kug et al., 2009;Huang et al., 2016) EP events.

Seasonal Wet/Dry Patterns and CP El Niño
Next, we focus on very high and low seasonal precipitation along with the so far less intensively studied CP El Niño (cf. Figures 3B,D,F). We first discuss those regions that display significant event coincidence rates for EP El Niños (see above) but not for the corresponding CP periods. Notably, very dry seasons over tropical South America that frequently co-occur together with EP El Niños are markedly less prominent for CP El Niños during DJF ( Figure 3D) and do not display any significant event coincidence rates at all in SON (Figure 3B). The latter also holds true for very wet seasons along the western coast of Central and North America that have been observed for EP El Niños. In the same manner, the wet SON patterns over southern China, the Mediterranean, and East Africa frequently co-occurring with EP El Niños in SON do not exhibit significant coincidence rates with CP El Niños (Figure 3B). We further note, that the large-scale dry events over Indonesia observed along with EP El Niños during SON become less spatially coherent for CP El Niños ( Figure 3B). For MAM, the wet patterns over the northern Iberian peninsula that are observed for EP El Niños cannot be identified for CP El Niños (Figure 3F).
While the aforementioned observations indicate decreased or weakened impacts of CP El Niños in comparison with the EP flavor, we also observe new additional patterns of significant event coincidence rates that are not present during EP El Niños but emerge only along with CP periods. Most notably, very dry conditions become more likely along Australia's east coast during SON (Figure 3B). Such hydrometeorological anomalies could thus result in severe impacts on river ecosystems and agriculture in that region (Leblanc et al., 2012). In particular, marked drought phases in Eastern Australia are likely to cause a cascade of low river inflows, general water scarcity and large scale floodplain forest mortality as well as an increase of toxicity in the surrounding lakes (Leblanc et al., 2012). These natural hazards can have substantial effects on agricultural production in terms of a severe reduction in irrigated crop yields (van Dijk et al., 2013). In addition to reduced rainfall responses, significant event coincidence rates with very wet conditions are found over southern Chile pointing toward increased rainfall during CP El Niños as compared to their canonical counterparts. In DJF months coinciding with CP El Niño periods, we observe the emergence of new wet patterns over Central Asia well as a dry pattern over the north of Peru and Ecuador ( Figure 3D). Finally, we observe a pronounced dry pattern over Southeast Africa in MAM ( Figure 3F).
Generally, we note lower event coincidence rates between seasonal wet/dry conditions and CP El Niños as compared to EP periods (compare Figures 3A,C,E and Figures 3B,D,F). This might be partly explained by the larger of number of 17 CP events as compared to 7 EP events over the considered study period which is consistent with previously reported frequencies of the two flavors (Hendon et al., 2009;Graf and Zanchettin, 2012;Preethi et al., 2015;Wiedermann et al., 2016). In addition, recent studies suggest that CP El Niños might be further discriminated into two subtypes based on their specific impacts on Pacific rainfall and the modes of the Indian Ocean dipole Wang, 2013, 2014;Wang et al., 2018). Along those lines, future work should therefore investigate whether the comparatively lower significant event coincidence rates between CP El Niños and strong/weak seasonal precipitation can indeed also be attributed to the presence or absence of any of these two possible subtypes.

Seasonal Wet/Dry Patterns and EP/CP La Niña
Ultimately, we perform the same analysis as above for La Niña periods. For the EP, i.e., canonical, La Niña phases (Figures 4A,C,E) we again find various patterns that have already been reported in previous studies. Specifically, during SON coinciding with EP La Niñas (Figure 4A), we recapture an increased probability for very wet conditions over Australia and Indonesia (Arblaster et al., 2002) and exceptionally dry conditions in southern Europe (Pozo-Vázquez et al., 2005) and the south of Brazil and Uruguay (Ropelewski and Halpert, 1996). We further observe significant event coincidence rates for dry seasons in the Middle East contrasted by more intense than normal rainfall over central Europe. For DJF, our analysis confirms previous findings of very wet conditions over South Africa and dry episodes over West Africa (Nicholson and Selato, 2000) (Figure 4C), the latter having been previously linked to potential agricultural losses (Karpouzoglou and Barron, 2014) and considerable health risks in that specific area (Rataj et al., 2016). We further observe a prominent seasonal precipitation dipole with dry conditions over Mexico and elevated rainfall over southwestern Canada. The latter has become an important aspect of local water resource management (Lute and Abatzoglou, 2014), but together with increasing air temperature and more frequent storms also poses the threat of landslides in corresponding coastal areas (Guthrie et al., 2010). For MAM seasons associated with EP La Niñas, we observe a tendency towards strong rainfall over the Amazon (Rogers, 1988) and parts of Northern Australia (Arblaster et al., 2002) and the Philippines (Figure 4E).
In contrast to the wide-spread spatially coherent wet/dry anomaly patterns observed for EP La Niña periods, we find much fewer spatially extended structures along with CP La Niñas (Figures 4B,D,F). Most prominently, we recover previously reported wet conditions over parts of Australia in SON (Arblaster et al., 2002;Cai and Cowan, 2009) (Figure 4B). Additionally, we uncover strongly reduced rainfall over Florida in DJF (Figure 4D), and over the United Kingdom, Ireland and the west of Kazakhstan in MAM ( Figure 4F). In summary, we observe that CP La Niñas co-occur with less spatially coherent precipitation responses as compared to their canonical counterparts. We also note that event coincidence rates between La Niña and seasonal wet/dry precipitation signals are quantitatively more similar across EP and CP periods than El Niño (cf. Figures 4A-F), which could be partly explained by the comparatively lower difference between the total number of 8 EP and 12 CP La Niña events. In addition, future work should again investigate the potential for further discriminating CP La Niñas into two distinct subtypes in analogy to El Niño Wang, 2013, 2014;Wang et al., 2018). Such an approach would possibly allow us to attribute significant event coincidence rates between CP La Niña and strong/weak seasonal precipitation signals to the presence or absence of a certain event type and help explain the different magnitude in significant event coincidence rates observed in Figure 4.

DISCUSSION AND CONCLUSIONS
We have carried out a detailed analysis of ENSO imprints in global patterns of very wet/dry seasons over land. Specifically, we distinguished both, El Niño and La Niña, into two distinct flavors (East Pacific and Central Pacific) by utilizing a classification based on an extensive literature review paired with the assessment of a recently proposed complex network-based index (Wiedermann et al., 2016). From this classification, we obtained event series describing the occurrence-times of the four distinct types of ENSO. Strong seasonal precipitation anomalies have been obtained from the globally gridded GPCC rainfall data set by identifying seasons with precipitation sums above the empirical 80th (below the empirical 20th) percentile of all values for a given grid point as very wet (very dry) periods. This definition follows the spirit of the 3-month aggregate standardized precipitation index (SPI-3), which just provides a monotonic rescaling of seasonal precipitation sums according to the local distributional characteristics of seasonal rainfall. Accordingly, our results can be interpreted in terms of seasonal drought characteristics. Modifying the considered thresholds for defining very wet/dry seasons within reasonable ranges did not qualitatively alter the results reported in this work (see Supplementary Material).
We have then employed event coincidence analysis Siegmund et al., 2016a) to identify grid points with significant event coincidence rates between different types of ENSO phases and very high or low seasonal precipitation sums. Our analysis confirmed that previously observed interrelationships based on linear correlation or composite analysis in many cases also apply to the timing of events corresponding to the tails of the probability distributions of seasonal precipitation sums. In addition, we identified further patterns of very wet/dry conditions with elevated probabilities alongside different ENSO types, which have to our best knowledge not been described so far. These include increased rainfall over Central Asia during CP El Niños contrasted by rainfall reduction over the same area during CP La Niñas (compare also Figures 5B,D) which implies that for such cases standard (linear) statistical methods seem to be insufficient to unveil the underlying interrelations of events. Moreover, our present analysis demonstrates that even though a general linear relationship between ENSO and precipitation might be relatively weak or even absent for a given region, dry periods or very wet seasons can still be possible consequences of the presence of certain ENSO phases. At the same time we observe that several previously reported links between ENSO and global precipitation patterns mostly apply to EP ENSO flavors as far as particularly wet/dry seasons are concerned.
We found that the CP flavor of El Niño (and in parts also that of La Niña) significantly, i.e., non-randomly, cooccurs alongside wet/dry seasonal precipitation signals at lower event coincidence rates which might in parts be explained by the more frequent occurrence (Hendon et al., 2009;Graf and Zanchettin, 2012;Preethi et al., 2015;Wiedermann et al., 2016) compared to its EP counterpart. This finding also suggests that the absolute number of event coincidences [i.e., the numerator in Equation (2)] between EP or CP ENSO periods with wet/dry seasonal precipitation per grid-cell might be of rather comparable size over the period of study, an effect that should be investigated more thoroughly in future research. In addition, recent studies suggest that the CP El Niño can be further discriminated into two types with distinct imprints on precipitation signals especially in the Indo-Pacific area Wang, 2013, 2014;Wang et al., 2018). Future work should hence investigate if the comparatively lower significant event coincidence rates between CP El Niño and wet/dry seasonal conditions can potentially be attributed to the presence or absence of a further subtype of CP El Niño. Since El Niño and La Niña show a large degree of symmetry in terms of their frequency and potential for discrimination into different types (Ashok et al., 2007;Kao and Yu, 2009;Hidayat et al., 2018), a similar analysis could also be performed for La Niña as we found lower, yet significant, event coincidence rates for the CP flavor of ENSO's negative phase as well. Along those lines, our analysis provides a complementary impact oriented view on the recently raised question whether it is actually (statistically and/or dynamically) meaningful to distinguish La Niña into two types in a similar fashion as El Niño (Kug and Ham, 2011). While some previous studies advocated for such a distinction (Ashok and Yamagata, 2009;Kao and Yu, 2009), others argued that based on correlation analyses between La Niña related SST patterns there is a lack of evidence for the existence of two distinct types (Kug et al., 2009;Ren and Jin, 2011). Contributing to this discussion, Wiedermann et al. (2016) already demonstrated that according to the transitivity index of ENSO's global teleconnections it is indeed meaningful to provide a discrimination of La Niña periods into two flavors in close analogy to El Niño. Complementing recent findings on differential hydrometeorological conditions alongside potentially distinct La Niña flavors (Shinoda et al., 2011;Magee et al., 2017;Song et al., 2017;Hidayat et al., 2018), our results further demonstrate that seasonal wet/dry patterns accompanied by EP La Niñas are generally more likely to arise in a spatially coherent way than such observed along with CP periods (Figure 4). The same finding also applies to El Niño periods, which highlights that there exists a certain symmetry between both types of ENSO phases not only in the spatial SST anomaly patterns of El Niño and La Niña themselves, but also with respect to their effects on global precipitation patterns. Thus, from an impact oriented point of view, our work provides further evidence in favor of a distinction between two flavors of La Niña indicated by the presence or absence of very wet/dry regional conditions along with either of the two types of ENSO conditions. In other words, from the viewpoint of event based statistics (and therefore not necessarily in agreement with results based on linear correlations) it appears not only reasonable, but actually relevant to discriminate La Niña into two types in a similar way as common for El Niño periods.
In conclusion, our analysis provides a detailed and global overview on the large-scale differential imprints of different ENSO phases and flavors in the emergence of very wet/dry seasons. All findings reported in this work are ultimately summarized in Figure 5, which highlights schematically the main regions where the four different types of ENSO show large scale patterns of significant event coincidence rates with very high/low seasonal precipitation. Especially with respect to the CP flavors we find a variety of regions [such as (northern) Australia or southern Africa] where El Niño and La Niña show opposite impacts in terms of wet/dry seasons (Figures 5B,D). In addition, we observe that both flavors of one ENSO phase can also display similar impacts over the same regions, such as reduced precipitation over Australia or south-east Asia for EP and CP El Niños (Figures 5A,B) or enhanced precipitation over Australia or northern South America for EP and CP La Niñas (Figures 5C,D). We further uncover unique signatures that are only observed for a single type, e.g., reduced precipitation over northern India and the Middle East for EP (and not CP) La Niñas ( Figure 5C). The general trend toward less significant clusters of large-scale coincidences between ENSO and very wet/dry seasons along with CP periods is observable for both, El Niño and La Niña, as well (c.f. Figures 5A,B and Figures 5C,D). Our results thus demonstrate that a thorough discrimination of ENSO can be crucial for properly anticipating strong regional and seasonal hydrometeorological anomalies since its differential impacts may not only vary depending on the presence of El Niño and La Niña conditions but are additionally modulated by their respective EP and CP flavors.

PERSPECTIVES AND OUTLOOK
Future work should further apply the statistical concepts used in this work to also study ENSO related imprints on other climate variables (e.g., surface air temperature) as well as corresponding effects on socioeconomically and ecologically relevant observables like agricultural yields or water availability. This also implies that if the specific flavor of a developing El Niño or La Niña was to be detected early enough, possible threats like droughts or elevated flood risks, as well as their ecological and socio-economic consequences, could be better anticipated.
Since reliable predictions of EP or CP types of ENSO are scarce, the framework surrounding our present analysis has the potential to serve as an indirect early classification scheme of such periods based on observed globally distributed impacts. Specifically, the observation of early strong/weak climate signals either in SON or the presently ignored boreal summer (June-August, JJA) season of the onset year could be employed to assess the likelihood of approaching an ENSO phase with a specific flavor. Here, ECA could be used to analyze probabilities for specific types of seasonal climate anomalies-that could arise through a variety of conditions and drivers-to be preceded by certain ENSO periods and flavors (using so-called precursor event coincidence rates; Donges et al., 2016, see section 2.4). Such assessments would complement the estimation of trigger coincidence rates  at which a specific ENSO flavor is followed by a certain precipitation response as studied in our present work. The corresponding impact-based classification of ENSO periods could then be used to systematically estimate probabilities of potentially upcoming later (DJF, MAM, or even JJA of the year following the onset of an ENSO period) strong/weak precipitation seasons. In this context, we recall that for the case of instantaneous co-occurrences as studied here, precursor and trigger event coincidence rates solely differ by which of the two series (ENSO or precipitation events) is considered as the reference, i.e., the denominator in Equation (2). Since (sea or air) temperature anomalies are known to exhibit considerable persistence, it may therefore be another relevant extension of the present work to study non-instantaneous (i.e., multi-year) statistical linkages between ENSO episodes and wet/dry patterns, in which case trigger and precursor event coincidence rates become distinct statistical concepts.
Taken together, a corresponding analysis for early classification of ENSO flavors would need to be accompanied by a thorough review of recent approaches for ENSO prediction and attribution, a proper selection of climate variables beyond precipitation that potentially exhibit a mechanistic link with ENSO variability (Di Capua et al., 2019), and an associated estimation of appropriate temporal scales. Ultimately, such a predictive analysis would require a reliable validation of the obtained results by means of dividing the considered study period into training and test intervals, potentially accompanied by extending the analysis at least over the complete twentieth century or even beyond to ensure a sufficient number of events. Even though such an endeavor is clearly beyond the scope of the present work, it remains as a proposal for a potentially important subject of future research with the potential for an early attribution of EP and CP flavors to developing El Niños and La Niñas.
Finally, a more detailed intercomparison between the results obtained from traditional linear and complementary event based statistics could prove useful in assessing which regions are most affected in terms of strong climate responses to the presence of any combination of ENSO phase and flavor. Specifically, and in order to not only assess timing and regional extent of coincidences between ENSO events and precipitation signals, future work should systematically combine ECA with complementary techniques such as composite analysis to also study quantitatively the (relative) changes in magnitude of enhanced/reduced precipitation occurring alongside EP or CP flavours of ENSO. While the current work has solely focused on observational data, it should also be investigated if the reported co-occurrence patterns can be equally observed in historical simulations of state of the art coupled climate models [e.g., those contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP 6)], with any difference pointing to potentially insufficiently represented key processes in the models, thereby contributing to the identification of such processes and, hence, future model improvements. Ultimately, given that the frequency and magnitude of different ENSO phases might be markedly affected by global climate change (Yeh and Kirtman, 2007;Stevenson, 2012;Cai et al., 2014Cai et al., , 2015, it appears additionally useful to further apply the presented framework to future climate projections in order to assess possible changes in the spatial extent and frequency of ENSO related extreme events.

AUTHOR CONTRIBUTIONS
All authors designed the study. MW and JS analyzed the data and wrote the manuscript with input from all authors. JD assisted with the analysis. RD supervised the study. MW, JD, and RD substantively revised the work.