Forecasting the Future Temperature Using a Downscaling Method by LARS-WG Stochastic Weather Generator at the Local Site of Phitsanulok Province, Thailand

The study evaluates the effect of climate change on temperature, which is one of the most important variables in water resources management and irrigation scheduling. Climate prediction is necessary in the agricultural and hydrological analysis. This study proposed an approach to the application of the Long Ashton Research Station Weather Generator (LARS-WG) in Coupled Model Inter-comparison Project Phase 5 (CMIP5) under EC-Earth and MPI-ESM-MR. The first step is model calibration, where the observed dataset is analyzed statistically. In the second stage, the synthetic data and observed data are checked for Kolmogorov-Smirnov and the means and standard deviations. In order to evaluate the response of temperature under future warmer climate trends, the approach was assessed using data series. These parameters consisted of the minimum and maximum temperature at the Phitsanulok Meteorological Station (WMO Index 48378) and RCP4.5 climate change scenario for the base period as well as for 2021-2040 (the near future), 2041-2060 (the medium future) and 2061-2080 (the far future). The results of the numerical applications indicated that the linkage between the observed data spatially downscaled from LARS-WG simulations with the historical one of the locations during the baseline period had a very good accuracy. It was also found that the future climate change of temperature contributed to higher change. The mean of minimum temperature in the baseline year was 23.13 ̊C while the mean of minimum temperature in the projection period for 2021-2040, 2041-2060 and 2061-2080 is expected to be 24.09 (+4.18%), 24.49 (+5.94%) and 24.82 (+7.36%) ̊C, and 24.12 (+4.32%), 24.82 (+7.36%) and 25.08 (+8.48%) ̊C for the EC-Earth and MPI-ESM-MR, respectively. While, the mean of maximum temperature in the baseline year was 33.41 ̊C, the How to cite this paper: Punyawansiri, S. and Kwanyuen, B. (2020) Forecasting the Future Temperature Using a Downscaling Method by LARS-WG Stochastic Weather Generator at the Local Site of Phitsanulok Province, Thailand. Atmospheric and Climate Sciences, 10, 538-552. https://doi.org/10.4236/acs.2020.104028 Received: August 12, 2020 Accepted: October 9, 2020 Published: October 12, 2020 Copyright © 2020 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access S. Punyawansiri, B. Kwanyuen DOI: 10.4236/acs.2020.104028 539 Atmospheric and Climate Sciences maximum temperatures are projected to increase at 34.47 (+3.19%), 34.88 (+4.43%) and 35.21 (+5.40%) ̊C, and 34.53 (+3.36%), 35.19 (+5.34%) and 35.30 (+5.67%) ̊C, respectively. Furthermore, the future local surface temperatures from the MPI-ESM-MR project tended to be higher than EC-Earth. In conclusion, the study results indicate that in coming three time periods, the minimum and maximum temperature increase is expected in Phitsanulok province, Thailand.


Introduction
According to existing evidence as well as the fluctuating influences of greenhouse gases and human activities, climate is a complicated system that is continuously evolving. Appropriate approaches for decreasing the fatal effects of climate change phenomena are possible by forecasting and examining the rate of change in meteorological considerations. Employing modelled climatic variables by General Circulation Models or Global Climate Models (GCMs) is one of the most operative methods for assessing the influences of climate change occurrences [1] [2]. These simulations can forecast the changes resulting from the impact of greenhouse gasses on weather conditions. The output of these models should be reduced when they show low spatial resolution power. Awareness of potential climate change leads to the adoption of mitigation strategies, as these effects have implications for areas of public interest including agriculture [3] [4], hydrological [5] [6] and water management [7] [8].
The IPCC describes GCMs as a computational (quantitative) description of the climate system on the basis of the physical, chemical and biological characteristics of its elements, their correlations and feedback mechanisms [9]. GCMs, centered on statistical models of the atmosphere, ocean and land surface systems, are thought to be the only reliable instruments presently offered to model the reaction of the global climate system to rising greenhouse gas levels. Because of the restricted depiction of scale atmospheric processes, topography and land-sea distribution in GCMs, however, the direct use of output from GCMs is frequently insufficient [10].
Potential climate forecasts of GCMs in various conditions are typically used to establish climate mitigation approaches. Current GCMs still have insufficient ability to predict complex and local climate characteristics or provide accurate information on temperature and rainfall, which is the key input for hydrological effect evaluation criteria. Additionally, the contributions given by GCMs are too imprecise to be used by these hydrological models as they involve details at finer scales. Downscaled GCM results are typically used to present the fine-resolution or point-scale information needed for impact simulations [11].

The Study Area
The study site was Phitsanulok Meteorological Station (PMS) which is located at latitude 16.783˚N, longitude 100.267˚E and altitude 44 meters above mean sea level, as shown in Figure 1. This area supplies much of the agricultural produce of Thailand, and is under intense demand pressures regarding climate change for agricultural and irrigation system.

LARS-WG Stochastic Weather Generator
LARS-WG is a stochastic weather generator which could be utilized to predict weather data at a single station under both present and potential atmospheric conditions. The data comprises the regular time series for a number of environmental factors, including rainfall, maximum temperature and minimum temperature, and solar radiation. LARS-WG employs a semi-empirical distribution (SED) known as the cumulative probability distribution function (CDF) to estimate the confidence interval for minimum and maximum temper- EPM is a distribution histogram of 23 distinct intervals (a i-1 , a i ) where a i-1 < a i and gives a more precise picture of the observable distribution compared to 10 used in the earlier models. By disrupting the distribution parameters for a site with projected climate change from global or regional climate simulations, a regular environment model for this site may be created and used in tandem with a process-based impact assessment model. Generally, the approach for producing synthetic weather data can be divided into three different steps: model initialization, model validation, and scenario creation, as described and briefly defined below [19]. In the case of the minimum and maximum temperature two values close to 0 and 1 are used to interpret for extremely low and high temperatures, i.e. p 2 = 0.01, p 3 = 0.02, p n-1 = 0.99 and p n-2 = 0.98. Both values of pi (0 < I < n). For the latest edition of LARS-WG 6.0, the low and high temperatures for dry and rainy days are estimated by SEDs measured for each month [20]. The contribution to LARS-WG is a collection of seasonal data (rainfall, low temperature and high temperature) from the base time  and site data (latitude, longitude and altitude). The perceived daily weather data at a given site was employed to establish a set of parameters for the probability distributions of weather variables after input data preparation and quality control. Such criteria are used to produce a simulated weather sequence of variable lengths by randomly choosing values from the suitable distributions, with the same statistical properties as the original data observed, but varying on a day-to-day level. LARS-WG determines wet days from dry days depending on whether or not rainfall is zero. Alternating wet and dry series estimated by semi-empirical probability distributions are used to model the existence of rainfall. The statistical properties of the actual and simulated weather data during the model testing were evaluated to assess whether there were any statistically important variations using the Chi-square strength test (Kolmogorov-Smirnov, KS) and the means and standard deviations using the t and F measures, respectively. By varying the constraints of LARS-WG, number of years and seed numbers can be conducted. The LARS-WG baseline constraints, which are analyzed from the perceived weather for a baseline period  are regulated by the delta-changes for the future period and the emissions projected by a GCM for each climatic variable covering the site for creation of climate scenarios at the site for a particular time in the future with a designated emissions scenario.
In this study, the local-scale climate scenarios based on RCP4.5 scenarios simulated by the selected two GCMs are generated for the time periods of 2021-2040, 2041-2060 and 2061-2080 to predict the future change of temperature in PMS. The coefficient of determination for a linear regression model is the quotient of the variances of the generated data and observed data values [21]. The coefficient of determination is computed according to the Equation (2) below: The coefficient of determination (R 2 ), represents the proportion of the total sample variation in y (measured by the sum of squares of deviations of the sample y values about their mean y) which is explained by the linear relationship between generated data and observed data.

Global Climate Models (GCMs)
The new version of the LARS-WG 6.0 was applied for this study that incorporated predictions from two GCMs, which were used in the IPCC's Fifth Assessment Report (AR5) based on Emissions Scenarios RCP4.5 (medium greenhouse

Calibration and Validation of LARS-WG
The "Site Analysis" function in the LARS-WG model with two data sets, comprising observed station and GCMs data, were utilized to conduct the calibration and validation. The Kolmogorov-Smirnov test, F-test and t-test were operated to assess the functioning of the weather generator during calibration and validation. To assess the similarity of distributions for daily maximum and minimum temperature, the KS test was used. The equivalence of 30-day mean of daily maximum temperature (Tmax) and the 30-day mean of daily minimum temperature (Tmin) were assessed by the F-test. Both experiments measure a p-value that is used to support or dismiss the hypothesis that the two data sets (detected and created) may have been obtained from the same distribution at a 5 percent level of significance. The coefficient of determinant (R 2 ) was used to check the performance. Tables 2-6 show the estimation of suitability for LARS-WG function in replicating for PMS. The model executed very well in fitting for the two datasets, as shown by the KS test. LARS-WG performance was effective in simulating both minimum and maximum temperature for both data sets, as revealed by the assessment results. Figures 2-5 show the relationship linking the 30-day mean and standard deviation of temperature for the two data sets employed in the study. For measuring the 30-day means of Tmin and Tmax statistics, the outcomes displayed excellent performance by LARS-WG. Using LARS-WG, the mean periodic totals for the minimum and maximum temperatures were formed efficiently, which showed that the temperature can be analyzed using the daily time series. The findings of the experimental analyses showed that there were no major variations in the calculated mean of the measured average temperature relative to the estimates. LARS-WG precisely created the 30-day means of minimum and maximum temperature values, offering a coefficient of determinants equivalent to units. As shown in Table 2, the R 2 for the mean 30-day minimum and maximum temperature had a robust linear correlation linking the experimental, created and artificial data.     Compared to those calculated from the experimental time series, the means and variances of day-to-day artificial weather data were assumed to be insignificantly distinctive. It is likewise critical that the artificial weather series adopt a confidence interval that is not significantly varied from the actual time series. As seen in Figure 4, Figure 5, the 30-day temperature standard deviations for the created values were miscalculated for both data sets (perceived and produced).
The evaluation results exhibit that LARS-WG execution in reproducing the distributions of the day-to-day temperature in all months is ideal, as seen from the results in Table 3 and Table 4.

Future Scenarios of Minimum and Maximum Temperature
The EC-Earth and MPI-ESM-MR GCMs and RCP4.5 scenarios in LARS-WG version 6.0 were used in this study to generate future climate scenarios to better deal with un-certainties. In Figure 6, Figure 7, graphs of the temporal evolution of the indicator Tmin and Tmax in the PMS part of Phitsanulok Province of Thailand are presented. Table 5, Table 6

Conclusion
Using a LARS-WG 6.0 stochastic weather generator, this study replicated the meteoro-logical constraints from Phitsanulok meteorological station (minimum and maximum temperature). The findings revealed that these criteria were modelled with very high precision. LARS-WG may be used to produce climate mod- likely to be superior to EC-Earth. It was concluded from research that the LARS-WG model had the reasonable ability to imitate the day-to-day and 30-day temperature constraints and could be used as an operative tool for integrating climate change impacts into viable improvement. The work presented in this paper strongly supports the use of multi-GCMs model downscaling for providing the required data for agricultural and hydrological analysis. Additionally, we found that one of the many reasons for the large differences of the future minimum and maximum temperature is the CMIP5 models' resolution. Finally, the LARS-WG 6.0 provides a means of synthesizing long-term climate data that can then be used in impact of climate change assessment. Evaluate the performance of the generators for multi-GCMs and climate models are required to provide perspective into downscaling model vulnerabilities in order to establish alternate mitigation and adaptation approaches in climate change.