Interdecadal Variability in Myanmar Rainfall in the Monsoon Season (May–October) Using Eigen Methods

: In this study, we investigated the interdecadal variability in monsoon rainfall in the Myanmar region. The gauge-based gridded rainfall dataset of the Global Precipitation Climatology Centre (GPCC) and Climatic Research Unit version TS4.0 (CRU TS4.0) were used (1950–2019) to investigate the interdecadal variability in summer monsoon rainfall using empirical orthogonal function (EOF), singular value decomposition (SVD), and correlation approaches. The results reveal relatively negative rainfall anomalies during the 1980s, 1990s, and 2000s, whereas strong positive rainfall anomalies were identiﬁed for the 1970s and 2010s. The dominant spatial variability mode showed a dipole pattern with a total variance of 47%. The power spectra of the principal component (PC) from EOF revealed a signiﬁcant peak during decadal timescales (20–30 years). The Myanmar summer monsoon rainfall positively correlated with Atlantic multidecadal oscillation (AMO) and negatively correlated with Paciﬁc decadal oscillation (PDO). The results reveal that extreme monsoon rainfall (ﬂood) events occurred during the negative phase of the PDO and below-average rainfall (drought) occurred during the positive phase of the PDO. The cold phase (warm phase) of AMO was generally associated with negative (positive) decadal monsoon rainfall. The ﬁrst SVD mode indicated the Myanmar rainfall pattern associated with the cold and warm phase of the PDO and AMO, suggesting that enhanced rainfall for about 53% of the square covariance fraction was related to heavy rain over the study region except for the central and eastern parts. The second SVD mode demonstrated warm sea surface temperature (SST) in the eastern equatorial Paciﬁc (El Niño pattern) and cold SST in the North Atlantic Ocean, implying a rainfall deﬁcit of about 33% of the square covariance fraction, which could be associated with dry El Niño conditions (drought). The third SVD revealed that cold SSTs in the central and eastern equatorial Paciﬁc (La Niña pattern) caused enhance rainfall with a 6.7% square covariance fraction related to ﬂood conditions. Thus, the extra-subtropical phenomena may affect the average summer monsoon trends over Myanmar by enhancing the cross-equatorial moisture trajectories into the North Atlantic Ocean.


Introduction
Rainfall variability is among the basic indicators of climate and water cycle changes in a region [1]. Anomalous changes in rainfall dynamics may cause hydro-meteorological hazards such as flood, drought, and storms [2,3], thus ultimately resulting in loss of human

Data
We used monthly gauge-based gridded rainfall datasets of the Global Rainfall Climatology Centre (GPCC) and Climatic Research Unit version TS4.0 (CRU TS4.0) during the period of 1950-2019 [39][40][41][42]. The gridded climate datasets were collected with a spatial resolution of 0.5° × 0.5° over the study region. Several institutes have developed gridded climate datasets during the past few decades, which are frequently used in hydro-climatic assessments [40,43,44]. Among others, the data products of CRUTS4.0 and GPCC are commonly used due to a longer temporal span [39,41]. The in situ observation data of 35 meteorological stations were acquired from the Department of Meteorology and Hydrology (DMH), Myanmar. The long-term in situ observation data were further used after validation with observed summer monsoon rainfall (1950-2019), which is highly correlated with GPCC gridded data, to investigate the interdecadal variability.
The monthly data of the Pacific decadal oscillation (PDO) Index from May-October for the period 1950 to 2010 were used in this study [22,45]. The AMO index is defined as a local estimate over the entire North Atlantic low-resolution annual SST anomalies after removing any direct trend [46]. The North Pacific (NP) index, North Atlantic Oscillation (NAO) index, and the Arctic Oscillation (AO) index data were also used in this study [47].

Data
We used monthly gauge-based gridded rainfall datasets of the Global Rainfall Climatology Centre (GPCC) and Climatic Research Unit version TS4.0 (CRU TS4.0) during the period of 1950-2019 [39][40][41][42]. The gridded climate datasets were collected with a spatial resolution of 0.5 • × 0.5 • over the study region. Several institutes have developed gridded climate datasets during the past few decades, which are frequently used in hydro-climatic assessments [40,43,44]. Among others, the data products of CRUTS4.0 and GPCC are commonly used due to a longer temporal span [39,41]. The in situ observation data of 35 meteorological stations were acquired from the Department of Meteorology and Hydrology (DMH), Myanmar. The long-term in situ observation data were further used after validation with observed summer monsoon rainfall (1950-2019), which is highly correlated with GPCC gridded data, to investigate the interdecadal variability.
The monthly data of the Pacific decadal oscillation (PDO) Index from May-October for the period 1950 to 2010 were used in this study [22,45]. The AMO index is defined as a local estimate over the entire North Atlantic low-resolution annual SST anomalies after removing any direct trend [46]. The North Pacific (NP) index, North Atlantic Oscillation (NAO) index, and the Arctic Oscillation (AO) index data were also used in this study [47]. Large-scale atmospheric circulation, such as air temperature and relative humidity, was collected from ERA5 the 5 th generation reanalysis of the European Centre for Medium-Range Weather Forecasts (ECMWF) to determine the related breakouts and their moisture transport [48,49]. The extended version of the Reproduced Ocean Surface Temperature (ERSST, version3b), published in the National Maritime and Barometrical Administration/National Climatic Information Center, was used for seawater heating [50]. Analytical data from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) [51], such as speed, zonal and meridional wind, vertical speed (0.5 • × 0.5 • ), were used to explore long-term wind zones.

Methods
The EOF is a statistical method broadly used to minimize the multidimensionality of complex climate information and recognize the most imperative physical modes with the least chance of data being misplaced [52,53]. Principal component analysis (PCA) is used for the assessment of atmospheric information based on the EOF method [54]. The technique describes the variance and covariance of the data, describing a few modes of variability. The modes that explain the largest percentage of the original variability are considered significant. The modes can be represented by orthogonal spatial patterns (eigenvectors) and corresponding time series (principal components). The EOF technique has been successfully used in Myanmar to investigate the interannual variability in summer monsoon rainfall [25]. With regard to the temporal behavior of the EOF mode, a 10-year running average was applied to obtain the long-term variability. The summer monsoon rainfall data were normalized to prevent areas and seasons of maximum variance from dominating the eigenvectors [55]. The standardized rainfall values were computed for all the years from the long-term mean, yearly mean, and the standard deviation using Equation (1): where Z represents the standardized departure, X is the long-term mean value, and S d is the standard deviation from the mean. The Z value provides immediate information about the significance of a specific deviation from the mean [25,56]. The orthogonal function of EOF is defined as: where z(x, y, t) denotes the function of space (x, y) and time (t), while EOF(x, y) represents the spatial structure in relation to the temporal variation in Z.
A power spectrum was analyzed to evaluate the multiscale temporal variations resulting from the wavelet transform employed in the long-term rainfall dataset [57][58][59]. The decadal spectra were used to estimate the changes that occur in the rainfall pattern on a decadal basis. Prior to performing the power spectrum analysis from the results of the EOF, PCA was used to verify the significance of the scales of interdecadal variation in this study. The red noise spectrum of the decadal spectrum was computed [54]. The confidence limits related to the read noise spectrum were calculated with 90% confidence levels.
Composites include classifying and averaging one or more category of the variable fields according to their relation with the prevailing situations [60]. Composites results are used to produce hypotheses for patterns related to individual scenarios [53,61]. The goal of the composite analysis was to produce the wet and dry years separately. It is primarily used in China and Myanmar to detect the circulation anomalies connected with wet/dry events [25,45,62].
Horizontal moisture transfer rate/flux convergence (MFC), known as moisture convergence, was identified for wet and dry years using the vector character. The MFC results Water 2021, 13, 729 6 of 27 from the conservation of water vapor in pressure coordinates, which has been discussed in detail previous [63]. MFC can be calculated as: where ∇ =î ∂ ∂x +ĵ ∂ ∂y , q is the specific humidity, and V h = (u, v). The nonparametric Mann-Kendall test (MK) test was used to detect the trend in monsoon rainfall [64,65]. The MK test is widely used to evaluate trends in rainfall, agrometeorological, and hydrological time series [66][67][68][69][70][71][72]. This trend test's limitations are normally associated with its own null hypothesis, which assumes that the data are independently and identically distributed. However, the null hypothesis is commonly taken as evidence of the null hypothesis, which is commonly used as evidence of a trend in a given agrometeorological and/or hydrological time series. The MK test investigates the sequential change point to emphasize an abrupt change significant at the 95% confidence level. This approach has been used successfully worldwide [70,73,74] to detect climate variables' change points. For certain datasets consisting of x values with a sample size Sen slope (SS), the MK calculation begins by estimating the S statistic as: As presented in Mann (1945) and Kendall (1975), when SS ≥ 8, the distribution of S approaches the Gaussian form with mean E(S) = o and variance V (S) is calculated using: where ti is the number of ties of length m. The S statistic is standardized Z, as shown in Equation (6), and its significance can be estimated from the normal cumulative distribution function.
Singular value decomposition (SVD) analysis is used in two techniques: meteorological and oceanographic data analysis. SVD is applied between two identical datasets of two jointly analyzed fields to identify pairs of the coupled spatiotemporal variations. Bretherton et al., (1992) [75] and Wallace et al., (1992) [76] provided a brief introduction to the SVD method as a fundamental matrix operation that can be considered an extension of rectangular matrices of the diagonalization of a square symmetric matrix. In the calculation, each pair defines the covariance percentage between the two joint fields, allowing the extraction of the dominant modes of pair covariability between the two analyzed fields. Moreover, the SVD of the cross-covariance matrix identifies from two data fields pairs of spatial patterns that explain as much as possible of the mean-squared temporal variance between the two fields. This method is widely used [25] because it works on common and ambiguous data sets. In this work, SVD was used to find the maximum covariance or correlation between the summer monsoon rainfall over Myanmar and the SST over the Indian Ocean and Pacific Ocean.
After this, we used ensemble empirical mode decomposition (EEMD), which is an adaptive time-frequency data analysis technique developed by Wu et al., (2007) [77] and Wu and Huang (2009) [78]. It is used to resolve problems of mixed-mode and false components in the empirical mode decomposition (EMD) decomposition process. This method is a major modification of the original EMD method, emphasizing the adaptiveness and temporal locality of the data decomposition [79]. The EMD method can decompose any complicated data series into a few amplitudes of frequency-modulated oscillatory components, called intrinsic mode functions (IMFs), of distinct timescales [80]. EMD is based on the local characteristic time scale of a signal, and decomposes the complicated signal into a number of IMFs. The main steps of the EEMD method employed in this study were similar to those outlined by Qian (2009) [81,82]. In this study, EEMD was used to decompose the summer monsoon rainfall and PDO index. The climatological annual cycle of the monthly mean monsoon rainfall over Myanmar (1950-2019) was estimated based on: in situ observation, CRUTS4.0, and GPCC datasets; the results are shown in Figure 2. The error bars shows the monthly standard deviation of climatological annual cycle with in each active and break spells. The horizontal solid lines show the monthly standard deviations of all datasets with blue, grey, and orange bars. Interestingly, the significant number of active spells (blue and orange bars) show positive rainfall pattern from June to August whereas the break spells show relatively weak rainfall pattern during September to October. Table 1 reveals the amount of average monthly summer monsoon rainfall (May-October) from in situ observations, CRUTS4.0, and GPCC for the period of 1950-2019. The country experienced a relatively high unimodal rainfall pattern in June, July, and August, observed in all datasets. The average monthly rainfall from in situ observations was 480 mm in June, 520 mm in July, and 524 mm in August. From GPCC, these amounts were 478, 519, and 519 mm, respectively. The rainfall from CRUTS4.0 exhibited a weaker pattern in these three months (444, 466, and 476 mm, respectively). The higher rainfall agrees with the northern hemisphere summer and monsoon season. The rainfall season in Myanmar extends from May to October, and the rainfall season coincides with summer monsoon onset in the country, which generally occurs in the third week of May. Similar monsoon onset periods were observed by  in a simulation study using Regional Climate Model (RegCM3). Notably, the GPCC dataset produced more accurate results compared to the in situ observations. CRUTS4.0 reported less rainfall compared to the in situ observations and GPCC. We found a strong correlation (0.94) between the in situ observations and GPCC datasets, but a relatively weak correlation (0.52) between the in situ observations and the CRUTS4.0 dataset ( Table 2). The root mean square error (RMSE) for GPCC (14.24) for summer rainfall was lower than that of the CRUTS4.0 dataset (45.29) (Table 2). Therefore, GPCC was thus used in most analyses since the data are continuous and can be extrapolated to cover all parts of the country to identify longer time periods. Roy et al., (2011) studied the interannual variability in summer monsoon rainfall in relation to Indian Ocean Dipole (IOD) and ENSO using the GPCC dataset. Therefore, the use of GPCC data provided better results, due to extended longer periods, which were more complete as compared to the observed data.      Both datasets showed relatively similar spatial rainfall patterns over the country (Figure 3). Maximal rainfall was received along the western coast (Rakhine State), and southern (Mon, Kayin States and Tanintharyi Region) and northern tips (Kachin State) of Myanmar. This aligns with the recent extreme events (floods, strong winds, and landslides) that have occur in these areas due to the strong monsoon rainfall. Along the coastal area (west and south), regions close to the BoB, Andaman Sea, and northern tip (Kachin State) occupy the highest mountain, called Hkakabo Razi. The central (dry zone) and eastern areas receive less rainfall than other places in the region. The central dry zone receives the lowest amount of rainfall compared to rest of the study area ( Figure 3). The observed low rainfall in the central dry zone is attributed to the topography: the area is situated between two mountain ranges: the Shan Plateau to the east and Rakhine Yoma Mountain to the west. Kumar et al., (2006) indicated that increased rains on the Myanmar coast contribute more than total rainfall in the West Ghats, associated to the central and northern regions.
Water 2021, 13, x FOR PEER REVIEW 9 of 28 with the recent extreme events (floods, strong winds, and landslides) that have occur in these areas due to the strong monsoon rainfall. Along the coastal area (west and south), regions close to the BoB, Andaman Sea, and northern tip (Kachin State) occupy the highest mountain, called Hkakabo Razi. The central (dry zone) and eastern areas receive less rainfall than other places in the region. The central dry zone receives the lowest amount of rainfall compared to rest of the study area ( Figure 3). The observed low rainfall in the central dry zone is attributed to the topography: the area is situated between two mountain ranges: the Shan Plateau to the east and Rakhine Yoma Mountain to the west. Kumar et al. (2006) indicated that increased rains on the Myanmar coast contribute more than total rainfall in the West Ghats, associated to the central and northern regions.

Spatial and Temporal Interdecadal Variability in Summer Monsoon Rainfall.
The change in decadal mean May-October rainfall from the long-term mean rainfall over Myanmar was estimated using GPCC data ( Figure 4). In the six decades, the results show strong negative anomalies during the periods of 1971-1980, 1981-1990, and 1991-2000
Water 2021, 13, x FOR PEER REVIEW 10 of 28 reported that the occurrence of extreme events (i.e., drought and flood) is significantly associated with large-scale atmospheric circulations of the Pacific and Indian Oceans, particularly during the positive phase of the ENSO. The overall strong and positive decadal rainfall anomalies occurred in La Niña years over the target region. Table 3 reveals that the observed extreme rainfall between 1950 and 2019 occurred in the first three decades for the west coast (Rakhine State), for two decades along the south coast (Tanintharyi Region), and for one decade in the Upper Sagaing Region. This result agrees with the spatial decadal monsoon rainfall in Figure 4.  From Figure 4 and Table 3, the decadal rainfall fluctuated during the period of 1950-2019, with a positive skewness and increased frequency in recent decades. This suggests an increase in the extreme rainfall events with increased frequency in all study regions, especially in the western and southern regions. This may be associated with increased  From Figure 4 and Table 3, the decadal rainfall fluctuated during the period of 1950-2019, with a positive skewness and increased frequency in recent decades. This suggests an increase in the extreme rainfall events with increased frequency in all study regions, especially in the western and southern regions. This may be associated with increased global-warming-induced changes in monsoon and circulation patterns in the regions; similar rainfall patterns were previously reported [83].
The results for various regions in Myanmar revealed that the west, south, and north regions received the most rainfall, while the east and central regions received relatively less ( Figure 5). Figure 5 shows that an overall increasing trend in rainfall occurred over all regions, which corresponds to previously reported findings [18,84,85]; however, the annual rainfall variability on the regional scale differed. In the northern part, the rainfall variability appeared to be increasing, with an amount of 500 mm. The peak rainfall years (dry years) were recorded during 1992, 1997, and 2010 (1982, 1991, 2005, and 2009). In the eastern part, monsoon rainfall slightly increased, as depicted by a rainfall amount of 250 mm with dry and wet years appearing to be the possible reason for balancing the overall trend in the region. In this region, the highest (lowest) monsoon rainfall was observed in 1988 and 2013 (1987 and 2009). In the central part, the monsoon rainfall exhibited an increasing trend at the rate of 300 mm per year. The analysis showed that 1960, 1997, 2010, and 2014 (1987, 1991, 1999 and 2009) were the peak wet (dry) years in the northern and central parts. Our results concur with the findings of Zaw et al., (2021), who found that 2-4 year high-frequency periodicities from spectral peaks, and predominant regions of high spatial correlations indicated the summer rainfall in Myanmar is associated with large-scale atmospheric circulations, mainly linked with the ENSO events due to SST deviations in the tropical Pacific Ocean. The variations in the western region's dry condition were observed in the late 1960s and 1980s ( Figure 5d) and in the late 2000s in southern region ( Figure 5e). In the western part, the monsoon rainfall experienced an increasing trend initially with an amount of 250 mm, but a decreasing pattern in the last few decades. In the southern part, the opposite pattern to the western region was observed, recording the highest rainfall (530 mm) over last decades, being the highest amongst all regions. The maximum (minimum) rainfall was recorded in 1960, 1980, 1998, 2006, and 2011 (1987, 1991, 1999, 2002, 2004, 2005, and 2014) in the region.
global-warming-induced changes in monsoon and circulation patterns in the regions; similar rainfall patterns were previously reported [83].
The results for various regions in Myanmar revealed that the west, south, and north regions received the most rainfall, while the east and central regions received relatively less ( Figure 5). Figure 5 shows that an overall increasing trend in rainfall occurred over all regions, which corresponds to previously reported findings [18,84,85]; however, the annual rainfall variability on the regional scale differed. In the northern part, the rainfall variability appeared to be increasing, with an amount of 500 mm. The peak rainfall years (dry years) were recorded during 1992, 1997, and 2010 (1982, 1991, 2005, and 2009). In the eastern part, monsoon rainfall slightly increased, as depicted by a rainfall amount of 250 mm with dry and wet years appearing to be the possible reason for balancing the overall trend in the region. In this region, the highest (lowest) monsoon rainfall was observed in 1988 and 2013 (1987 and 2009). In the central part, the monsoon rainfall exhibited an increasing trend at the rate of 300 mm per year. The analysis showed that 1960, 1997, 2010, and 2014 (1987, 1991, 1999 and 2009) were the peak wet (dry) years in the northern and central parts. Our results concur with the findings of Zaw et al. (2021), who found that 2-4 year high-frequency periodicities from spectral peaks, and predominant regions of high spatial correlations indicated the summer rainfall in Myanmar is associated with largescale atmospheric circulations, mainly linked with the ENSO events due to SST deviations in the tropical Pacific Ocean. The variations in the western region's dry condition were observed in the late 1960s and 1980s ( Figure 5d) and in the late 2000s in southern region ( Figure 5e). In the western part, the monsoon rainfall experienced an increasing trend initially with an amount of 250 mm, but a decreasing pattern in the last few decades. In the southern part, the opposite pattern to the western region was observed, recording the highest rainfall (530 mm) over last decades, being the highest amongst all regions. The maximum (minimum) rainfall was recorded in 1960, 1980, 1998, 2006, and 2011 (1987, 1991, 1999, 2002, 2004, 2005, and 2014) in the region.  The time series of interdecadal variability in rainfall over the whole region is shown in Figure 6. The mean rainfall was observed to increase in the late 1950s, reaching its peak around 1965, and subsequently decreased during the mid-1980s before increasing again in the early 1990s. The second climax was observed in the early 2000s; however, the rainfall decreased in the recent decade. We classified wet and dry years based on the standardized rainfall anomaly as ≥1 and ≤-1 for wet and dry years, respectively. The same approach was used effectively by other authors (Chen and Huang, (2017)  The results from the sequential Mann-Kendall test statistic were calculated for interdecadal rainfall variability over the region to detect the statistically significant turning points in the decadal trend in rainfall (Figure 7). An abrupt interdecadal change was noted during the 1970s; however, there was a significant reduction (at the 95% confidence level) in rainfall between the 1980s and 2000s (Figure 7), which agrees with spatial map shown in Figure 4.
in Figure 6. The mean rainfall was observed to increase in the late 1950s, reaching its peak around 1965, and subsequently decreased during the mid-1980s before increasing again in the early 1990s. The second climax was observed in the early 2000s; however, the rainfall decreased in the recent decade. We classified wet and dry years based on the standardized rainfall anomaly as ≥1 and ≤-1 for wet and dry years, respectively. The same approach was used effectively by other authors (Chen and Huang, (2017) in China, Lone et al. (2019) in India, and  in Myanmar). The wet periods were 1962-1967, 1998, 2003, and 2004; typical dry years were observed during the periods of 1978-1986 and 1988-1990. The results from the sequential Mann-Kendall test statistic were calculated for interdecadal rainfall variability over the region to detect the statistically significant turning points in the decadal trend in rainfall (Figure 7). An abrupt interdecadal change was noted during the 1970s; however, there was a significant reduction (at the 95% confidence level) in rainfall between the 1980s and 2000s (Figure 7), which agrees with spatial map shown in Figure 4.    in Figure 6. The mean rainfall was observed to increase in the late 1950s, reaching its peak around 1965, and subsequently decreased during the mid-1980s before increasing again in the early 1990s. The second climax was observed in the early 2000s; however, the rainfall decreased in the recent decade. We classified wet and dry years based on the standardized rainfall anomaly as ≥1 and ≤-1 for wet and dry years, respectively. The same approach was used effectively by other authors (Chen and Huang, (2017) in China, Lone et al. (2019) in India, and  in Myanmar). The wet periods were 1962-1967, 1998, 2003, and 2004; typical dry years were observed during the periods of 1978-1986 and 1988-1990. The results from the sequential Mann-Kendall test statistic were calculated for interdecadal rainfall variability over the region to detect the statistically significant turning points in the decadal trend in rainfall (Figure 7). An abrupt interdecadal change was noted during the 1970s; however, there was a significant reduction (at the 95% confidence level) in rainfall between the 1980s and 2000s (Figure 7), which agrees with spatial map shown in Figure 4.     6 and 7 show that regional variability in monsoon rainfall over Myanmar was obvious during the 1950-2019 period. The magnitude of the rainfall variability could possibly be associated with the geographical location of the region and monsoon circulation patterns. The temporal trend in the Asian and Indian monsoon system showed a strong variation on different time scales, including monthly, inter-annual, intra-seasonal, and annual time scale [18,84,85]. Such large-scale variation due to different factors could induce changes in monsoon rainfall over the diverse study regions [21,86,87].

Interdecadal Rainfall Variability
The EOFs were examined with the first three eigenvectors of the summer monsoon rainfall in Myanmar for the decadal 10 year running mean ( Figure 8). In Myanmar, EOF approaches have been used to detect the interannual variability in summer monsoon rainfall [31,83]. These methods are suitable for finding the total variance in monsoon rainfall in the region. The results reveal that the first EOF (EOF1), explaining 47% of the total variance, showed a dipole spatial pattern with positive loading on the western coast (Rakhine region) and negative loading in northwestern region (Sagaing Region) over the target region (Figure 8a, top and bottom), and explained the corresponding PC-1, which clearly exhibited interdecadal variability. EOF2, explaining 18% of the total variance, exhibited a dipole pattern, with the opposite pattern to EOF1 (Figure 8b, top and bottom), revealing that the PC time series is related to EOF2 (PC2). The third EOF (EOF3), explaining 11% of the total variance, showed strong positive loading in the southern region, and the time series of the consistent PC is connected to EOF3 (PC3), as shown in Figure 8c. To validate the significance of the scales of interdecadal variation, we used the power spectrum from the results EOF: the PC for 53 years was analyzed ( Figure 9). PC1-PC3 exhibited a significant spectral peak detected at nearly 20-30 year time scales.
Water 2021, 13, x FOR PEER REVIEW 13 of 28 patterns. The temporal trend in the Asian and Indian monsoon system showed a strong variation on different time scales, including monthly, inter-annual, intra-seasonal, and annual time scale [18,84,85]. Such large-scale variation due to different factors could induce changes in monsoon rainfall over the diverse study regions [21,86,87].

Interdecadal Rainfall Variability
The EOFs were examined with the first three eigenvectors of the summer monsoon rainfall in Myanmar for the decadal 10 year running mean (Figure 8). In Myanmar, EOF approaches have been used to detect the interannual variability in summer monsoon rainfall [31,83]. These methods are suitable for finding the total variance in monsoon rainfall in the region. The results reveal that the first EOF (EOF1), explaining 47% of the total variance, showed a dipole spatial pattern with positive loading on the western coast (Rakhine region) and negative loading in northwestern region (Sagaing Region) over the target region (Figure 8a, top and bottom), and explained the corresponding PC-1, which clearly exhibited interdecadal variability. EOF2, explaining 18% of the total variance, exhibited a dipole pattern, with the opposite pattern to EOF1 (Figure 8b, top and bottom), revealing that the PC time series is related to EOF2 (PC2). The third EOF (EOF3), explaining 11% of the total variance, showed strong positive loading in the southern region, and the time series of the consistent PC is connected to EOF3 (PC3), as shown in Figure 8c. To validate the significance of the scales of interdecadal variation, we used the power spectrum from the results EOF: the PC for 53 years was analyzed ( Figure 9). PC1-PC3 exhibited a significant spectral peak detected at nearly 20-30 year time scales.   From Figures 8 and 9, we inferred that the interdecadal rainfall pattern over Myanmar has significantly changed over the study period, especially during EOF2 and EOF3. PC1 and PC2 experienced an increasing spectrum over the target region during the period 1950-2019. The humid regions of the country appeared to experience decreasing rainfall, as can be seen in Figure 8c, while the arid regions in the northwestern area experienced an increase in monsoon rainfall. The humid central region appeared to be the core monsoon region of Myanmar; thus, this could imply a shift toward the northwestern parts of the country. Further studies with more data and modeling approaches are needed to validate the current hypotheses.

Large-Scale Atmospheric Circulation Pattern
Rainfall variability is influenced by atmospheric circulation, geographical location, and low latitude in a country like Myanmar. In this section, the general circulation anomalies over low latitudes are examined as a major issue. Circulation patterns on an interdecadal scale are analyzed over the region. To examine the interdecadal variability in the country's rainfall related to the typical circulation patterns on the interdecadal time scale, we employed a low level of horizontal and vertical wind anomalies of 850 hPa for the wet and dry periods ( Figure 10). The results show a significant negative westerly or southwesterly winds anomaly from the BoB and Andaman Sea over almost the whole region except the eastern and southern parts during wet years, as shown in Figure 10a. During dry years (Figure 10a), significant negative northeasterly and easterly wind anomalies prevailed in the region. The southwesterly and BoB appeared, which are assumed to be the prime sources of oceanic water transport to continental land masses for rainfall during monsoon seasons ( Figure S2). A similar pattern was observed for wind anomalies, further explaining dry conditions and elevated air temperature due to less rainfall. The continental land mass of South East Asia (SEA) experiences a strong anticyclonic pattern, which (in the northern hemisphere) is associated with decreased wind and rainfall. The anticyclonic pattern appeared to be pushing the cyclonic activity toward the tropical region and northern parts of the BoB by pushing the south easterlies to the peninsular tip of India ( Figure S2). From Figures 8 and 9, we inferred that the interdecadal rainfall pattern over Myanmar has significantly changed over the study period, especially during EOF2 and EOF3. PC1 and PC2 experienced an increasing spectrum over the target region during the period 1950-2019. The humid regions of the country appeared to experience decreasing rainfall, as can be seen in Figure 8c, while the arid regions in the northwestern area experienced an increase in monsoon rainfall. The humid central region appeared to be the core monsoon region of Myanmar; thus, this could imply a shift toward the northwestern parts of the country. Further studies with more data and modeling approaches are needed to validate the current hypotheses.

Large-Scale Atmospheric Circulation Pattern
Rainfall variability is influenced by atmospheric circulation, geographical location, and low latitude in a country like Myanmar. In this section, the general circulation anomalies over low latitudes are examined as a major issue. Circulation patterns on an interdecadal scale are analyzed over the region. To examine the interdecadal variability in the country's rainfall related to the typical circulation patterns on the interdecadal time scale, we employed a low level of horizontal and vertical wind anomalies of 850 hPa for the wet and dry periods ( Figure 10). The results show a significant negative westerly or southwesterly winds anomaly from the BoB and Andaman Sea over almost the whole region except the eastern and southern parts during wet years, as shown in Figure 10a. During dry years (Figure 10a), significant negative northeasterly and easterly wind anomalies prevailed in the region. The southwesterly and BoB appeared, which are assumed to be the prime sources of oceanic water transport to continental land masses for rainfall during monsoon seasons ( Figure S2). A similar pattern was observed for wind anomalies, further explaining dry conditions and elevated air temperature due to less rainfall. The continental land mass of South East Asia (SEA) experiences a strong anticyclonic pattern, which (in the northern hemisphere) is associated with decreased wind and rainfall. The anticyclonic pattern appeared to be pushing the cyclonic activity toward the tropical region and northern parts of the BoB by pushing the south easterlies to the peninsular tip of India ( Figure S2). Water 2021, 13, x FOR PEER REVIEW 15 of 28 During wet years (Figure 10b), the subtropical and Somali jet streams appeared to be pushing the oceanic moisture-rich winds toward the continental land mass. The southern region received northeasterly and easterly winds from the South China Sea, which is the core region experiencing monsoon rainfall. However, the anomalous cloud liquid water contents were obvious over the Western SEA, which includes Myanmar, with enhanced winds speed and reduced air temperature, implying increased rainfall. The anticyclonic pattern appeared to move to the BoB region strengthening the intensity of sub-tropical jets, pushing them further into continental land masses.
To assess the rainfall variability and its relationship with regional-scale atmospheric circulation, oceanic temperature, and monsoon indices, a composite analysis was performed. To do so, the regional-scale monthly and annual rainfall standardized index was developed to the inter-annual scale.
To further clarify the atmospheric circulation and its possible association with monsoon rainfall, we assessed the pressure vertical velocity (ω) for the wet and dry periods along latitude 15°, 20°, and 25° N, as shown in Figure 11. The results reveal that upward (rising) motions over the region were significant positive anomalies at 15° and 20° N (Figure 11a,c,e). Sinking motion occurred over the region during dry years, characterized by the significant positive anomaly at 15° and 20° N (Figure 11b,d,f). The difference in wet minus dry period moisture transport at 850 hPa revealed that anomalous moisture convergence (positive anomalies) emerge in the central and northern regions where moisture comes from BoB and the Gulf of Thailand ( Figure S1). We also identified a relationship between the distribution of air space at lower and higher levels and the flow of moisture transport. An unfamiliar airfield at low levels aligned with the location of a large junction of the complete flow of moisture in a stormy area. Figure S2 shows the distribution of field streams for the diversity of summer varieties for wet drying seasons at 850 hPa. The cyclonic broadcast in Southern Myanmar during the period 1950-2019 occurred in the southwest, with its center around 10° N (Taninthary region). Therefore, more water vapor entered the region, especially since the convergence center of the corresponding humidity was associated with a clear climate zone that cuts the region, increasing rainfall in that region. In these cases, water vapor from the BoB and the Andaman Sea was transported into the land.  (1962-1967, 1998, and 2003-2004) and (b) dry years (1978-1986 and 1988-1990). The red (blue) color indicates the significant positive (negative) anomalies over the target region.
During wet years (Figure 10b), the subtropical and Somali jet streams appeared to be pushing the oceanic moisture-rich winds toward the continental land mass. The southern region received northeasterly and easterly winds from the South China Sea, which is the core region experiencing monsoon rainfall. However, the anomalous cloud liquid water contents were obvious over the Western SEA, which includes Myanmar, with enhanced winds speed and reduced air temperature, implying increased rainfall. The anticyclonic pattern appeared to move to the BoB region strengthening the intensity of sub-tropical jets, pushing them further into continental land masses.
To assess the rainfall variability and its relationship with regional-scale atmospheric circulation, oceanic temperature, and monsoon indices, a composite analysis was performed. To do so, the regional-scale monthly and annual rainfall standardized index was developed to the inter-annual scale.
To further clarify the atmospheric circulation and its possible association with monsoon rainfall, we assessed the pressure vertical velocity (ω) for the wet and dry periods along latitude 15 • , 20 • , and 25 • N, as shown in Figure 11. The results reveal that upward (rising) motions over the region were significant positive anomalies at 15 • and 20 • N (Figure 11a,c,e). Sinking motion occurred over the region during dry years, characterized by the significant positive anomaly at 15 • and 20 • N (Figure 11b,d,f). The difference in wet minus dry period moisture transport at 850 hPa revealed that anomalous moisture convergence (positive anomalies) emerge in the central and northern regions where moisture comes from BoB and the Gulf of Thailand ( Figure S1). We also identified a relationship between the distribution of air space at lower and higher levels and the flow of moisture transport. An unfamiliar airfield at low levels aligned with the location of a large junction of the complete flow of moisture in a stormy area. Figure S2 shows the distribution of field streams for the diversity of summer varieties for wet drying seasons at 850 hPa. The cyclonic broadcast in Southern Myanmar during the period 1950-2019 occurred in the southwest, with its center around 10 • N (Taninthary region). Therefore, more water vapor entered the region, especially since the convergence center of the corresponding humidity was associated with a clear climate zone that cuts the region, increasing rainfall in that region. In these cases, water vapor from the BoB and the Andaman Sea was transported into the land.

Role of Oceans
From the composite analysis, we elucidated that large-scale atmospheric circulation had less (not-significant) influence on the summer monsoon rainfall variability over Myanmar. The observed relationship was rather weaker and could not be associated with rainfall variability over the study region. However, the composite analysis revealed that the anomalous wind pattern led the rainfall variability due to the regional rainfall dependence on moisture/water vapor. The trend in rainfall was obvious but the driver of the trend appeared to be complex; a simple mean trend may not provide enough explanation since Figure 11. Pressure vertical velocity (ω) (×10 −3 Pa s −1 ) over Myanmar during wet years (1962-1967, 1998, and 2003-2004) and dry years (1978-1986 and 1988 -1990

Role of Oceans
From the composite analysis, we elucidated that large-scale atmospheric circulation had less (not-significant) influence on the summer monsoon rainfall variability over Myanmar. The observed relationship was rather weaker and could not be associated with rainfall variability over the study region. However, the composite analysis revealed that the anomalous wind pattern led the rainfall variability due to the regional rainfall dependence on moisture/water vapor. The trend in rainfall was obvious but the driver of the trend appeared to be complex; a simple mean trend may not provide enough explanation since sub-seasonal variability may balance each other, which may smooth or even remove the trend. Figure 12 reveals a correlation map between the mean May-October SST (over the Indian Ocean and the Pacific Ocean) and mean May-October rainfall over Myanmar. A significant correlation (hatched regions) is displayed at the 95% confidence level in Figure 12. The results show a strong negative correlation between rainfall over Myanmar and Indian Ocean SST, which means that more (less) rainfall was recorded over Myanmar and with cooling (warming) over the Indian Ocean. The results further indicate that the air temperature in the northwest part of the SEA was obviously higher than in the southeast part of the region, which may affect the indicator of regional extreme temperature. Interestingly, we detected increasing trends in warm temperature events in the northern part of the region, in agreement with previous studies in the target regions [81,88]. Kreft and Eckstein (2013b) examined the monsoon rainfall over Myanmar and found a negative correlation with the positive phase of the Indian Ocean dipole (IOD).
Water 2021, 13, x FOR PEER REVIEW 17 of 28 sub-seasonal variability may balance each other, which may smooth or even remove the trend. Figure 12 reveals a correlation map between the mean May-October SST (over the Indian Ocean and the Pacific Ocean) and mean May-October rainfall over Myanmar. A significant correlation (hatched regions) is displayed at the 95% confidence level in Figure 12. The results show a strong negative correlation between rainfall over Myanmar and Indian Ocean SST, which means that more (less) rainfall was recorded over Myanmar and with cooling (warming) over the Indian Ocean. The results further indicate that the air temperature in the northwest part of the SEA was obviously higher than in the southeast part of the region, which may affect the indicator of regional extreme temperature. Interestingly, we detected increasing trends in warm temperature events in the northern part of the region, in agreement with previous studies in the target regions [81,88]. Kreft   This instability in the increase/decrease in SST and maximum/minimum rainfall may be associated with the increasing (decreasing) trend in temperature over the target region, which might affect the extreme temperature occurrences in the study area. The increase in frequency and magnitude of daily extreme warm temperature and decrease in extreme cold temperature may occur on the global scope, resulting in the increase in length, frequency, and/or intensity of warm periods or heat waves for most land areas [4,33]. Strong winds can transport water vapor from the ocean to nearby coastal areas, which can affect regional temperatures and ultimately affect the cooling trends, and vice versa [4,5].
The interdecadal variability in the 10 year running mean in Myanmar rainfall anomaly was computed using standardized deviation, PDO, and AMO indices, with the results provided in Figure 13. The effect of positive PDO was low rainfall (dry/drought), and that of negative PDO was more rainfall (wet/flood) in the region (Figure 13). The PDO index was statistically significant with a strong negative correlation with rainfall (-0.79), and the AMO index correlation coefficient explained the positive correlation (0.53), as shown in Table 4. The correlation coefficients between the Myanmar rainfall and various climate indices for the interannual and interdecadal components are presented in Table 4. The correlation coefficients between PDO and AMO indices and Myanmar rainfall for both This instability in the increase/decrease in SST and maximum/minimum rainfall may be associated with the increasing (decreasing) trend in temperature over the target region, which might affect the extreme temperature occurrences in the study area. The increase in frequency and magnitude of daily extreme warm temperature and decrease in extreme cold temperature may occur on the global scope, resulting in the increase in length, frequency, and/or intensity of warm periods or heat waves for most land areas [4,33]. Strong winds can transport water vapor from the ocean to nearby coastal areas, which can affect regional temperatures and ultimately affect the cooling trends, and vice versa [4,5].
The interdecadal variability in the 10 year running mean in Myanmar rainfall anomaly was computed using standardized deviation, PDO, and AMO indices, with the results provided in Figure 14. The effect of positive PDO was low rainfall (dry/drought), and that of negative PDO was more rainfall (wet/flood) in the region (Figure 14). The PDO index was statistically significant with a strong negative correlation with rainfall (-0.79), and the AMO index correlation coefficient explained the positive correlation (0.53), as shown in Table 4. The correlation coefficients between the Myanmar rainfall and various climate indices for the interannual and interdecadal components are presented in Table 4. The correlation coefficients between PDO and AMO indices and Myanmar rainfall for both the interannual and interdecadal components were significant at the 95% confidence level. The Myanmar rainfall showed a highly negative correlation with the PDO index and a positive correlation with the AMO index, particularly for the interdecadal component. The cumulative role of these large-scale indices, especially during monsoon season, is perceived as affecting the SEA monsoon circulation pattern and thus resulting in dry/wet conditions during El Niño/La Niña years. The findings of the current study showed no obvious or significant association between ENSO, PDO, and monsoon onset timing with rainfall variability over Myanmar. The cold phase of the AMO is generally associated with negative Myanmar rainfall anomalies, which is lower than normal rainfall. The warming or positive phase PDO may decrease rainfall (drought) conditions, and cooling or negative-phase PDO may increase rainfall (flood) conditions over Myanmar. Moreover, there were nine extreme years (1952, 1959, 1961, 1965, 1970, 1973, 1974, 1999, and 2001), as shown in Figure 13. The results of extreme rainfall measures during the negative phases of the PDO show that no extreme events occurred during the positive phase of the PDO (Figure 13).
Various studies of linearly coupled dominating patterns between the global SST and rainfall variations in Myanmar related to the interdecadal variability have been conducted using the SVD method [4,31,37,81]. This method has successfully been applied in the region to identify the interannual variability in summer monsoon rainfall [89]. Figure 15 indicates the heterogeneous correlation of the first three SVD modes amongst global SST and Myanmar rainfall. The SVD of the SST modes closely represents the long-term global warming pattern (Figure 15a,e), with a correlation coefficient of 0.94 between the time series of its principal components and the low-pass-filtered global SST anomaly. The southern and eastern regions of Myanmar and regions of the mid-west extending up to Indian Ocean are strongly correlated with SST in the tropical and northern Pacific (Figure 15a,b). We found subtle differences in the global spatial structure of the tropical Pacific SST correlation pattern in this period (Figure 15c). Although it is unclear if the differences in the global SST spatial patterns for the three modes regional rainfall are statistically significant, they do provide some insights. The findings reported using heterogeneous correlation patterns [11] as the first SVD of SST mode 1 were mainly associated with the cold phase of PDO in the Pacific Ocean, north of 20 • N, and the warm phase of AMO in the North Atlantic Ocean (Figure 15a). The rainfall mode 1 explains much of the rainfall over the regions except for the central core and eastern regions (Figure 15b). The squared covariance fraction of this mode is about 53%, and the correlation coefficient between the two fields is 0.94. SVD2 of SST mode 2 demonstrated a warm SST in the eastern equatorial Pacific SST, like the El Niño SST pattern in the Pacific Ocean, and a cold SST in the North Atlantic Ocean (Figure 15c). The result show that relatively low rainfall over the region contributed to rainfall mode 2 ( Figure 15d). Thus, the El Niño SST pattern in the Pacific Ocean and cold phase of the AMO were associated with dry/drought conditions over the region. The time series of the coefficients for the three SVD modes of global ocean SST and Myanmar rainfall are presented in Figure 16. The squared covariance fraction of SVD mode 2 was detected at a rate of 33%, and the correlation coefficient between the two fields' time coefficients is around 0.94. SST mode 3 explains cold SSTs in the central and eastern equatorial Pacific SSTs as a La Niña SST pattern (Figure 15e). Rainfall mode 3 exhibits enhanced rainfall over the region (Figure 15f). SVD3 of the squared covariance fraction of this mode was detected with a 6.7% rate, and the correlation coefficient between the two fields is 0.93 (Table 5). Therefore, we inferred that the role of the oceanic circulation is important for decadal Myanmar rainfall.    [4,31,37,81]. This method has successfully been applied in the region to identify the interannual variability in summer monsoon rainfall [89]. Figure 15 indicates the heterogeneous correlation of the first three SVD modes amongst global SST and Myanmar rainfall. The SVD of the SST modes closely represents the long-term global warming pattern (Figure 15a,e), with a correlation coefficient of 0.94 between the time series of its principal components and the low-pass-filtered global SST anomaly. The southern and eastern regions of Myanmar and regions of the mid-west extending up to Indian Ocean are strongly correlated with SST in the tropical and northern Pacific (Figure  15a,b). We found subtle differences in the global spatial structure of the tropical Pacific SST correlation pattern in this period (Figure 15c). Although it is unclear if the differences in the global SST spatial patterns for the three modes regional rainfall are statistically significant, they do provide some insights. The findings reported using heterogeneous correlation patterns [11] as the first SVD of SST mode 1 were mainly associated with the cold phase of PDO in the Pacific Ocean, north of 20° N, and the warm phase of AMO in the the interannual and interdecadal components were significant at the 95% confidence level.
The Myanmar rainfall showed a highly negative correlation with the PDO index and a positive correlation with the AMO index, particularly for the interdecadal component. The cumulative role of these large-scale indices, especially during monsoon season, is perceived as affecting the SEA monsoon circulation pattern and thus resulting in dry/wet conditions during El Niño/La Niña years. The findings of the current study showed no obvious or significant association between ENSO, PDO, and monsoon onset timing with rainfall variability over Myanmar. The cold phase of the AMO is generally associated with negative Myanmar rainfall anomalies, which is lower than normal rainfall. The warming or positive phase PDO may decrease rainfall (drought) conditions, and cooling or negative-phase PDO may increase rainfall (flood) conditions over Myanmar. Moreover, there were nine extreme years (1952, 1959, 1961, 1965, 1970, 1973, 1974, 1999, and 2001), as shown in Figure 14. The results of extreme rainfall measures during the negative phases of the PDO show that no extreme events occurred during the positive phase of the PDO ( Figure 14).

Discussion
In this study, we evaluated the interdecadal variability in the summer monsoon rainfall (May-October) in Myanmar (1950-2019) using EOF, SVD, EEMD, and correlation analyses based on 10 year running means. This study is the first of its kind to use a large number of statistical approaches to investigate the interdecadal variability in monsoon rainfall, which is beneficial for the scientific community and climate change adaptation in Myanmar. The results reveal that the dominant summer monsoon regions (south, central, and east) observed dry conditions during 1998 and 2002.  described the lowest rainfall over Myanmar since 1997, and the most intensive heat wave (60%) occurred in 1998 during an ENSO year. Webster and Yang (1992) [90] showed that 1998 was the driest year during the period of 1950-2019, with extreme temperature recorded in May 1998 (43.6 °C). Our results agree with the findings of previous studies [27], reporting similar results over the target regions. We used EOF analysis with the 10 year running means of summer monsoon rainfall, and the results show that the first leading EOF pattern explained 47% of the variance. The maximum loading occurs along the western coast and

Discussion
In this study, we evaluated the interdecadal variability in the summer monsoon rainfall (May-October) in Myanmar (1950-2019) using EOF, SVD, EEMD, and correlation analyses based on 10 year running means. This study is the first of its kind to use a large number of statistical approaches to investigate the interdecadal variability in monsoon rainfall, which is beneficial for the scientific community and climate change adaptation in Myanmar. The results reveal that the dominant summer monsoon regions (south, central, and east) observed dry conditions during 1998 and 2002. Sein et al., (2015) described the lowest rainfall over Myanmar since 1997, and the most intensive heat wave (60%) occurred in 1998 during an ENSO year. Webster and Yang (1992) [90] showed that 1998 was the driest year during the period of 1950-2019, with extreme temperature recorded in May 1998 (43.6 • C). Our results agree with the findings of previous studies [27], reporting similar results over the target regions. We used EOF analysis with the 10 year running means of summer monsoon rainfall, and the results show that the first leading EOF pattern explained 47% of the variance. The maximum loading occurs along the western coast and the minimum in the northwest region, while the corresponding PC1 exhibits a dipole pattern. Our results support the findings of Zhang et al., (2014), who reported similar results over South China using a drought index. The circulation patterns on an interdecadal scale 2343 analyzed using the composites (dry and wet years) in this study. Interdecadal variability in circulation mode consists of upward (rising) motion over the region, significantly positive anomaly over central Myanmar. The southwest moisture flux comes from the north Indian Ocean across the entire west coast of Myanmar.
Reckien and Petkova (2019) [91] studied the first component of EOF of the vertically integrated moisture transport vector in summer, and strong northward moisture transport from the tropical Indian and West Pacific Oceans dominated the whole East Asian region. In this study, PDO index values were found to be significantly related to Myanmar summer monsoon rainfall (  . They found that the overall negative relationship between ENSO and rainfall patterns had a strong positive association across Myanmar. The role of ENSO in the resulting rainfall patterns during El Niño years was stronger during the cold phase of the PDO compared with during La Niña years. The difference between rainfall during the different phases of ENSO was superimposed with the warm and cold periods of the PDO, showing overall higher average annual rainfall during cold PDO episodes [92]. The results reveal that the Indian monsoon is more vulnerable to drought when El Niño events occur during warm phases of the Pacific interdecadal variability. Conversely, wet monsoons are more likely to prevail when La Niña events coincide during cold phases of the Pacific interdecadal variability. Based on the EEMD analysis, highly extreme summer monsoon rainfall (i.e., 1952, 1959, 1961, 1965, 1970, 1973, 1974, 1999, and 2001) occurs during the negative phases of the PDO; no extreme events occurred during the positive phase of the PDO. Our results support the findings of Chen and Huang (2017), who described that between 1951 and 2012, there were five extreme events that had a strong correlation with the negative phases of the PDO. The AMO index was significantly associated with Myanmar summer monsoon rainfall recorded, at 0.54. Here, to determine the interdecadal summer monsoon rainfall modes and SST, a further SVD analysis was performed. The results of SVD mode 1 explain that the positive phase of the AMO receives higher summer monsoon rainfall, except in central and east regions in Myanmar. The positive phase of the AMO significantly correlated with higher summer monsoon rainfall and negatively correlated with lower summer monsoon rainfall over Myanmar, which is in agreement with .
Our findings do not negate the potential impact of large-scale oceanic indices but rather postulate the hypothesis that these controls need to be reviewed with clear links to summer monsoon rainfall variability over Myanmar. Another possible reason could be the quality of station data, since they are recorded for synoptic-scale use and climatological processes usually occur on longer time scale; thus, the data may lack such large-scale signals. These are some of the points that need to be thoroughly investigated in future studies, including the frequency and intensity of extreme rainfall events and their role as drivers of regional rainfall variability.

Conclusions
In this study, we attempted to explore the spatiotemporal trend and variability in monsoon rainfall during the period of 1950-2019 over Myanmar. Using statistical and composite analyses, we found significant variation in rainfall on seasonal, annual, and interdecadal scales. The main conclusions based on the findings are as follows: (1) The interdecadal variation over the target region increased in the late 1950s, reaching its peak around 1965, and subsequently decreased in the mid-1980s before increasing again in the early 1990s. An abrupt rainfall shift was recorded during the 1970s and a significant reduction (at the 95% confidence level) in rainfall was evident between the 1980s and the 2000s. (2) In terms of rainfall trend, the regional variability in monsoon rainfall over Myanmar was obvious during the 1950-2019 period. The magnitude of the rainfall variability could possibly be associated with the geographical location of the region and monsoon circulation patterns. The temporal trend over the Asian and Indian monsoon systems showed strong variation on different time scales, including monthly, inter-annual, intra-seasonal, and annual time scales. Such large-scale variation caused by different factors could induce changes in monsoon rainfall over the diverse study regions. (3) The widely used EOF method showed that the interdecadal rainfall pattern over Myanmar significantly changed over the study period, especially during EOF2 and EOF3. PC1 and PC2 experienced an increasing trend over the target region during the period 1950-2019. (4) The results further demonstrate the significant negative westerly and southwesterly wind anomalies from the BoB and Andaman Sea over almost the whole region except for the eastern and southern regions during wet years. The moisture transport at 850 hPa revealed that the anomalous moisture convergence (positive anomalies) emerged in the central and northern region, transporting moisture from the Western BoB and Eastern Gulf of Thailand. This cyclonic circulation may cause more water vapor from the BoB and Andaman Sea to be transported into the country. In addition, warm or positive PDO phase was associated with low rainfall in Myanmar, and cooling or negative PDO phase can increase rainfall over Myanmar. The AMO index had a mean average value of 0.53 for the same period. In addition, the effects of excessive rainfall occurred in the negative stages of the PDO, and no serious adverse events occurred in the positive phase of the PDO. The first SVD of SST mode 1 was mainly linked to the PDO cold phase in the Pacific Ocean above the North Pacific Ocean (20 • N) and the warm AMO phase in the North Atlantic Ocean. (5) The SST patterns showed that the warm SST of the eastern equatorial Pacific SST is a pattern of the El Niño SST in the Pacific Ocean and the cold SST in the North Atlantic Ocean. The fractional covariance fraction of SVD2 was recorded as 33%, whereas the coefficient of correlation between the coefficients for the two-phase period was 0.94. Therefore, we concluded that the findings of the present study provide valuable results regarding the interdecadal variability in Myanmar summer monsoon rainfall, which is shown as a negative correlation with the PDO and a positive association with the AMO. In addition, Myanmar's decadal summer monsoon rainfall and relationships with the PDO and AMO indexes explain the flood and dry conditions (drought) in the region.

Informed Consent Statement: Not Applicable.
Data Availability Statement: Not Applicable.