Projected future changes in rainfall in Southeast Asia based on CORDEX–SEA multi-model simulations

This paper examines the projected changes in rainfall in Southeast Asia (SEA) in the twenty-first century based on the multi-model simulations of the Southeast Asia Regional Climate Downscaling/Coordinated Regional Climate Downscaling Experiment–Southeast Asia (SEACLID/CORDEX–SEA). A total of 11 General Circulation Models (GCMs) have been downscaled using 7 Regional Climate Models (RCMs) to a resolution of 25 km × 25 km over the SEA domain (89.5° E–146.5° E, 14.8° S–27.0° N) for two different representative concentration pathways (RCP) scenarios, RCP4.5 and RCP8.5. The 1976–2005 period is considered as the historical period for evaluating the changes in seasonal precipitation of December–January–February (DJF) and June–July–August (JJA) over future periods of the early (2011–2040), mid (2041–2070) and late twenty-first century (2071–2099). The ensemble mean shows a good reproduction of the SEA climatological mean spatial precipitation pattern with systematic wet biases, which originated largely from simulations using the RegCM4 model. Increases in mean rainfall (10–20%) are projected throughout the twenty-first century over Indochina and eastern Philippines during DJF while a drying tendency prevails over the Maritime Continent. For JJA, projections of both RCPs indicate reductions in mean rainfall (10–30%) over the Maritime Continent, particularly over the Indonesian region by mid and late twenty-first century. However, examination of individual member responses shows prominent inter-model variations, reflecting uncertainty in the projections.


Introduction
The Southeast Asia (SEA) region, which comprises least developed and developing countries, has been frequently impacted by climate-related disasters such as floods and droughts (Hijioka et al. 2014). With high exposure and low resilience, the SEA region is vulnerable to climate change impacts as the intensity and frequency of extreme events may increase in the future (Ge et al. 2019;Kang et al. 2019;Manomaiphiboon et al. 2013;Ngo-Duc et al. 2014;Raghavan et al. 2017). In particular, a major shift in seasonal rainfall patterns, as well as increases in occurrences of floods and drought would have huge implications on livelihoods, infrastructures, agricultural outputs and food security. For countries to mitigate and adapt to the impacts of future climate change and to increase resilience, policymakers and practitioners must have access to robust information on how future climate change will affect critical sectors in their countries. However, the Working Group II of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report indicated significant knowledge gap on impacts of future climate change on critical sectors in SEA (Hijioka et al. 2014), which could partly be attributed to the unavailability of high-resolution multi-model climate simulations in the region.

3
While future climate projections at global scales have existed from the large collection of simulations using general circulation models (GCMs) under the Coupled Model Intercomparison Projects (CMIPs) (Meehl et al. 2000), the coarse resolution of GCMs (typically at 100-300 km) rendered their applications in climate change impact assessments at local scales less practical (Giorgi et al. 2009;Giorgi and Mearns 1991). Regional climate downscaling, either statistically or dynamically, offers a solution in which information at the local and regional scales can be generated from coarse GCM outputs (Christensen et al. 2007;Fowler et al. 2007;Giorgi et al. 2009;Giorgi and Mearns 1991;Hewitson and Crane 1996;McGregor 1997;McGregor et al. 1993;Wilby and Wigley 1997). While statistical downscaling (SD) does not require large computing resources, its application is dependent on the availability of long time-series of reliable and good quality observed data (Murphy 1999;Tang et al. 2016). Since there is limited access to observed data with sufficiently long temporal coverage and good quality in SEA, the application of SD in this region can be less practical. Furthermore, in the SD approach, the dynamics and climate feedbacks at small scales are not taken into account, and statistical relationships between large and small scales in the future might be different from those in the past and present. On the other hand, dynamically downscaling (DD) can be difficult and time-consuming to implement because of the huge computing resources required to have multiple GCMs, regional climate models (RCMs) and greenhouse gas (GHG) concentration scenarios to generate a range of plausible regional climate change projections (Giorgi et al. 2009;Giorgi and Mearns 1991;Valle et al. 2009).
Modulated by the Asian-Australian monsoon system, the climate of SEA is considered very complex. It is largely characterized by the asymmetric seasonal march of maximum convection from the northern parts of the region during the summer monsoon to the southern parts of SEA during the winter monsoon (Chang et al. 2005;Robertson et al. 2011;Waliser and Gautier 1993). Hence, the maximum rainfall in Indochina occurs during boreal summer while most areas in the Maritime Continent experience maximum rainfall during boreal winter (Chang et al. 2005;Waliser and Gautier 1993). However, with complex terrains where many islands of different sizes are interspersed among regional seas, the regional and local climates of SEA are considered unique in its spatial and temporal variability. In addition, the SEA climate also shows intraseasonal to interannual variations concurrent with the El Niño-Southern Oscillation (ENSO) (e.g. Wang et al. 2003;Tangang and Juneng 2004;Juneng andTangang 2005, Phan et al. 2009;Salimun et al. 2014;Villafuerte and Matsumoto 2015;Tangang et al. 2017;Limsakul et al. 2017;Supari et al. 2018) and Madden-Julian Oscillation (e.g. Tangang et al. 2008;Hidayat and Kizu 2010;van der Linden et al. 2016;Birch et al. 2016;Kim et al. 2017;Xavier et al. 2020). Hence, simulating the past and projecting future climate over the region can be challenging.
In contrast to developed regions such as Europe and North America, the number of regional climate simulations carried out in SEA, which can be reliably used for impact assessments, is limited. Most studies have been limited by the use of a single GCM and RCM, except a few recent ones where multiple GCMs have been downscaled using a single RCM, focusing on some parts of SEA. For example, Kang et al. (2019) used a single RCM to downscale three GCMs over a model domain that covered the Maritime Continent but excluded mainland Southeast Asia. On the other hand, in Rahmat et al. (2014), ECHAM5 and several perturbed Hadley Centre GCMs have been dynamically downscaled over the SEA region. Some countries in SEA have also carried out future climate projection studies, e.g. Vietnam (Katzfey et al. 2014;Ngo-Duc et al. 2014), Philippines (Villafuerete et al. 2019), Malaysia (NAHRIM 2006), and Indonesia (McGregor et al., 2016. On a more encouraging note for SEA was the recent completion of high-resolution multi-model regional climate simulations of the Southeast Asia Regional Climate Downscaling/Coordinated Regional Climate Downscaling Experiment-Southeast Asia (SEACLID/CORDEX-SEA, hereafter referred to as CORDEX-SEA) (Tangang et al. , 2019. Launched in November 2013, the COR-DEX-SEA project was carried out by a consortium of many institutions from within and outside the region, and represents the most comprehensive set of simulations over SEA thus far (https ://www.ukm.edu.my/seacl id-corde x). Based on these simulations, several papers have already been published recently. Tangang et al. (2019) and Trinh-Tuan et al. (2019) used a subset of the CORDEX-SEA simulations for assessment over Thailand and Vietnam, respectively. Herrmann et al. (2020) used one model member of the CORDEX-SEA simulations in investigating the projected changes of sea surface wind in the region. Ge et al. (2019) used three members of the CORDEX-SEA simulations in analysing extreme precipitation over SEA related to 1.5 and 2.0 °C warming. Tangang et al. (2018) also investigated extreme precipitation over SEA under global warming of 2 °C based on a subset of the COR-DEX-SEA simulations. A recent study by Supari et al. (2020) evaluated changes in extreme precipitation over SEA region in the twenty-first century. Using available members of the CORDEX-SEA simulations, the present study aims to analyse the performance of models in simulating December-January-February (DJF) and June-July-August (JJA) mean rainfall in the historical period and evaluate future changes. This paper should serve as a standard reference for CORDEX simulations over the SEA region.
1 3 2 Models, simulation setups and observed data Table 1 provides a list of the simulations conducted in CORDEX-SEA with 7 RCMs and 11 driving CMIP Phase 5 (CMIP5) GCMs and considered as an "ensemble of opportunity" (Daron et al. 2018). A description of the RCMs used is provided in Table 2. Most of the GCMs have been evaluated on their ability to reproduce the monsoonal circulations and climatology over this region in previous studies McGregor et al. 2016;McSweeney et al. 2015;Nguyen and McGregor 2009;Siew et al. 2014). Figure 1 -2040, 2041-2070, and 2071-2099 as early, mid, and late twenty-first century projection periods, respectively. The seven RCMs used in this study utilized various parameterization schemes (Table 2). For RegCM4, WRF and RCA4, brief descriptions of parameterization schemes employed, model configuration and setup were provided in Tangang et al. (2019). The details of the ROM implementation can be found in Sein et al. (2015), while the NHRCM details are given in Cruz and Sasaki (2017). Within the PRE-CIS regional climate modelling system (Jones et al. 2004), the HadRM3P RCM (Massey et al. 2015) is a hydrostatic, 19 hybrid sigma level model, which utilises a rotated north pole coordinate system at 25 km × 25 km horizontal resolution. HadRM3P dynamically downscales HadGEM2-ES by gradually relaxing data saved at six hourly instantaneous frequency for surface pressure as well as three-dimensional horizontal winds, moisture and temperature. HadRM3P is forced with surface boundary conditions (skin temperature and sea ice fraction) from the driving HadGEM2-ES as well as atmospheric boundary conditions for ozone concentration, sulphur dioxide and dimethyl sulphide and volcanic aerosol.
Simulations using CCAM, a stretched-grid GCM, were performed in two stages. First, three quasi-uniform simulations with 100 km resolution are performed from 1971 to 2071 driven by the sea-ice and bias-corrected monthly  Fritsch (1990, 1993) Gregory (1999) Updated version of mass-flux closure McGregor (2003) Fritsch (1990,1993); configured based on Cruz et al. (2016) Tiedtke (1989) with closure for deep convection modified according to Nordeng (1994) Radiation scheme Kiehl et al. (1996) Community Atmos-  (2015) Rotstayn (1997) Lin et al. (1983), Murakami (1990) Lohmann andRoeckner (1996) Land surface scheme Biosphere-Atmos-   is regarded as the "present-day" 25 km simulations for all the runs, as the other GCMs will produce similar climatologies for this period as a result of the biascorrection of their SSTs. The CCAM regional climate simulations used a quasi-uniform C96 conformal-cubic global grid (6 × 96 × 96 grid points). The 25 km simulations used a C120 grid with Schmidt (1977) transformation having a stretching factor of 3.2. Details of the CCAM dynamical formulation and downscaling strategy are provided by McGregor (2005aMcGregor ( , 2015. The cumulus convection scheme is an updated version of the mass-flux closure described by McGregor (2003) and includes entrainment, detrainment and downdrafts. A description of the other model parameterizations is provided by McGregor et al. (2016).
In this study, model performances have been evaluated for the historical period   The statistical significance of biases and projected changes have been calculated using thresholds based on the Monte Carlo permutation test (Baez and Tweed 2013). The 90% significance level of biases and projected changes above random noise are indicated in the maps as forward slashes. Furthermore, the model agreement in the sign of change (henceforth referred as robustness) is also determined, wherein the minimum number of models required to agree was determined using a binomial probability function (Vautard et al. 2014). This model agreement in the sign of changes at 75% level of confidence is indicated in the maps as backslashes. Outputs of driving GCMs are also included in the analysis for comparison with RCM simulations.

Monsoon circulations
The atmospheric circulation over Southeast Asia is largely modulated by the Asian-Australian monsoon, in which the migration of the Intertropical Convergence Zone (ITCZ) and monsoon trough influences the rainfall distribution over the region (Chang et al. 2005;Robertson et al. 2011;Waliser and Gautier 1993). Hence, models must be able to simulate this regional circulation to correctly capture the rainfall distribution. In addition, knowledge of future changes in regional circulation is relevant in explaining changes in precipitation. It is thus necessary to evaluate the performance of the models in simulating the monsoonal circulation, which could be represented by the low-level circulation during DJF and JJA for the winter and summer monsoon circulations over Southeast Asia, respectively (Chang et al. 2005).
The RCM ensemble mean 850 hPa winds for DJF and JJA are shown in Fig. 2. Due to the unavailability of wind data from some simulations at the time of analysis, the ensemble consisted only the first nine members in Table 1, excluding HadGEM2-AO(WRF), HadGEM-ES(PRECIS), and MRI-AGCM(NHRCM). For comparison, the climatological mean of ERA-Interim and ensemble mean of the forcing GCMs are also shown in Fig. 2, as well as the corresponding biases of RCM and GCM relative to ERA-Interim. Despite notable magnitude of biases, both GCM and RCM ensemble means display some consistency, and capture the general pattern of the low-level circulation of ERA-Interim for both seasons. During DJF, the easterly component of the winds in the RCM tends to be stronger over mainland SEA and the Bay of Bengal compared to that of ERA-Interim. RCM wind speeds are also overestimated over Java and an area north of Papua. Similar biases are evident in the GCM. Relatively large biases can be seen during JJA with stronger simulated monsoonal wind along the latitudinal band of 18° N in RCM compared to ERA-Interim. Biases are also notable over the Maritime Continent and western Pacific. Despite the consistency of both RCM and GCM ensemble means to ERA-Interim, inter-model variations within RCMs and GCMs exist (Supplementary Figures Fig. S1 to Fig. S4). The differences in RCMs can largely be traced to the forcing GCMs. However, in some simulations notable modifications can be made by the RCM. For example, in HadGEM2-ES (RegCM4) during DJF, strong southerly winds can be found over eastern Indian Ocean and west of Sumatra in RCM but absent in HadGEM2-ES (Fig. S3 vs. Fig. S4).
In addition to the low-level circulation, it is also useful to evaluate how the models simulate the subsidence field 1 3 as rainfall and subsidence can be tightly linked (e.g. Luis and Pandey 2003;Chang et al. 2005;Hill et al. 2018). However, since the vertical velocity is not available for some models, the divergence at 850 hPa was used as the divergence (convergence) corresponds to subsidence (rising motion) (e.g. Tangang et al. 2019). Figure 3 shows the ERA-Interim and RCM ensemble mean divergence at 850 hPa for DJF and JJA. The divergence fields of GCM have similar patterns with the RCM (not shown). Despite the existence of inter-model variations among the experiments (not shown), the RCM ensemble mean divergence captured the general patterns of those of ERA-Interim well. In terms of consistency, models tend to show high agreement in the sign of biases over the ocean than over land (Fig. 3).
During DJF, both ERA-Interim and the RCM ensemble mean indicate mostly convergence (rising motion) over the Indonesian region centered at 5° S, and divergence (subsidence) over the northern part of the domain centered at 18° N, especially over mainland SEA. The area of convergence over Indonesia corresponds well with the location of the ITCZ during winter monsoon (Waliser and Gautier 1993). Interestingly, despite overall divergent low-level flow in mainland SEA, low-level wind convergence can be seen over the northwest corner of the domain covering the north of Myanmar, Thailand and Laos. Convergence over an area over the southeast of Vietnamese coast can also be found in the ERA-Interim and this is well simulated by the RCM ensemble mean. The RCM ensemble mean also approximates the divergences shown in ERA-Interim during JJA, in which the pattern is roughly opposite to that of DJF. However, some differences at sub-regional/local scales exist. For example, over Peninsular Malaysia, divergence (subsidence) is indicated in ERA-Interim. However, the area is dominated by convergence (rising motion) in the RCM ensemble mean, creating large negative biases in divergence fields. Figure 4 shows the zonally averaged annual cycle of precipitation of GPCC, and the ensemble means of RCMs and GCMs, averaged over the historical period. The "A" shape of GPCC represents the rainfall temporal-spatial distribution over SEA, which is modulated by the Asian-Australian monsoon cycle. During boreal winter, the ITCZ is located around 3-5°S with relatively higher rainfall over these latitudinal bands (Chang et al. 2005;Robertson et al. 2011). Due to the asymmetrical nature of the Asian-Australian monsoon system, the ITCZ is located around 15-20° N during boreal summer (Chang et al. 2005). Hence, the "A" shape represents the ITCZ migration between these two latitudinal bands. Despite inter-model variations (not shown), which is largely reflected in rainfall intensity, all RCMs and GCMs generally simulated the observed "A" pattern. However, the rainfall intensity in RCM ensemble appears higher compared to GCM. For the RCM simulations, RegCM4 and WRF simulations tended to have wet biases while RCA4, PRECIS, ROM and NHRCM produced dry biases, relative to GPCC (not shown).  Figure 5 shows the seasonal climatological rainfall of DJF and JJA of GPCC, the ensemble means of RCMs and GCMs, and their corresponding biases relative to GPCC, averaged over the historical period. Consistent with Fig. 4, the observed spatial distribution of mean rainfall of GPCC shows the north-south gradient, reflecting the Asian-Australian monsoon modulation (Chang et al. 2005;Robertson et al. 2011;Waliser and Gautier 1993). During DJF low (high) rainfall is depicted over Indochina (Maritime Continent). The high rainfall over the Maritime Continent is Fig. 2 The comparison between ERA-Interim 850 hPa winds (first column) and RCM ensemble mean (middle column), and GCM ensemble mean (third column) during historical period. The respec-tive biases are shown in the last two columns. The first (second) row represents DJF (JJA) seasonal mean. The RCM and GCM ensembles mean were averaged from 9 and 8 members, respectively Fig. 3 The seasonal divergence field of ERA-Interim (first column), RCM ensemble mean (second column) and biases (third column). The ensemble mean was based on 9 members. The first (second) row represents DJF (JJA) seasonal mean. Backslashes indicate model agreement in the sign of biases at 75% confidence level 1 3 associated with the monsoonal winds that transport moisture from the South China Sea (Fig. 2) and enhanced convergence (Fig. 3). However, spatial variability can be prominent in some areas in the Maritime Continent, e.g. over Borneo, southern Sumatra and Java (Fig. 5). Over central-western Borneo, high rainfall during DJF is modulated by the existence of the synoptic-scale Borneo vortex (Chen et al. 2013a, b;Tangang et al. 2012). Over Java and southern Sumatra, high rainfall usually occurs in February when the crossequatorial wind strengthens and the ITCZ is located over its southernmost location (e.g. Chang et al. 2005;Tangang et al. 2008). During JJA, rainfall is high over Indochina especially along the coast of Myanmar due to the boreal summer monsoonal winds that transport moisture from the Indian Ocean and Bay of Bengal (Chang et al. 2005). Due to the topographic effect, high rainfall can also be seen in the mountainous region over Laos, and north and central regions of Vietnam.

Seasonal precipitation
Despite the complexity in the spatial distribution of seasonal rainfall, both the RCM and GCM ensemble means appear to well reproduce the main feature, i.e., the north-south gradient in the rainfall pattern and intensity. However, unlike the RCM, the GCM was unable to simulate the fine details of the rainfall distribution. Features of high rainfall over central-western Borneo, Java and southern Sumatra during DJF, and over Laos and northern Vietnam during JJA, were missing in the GCM ensemble mean. In contrast, these features were well simulated by the RCM, implying added values of RCM over GCM. The RCM bias map for DJF over central-western Borneo, Peninsular Malaysia, Sumatra and Java also indicate near zero bias as opposed to the underestimation of − 4 to − 2 mm day −1 in the GCM ensemble mean. During JJA, bias reduction is also featured in the RCM over mountainous areas of Laos, northern and central Vietnam, implying RCM's potential improvement over areas with complex topography (Fig. 5).
For further discussion on inter-model variations and consistency between RCMs and GCMs, biases maps for each RCM and GCM, relative to GPCC, are shown in the Supplementary Figures for DJF (Fig. S5; Fig. S6) and JJA ( Fig. S7; Fig. S8). In Fig. S6 and Fig. S8, the biases of MRI-AGCM was not included due to the unavailability of the data during the time of the analysis. Generally, the RCM biases appear consistent to some driving GCMs especially for RCA4, ROM and PRECIS. Noticeably however, RegCM4 and WRF tended to overestimate the mean rainfall during DJF, e.g., wet biases are featured over Peninsular Malaysia, Borneo, Sumatra and Java in CNRM-CM5 (RegCM4) as opposed to dry biases in CNRM-CM5 GCM. This is also the case for EC-EARTH (RegCM4), MPI-ESM-MR (RegCM4) and HadGEM2-ES (RegCM4). The HadGEM2-AO (WRF) show wet biases over northern Borneo and Vietnam as opposed to dry biases in the HadGEM2-AO GCM. The tendency for overestimation in RegCM4 simulations is also indicated during JJA ( Fig. S7;  Fig. S8). However, not all RegCM4 simulations produced wet biases and opposite to the driving GCMs. The biases in CSIRO-Mk3.6.0 (RegCM4) appear consistent with that of CSIRO-Mk3.6.0 GCM.
A detailed intercomparison of the RCM simulations is shown in the Taylor Diagram (Taylor 2001) in Fig. 6, summarising the performances of all members during the historical period against GPCC as a reference observation. Also shown are four other gridded observational products, i.e. CRU, APHRODITE, CHIRPS and TRMM. Despite differences among gridded observations, as previously highlighted in Juneng et al. (2016), Fig. 6 shows that these observation products are relatively close to each other, such that replacing GPCC with either TRMM, CHIRPS or APHRODITE may produce similar outcomes. However, since TRMM covers only 8 years (1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005) of the historical period, GPCC was chosen to provide a stable climatology of 30 years. The Fig. 4 Zonally-averaged annual cycle of precipitation of GPCC and the ensemble means of RCM and GCM, which were based on 12 and 8 members, respectively 1 3 use of GPCC as reference data for this study was in line with that of Tangang et al. (2019). Figure 6 depicts inter-model variations where most models have root-mean-square difference (RMSD) greater than 1 for both seasons except for CCAM and MPI-ESM-LR (ROM). In terms of spatial correlation, the performances seem dependent on the season where most models registered correlations of 0.5-0.8 for DJF, and 0.4-0.7 for JJA. However, the correlation values of some models in JJA, e.g., EC-EARTH (RegCM4) and CSIRO-Mk3-6-0 (RegCM4), were lower than 0.4. For both seasons, the ensemble mean performed best with correlation values of 0.8 and 0.7 during DJF and JJA, respectively.
In CORDEX-SEA, one of the main interests was to evaluate how well the simulations at 25 km resolution capture the climatology of complex areas within the region. Figure 7 provides such information in terms of annual cycles of rainfall over the 20 sub-regions of interest (see Fig. 1). The rationale for selecting these 20 sub-regions was described in Juneng et al. (2016). In sub-regions covering many islands, e.g., R1, R2, R7, and R10, it would be important to analyse if simulations at 25 km resolution provide reasonable approximation. Figure 7 also shows the annual cycles of GPCC rainfall and the range of values from the RCM simulations. Generally, the RCMs overestimated rainfall in most sub-regions, which was likely caused by the wet biases in the RegCM4 simulations ( Fig. S5; Fig. S7). In most sub-regions, the GPCC means were close to the lower bounds of the RCM simulations, implying the relatively large contribution of the RegCM4 simulations.
Over the sub-regions of R15-R20 of mainland SEA, the ensemble mean performs extremely well especially over R18-R20. For these sub-regions, the annual cycles are characterized by a peak in rainfall during boreal summer and minimal rainfall during boreal winter (Tangang et al. 2019). However, the performance of an individual model can be different from that of the ensemble mean as indicated by the large difference between the lowest and the highest values, especially during JJA, which represent errors of RCMs and GCMs. Hence, this reiterates the requirement of having multi-model simulations and using the ensemble mean rather than an individual model (Valle et al. 2009).
In the Philippines (R1 and R2), relatively large wet biases are shown, particularly during April to October (Fig. 7), possibly due to the wet bias in the RegCM4 simulations. Other factors that may contribute include the inadequacy of the model's 25 km × 25 km resolution to resolve the complex topography and island features, as well as difficulty in adequately simulating the monsoonal circulation and tropical cyclones, which both influence seasonal rainfall (Cruz and Sasaki 2017;Juneng et al. 2016). Over Borneo Island (R3, R4 and R5), the ensemble mean approximates the annual cycles reasonably well although a slight overestimation can be seen especially in R3 during February to April. Large inter-model variation can be seen throughout the year.
The overestimation by the ensemble mean is evident in the sub-regions of Sulawesi (R6), Maluku and West Papua (R7), Papua (R8) and Nusa Tenggara (R10). These subregions have been shown to have multiple climate regimes (Aldrian and Dwi Susanto 2003) and large rainfall variability  (8) members. Hatching indicates differences (or biases) that are significant at 95% above random noise 1 3 (Aldrian et al. 2004). In addition to the RegCM4 tendency for wet simulations, the relatively large biases over these sub-regions could also be due to the models' inability to capture the precipitation processes due to the existence of many islands. Furthermore, a gridded dataset such as GPCC may not adequately represent rainfall in these sub-regions. On the other hand, the ensemble mean approximated the annual cycle much better over the southern part of Papua (R9), Java (R11) and the southern part of Sumatra (R12), with similar performances over northern Sumatra and western Peninsular Malaysia (R14). However, over eastern Peninsular Malaysia (R13), rainfall from all models, including the ensemble mean, tend to peak in October instead of November and December, when the maximum rainfall usually occurs due to cold surges and the Borneo vortex (Tangang et al. 2008(Tangang et al. , 2012Chen et al. 2013a, b;Loh et al. 2016). Hence, such rainfall deficit could be due to the models' inability to correctly simulate the complex circulations and precipitation dynamics over this region. However, the analysis of this issue requires a detailed analysis of the circulation in each model, and is beyond the scope of this paper.

Future projections
3.2.1 Regional circulation Figure 8 shows the projected changes in the RCM ensemble mean 850 hPa seasonal winds for the early, mid and late twenty-first century periods for both RCP4.5 and RCP8.5. For RCP4.5, the ensemble means of early, mid and late twenty-first century periods were averaged from 8, 8 and 7 members, respectively. For the RCP8.5, the corresponding number of ensemble members were 8, 11 and 7. Gradual changes over time are depicted for both RCPs and seasons with prominent changes by the end of the twenty-first century, especially for RCP8.5. During DJF, the directions of anomalous winds over north of the equator are largely opposite to the climatological winds (Fig. 2), implying a slight weakening of the winter monsoonal winds. Similarly, over Java Sea and south Indian Ocean, anomalous southeasterlies prevail implying weakened cross-equatorial and monsoonal winds during future periods, especially in the late twentyfirst century. On the other hand, anomalous westerlies are depicted over the eastern part of the domain. During JJA, stronger monsoonal winds can be seen over Indochina and northern South China Sea whereas winds are weaker over Java Sea. Overall, the projected low level winds indicated low level divergent flow over the Maritime Continent, largely consistent with the low level divergence shown in Fig. 9. Generally, enhanced divergence (convergence) over the Maritime Continent (Indochina) is projected in the future, especially during JJA (DJF).

Seasonal precipitation
The projected changes in seasonal mean precipitation from the RCM and GCM ensemble means, relative to the historical period 1976-2005 for the early, mid and late twenty-first century under RCP4.5 and RCP8.5 scenarios, are shown in Figs. 10 and 11, respectively. For RCP4.5, the RCM ensemble means of early, mid and late twentyfirst century periods were averaged from 9 members while for the RCP8.5, the number of ensemble members were 11, 13 and 11 for early, mid and late twenty-first century, respectively. For all periods and RCPs, the GCM ensemble means were based on 8 members. The projected changes in GCM and RCM simulated mean rainfall are generally similar, although notable differences can be seen. Significant and robust increases of rainfall (10-20%) are depicted during DJF over northern Thailand, Laos and northern Vietnam. The projected changes over Thailand is consistent with Tangang et al. (2019). Over northern Myanmar, while minimal changes are projected during early twentyfirst century for both RCPs and models, significant and robust increases (10-20%) can be found by mid and late twenty-first century. However, during DJF over Cambodia and southern Vietnam, no significant and robust changes are projected except over Cambodia during mid twentyfirst century for RCP8.5. The projected changes during mid twenty-first century over Vietnam, Cambodia, Laos and eastern Thailand are generally consistent with that of the national report of Vietnam on future projected climate (Katzfey et al. 2014). In the Philippines, the projected changes show a slight wetting (drying) tendency in DJF (JJA) (Fig. 10), which is fairly consistent with the projected changes in Villafuerete et al. (2019) for mid twenty-first century under RCP8.5. In the GCM however, a slightly wetter condition is projected over the Philippines in JJA (Fig. 11).
While consistency between GCM and RCM is high during DJF for RCP4.5 in Indochina, some differences are depicted over the Maritime Continent, especially during mid and late twenty-first century. Significant and robust increases are projected over Sumatra, Kalimantan and Papua by the GCM in DJF but largely no changes in the RCM projection. This is also the case for RCP8.5 during late twenty-first century. Such inconsistencies can be caused by the modifications of the sign of change in the RCM compared to the driving GCMs. As shown in the Supplementary Figures of Fig. S9 and Fig. S10, drier condition is projected in some RCMs in contrast to the wetter condition in GCMs. For example, during DJF of late twenty-first century under RCP8.5, CNRM-CM5 (RegCM4) and MPI-ESM-MR (RegCM4) projected dry condition over Sumatra and Kalimantan (Fig. S9) whereas the CNRM-CM5 and MPI-ESM-MR GCMs indicated the opposite future change (Fig. S10). Such inconsistency between RCM and GCM is not uncommon and has also been reported in other studies (e.g., Dosio and Panitz 2016;Fernández et al. 2019).
There are some notable differences between projected rainfall in RCM and GCM for JJA under both RCPs (Figs. 10, 11). While both RCM and GCM show agreement on significant and robust reduction in rainfall (10-30%) over Indonesian region, disagreement can be found in other areas. These include Indochina, northern Borneo, Peninsular Malaysia, Papua and the Philippines where mean rainfall increases are projected in the GCM but generally drier condition is indicated in RCM. Such inconsistencies can be attributed to the switching of the sign of changes in RCM 1 3 compared to GCM (Fig. S11 and Fig. S12). Furthermore, it is not only RegCM4 that tended to switch the sign of change but RCA4 as well. However, over Myanmar for both RCPs and over northern Borneo, northwestern Thailand and Laos for RCP8.5, both GCM and RCM projections appeared consistent. While the projection of reduced rainfall over eastern Vietnam during mid twenty-first century for RCP8.5 is consistent with Katzfey et al. (2014), other areas indicated opposite changes.

Discussion
Due to the complex topography, archipelagic features and oceanic influences, regional climate simulation over Southeast Asia can be considered a very challenging task (e.g., Juneng et al. 2016;Ngo-Duc et al. 2017;Ratna et al. 2017;Chung et al. 2018;Kang et al. 2019). Error and uncertainty in the simulations are difficult to reduce and can be attributed to the shortcomings in the RCMs (e.g., physical parameterizations, resolutions, etc.) and GCMs (as boundary forcing). A multi-model approach is used in the Fig. 8 The RCM projected changes in low-level wind for early (left column), mid (middle column) and late twenty-first century (right column). The first and second rows (third and fourth rows) indicate the projected changes in DJF (JJA) for RCP4.5 and RCP8.5 1 3 CORDEX-SEA simulations at 25 km horizontal resolution using 7 RCMs and 11 GCMs, which allowed the uncertainty associated with the RCMs and GCMs to be sampled. The circulations (both low level winds and divergence field) and rainfall climatology were reasonably approximated by the ensemble means of RCMs and GCMs. Furthermore, the ensemble mean of RCMs showed added values by simulating some finer details in the rainfall spatial patterns. Nevertheless, consistent with Juneng et al. (2016), RegCM4 persistently produced wet simulations, contributing to the overall wet bias in the ensemble mean. In future work, RegCM4 may require further tuning and sensitivity analysis to reduce the wet biases. While the wet biases in the RegCM4 simulations contributed to the large inter-model variations when evaluated at sub-regions, a higher resolution than 25 km may be required in future simulations to resolve the complex topography of many islands. Some synoptic regional circulations appear to be captured by the RCMs, which was missing in the GCMs, e.g., the Borneo vortex that brings high rainfall over Borneo (Chen et al. 2013a). However, the simulated rainfall over east of Peninsular Malaysia peak in October, not during November and December. This appears to indicate that the rainfall delivery mechanism, i.e., through cold surges and easterly waves (Tangang et al. 2012;Chen et al. 2013a), may not be well simulated in both GCM and RCM. The projected changes in low-level winds and divergence fields indicated enhanced subsidence (rising motion) over the Maritime Continent (Indochina) especially during boreal summer (winter). Projected increased rainfall of 10-20% over Indochina during DJF throughout the rest of the twentyfirst century, which is consistent in both RCM and GCM, would bring positive effects to various sectors including water resources, agriculture, power-generation etc., given the very low rainfall over this region during this season. At the same time, a reduction of 10-20% in mean rainfall during JJA over Cambodia, Vietnam and eastern Thailand could have serious implications especially for future agriculture and food production, which may have worldwide implications as these areas are currently considered as food-basket regions (Keskinen et al. 2010).
The most striking changes are the significant and robust mean rainfall reductions of 10-30% over Indonesia during JJA, especially over Sumatra and Kalimantan by the middle and end of the twenty-first century under RCP4.5 and RCP8.5. This tendency for enhanced drying is consistent in both RCM and GCM simulations, which has also been highlighted in other studies (e.g. Kang et al. 2019;Giorgi et al. 2019;IPCC 2013), and can be associated with enhanced divergence and subsidence over the Maritime Continent (Fig. 9). While the causes may require further analysis, previous studies seem to indicate the enhanced subsidence can be associated with the deep tropical squeeze resulting from the equatorward contraction of the rising branch of the Hadley Circulation as the climate continues to warm (Fu 2015).  Fig. 9 except for precipitation, expressed as a percentage anomaly (%) relative to the mean in the historical period. Forward slashes indicate the changes that are significant at 90% level above random noise while backslashes show model agreement in the sign of change at 75% level of confidence 1 3 A recent study by Supari et al. (2020) indicated that the consecutive dry days (CDD) over Indonesia, especially Sumatra and Kalimantan, are projected to increase by more than 50% for JJA and September-October-November (SON). Such a drying condition is experienced only inter-annually during El Niño events in the current climate (Juneng and Tangang 2005;Lestari et al. 2016;Supari et al. 2018), which usually provides favorable conditions for widespread peatland and forest fires, causing haze episodes that have significant economic impacts in Indonesia and surrounding countries (Tangang et al. 2010). With the likelihood of yearly occurrence of drought in future warmer periods, instead of inter-annually under El Niño conditions (Supari et al. 2018), the impacts of forest fires can be expected to worsen. Furthermore, the increasing likelihood for more frequent extreme El Niño events in future warmer periods (Cai et al. 2014) will exacerbate future challenges in mitigating drought and forest fires in Indonesia.
Information on uncertainty and ranges of projected changes in mean rainfall can be useful especially for users at the indicated 20 sub-regions. This information is provided in the Supplementary Figures of Fig. S13 and Fig. S14, showing plots of the projected surface temperature changes (in °C) versus the projected rainfall changes (in %) for all models including the ensemble mean. Only 11 ensemble members were available for the end of the twenty-first century projection (see Table 1). Although the ensemble mean shows no significant future changes in many sub-regions, the range of projections can be large. This represents the uncertainty in the projections, which is attributed to inter-model variations in both RCM and GCM. The users are recommended to consider these uncertainties when using CORDEX-SEA simulation products. The upper and lower limits of the projected changes in mean rainfall should be considered when estimating possible impacts of future climate change in those subregions (e.g. Daron et al. 2018). However, the use of individual model outputs for driving impact models may require bias-adjustment procedures (Ngai et al. 2017). In addition, large biases in some sub-regions within the archipelagic Maritime Continent (Fig. 7), tendency for wet biases for RegCM4 simulations (Fig. S5 to Fig. S8), and inconsistency on the sign of change with the forcing GCMs (Fig. S9 to Fig. S12) justify the need to fine tune the simulations in future undertaking. Furthermore, the use of recent CMIP Phase 6 (CMIP6) GCM simulations (Eyring et al. 2016) and higher spatial resolutions (< 25 km) in climate simulations may reduce the biases and uncertainty.

Conclusion
With 7 RCMs and 11 forcing GCMs, the regional climate simulations of CORDEX-SEA represent the most comprehensive set of high-resolution regional climate simulations carried out over the Southeast Asia (SEA) region thus far. This paper examines the performance of these models in the historical period  and assesses future climate projection over SEA for early (2011-2040), mid (2041-2070) and late twenty-first century (2071-2099) periods under both RCP4.5 and RCP8.5 scenarios. With reasonably large number of ensemble members, which come from various RCMs, GCMs and constrained by different parameterization schemes, the set of simulations of CORDEX-SEA should provide robust estimates of future climate in this region, not just in terms of mean climate but also associated uncertainties. The following are key findings of the analysis: (i) Seasonal low-level circulations of ERA-Interim during the historical period is captured reasonably well by the RCM ensemble mean. This is best reflected in both the 850 hPa wind and divergence fields, where the monsoonal features of both fields have been simulated well by the models. However, biases exist, reflecting the inter-model variations among RCMs and GCM forcings. (ii) The model ensemble mean captures well the spatial -temporal mean rainfall distribution over Southeast Asia, reflecting the models' ability to simulate monsoonal rainfall migration in the region associated with the Asian-Australian monsoon system. Compared with GCMs, RCMs showed added value, but significant and mostly wet biases reflect inter-model variations and overestimation in some models, primarily contributed by RegCM4. The ensemble mean simulated reasonably well the annual precipitation cycle of GPCC over the 20 sub-regions, especially over Indochina, except over some sub-regions comprising many islands in the Maritime Continent. (iii) Projected low-level circulations show low-level divergent flow over the Maritime Continent, causing enhanced subsidence in this region especially during JJA, and rising motion over Indochina. (iv) Enhanced mean rainfall of 10-20% are projected over Indochina during DJF for both RCPs throughout the twenty-first century. During JJA, enhanced drying condition and mean rainfall reductions of 10-30% are projected over the Indonesian region.
Overall, in addition to previous but limited studies (Kang et al. 2019;Rahmat et al. 2014;Katzfey et al. 2014;Villafuerete et al. 2019 and other initiatives at national levels), this study advances the scientific understanding of future climate change over the SEA region. Future studies should focus on reducing the biases in models, and to analyse whether individual models capture the complex and multiscale circulations over the region. In order to facilitate future studies and data dissemination, the model outputs of COR-DEX-SEA are archived in a dedicated data portal called the Southeast Asia Regional Climate Change Information System (SARCCIS), hosted by the Ramkhamhaeng University Center for Regional Climate Change and Renewable Energy (RU-CORE) in Bangkok, Thailand (https ://www.rucor e.ru. ac.th/SARCC IS). This data archive is linked to the Earth System Grid Federation (ESGF).

3
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/. 1 3