A Proposal to Evaluate Drought Characteristics Using Multiple Climate Models for Multiple Timescales

This study presents a method to investigate meteorological drought characteristics using multiple climate models for multiple timescales under two representative concentration pathway (RCP) scenarios, RCP4.5 and RCP8.5, during 2021–2050. The methods of delta change factor, unequal weights, standardized precipitation index, Mann–Kendall and Sen’s slope are proposed and applied with the main purpose of reducing uncertainty in climate projections and detection of the projection trends in meteorological drought. Climate simulations of three regional climate models driven by four global climate models are used to estimate weights for each run on the basic of rank sum. The reliability is then assessed by comparing a weighted ensemble climate output with observations during 1989–2008. Timescales of 1, 3, 6, 9, 12, and 24 months are considered to calculate the standardized precipitation index, taking the Vu Gia-Thu Bon (VG-TB) as a pilot basin. The results show efficient precipitation simulations using unequal weights. In the same timescales, the occurrence of moderately wet events is smaller than that of moderately dry events under the RCP4.5 scenario during 2021–2050. Events classified as “extremely wet”, “extremely dry”, “very wet” and “severely dry” are expected to rarely occur under the RCP8.5 scenario.


Introduction
Drought is a natural hazard related to a deficiency of precipitation for an extended period that results in water shortage for some activities or for some economic sectors [1].The meanings of "drought" depend on different perspectives of stakeholders from farmers to meteorologists [2].Commonly, according to the studies of droughts [3][4][5], the concept of drought is clustered into four types consisting of meteorological, agricultural, hydrological and socioeconomic types.The Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report [6] emphasizes that the world indeed has become more drought-prone during the past 25 years.Drought-affected areas will likely increase in frequency and severity, with implications for sustainable development (e.g., agriculture and forestry production or land degradation).Observed changes in characteristics of droughts (i.e., more intense and longer duration droughts) are widely documented for a variety of regional and ocean basin scales since the 1970s with the emphasis on tropics and substropics [6].In comparison to the Medieval Climate Anomaly (950-1250), more megadroughts appeared in monsoon Asia and wetter conditions became dominant in arid Central Asia and the South American monsoon region during the Little Ice Age (1450-1850) [7].Over a global scale, it is observed that the intensity and/or duration of droughts likely increase in the Mediterranean and West Africa and decrease in central North America and north-west Australia [6].Increased drying is directly linked to higher temperatures and decreased precipitation.It is noteworthy that the palaeoclimate records show that droughts prolonging with a scale of decades or longer have been very likely a repetitive feature of Climate 2018, 6, 79 3 of 16 from multiple models perform better than the single models.Furthermore, the projection errors can be further eliminated by using "unequal weighting".To get this, however, single model skills and relative contributions of the joint model error are required.More importantly, Timothy et al. [20] used a statistical test to define whether an ensembles of multi models with "unequal weighting" is significantly better than without "unequal weighting".The study showed that a value for the relatively small global fraction is illustrated with the method of "unequal weighting".Sanderson et al. [18] suggested a weighting scheme to eliminate some aspects of model codependency in the ensemble.Also, a weighting strategy for an ensemble of CMIP5 was presented by Sanderson et al. [21] in the fourth National Climate Assessment.In general, these studies use a distance metric of models to observations and the distance metric of a pair of models i and j, and a relationship to convert those into a weight.The equal weighting is, remarkably, often used to develop the global ensemble scenarios as a safer and more transparent way to combine models [1], but unequal weights can be better for some areas of the global [20].In most cases of existing ensemble members, the application of any kind of weighting to ensemble variance is mostly discarded, but only considered the weighted mean [1,22,23].This can ignore the intermodel relationships and unexpected values, as extreme weather variables can potentially be more sensitive to changes in the variance [24].In addition, ensemble members typically come from the same model.
In the context of changing climate, studies in drought events at scales of region and basin are valuable for understanding their evolution and impacts on a wide range of fields (e.g., agriculture, socio-economic, environment and natural resources) that occur over certain areas.In this sense, the present study aims to evaluate the wet and drought events in the 21stcentury using multiple climate models for multiple timescales.The precipitation projections from multiple regional climate models driven by multiple global climate models are separately corrected using the method of delta change factor.In addition, an unequal weight method is proposed in the expectation of a better performance of multiple climate projections of precipitation at basin scale in Vietnam as a case study.Weights for each climate simulation are calculated on the basic of rank sum metric of each climate simulation.The rank sum is defined from statistical indices of each climate simulation in comparison with observation.In comparison with the existing methods mentioned above, this approach can measure not only the absolute performance of each model, but the performance compared with the other models in the ensemble with its ranks.The methods of the non-parametric Mann-Kendall (MK) test [25,26] and Sen's slope [27] are then applied to detect the projection trends in meteorological drought for multiple timescales at a significance level of 0.05.The reason for this is that the MK test is widely applied [28][29][30] with advantages of a rank correlation without any request of a particular distribution of data and not affected by the data errors and outliers.

Description of the Case Study Area: Vu Gia-Thu Bon Basin
A plot basin, Vu Gia-Thu Bon (VG-TB), is selected in this study.It is located in central Vietnam, elongating from 16 • 55 through 14 • 55 and from 107 • 15 through 108 • 24 and covers a total of area of approximately 12,577 km 2 .The VG-TB basin is surrounded by two main provincial administrative territories Quang Nam and Da Nang.The basin is characterized by a steep topography and the altitude ranging from 0 m at the coast to 2567 m in elevation above sea level (m.a.s.l) in the west (Figure 1).

Observational Data
The monthly precipitation records are obtained from the Vietnam HydroMeteorological Data Center of the Ministry of National Resources and Environment of Vietnam (MONRE).They are aggregated from the daily series of data.There are two national rain gauge stations (i.e., Danang and Tramy).Other stations including Ainghia, Camle, Giaothuy, Caulau, Hien, Hiepduc, Hoian, Khamduc, Nongson, Queson, Thanhmy and Tienphuoc are popular rain gauge stations which operate manually on the base of volunteers.The location of these stations is displayed on Figure 1.The data is available from 1986-2015.

Gridded Data
The precipitation products are from different assembliese of regional models: (1) The Regional Climate Model version 4 (RegCM4), developed by the International Centre for Theoretical Physics (ICTP).The dynamical structure of RegCM4 firstly developed at the National Center of Atmosphere Research (NCAR) and Penn State University (PSU) for a hydrostatic version of the Meso-scale Model (MM5).A detailed description of RegCM4 can be found in Giorgi et al. [31].The model output of HadGEM2-AO produced by the National Institute of Meteorological Research (NIMR)/Korea Meteorological Administration (KMA) are used as an initial and boundary conditions, referred to REG/HadGEM.Details of HadGEM2-AO are given by Collins et al. [32]; (2) The model of SNU-MM5(Seoul National University Meso-scale Model version 5) [33] is based on a hydrostatic version of the Meso-scale Model and the community land model version 3 (CLM3).The future climatic projections are produced with the HadGEM2-AO, referred to SNU/HadGEM; (3) A regional spectral model, which is also known as Regional Model Program (RMP) of the Global/Regional Integrated Model System (GRIMs) [34]is used in this study.The dynamic frame of RMP is rooted in the National Center for Environmental Prediction (NCEP) RSM.More detailed information about the GRIMs-RMP is provided by Hong et al. [34].This model is driven by the HadGEM2-AO, referred to RSM/HadGEM; (4) The RegCM4 is forced by the model of Max Planck Institute for Meteorology Earth System Model MR (MPI-ESM-MR), which has an ocean horizontal resolution of 0.4° × 0.4° and atmosphere horizontal resolution of 1.9° × 1.9°.It is written under a short symbol of REG/MPI; (5)

Observational Data
The monthly precipitation records are obtained from the Vietnam HydroMeteorological Data Center of the Ministry of National Resources and Environment of Vietnam (MONRE).They are aggregated from the daily series of data.There are two national rain gauge stations (i.e., Danang and Tramy).Other stations including Ainghia, Camle, Giaothuy, Caulau, Hien, Hiepduc, Hoian, Khamduc, Nongson, Queson, Thanhmy and Tienphuoc are popular rain gauge stations which operate manually on the base of volunteers.The location of these stations is displayed on Figure 1.The data is available from 1986-2015.

Gridded Data
The precipitation products are from different assembliese of regional models: (1) The Regional Climate Model version 4 (RegCM4), developed by the International Centre for Theoretical Physics (ICTP).The dynamical structure of RegCM4 firstly developed at the National Center of Atmosphere Research (NCAR) and Penn State University (PSU) for a hydrostatic version of the Meso-scale Model (MM5).A detailed description of RegCM4 can be found in Giorgi et al. [31].The model output of HadGEM2-AO produced by the National Institute of Meteorological Research (NIMR)/Korea Meteorological Administration (KMA) are used as an initial and boundary conditions, referred to REG/HadGEM.Details of HadGEM2-AO are given by Collins et al. [32]; (2) The model of SNU-MM5(Seoul National University Meso-scale Model version 5) [33] is based on a hydrostatic version of the Meso-scale Model and the community land model version 3 (CLM3).The future climatic projections are produced with the HadGEM2-AO, referred to SNU/HadGEM; (3) A regional spectral model, which is also known as Regional Model Program (RMP) of the Global/Regional Integrated Model System (GRIMs) [34] is used in this study.The dynamic frame of RMP is rooted in the National Center for Environmental Prediction (NCEP) RSM.More detailed information about the GRIMs-RMP is provided by Hong et al. [34].This model is driven by the HadGEM2-AO, referred to RSM/HadGEM; (4) The RegCM4 is forced by the model of Max Planck Institute for Meteorology Earth System Model MR (MPI-ESM-MR), which has an ocean horizontal resolution of 0.4 • × 0.4 • and atmosphere horizontal resolution of 1.9 • × 1.9 • .It is written under a short symbol of REG/MPI; (5) The RegCM4 is forced by the model of Institut Pierre Simon Laplace CM5A-LR (IPSL-CM5A-LR), which is the low-resolution version of the IPSL-CM5A Earth system model.It has a horizontal resolution of 1.875 • × 3.75 • with 39 vertical level for the atmosphere and about 2 • (with a meridional increased resolution of 0.5 • near the equator) and with 31 vertical levels for the ocean.In this study, it implies to REG/IPSL; (6) The RegCM4 is forced by the model of Irish Centre for High-End Computing European community Earth-System (ICHEC-EC-EARTH), which is a new Earth system model on the basic of the operational seasonal forecast system of the European Centre for Medium-Range Weather Forecasts (ECMWF).This case is written with a short symbol of the REG/ICHEC.More importantly, all simulations and projections of climatic are run under two IPCC RCP4.5/8.5 scenarios.The RCP4.5 is a stabilization scenario where total radiative forcing is stabilized before 2100 by employment of a range of technologies and strategies for reducing greenhouse gas emissions.Meanwhile, the RCP8.5 is characterized by increasing greenhouse gas emissions over time representative for scenarios in the literature leading to high greenhouse gas concentration.Table 1 lists the name of the models and the number of runs of historical control and RCPs as well as the considered periods.A total of 18 climatic simulations and projections are considered in this study as shown in Table 1.

Delta Change Factor
The raw climate simulations, especially for precipitation time series, are highly biased as mentioned in many studies [6,35,36].Thus, an additional post-processing step (e.g., bias correction of climatic variables) is a standard procedure for related climate change studies.In this study, the delta change method is adopted due to its simple and common use as described in Olsso et al. [37], Lenderink et al. [38], Teutschbein and Seibert [36] and Maraun [35].This approach does not adjust the output of climate models, but uses observations and the change signal of regional climate models forced by global climate models to generate future data.Also, this approach is to avoid considerable variability in the day-to-day change signals and changes in extremes are linearly scaled with changes in the mean.The core of this method is that the historical observations are transformed into future projections using monthly average correction factors that derived from the regional climate models forced by global climate model outputs for the baseline and future climate.It is expressed by the equation: where P GCM.f is the monthly precipitation from the future climate.P GCM.b is the monthly precipitation from the baseline climate.P OBS.f is future projections and P OBS.b is historical observations.With this approach, the correlation structure of downscaled future data in spatio-temporal terms are physically reasonable because it reflects observed conditions.The drawback of this approach, however, is that the future and baseline scenarios just differ in terms of their means and intensity, but other statistics of the data (e.g., skewness or structure of dry and wet days) are mostly unchanged.Also, the sample is limited to the length of the observed record.It should be noted that for a near-term future (2021-2050), the changes in the dry/wet days are probably trivial.Sun et al. [39] used the lag 1 autocorrelation to investigate the stationary land-based gridded annual precipitation .The results showed a stationary annual precipitation over an area of about 80% of the global land surface.Wilks and Wilby [40] show that calculating the autocorrelation of non-zero precipitation amounts is essential at time step of hourly (or sub-hourly) rather than a daily time step.Meanwhile, the autocorrelation between successive nonzero precipitation amounts is usually of little practical importance and quite small.At a time step of monthly precipitation, thus, a stationary assumption in the temporal correlation is made in this study.In other words, rescaling the precipitation time series observed during the baseline could lead to not realistic results when future scenarios that preserve the observed autocorrelation in time are not considered.

Unequal Weights
In this study, a method is suggested to estimate the unequal weights based on the rank sum from each climate simulation run for the analytical hierarchy process.It is called unequal weights.The rank sum is calculated using the statistical indices.The statistical indices are Nash-Sutcliffe efficiency with logarithmic value (ln(Nash)), root-mean-square error (RMSE) and Nash-Sutcliffe efficiency (Nash) for the studied domain.They are basically quantified by measuring the difference between the observed data and the outputs of regional climate models forced by lateral and surface boundary conditions from the European Centre for Medium-Range Weather Forecasts at the monthly scale.The Nash-Sutcliffe efficiency with logarithmic values ln(Nash) is selected because it can be added to expect a better quantification of the performance in different conditions (maximum and minimum values) [41].The other remaining value is widely applied in hydrometeorological fields.The formulations used to compute the goodness-of-fit statistical indicators for each climate simulation run are presented as follows: where n denotes number of months, F i and O i represent simulated and observed monthly data, respectively.The method of unequal weights is briefly expressed with the following steps with an assumption of N climate simulations: (1) Calculation of the statistical indices on the basic of the historical observations and climate simulations from regional climate models forced by the reanalysis data of the European Centre for Medium-Range Weather Forecasts during 1989-2008.Each climate simulation receives a rank from 1 to N depending on the levels of perfect score for each statistical index, starting with the best as 1 and the worst is N.As an example, if the RMSE index of the ith climate model has the best score (the perfect score of RMSE is zero), the received rank is 1.Then, an ensemble rank order (r) as an integer number is calculated from the average of the ranks they span for each climate simulation.(2) Calculation of rank sum for each climate simulation by N-r + 1 with N is the number of climate simulations.(3) Establishing a reciprocal matrix between sets of models a ij = 1/a ji with i,j ranging from 1 corresponding to the best climate simulation which has the largest rank sum to N and a ij = 1 as i = j.a ij is determined by the difference of rank sum between sets of models plus 1. (4) Estimation of weights matrix w ij = a ij /∑ N 1 a ij (i,j = 1, N) (5) Estimation of each weight for each climate simulation w i = ∑ N j=1 w ji (i = 1, N) with ∑ N i=1 w i = 1.

Standardized Precipitation Index (SPI)
In this study, the SPI is constructed for multiple timescales ranging from 1 month to 24 months.The calculation of the SPI is separately applied for each month on the basis of the Gamma distribution.The reason for this is that the Gamma distribution is the most frequently used and fits well with daily precipitation in different studies across Vietnam [42].The alpha and beta parameters of the Gamma probability density function are estimated for each station and for the required multiple timescales.They are used to calculate the cumulative probability distribution of accumulated precipitation.The maximum likelihood solutions [43] and Thom's study [44] are applied to optimally estimate alpha (α) and beta (β).It is especially noteworthy that the reference period adopted to compute the best-fit parameters for the gamma distribution spans 30 years .
A probability transformation is then applied to transform monthly precipitation to a standard normal distribution with a zero mean and standard deviation of one to yield SPI values by preserving probabilities [45].Figure 2 shows the fitness of the SPI data [46].
Climate 2018, 6, x FOR PEER REVIEW 7 of 16 order (r) as an integer number is calculated from the average of the ranks they span for each climate simulation.(2) Calculation of rank sum for each climate simulation by N-r + 1 with N is the number of climate simulations.
(3) Establishing a reciprocal matrix between sets of models aij = 1/aji with i,j ranging from 1 corresponding to the best climate simulation which has the largest rank sum to N and aij = 1 as i = j.aij is determined by the difference of rank sum between sets of models plus 1. (4) Estimation of weights matrix wij = aij/∑ a (i,j = 1, N) (5) Estimation of each weight for each climate simulation wi = ∑ w (i = 1, N) with ∑ w = 1 .

Standardized Precipitation Index (SPI)
In this study, the SPI is constructed for multiple timescales ranging from 1 month to 24 months.The calculation of the SPI is separately applied for each month on the basis of the Gamma distribution.The reason for this is that the Gamma distribution is the most frequently used and fits well with daily precipitation in different studies across Vietnam [42].The alpha and beta parameters of the Gamma probability density function are estimated for each station and for the required multiple timescales.They are used to calculate the cumulative probability distribution of accumulated precipitation.The maximum likelihood solutions [43] and Thom's study [44] are applied to optimally estimate alpha (α) and beta (β).It is especially noteworthy that the reference period adopted to compute the best-fit parameters for the gamma distribution spans 30 years .
A probability transformation is then applied to transform monthly precipitation to a standard normal distribution with a zero mean and standard deviation of one to yield SPI values by preserving probabilities [45].Figure 2 shows the fitness of the SPI data [46].The non-parametric Mann-Kendall (MK) test [25,26] is widely applied to detect the possible trends in many countries such as in China [47], in Serbia [48] in Brazil [49], in Canada [50].In this study, the MK test statistic, S, is applied and briefly represented by: where n is the number of data points; x j and x k are the data values in time series j and k respectively, and In this test, the null hypothesis (H o ) assumes that there is no trend in meteorological droughts over time; the alternative hypothesis (H 1 ) assumes that there is an upward or downward trend over time.The mathematical equations for calculating Var(S) and standardized test statistics Z are presented in previous studies [25,26,47,51,52].An upward, downward, or no trend will be assessed at α significance level of 0.05.The computed probability is greater than the specific significance level α (H o is rejected); the increasing trend responds to a positive value of Z and a negative value of Z indicates a decreasing trend.There is no trend if the computed probability is less than the level of significance (H o is accepted).At the α significance level of 0.05, the null hypothesis of no trend is rejected if |Z MK | > 1.96.

The Sen's Slope Estimator
In order to get the magnitude of a consistent trend (Q), the Sen's non-parametric method [19] is used.It is estimated by the slope of all data pairs (N = n(n − 1)/2) as the following formula: where Q is slope between data points x j and x k , x j and x k are the data values at time j and k (j > k) respectively, j is time after time k.The Sen's is computed by the median slope as: if N is odd where N is the number of calculated slopes.The confidence limits for the nonparameter slope estimator are estimated by the well-known studies [27,29,53,54].

Calculation of Weights for Each Climate Model
As the first step, the delta change factor is applied in this study.Then, the unequal weights are added to define weights for each climate simulations as mentioned in Table 1over the whole VG-TB.Table 2 shows the statistical indices for these climate simulations.For RMSE index, the perfect score is zero with a wide range (0 ≤ RMSE <∞).While the perfect score of Nash and ln(Nash) indices is 1 with a wide span (−∞ ≤ Nash < 1; −∞ ≤ ln(Nash) < 1).Based on these score, ranks (i.e., ensemble rank, rank sum) are calculated and presented in Table 3.It is observed that REG/IPSL gives the best score of ln(Nash), RMSE and Nash.This means that REG/IPSL simulation is well matched to observations.As opposed to this, RSM/HadGEM gives the worst score of ln(Nash), RMSE and Nash.Although REG/IPSL is the best, a multiple climate model is always highlighted in the climate simulations and projections due to reducing uncertainty [7].The reason for this is that the uncertainty comes from lots of factors such as convective cloud parameterization, greenhouse gases scenarios or land surface parameterization.The best rank sum belongs to the model REG/IPSL, followed by REG/ICHEC.In contrary to this, the worst rank sum belongs to the model RSM/HadGEM as shown in Table 3.A reciprocal matrix between sets of models is created based on the rank sum of climate models.As shown in Table 4, it is illustrated that REG/IPSL is more important than REG/ICHEC and much more important than REG/MPI.Meanwhile, REG/ICHEC is more important than REG/MPI and much important than SNU/HadGEM.The least importance within all considered climate simulations is RSM/HadGEM, followed by REG/HadGEM.After establishing a weights matrix, weights are estimated for each climate model based on the 4th and 5th steps in Section 2.3.2.Weights are 0.344 (REG/IPSL), 0.254 (REG/ICHEC), 0.177 (REG/MPI), 0.116 (SNU/HadGEM), 0.07 (REG/HadGEM) and 0.038 (RSM/HadGEM).The results from this method are against the results from the method of equal weights.Both of them are compared with observations using the statistical indices of mean absolute error (MAE) and RMSE.The pair of equal weights and observations gives results of 99.3 for MAE and 164 for RMSE, meanwhile a pair of unequal weights and observations gives results of 95.4 and 155 for MAE and RMSE, respectively.
The perfect score of MAE and RMSE is zero.In other words, the minimum of the RMSE and MAE is obtained when simulated time series perfectly matches observed time series.As mentioned in the delta change factor approach, an assumption of the delta factor is equal to one (i.e., there is no variation) for each climate runs, and the condition of minimum of the RMSE and MAE is obtained.
Climate models with convective precipitation schemes, however, often tend to produce higher precipitation [6,55] compared with the observed data.Moreover, incomplete understanding of natural and its representation is presented within the climate models.Thus, the delta factor assumed by one is unreal, especially in Vietnam where the rain regime is strongly dominated by convective processes.The calculated statistics of MAE and RMSE indicate a better implementation of unequal weights.Climate simulations are closely fit to observations with the unequal weights method as presented in Figure 3.
in the delta change factor approach, an assumption of the delta factor is equal to one (i.e., there is no variation) for each climate runs, and the condition of minimum of the RMSE and MAE is obtained.Climate models with convective precipitation schemes, however, often tend to produce higher precipitation [6,55] compared with the observed data.Moreover, incomplete understanding of natural and its representation is presented within the climate models.Thus, the delta factor assumed by one is unreal, especially in Vietnam where the rain regime is strongly dominated by convective processes.The calculated statistics of MAE and RMSE indicate a better implementation of unequal weights.Climate simulations are closely fit to observations with the unequal weights method as presented in Figure 3.

Calculation of SPI
The weights estimated from Section 3.1 are used to calculate the SPI for 14 stations over the whole VG-TB.It is noted that the weights are applied for both scenarios of RCP4.5 and RCP8.5 during 2021-2050 with an assumption of no change in future.The SPI is designed to project the precipitation deficit for multiple timescales with a weighted ensemble of multiple climate models in future.SPI values are computed for the timescales of 1-(SPI-1), 3-(SPI-3), 6-(SPI-6), 9-(SPI-9), 12-(SPI-12) and 24-(SPI-24) month for 14 stations.For instance, the temporal variation of drought for the timescales of 1 to 24months under RCP4.5 (Figure 4) and RCP8.5 (Figure 5) at Danang station from 2021 to 2050 is displayed.
According to the classification of SPI values [56], generally, there are no events classified as "extremely wet" and "extremely dry" for the Danang station under both scenarios of RCP4.5 and RCP8.5 during 2021-2050 for considered multiple timescales.Under these scenarios, events classified as "near normal" hold more than 90% out of considered months for all timescales of 14 stations.In particular, events classified as "extremely wet" and "extremely dry" are not found out for the all stations under both scenarios of RCP4.5 and RCP8.5 during 2021-2050 for SPI-1, SPI-3, SPI-6, SPI-9, SPI-12 and SPI-24 indices.More importantly, under RCP4.5 scenario, events classified as "very wet" and "severely dry" are lower than 0.6% (mainly in stations located in the east of basin

Calculation of SPI
The weights estimated from Section 3.1 are used to calculate the SPI for 14 stations over the whole VG-TB.It is noted that the weights are applied for both scenarios of RCP4.5 and RCP8.5 during 2021-2050 with an assumption of no change in future.The SPI is designed to project the precipitation deficit for multiple timescales with a weighted ensemble of multiple climate models in future.SPI values are computed for the timescales of 1-(SPI-1), 3-(SPI-3), 6-(SPI-6), 9-(SPI-9), 12-(SPI-12) and 24-(SPI-24) month for 14 stations.For instance, the temporal variation of drought for the timescales of 1 to 24months under RCP4.5 (Figure 4) and RCP8.5 (Figure 5) at Danang station from 2021 to 2050 is displayed.
According to the classification of SPI values [56], generally, there are no events classified as "extremely wet" and "extremely dry" for the Danang station under both scenarios of RCP4.5 and RCP8.5 during 2021-2050 for considered multiple timescales.Under these scenarios, events classified as "near normal" hold more than 90% out of considered months for all timescales of 14 stations.In particular, events classified as "extremely wet" and "extremely dry" are not found out for the all stations under both scenarios of RCP4.5 and RCP8.5 during 2021-2050 for SPI-1, SPI-3, SPI-6, SPI-9, SPI-12 and SPI-24 indices.More importantly, under RCP4.5 scenario, events classified as "very wet" and "severely dry" are lower than 0.6% (mainly in stations located in the east of basin as Caulau) (see Supplementary Material).Meanwhile there is a little change (less than 2.8%) in events classified as "very wet", "moderately wet", "moderately dry" and "severely dry" for SPI-1 for Ainghia, Camle, Giaothuy, Danang and Hien stations under the RCP8.5 scenario.Under this scenario, the remaining stations have no change in drought and wet events for all considered multiple time scales (see Supplementary Material).This indicates only little changes in drought and wet events in the eastern domain of basin under RCP8.5 for considered time scales."Moderately wet"events for the SPI-1 range from 1.4% (Hien, Hoian and Khamduc stations) to 2.8% (Danang station) under the RCP4.5 scenario.Hence, an emphasis of events classified as "moderately wet" and "moderately dry" for multiple timescales over the whole basin under the RCP4.5 scenario during 2021-2050 is documented in Table 5.
Climate 2018, 6, x FOR PEER REVIEW 11 of 16 as Caulau) (see Supplementary Material).Meanwhile there is a little change (less than 2.8%) in events classified as "very wet", "moderately wet", "moderately dry" and "severely dry" for SPI-1 for Ainghia, Camle, Giaothuy, Danang and Hien stations under the RCP8.5 scenario.Under this scenario, the remaining stations have no change in drought and wet events for all considered multiple time scales (see Supplementary Material).This indicates only little changes in drought and wet events in the eastern domain of basin under RCP8.5 for considered time scales."Moderately wet"events for the SPI-1 range from 1.4% (Hien, Hoian and Khamduc stations) to 2.8% (Danang station) under the RCP4.5 scenario.Hence, an emphasis of events classified as "moderately wet" and "moderately dry" for multiple timescales over the whole basin under the RCP4.5 scenario during 2021-2050 is documented in Table 5.As shown in Table 5, the maximum of moderately dry event is 6.5% for SPI-24, followed by 6.3% and 6.0% for SPI-12 and SPI-9, respectively.The median of moderately dry is projected to get the maximum value of 4.2% for SPI-24, followed by 4.0% and 3.7% for SPI-9 and SPI-12, respectively.Importantly, the maximum of moderately wet event is 5.9% for SPI-24, followed by 3.7% and 3.4% for SPI-12 and SPI-9, respectively.It is noteworthy that the values reported from Tables 5 and S1-S6 are assigned by the percentages of drought and wet events under considered multiple timescales.In other words, these values are defined as the percentage of number of months suffering wet and drought events among all months during a 30-year period under considered scenarios.Under the RCP4.5 scenario, in general, timescales are longer, intensity of moderately dry is stronger.In the same timescale, occurrence of moderately wet events is smaller than that of moderately dry under RCP4.5 scenario during 2021-2050.As shown in Table 5, the maximum of moderately dry event is 6.5% for SPI-24, followed by 6.3% and 6.0% for SPI-12 and SPI-9, respectively.The median of moderately dry is projected to get the maximum value of 4.2% for SPI-24, followed by 4.0% and 3.7% for SPI-9 and SPI-12, respectively.Importantly, the maximum of moderately wet event is 5.9% for SPI-24, followed by 3.7% and 3.4% for SPI-12 and SPI-9, respectively.It is noteworthy that the values reported from Table 5 and Tables S1-S6 are assigned by the percentages of drought and wet events under considered multiple timescales.In other words, these values are defined as the percentage of number of months suffering wet and drought events among all months during a 30-year period under considered scenarios.Under the RCP4.5 scenario, in general, timescales are longer, intensity of moderately dry is stronger.In the same timescale, occurrence of moderately wet events is smaller than that of moderately dry under RCP4.5 scenario during 2021-2050.

Projection Trends in SPI
Projection trends in SPI for multiple timescales are considered with an emphasis of drought events under the RCP4.5 scenario.The purpose of this is to detect any trends or changes in drought events.In addition, as clarified in the previous section, just little, even no wet and drought events of extremely wet, extremely dry, very wet and severely dry are found out under the RCP8.5 scenarios.During 2021-2050, out of 14 stations, only five show significant upward trends in November for the SPI-1 index as presented in Figure 6b.There is a similar of upward trends for SPI-3 and SPI-6 in November.Importantly, no trends are detected in November for SPI-9, SPI-12 and SPI-24.Besides that, it is noteworthy that a poor signal of upward trends in SPI-9, SPI-12 and SPI-24 indices is illustrated for months in a year.It indicates an importance of short-term, plans more than long-term plans in related aspects such as water resources planning.
In February, insignificant downward trends are mostly observed over the whole basin (Figure 6a).In the months of March, May, July, October and December, more than half of total considered stations are defined with no trends in SPI-1 index.No trends in SPI-3 index are calculated for almost all the months in a year, except for January, April and November.Downward trends in SPI-1 (in the months of February, April and June), SPI-3 (in April) and SPI-6 (in the months of July and August)

Projection Trends in SPI
Projection trends in SPI for multiple timescales are considered with an emphasis of drought events under the RCP4.5 scenario.The purpose of this is to detect any trends or changes in drought events.In addition, as clarified in the previous section, just little, even no wet and drought events of extremely wet, extremely dry, very wet and severely dry are found out under the RCP8.5 scenarios.During 2021-2050, out of 14 stations, only five show significant upward trends in November for the SPI-1 index as presented in Figure 6b.There is a similar of upward trends for SPI-3 and SPI-6 in November.Importantly, no trends are detected in November for SPI-9, SPI-12 and SPI-24.Besides that, it is noteworthy that a poor signal of upward trends in SPI-9, SPI-12 and SPI-24 indices is illustrated for months in a year.It indicates an importance of short-term, plans more than long-term plans in related aspects such as water resources planning.
In February, insignificant downward trends are mostly observed over the whole basin (Figure 6a).In the months of March, May, July, October and December, more than half of total considered stations are defined with no trends in SPI-1 index.No trends in SPI-3 index are calculated for almost all the months in a year, except for January, April and November.Downward trends in SPI-1 (in the months of February, April and June), SPI-3 (in April) and SPI-6 (in the months of July and August) indices are clearly observed.More importantly, the intensity of drought events in the north and east of the basin is generally stronger than that in other areas.

Conclusions
This study presents a proposed approach to evaluate drought characteristics using multiple climate models for multiple timescales under a context of global warming.The method of unequal weights is proposed to reduce uncertainty of climate projections.The Mann-Kendall and Sen's slope methods are then applied to find the trends in drought characteristics for multiple timescales of 1, 3, 6, 9, 12, and 24 months with six climate projections.The unequal weights are proposed on the basic of rank sum of each climate model.Vu Gia-Thu Bon basin located in central Vietnam is selected to implement this study as a pilot basin.The major findings of the present study can be exposed as an impulse of the precipitation simulations using unequal weights for multiple climate models.Under the scenarios of RCP4.5 and RCP8.5, there are no events classified as "extremely wet" and "extremely dry" for the all stations during 2021-2050 for SPI-1, SPI-3, SPI-6, SPI-9, SPI-12 and SPI-24 indices.A higher magnitude of the drought conditions in the north and east of the basin compared toother areas is found out.Under the RCP4.5 scenario, more importantly, timescales are longer, and the intensity of moderately dry is stronger.Moreover, the occurrence of moderately wet events is smaller than that of moderately dry under the RCP4.5 scenario during 2021-2050 in the same timescales.Under the RCP8.5 scenario, events classified as "extremely wet", "extremely dry", "very wet" and "severely dry" are expected to rarely occur.

Conclusions
This study presents a proposed approach to evaluate drought characteristics using multiple climate models for multiple timescales under a context of global warming.The method of unequal weights is proposed to reduce uncertainty of climate projections.The Mann-Kendall and Sen's slope methods are then applied to find the trends in drought characteristics for multiple timescales of 1, 3, 6, 9, 12, and 24 months with six climate projections.The unequal weights are proposed on the basic of rank sum of each climate model.Vu Gia-Thu Bon basin located in central Vietnam is selected to implement this study as a pilot basin.The major findings of the present study can be exposed as an impulse of the precipitation simulations using unequal weights for multiple climate models.Under the scenarios of RCP4.5 and RCP8.5, there are no events classified as "extremely wet" and "extremely dry" for the all stations during 2021-2050 for SPI-1, SPI-3, SPI-6, SPI-9, SPI-12 and SPI-24 indices.A higher magnitude of the drought conditions in the north and east of the basin compared toother areas is found out.Under the RCP4.5 scenario, more importantly, timescales are longer, and the intensity of moderately dry is stronger.Moreover, the occurrence of moderately wet events is smaller than that of moderately dry under the RCP4.5 scenario during 2021-2050 in the same timescales.Under the RCP8.5 scenario, events classified as "extremely wet", "extremely dry", "very wet" and "severely dry" are expected to rarely occur.

Figure 2 .
Figure 2. The probability transformation from fitted gamma distribution to the standard normal distribution.

Figure 2 .
Figure 2. The probability transformation from fitted gamma distribution to the standard normal distribution.

Figure 3 .
Figure 3. Scatter plot ofthe precipitation simulations with equal and unequal weights minus observations over the whole VG-TB during 1989-2008 (mm/month).

Figure 3 .
Figure 3. Scatter plot ofthe precipitation simulations with equal and unequal weights minus observations over the whole VG-TB during 1989-2008 (mm/month).

Climate 2018, 6 ,
x FOR PEER REVIEW 13 of 16 indices are clearly observed.More importantly, the intensity of drought events in the north and east of the basin is generally stronger than that in other areas.

Figure 6 .
Figure 6.Spatial distribution of Sen's slope in February (a) and in November (b) for SPI-1 index.

Figure 6 .
Figure 6.Spatial distribution of Sen's slope in February (a) and in November (b) for SPI-1 index.

Table 1 .
Climate models and number of runs.

Table 2 .
Statistical indices for multiple climate models.

Table 3 .
Ranks of climate models.

Table 4 .
A reciprocal matrix between sets of models. REG/