Dual response of Arabian Sea cyclones and strength of Indian monsoon to Southern Atlantic Ocean

Variability and trends of the south Asian monsoon at different time scales makes the region susceptible to climate-related natural disasters such as droughts and floods. Because of its importance, different studies have examined the climatic factors responsible for the recent changes in monsoon strength. Here, using observations and climate model experiments we show that monsoon strength is driven by the variations of south Atlantic Ocean sea surface temperature (SASST). The mechanism by which SASST is modulating the monsoon could be explained through the classical Matsuno-Gill response, leading to changes in the characteristics of vertical wind shear in the Arabian Sea. The decline in the vertical wind shear to the warming of SASST is associated with anomalous lower (upper)-level easterlies (westerlies). This further leads to a strong increase in the frequency of the Arabian Sea cyclones; and also prohibits the transport of moisture to the Indian landmass, which eventually reduces the strength of monsoon. The conditions in the SASST which drove these responses are aggravated by greenhouse gas emission, revealing the prominent role played by anthropogenic warming. If, with proper mitigation, these emissions are not prevented, further increases in the SASST is expected to result in increased Arabian sea cyclones and reduced monsoon strength.


Introduction
The South Asian summer monsoon rainfall, which occurs during June to September (JJAS), has a significant impact on nearly one-half of the world's population. Variation in strength of the monsoon, due to both natural and anthropogenic influences, has a tremendous impact towards increasing the natural disaster risks (Sharma et al. 2020;Vittal et al. 2020a). Considering this, there has been extensive work on understanding the past and future changes in the monsoon using both observed data and climate model simulationsyet many aspects of the monsoon remain poorly understood, leading to uncertainties in future predictions (Shashikanth et al. 2018). Studies show that the south Asian monsoon rainfall (Kitoh et al. 2013), Hadley circulation over the Indian ocean (Roxy et al. 2015) and Walker circulations over the Indo-pacific region (Vecchi et al. 2006) have weakened post-1950. To explain this, multiple hypotheses have been proposed in literature, viz., increase in aerosol emissions (Bollasina et al. 2011), Western Indian Ocean warming (Roxy et al. 2015) and changes in land use land cover (Paul et al. 2016).
The hypothesis put forward for explaining the declining monsoons are mostly related to weakening of the monsoon south-westerly winds. However, along with the low-level westerlies, the formation of depressions over Bay of Bengal (BoB), which transport enormous amounts of moisture, also contribute significantly to the mean monsoon rainfall during summer (Vishnu et al. 2016). The decreasing trend in the frequency of monsoon depressions in the BoB is related to anomalous moisture convergence over the western equatorial Indian Ocean caused by an increase in sea surface temperature (Vishnu et al. 2016). Another reason for the reduction in monsoon depressions over the BoB may be linked with the episodic surges of moisture transport from the Arabian Sea (AS) (Roxy et al. 2017). Precipitation usually exceeds evaporation in BoB region during summer season and this is maintained by a substantial moisture influx from the AS (Mooley 1973). However, during recent past, the moisture influx from the AS has been reducing significantly which in turn reduced the genesis of monsoon depressions (Roxy et al. 2017). Moreover, while analysing an increase in extreme precipitation over Central India, Roxy et al. (2017) identified that the primary source of moisture is from the Arabian Sea. They quantify that around 36% of total moisture comes from the AS, while 26% and 9% arise from the BoB and the central Indian Ocean, respectively.
Thus, monsoon depressions in the BoB are reducing and the dynamics of the Arabian Sea has played a major role in doing so. In addition, the moisture supply from the AS to the Indian rainfall is higher when compared to BoB. However, the climate system over AS has changed significantly during recent past. For example,  and Murakami et al. (2017) showed that the frequency of both pre-(May-June) and post-monsoon (October-December) AS tropical cyclones increased during last decade due to anthropogenic warming. In such a scenario, AS cyclones could inflict some 'collateral damage' on the progress of monsoon over mainland India and the monsoon may remain weak until the storm dissipates, specifically during the month of June when the onset of monsoon happen for most regions of India. Although these cyclones can increase the intensities of extreme precipitation over India (as reported in Roxy et al. 2017), it may weaken the mean monsoons, a mechanism that has not been explored in the previous studies. A recent study by Zhang and Villarini (2019), however, showed that the forcing of the south Atlantic Ocean surface temperature (SASST) can drive changes in the AS cyclones. Nonetheless, the study is conducted during the pre-monsoon season. Moreover, the SASST also significantly alters the monsoon rainfall strengths over India during the monsoon season (June-September) by the Matsuno-Gill responses (Kucharski et al. 2007(Kucharski et al. , 2008(Kucharski et al. , 2009Pottapinjara et al. 2014;Sabeerali et al. 2019a). The responses of large-scale circulation to SASST forcing feature a strong upper-level cyclone and low-level anticyclone in the Indian ocean basin, including AS and BoB (Kucharski et al. 2009;Zhang and Villarini 2019), which could lead to noticeable impacts on vertical wind shear over the basin, which is crucial to cyclonic activity. Therefore, the main objective of the present study is to examine whether and how the SASST can modulate the frequency of the AS cyclones using both the observations and climate model perturbation simulations, specifically during the monsoon season, and also explore changing dynamics of the Indian summer monsoon in response to SASST anomalous warming. We then extend our analysis with the use of suite of climate model simulations from the Coupled Model Intercomparison Project phase 5 (CMIP5) to understand how the SASST would change at the end of this century, and underpin the role of the anthropogenic warming in exacerbating such changes.

Data and methods
We procure a long-term precipitation data for Asian monsoon region from the Climate Research Unit (CRU) (Mitchell and Jones 2005) for the period 1951-2017; with a spatial resolution of 0.5° at a monthly temporal scale. The monthly sea surface temperature (SST) data for the period 1870-2018 are obtained from HadISST dataset (Rayner et al. 2014) provided by the Met Office Hadley Centre. Details of cyclone tracks and its intensity for the AS for the period 1979-2010 is obtained from India Meteorological Department (IMD), which is a nodal government agency that provides weather services related to cyclones in India that classifies the pressure system from depression to very severe Cyclonic storms. Further, the large-scale atmospheric information, such as low-and high-level wind, is obtained from the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) reanalysis project (NCEP-NCAR-1) (Kalnay et al. 1996) for the period 1951-2017 available at 2.5° resolution. To examine the characteristics of the mechanism responsible for the changing characteristics of monsoon strength and AS cyclones, we procure the state-of-the-art global climate model simulations from the Coupled Model Intercomparison Project phase 5 (CMIP5, (Taylor et al. 2012); detailed description is provided in Table 1. To quantify the effect of human activities in the past , three types of monthly forcings are analysed: (1) natural-only (HistNat), (2) green-house emission-only (HistGhG), and historical (Hist). Here, the HistNat simulations contain only the effects of natural forcing (e.g., changes in solar radiation, volcanic eruptions), and the HistGhG consist of only the anthropogenic (i.e., GHG concentrations) effects. Hist, on the other hand, considers both natural and anthropogenic effects. To assess possible future climate scenarios, we procure two Representative Concentration Pathway (RCP) scenarios (4.5 and 8.5), which are available for 2006-2100. where b is determined using the linear regression between SST at a respective grid and Niño 3.4 index.

Moisture transport analysis
The amount of moisture transported to the Indian continent, which is surrounded by oceans, is directly linked to the rainfall activity over the region (Levine and Turner 2012;Pathak et al. 2017;Roxy et al. 2017). Thus, we analyse the condition of moisture transport during our identified largescale atmospheric phenomena in terms of integrated water vapour transport (IVT), the key link between source and sink regions of moisture (Lavers et al. 2015). We use specific humidity and the meridional and zonal wind components from the surface to mid-troposphere pressure levels to calculate the magnitude of IVT (Lavers et al. 2012), where q is the specific humidity (kg kg −1 ), u and v are the zonal and meridional wind components (ms −1 ), respectively, g is the acceleration due to gravity (ms −2 ), and dp is the pressure level difference.

Dynamic recycling model
The relative contribution of evaporated water from different ocean sources to the precipitation over the Indian subcontinent (sink) is also analysed using Dynamic Recycling Model [see Martinez and Dominguez (2014) and Pathak et al. (2017) for detailed information]. In this model, the conservation of the total column water vapour is solved by using a classical Lagrangian approach. We use six-hourly fields of vertical integral of eastward and northward water vapor flux to obtain the moisture-weighted wind component in two horizontal directions [see Martinez and Dominguez (2014) and Pathak et al. (2017)] to study the atmospheric moisture transport.

Climate model perturbation experiments
To perform the perturbation experiments, we use atmospheric general circulation model developed by International Centre for Theoretical Physics (ICTP AGCM) (Molteni 2003;Kucharski et al. 2013;Vittal et al. 2020c). Here, we use the latest version (v41) of this model, which consists of 8 vertical levels along-with the horizontal spectral truncation set at T30 ( ∼ 3.75° spatial resolution). This atmospheric model consists of physically based parameterizations of large-scale condensation, shallow and deep convection, short-wave and long-wave radiation, surface fluxes of momentum, heat and moisture, and vertical diffusion. To evaluate the impacts of the SST anomaly on Asian monsoon and large-scale circulation, we performed experiments which are integrated for 30 years. Initially, we start the experiment by prescribing the seasonal climatology of SST based on HadISST (CLIM). Further, to evaluate the response of SST anomaly, which could influence the changing characteristics of both the monsoon strength and AS cyclone, the observed composite anomaly of SASST during positive and negative phase for summer monsoon season are superimposed on the seasonal climatology, keeping SST identical to the CLIM experiment for other months. The subtraction of CLIM experiment with the SST anomaly forced experiment then provides the response of large-scale atmospheric circulation, which is associated to the monsoon and AS cyclone, to the SST anomaly.

Results and discussions
Before we do a comprehensive analysis, it is important to identify large-scale atmospheric circulation patterns which could affect both the monsoon strength and frequency of AS cyclone. For this purpose, we first analyse the association between rainfall over the Indian core monsoon zone and sea surface temperature (SST) during JJAS (Fig. 1a). We identify a prominent negative correlation in the AS, equatorial Pacific Ocean, which also encompasses the ENSO region, and southern Atlantic Ocean SST (SASST). Provided that ENSO exacerbates the strong influence on the surrounding regions, we then isolate the ENSO-independent contribution of SASST on the Indian monsoon (Eq. 1). We notice that negative association between SASST (20° S-1° S, 30° W-10° E) and the Indian monsoon became more evident (Fig. 1b). To further investigate this association, we spatially aggregate SASST to form SASST index, with ENSO removed, and conduct a grid-wise correlation with monsoon rainfall over the Indian mainland. Figure 1c indeed acknowledges the negative relationship over India, with the highest signal noticed in Central India and the core monsoon zone. It is hence evident from this initial analysis that the SASST influences the Indian monsoon. This observation has been reported by previous studies (Kucharski et al. 2007(Kucharski et al. , 2009Pottapinjara et al. 2016;Yadav et al. 2018). SASST mode is similar to ENSO with high SST anomalies in eastern equatorial Atlantic and Gulf of Guinea (GG), however, with a weaker signal and mostly short-lived (Lübbecke et al. 2010(Lübbecke et al. , 2018Burls et al. 2012). Kucharski et al. (2009), Pottapinjara et al. (2014 and Sabeerali et al. (2019a) showed that SASST effects the Indian monsoon through a simple Matsuno-Gill-type response; and this is responsible for the formation of monsoon depressions in the BOB. The warming of SASST causes fewer monsoon depressions, due to reduction in low-level cyclonic vorticity and midtropospheric humidity, which eventually leads to reduced monsoon over India. It is also clear from our analysis that SASST alters the monsoon strength.
With the use of NCEP reanalysis data, we regress the lower (850 hpa) and upper (200 hpa) level wind fields onto SASST from 1951 to 2017 (Fig. 2a, b). The result depicts a strong easterly (westerly) at 850 (200) hpa level, leading to a significant reduction in the wind shear (defined as change in wind speed with respect to 200 and 850 hpa pressure levels) over the AS and BoB. SASST anomaly plot further shows that in recent times (post-1990) they are predominately characterised by a higher number of positive anomaly years (Fig. 2c). From this plot, we extract the information pertaining to the positive and negative phases of SASST. During the positive phase, we do notice easterlies at 850 hpa and westerlies at 200 hpa level, thus leading to the reduction  in the vertical wind shear (Fig. 2d, e) over AS. This situation could provide a conducive environment for the formation of cyclones over AS Wang et al. 2012;Rajeevan et al. 2013). Although, these mentioned studies provide a clear information pertaining to the increased frequency of AS cyclone, the underlying physical mechanism was not clear, which we try to identify in the present study. During the negative phase, however, we notice an increased vertical wind shear with strong westerlies (easterlies) at 850 (200) hpa level (Fig. 2f, g).
The physical mechanism by which the heating anomalies of SASST interconnect to the Indian ocean basin (including Arabian sea) could be explained by a Gill-Matsuno-type quadrupole response (Gill 1980), which is developed by an eddy stream function with kelvin and Rossby waves transporting the signal respectively to the east and to the west (Kucharski et al. 2009). Complete description of the mechanism of SASST on Indian ocean region can be found in Kucharski et al. (2009) and Barimalala et al. (2012). Nonetheless, here, we provide a brief discussion of the same. Warm anomalous SASST signal provides a favourable condition for a rising motion, upperlevel divergence and a minimum velocity potential over equatorial Atlantic and adjoining African region. This condition is balanced by an increased upper-level convergence over Indian ocean and central-western Pacific regions, thereby causing a higher values of velocity potential (Kucharski et al. 2009). However, on the lower levels, opposite anomalies will be formed, indicating a baroclinic response to SASST. Further, in upper-level velocity potential gradient, an upper-level (lower-level) cyclonic (anticyclonic) circulation will be developed, particularly over sub-tropical region, because of Sverdrup balance (Chen 2003;Kucharski et al. 2011). This process further leads to a high-pressure anomaly at the surface over Arabian sea and Indian mainland, triggering anomalous low-level divergence and a weakening of summer monsoon strength over Indian core monsoon zone. Further, at upper levels (low-levels), a warm SASST anomaly leads to westerly (easterly) wind anomalies in Indian Ocean region, including Arabian sea, causing the changes in wind shear characteristics and also does influences a cyclone characteristic over the basin. Additionally, it has also been observed that anomalous warm SASST weakens the time-mean Walker circulation, and it may further contribute to strengthening the teleconnection between SASST and atmospheric circulation over Indian basin (Kucharski et al. 2007(Kucharski et al. , 2008. Figure 3 depicts a schematic representation which encapsulate the principal mechanisms associated with an anomalous warming of SASST. Overall the analysis from observations showed the influence of SASST in changing large-scale atmospheric characteristics over AS. To comprehend the results procured from observations, we further extend our analysis with the use of climate perturbation experiment. Here, we use ICTP AGCM, which could successfully demonstrate the large and persistent Atlantic SST anomalies forcing the low-level convergent flow, bringing moist air during the African and Indian summer monsoon and modulating the summer climate over these regions, consistent with the observed patterns (Kucharski et al. 2007(Kucharski et al. , 2009Feudale and Kucharski 2013;Zhang and Villarini 2019). In particular ICTP AGCM could reasonably reproduce the summer monsoon features over India (Fig. 4). Indian summer monsoon season is characterised by the strengthening of south-westerly winds that brings moisture from ocean to Indian main land (Fig. 4a), which can also be seen in ICTP AGCM simulations (Fig. 4b). We further consider the spatial average of rainfall to compute the climatology for all the months. The observed climatology shows a gradual increase of rainfall from June, after onset, reaching peak during July-August and then decrease towards September, followed by a monsoon withdrawal (Fig. 4c). We find that ICTP AGCM could also capture this observed pattern of climatology over India. Considering the fact that the ability of ICTP AGCM in simulating the monsoon characteristics and also reproducing the climatic response of Atlantic SST warming over the Africa and Indian regions, we therefore select this model in our study. Fig. 3 Schematic representation of the mechanism of SASST teleconnection for an anomalous positive SST and showing its influence on Arabian sea and Indian region, which is adapted from Barimalala et al. (2012). Shading represents the anomalies of surface pressure during the positive phase of SASST, whereas winds at surface level are shown in vectors With respect to ITCP AGCM perturbation study, we conduct two sets of experiments by superimposing positive (Fig. 5a) and negative anomaly (Fig. 5b) of SST to the seasonal climatology for JJAS months in the south Atlantic region. The positive SASST forcing resulted in a decrease in vertical wind shear over AS (Fig. 5c) due to anomalous easterlies (westerlies) at 850 (200) hpa; and this bears a strong agreement with the observations (Fig. 2). The negative SASST forcing, on the other hand, responded opposite to that of positive SASST forcing (Fig. 5d). Thus, these results further confirm the effect of SASST in changing the characteristics of lower (upper) level winds, further leading to changes in vertical wind shear. Analysing the premonsoon cyclones over AS, Zhang and Villarini (2019) found a significant positive correlation with the warming of the south Atlantic Ocean and the frequency of cyclones. Similar to our study, they also identified a reduced vertical wind shear to the warming of SASST. Considering this, and also from previous studies Rajeevan et al. 2013), which reported a favourable condition of reduced vertical shear to the formation of cyclone genesis, we further perform the analysis pertaining to the frequencies of monsoon (JJAS) cyclones over AS during SASST phases. Figure 6a shows the total number of AS cyclones for the period 1979-2010. Many of these cyclones made a landfall over middle-east and Indian subcontinent regions. Additionally, we note that most of these cyclones were formed near the west coast of India. To underpin the role of SASST, we show AS cyclonic tracks during the positive (Fig. 6b) and negative phase (Fig. 6c). It is worth mentioning that frequencies were significantly higher during the positive phase than that of the negative phase of SASST. It should be noted that tropical cyclones over AS usually forms during pre-and post-summer monsoon seasons. This is due to the fact that during summer monsoon season, atmospheric monsoon circulation and associated  (Gray 1968;Evan and Camargo 2011;). In our case, shown in Fig. 6, most of the pressure systems formed during the month of June (transition from pre-monsoon to monsoon season) and September (transition from monsoon to post-monsoon season) of the summer monsoon months. Having said that, the cyclonic storms which form during the month of June, does delay the monsoon onset and also impede the propagation of the monsoon within the subcontinent (Murakami et al. 2017). Thus, it's impact on the billion-strong population residing in India cannot be neglected. Our results based on both observations and climate model perturbation experiments perhaps provide a prominent evidence on the on the role of SASST in influencing the frequency of AS cyclones during JJAS season.  The tropical cyclones (including AS region) and its fragments could be one of the significant drivers of precipitation extremes across different parts of the mid-and high-latitude regions. While, the contribution of cyclones to overall rainfall during JJAS is negligible (Khouakhi et al. 2017) in south Asia, its contribution towards extreme precipitation can hardly be overemphasized (Lau et al. 2008;Knight and Davis 2009;Prat and Nelson 2016). A recent study by Scoccimarro et al. (2020) identified the prominent role of tropical cyclones in significantly reducing the mean rainfall over tropical domain of Maritime Continent. They further add that this reduction in mean rainfall during active cyclones is due to net reduction in moisture transport to the Maritime Continent influenced by the cyclone-associated circulations.
On the Indian continent however, the moisture is primarily transported from the AS, which is approximately around 36% of the total moisture transported by the surrounding oceanic regions (Roxy et al. 2017). As we have already shown here that SASST increases the possibility of cyclones over AS, it is thus important to analyse the moisture transport conditions from the AS during the SASST phases. Figure 7a shows JJAS climatology of moisture transport (in terms of IVT, based on Eq. 3) from 1951 to 2017. Majority of moisture is transported from the AS, western Indian Ocean and BOB-of which AS region appears to be the dominant one agreeing with Roxy et al. (2017) and Pathak et al. (2017). Along with increased IVT, we also notice a strengthening of westerlies at 850 hpa which enables the transport of moisture from source (AS) to sink (Indian continent). Further extension of analysis pertaining to the positive phase of SASST (Fig. 7b) depicts a weakening of lower wind westerlies together with declined eastward water flow into the Indian mainland. During the negative phase of SASST, prohibition of eastward water flow into the Indian mainland is nullified along with strengthening of westerlies at 850 hpa (Fig. 7c). We further determine the moisture contributions from different Oceanic sources to the Indian mainland using a dynamic recycling model [refer Pathak et al. (2014) and Pathak et al. (2017) for more details]. Here, we consider two main oceanic sources viz., western Indian ocean (WIO, also encompasses AS) and central Indian Ocean (CIO), which are the important moisture sources for Indian subcontinent , and also noticed in Fig. 7a. With the use of this model, we compare the difference in mean JJAS moisture contribution originated from the WIO (Fig. 7d) and the CIO (Fig. 7e) during negative phase and positive phase of SASST. A weakening of moisture transport from both of the sources (i.e. WIO and CIO) is observed during positive phase of SASST. The weakening of moisture contribution from WIO and CIO during positive phase is linked to weakening of vertical wind shear and hence weakening of monsoonal current. With the decline of moisture transport during the positive phase, the monsoon rainfall is seen to reduce (Fig. 7f)-the opposite can be noticed during the negative SASST phase (Fig. 7g). The spatially aggregated monsoon rainfall over the core monsoon zone is significantly lower during the positive phase compared to the negative phase of SASST (Fig. 7h). These analyses perhaps provide an evidence that SASST affects the monsoon rainfall over India along with modulating the genesis of cyclones over AS.

Summary and conclusions
Our study is an extension of previous studies (Pottapinjara et al. 2014(Pottapinjara et al. , 2016Sabeerali et al. 2019a), which discussed the relationship between Atlantic zonal mode/SASST and frequency of monsoon depression in BoB and its consequences in reducing the strength of the monsoon. Nonetheless, these studies overlooked the influence of SASST over AS basin, whose moisture contribution to the Indian mainland during summer monsoon season is relatively more than the BoB (Roxy et al. 2017). Furthermore, Pathak et al. (2017) showed that contribution of the WIO (which includes the AS) to Indian land mass during initial phase of monsoon is much larger (35-40%) than the moisture contribution from the BoB. This perhaps indicates that understanding the dynamics of Arabian Sea is more important, which were lacking in the literature. Based on the analysis presented in this study, we underpin the significant role played by SASST on the variability of Indian monsoon and on the frequencies of cyclone over AS; and its associated physical mechanisms. Reduced vertical wind shear due to warming of south Atlantic Ocean can be explained by a classical Matsuno-Gill response, which is consistent both with observations and climate model perturbation experiments. This reduced vertical wind shear aids in the development of cyclone genesis over AS and also prohibits the transport of moisture to the Indian mainland, eventually causing weakening of monsoon circulation-leading to dry condition over India. Our analysis provides sufficient evidence to use the SST anomalies in the south Atlantic Ocean as a potential predictor for both AS cyclones and Indian monsoon strength at seasonal time scales, which is misrepresented by most forecast models (Sabeerali et al. 2019b).
SST over Atlantic Ocean has significantly warmed during recent past (Saunders and Lea 2008). Given a significant role of SASST on AS cyclones and JJAS Indian rainfall, we went one step further and examined the relative role played by the human-induced climate change in driving the increasing trend of SST-with the use of state-of-the-art climate model simulation from CMIP5 (details provided in Table 1). Figure 8a shows the similarity between Hist and HistGhG experiments, especially during the late 20th century-the period during which the effect of anthropogenic warming was apparent (Diffenbaugh et al. 2015). Moreover, decadal variability of Hist simulations are consistent with observed SASST patterns, specifically during post 1960s (inset of Fig. 8a). Further, we perform relative importance analysis [detailed in Johnson and LeBreton (2004), Grömping (2007) and Vittal et al. (2020c)] to quantify the contribution of HistGHG towards changes in SASST. We notice a significant contribution from HistGHG towards demonstrating the observed SST pattern compared to the HistNAT during post 1960 (Fig. 8b). If the GhG emission is not controlled with proper mitigation measures, SASST will increase strongly towards the end of the 21th century as shown by RCP 8.5 and 4.5 scenarios (Fig. 8a).
As the warming of the Atlantic Ocean due to increased GhG emission is apparent ), more frequent cyclones in AS and weakening of monsoon over India can be expected in near future. One way to tackle this issue is by adhering to recent Paris Agreement, which urges all nations to restrict emission such that the warming level does not exceed 2 °C compared to pre-industrial levels and to pursue the limit of 1.5 °C (Vicedo-cabrera et al. 2018). Studies reported a wide range of benefits, such as reduced tropical cyclones (Wen et al. 2019) to increase in mean precipitation over India (Maharana et al. 2020), if the warming limit is set to 1.5 °C based on the Paris agreement. different organizations for making the data available, which includes CMIP5 data from the World Climate Research Programme's Working Group on Coupled Modelling; CRU, HadISST and NCEP reanalysis database. Tropical cyclone track data for the Arabian Sea in the form of eAtlas from India Meteorological Department. The cyclone e-Atlas is produced by India Metorological Department, Cyclone Warning Research Centre (CWRC) and Regional Meteorological Centre (RMC), Chennai. VH acknowledges Dr. Wei Zhang, Univeristy of Iowa, for his contribution in knowledge sharing. Many useful discussions with him helped VH to generate and address the research problem which is presented in this article. VH also acknowledges Prof. Mat Collins, University of Exeter, for his helpful suggestions and comments. Data analysis was conducted at the High-Performance Computing (HPC) Cluster EVE, a joint effort of both the Helmholtz Centre for Environmental Research-UFZ and the German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig. Finally, we would like to thank two anonymous reviewers and the Editor for their constructive comments, which improved the quality of the manuscript.
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/. Fig. 8 a Yearly variation of SASST during JJAS season based on observational data-set and the ensemble (N = 10) of climate model simulations from CMIP5 under different experimental scenarios: natural only historical (HistNat), green-house emission-only (Hist-GhG), all-forcings historical (Hist), and two future RCPs (4.5, and 8.5). Thick solid lines show the multimodel means, and filled areas represent 95% confidence intervals based on the sampling distribution of mean. The inset plot in a represents variations of SASST, but with 10-year moving average. Relative importance analysis of HistGhG and HistNaT simulations on observed SASST for post-1960 (b). The bar plot shows the contribution of each simulations, whereas error bar represents the bootstrapped 95% confidence for 1000 re-samples