Multi-Decadal Trend and Decadal Variability of the Regional Sea Level over the Indian Ocean since the 1960 s : Roles of Climate Modes and External Forcing

Previous studies suggest that anthropogenic warming has affected the multi-decadal trend patterns of sea level over the Indian Ocean (IO). This effect, however, has not been quantified. Using observational datasets combined with large ensemble experiments from two climate models, this paper assesses the effects of natural internal variability versus external forcing on the observed, multi-decadal trend pattern and the decadal sea level anomaly (SLA) of the IO since the 1960s. Because the global mean sea level rise (SLR), which results largely from external forcing, has been removed before the examination, the paper focuses on the regionally uneven distribution of trend and SLA. The impacts of climate modes are quantified using a Bayesian Dynamic Linear Model. For the regional trend pattern of 1958–2005, the effects of internal variability dominate external forcing. Over the Seychelles area where sea-level variations obtain the maximum, internal variability (external forcing) contributes 81% (19 ± 2.4%) of the observed trend. For decadal SLA, internal variability is the predominant cause, with a standard deviation (STD) ratio of externally forced/observed SLA being 18 ± 17% over Seychelles and 17 ± 11% near the Indonesian Throughflow (ITF) area. Climate modes account for most observed SLA during boreal winter, with the total effects of decadal ENSO, Indian Ocean Dipole (IOD), and monsoon accounting for 78–86% of the observed STD near the Seychelles region, ITF area, and coasts of Sumatra and the Bay of Bengal. During summer, climate modes explain 95% of observed STD near the ITF but only 58–67% in other regions. Decadal ENSO dominates the SLA in the south tropical IO for both seasons and near the coasts of Sumatra and the Bay during winter. Decadal IOD and monsoon, however, control the coastal SLA during summer. Remote and local winds over the IO are the main drivers for decadal SLA, while the Pacific influence via the ITF is strong mainly in the southeast basin.


Introduction
Many low-lying coastal areas and island nations of the Indian Ocean (IO) are highly vulnerable to sea level rise (SLR) as they harbor one-third of the world's populations that are heavily populated by developing countries, which have poor disaster management resources [1].The rapid population growth will further increase their vulnerability in the future [2].Consequently, there is a strong societal demand for improved projections of future SLR and decadal predictions of sea-level variations over the IO region [1,3,4].Yet, our understanding of IO regional, decadal sea-level variations and their major drivers is very limited, whereas this understanding is instrumental for reliable decadal predictions and uncertainty quantifications of future SLR projections.This paper aims to assess the effects of natural internal variability, particularly climate modes, and external forcing (natural plus anthropogenic) on the observed IO decadal sea level variability and multi-decadal trend pattern, with global mean SLR (due largely to external forcing) being removed.Here, "decadal variability" is collectively referred to as variations from one to a few decades.

Background
Regional sea level variability can be induced by various factors: changes in atmospheric and oceanic circulations, large-scale deformation of ocean basins, the variation of Earth's gravity field, and local land movement [4][5][6].Among these factors, changes in atmospheric and oceanic circulations are shown to be the major cause for basin-wide spatial patterns of contemporary decadal sea-level variations [5].Over the tropical IO, multi-decadal trends of sea level since the 1960s exhibit a distinct spatial pattern, with the falling sea level occurring in the southwest basin where the thermocline ridge and the Seychelles Islands are located [7][8][9][10] but rising sea level elsewhere [11][12][13].Experiments using ocean general circulation models (OGCMs) show that the changing surface wind is the major force for this basin-wide sea level pattern [11], and variations of the Indonesian Throughflow (ITF) may have also contributed to the thermocline cooling [14] and, thus, the sea level fall in the tropical southwest IO [15].While the multi-decadal trend of the surface wind may be largely associated with natural internal variability, experiments using an atmospheric GCM suggest that the changing wind could also be partly forced by IO warming [11], which was attributed primarily to anthropogenic greenhouse gas (GHG) forcing [16,17].What are the relative roles of natural internal climate variability versus external forcing (for example, solar, volcano, GHG) in causing the multi-decadal trend pattern of sea level over the IO?This is an important science question remaining to be answered.
Overlying the multi-decadal trend, satellite observations show large decadal variations of basin-wide sea level patterns in recent decades, with reversing trends between 1993-2000 and 2000-2006 [18].Over the North IO, the sea level experienced basin-wide falls from 1993-2003 but sharp rises from 2004-2013 [19,20].Both observational analyses and OGCM experiments suggest that the surface wind over the IO is the primary cause for the observed decadal sea level variability [19][20][21][22], with the ITF making significant contributions primarily in the southeast tropical IO [21,23].Thermosteric sea level is the primary contributor [22,24,25], with the halosteric sea level having apparent influence in some regions [22], such as the tropical southeast IO [26].While surface wind associated with natural internal variability is suggested to be the major cause for decadal sea level variability [5,27], the effect of external forcing on regional, decadal sea-level variations over the IO remains unknown.
To what extent is the IO decadal sea level variability associated with internal climate modes?To answer this question, Li and Han performed a suite of OGCM experiments to assess the forcing by surface wind stress associated with climate modes, stochastic wind, and buoyancy flux from 1948-2012 [28].To extract the wind stress signals associated with climate modes, they first performed multiple linear regression onto the monthly Multivariate ENSO Index (MEI) and IO dipole (IOD) mode index (DMI) [29,30], using NCEP reanalysis wind and HadISST sea surface temperature (SST) data.Then, they removed the ENSO-and IOD-associated signals from the surface wind and SST and applied singular value decomposition (SVD) to obtain the wind-SST covariance for the residuals.The sum of the wind stress anomalies associated with ENSO, IOD, and the leading 7 SVD modes of the residuals represent the total effect of all climate modes.Their results suggested that the decadal sea level variability north of 20 • S is primarily induced by wind stress associated with climate modes, with the largest variability occurring in the thermocline ridge region.
To assess the effect of an individual climate mode-the Pacific Decadal Oscillation (PDO)-on Indo-Pacific decadal sea level variability, Frankcombe et al. obtained the multiple linear regression patterns by regressing the reanalysis and OGCM-simulated sea level onto the PDO index [31].Their results showed that the PDO is moderately correlated with decadal sea level variability in some regions of the tropical IO, and the regression coefficients are much weaker than those of ENSO and IOD on interannual timescales, and the association of sea level with the Southern Annular Mode (SAM) is weak.
Note that both these studies [28,31] used the multiple linear regression approach, which assumes stationary relationships (that is, constant regression coefficients) between climate modes and their associated fields (for example, surface winds and sea level).In the real climate system, however, their relationships are often non-stationary and evolve with time.For instance, it has been shown that the influence of El Niño on the subtropical northwest Pacific climate has significantly increased after the 1970s [32]; decadal variability of Indian summer monsoon convection is more important since the 1990s compared to prior decades in determining the decadal variability of the Indian Walker Cell [33].While strong El Niño events generally induce a weaker Indian summer monsoon and thus, severe drought, some strong El Niño events correspond to normal monsoons [34].The conventional multiple linear regression cannot capture these changing relationships, which is a limitation in its application (see Section 2.4 for details).Moreover, the relative importance of climate modes in determining decadal sea-level variations in different regions of the IO is unknown, and the seasonality of their impact remains unclear [28].Yet, understanding the impacts of climate modes has important implications for climate predictions, because internal climate modes with physical bases, such as IOD and ENSO, possess predictability and can be predicted with positive skills two or three seasons in advance [35,36], even though no decadal-scale predictability of ENSO behavior was found in a long-control simulation of a global climate model [37].

Present Research
The goal of this research is to quantify the effects of internal climate modes and external forcing on the observed IO regional, uneven distribution of the sea level trend and decadal variability (with global mean SLR removed).To achieve the goal, we first analyze the observed patterns of multi-decadal trend and decadal sea level variability (departure from the global mean), and then extract the imprints of internal climate modes in the observed sea level and surface winds using a relatively new approach, the Bayesian Dynamic Linear Model (DLM) [38].The Bayesian DLM can capture the non-stationary impacts of climate modes, which cannot be achieved by the conventional multiple linear regression, referred to as static linear model (SLM) in previous studies [33,38].Comparisons between Bayesian DLM and SLM in assessing the climate modes' impacts on decadal variability of Indian and Pacific Walker Cells reveal that the DLM considerably outperforms SLM in simulating the Walker Cells' variability, and has captured the non-stationary nature of climate modes' impacts over different decades [33].To examine the effects of external forcing, we analyze the results from large ensemble (LE) experiments of two global climate models and compare them with observational datasets (see Section 2 for more details) because the influence of external forcing is difficult to quantify using observational analysis alone.Our major period of interest is from the 1960s to the present when more reliable ocean-atmosphere datasets are available.Results from this research will contribute to a better understanding of decadal predictions and uncertainties of the externally forced, regional sea level signals over the IO and, therefore, benefit a large fraction of the region's populations.
Here, we focus on decadal-scale variability for two reasons.First, decadal prediction emerges as a pressing priority in climate research due to strong society demands for planning purposes [39][40][41][42][43][44][45].Second, empirical approaches, including the Bayesian DLM, build simultaneous relationships between climate mode indices and sea level variability (Section 2.4), indicating that variations of climate modes are in quasi-equilibrium with variations of sea level.The physical bases for their relationships are that climate modes affect sea level through affecting surface forcing fields (for example, surface wind, precipitation) [28][29][30]36].While surface wind and rainfall over the IO can respond quickly (a few days) to SST associated with climate modes (for example, ENSO) [29, 30,36], variations of sea level are not always in equilibrium with surface wind forcing due to oceanic wave propagations.On decadal timescales, surface wind and sea level over the tropical IO basin are largely in quasi-equilibrium, but this is untrue at inter-annual and shorter timescales (see Section 2.5 for further discussion).

Datasets: Sea Level and Surface Wind
To document the observed multi-decadal trend and decadal variability of IO regional sea level, we analyze monthly sea level data from the ECMWF Ocean Reanalysis System 4 (ORAS4) available for 1958-2015 [46], and from the Simple Ocean Data Assimilation (SODA version 2.2.4) reanalysis available for 1871-2010 [47].Note that the OGCMs for producing these reanalysis products made the Boussinesq approximation and, thus, implicitly assume the total ocean volume is conserved, even though the global mean temperature is rising and continental ice is melting.Consequently, the global mean SLR is excluded and the regional sea-level variations (both trend pattern and decadal variability) analyzed here are departures from the global mean (that is, the dynamical sea level).The ORAS4 data, however, assimilated the satellite-observed sea level [46], yielding a global mean rising trend since 1993 [48].Variances of decadal sea level anomaly (SLA) from satellite altimetry are reasonably reproduced by both reanalysis data over the IO, even though there are quantitative differences in some regions [49].Since the thermosteric sea level is the dominant component of basin-wide decadal SLA (Section 1), the upper 700 m thermosteric sea level data from the world ocean atlas 2013 (WOA13) for 1955-2015 and from Ishii and Kimoto for 1945-2014 [50,51] are also analyzed.Finally, even though the record is short, the multiple-satellite merged monthly SLA data from the French Archiving, Validation, and Interpolation of Satellite Oceanographic Data (AVISO) project [52,53] for 1993-2015 are also examined.
To help understand the effects of climate-modes-associated surface winds on IO decadal SLA, we also examine the monthly surface winds from the Japan Meteorological Agency 55-yr Reanalysis (JRA55) from 1958-2015 [54].Han et al. compared various reanalysis wind products, and showed that the decadal variability of the Indian Walker Cell, represented by zonal surface wind of the equatorial IO since the 1960s, and its seasonality are robust to cross-dataset differences, and they agree with the zonal sea level gradient detected by ORAS4 and upper 700 m thermosteric sea level [33].The decadal variability and seasonality of the Indian Walker Cell detected by the JRA55 wind also agree with the cross-calibrated multiplatform (CCMP) satellite winds version 1 [55] during their overlapping period [33].The CCMP monthly winds since 1987, satellite-observed monthly Outgoing Longwave Radiation (OLR)-a proxy of tropical deep convection-available since 1974 [56], and Global Precipitation Climatology Project (GPCP) monthly precipitation available since 1979 [57] are also analyzed.
To isolate decadal signals, a Lanczos lowpass filter [58] with half power at the 8 year period is used.The response curve of the 8 year lowpass filter retains ~90% of the amplitude at the 10 year period and almost full amplitudes for longer periods, which essentially retains the decadal variability signals with periods of 10 year and longer.Similar lowpass filters have been successfully applied to investigate decadal variability by many existing studies [31,33,[59][60][61].

Large Ensemble (LE) Experiment Results from Two Climate Models
To assess the contribution of external forcing on the regional sea level trend and decadal SLA over the tropical IO, we utilize the recently available LE simulations with all (both natural and anthropogenic) forcing from two global climate models: the National Center for Atmospheric Research (NCAR) community earth system model version 1 (CESM1) [62] and the Max-Planck-Institute for meteorology (MPI, Hamburg, Germany) earth system model, an updated version of the same model used in CMIP5 [63].Monthly data from the 40-member ensemble of CESM1 simulations are available for 1920-2005, and that from the 100-member ensemble of MPI are available for 1850-2005.Their monthly dynamical sea level data are analyzed and the model drifts are removed before the analyses.The 40-and 100-member ensembles are used to assess the effect of external forcing.
Note that the cross-model differences of regional sea level variability are large between the two LE model results, not only for multi-decadal trends but also for decadal SLAs.Consequently, we focus only on analyzing and comparing their SLA magnitudes (for example, variances) rather than their phases, when we estimate the effects of external forcing.Note also that we use the LE experiments only to assess the contribution of external forcing relative to the observed SLAs.Examining and comparing the natural internal variability versus externally forced SLAs in each climate model is beyond the scope of this paper.

Datasets for Climate Modes and Other Climate Indices
To understand the IO decadal sea level variability associated with internal climate modes, we examine the relationships between sea level and ENSO and IOD, the two dominant inter-annual climate modes over the Indo-Pacific basin with evident decadal modulation [15,27,64,65].Note that decadal variability of ENSO, represented by the 8 year low-passed Nino3.4 index and MEI, is highly correlated with the PDO and Inter-decadal Pacific Oscillation (IPO) indices, with Nino3.4-IPO and IPO-PDO correlations both being ~0.88 from 1900-2008 and the IPO-MEI correlation being 0.89 from 1900-2001 [61,66].Following Han et al. (2017) [33], we use the 8 year low-passed Nino3.4 index to represent the decadal variability of the ENSO or IPO.Hereafter, we refer it to as decadal ENSO.For the IOD, its decadal variability is represented by the 8 year low-passed DMI.Here, we calculate the DMI for each month, whereas the IOD is active only during summer-to-fall and peaks in Sep-Nov [33].Consequently, the "DMI" here essentially represents the east-west SST gradient during January-May.For convenience, we still refer it to as the IOD.Indices of ENSO and IOD are calculated using monthly HadISST data [67] available since 1870, and then the 8 year low-passed Lanczos filter is applied.
In addition to ENSO and IOD, the decadal variability of monsoon winds over the IO can drive the large variability of IO sea level.Since the strengths of monsoon winds are tied to that of monsoon rainfall over India during boreal summer, the Indian Institute of Tropical Meteorology homogeneous monthly All India Rainfall dataset available from 1871-2014 [68] is used to serve as the monsoon index.Note that the leading Empirical Orthogonal Function (EOF) mode of IO SST (IOSST) is a basin-wide warming/cooling pattern at both inter-annual and decadal timescales [15,61,69].Because the basin-scale IOSST is primarily induced by ENSO [61,69], its influence is largely imprinted in the ENSO index.Since the 1980s, however, the relationship between decadal ENSO and IOSST has reversed [15,61] due to volcanic eruptions and anthropogenic warming [70,71].As a result, decadal IOSST has had an active impact on the Pacific trades and sea level during recent decades [61,72].In this paper, the effect of decadal IOSST, the part that is independent of decadal ENSO, will also be examined using Bayesian DLM (see below).

The Bayesian DLM and Partial DLM
Since the Bayesian DLM and partial DLM have been discussed in detail in Han et al. (2017) [33], here, we only discuss the aspects essential to this paper.The DLM consists of two equations: an "observation equation" shown by Equation (1), and a "state equation" represented by Equation (2): The observation Equation (1) builds the relationship between predictand Y and predictors (X 1 , . . .X M ), where M is the number of predictors, and the state Equation (2) controls the dynamical evolution of coefficients b i (i = 0, 1, 2, . . .M). Equation (2) states that the predictive distribution of b i at each time step t (that is, posterior) is updated based on its previous step t − 1 distribution (that is, prior) and the probability of observations Y conditional on b i at time t (that is, the likelihood) using Bayes theorem [38].Coefficient b i (t) is obtained by applying Kalman filtering and smoothing, and it measures the "dynamical" (time-evolving) impacts of X i on Y, which overcomes the limitation of the constant coefficients of the conventional multiple linear regression.In the above, term b 0 (t) represents a time-varying intercept whose variability is unexplained by the predictors X i , while the b i terms represent the non-stationary influence of X i on Y; ε(t) and w i (t) are independent white noises or errors, distributed normally with a mean of 0 and variances of V(t) and W i (t).The choices and calculations for W i (t) and V(t) follow those in Reference [38].For more details of the Bayesian DLM, please see References [33,38].
Here, predictors X i represent the climate mode indices, and predictand Y represents the SLA.Since the decadal variability of ENSO and IOD are moderately correlated, yet the DLM requires the predictors to be independent, we apply the partial Bayesian DLM to remove the ENSO influence from both the DMI and SLA before we examine the IOD effect.Specifically, using ENSO as the predictor, Equation (1) becomes Then, we apply the Bayesian DLM to the residual, SLA = SLA(t) − b 1 (t)ENSO(t), using DMI and monsoon index (with ENSO signal removed) as predictors: Note that in Equation ( 4), the monsoon effect on IOD is also removed from the DMI.In addition, the decadal IOSST, which is induced by the external forcing since the 1980s [71], may also influence IO convection, winds and, thus, sea level.To consider this effect, we perform partial Bayesian DLM to the SLA residual, In Equation (5), residual , and the IOSST index with the decadal ENSO effect removed may, to some degree, represent external forcing.The time series of the decadal variations of these climate indices are shown in Figure 1.
Note that in Equations ( 3)-( 5), we removed the ENSO signals from both the DMI and the monsoon index, as well as monsoon signal from the DMI, to ensure the independence of the predictors.This procedure implies that ENSO affects the IOD and monsoon, and the monsoon affects the IOD, but not the other way around.The caveat is that the IOD and monsoon are integral components of the tropical biennial oscillation [73,74] and that the IOD can also have an active impact on ENSO [75], even though ENSO does have a significant influence on the IOD, particularly during northern winter [36].Consequently, part of the IOD impact on SLA has been imprinted in the effect of ENSO and the monsoon in Equations ( 3) and ( 4), which may result in the overestimation of decadal ENSO and monsoon influences and the underestimation of the decadal IOD effect.
In Bayesian DLM, the SLA is determined not only by the magnitudes of climate indices, but also by their corresponding time-varying coefficients b i (t), as shown by Equations (1)- (5).Physically, this indicates that not all strong climatic events can produce strong SLA.For instance, the strong 1997-1998 eastern Pacific El Niño is associated with a normal Indian summer monsoon, contrasting the drought year that an El Niño often produces.This is because it is the central Pacific El Niño that is more effective in producing subsidence and, therefore, drought over India, compared to the eastern Pacific El Niño [34].For sea level, the eastern and central Pacific El Niño events-with both being represented by the Nino3.4index in Bayesian DLM-may have different impacts on the IO regional SLA even though their Nino3.4indices are the same, because their differences in convection patterns may drive different wind anomalies over the IO.This changing relationship can be captured by b i (t) Climate 2018, 6, 0 7 of 36 using Bayes theorem in Equation ( 2), whereas it cannot be captured by the constant b i of SLM.In SLM, the SLA is only determined by the magnitude of the Nino3.4index.The Bayesian DLM approach has been recently applied to climate studies to understand the non-stationary effects of ENSO and IOD on the inter-annual variability of the Indian Summer Monsoon [76] and Indonesian rainfall [77], and it has also been used to explore the climate modes' impacts on decadal variability of the Indian and Pacific Walker Cells [33].

Choice of Analysis Region
Our analyses focus on the tropical IO north of 20° S, because at 20° S, the oceanic adjustment time, which is the time it takes Rossby waves to cross the basin, is ~5 year for the first and ~11 year for the second baroclinic modes, the two gravest modes of oceanic response to surface forcing.Equatorward of this latitude, the adjustment times are shorter and decrease with the decrease of latitude [78][79][80].This choice ensures that on decadal and longer timescales, the Bayesian DLM can be applied with zero-lag between the climate mode indices and sea level, as shown by Equations ( 3) and (4).This is because Rossby waves associated with the first two gravest baroclinic modes, together with the barotropic mode which is much faster, have propagated across the tropical IO basin within

Choice of Analysis Region
Our analyses focus on the tropical IO north of 20 • S, because at 20 • S, the oceanic adjustment time, which is the time it takes Rossby waves to cross the basin, is ~5 year for the first and ~11 year for the second baroclinic modes, the two gravest modes of oceanic response to surface forcing.Equatorward of this latitude, the adjustment times are shorter and decrease with the decrease of latitude [78][79][80].This choice ensures that on decadal and longer timescales, the Bayesian DLM can be applied with zero-lag between the climate mode indices and sea level, as shown by Equations ( 3) and (4).This is because Rossby waves associated with the first two gravest baroclinic modes, together with the barotropic mode which is much faster, have propagated across the tropical IO basin within ~10 years' time.Consequently, variations of surface forcing fields (for example, wind) associated with climate modes are approximately in equilibrium with variations of sea level at timescales of 10 year and longer.

Results
In this section, we first examine the regional, multi-decadal trend pattern of the sea level over the IO from 1958-2005 when both observational datasets and climate model results are available, and assess the contributions of external forcing versus internal variability to the observed regional trend (Section 3.1).Here, we focus on understanding the "regional distribution of trend" after the global mean SLR, which results largely from anthropogenic warming, being removed.Then, we investigate the decadal variations of the observed sea level (with both global mean and regional trend removed), assess the relative magnitudes of the natural internal variability versus external forcing, and extract the part that is associated with internal climate modes (Section 3.2).Finally, we discuss the importance of winds in driving the regional, decadal sea level patterns (Section 3.3).

Multi-Decadal Trends
The ORAS4 and SODA products show overall consistent spatial patterns of sea level trend from 1958-2005, with the maximum sea-level fall located in the southwestern tropical basin Seychelles Islands region, the weaker sea level fall in the Arabian Sea, and the SLR in the eastern equatorial basin and the Bay of Bengal (Figure 2a,b) where freshwater fluxes due to strong precipitation and the Bay of Bengal river's discharges are large [81].These patterns resemble those of earlier studies from the standalone ocean model experiments, reanalysis data, and reconstructed sea level for a similar period of time [11][12][13]82].Given that global mean SLR has been removed from each dataset, our results suggest that the rates of SLR near the coasts of Sumatra, Java, and the Bay of Bengal from 1958-2005 are considerably larger than the global mean.Note that the vertical land motions detected at some locations near the coasts and the Islands of the IO for the past 20 year [83] may have also affected the sea level trends there.
Qualitative and quantitative differences, however, exist between ORAS4 and SODA.One apparent disparity is in the southeastern basin near the ITF region, where ORAS4 exhibits SLR whereas SODA shows sea level fall.Indeed, the upper 700 m thermosteric sea level data from both WOA13 and Ishii show falling trends in this area, and also in the Bay of Bengal (Figure 2c,d) where ORAS4 and, to a lesser extent, SODA detect rising trends (Figure 2a,b).Recall that ORAS4 assimilated satellite altimeter data since 1993 but SODA did not (see Section 2.1), and the thermosteric sea level excludes the salinity effect.The disagreements between ORAS4 and SODA (also thermosteric sea level) in the ITF and Bay of Bengal regions suggest that the halosteric sea level associated with an intensified hydrological cycle may have contributed to the multi-decadal SLR in these areas since the 1950s [84-88].Quantitatively, there are also significant differences.For instance, the falling rate averaged over the Seychelles Islands (boxed region of Figure 2a) is −1.74 ± 0.12 mm year −1 in ORAS4 but −0.73 ± 0.12 mm year −1 in SODA, which is ~42% that of ORAS4.Various reasons may contribute to the qualitative and quantitative discrepancies: for example, differences in their OGCMs and surface forcing fields, lack of in situ data that constrain the models in the Seychelles region, and different datasets (for example, altimeter data) assimilated into their models.
Indeed, the trend patterns of ORAS4 from 1993-2010 bear the most resemblance with those of the satellite data (Figure 3a,b), with a pattern correlation of r = 0.84 (>99% significance).By contrast, the SODA data show a much lower pattern correlation of r = 0.48 (Figure 3c).The upper 700 m-thermosteric sea level has a pattern correlation of 0.48 for AVISO/WOA13 and 0.73 for AVISO/Ishii (Figure 3d,e).Figure 3 suggests that ORAS4 data most accurately represent the dynamical sea level variability over the IO during the satellite era compared to other products; the good agreements between Ishii, AVISO, and ORAS4 data (Figure 3a,b,e) also suggests that the upper-ocean thermosteric sea level is the dominant component for decadal sea level variability, with halosteric sea level being less important even in the eastern basin where the salinity effect is suggested to be important for the multi-decadal trend of SLR (Figure 2).Nidheesh et al. [49] examined the decadal SLAs over the Indo-Pacific basin since the 1960s using 10 sea level products from in-situ based, reanalysis, and reconstructed datasets, and their results show significant differences over the IO basin.It is difficult to determine which dataset is the best prior to the satellite era because reanalyses and reconstructions suffer from different issues.For of sea level are either fixed to those of the satellite era [89] or assumed to be cyclostationary [82], and their temporal variability is constrained by tide gauge observations before the satellite era [60,82,89].There are, however, only a dozen tide gauge stations with record lengths greater than 50 year around the IO rim [90] and no data are available in the ocean interior.Since our research focus is to understand the regional patterns of sea-level variations, we use reanalysis products because they assimilated available observations throughout the basin.Nidheesh et al. [49] examined the decadal SLAs over the Indo-Pacific basin since the 1960s using 10 sea level products from in-situ based, reanalysis, and reconstructed datasets, and their results show significant differences over the IO basin.It is difficult to determine which dataset is the best prior to the satellite era because reanalyses and reconstructions suffer from different issues.For reanalyses, large errors may arise from uncertainties in the model forcing fields, model deficiency, and lack of long-term in situ records to constrain the models.For reconstructions, the spatial patterns of sea level are either fixed to those of the satellite era [89] or assumed to be cyclostationary [82], and their temporal variability is constrained by tide gauge observations before the satellite era [60,82,89].There are, however, only a dozen tide gauge stations with record lengths greater than 50 year around the IO rim [90] and no data are available in the ocean interior.Since our research focus is to understand the regional patterns of sea-level variations, we use reanalysis products because they assimilated available observations throughout the basin.
While the sea level trends from ORAS4 and SODA include the effects of both natural internal variability and external forcing, the 100-member ensemble average of MPI model (with global SLR removed) isolates the effect of external forcing on the regional distribution of sea-level variations.Its trend patterns resemble that of ORAS4 in most regions of the IO, such as the evident sea level fall in the southwest basin, a weaker fall in the Arabian Sea, and a rise in the southeast basin near the ITF region (compare Figure 2a,e).An exception is the eastern equatorial basin and the Bay of Bengal, where the MPI model produces a sea level fall instead of the SLR shown in ORAS4 and SODA (Figure 2; the right column of Table 1).This difference may reflect the true effect of external forcing, but it may also arise from the common climate model biases associated with convection and winds, which resulted in an underestimated SLR there in the eastern IO [91][92][93].Based on the 100-member ensemble mean of the MPI model, external forcing contributes −0.29 ± 0.04 mm year −1 to the falling trend of the thermocline ridge (Figure 2e; Table 1), which is ~17 ± 2.3% of that from ORAS4 and ~40 ± 5.5% from SODA.
Table 1.The linear trends (mm year −1 ) of regional sea level with global sea level rise (SLR) removed (that is, departure from the global mean) averaged in the Seychelles area and in the eastern equatorial basin (boxed regions of Figure 2) for the period of 1958-2005 from observations and climate model simulations.Observations are based on monthly ECMWF Ocean Reanalysis System 4 (ORAS4) and Simple Ocean Data Assimilation (SODA) reanalysis products, with all trend values exceeding 95% significance and errors being the slope uncertainties of linear fitting.The model simulations are from the Max-Planck-Institute of meteorology (MPI) 100-member ensemble, National Center for Atmospheric Research (NCAR) community earth system model version 1 (CESM1) 40-member ensemble, and their total of 140-member ensemble; their trend values are obtained from the 100-member ensemble averages of external forcing effect, with each of the 100 members being an ensemble mean of randomly chosen N members using the Monte Carlo method; uncertainties indicate 1.65 standard deviations of the external forcing spread, as discussed in detail in the text.

West: Tropical South IO (Seychelles)
(50 Compared to the MPI 100-member ensemble average, the 40-member ensemble mean of CESM1 shows some similarities in sea level patterns; but the falling trend extends farther to the south, yielding a ~0 trend averaged in the southwest basin boxed region over Seychelles.If we move the box southward by 10 • , the averaged value is −0.21 mm year −1 , which is close to the −0.29 mm year −1 from the MPI model (Figure 2e,f).Randomly chosen 40-member ensemble means of MPI show similar sea level trend patterns to the 100-member average, even though there are still some quantitative differences in some areas (compare Figure 2e,g,h).Over the north IO, the two models show overall similar patterns but with significant quantitative differences.While both model simulations yielded a sea level fall in the eastern equatorial basin, CESM1 produces a SLR but MPI generates a weak sea level fall in the Arabian Sea and the western Bay of Bengal.CESM1 shows falling sea level over the entire eastern IO, including the ITF region where ORAS4 and MPI model show rising sea level.This pattern resembles the bias of an overly strong east-west dipole pattern in the predecessor of the CESM1 model [91][92][93].
The above results indicate that the 40-member ensemble average can reasonably damp the natural internal variability in the multi-decadal (1958-2005) regional sea level trend, and the differences between the MPI 100-member and CESM1 40-member ensemble averages arise mainly from the differences between the two models.Since the observed IO warming in the past few decades can induce surface winds that likely drive a similar sea level pattern to Figure 2e in most regions of the tropical IO [11], and the warming has been attributed primarily to anthropogenic forcing [16,17], there is reason to believe that the MPI model may have produced more reasonable externally-forced sea level trend patterns than CESM1.Indeed, apparent differences between MPI and CESM sea level trends exist over the entire Indo-Pacific basin, with opposite signs occurring in some regions (for example, the western tropical Pacific; figure not shown).Their large disparity suggests that the two models have very different responses to the same external forcing, yielding different or even opposite sea level trends in some regions.
Using the 140-member ensemble mean of the two models (100 MPI and 40 CESM1) to represent external forcing (−0.23 ± 0.03 mm year −1 ), and the average of ORAS4 and SODA to represent observation (−1.23 ± 0.12 mm year −1 ), external forcing contributes 19 ± 2.4% of the observed falling trend over the Seychelles area (Model-mean versus Obs-mean in Table 1; see Supplementary Figure S1 for bar graph).Using the MPI 100-member average to represent external forcing, the ratio of the forced/observed trend value becomes 24 ± 3.3% over the Seychelles area.Since the climate models produce a falling trend in the eastern equatorial IO where observations detect a rising trend, the external forcing contributes less to the observed east-west contrasting pattern (Figure 2; Table 1) than to the Seychelles falling trend.Be aware that the mean of the two reanalysis products may not adequately represent the ensemble mean of the multiple reanalysis datasets [94] and the true observational values, giving rise to a larger uncertainty in the estimated ratio for external forcing.
Note that here, we use the Monte Carlo method to calculate the uncertainty of external forcing because the conventional approach of directly calculating the ensemble member spread (represented by 1.65 standard deviation (STD)) used in Slangen et al. and Meyssignac et al. [95,96] does not apply to our case.This is because after global SLR is removed, the ensemble average-which assesses the effect of external forcing-has a weak magnitude, and this weak signal is overwhelmed by the large-amplitude internal variability in each ensemble member.The trend spread (1.65STD) is ~2.5 times that of the trend signal of the ensemble mean for 1958-2005, and the spread is way larger than the signal for decadal SLA (Supplementary Figure S2).Physically, the huge spread demonstrates that natural internal variability dominates over external forcing and it cannot be used to represent the "uncertainty range" of external forcing here, even though it worked for the situation when global SLR was retained [95,96].To properly assess the uncertainty of the external forcing effect, we first randomly choose N members from one model ensemble, and then obtain the N-member ensemble average to represent "external forcing effect".Using this procedure, we obtain a set of 100 members for external forcing, with each being a N-member ensemble mean.Then we calculate the 100-member average and spread (1.65STD) to measure the effect of external forcing and uncertainty.Size N is chosen to ensure that the N-member ensemble mean can reasonably damp the internal variability.Based on our test and analysis (Supplementary Figure S3), choosing N as ~70% of the total members will reduce the uncertainty of the externally forced decadal SLA to below the signal values in both the western and eastern basins.Consequently, we use N = 98, 70, and 28 for the total model (MPI + CESM1), MPI, and CESM1, respectively.For the 1958-2005 trend, using N = 40 yields reasonable results (not shown).For consistency, however, we use the same N values for the multi-decadal trend (Table 1) and decadal SLA (discussed below) to assess the effects of external forcing and their uncertainties.

Natural Internal Variability Versus External Forcing
To understand in which areas the observed decadal (8 year lowpass filtered) SLAs obtain their largest amplitudes, we calculate their STD maps using the annual mean SLAs from the ORAS4 and SODA data (Figure 4a-c).Albeit with quantitative differences between the two reanalysis products, decadal SLAs reach their maximum magnitudes over the Seychelles region (Reg. 1 of Figure 4a-c) where the multi-decadal trend is also the strongest (see Figure 2), near the ITF area north of Australia (Reg. 2), and to a lesser degree, along the perimeter of the Bay of Bengal for both datasets.The maximum decadal SLA in the Seychelles region is consistent with the result of Li and Han (2015) [28].The STD of SLA averaged over Reg. 1 is 15.5 mm for ORAS4 and 9.8 mm for SODA, and the correlation between the two is r = 0.46 (below 90% significance; Figure 4g).For Reg. 2, the STD is 14.9 mm for ORAS4 and 17.1 mm for SODA, and their correlation coefficient is r = 0.72 (>95% significance).The less agreement over Reg. 1 further supports our earlier argument regarding the lack of observations in this area for constraining the models, yielding larger uncertainties of the detected decadal SLA in the Seychelles region.(Reg.2), and to a lesser degree, along the perimeter of the Bay of Bengal for both datasets.The maximum decadal SLA in the Seychelles region is consistent with the result of Li and Han (2015) [28].
The STD of SLA averaged over Reg. 1 is 15.5 mm for ORAS4 and 9.8 mm for SODA, and the correlation between the two is r = 0.46 (below 90% significance; Figure 4g).For Reg. 2, the STD is 14.9 mm for ORAS4 and 17.1 mm for SODA, and their correlation coefficient is r = 0.72 (>95% significance).The less agreement over Reg. 1 further supports our earlier argument regarding the lack of observations in this area for constraining the models, yielding larger uncertainties of the detected decadal SLA in the Seychelles region.External forcing also induces its largest decadal SLA over the Seychelles Islands, near the ITF and in the Bay of Bengal (Figure 4d-f), the regions where strong variability occurs in the reanalysis data (Figure 4a-c).The largest magnitude over Seychelles suggests that the shallow, mean thermocline there is prone to the influence of external forcing and internal climate variability.Comparing the MPI 100and 40-member ensemble averages, we see apparently weaker magnitudes in the 100-member mean (compare their STD values in Figure 4e-h).The STD of the 100-member averaged SLA is 2.58 mm in Reg. 1, which is ~17% of the STD of ORAS4 and ~26% of SODA.For Reg. 2, the 100-member mean SLA STD is 2.62 mm, which is ~18% of the STD of ORAS4 and ~15% of SODA.Compared to the MPI model, CESM1 produces larger SLAs (compare Figure 4d-h; also see Figure 5b,c,e,f).The correlation between the 40-member ensemble averages from the MPI and CESM1 is −0.39 for Reg. 1 and 0.2 for Reg. 2, and neither correlation coefficient is statistically significant at 90%.The large differences between the two models suggest that decadal SLAs induced by the external forcing of the climate system are highly sensitive to the models, and it is unknown which model produces the correct response in nature.For this reason, we cannot simply subtract the effect of external forcing from observations to obtain the effect of internal variability.External forcing also induces its largest decadal SLA over the Seychelles Islands, near the ITF and in the Bay of Bengal (Figure 4d-f), the regions where strong variability occurs in the reanalysis data (Figure 4a-c).The largest magnitude over Seychelles suggests that the shallow, mean thermocline there is prone to the influence of external forcing and internal climate variability.Comparing the MPI 100-and 40-member ensemble averages, we see apparently weaker magnitudes in the 100-member mean (compare their STD values in Figure 4e-h).The STD of the 100-member averaged SLA is 2.58 mm in Reg. 1, which is ~ 17% of the STD of ORAS4 and ~26% of SODA.For Reg. 2, the 100-member mean SLA STD is 2.62 mm, which is ~18% of the STD of ORAS4 and ~15% of SODA.Compared to the MPI model, CESM1 produces larger SLAs (compare Figure 4d-h; also see Figure 5b,c,e,f).The correlation between the 40-member ensemble averages from the MPI and CESM1 is −0.39 for Reg. 1 and 0.2 for Reg. 2, and neither correlation coefficient is statistically significant at 90%.The large differences between the two models suggest that decadal SLAs induced by the external forcing of the climate system are highly sensitive to the models, and it is unknown which model produces the correct response in nature.For this reason, we cannot simply subtract the effect of external forcing from observations to obtain the effect of internal variability.Using the 1962-2001 STD of the ORAS4 and SODA average to represent observation (STDo), and the 1962-2001 STD of the 100-member mean of the external forcing ensemble (see detailed descriptions above) from the total MPI + CESM1 to represent external forcing (STDem), we obtain their ratio, STDem/STDo, as 18 ± 17% in the Seychelles region and 17 ± 11% near the ITF area (Figure 5a,d).The uncertainty is the ratio between 1.65STD of the 100-member spread (STDes) and the observation, 1.65STDes/STDo, which is 17% in the Seychelles region and 11% near the ITF.These results clearly demonstrate that natural internal variability is the predominant cause for the observed decadal SLA in the two regions where both observed and externally forced signals reach their maxima (Figure 4).

Effects of Climate Modes
Given that ORAS4 data more accurately represent the decadal sea level trends than SODA during the satellite era of 1993-2010 (Figure 3) and that the dataset has a longer record, hereafter, we use ORAS4 data to assess the decadal SLAs associated with climate modes.
(a) Contribution in regions of large SLA.Evidently, climate modes plus IOSST, which is independent of decadal ENSO and partly represents external forcing (Section 2.4), can explain most of the observed variances of the decadal SLAs based on annual mean data, and climate modes (ENSO + IOD + monsoon) are the main contributor in all key regions where the observed SLA STDs are large (Figure 6a-c), with IOSST (Figure 6h) having identifiable contributions over Seychelles (Reg.1), the Sumatra coasts (Reg.3), and the Bay of Bengal coasts (Reg.4).Each individual mode is associated with distinct SLA patterns (Figure 6e-g), with ENSO contributing the most in the south IO over Seychelles and near the ITF (Reg.2), IOD over Seychelles, and monsoon over Seychelles, the Sumatra coasts, and Bay of Bengal coasts.It appears that monsoon variability dominates decadal SLA along the coasts of Sumatra and the Bay, with ENSO and IOD having weak contributions there, even though the east antinode of the IOD covers a large part of Reg. 3. As we shall see below, this contradiction arises partly from the strong seasonality of climate modes and, thus, the seasonality of their impacts.
To reveal the seasonality of the SLA associated with climate modes, we calculated the decadal SLA based on the boreal winter (November-April) and summer (May-October) mean values.These months are chosen for the following reasons: (1) ENSO reaches its peak amplitude during boreal winter, and Indian summer monsoon rainfall obtains its maximum during boreal summer, that is, May-October; (2) For the IOD, about one third of events are associated with ENSO and two thirds are independent of ENSO; the independent IOD events usually develop in June and quickly weaken after October [36]; and (3) November-April and May-October have been characterized as the wet and dry seasons of Indonesia that resides in the region of Indo-Pacific warm pool convection [77], and were used to document the seasonality of the relationships between decadal variability of the Indo-Pacific Walker Cells and climate modes [33].
During November-April, climate modes can explain most of the SLAs in all key regions, but IOSST also has non-negligible contributions (Figure 7a-c,h), similar to the annual mean (Figure 6).The residual, the portion unexplained by climate modes and IOSST, is weak in all other areas but Reg. 1, where both the observed SLA and data uncertainty are the largest, as discussed earlier.While the decadal ENSO has its largest influence in the south IO as the annual mean (Regs.1-2, Figure 7e), it also exerts a strong influence on SLA along the Sumatra and Bay of Bengal coasts (Regs.[3][4].The effect of the monsoon is comparable to ENSO in all regions except for the ITF area (Figure 7g).The IOD has a significant contribution only over Seychelles.The weak influence of the decadal IOD near the Sumatra coast is likely because we removed the ENSO signals from the IOD index in Equation (4).Consequently, the impact of the ENSO-associated IOD on SLA is imprinted in the effect of ENSO; this procedure may overestimate the impact of ENSO but underestimate the effect of IOD, a caveat that was pointed out in Section 2.4.Similar to winter, climate modes and IOSST (ENSO + IOD + monsoon + IOSST) explain most of the observed SLA variances in all regions during summer.Different from winter, climate modes contribute less to the observed variances near the Sumatra, Java, and Bay of Bengal coasts, and the errors are apparently larger in these regions compared to winter (compare Figures 7a-d and 8a-d).
In addition, during May-October, the strong impacts of ENSO are confined to the south IO (Regs. 1 and 2), with little influence on the SLAs around the Bay of Bengal rim and west coasts of Sumatra (Figure 8e).The IOD, on the other hand, has the strongest influence on SLAs near the Sumatra and Bay of Bengal coasts, compared to the impacts of ENSO and the monsoon (Figure 8e-g).The large IOD impact on the coastal areas of the eastern IO results from the fact that most of the IOD events that are independent of ENSO develop in early summer and quickly weaken after October [36].The Similar to winter, climate modes and IOSST (ENSO + IOD + monsoon + IOSST) explain most of the observed SLA variances in all regions during summer.Different from winter, climate modes contribute less to the observed variances near the Sumatra, Java, and Bay of Bengal coasts, and the errors are apparently larger in these regions compared to winter (compare Figures 7a-d and 8a-d).
In addition, during May-October, the strong impacts of ENSO are confined to the south IO (Regs. 1 and 2), with little influence on the SLAs around the Bay of Bengal rim and west coasts of Sumatra (Figure 8e).The IOD, on the other hand, has the strongest influence on SLAs near the Sumatra and Bay of Bengal coasts, compared to the impacts of ENSO and the monsoon (Figure 8e-g).The large IOD impact on the coastal areas of the eastern IO results from the fact that most of the IOD events that are independent of ENSO develop in early summer and quickly weaken after October [36].The IOD, together with monsoon and IOSST, is responsible for a large fraction of the observed SLAs near the east IO coasts (Figure 8f-h).

(b) Contribution to basin-wide EOF patterns.
To reveal the dominant spatial patterns of decadal SLA, we carry out an EOF analysis for both winter and summer from 1962-2011.During winter, the leading EOF (EOF1; explaining 38% variance) of ORAS4 SLA shows the falling sea level in the central and western basin with a maximum magnitude occurring in the Seychelles region, but rising sea level in the eastern basin with large amplitudes occurring near the ITF and coastal areas of the north IO (Figure 9a).This basin-wide pattern exhibits significant decadal variability, as shown by its principal component (PC1; the black curve of Figure 9d).Consistent with the STD analysis (Figure 7), the EOF1 of SLA from the sum of climate modes plus IOSST can reasonably represent the observed EOF1 pattern and its temporal evolution (Figure 9a,b), with the correlation between their PCs being 0.94 (Figure 9d); the climate modes account for most of the variability (Figure 9c,d), with decadal ENSO and monsoon (to a lesser degree) having large contributions, and the influence of decadal IOD being weak near the eastern IO coasts.During summer, similar EOF1 patterns are found from observations (Figures 9a and 10a), and the total effect of the climate modes plus IOSST can still explain a large fraction of the basin-wide pattern.Similar to the STD analysis (Figure 8), climate modes plus IOSST explain less SLA variance in the eastern IO coastal regions compared to winter (Figure 10a-c) and the correlations between the observed and DLM simulated PC1s are 0.72 for climate modes plus IOSST and 0.50 for climate modes, which are significantly lower than the 0.94 and 0.84 of winter.The effect of decadal ENSO is confined to the south IO (Figure 10e), and decadal IOD is important for causing the SLAs along the coasts of During summer, similar EOF1 patterns are found from observations (Figures 9a and 10a), and the total effect of the climate modes plus IOSST can still explain a large fraction of the basin-wide pattern.Similar to the STD analysis (Figure 8), climate modes plus IOSST explain less SLA variance in the eastern IO coastal regions compared to winter (Figure 10a-c) and the correlations between the observed and DLM simulated PC1s are 0.72 for climate modes plus IOSST and 0.50 for climate modes, which are significantly lower than the 0.94 and 0.84 of winter.The effect of decadal ENSO is confined to the south IO (Figure 10e), and decadal IOD is important for causing the SLAs along the coasts of the Sumatra and Bay of Bengal (Figure 10f).Different from the STD analysis, monsoon variability plays a larger role than IOD in determining the observed EOF1 pattern and its temporal evolution (Figure 10g,h).The dominance of decadal IOD in causing the regional SLAs near the coasts of Sumatra and Bay of Bengal (Figure 8) is not reflected in the EOF1s and PC1s, although the EOF2 of the IOD-associated SLA does show the strongest variability along the eastern IO coasts compared to ENSO and monsoon (not shown).Indeed, the EOF1 of climate modes essentially fails to capture the strong SLA along the coasts of Sumatra and Bay of Bengal (compare Figure 10a,c), even though the effect of IOSST improves this situation (Figure 10b).These differences arise from the fact that the basin-wide EOF1 patterns may not properly capture the regional/local features of the climate modes' impacts shown in the STD maps.(c) Time series analysis.To quantify the temporal evolution of the climate modes' impacts on regional SLA, we analyze the time series of decadal SLA averaged for each of the four key areas shown in Figures 6-10 where the SLA STDs are large.During winter, the results of Bayesian DLM show that the STD ratios between climate-modes-associated and observed SLAs (thin and thick solid curves of Figure 11) are 78%, 86%, 84%, and 80%, and their correlations are 0.80, 0.89, 0.87, and 0.85 (>95%) in Reg. 1, 2, 3, and 4, respectively.Adding the IOSST effect to the climate modes (green and (c) Time series analysis.To quantify the temporal evolution of the climate modes' impacts on regional SLA, we analyze the time series of decadal SLA averaged for each of the four key areas shown in Figures 6-10 where the SLA STDs are large.During winter, the results of Bayesian DLM show that the STD ratios between climate-modes-associated and observed SLAs (thin and thick solid curves of Figure 11) are 78%, 86%, 84%, and 80%, and their correlations are 0.80, 0.89, 0.87, and 0.85 (>95%) in Reg. 1, 2, 3, and 4, respectively.Adding the IOSST effect to the climate modes (green and black solid curves of Figure 9a-d), the STD ratios increased to 83%, 92%, 93%, and 91% and the correlations to 0.87, 0.96, 0.95, and 0.94, indicating that the external forcing has improved the STD ratios by 5-11% in Regs.1-4.The 5-6% values in Regs.1-2 are smaller than the external forcing of 17-18% estimated from the total of MPI + CESM1 relative to the average of the ORAS4 and SODA data (Figure 5; Section 3.2.1).This is expected, given that the IOSST may only partly reflect the effect of the external forcing since the 1980s (Section 2.4).Decadal ENSO is the largest contributor in all regions (the right column of Figure 11) and it is the controlling factor for the ITF area (Figure 11d).The effects of the decadal monsoon are also strong in all regions except for near the ITF, but their values are systematically weaker than that of decadal ENSO.The influence of IOD is relatively weak, with the largest effect (STD ratio of ~30% to observation) occurring over the Seychelles region.
During summer (Figure 12), the climate modes explain a lower percentage variance of the observed SLA in all regions (Reg. 1, 3, and 4) except for the ITF area (Reg.2).Over Seychelles, the Sumatra coasts, and the Bay of Bengal coasts, the STD ratios between the climate modes-associated and observed SLAs are 67%, 58%, and 63% with r = 0.76, 0.86, and 0.86, respectively.Adding IOSST increases the STD ratios to 73%, 75%, and 82% and correlation to 0.81, 0.93, and 0.95.Near the ITF, however, the climate modes contribute to 95% of the observed STD with r = 0.95, which are higher than their corresponding winter values.
While ENSO remains the largest contributor to decadal SLA over Seychelles, the Indian summer monsoon and, to a lesser extent, the IOD also have significant contributions in this region (Figure 12b).Indeed, the effect of decadal monsoon is comparable to that of decadal ENSO in the past decade or so, with an SLA STD being 9.5 mm due to the monsoon and 9.7 mm due to ENSO from 1997-2011.Consistent with Figure 8, decadal IOD is the largest contributor near the coasts of Sumatra and the Bay of Bengal overall from 1962-2011, and its STD is significantly larger than that of the decadal monsoon, with decadal ENSO being the least important (Figure 12f,h).Note, however, that decadal IOD has its strongest contribution from the late 1980s to early 1990s, and its influence is rather weak during other times (blue curves of Figure 12f,h).In contrast, the effect of the monsoon is the largest in the past decade or so, with SLA STD being 8.8 mm and 7.8 mm in Reg. 3 and Reg. 4, compared to the IOD contribution of 2.3 mm and 2.1 mm for 1997-2011.The increased influence of the monsoon on coastal SLAs in the past decade will be further discussed in Section 3.3.As for winter, decadal ENSO plays a deterministic role in causing SLA near the ITF region, with the monsoon also being important in the recent decade (Figure 12d).The effect of the IOD and IOSST are weak.During summer (Figure 12), the climate modes explain a lower percentage variance of the observed SLA in all regions (Reg. 1, 3, and 4) except for the ITF area (Reg.2).Over Seychelles, the Sumatra coasts, and the Bay of Bengal coasts, the STD ratios between the climate modes-associated and observed SLAs are 67%, 58%, and 63% with r = 0.76, 0.86, and 0.86, respectively.Adding IOSST  While ENSO remains the largest contributor to decadal SLA over Seychelles, the Indian summer monsoon and, to a lesser extent, the IOD also have significant contributions in this region (Figure 12b).Indeed, the effect of decadal monsoon is comparable to that of decadal ENSO in the past decade

Influences of Winds Associated with Climate Modes
To understand the effects of winds in driving the regional distributions of decadal sea level variability, we carry out composite analyses for decadal SLA, together with their associated surface wind stress and Ekman pumping velocity, defined as ρ f , where τ x and τ y are zonal and meridional wind stress, ρ = 1000 kg m −3 is the mean water density, and f is the Coriolis parameter.The composites are calculated for both the winter and summer seasons (Figures 13 and 14).For the November-April winter season, the years for the composite are chosen when the sea level fall (that is, the negative SLA) over Seychelles exceeds −0.5 STD in the ORAS4 data (Reg.1, Figure 11a); for the May-October summer season, the years of the composite are chosen when the sea level rise (that is, positive SLA) adjacent to the Bay of Bengal coasts exceeds +0.5 STD (Reg.4, Figure 12g).The criteria of the −1 STD are also tested and they yield similar results except for larger magnitudes for each field (not shown).
the IOD contribution of 2.3 mm and 2.1 mm for 1997-2011.The increased influence of the monsoon on coastal SLAs in the past decade will be further discussed in Section 3.3.As for winter, decadal ENSO plays a deterministic role in causing SLA near the ITF region, with the monsoon also being important in the recent decade (Figure 12d).The effect of the IOD and IOSST are weak.

Influences of Winds Associated with Climate Modes
To understand the effects of winds in driving the regional distributions of decadal sea level variability, we carry out composite analyses for decadal SLA, together with their associated surface wind stress and Ekman pumping velocity, defined as  = − , where  and  are zonal and meridional wind stress,  = 1000 kg m −3 is the mean water density, and  is the Coriolis parameter.The composites are calculated for both the winter and summer seasons (Figures 13 and  14).For the November-April winter season, the years for the composite are chosen when the sea level fall (that is, the negative SLA) over Seychelles exceeds −0.5 STD in the ORAS4 data (Reg.1, Figure 11a); for the May-October summer season, the years of the composite are chosen when the sea level rise (that is, positive SLA) adjacent to the Bay of Bengal coasts exceeds +0.5 STD (Reg.4, Figure 12g).The criteria of the −1 STD are also tested and they yield similar results except for larger magnitudes for each field (not shown).The vital role played by equatorial winds, together with the increased importance of the monsoon in the past decade, agrees well with the results of Han et al. (2017) [33], who showed that decadal variability of the Indian Walker Cell, measured by its surface branch of the zonal wind in the equatorial IO, was determined primarily by decadal variations of the Indian summer monsoon since the late 1990s [33].A thorough understanding of the non-stationary impacts of climate modes and their causes on the IO regional, decadal SLA is beyond the scope of this paper, but it will be investigated in our future research.During winter, the composite SLAs show a basin-scale pattern that resembles EOF1 (Figures 9a  and 13a).The positive Ekman pumping velocity w e in the south tropical IO (white contours of Figure 13a) drives surface mass divergence, lowering sea level and enhancing the upwelling of colder subsurface water into the surface.The colder water is denser and, thus, further lowers the sea level.Dynamically, the decadal SLAs are associated with the zonal integral of w e and on decadal timescales, the sea level variability of the tropical IO is in quasi-equilibrium with the wind forcing (Section 2.5).As a result, the amplitude of the sea level fall intensifies toward the west when w e keeps the same positive sign, attaining its maximum over Seychelles, but weakens near the western boundary where w e changes from positive to negative (black contours).
Near the equator, westerly wind anomalies cause surface Ekman convergence and raise the sea level.The high sea level signals propagate eastward as the equatorial Kelvin waves and then poleward along the eastern boundary as coastal Kelvin waves, increasing the sea level in the eastern equatorial basin and along the coasts of Sumatra, Java, and the Bay of Bengal (Figure 13a, color contours).Subsequently, the positive SLAs along the eastern boundary radiate westward as Rossby waves, weakened by the negative SLAs induced by the positive w e in the basin interior (for example, white contours over the Bay and in the south IO).Along the Sumatra and Java coasts, the negative w e induces Ekman convergence, which further increases the sea level in these regions.Over the Bay of Bengal, however, the local alongshore winds and w e are weak near the coasts; therefore, the coastal SLAs result mainly from the remote, equatorial wind forcing.The transient wave processes over the equatorial and north IO occur at timescales shorter than eighteen months [80,97] (Section 2.5).Consequently, the decadal SLAs are in equilibrium with decadal winds.Friction and the topographic beta effect associated with continental shelf, however, will trap part of the Rossby wave energy to the eastern boundary [98], contributing to the observed decadal SLAs near the coastal regions.The effects of topographic trapping and friction cannot be quantified by observational analyses, and high-resolution model experiments are needed to quantify these effects.
To extract the surface winds associated with climate modes, we apply the Bayesian DLM Equations ( 3)-( 5) to the 8 year low-passed JRA55 winds.Like sea level, the decadal variations of surface winds also result mainly from climate modes (Figure 13a-d), with ENSO being the primary contributor during winter (Figure 13e).The influences of winds associated with monsoon and IOD are also evident, except that the IOD contributes little to the coastal SLAs in the Bay of Bengal, and the effect of IOSST is weak (Figure 13f-h).
Similar to winter, the composite analyses of summer show the deterministic role of surface winds in driving the IO regional SLAs, with the positive w e accounting for the negative SLAs over Seychelles, and equatorial winds inducing increased sea level in the eastern basin (Figure 14a).The local w e also enhances the SLA along the coasts of Sumatra and Java.Different from winter, the local alongshore wind and negative w e near the eastern boundary of the Bay also enhance the positive SLAs around the perimeter of the Bay.Climate modes plus IOSST explain a large fraction of the SLA and wind (Figure 14a-c); the residual, however, is much larger than that of winter (Figures 13d and 14d).The IOD-associated wind is the major cause for the coastal SLAs of Sumatra and the Bay of Bengal, and the effect of Indian summer monsoon wind is also evident but apparently weaker in the composite (Figure 14f,g) compared to its effect shown in the time series (Figure 12).This is because during the composite years of 1993 and 1999-2001, the decadal SLAs induced by monsoon were weak (Figure 12f,h), even though the monsoon convection plays a dominant role during the past decade.The influence of ENSO is relatively weak (Figure 14e), but the effect of IOSST is stronger compared to winter (Figures 13h and 14h).
The vital role played by equatorial winds, together with the increased importance of the monsoon in the past decade, agrees well with the results of Han et al. (2017) [33], who showed that decadal variability of the Indian Walker Cell, measured by its surface branch of the zonal wind in the equatorial IO, was determined primarily by decadal variations of the Indian summer monsoon since the late 1990s [33].A thorough understanding of the non-stationary impacts of climate modes and their causes on the IO regional, decadal SLA is beyond the scope of this paper, but it will be investigated in our future research.Near the ITF area, w e anomalies are positive during both winter and summer, yet decadal SLAs are strong and positive (Figures 13a and 14a), suggesting that the w e anomalies are not the cause of the positive SLAs.While the strong variability of zonal equatorial winds and local alongshore winds near the northwest Australian coast may contribute to the decadal SLAs in Reg. 2 via coastal Kelvin and westward radiating Rossby waves, the remote forcing from the Pacific via the ITF can have a strong influence in this region (see Section 1).This is supported by the vital role of the decadal ENSO in determining the SLAs in this area for both seasons (Figures 11d and 12d).
To further confirm the composite SLA and surface wind patterns obtained from ORAS4 and JRA55 reanalyses, we compare them with satellite-observed AVISO SLA and CCMP1 winds for their overlapping period of 1993-2011 (Figure 15).First, the composite SLA and surface winds using ORAS4 and JRA55 data for 1998-2002 (Figure 15a), the five years when the November-April mean decadal SLAs exceeded the −0.5STD in Reg. 1 within the 1993-2011 period and resembled those of 1962-2011 (Figure 13a), albeit with larger SLA amplitudes for the shorter 1993-2011 period (note the scale changes in the color bars of Figures 13a and 15a).The 1998-2002 composites based on the 8 year low-passed SLAs are similar to those based on the unfiltered fields, and they are also similar to the satellite-observed SLAs and surface winds (Figure 15a-c).The surface wind anomalies are tied to the anomalous tropical deep convection and precipitation, converging to the increased convection (negative OLR) and precipitation over the Indo-Pacific warm pool and diverging from the reduced convection (positive OLR) and precipitation in the tropical southwest basin and northern IO (color and line contours of Figure 15d).Near the ITF area, we anomalies are positive during both winter and summer, yet decadal SLAs are strong and positive (Figures 13a and 14a), suggesting that the we anomalies are not the cause of the positive SLAs.While the strong variability of zonal equatorial winds and local alongshore winds near the northwest Australian coast may contribute to the decadal SLAs in Reg. 2 via coastal Kelvin and westward radiating Rossby waves, the remote forcing from the Pacific via the ITF can have a strong influence in this region (see Section 1).This is supported by the vital role of the decadal ENSO in determining the SLAs in this area for both seasons (Figures 11d and 12d).
To further confirm the composite SLA and surface wind patterns obtained from ORAS4 and JRA55 reanalyses, we compare them with satellite-observed AVISO SLA and CCMP1 winds for their overlapping period of 1993-2011 (Figure 15).First, the composite SLA and surface winds using ORAS4 and JRA55 data for 1998-2002 (Figure 15a), the five years when the November-April mean decadal SLAs exceeded the −0.5STD in Reg. 1 within the 1993-2011 period and resembled those of 1962-2011 (Figure 13a), albeit with larger SLA amplitudes for the shorter 1993-2011 period (note the scale changes in the color bars of Figures 13a and 15a).The 1998-2002 composites based on the 8 year low-passed SLAs are similar to those based on the unfiltered fields, and they are also similar to the satellite-observed SLAs and surface winds (Figure 15a-c).The surface wind anomalies are tied to the anomalous tropical deep convection and precipitation, converging to the increased convection (negative OLR) and precipitation over the Indo-Pacific warm pool and diverging from the reduced convection (positive OLR) and precipitation in the tropical southwest basin and northern IO (color and line contours of Figure 15d).Similarly, the filtered and unfiltered summertime composites of SLA and surface wind anomalies from reanalysis data, based on the four years of 1993 and 1999-2001 when ORAS4 SLA in Reg. 4 exceeded +0.5STD (Figure 12d), agree well with the satellite observations (compare Figures 14a and 15e-g), and the surface wind patterns are closely tied to anomalous convection and precipitation, converging to (diverging from) the enhanced (reduced) convection and precipitation areas (Figure 15h).The increased convection and precipitation in the eastern equatorial basin and reduced convection and precipitation in the west agree with the negative IOD events, which are shown to be the dominant cause for the summertime increased SLAs along the coasts of Sumatra and the Bay of Bengal (Figure 14).Over the Indian subcontinent, convection and precipitation are also enhanced, indicating a stronger summer monsoon situation.However, whether or not the negative decadal IOD events always accompany enhanced decadal convection of Indian summer monsoon requires further investigation.The data record is too short to make a definitive conclusion, and the four-year composite may still contain some inter-annual variability.
While the increased convection in the warm pool reduces the downward shortwave flux and, thus, decreases the sea level through thermosteric effect, the increased precipitation tends to increase the sea level via halosteric effect.The opposing influence of shortwave flux and precipitation may partly account for the compensating effect of the thermal and halosteric sea level on regional scales discussed in Stammer et al. [5].This compensation weakens the contribution of buoyancy flux to regional decadal SLA, which may further contribute to the dominance of winds in causing regional, decadal sea level variability over the IO.

Summary and Discussions
By analyzing observational datasets and the results from large ensemble experiments of two global climate models, we assessed the influences of natural internal variability versus external (natural plus anthropogenic) forcing on the observed multi-decadal trend patterns and decadal SLA in the tropical IO north of 20° S since the 1960s.Here, we focused on understanding the "regional trend pattern" and decadal SLA after removing the global mean SLR.Consequently, the significant effect of external forcing on the global mean SLR was removed before our analysis.Using Bayesian DLM, we quantified the influences of climate modes (that is, ENSO, IOD, and monsoon; Figure 1) on the IO regional, decadal SLA.The Bayesian DLM can capture the non-stationary impacts of climate modes on the regional sea level and, thus, it consistently outperforms the conventional multiple linear Similarly, the filtered and unfiltered summertime composites of SLA and surface wind anomalies from reanalysis data, based on the four years of 1993 and 1999-2001 when ORAS4 SLA in Reg. 4 exceeded +0.5STD (Figure 12d), agree well with the satellite observations (compare Figures 14a and  15e-g), and the surface wind patterns are closely tied to anomalous convection and precipitation, converging to (diverging from) the enhanced (reduced) convection and precipitation areas (Figure 15h).The increased convection and precipitation in the eastern equatorial basin and reduced convection and precipitation in the west agree with the negative IOD events, which are shown to be the dominant cause for the summertime increased SLAs along the coasts of Sumatra and the Bay of Bengal (Figure 14).Over the Indian subcontinent, convection and precipitation are also enhanced, indicating a stronger summer monsoon situation.However, whether or not the negative decadal IOD events always accompany enhanced decadal convection of Indian summer monsoon requires further investigation.The data record is too short to make a definitive conclusion, and the four-year composite may still contain some inter-annual variability.
While the increased convection in the warm pool reduces the downward shortwave flux and, thus, decreases the sea level through thermosteric effect, the increased precipitation tends to increase the sea level via halosteric effect.The opposing influence of shortwave flux and precipitation may partly account for the compensating effect of the thermal and halosteric sea level on regional scales discussed in Stammer et al. [5].This compensation weakens the contribution of buoyancy flux to regional decadal SLA, which may further contribute to the dominance of winds in causing regional, decadal sea level variability over the IO.

Summary and Discussions
By analyzing observational datasets and the results from large ensemble experiments of two global climate models, we assessed the influences of natural internal variability versus external (natural plus anthropogenic) forcing on the observed multi-decadal trend patterns and decadal SLA in the tropical IO north of 20 • S since the 1960s.Here, we focused on understanding the "regional trend pattern" and decadal SLA after removing the global mean SLR.Consequently, the significant effect of external forcing on the global mean SLR was removed before our analysis.Using Bayesian DLM, we quantified the influences of climate modes (that is, ENSO, IOD, and monsoon; Figure 1) on the IO regional, decadal SLA.The Bayesian DLM can capture the non-stationary impacts of climate modes on the regional sea level and, thus, it consistently outperforms the conventional multiple linear regression model in simulating the observations and in estimating the impacts of each climate mode [33] (Section 2.4).Our major findings are summarized below.
(1) For the 48-yr trend pattern of the sea level (with global SLR removed) from 1958-2005, natural internal variability dominates the external forcing in causing the spatially uneven sea level change patterns over the IO, with natural variability (external forcing) accounting for ~81% (19 ± 2.4%) of the observed falling trend over the Seychelles Islands where both the observed and externally forced trends obtain their maxima (Figure 2).These percentage values are obtained using the average of ORAS4 and SODA to represent the observation and the 100-member mean of the external forcing ensemble, with each member being the ensemble average of 98 members (70% of 140) randomly chosen from the 140 members of the two climate models (100 from MPI and 40 from CESM1), to represent external forcing (Table 1; see Section 3.1 for detailed descriptions on the Monte Carlo method of assessing external forcing and uncertainty).This is because large quantitative or even qualitative differences exist in some regions between ORAS4 and SODA data and between the MPI and CESM1 results (Figures 2-4), even though both reanalysis data produce similar basin-wide patterns.The ORAS4 product agrees the best with the AVISO sea level since 1993 (pattern correlation r = 0.84) compared to SODA (r = 0.48) and the thermosteric sea level (Figure 3), likely due to the ORAS4 assimilated satellite-observed sea level data.The 100-member ensemble mean agrees well with the subsamples of the 40-member ensemble means of the MPI model, even though some quantitative differences remain in some regions.This suggests that 40-member ensemble averages can reasonably damp internal variability in the multi-decadal trend, and the large differences between the MPI and CESM1 simulations result from the model dependence of the sea level response to the same external forcing.(2) For decadal variability, represented by 8 year low-passed SLA from 1962-2011 (both global SLR and regional trend removed), the largest contribution of external forcing (measured by the 1.65STD of the 100-member mean of the external forcing ensemble) is 18 ± 17% of the observed SLA STD (ORAS4 and SODA average) in the Seychelles region and 17 ± 11% in the southeastern basin near the ITF (Figures 4 and 5).Consequently, the observed decadal SLA results predominantly from natural internal variability.While the variances of the decadal SLA from the MPI and CESM1 reach large magnitudes in similar areas, the time series of their LE averaged decadal SLAs for the same areas differ considerably in amplitude and phase, showing weak or negative correlations.These results suggest that the decadal predictions and future projections of the regional sea level can be strongly model-dependent because they can respond very differently to the same external forcing on regional scales.Efforts should be made to understand the causes for the strong model-dependence and identify their biases, with a hope to improve the physical representations of processes in climate models, and ultimately improve their simulations, predictions, and future projections.(3) Decadal SLAs exhibit distinct basin-wide patterns, with the sea level fall in the southwest tropical IO accompanying the sea level rise in the eastern basin (Figures 9 and 10).The strong decadal variability overlies on the multi-decadal falling trend in the western basin, but on the rising trend in the eastern basin, particularly in the coastal regions of Sumatra and the Bay of Bengal.This can increase the stress of SLR for coastal countries such as Bangladesh, whose elevation is near the sea level.Climate modes (the sum of decadal ENSO, IOD, and monsoon) explain a large fraction of the observed decadal sea level variability, regardless of whether the decadal SLAs are based on the annual (12-month) mean winter (November-April) or summer (May-October) seasons, and the impacts of the individual modes exhibit apparent seasonality (Figures 6-12).
During November-April (Figure 11), the STD ratios between the climate-modes-associated and observed decadal SLAs are ~80-86% and their correlation coefficients are r = 0.80-0.89(>95% significance) in the four key regions where the decadal SLAs are large, that is, the Seychelles Islands (Reg.1), the ITF area (Reg.2), west coasts of Sumatra (Reg.3), and the east and north coasts of the Bay of Bengal (Reg.4).The smallest impact occurs over Seychelles (78% and r = 0.80) and the largest impact occurs near the ITF (86% and r = 0.89).Adding the IOSST, which represents part of the external forcing, increases the STD ratios to >90% and r to ≥0.94 in all areas but the Seychelles, where the ratio increases to 83% and r to 0.87 (Figure 11).The decadal ENSO is the largest contributor to the basin-wide pattern and individual areas in both the basin interior and along the coasts (Figures 7 and 9) and it controls the SLA near the ITF area.Decadal variability of the monsoon also has evident effects in all regions except for the ITF, and the influence of the decadal IOD is relatively weak, with its strongest effect occurring over Seychelles.The caveat for the weak IOD influence in the eastern IO coast during boreal winter is that part of the IOD impact is embedded in the effect of ENSO because we removed the ENSO signals from the DMI to ensure the independence of predictors in the DLM equations (Section 2.4).The removal of ENSO from the DMI implies that the ENSO affects IOD but IOD does not affect ENSO.Yet existing studies have suggested the active impact of IOD on ENSO.This procedure may, to some degree, overestimate the ENSO effect and underestimate the IOD effect (Section 3.2.2).
During May-October, the climate modes explain a lower percentage of the observed SLA variance compared to winter, both for the basin-wide EOF patterns and for the individual regions over Seychelles, Sumatra coasts, and Bay of Bengal coasts where the STD ratios between climate-modes-associated and observed SLAs are 67%, 58%, and 63% and correlations are r = 0.76, 0.86, and 0.86, respectively (Figures 8, 10 and 12).Adding IOSST increases the STD ratios to 73%, 75%, and 82% and the correlation to 0.81, 0.93, and 0.95.An exception is for the ITF area where the total effects of the climate modes, which are higher than their corresponding winter values, explain 95% of the observed STD with r = 0.95.Decadal ENSO is the determining factor for the ITF area, with the decadal monsoon having a comparable contribution in the past decade or so (1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011).Over the Seychelles area, decadal ENSO is the largest contributor, with decadal monsoon also having considerable contributions particularly in the past decade.The effect of decadal IOD is also non-negligible.Along the coasts of Sumatra and the Bay, decadal IOD is the largest contributor overall from 1962-2011 and is the controlling factor from the end of the 1980s to early 1990s, while decadal monsoon is the largest contributor in the past decade.(4) Decadal variability of the surface winds over the IO plays a vital role in causing decadal SLAs over most regions of the IO, except for the ITF area where the remote forcing from the Pacific via the ITF is also important.While the Ekman pumping velocity is the major cause for the large SLAs over the Seychelles Islands, the zonal equatorial wind anomalies associated with the Walker Circulation primarily determines the SLAs along the coasts of Sumatra, Java, and the Bay of Bengal for both the winter and summer seasons.The local Ekman pumping velocity near the Sumatra and Java coasts, together with the alongshore winds and Ekman pumping velocity in the eastern Bay during summer, enhance the SLA amplitudes in the coastal regions.These results, which are based on reanalysis data, are supported by the AVISO satellite sea level and CCMP winds for their overlapping periods during recent decades.In addition, satellite-observed convection (that is, OLR) and GPCP precipitation further confirm the wind variations, and they also suggest a compensating effect of thermosteric and halosteric sea level.This is because in regions of increased convection and precipitation, such as the western equatorial IO associated with +IOD, the reduction in the shortwave flux and upper ocean temperature act to decrease thermosteric sea level; meanwhile, the increased precipitation increases the halosteric sea level.This compensation weakens the total effect of buoyancy forcing on regional, decadal SLA.Topographic trapping of the Rossby wave energy by the continental shelf along the west Australian coast [98], together with frictional damping of coastal Kelvin and Rossby waves will also affect coastal SLA.These coastal effects will be investigated in our future research using high-resolution numerical models.
Even though the record length of ~50 year is relatively short, our results from the observational analysis using Bayesian DLM may, to a certain degree, serve as a benchmark for verifying the skills of global climate models in simulating climate modes and their associated decadal sea level patterns over the Indian Ocean.Note that the positive skill of decadal predictions in the first few years results from the observation-based initialization [41,42], which carries climate mode signals that evolve in time via air-sea coupling and interactions among different components of the climate system.Model biases and mean state drift, however, significantly limit the prediction of climate modes' evolution.For this reason, understanding the past sea level behaviors associated with climate modes, as is done in this paper, and using them to infer the near-future evolutions of regional sea level using advanced empirical methods, are as important as the efforts of improving the simulation skills of climate models.The two approaches should be explored in parallel to satisfy the short-and long-term societal needs of sea level predictions on regional scales.

Climate 2018, 6 , 35 Figure 1 .
Figure 1.(a) The 8 year low-passed climate indices based on the annual mean Nino3.4 index, Indian Ocean Dipole (IOD) Mode Index (DMI), and the annual total of all India monsoon rainfall (mm) index, and the Indian Ocean SST (IOSST) index; (b) same as (a) but based on winter (November-April) data; (c) same as (a) but based on summer (May-October) data.The 1962-2011 linear trend is removed and each index is normalized by its standard deviation (STD).Panels (d-f) are the same as (a-c), except for removing the El Niño-Southern Oscillation (ENSO) effect on the DMI, monsoon, and IOSST, and removing the monsoon effect on the DMI, using Bayesian dynamical linear model (DLM; see Section 2.4).Consequently, the predictors shown in (d-f) are independent of each other.

Figure 1 .
Figure 1.(a) The 8 year low-passed climate indices based on the annual mean Nino3.4 index, Indian Ocean Dipole (IOD) Mode Index (DMI), and the annual total of all India monsoon rainfall (mm) index, and the Indian Ocean SST (IOSST) index; (b) same as (a) but based on winter (November-April) data; (c) same as (a) but based on summer (May-October) data.The 1962-2011 linear trend is removed and each index is normalized by its standard deviation (STD).Panels (d-f) are the same as (a-c), except for removing the El Niño-Southern Oscillation (ENSO) effect on the DMI, monsoon, and IOSST, and removing the monsoon effect on the DMI, using Bayesian dynamical linear model (DLM; see Section 2.4).Consequently, the predictors shown in (d-f) are independent of each other.

Figure 2 .
Figure 2. The linear trends of regional sea level over the tropical Indian Ocean from 1958-2005.Monthly data are used, with the 1958-2005 monthly climatology and global mean sea level rise (SLR) time series removed from each dataset before the trend calculation.(a) The ECMWF Ocean Reanalysis System 4 (ORAS4) reanalysis; (b) the Simple Ocean Data Assimilation (SODA) reanalysis; (c) the world ocean atlas 2013 (WOA13) upper 700 m thermosteric sea level; (d) Ishii upper 700 m thermosteric sea level; (e) the Max-Planck-Institute of meteorology (MPI) 100-member ensemble mean dynamic sea level; (f) the National Center for Atmospheric Research (NCAR) community earth system model version 1 (CESM1) 40-member ensemble mean dynamic sea level; (g) MPI ensemble mean for members 1-40; (h) MPI ensemble mean for members 41-80.The two boxes show the maximum sea level fall and rise areas in both the ORAS4 and SODA dynamic sea level, with Reg. 1 located in the thermocline ridge of the southwest tropical basin (50° E-85° E, 17° S-5° S) and Reg. 2, eastern equatorial Indian Ocean (90° E-103° E, 7° S-7° N).The number in each box shows the trend value of SLA averaged in the box.Trend values exceeding 95% significance are shown in color contours, and those below 95% are shown in white.

Figure 2 .
Figure 2. The linear trends of regional sea level over the tropical Indian Ocean from 1958-2005.Monthly data are used, with the 1958-2005 monthly climatology and global mean sea level rise (SLR) time series removed from each dataset before the trend calculation.(a) The ECMWF Ocean Reanalysis System 4 (ORAS4) reanalysis; (b) the Simple Ocean Data Assimilation (SODA) reanalysis; (c) the world ocean atlas 2013 (WOA13) upper 700 m thermosteric sea level; (d) Ishii upper 700 m thermosteric sea level; (e) the Max-Planck-Institute of meteorology (MPI) 100-member ensemble mean dynamic sea level; (f) the National Center for Atmospheric Research (NCAR) community earth system model version 1 (CESM1) 40-member ensemble mean dynamic sea level; (g) MPI ensemble mean for members 1-40; (h) MPI ensemble mean for members 41-80.The two boxes show the maximum sea level fall and rise areas in both the ORAS4 and SODA dynamic sea level, with Reg. 1 located in the thermocline ridge of the southwest tropical basin (50 • E-85 • E, 17 • S-5 • S) and Reg. 2, eastern equatorial Indian Ocean (90 • E-103 • E, 7 • S-7 • N).The number in each box shows the trend value of SLA averaged in the box.Trend values exceeding 95% significance are shown in color contours, and those below 95% are shown in white.

Figure 3 .
Figure 3.The linear trend map of sea level (mm year −1 ; global mean SLR removed) over the Indian Ocean from 1993-2010, when all datasets are available.(a) The French Archiving, Validation, and Interpolation of Satellite Oceanographic Data (AVISO) satellite altimeter data; (b) ORAS4 reanalysis; (c) SODA reanalysis; (d) WOA13 upper 700 m thermosteric sea level; (e) Ishii upper 700 m thermosteric sea level.Spatial pattern correlation coefficient between AVISO and each product is also plotted, with r = 0.84 for AVISO/ORAS4, 0.48 for AVISO/SODA, 0.48 for AVISO/WOA13, and 0.73 for AVISO/Ishii, with all coefficients exceeding 99% significance.

Figure 3 .
Figure 3.The linear trend map of sea level (mm year −1 ; global mean SLR removed) over the Indian Ocean from 1993-2010, when all datasets are available.(a) The French Archiving, Validation, and Interpolation of Satellite Oceanographic Data (AVISO) satellite altimeter data; (b) ORAS4 reanalysis; (c) SODA reanalysis; (d) WOA13 upper 700 m thermosteric sea level; (e) Ishii upper 700 m thermosteric sea level.Spatial pattern correlation coefficient between AVISO and each product is also plotted, with r = 0.84 for AVISO/ORAS4, 0.48 for AVISO/SODA, 0.48 for AVISO/WOA13, and 0.73 for AVISO/Ishii, with all coefficients exceeding 99% significance.

Figure 4 .
Figure 4. (a)The standard deviation (STD) map of decadal SLA (with global mean SLR and regional trend removed) from ORAS4 reanalysis over the tropical Indian Ocean for the 1962-2011 period.First, the global SLR was removed from the monthly sea level data for 1958-2015; then monthly climatology and the linear trend at each grid point were removed; finally, a Lanczos 8 year lowpass filter was applied; the first and last four years of data(1958-1961 and 2012-2015)  are excluded to remove the endpoint effect of the filter; the two boxes show the maximum STD regions, with Reg. 1 representing the thermocline ridge (50° E-85° E, 17° S-5° S) shown in Figures2 and 4, and Reg. 2, the near ITF region (105° E-120° E, 20° S-12° S); (b) same as (a) but for the ORAS4 for the 1962-2001 period; (c) same as (b) but for the SODA reanalysis data from 1962-2001; (d) same as (b) but for the CESM1 40member ensemble mean from 1962-2001; (e) same as (d) but for the MPI 100-member ensemble mean from 1962-2001; (f) same as (e) but for the MPI 40-member ensemble mean; (g) the time series of decadal SLA averaged over Reg.1; (h) same as (g) but for the decadal SLA averaged over Reg. 2. Units: mm.

Figure 4 .
Figure 4. (a)The standard deviation (STD) map of decadal SLA (with global mean SLR and regional trend removed) from ORAS4 reanalysis over the tropical Indian Ocean for the 1962-2011 period.First, the global SLR was removed from the monthly sea level data for 1958-2015; then monthly climatology and the linear trend at each grid point were removed; finally, a Lanczos 8 year lowpass filter was applied; the first and last four years of data(1958-1961 and 2012-2015)  are excluded to remove the endpoint effect of the filter; the two boxes show the maximum STD regions, with Reg. 1 representing the thermocline ridge (50 • E-85 • E, 17 • S-5 • S) shown in Figures2 and 4, and Reg. 2, the near ITF region (105 • E-120 • E, 20 • S-12 • S); (b) same as (a) but for the ORAS4 for the 1962-2001 period; (c) same as (b) but for the SODA reanalysis data from 1962-2001; (d) same as (b) but for the CESM1 40-member ensemble mean from 1962-2001; (e) same as (d) but for the MPI 100-member ensemble mean from 1962-2001; (f) same as (e) but for the MPI 40-member ensemble mean; (g) the time series of decadal SLA averaged over Reg.1; (h) same as (g) but for the decadal SLA averaged over Reg. 2. Units: mm.

Figure 5 .
Figure 5.The time series of the 8 year lowpass filtered decadal SLA (mm) averaged over the Seychelles Island area in the west (Reg. 1 of Figure 4) and near the ITF area in the east (Reg. 2 of Figure 4) from 1962-2001.(a) The time series of the 100 members that represent external forcing, with each member being the ensemble average of the randomly chosen N members from the total MPI and CESM (140 members) averaged in Reg. 1 (thin curves), the effect of the external forcing represented by the 100member ensemble mean (thick black curve), and the observation represented by the average of ORAS4 and SODA (red curve); the 18% ± 17% are the STD ratios between the external-forced and observed SLA and its uncertainty (1.65STD of spread, shown by the gray shading); (b,c) are the same as (a) but for the MPI 100 members and the CESM 40 members, respectively; (d-f) are the same as (ac) except that they are averaged in Reg. 2. See the last paragraph of Section 3.1 for a detailed discussion on the Monte Carlo method to obtain the 100 members of external forcing and the choice of N.

Figure 5 .
Figure 5.The time series of the 8 year lowpass filtered decadal SLA (mm) averaged over the Seychelles Island area in the west (Reg. 1 of Figure 4) and near the ITF area in the east (Reg. 2 of Figure 4) from 1962-2001.(a) The time series of the 100 members that represent external forcing, with each member being the ensemble average of the randomly chosen N members from the total MPI and CESM (140 members) averaged in Reg. 1 (thin curves), the effect of the external forcing represented by the 100-member ensemble mean (thick black curve), and the observation represented by the average of ORAS4 and SODA (red curve); the 18 ± 17% are the STD ratios between the external-forced and observed SLA and its uncertainty (1.65STD of spread, shown by the gray shading); (b,c) are the same as (a) but for the MPI 100 members and the CESM 40 members, respectively; (d-f) are the same as (a-c) except that they are averaged in Reg. 2. See the last paragraph of Section 3.1 for a detailed discussion on the Monte Carlo method to obtain the 100 members of external forcing and the choice of N.

Figure 6 .
Figure 6.(a) The standard deviation map of ORAS4 8 year low-passed SLA for the 1962-2011 period based on the yearly mean (all 12 months) sea level data; it is the same as Figure 4a except for adding boxed region R3 in the eastern equatorial basin (90° E-102° E, 5° S-5° N) and boxed region R4 in the Bay of Bengal coast.R4 consists of an area east of 92° E but north of 10° N and an area west of 92° E but north of 18° N; (b) same as (a) but for SLA associated with all climate modes (ENSO + IOD + monsoon) plus the effect of IOSST, which to some extent represents external forcing; (c) same as (b) but only for climate modes (ENSO + IOD + monsoon), and it equals the IOSST effect subtracted from panel (b); (d) same as (b) but for the residual, difference for (a,b), measuring the effects that cannot be explained by the climate modes and IOSST; (e-h) are the same as (b) but for the individual effect of decadal ENSO, IOD, monsoon, and IOSST.Units: mm.The boxed regions have large variances and are used to calculate the time series later.

Figure 6 .
Figure 6.(a) The standard deviation map of ORAS4 8 year low-passed SLA for the 1962-2011 period based on the yearly mean (all 12 months) sea level data; it is the same as Figure 4a except for adding boxed region R3 in the eastern equatorial basin (90 • E-102 • E, 5 • S-5 • N) and boxed region R4 in the Bay of Bengal coast.R4 consists of an area east of 92 • E but north of 10 • N and an area west of 92 • E but north of 18 • N; (b) same as (a) but for SLA associated with all climate modes (ENSO + IOD + monsoon) plus the effect of IOSST, which to some extent represents external forcing; (c) same as (b) but only for climate modes (ENSO + IOD + monsoon), and it equals the IOSST effect subtracted from panel (b); (d) same as (b) but for the residual, difference for (a,b), measuring the effects that cannot be explained by the climate modes and IOSST; (e-h) are the same as (b) but for the individual effect of decadal ENSO, IOD, monsoon, and IOSST.Units: mm.The boxed regions have large variances and are used to calculate the time series later.

Figure 7 .
Figure 7. (a) The standard deviation map of the ORAS4 8 year low-passed SLA for the 1962-2011 period of the northern hemisphere winter season, based on the November-April averaged data of each year; (b) same as (a) but for SLA associated with all climate modes (ENSO + IOD + monsoon) plus the effect of IOSST which, to some extent, represents external forcing; (c) same as (b) but only for climate modes (ENSO + IOD + monsoon) and it equals the IOSST effect subtracted from panel (b); (d) same as (b) but for the residual, difference for (a,b), measuring the effects that cannot be explained by the climate modes and IOSST; (e-h) are the same as (b) but for the individual effect of decadal ENSO, IOD, monsoon, and IOSST.Units: mm.The figure is the same as Figure 6 except for northern winter season.

Figure 7 .
Figure 7. (a) The standard deviation map of the ORAS4 8 year low-passed SLA for the 1962-2011 period of the northern hemisphere winter season, based on the November-April averaged data of each year; (b) same as (a) but for SLA associated with all climate modes (ENSO + IOD + monsoon) plus the effect of IOSST which, to some extent, represents external forcing; (c) same as (b) but only for climate modes (ENSO + IOD + monsoon) and it equals the IOSST effect subtracted from panel (b); (d) same as (b) but for the residual, difference for (a,b), measuring the effects that cannot be explained by the climate modes and IOSST; (e-h) are the same as (b) but for the individual effect of decadal ENSO, IOD, monsoon, and IOSST.Units: mm.The figure is the same as Figure 6 except for northern winter season.

Climate 2018, 6 , 35 Figure 8 .
Figure 8.(a) The standard deviation map of the ORAS4 8 year low-passed SLA for the 1962-2011 period of the northern hemisphere summer season, based on the May-October averaged data of each year; (b) same as (a) but for SLA associated with all climate modes (ENSO + IOD + monsoon) plus the effect of IOSST which, to some extent, represents external forcing; (c) same as (b) but only for climate modes (ENSO + IOD + monsoon) and it equals the IOSST effect subtracted from panel (b); (d) same as (b) but for the residual, difference for (a,b), measuring the effects that cannot be explained by the climate modes and IOSST; (e-h) are the same as (b) but for the individual effect of decadal ENSO, IOD, monsoon, and IOSST.Units: mm.The figure is the same as Figure 6 except for the northern summer.

Figure 8 .
Figure 8.(a) The standard deviation map of the ORAS4 8 year low-passed SLA for the 1962-2011 period of the northern hemisphere summer season, based on the May-October averaged data of each year; (b) same as (a) but for SLA associated with all climate modes (ENSO + IOD + monsoon) plus the effect of IOSST which, to some extent, represents external forcing; (c) same as (b) but only for climate modes (ENSO + IOD + monsoon) and it equals the IOSST effect subtracted from panel (b); (d) same as (b) but for the residual, difference for (a,b), measuring the effects that cannot be explained by the climate modes and IOSST; (e-h) are the same as (b) but for the individual effect of decadal ENSO, IOD, monsoon, and IOSST.Units: mm.The figure is the same as Figure 6 except for the northern summer.

Climate 2018, 6 , 35 Figure 9 .
Figure 9.The leading Empirical Orthogonal Function (EOF1) of the 8 year low-passed November-April mean SLA over the tropical Indian Ocean for the period of 1962-2011 from (a) ORAS4 data; (b) Bayesian DLM simulation from the sum of all modes and IOSST (ENSO + IOD + monsoon + IOSST); and (c) DLM simulation from all modes (ENSO + IOD + monsoon); (d) the principal component of EOF1 (PC1) for ORAS4 (black), all modes + IOSST (blue), and all modes (red); (e-g) are the same as (c) except for the DLM simulation from ENSO, IOD, and monsoon, respectively; (h) Same as (d) except for the PC1 of ENSO (black), IOD (red), and monsoon (blue).The correlation coefficients shown in (d,h) are with the observed ORAS4 PC1 (the black curve of panel (d)).

Figure 9 .
Figure 9.The leading Empirical Orthogonal Function (EOF1) of the 8 year low-passed November-April mean SLA over the tropical Indian Ocean for the period of 1962-2011 from (a) ORAS4 data; (b) Bayesian DLM simulation from the sum of all modes and IOSST (ENSO + IOD + monsoon + IOSST); and (c) DLM simulation from all modes (ENSO + IOD + monsoon); (d) the principal component of EOF1 (PC1) for ORAS4 (black), all modes + IOSST (blue), and all modes (red); (e-g) are the same as (c) except for the DLM simulation from ENSO, IOD, and monsoon, respectively; (h) Same as (d) except for the PC1 of ENSO (black), IOD (red), and monsoon (blue).The correlation coefficients shown in (d,h) are with the observed ORAS4 PC1 (the black curve of panel (d)).

35 Figure 10 .
Figure 10.The EOF1 of the 8 year low-passed May-October mean SLA over the tropical Indian Ocean for the period of 1962-2011 from (a) ORAS4 data; (b) DLM simulation from the sum of all modes and IOSST (ENSO + IOD + monsoon + IOSST); and (c) DLM simulation from all modes (ENSO + IOD + monsoon); (d) PC1 for ORAS4 (black), all modes + IOSST (blue) and all modes (red); (e-g) are the same as (c) except for the DLM simulation from ENSO, IOD, and monsoon, respectively; (h) same as (d) except for the PC1 of ENSO (black), IOD (red), and monsoon (blue).The correlation coefficients shown in (d,h) are with the observed ORAS4 PC1 (the black curve of panel (d)).This figure is the same as Figure 9 except for the northern summer (May-October mean).

Figure 10 .
Figure 10.The EOF1 of the 8 year low-passed May-October mean SLA over the tropical Indian Ocean for the period of 1962-2011 from (a) ORAS4 data; (b) DLM simulation from the sum of all modes and IOSST (ENSO + IOD + monsoon + IOSST); and (c) DLM simulation from all modes (ENSO + IOD + monsoon); (d) PC1 for ORAS4 (black), all modes + IOSST (blue) and all modes (red); (e-g) are the same as (c) except for the DLM simulation from ENSO, IOD, and monsoon, respectively; (h) same as (d) except for the PC1 of ENSO (black), IOD (red), and monsoon (blue).The correlation coefficients shown in (d,h) are with the observed ORAS4 PC1 (the black curve of panel (d)).This figure is the same as Figure 9 except for the northern summer (May-October mean).

Figure 11 .
Figure 11.(a) The November-April averaged time series of the 8 year low-passed SLA averaged in the thermocline ridge region (Reg. 1 of Figures 6-10) from 1962-2011 from ORAS4 (thick solid black), the sum of climate modes' effects (ENSO + IOD + monsoon; thin solid black), all modes plus IOSST (green), and residual (black dotted), which is (ORAS4 − (all modes + IOSST)); linear trend for 1962-2011 is removed from each curve; (b) same as (a) except for contributions from all modes (black), decadal ENSO (red), decadal IOD (blue), and decadal monsoon (purple); panels (c-h) are the same as (a,b) except for the ITF area (Reg. 2 of Figures 6-10), eastern equatorial basin (Reg.3 of Figures 6-10), and the Bay of Bengal coasts (Reg.4 of Figures 6-10), respectively.The correlation coefficients r in each panel of the left column are between the observed and DLM simulated SLA from climate modes (black) and climate modes + IOSST (green), and those in the right column show the r between each climate mode and the sum of all climate modes.The two horizontal lines in (a) show the ±0.5 STD of the ORAS4 observed SLA, which serve as a reference for the composites of Section 3.3 below.

Figure 11 .
Figure 11.(a) The November-April averaged time series of the 8 year low-passed SLA averaged in the thermocline ridge region (Reg. 1 of Figures 6-10) from 1962-2011 from ORAS4 (thick solid black), the sum of climate modes' effects (ENSO + IOD + monsoon; thin solid black), all modes plus IOSST (green), and residual (black dotted), which is (ORAS4 − (all modes + IOSST)); linear trend for 1962-2011 is removed from each curve; (b) same as (a) except for contributions from all modes (black), decadal ENSO (red), decadal IOD (blue), and decadal monsoon (purple); panels (c-h) are the same as (a,b) except for the ITF area (Reg. 2 of Figures 6-10), eastern equatorial basin (Reg.3 of Figures 6-10), and the Bay of Bengal coasts (Reg.4 of Figures 6-10), respectively.The correlation coefficients r in each panel of the left column are between the observed and DLM simulated SLA from climate modes (black) and climate modes + IOSST (green), and those in the right column show the r between each climate mode and the sum of all climate modes.The two horizontal lines in (a) show the ±0.5 STD of the ORAS4 observed SLA, which serve as a reference for the composites of Section 3.3 below.

Figure 12 .
Figure 12.(a) The May-October averaged time series of the 8 year low-passed SLA averaged in Reg. 1 of Figures 6-10 from 1962-2011 from ORAS4 (thick solid black), the sum of climate modes' effects (ENSO + IOD + monsoon; thin solid black), all modes plus IOSST (green), and residual (black dotted), which is (ORAS4 − (all modes + IOSST)]; the linear trend for 1962-2011 is removed from each curve; (b) same as (a) except for the contributions from all modes (black), decadal ENSO (red), decadal IOD (blue), and decadal monsoon (purple); panels (c-h) are the same as (a,b) except for the ITF area (Reg. 2 of Figures 6-10), eastern equatorial basin (Reg.3 of Figures 6-10), and Bay of Bengal coasts (Reg.4 of Figures 6-10), respectively.The correlation coefficients r in each panel of the left column are between the observed and DLM simulated SLA from climate modes (black) and climate modes + IOSST (green), and those in the right column show the r between each climate mode and the sum of all climate modes.The two horizontal lines in (g) show the ±0.5 STD of ORAS4 observed SLA, which serve as a reference for the composites of Section 3.3.The figure is the same as Figure 11 except for the northern summer.

Figure 12 .
Figure 12.(a) The May-October averaged time series of the 8 year low-passed SLA averaged in Reg. 1 of Figures 6-10 from 1962-2011 from ORAS4 (thick solid black), the sum of climate modes' effects (ENSO + IOD + monsoon; thin solid black), all modes plus IOSST (green), and residual (black dotted), which is (ORAS4 − (all modes + IOSST)]; the linear trend for 1962-2011 is removed from each curve; (b) same as (a) except for the contributions from all modes (black), decadal ENSO (red), decadal IOD (blue), and decadal monsoon (purple); panels (c-h) are the same as (a,b) except for the ITF area (Reg. 2 of Figures 6-10), eastern equatorial basin (Reg.3 of Figures 6-10), and Bay of Bengal coasts (Reg.4 of Figures 6-10), respectively.The correlation coefficients r in each panel of the left column are between the observed and DLM simulated SLA from climate modes (black) and climate modes + IOSST (green), and those in the right column show the r between each climate mode and the sum of all climate modes.The two horizontal lines in (g) show the ±0.5 STD of ORAS4 observed SLA, which serve as a reference for the composites of Section 3.3.The figure is the same as Figure 11 except for the northern summer.

Figure 13 .
Figure 13.The composite analyses for the November-April mean, 8 year low-passed decadal SLA of November-April mean (color contours), the Japan Meteorological Agency 55-year Reanalysis (JRA55) surface wind stress (arrows; units: ×0.2 dyn cm −2 ), and the Ekman pumping velocity,  (line contours, positive (negative) values are in white (black), beginning at 1 (−1) with an interval of 2; zero contours are omitted; units: ×10 −7 m s −1 ), based on the years when the magnitudes of the ORAS4 decadal SLA in the thermocline ridge area (Reg. 1 of Figures 6-10) exceeded the −0.5 STD line shown in Figure 11a for: (a) the total SLA and wind anomalies; (b) the sum of the DLM results from all modes plus IOSST; (c) the sum of the DLM results from all modes; (d) the residual, which is panel (a) minus panel (b); (e-h) are the same as (c) but for the individual mode ENSO, IOD, monsoon, and IOSST, respectively.Note that the  values within 5° S-5° N are not shown.

Figure 13 .
Figure 13.The composite analyses for the November-April mean, 8 year low-passed decadal SLA of November-April mean (color contours), the Japan Meteorological Agency 55-year Reanalysis (JRA55) surface wind stress (arrows; units: ×0.2 dyn cm −2 ), and the Ekman pumping velocity, w e (line contours, positive (negative) values are in white (black), beginning at 1 (−1) with an interval of 2; zero contours are omitted; units: × 10 −7 m s −1 ), based on the years when the magnitudes of the ORAS4 decadal SLA in the thermocline ridge area (Reg. 1 of Figures 6-10) exceeded the −0.5 STD line shown in Figure 11a for: (a) the total SLA and wind anomalies; (b) the sum of the DLM results from all modes plus IOSST; (c) the sum of the DLM results from all modes; (d) the residual, which is panel (a) minus panel (b); (e-h) are the same as (c) but for the individual mode ENSO, IOD, monsoon, and IOSST, respectively.Note that the w e values within 5 • S-5 • N are not shown.

Figure 14 .
Figure 14.The composite analyses for the May-October mean, 8 year low-passed decadal SLA (color contours), the JRA55 surface wind stress (arrows; units: ×0.2 dyn cm −2 ), and the Ekman pumping velocity,  (line contours, positive (negative) values are in white (black), beginning at 1 (−1) with an interval of 2; zero contours are omitted; units: ×10 −7 m s −1 ), based on the years when the SLA values in the eastern Bay of Bengal (Reg. 4 of Figures 6-10) exceeded the +0.5 STD line shown in Figure 12g for: (a) the total SLA and wind anomalies; (b) the sum of the DLM results from all modes plus IOSST; (c) the sum of the DLM results from all modes; (d) the Residual, which is panel (a) minus panel (b); (e-h) are the same as (c) but for the individual mode ENSO, IOD, monsoon, and IOSST, respectively.The  values within 5° S-5° N are omitted.

Figure 14 .
Figure 14.The composite analyses for the May-October mean, 8 year low-passed decadal SLA (color contours), the JRA55 surface wind stress (arrows; units: ×0.2 dyn cm −2 ), and the Ekman pumping velocity, w e (line contours, positive (negative) values are in white (black), beginning at 1 (−1) with an interval of 2; zero contours are omitted; units: × 10 −7 m s −1 ), based on the years when the SLA values in the eastern Bay of Bengal (Reg. 4 of Figures 6-10) exceeded the +0.5 STD line shown in Figure 12g for: (a) the total SLA and wind anomalies; (b) the sum of the DLM results from all modes plus IOSST; (c) the sum of the DLM results from all modes; (d) the Residual, which is panel (a) minus panel (b); (e-h) are the same as (c) but for the individual mode ENSO, IOD, monsoon, and IOSST, respectively.The w e values within 5 • S-5 • N are omitted.

Figure 15 .
Figure 15.(a) The composite for the November-April mean, 8 year low-passed ORAS4 SLA (color contours), and the JRA55 surface wind stress (arrows; units: ×0.2 dyn cm −2 ) for years 1998-2002 when the magnitudes of ORAS4 SLA of Reg. 1 (Figures 6-10) exceeded −0.5 STD (see Figure 11a); it is the same as Figure 13a except for the years within 1993-2011 when all the datasets are available; linear trend for 1993-2011 is removed from each field; (b) same as (a) but with no lowpass filtering; (c) same as (b) but for the AVISO SLA and cross-calibrated multiplatform (CCMP) satellite winds with no filtering; (d) same as (c) but with color contours showing the Outgoing Longwave Radiation (OLR) anomaly (white (red) line contours showing the Global Precipitation Climatology Project (GPCP) negative (positive) precipitation anomaly, beginning at ±0.2 mm day −1 with an interval of 0.4; the zero contour suppressed); the OLR and precipitation anomalies are shown over both ocean and land; panels (e-h) are the same as (a-d) except for the May-October mean composites for years 1993 and 1999-2001 when the magnitudes of ORAS4 decadal SLA of Reg. 4 exceeded +0.5STD (see Figure 12g).

Figure 15 .
Figure 15.(a) The composite for the November-April mean, 8 year low-passed ORAS4 SLA (color contours), and the JRA55 surface wind stress (arrows; units: ×0.2 dyn cm −2 ) for years 1998-2002 when the magnitudes of ORAS4 SLA of Reg. 1 (Figures 6-10) exceeded −0.5 STD (see Figure 11a); it is the same as Figure 13a except for the years within 1993-2011 when all the datasets are available; linear trend for 1993-2011 is removed from each field; (b) same as (a) but with no lowpass filtering; (c) same as (b) but for the AVISO SLA and cross-calibrated multiplatform (CCMP) satellite winds with no filtering; (d) same as (c) but with color contours showing the Outgoing Longwave Radiation (OLR) anomaly (white (red) line contours showing the Global Precipitation Climatology Project (GPCP) negative (positive) precipitation anomaly, beginning at ±0.2 mm day −1 with an interval of 0.4; the zero contour suppressed); the OLR and precipitation anomalies are shown over both ocean and land; panels (e-h) are the same as (a-d) except for the May-October mean composites for years 1993 and 1999-2001 when the magnitudes of ORAS4 decadal SLA of Reg. 4 exceeded +0.5STD (see Figure 12g).