Forecasting the Future Temperature Using a Downscaling Method by LARS-WG Stochastic Weather Generator at the Local Site of Phitsanulok Province, Thailand ()
1. 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].
A computer algorithm that incorporates weather and climate records to generate a long series of synthetic daily weather data of continuous integration for a location due to the statistical features of the studied weather data for that location is a stochastic weather generator. One of the approaches for advancing local scale future climate situations from large-scale climate changes modelled by GCMs is stochastic weather generators [12]. Recently, establishing the prospective influence on hydrological applications has been carried out in climate change and variability studies using simulated weather data created by weather generators [13] [14].
These selected weather generators have been used in different climate change impact studies, including studies with LARS-WG weather generators [15]. The current version 6.0 of LARS-WG incorporates climate projections from the CMIP5 ensemble used in the IPCC Fifth Assessment Report (AR5) [16]. Previous version of LARS-WG 5.5 which climate projections from the CMIP3 provided the scientific base for the Fourth Assessment Report (AR4) of IPCC published in 2007 [17]. The LARS-WG 6.0 has been well validated in diverse climates around the world. The LARS-WG version 6.0 can be download from https://www.rothamsted.ac.uk/, the link for downloading LARS-WG 6.0 is given in https://sites.google.com/view/lars-wg/.
The objectives of this study were to use the new version of the LARS-WG in predicting climate change for the periods 2021-2040 (the near future), 2041-2060, (the medium future), and 2061-2080 (the far future) according to the projections of EC-Earth and MPI-ESM-MR GCMs for RCP4.5 scenarios (medium greenhouse gas emission scenario) at the local site of Phitsanulok province, Northern of Thailand.
2. Data and Methods
2.1. 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.
The site is influenced by the tropical region which has two seasons: the dry season from November to April, and the rainy season from May to October. The average minimum and maximum temperature were 22.73˚C, 33.26˚C in 1961-2019.
Historical climate data for daily precipitation, the minimum and maximum temperature were obtained from PMS (WMO Index 48378) for the period 1961-2019. Only data for the period 1980-2010 (the baseline year) were selected for analysis as required for the calibration and validation of the weather generator. The weather data did not include solar radiation. The climate data of 3.17% were missing from the observed station data. Missing data in the observed dataset may alter the statistics, but the long-term climate series will still be generated, missing data values should be coded −99.0 of the layout of a typical file of weather data for use in LARS-WG [18]. The generator outputs result of statistical tests comparing the observed and generated data.
2.2. 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 temperatures.
. (1)
EPM is a distribution histogram of 23 distinct intervals (ai-1, ai) where ai-1 < ai 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. p2 = 0.01, p3 = 0.02, pn-1 = 0.99 and pn-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 (1980-2010) 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 (1980-2010) 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:
(2)
The coefficient of determination (R2), 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.
2.3. 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 gas emission scenario) for three times the period as listed in Table 1. These two models differ in spatial resolution power, design institute, predictability of atmospheric variables, and predictability of oceanic variables [22] [23]. The simulation of GCMs assumes the year length is 365 calendar days, with global coverage and daily.
3. Results and Discussion
3.1. 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 (R2) 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 R2 for the mean 30-day minimum and maximum temperature had a robust linear correlation linking the experimental, created and artificial data.
Table 1. The global climate models from IPCC AR5.
Table 2. Coefficient of determination (R2) in calibration.
Table 3. KS-test of daily distributions for observed data.
Table 4. KS-test of daily distributions for generated data.
Table 5. Future mean of minimum temperature (˚C) under RCP4.5 scenarios.
Table 6. Future mean of maximum temperature (˚C) under RCP4.5 scenarios.
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.
Figure 2. Tmin monthly mean of observed versus generated data.
Figure 3. Tmax monthly mean of observed versus generated data.
Figure 4. Tmin monthly standard deviations of observed versus generated data.
Figure 5. Tmax monthly standard deviations of observed versus generated data.
3.2. 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 show the mean of Tmin and Tmax in base time (1990s) would be 23.12 and 33.41 (˚C). The Tmin would increase (˚C) by 0.97, 1.37 and 1.70 under EC-Earth RCP4.5 (2030s), RCP4.5 (2050s) and RCP4.5 (2070s), respectively. From the MPI-ESM-MR, results showed the increase (˚C) in the mean Tmin by 1.00 under RCP4.5 (2030s) after that an increase in the mean of Tmin by 1.70 and 1.96 under RCP4.5 (2050s) and RCP4.5 (2070s), respectively.
The mean of Tmax would increase (˚C) by 1.06, 1.48, and 1.80 under EC-Earth RCP4.5 (2030s), RCP4.5 (2050s) and RCP4.5 (2070s), respectively. While the MPI-ESM-MR, results showed the increase (˚C) in the mean of Tmax by 1.12, 1.78 and 1.89 under RCP4.5 (2030s), RCP4.5 (2050s) and RCP4.5 (2070s), as in Table 6. The changes in temperature are 0.97˚C - 1.80˚C (EC-Earth), 1.00˚C - 1.89˚C (MPI-ESM-MR) under RCP4.5 scenarios generated using the observed station data.
Figure 6. Baseline and future of minimum temperature (A) EC-Earth, (B) MPI-ESM-MR.
Figure 7. Baseline and future of maximum temperature (A) EC-Earth, (B) MPI-ESM-MR.
A step towards a widespread evaluation of uncertainty in future climate predictions is possible with the multi-ensembles method using various climate models and emissions circumstances. Nonetheless, the relative change factors (RCFs), resulting from the direct GCM output of the variables that precisely match the estimates, are used by the LARS-WG. Therefore, the reality that GCMs are not very accurate when replicating temperature means that the error caused by the GCM output for temperature extends the error of downscaling that causes the efficiency of LARS-WG to downscale temperature to demand more vigilance [24].
4. 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 models for present and future PMS simulations. LARS-WG projected the 30-day mean low and high temperatures, which are consistent with the similarity between the observable and produced mean values of 0.999 and 0.998 respectively. The low and high temperature for PMS would rise in the 2030s, 2050s and 2070s under RCP4.5, as revealed by the results. Climate signals for temperature values in the 2030s, 2050s, and 2070s would rise under perceived station data, as suggested by the ensemble mean for the EC-Earth and MPI-ESM-MR GCMs. Further, impending local surface temperatures from the MPI-ESM-MR project were 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.
Acknowledgements
The authors would like to thank the Agricultural Research Development Agency, Thailand for supporting the study. The researchers also wish to thank Prof. Mikhail A. Semenov for the provision of Long Ashton Research Station Weather Generator, the license used in this study, and the Thai Meteorological Department (TMD) for providing the weather data.