The variability of the hydroclimate over mainland Southeast Asia is strongly influenced by the El Niño–Southern Oscillation (ENSO), which has been linked to severe droughts and floods that profoundly influence human societies and ecosystems alike. Although the significance of ENSO is well understood, there are still limitations in the understanding of its effects on hydroclimate, particularly with regard to understanding the spatio-temporal characteristics and the long-term variation of its effects. Therefore we analysed the seasonal evolution and spatial variations in the effect of ENSO on precipitation over the period of 1980–2013 and the long-term variation in the ENSO teleconnection using tree-ring-derived Palmer drought severity indices (PDSIs) for the March–May season that span over the time period 1650–2004. The analyses provided an improved understanding of the seasonal evolution of the precipitation anomalies during ENSO events. The effects of ENSO were found to be most consistent and expressed over the largest areal extents during March–May of the year when the ENSO events decay. On a longer timescale, we found that ENSO has affected the region's March–May hydroclimate over the majority (95 %) of the 355-year study period and that during half (52 %) of the time ENSO caused a significant increase in hydroclimatic variability. The majority of the extremely wet and dry March–May seasons also occurred during ENSO events. However, considerable variability in ENSO's influence was revealed: the spatial pattern of precipitation anomalies varied between individual ENSO events, and the strength of ENSO's influence was found to vary through time. Given the high variability in ENSO teleconnection that we described and the limitations of the current understanding of the effects of ENSO, we suggest that the adaptation to ENSO-related extremes in hydroclimate over mainland Southeast Asia needs to recognise uncertainty as an inherent part of adaptation, must go beyond “predict and control”, and should seek adaptation opportunities widely within society.
Extremes or changes in the mean state of climate can result in great duress to societies, especially during periods of prolonged drought or flood. A well-known source of droughts and floods on a global scale is the ocean–atmosphere coupled phenomena El Niño–Southern Oscillation (ENSO) (Cane, 2005; Ward et al., 2014). ENSO is an evolving phenomenon (Trenberth and Shea, 1987), and it has become increasingly variable over recent decades (McGregor et al., 2013; Cai et al., 2014).
Over mainland Southeast Asia, henceforth MSEA, ENSO explains a large part of the inter-annual hydrological variability (Juneng and Tangang, 2005), and many of the recent severe droughts and floods occurred during ENSO events (see e.g. Räsänen and Kummu, 2013). Changes in hydroclimate variability is of great concern to the largely agrarian population of MSEA, as their livelihoods, economy, and food security are strongly dependent upon hydroclimatic conditions (MRC, 2010; Keskinen et al., 2010; ADB, 2016; Pech and Sunada, 2008). This dependency has triggered several studies that investigate the hydroclimatic variability and particularly the role of ENSO over MSEA.
Past research has shown that ENSO modulates precipitation and river flows over MSEA (Cook et al., 2012; Anchukaitis et al., 2016). Precipitation over MSEA is known to decrease during warm-phase (El Niño) events and increase during cool-phase (La Niña) events (Juneng and Tangang, 2005; Singhrattna et al., 2005b; Räsänen and Kummu, 2013; Kripalani and Kulkarni, 1997). The effects of El Niño on precipitation have been reported to evolve over Southeast Asia from south to north during development and decay phases of the events and particularly affect the southern parts of MSEA during the spring when the events decay (Räsänen and Kummu, 2013; Juneng and Tangang, 2005). ENSO's correlation with precipitation is known to be strongest in southern MSEA, weakening towards the north, and there are indications of opposite correlation between southern and northern areas (Kiem et al., 2005; Räsänen and Kummu, 2013; Zhang et al., 2007). These studies have contributed to the understanding on the effects of ENSO on precipitation over MSEA, and they provide either a high-resolution view over a smaller area or a coarse-resolution view over a larger area, but they do not provide a high-resolution view over the entire MSEA and its largest river basins, particularly on seasonal scales.
ENSO-related hydroclimatic anomalies over MSEA are known to vary through time. In general, during periods when hydrological conditions are below (above) average the effects of El Niño (La Niña) on precipitation are more severe (Kripalani and Kulkarni, 1997). Precipitation analyses over Thailand show that the connection between precipitation and ENSO has become stronger in the post-1980 period (Singhrattna et al., 2005b). Variation in the relationship between ENSO and hydroclimate is also found in the river flows. The analyses of the Mekong River show a stronger relationship between ENSO and river flow before the 1940s and after the late 1970s (Räsänen and Kummu, 2013; Darby et al., 2013). The changes in the relationship between ENSO and hydroclimate are linked to changes in ENSO's connection to different monsoon components. MSEA lies between the Indian summer monsoon (ISM) and western North Pacific summer monsoon (WNPSM) regions, and since the late 1970s the relationship between ENSO and WNPSM has strengthened while the relationship between ENSO and ISM has weakened (Wang et al., 2008; Hsu et al., 2014). These studies have shown temporal variations in the effects of ENSO in MSEA, albeit only over the last 100 years or so.
Xu et al. (2013) reconstructed the multivariate ENSO index using stable
isotopes of Oxygen (
Altogether, the body of research described above shows that the understanding of the linkage between ENSO and hydroclimate over MSEA has developed rapidly over recent years, but gaps exist, and there is a need to draw a more coherent picture. In this paper we focus on a research need consisting of combined analysis of three aspects: (1) high-spatial-resolution understanding of the seasonal evolution of correlation patterns between ENSO and precipitation, covering MSEA and its largest river basins; (2) spatial variation in precipitation anomaly patterns between individual ENSO events over MSEA; and (3) long-term temporal variation and stationarity of the ENSO teleconnection over MSEA. Advancing the knowledge of these three aspects would improve the scientific understanding of ENSO teleconnection and thus provide valuable information for adaptation to ENSO-related hydrological variability over MSEA and its largest river basins.
We therefore aim to analyse the instrumental and proxy records of hydroclimate over the region (Fig. 1). First we analyse instrumental records of precipitation over the period of 1980–2013 in order to investigate the seasonal evolution and spatial variation in the effect of ENSO on precipitation over MSEA, and second we analyse tree-ring-based proxy Palmer drought severity index data (PDSI; see Palmer, 1965) for the March–May season from two areas in MSEA that cover the time period of 1650–2004 to investigate the long-term variations in ENSO teleconnection. The methodology of using both precipitation and proxy PDSI data together aims to provide a more coherent view of the spatial and temporal variability in the effects of ENSO.
Map of the study area: mainland Southeast Asia. The spatial
variability of ENSO's influence was analysed using annual precipitation data
over the period of 1980–2013 with a focus on the area covering Myanmar,
Thailand, Lao PDR, Vietnam, and Cambodia and the largest river basins: the
Irrawaddy, Salween, Chao Phraya, Mekong, and Red rivers. The temporal
variability of ENSO's influence was analysed using proxy Palmer drought
severity index (PDSI) data for the March–May season over the period of
1650–2004 with a focus on two regions shown in the figure with rectangles
denoting the PDSI
The spatial and temporal analysis of ENSO's influence on hydroclimate is divided into two parts: analysis of seasonal precipitation over MSEA for the period 1980–2013, and analysis of proxy PDSI (for PDSI see Palmer, 1965) for the March–May season from two locations in MSEA over the period of 1650–2004.
The precipitation analysis aimed to improve our understanding of the spatial and temporal patterns of ENSO-related precipitation anomalies as well as our understanding of the strength of the relationship between ENSO and hydroclimate over the two proxy PDSI regions. The precipitation was analysed using precipitation data from the Global Precipitation and Climatology Centre (GPCC) (Schneider et al., 2015), the multivariate ENSO index (MEI) (Wolter and Timlin, 1993, 1998), and correlation analyses.
The analyses of proxy PDSI data aimed to improve the understanding of how the ENSO–hydroclimate teleconnection over MSEA has varied through time. Our analyses focus on the months of March–May, which span the transition period from dry to wet season, when the monsoon precipitation gradually starts (Adamson and Bird, 2010). March–May is also the beginning of the sowing season of rainfed rice in many areas (see e.g. Sawano et al., 2008), and the conditions of the early monsoon affect the transplanting of rice and thus the productivity of the crops (Fukai et al., 1998). The analyses of proxy PDSI data were based on two tree ring reconstructions from southern and northern Vietnam (Sano et al., 2008; Buckley et al., 2010; Cook et al., 2010), the unified ENSO proxy (McGregor et al., 2010), and correlation and wavelet methods (e.g. Torrence and Compo, 1998). In addition we analysed the co-occurrence of extremely dry and wet March–May seasons with ENSO events.
The seasonal precipitation analysis was based on GPCC v.7 data
(Schneider et al., 2015), which is an observation-based gridded
climatological data set with temporal coverage of 1901–2013 and spatial
resolution of 0.5
Description of the data sets used in the analyses of this study.
We also considered precipitation data from the Climate Research Unit (CRU TS v.3.21; Harris et al., 2014) and from Asian Precipitation – Highly Resolved Observational Data Integration Towards Evaluation (APHRODITE; Yatagai et al., 2009, 2012) precipitation data for the analyses, but comparisons suggested that GPCC v.7 was the most suitable. CRU TS v.3.21 had major gaps in stations in the region of Myanmar, and APHRODITE covers only a time period until 2007 and therefore does not capture the most recent influential ENSO events. The comparison of GPCC v.7 and APHRODITE over their common period provided very similar results.
First, the seasonal evolution of ENSO-related precipitation patterns was
analysed. ENSO events are generally 2-year phenomena that start to develop
in spring, mature late in the same year or early the next year, and decay in the
following summer. Therefore the precipitation was aggregated into JJA(0),
SON(0), DJF(0/1), MAM(1), JJA(1), and SON(1) seasonal sums and correlated with
the time series of the January–February–March value of MEI
(NOAA, 2015) from the second year of each ENSO event
(MEI
Second, we analysed the seasonal precipitation anomalies for each ENSO event and for each season over MSEA. Anomalies were calculated as deviations from the 1980–2013 average precipitation and reported as percentages. This yielded seasonal precipitation anomaly maps of all El Niño and La Niña events for the period of 1980–2013. In addition, we analysed the precipitation anomalies in more detail at the locations of proxy PDSI data in order to understand how strongly the hydroclimate at those locations is related to ENSO. This helps to assess how well the PDSI proxies are suited for analysing long-term ENSO teleconnection.
The temporal variability of ENSO's teleconnection to MSEA was analysed using
two tree-ring-based PDSI reconstructions developed by Sano et al. (2008)
and Buckley et al. (2010), for northern and southern Vietnam,
respectively. These two reconstructions marked the first two successful
calibration-verification model schemes from tropical tree rings, both from
the long-lived Vietnamese cypress (
We used the unified ENSO proxy (UEP), an index based on the 10 most
commonly used ENSO proxies that was published by McGregor et
al. (2010), to describe ENSO behaviour over the period 1650–2005. The
original UEP is annual data and covers the time period from 1650 to 1977. We
extended the UEP up to the year 2004 by using MEI in order to match the time
period of the PDSI data. To do so, we scaled the UEP variance to match the
variance of MEI (UEP
The PDSI
The applied wavelet methods included the computation of the wavelet power spectrum of a single time series, as well as the cross-wavelet power spectrum and wavelet coherence spectrum of two time series together. The computations were done using the WaveletComp R package developed by Rösch and Schmidbauer (2014). The wavelet power spectrum shows the time series in time-frequency space, which allows the examination of variations and their power with respect to their frequency and occurrence in time, while the cross-wavelet power spectrum shows where the variations of two time series have high common power in the time-frequency space. The wavelet coherence spectrum shows the coherence (i.e. localised correlation) between the two time series in time-frequency space, while the cross-wavelet power spectrum and the wavelet coherence spectrum also show the phase relationship between the two time series. In the case of correlated phenomena, the phase relationship is expected to be consistent in time. A more complete treatment of the wavelet methods can be found in Torrence and Compo (1998) and Grinsted et al. (2004).
The wavelet methods were used to identify temporal variability in the strength of ENSO's influence on hydroclimate over MSEA. Two categories were used for this identification: (i) primary ENSO-related variance and (ii) secondary ENSO-related variance in the hydroclimate of MSEA. These periods were defined according to regions in wavelet power, cross-wavelet power, and coherence spectrum that were overlapping in time-frequency space and fulfilled specific criteria. The specific criteria are explained in detail in Table 2. The major difference between the two categories is that in the former the increase of the wavelet power is statistically significant. Non-significant ENSO-related increases in wavelet power are also analysed as they reveal periods that still do have statistical relationship between ENSO and hydroclimate and provide an indication of the variations in the strength of ENSO teleconnection over MSEA. The wavelet analyses focused on periodicities from 2 to 10 years as they represent the frequencies of inter-annual ENSO variability. The statistical significance of the wavelet power and coherency was tested against white noise at the 5 % significance level.
The identification criteria for periods with ENSO-related variance in March–May hydroclimate. Two types of variance periods were identified from unified ENSO proxy and PDSI proxy data: primary ENSO-related variance and secondary ENSO-related variance in the hydroclimate. These periods were defined according to regions in wavelet power spectrum (WP), cross-wavelet power (CWP), and coherence spectrum (WC) that were overlapping in time-frequency space and fulfilled the criteria in the table. Variance period refers to period when ENSO had increased influence on the March–May hydroclimate in mainland Southeast Asia.
In addition to wavelet analysis, we employed a variance analysis of the PDSI with an 11-year moving window in order to identify periods with high inter-annual variability in the time domain. This process also enabled us to see how well these periods correspond with the high-variability periods identified from wavelet analysis. We chose 11 years in order to capture the band of inter-annual variability without the decadal variability.
The co-occurrence of extremely dry and wet years with ENSO events was based on the Gergis and Fowler (2009) multi-proxy ENSO event reconstruction over the period of 1525–2002. The extreme years were defined from PDSI data using 5th and 95th percentiles, which meant that 10 % of all years of PDSI data were defined as extreme. The co-occurrence of extreme years with warm- and cool-phase ENSO events was then identified by comparing the multi-proxy ENSO event reconstruction and extreme PDSI values.
The seasonal correlation analysis of precipitation and MEI
Map of correlation of January–February–March values of
multivariate ENSO index (MEI
Taking a closer look at these patterns, the negative correlations are observed during SON(0) in the southern coastal regions of MSEA in the west in Thailand and Myanmar and in the east in southern Vietnam and Cambodia. In DJF(0/1) the areas of negative correlation are pushed further south by areas of positive correlations. In MAM(1) the negative correlations are widespread and cover most of the study area, except northern Myanmar and parts of China. In JJA(1) the areas of negative correlations are observed mainly in western Thailand and in southern Myanmar, and in SON(1) the negative correlations have more or less disappeared. Another interesting feature is the statistically significant positive correlation during the JJA(0) season in southern Myanmar, southern Lao PDR, and northern Cambodia, separated by an area of negative correlation in Thailand.
The analysis of precipitation anomalies shows spatially varying anomaly patterns between ENSO events. This can be observed in Fig. 3, which shows the MAM(1) precipitation anomalies of eight El Niño and four La Niña events during the period of 1980–2013 (see also Figs. S1 and S2 in the Supplement, which show precipitation anomalies for all seasons for the same El Niño and La Niña events as in Fig. 3).
March–May precipitation anomalies (%) during the second year
(MAM(1)) of
In the case of the El Niño events of 1982–1983, 1986–1987, 1991–1992, 1994–1995, 1997–1998, and 2009–2010 (Fig. 3a–h), the MAM(1) precipitation anomalies are widely negative in large parts of the study area. During the El Niño event of 2002–2003 the negative precipitation anomalies are smaller in magnitude, and positive anomalies are observed in some regions, for example in southern Myanmar and at the border between southern Lao PDR and western Thailand. During the El Niño event of 2006–2007 the precipitation anomalies are mainly positive and thus inconsistent with other El Niño events.
In the case of La Niña events there is greater inconsistency in spatial patterns of MAM(1) precipitation anomalies than in the case of El Niño. During the 1998–1999 La Niña event, the MAM(1) precipitation anomalies are largely positive and cover Thailand, Cambodia, southern Lao PDR, southern Vietnam, and large parts of Myanmar. During the 1988–1989 event, the positive precipitation anomalies are confined to the eastern part of the study area in Vietnam, in 2007–2008 the precipitation anomalies are smaller but more widespread and can be seen particularly in Cambodia and eastern Thailand, and in the 2010–2011 event the positive precipitation anomalies are mainly in Myanmar and in western Thailand.
The time series analysis of MAM(1) precipitation for the areas of
PDSI
ENSO events (NOAA, 2015) and March–May precipitation
anomalies in the areas of PDSI
During MAM(1
The precipitation analyses provided a good understanding of the hydroclimate
and its relationship to ENSO in the areas of PDSI
The correlation analysis between ENSO
Correlations between ENSO
Periods with evidence of ENSO teleconnection in March–May
hydroclimate of mainland Southeast Asia over the period of 1650–2004.
“Correlation periods” refer to periods with statistically significant
correlation in moving-window correlation analysis (Fig. 4), and
“periods with primary and secondary ENSO-related variance in
hydroclimate” refer to periods when ENSO had stronger influence on
hydroclimate according to wavelet analyses (Figs. 5–6). Statistically
significant periods (
The wavelet analyses show a connection between ENSO and the hydroclimate of
the region throughout the study period (Figs. 5–6). The connection can be
observed as a relatively consistent temporal distribution of statistically
significant areas in the wavelet coherence spectrum of ENSO
Wavelet analysis of ENSO and PDSI
Wavelet analysis of ENSO and PDSI
The wavelet analyses of PDSI
The wavelet analyses of PDSI
Periods with evidence of ENSO-related hydrological variability in
March–May hydroclimate of mainland Southeast Asia over the period of
1650–2004. The periods with statistically significant correlation between
the time series of ENSO
The wavelet analyses also reveal that increased variance in ENSO does not
always result in increased hydroclimatic variance over MSEA. For example,
the statistically significant increases in wavelet power of ENSO
The analysis of extreme PDSI values in Figs. 5e and 6e shows that the
majority of the most-extreme dry and wet MAM seasons occurred during ENSO
events, particularly in the region of PDSI
When the results of the moving-window correlation analyses and the wavelet
analyses of both PDSI
The research approach used, based on a combination of precipitation data and proxy PDSI data derived from tree ring records, provides a more uniform and coherent picture of the spatio-temporal effects of ENSO over MSEA and its largest river basins. The analysis of precipitation data showed how the precipitation anomalies evolve in time during ENSO events and how they vary in space between individual ENSO events. The analysis of proxy PDSI data in turn showed not only how the effects of ENSO have varied for the monsoon transition period (March–May) over a longer timescale but also how the effects have varied spatially between northern and southern areas of MSEA. In the following sections we further discuss important aspects of the methodology, state our contributions and compare our findings with past research, and suggest directions for future work as well as for adaptation to ENSO-related hydroclimatic anomalies.
The analysis of the long-term ENSO-hydroclimate relationship using two
methods (moving-window correlation and wavelets) and two hydroclimate
proxies derived from tree rings (PDSI
However, there are certain limitations in the above approach in providing
exact years for the periods with connection between ENSO and hydroclimate.
First, the proxy PDSI analyses focused only on the MAM season. However, this
season is deemed to be appropriate for detecting an ENSO signal over MSEA,
as our analyses revealed that the correlation between ENSO and precipitation
over MSEA was strongest and statistically significant over the largest area
of MSEA during the MAM season compared to other seasons (Fig. 2). Moreover,
the proxy PDSI data are most accurate for the MAM season: the tree ring data
have the strongest correlation with instrumental PDSI data and provide the best
verification results for the MAM season (Sano et al., 2008; Buckley et
al., 2010). As argued in the Methodology section, the MAM season is also
hydrologically important. For example, our analyses showed that, in the area
of PDSI
Further, the moving-window correlation was based on a window size of 21 years, resulting in ambiguity in the dating of the statistically significant periods. Third, the visual interpretation of the wavelet images involves a certain amount of subjectivity when multiple images are compared simultaneously. For example, subjective judgement was needed when the statistically significant areas in wavelet power, cross-wavelet power, and coherence spectrum images were of different size and not perfectly overlapping and when the phase arrows varied slightly from the expected direction. In order to minimise the errors from subjectivity, clear rules for consistent interpretation were developed and followed (see Sect. 2). Fourth, the size of statistically significant areas in wavelet images depended on parameters of the wavelet analysis. For example, the choice of statistical significance testing method affected the size of the statistically significant areas, which may change the timing and duration of any such identified ENSO periods with so few years. Last, it is likely that the approach used was not able to capture all individual ENSO events that resulted in hydroclimate anomalies. Despite these limitations, the results are based upon standard methods in time series analysis and are therefore considered to be reliable estimates of ENSO-related hydrological variability.
The past research provides a view on the general influence of ENSO on precipitation over MSEA, as discussed in the Introduction section. The El Niño (La Niña) events result in drier (wetter) conditions over MSEA (Juneng and Tangang, 2005; Kripalani and Kulkarni, 1997) and at more local scales in Thailand (Singhrattna et al., 2005b) and in the Mekong River basin (Räsänen and Kummu, 2013). These studies also suggest stronger correlation between ENSO and hydroclimate in central and southern parts of MSEA. A transition of the influence of ENSO to opposite sign from south to north is also reported for the Mekong River basin (Räsänen and Kummu, 2013) and is supported by studies focusing on the upper reaches of the Mekong and Yangtze River basins (Kiem et al., 2005; Zhang et al., 2007). The precipitation anomalies are also shown to evolve northeastward during El Niño events from southern parts of Southeast Asia to MSEA (Juneng and Tangang, 2005).
The current research confirms these past findings on the effects of ENSO on precipitation and expands the existing knowledge in three aspects: first, by providing more detailed and informative description of the seasonal evolution of the effects of ENSO in MSEA and by showing this evolution in more northern areas (compared to Juneng and Tangang, 2005) (Fig. 2); second, by showing the areas and the season when the transition of the influence of ENSO to opposite sign occurs (Fig. 2c); and third, by showing how the spatial patterns of precipitation anomalies have varied between individual ENSO events over MSEA (Fig. 3, Figs. S1 and S2). Through these contributions the current research provides a more accurate and uniform picture of the spatio-temporal effects of ENSO on precipitation and thus allows a more detailed comparison of effects of ENSO between different regions and seasons of MSEA and its largest river basins.
The long-term relationship of ENSO and hydroclimate in MSEA has been shown to exist at centennial scales by several studies (Xu et al., 2013; Buckley et al., 2007, 2010; Sano et al., 2008), but the variation of the relationship of ENSO and hydroclimate has been studied only over the past 100 years or so. Studies conducted in Thailand (Singhrattna et al., 2005b) and the Mekong River basin (Räsänen and Kummu, 2013; Räsänen et al., 2013; Darby et al., 2013) report that the most recent periods of stronger relationship between ENSO and hydroclimate occurred during the beginning of the 20th century and lasted until the 1940s, while the second period began around the 1960s–1980s.
The current research agrees with the findings from Thailand and the Mekong
River basin and suggests a period of weaker relationship between ENSO and
hydroclimate during the 1930s–1950s in the southern parts of the study area
(PDSI
It is worthwhile to further highlight that the article's demonstration of the strong inverse relationship between the reconstructed drought metric PDSI and ENSO fit within a broader context of studies demonstrating the importance of ENSO. The tree ring studies used here (Sano et al., 2008; Buckley et al., 2010) focus on northern and southern Vietnam. Buckley et al. (2014) expand upon this discussion by using tree ring records from all across monsoon Asia and North America, illustrating that the dominant mode of climate variability across both sides of the Pacific is driven by ENSO-like variability, particularly at decadal scales (i.e. the Interdecadal Pacific Oscillation, or IPO – see Meehl and Hu (2006) and Buckley et al. (2010) for further details). Indeed, other tree ring sites from Thailand (Buckley et al., 2007) and Myanmar (D'Arrigo et al., 2011) confirm the strength of this relationship in these regions as well.
The findings of the current paper indicate considerable uncertainties in the
effects of ENSO on hydroclimate and how this relationship develops through
time. For example, clear patterns were found in the seasonal evolution of
precipitation anomalies during ENSO events, but at the same time the
precipitation anomalies and their spatio-temporal patterns were found to
vary considerably between ENSO events. This leads to two potentially useful
research directions related to ENSO. The first research direction would
explore the physical characteristics (e.g. sea surface temperature, air
pressure, wind, and moisture fluxes patterns) of each ENSO event and how they
translate into anomalies of MSEA hydroclimate. For example, it is
hypothesised that the placement of the descending limb of the Walker
circulation could affect the ENSO teleconnection over MSEA
(Singhrattna et al., 2005b). The second direction could be to
explore how other climatic and oceanic phenomena interact with the ENSO
teleconnection over MSEA. For example, it is known that the Indian Ocean
Dipole (IOD) affects the hydroclimate over MSEA (Darby et al., 2013), and
there are good indications of the effect of the Pacific Decadal Oscillation
(PDO) (Delgado et al., 2012). In addition, the current
research discovered statistically significant positive correlation between
precipitation and MEI
The findings of this study also provide perspectives for adaptation to extremes in hydroclimate. The findings suggest some degree of statistical predictability of ENSO-related anomalies in hydroclimate, but at the same time the findings revealed large variation and thus uncertainties in the effects of ENSO over MSEA. It is well known that statistical approaches can have limitations when it comes to predicting extreme events (see e.g. Nassim, 2010). Thus, given the high variability in the effects of ENSO, limitations in the current knowledge, and statistical approaches we suggest exploration of adaptation approaches that embrace uncertainty and complexity and seek adaptation opportunities in multiple sectors and levels of society (see e.g. Resilience concept: Walker et al., 2004, 2013) while considering ongoing anthropogenic environmental changes (Keskinen et al., 2010; Lauri et al., 2012; Pech and Sunada, 2008). For example, adaptation only through engineering solutions is likely to aggravate already-existing challenges (e.g. Baran and Myschowoda, 2009). The suggested adaptation approaches could further benefit from analysis of the societal impacts of the identified historical events and the coping mechanisms used to deal with them in the past (Nuorteva et al., 2010; Buckley et al., 2010).
Hydroclimate variability affects various economic activities, local livelihoods, and food security across MSEA. This research sought to improve our understanding of the hydroclimate variability by investigating the spatial and temporal variability of MSEA's ENSO over the period of 1650–2013. The investigations were based on analyses of gridded seasonal precipitation data (1980–2013), proxy Palmer drought severity index values for the March–May season, and proxy ENSO data (1650–2004).
These analyses provide a more accurate and uniform picture of the spatio-temporal effects of ENSO on precipitation and improve our understanding of the long-term ENSO teleconnection and its variability over MSEA. The research reveals new information on the seasonal evolution of the effects of ENSO over MSEA, and it shows how the spatial patterns of the effects of ENSO vary between individual events. On a longer timescale, the strength of the effects of ENSO on hydroclimate of the March–April season (the important monsoon transition season with most-widespread ENSO effects in MSEA) was shown to vary between periods of weaker and stronger effects. Altogether our findings reinforce the significance of ENSO over MSEA, but they also expand the past knowledge by describing the high degree of variability and non-stationarity in the effects of ENSO. This described variability implies challenges for understanding and predicting the effects of ENSO over MSEA into the future.
Given the high impact and variability of ENSO, and limitations in the current knowledge and predictive skill, adaptive approaches for mitigating the negative effects of ENSO are recommended. Adaptation should embrace uncertainty, seek adaptation opportunities within multiple sectors and levels of society, and consider climate-related adaptation as part of broader adaptation to ongoing social and environmental changes. Forecasting- and engineering-based approaches alone are likely to be inadequate and will risk creating further challenges.
All data used in this study is publicly available. The precipitation data
(GPCC v.7; Schneider et al., 2015) are is available at Deutscher Wetterdiest:
Timo A. Räsänen and Ville Lindgren received funding from Maa- ja vesitekniikan tuki ry, Joseph H. A. Guillaume from the Academy of Finland-funded project NexusAsia (grant no. 269901), Brendan. M. Buckley from NSF grants GEO 09-08971 and AGS 130-3976, and Matti Kummu from the Academy of Finland-funded project SCART (267463) and the Emil Aaltonen Foundation-funded project “eat-less-water”. Lamont contribution number 8055. Edited by: D. Fleitmann Reviewed by: two anonymous referees