Changes in Climate Extremes over North Thailand, 1960–2099

1Global Institute for Water Security, School of Environment and Sustainability, University of Saskatchewan, 11 Innovation Boulevard, Saskatoon, SK, Canada S7N 3H5 2Agricultural Systems and Engineering, School of Environment, Resources and Development, Asian Institute of Technology, P.O. Box 4, Klong Luang, Pathumthani 12120, Thailand 3School of Engineering and Technology, Asian Institute of Technology, P.O. Box 4, Klong Luang, Pathumthani 12120, Thailand


Introduction
Analyses of climate change, climate variability, and their extremes are becoming more important as they clearly affect the human society and are essential for exploration of ecological and societal changes.One of the motivations for such analyses is to establish if and how global warming [1] might have influenced the climate on a global as well as on regional scale.In the environment and on the economic activities, the impact of climate change is reflected by the change of frequency, duration, and intensity of climate extreme events.After the third assessment report of the Intergovernmental Panel on Climate Change (IPCC), projections and detection of climate change are no longer relegated to global averages and have expanded to include extremes [2].Various climate extreme events such as more recurrent and more intense droughts, floods, cyclones, and extreme rainfall events pose increasing and, oftentimes, limitless threats to environment, water resources, and agricultural production [3,4].
General Circulation Models (GCMs) provide a reasonable basis for estimating future climates; however, outputs from GCM cannot be applied directly in many environmental and water resources protection studies because of some limitations, which include poor representation and simulation of local-scale processes and coarse spatial resolution.Therefore, these issues can be addressed by downscaling approaches.First approach of downscaling (statistical downscaling) climate data is through the use of statistical regressions that establish the relationship between large-scale variables, such as the driving factors derived from GCMs, to local level climate conditions.In the second approach (dynamical downscaling) Regional Climate Models (RCMs) are embedded into GCMs.The RCM is essentially a numerical model in which GCMs are used to fix boundary conditions.So, applying a downscaling technique is recommended for regional-scale analysis in IPCC-AR4 [5].More recently, downscaling is widely applied in climatology for scenario construction and simulation/prediction of mean, minimum, and maximum air temperature and rainfall [6][7][8].
A number of investigations have attempted to identify observed and projected future changes in climate extremes on a global scale (e.g., [9][10][11][12]).A study on the variability in occurrence of extreme climate events at a global scale during the second half of the twentieth century [9] showed general increase in warm nights, decrease in the number of frost days, more frequent heavy rainfall events, less frequent cold temperature extremes, and decrease in the intra-annual extreme temperature range.Alexander et al. [10] showed that over 70% of global land area observed significant increase in the annual occurrence of warm nights and significant decrease in occurrence of cold nights.Projected changes in precipitation and dryness extremes are more ambiguous than that in temperature extremes at the end of 21st century [11].Another study by Donat et al. [12], a research to improve global gridded dataset of extreme indices (HadEX to HadEX2, [10]), showed widespread significant changes in temperature extremes with warming especially for those indices derived from minimum temperature.Precipitation also showed similar changes but the changes were much more spatially heterogeneous compared to temperature changes.
A few studies have examined trends in climate extremes, other than changes in mean values across the Southeast Asia.Hence, little is known about trends in extremes of temperature and rainfall in this region.Manton et al. [13] analysed climate extreme indices for historical period in the Southeast Asia and South Pacific.Regional studies across the Asia-Pacific region have shown significant increases in occurrences of annual number of hot days and warm nights and decreases in occurrences of annual number of cool days and cold nights over the past few decades (e.g., [13,14]).Results from Manton et al. [13] revealed significant decrease in the number of rainy days, significant increase in the proportion of total rainfall from extreme rainfall, significant increase in the minimum and maximum temperature, increase in the number of warm nights, and decrease in the number of cold nights and cool days at Nan Province of North Thailand.Chinvanno et al. [15] developed high-resolution climate scenarios using PRECIS (Providing Regional Climates for Impact Studies) RCM and showed that, during the observed period, most part of the North and Northeastern Thailand had 3-4 hot months over the year, while the Central Plain and the Southern Region had slightly longer summertime.Preliminary finding on MM5-RCM simulations of future projection by Kreasuwan [16] during 2010-2029 revealed that the average temperature in Thailand will increase slightly, and this trend will continue further in the future (during 2040-2059).In a study of Southeast Asia START (System for Analysis, Research and Training) Regional Centre (SEA START RC) [17], it was found that the night time temperature will get warmer than the daytime temperature by the 2045-2065 period.For the Northern mountainous and valley zone, annual average maximum and minimum temperature are expected to rise from 32.41 to 35.82 ∘ C and 20.43 to 24.82 ∘ C, respectively, when compared to the median value among the results from 8 GCMs.
North Thailand is a critical area for a number of regional sustainable developmental issues.The mountainous landscape has relatively high forest cover.It contains the upper reaches of most of the watersheds feeding into the Chao Phraya River system, including the largest one, the Ping River Basin (PRB, catchment area of about 35,000 km 2 ).This basin is an important wet-season source of water stored in the downstream Bhumibol Reservoir (located in Tak Province of North Thailand), which serves irrigation and commercial users during a lengthy dry season when drought may be problematic [18].The basin area was previously covered by subtropical forests that have been slowly converted to agricultural lands in the recent past due to a host of economic, social, and political drivers that have commonly been driving land-cover/land-use changes through montane mainland Southeast Asia [19].Greater Chiang Mai city area experienced a farm-to-city transition where built-up areas increased from 9% to over 33% between 1989 and 2009 [20][21][22].The effects of urbanization that include the construction of major ring roads and changes in the Ping River bank and floodplains are blamed for causing more serious floods in recent years [21,[23][24][25].Flooding and mudslides affected a large area in North Thailand, including Chiang Mai, Chiang Rai, Phayao, and Mae Hong Song Provinces [26].The 2011 Chao Phraya River flood that caused $45 billion damage in Thailand highlights the difficulty of water management [27,28].Despite today's popularization of irrigation, rainfed crop fields remain in a large part of North Thailand.Since this type of cultivation largely depends on weather conditions, it has potential vulnerability against future climate changes accompanied by global warming.Changes in climate extremes have relations and significant impact on all of the anthropogenic activities and natural hazards.These issues constituted the motivation for this research.The purpose of the study is not to find out the best climate model rather than identifying the past and future changes of climate extremes in the regional scale through different emission scenarios.Therefore, in this study, observed trends in temperature and rainfall extremes in North Thailand were investigated over the last 51 years.Also, to detect the possible changes of extremes up to the year 2100, assessments of regional climate projections were made, derived from PRECIS regional model and HadCM3 (Hadley Centre Coupled Model version 3) global model.
The paper is organized as follows.Section 2 presents an overview of the study area and data used in this research, followed by methodology in Section 3. In Section 4, results of various analyses are presented and finally conclusions of the study are given in Section 5.

Geography. Administratively as well as geographically
Thailand has four regions, namely, the Northeastern, North, Central Plain, and Southern.This study focused on four provinces of North Thailand, namely, Chiang Mai, Chiang Rai, Lamphun, and Tak (Figure 1).These provinces cover an important drainage area of the PRB and have a significant impact on country's agricultural activities.The Ping River, a major tributary of the Chao Praya River, plays a vital role in  the country's economy.The Bumibol Dam on the Ping River adds additional importance to this region.Geographically, this region is mountainous and hilly with narrow valleys and occasional intermountain depressions forming a highly complex relief that merges across the international borders into Myanmar, Laos, and Yunnan (China), situated between 14.56 ∘ -20.27 ∘ North and 97.20 ∘ -101.47 ∘ East.These northsouth oriented hill ridges are parallel from west to east and intersected by some major valeeys, particularly those near Chiang Mai and Chiang Rai provinces.It experiences a strong monsoon climate with heavy summer rains between mid-May or June and October but experiences dry weather in winter because of the northeast monsoon which is a main factor that controls the climate of this region.The northern part of the country usually experiences an extended period of warm weather because of its inland nature and tropical latitude zone [29].In this region, the average annual temperature ranges from 20 to 34 ∘ C and the average annual rainfall varies between 1000 and 1800 mm depending on the location.

Data Used.
Observed daily minimum temperature ( ∘ C), maximum temperature ( ∘ C), and total daily rainfall (mm) data were collected from Thailand Meteorological Department for the period of 1960 to 2010 for Chiang Rai, Chiang Mai, and Tak provinces (Table 1).In Lamphun Province recording of climate data started from 1981.Two sets of model data, each of which has two SRES scenarios (A2 and B2), were used in this study to the end of 21st century .The first set of data was collected from SEA START RC, where the climate data were downscaled using PRECIS RCM, which used ECHAM4 GCM data as an initial condition.Another set of climate data was produced by using a Statistical Downscaling tool (Statistical Downscaling Model (SDSM)) for the coarse resolution climate variables obtained from HadCM3 GCM.

Methodology
3.1.Downscaling Climate Data.The SDSM [30] calculates statistical relationships, based on multiple linear regression techniques, between large-scale (the predictors) and local (the predictand) climates.There are three implicit assumptions involved in statistical downscaling [31].Firstly, the predictors are variables of relevance and are realistically modeled by the host GCM.Secondly, the empirical relationship is also valid under altered climatic conditions.Thirdly, the employed predictors fully represent the climate change signal.These relationships are developed using observed weather data, assuming that these relationships will remain valid in the future, that can be used to obtain downscaled local information for some future time period by driving the relationships with GCM-derived predictors.The observed large-scale predictors have been derived from the National Centers for Environmental Prediction (NCEP) reanalysis dataset [32].This model becomes more accepted in recent years due to its applicability in a wide region and simplicity of establishing a relationship between predictor and predictand variables for future time zone [33].
Substantial differences may exist between projections from climate models; use of a single model does not provide results that adequately reflect the uncertainty inherent in regional climate scenario generation [34].This study incorporates output of two models for a reliable projection of future climate.The SDSM is used to downscale the scenarios (A2 and B2) of HadCM3 by using the predictor variables developed by NCEP.HadCM3 was chosen for this analysis because it is one of four sets of climate model (HadCM2, HadCM3, CGCM2, and CSIRO) output provided with SDSM.The SDSM has successfully been used in several other downscaling studies [35][36][37].Before forecasting the future climate variables, SDSM has been validated for the 1980-1989 period.The SEA START RC used ECHAM4 GCM (resolution 280 × 280 km 2 ; Max Plank Institute, Germany) for climate scenario generation and PRECIS RCM (developed at the Hadley Centre at the Meteorological Office, spatial resolution of 25 × 25 km 2 and emission scenarios B1, B2, A2, and A1FI) for downscaling the GCM data for specific region.The main reason to select PRECIS RCM in this research is that SEA START RC specially handles Southeast Asian region to predict more precise climate and freely provide data to researchers.Hence, A2 and B2 scenarios of PRECIS RCM data were collected from Southeast Asian START Regional Centre website (http://www.start.or.th/),where data are readily available to download.
A2 and B2 scenarios are the storylines which describe the development of different social, economic, technological, environmental, and policy dimensions.A2 is the future of the differentiated world, with emphasis on human wealth.Major underlying themes include strengthening regional cultural identities, with an emphasis on family values and local traditions, high population growth, and less concern for rapid economic development.B2 scenario emphasizes on sustainability and equity (Mixed Green Bag) and local solutions to economic, social, and environmental sustainability in the world [38].3.2.Extreme Indices.Extreme climatic events, through the prolonged hot and cold spells, have numerous impacts on a nation which include significant loss of life, agriculture, production, water resources, forestry, energy, and infrastructure.
It is essential to observe the climate extremes and continue to look for identification of changes in extremes as our resources and infrastructure are becoming more vulnerable to ruthless and extreme weather [39].All of the extreme indices (13 for temperature and 11 for rainfall, Table 2) are calculated from the daily observation of temperature and rainfall by using RClimDex software [40].These indices are recommended by the Climate Variability and Predictability (CLIVAR) Expert Team on Climate Change Detection, Monitoring, and Indices (ETCCDMI), adopted by the IPCC AR4.Data quality control is a necessary step before the analysis of climate variation, because erroneous outliers can impact seriously on trends [41].For this, each data point was examined based on neighboring data, as well as on the understanding of the climate of the region.The data were edited if the problem is obvious (e.g., negative precipitation, 195 degrees changed to 19.5, and maximum temperature less than minimum temperature), set to missing if it is clearly a problem with unknown solution, or kept if deemed probably valid.

Trend Analysis.
The trend of climate extremes presented in this paper is calculated by using the nonparametric Mann-Kendall (MK) test [39,40] that has been recommended by the World Meteorological Organization in assessing trends in environmental time series data [42].The test has been widely used in trend detection analysis for hydrological and meteorological data [12,[43][44][45][46][47].It is nearly as powerful as its parametric competitors.This method was chosen because it does not consider any distribution of the variable and some of the climate indices do not consider Gaussian distribution [12].To estimate the true slope of an existing trend (as change per year), Sen's nonparametric method [48] is used.Sen's slope method gives a robust estimation of the trend and is used in various trend analyses of climate data [10,12,42].Different significant levels (10%, 5%, 1%, and 0.1%) were used to consider a trend to be statistically significant.

Results and Discussion
In this section of the paper, results of various analyses ranging from trend analysis for observed period are presented first, followed by seasonal variation of climate parameters, model performance evaluation, and, finally, projected changes of climate parameters for the future period.Where necessary, link to other related studies is established in different subsections.For brevity, analyses of climate extremes results are presented for Chiang Rai Province and selected results for other provinces.Period (1960Period ( -2010)).Trends and statistical results of temperature extremes are presented in Table 3, where increasing trend of temperature is experienced in the observed period.Both maximum and minimum temperature indices (TXx, TXn, TNx, and TNn) increased with significant trend, with the fact that minimum temperature was increased more than maximum temperature.The increasing trend of summer days (SU35) and tropical nights (TR25) favored the increment of temperature in the study area except in Chiang Rai Province.Warm days and warm nights (TX90p and TN90p) also increased with very high significant trend whereas cool days (TX10p) and cool nights (TN10p) decreased [13,49] with the same trend along with mild to steep slope.Similar results were found by another study of Sharma and Babel [50] for Western Thailand.In the observed period, significant increases in warm extremes (TX90p, TN90p, SU35, TR25, and WSDI) were experienced in North Thailand where reverse trend was observed in case of cold extremes (TX10p, TN10p, and CSDI) which is also favored by the negative trend of diurnal temperature range (DTR) except in Lamphun Province (Figure 2).Zhang et al. [51] and Liu et al. [52]  All rainfall extremes showed insignificant trend except consecutive dry days (CDD) which was increasing with significant trend during the observed period of 1960-2010 (Table 3) in Chiang Rai and Tak provinces.Similar results were observed in another study of Sharma and Babel [50].Except Lamphun Province, negative trends of PRCPTOT and R10 indicate that the rainfall was decreasing (Figure 2) over the study area.This less rainfall distribution supports the negative trend of consecutive wet days (CWD).

Seasonal Variation in Observed Period.
During observed period, maximum temperature ( max ) showed a significant increasing trend in all seasons except spring which has an insignificant negative trend in all provinces.Minimum temperature ( min ) increased significantly in Chiang Rai Province in all seasons, which is followed by Tak, Chiang Mai, and Lamphun provinces.The Mann-Kendall test statistic  showed that  min has increasing trend more than  max .In mountainous provinces like Chiang Rai and Chiang Mai, the  The ratio of annual total rainfall to the number of wet days (≥1 mm) mm/day temperature is observed to be higher than other provinces.A similar trend was identified by Zongxing et al. [41].They also noticed the influence of topography on the regional pattern of temperature indices in Southwestern China; significant trends were on the Eastern Xizang Plateau and the Hengduan Mountains, while stations with insignificant trends were in the Sichuan Basin and on the Yunnan-Guizhou Plateau.
There was no significant trend of rainfall detected during the observed period except winter (downward trend) and spring (upward trend) in Chiang Rai and Chiang Mai provinces, respectively (Table 4).

Model Performance Evaluation.
Climate parameters can theoretically be projected by using appropriate GCM and RCM for intended place, but this does not imply that all the data will perfectly correspond to the real situation.Some data may be consistent with a particular time of month or year while other may not.To evaluate model performance, after screening all the 26 GCM variables in the SDSM model for significant correlation, 12 out of 26 variables were found significantly correlated ( < 0.05) to the observed temperature and rainfall data.Later on, these variables were used for the calibration of the model to develop a statistical relationship with the observed data and GCM output.For checking the consistency of observed and simulated maximum and minimum temperature and cumulative rainfall data, 1980s was taken as the base period.Results (Figure 3) revealed that none of the SRES of two models is fully unswerving with the observed data.However, based on the coefficient of determination, both PRECIS RCM and HadCM3 GCM have consistency with the base period.Similar performance was observed for HadCM3 GCM, which was downscaled using SDSM model in another study by Deb and Babel [53].The maximum temperature difference between the observed and simulated values was about ±0.5 ∘ C, which imply that two models are reliable to simulate future data.Both models showed higher accuracy in case of predicting the rainfall compared to the base period.The average difference was ±60 mm for all provinces except Tak which shows the maximum difference for HadCM3 GCM and minimum for PRECIS RCM.This performance evaluation is not to select the best model to simulate future scenarios rather than to + Moderately significant, * significant, * * highly significant, and * * * very highly significant.
Moderately significant = 90% level of confidence, significant = 95% level of confidence, highly significant = 99.9% level of confidence, and very highly significant = 99.99% level of confidence.provide more alternatives to know about the future climate and their extremes.plotted in Figure 4 for all provinces.None of the SRES is considered to have an outcome more likely than another but it is obvious that all demonstrate the continuation of the trends demonstrated by the historical data, which indicates increasing rates of temperature change for North Thailand.
There are some exceptions like when HadCM3 GCM predicts a flat trend of maximum temperature, which is below the observed average maximum temperature.Similarly, PRECIS RCM did the same for Talk Province.These exceptions are similar for minimum temperature.The ensemble averaged minimum temperature in the future period will increase by 0.87, 1.42, 0.93, and 0.49 ∘ C relative to the observed period in Chiang Rai, Chiang Mai, Lamphun, and Tak, respectively.For both models and scenarios considered herein, the average monthly maximum and minimum temperature for all provinces average increases projected for the 2011-2099 period are 0.99 and 0.93 ∘ C above the observed period (1960-2010) mean of 32.26 and 20.66 ∘ C, respectively.Erda et al. [54] reported that the temperature increase in China by the end of the twenty-first century may be between 3 and 4 ∘ C, simulated by A2 scenario of PRECIS model.In Southeast Asia, mean warming of temperature will be 2.5 ∘ C in 2080-2099 relative to 1980-1999 [55].According to IPCC ensemble-mean predictions, results estimated a net annual increase in temperature between 1.32 ∘ C (lowest emissions scenario B1) and 2.01 ∘ C (highest emissions scenario A1F1) for the 2040-2069 periods relative to the baseline 1961-1990 period [56].It is clear from the above discussion that all scenarios indicate higher future temperature for North Thailand than what would result from a continuation of the trend from the historical record.Multimodel and multiensembles temperature anomaly, regionally averaged, showed that the trend of increasing monthly minimum temperature is higher than that of monthly maximum temperature (Figure 5).As a result of higher interannual variability, prediction and identification of trends are less apparent for rainfall than temperature.However, each model and scenario evaluated project an increase in average annual rainfall of 9.65% for 2011-2099 compared to the observed average annual rainfall 1235 mm.Total annual rainfall projections for the range of SRES suggest small differences (Figure 6) for Chiang Rai.For the period 2011-2099, the HadCM3 GCM and PRECIS RCM models project increase in total annual rainfall by 20.95% and 18.46%, respectively, above the observed period mean of 1725 mm in Chiang Rai Province.PRECIS RCM projects increases in total rainfall but HadCM3 GCM projects decreases for Chiang Mai and Lamphun.This state is opposite for Tak Province where HadCM3 GCM projects increases in total rainfall.Future projection of monsoon rainfall is quite complicated due to the presence of aerosols (carbon aerosols) in the Asian monsoon [55].However, the rainfall anomaly, regionally averaged, showed a continuous increasing trend of annual rainfall (Figure 7).

4.5.
Projected Seasonal Variation.This section explains the ensemble average seasonal temperature and rainfall variation for 2020s, 2050s, and 2080s compared to the observed base period (1980s) (Tables 5 and 6).Daily maximum temperature in winter is projected to decrease during 2020s with respect to the base period.The daily maximum temperature projections of A2 and B2 scenarios of PRECIS RCM increase at higher rates in winter than that of the HadCM3 GCM and for the 2080s both scenarios of PRECIS RCM project daily maximum temperature increases by 2.33 and 1.40 ∘ C, respectively.Similarly, both models project an increasing daily maximum temperature for spring, summer, and autumn seasons except 2020s, where only PRECIS RCM showed a decline of daily maximum temperature.Daily minimum winter temperature is going to increase in 2020s followed by 2050s and 2080s compared to the base period 1980s.The maximum increment (4.99 ∘ C) is found for A2 scenario of PRECIS RCM in 2080s.In spring, PRECIS RCM projects an increase in daily minimum temperature in 2020s, 2050s, and 2080s regardless of scenarios but HadCM3 GCM projects the opposite.In summer and autumn, this minimum temperature is projected to increase in 2050s and 2080s by both scenarios of PRECIS and HadCM3 models.It is noteworthy that the increasing rate of daily minimum temperature is higher during winter and summer than that of daily maximum temperature, whereas this fact is opposite in case of spring and autumn.
Seasonal variations of annual rainfall in the future decades with respect to the base period 1980s are tabulated in Table 6.In winter, A2 and B2 scenarios of HadCM3 GCM project an increase in annual rainfall by 80% in 2080s, whereas PRECIS RCM's predictions are not consistent with this result.However, each scenario and model evaluated project a decrease in annual rainfall ranging from 13.84 to 61.10% in spring.However, in summer, the projection is opposite to spring where annual rainfall is increasing in 2020s, 2050s, and 2080s ranging from 19.83 to 46.55%.For autumn, the annual rainfall projections for the range of   SRES suggest big differences, because PRECIS RCM projects an increase in rainfall whereas HadCM3 GCM projects the opposite.Although both temperature and rainfall are rising in the future, there is no significant correlation found for temperature and rainfall in this region [57].
4.6.Trend Analysis for the Future Period.Figure 8 illustrates the two IPCC SRES for Chiang Rai Province as predicted by PRECIS RCM and HadCM3 GCM to see the trend of different warm and cold temperature extremes, developed by RClimDex, for the period of 2011-2099.Mann-Kendall test statistics  with Sen's slope are provided in Table 7 for the same province.Almost all the temperature rising extreme indices (TX90p, TN90p, TXx, TNx, SU35, TR25, and WSDI) are increasing (Table 7 and Figure 8) with highly (very highly) significant trend, and the temperature falling indices (TX10p, TN10p, TXn, TNn, and CSDI) are decreasing with same trend for both SRES A2 and B2.Also, the values of Mann-Kendall test statistic  and Sen's slope indicate that minimum temperature extremes (TNx and TNn) have stronger trend and steeper slope than maximum temperature extremes (TXx and TXn) and that the cold extremes (TX10p and TN10p) have weaker trend and milder slope than warm extremes (TX90p and TN90p).On an average, maximum and minimum temperature will increase by 1.56 and 2.24 ∘ C, respectively, in 2090s as compared to 2020s (Figure 8).In contrast to summer days (SU35) and tropical nights (TR25), tropical nights will increase far beyond summer nights with the same level of significance.Hence, minimum temperature will increase more than the maximum temperature which is supported by the very high significant negative trend of diurnal temperature range (DTR); that is, daynight temperature variation will be smaller in future.Daily minimum temperatures are projected to increase faster than daily maximum temperatures, leading to a decrease in diurnal temperature range in the world [55].
The summary of the statistical analysis of rainfall extreme indices with their level of significance is tabulated in Table 8 + Moderately significant, * significant, * * highly significant, * * * very highly significant, and ( ) = Sen's slope.
and their temporal trends are illustrated in Figure 9.The consecutive dry days (CDD) are indicating less rainfall in the study area.This extreme index is projected to decrease by an ensemble average of 1.38 days/year.The trend of annual total rainfall (PRCPTOT) reveals a very high significant increasing trend except B2 scenario of PRECIS RCM.This indicates that the study area will receive more rainfall in the future.Number of days with rainfall ≥10 mm (R10) and ≥20 mm (R20) are indicators of moderate and heavy rainfall.
Both of these indices are projected to increase significantly.However, R20 will increase with higher trend than that of R10.Two important measures of extreme rainfall are the percentage of total rainfall due to events above the 95th and 99th percentiles (R95p and R99p).Results revealed that both indices increase with the significant trend.The sign of the trend for the annual total rainfall is compared with the sign of the trends in the R99p index.At a station where the annual amount decreases, negative R99p trends are indicative of a large contribution of the extremes and vice versa [50].
The maximum event intensity indices (RX1day and RX5day) show a significant increasing trend, whereas all scenarios of both models show a very high significant trend for RX5day    + Moderately significant, * significant, * * highly significant, and * * * very highly significant.compared to that of RX1day.The simple daily intensity index (SDII) is the combined effect of the trend of PRCPTOT and consecutive wet days (CWD).All scenarios of both models predict an increase of SDII with nonsignificant trend except A2 scenario of HadCM3 GCM.Table 9 provides the summary of the statistical analyses of rainfall and temperature indices with their significance for Chiang Mai, Lamphun, and Tak provinces.With some exception especially for rainfall, the future projections for all temperature and rainfall indices are similar to the projection for Chiang Rai Province.

Conclusions
Our analysis enables a number of conclusions to be drawn as follows: (1) During observed 1960-2010 period, the entire temperature extremes showed a warming trend which is also supported by the statistical analyses.A warming trend is indicated by the summer days, tropical nights, Max  max , Min  max , Max  min , Min  min , cool nights, cool days, warm nights, warm days, diurnal temperature range, warm spell duration indicator, and cold spell duration indicator.The warming trend is higher in the southern part of the study area, that is, in Tak.
Trend and slope of night extremes (TNx, TNn, TN10, and TN90) appear to have been higher than those in day extremes (TXx, TXn, TX90, and TX10) resulting in the decrease of day-night temperature variation which is also indicated by diurnal temperature range (DTR).( 2) Considering all seasons (winter, spring, summer, and autumn), maximum temperature ( max ) and minimum temperature ( min ) have increased significantly in all seasons except spring which is insignificant over the whole study area.Statistical analyses indicated that the increasing trend of  min was higher than that of  max .(3) In future, the maximum temperature beyond 35 ∘ C (summer days) and minimum temperature beyond 25 ∘ C (tropical nights) will increase with the very high significant trend.For all SRES (A2 and B2) of both models considered in this study, the maximum and minimum temperature increases projected for the period of 2011-2099 were 0.99 and 0.93 ∘ C above the observed period.However, monthly maximum and minimum temperature will boost up by 1.56 and 2.24 ∘ C at the end of this century as compared to the 2020s.In future, all the day, night, cool, and warm extremes will change same as observed period.Diurnal temperature range will go down (0.64 ∘ C) notably over the study area.(4) During the observed period, annual total rainfall (PRCPTOT) and moderate rainfall (R10) indices were  decreased which is supported by the negative trend of consecutive wet days (CDD).However, none of the indices showed significant trend in the observed period and there was no significant trend found for seasonal rainfall also.Each model and scenario evaluated project an increase in average annual rainfall of 9.65% for 2011-2099 compared to the observed annual rainfall 1235 mm.
Topographic changes of the earth greatly influence the variation of climatic parameters and their extremes and anomalies.Hence, it is recommended to explore the changes of climate extremes and their anomalies over topography using a large number of stations.

Figure 1 :
Figure 1: Location map of the study area in Thailand.

Figure 2 :
Figure 2: Monthly indices of temperature and rainfall extremes for the observed period 1960 to 2010 averaged across the study area.The dotted line is a trend-line computed by least square fit and corresponding regression equation is showed for each index.

Figure 3 :
Figure 3: Comparison of observed and simulated cumulative rainfall (a), monthly maximum temperature (b), and monthly minimum temperature (c) for 1980-1989 of Chiang Rai.Results for the case of PRECIS RCM are shown in blue color and those for the case of HadCM3 GCM in red.

4. 4 .Figure 4 :Figure 5 :
Figure 4: Average monthly maximum and minimum temperature variation in the future for (a) Chiang Rai, (b) Chiang Mai, (c) Lamphun, and (d) Tak.Black solid thick line up to 2000s indicates the observed data.

Figure 6 :
Figure 6: Annual rainfall variation in the future for (a) Chiang Rai, (b) Chiang Mai, (c) Lamphun, and (d) Tak.Black solid thick line up to 2000s indicates the observed data.

Figure 7 :
Figure 7: Multimodel and multiensembles rainfall anomaly with respect to the control period 1970 to 2000.Results for the case of observed period (observed data) are shown in black and those for the case of future period (simulated data) in blue.The black solid is representing the 5-year moving average of anomalies.

Figure 8 :
Figure 8: Annual series of temperature extremes for Chiang Rai province as predicted by PRECIS RCM and HadCM3 GCM.

Figure 9 :
Figure 9: Annual series of rainfall extremes for Chiang Rai province as predicted by PRECIS RCM and HadCM3 GCM.

Table 1 :
Detail information for meteorological station.

Table 2 :
Definition of extreme temperature and rainfall indices used in this study.

Table 3 :
Observed trend of temperature and rainfall indices.

Table 4 :
Observed seasonal trend of temperature and rainfall.

Table 5 :
Average seasonal variation of daily maximum and minimum temperature in the future relative to 1980s of observed period.

Table 6 :
Average seasonal variation of annual rainfall in the future relative to 1980s of observed period.

Table 7 :
Trend of temperature indices for 2011-2099 period in Chiang Rai.

Table 8 :
Trend of rainfall indices for 2011-2099 period in Chiang Rai.

Table 9 :
Trend of temperature and rainfall indices for 2011-2099 period.