Fusion of Simulated and Observational Temperature Data in the Beijing-Tianjin-Hebei Region Based on High-Accuracy Surface Modeling

Temperature changes have a major impact on all aspects of human society and have attracted global attention. 'e scarcity of observation data and the inaccuracy of the models make obtaining accurate temperature distributions a challenge. 'is study introduces high-accuracy surface modeling (HASM) combined with temperature simulations from the Weather Research and Forecasting (WRF) model and temperature records from observation stations to investigate the spatiotemporal characteristics of temperature in the Beijing-Tianjin-Hebei region during the period of 1956–2005. Leave-one-out cross-validation is applied to verify the temperature fields before and after the fusion of the models.'e results indicate that theWRFmodel has a limited ability to simulate temperature conditions, but the overall deviation across the region is relatively large. 'e fusion results of the HASM decrease the mean absolute error (MAE) and the root-mean-square error (RMSE) by half in most instances, and the correlation between the fusion data and observation data is approximately 0.01–0.03 higher than that with theWRF simulation data. Based on the fusion data, obvious warming trends are observed during 1976–2005. In general, temperatures in spring, summer, and autumn increase rapidly from 1996 to 2005 but from 1976 to 1995 in winter. 'e substantial fluctuations in the interannual temperature during 1996–2005 in summer, autumn, and winter may be related to the frequent occurrence of extreme weather. 'e spatial distribution of temperature change differs in each season during 1956–1995. A dramatic increase in temperature occurs in the western part of the study area during 1996–2005 but with no seasonal difference.


Introduction
e impact of human activities on the natural environment has led to global and regional climate change, which has a major impact on all aspects of human society, such as the economy, agriculture, fishing industry, hydraulic engineering, and relevant policymaking [1][2][3][4].Temperature is the most direct meteorological factor for the characterization of climate change.
e Intergovernmental Panel on Climate Change (IPCC) fifth report [5] noted that the warming rate at the global scale reached 0.12 °C/10 yr from 1951 to 2012, which was accompanied by fewer cold days and nights and more warm days and nights.At the same time, extreme high/low-temperature weather and hydrological disasters, such as droughts and floods, occur frequently all over the world.
Human activities, such as urbanization and returning farmlands to forests, can greatly change the land cover, and greenhouse gas emissions from industrial activity, vehicles, and heating in winter can increase the concentration of greenhouse gas in the air; both of these activities drive climate change.ere are also regional differences in climate change due to differences in human aggregation and activity [6][7][8][9].China is a populous country with a fast-growing economy.From 1989 to 2018, China's gross domestic product (GDP) grew at an average annual rate of 9.55% [10], and the population also grew rapidly [11].In the Beijing-Tianjin-Hebei region, which is the political, economic, and cultural center of China and an urban agglomeration, land cover has changed dramatically in recent decades.Factors, such as urban expansion, farmland reduction, artificial vegetation, greenhouse gas emissions, aerosol changes, and the heat island effect, have large impacts on the climate in this region.
ese influences also affect socioeconomic factors such as people's livelihood and security [12,13].Studying the climate change characteristics of the Beijing-Tianjin-Hebei region also provides a basis for the formulation of relevant policies [14].
Many studies on spatiotemporal temperature patterns have been conducted in the Beijing-Tianjin-Hebei region.Yan et al. [15] integrated observation data from the Beijing Observatory Station and reported that the observed temperature data of the station were likely to be approximately 0.48 °C/10 yr.Li and Yao [16] performed an inhomogeneity analysis for observation data using the Multiple Analysis of Series for Homogenization (MASH) software package.Zhao et al. [17] used 327 observation stations to explore the temporal and spatial variations in temperature in the Beijing-Tianjin-Hebei region for 50 years and found that the annual temperature over 10 years largely changed in the southeast and that interdecadal variations of temperature were spatially stationary and varied from −0.05 °C to 0.05 °C.Because of the spatial discontinuity of the observed dataset, regional weather models have also been used to study climate change.For instance, Miao et al. [18] used the Weather Research and Forecasting (WRF) model to study the effects of urbanization on surface weather and the boundary layer structure and found that the urban canopy variables simulated by the model were in good agreement with the ground observations.Liu et al. [19] ran sensitivity tests on the regional weather model and suggested that heat waves enhanced urban warming intensity at night.Based on the Regional Climate Model version 3 (RegCM3), Zhang and Shi [20] studied the changes in precipitation, temperature, and certain related extreme weather events in northern China under the background of global warming.However, limited by the lack of understanding of atmospheric motion and human activities, biases remain in the results of the numerical weather models [21,22].e objective of this paper is to perform the highaccuracy surface modeling (HASM) to fuse the observation data from stations and the simulations from numerical weather models.At present, there is no study using a dynamical downscaling model as the input conditions of HASM to simulate the surface temperature.e previous studies performed statistical downscaling such as Bayesian statistical methods, linear regression methods, interpolation, and geographical weighted regression to produce the input conditions of HASM [17,[23][24][25].From the mechanism of temperature distribution, a dynamical downscaling model has more realistic and reliable results [26] and improves the accuracy of the final results of HASM.We use the WRF model to downscale the Community Earth System Model (CESM) climate data from 1956 to 2005 and apply these data to drive the HASM.en, the time continuous observation data of the Beijing-Tianjin-Hebei region during the study period are screened and used as the accuracy control conditions of HASM.By running HASM, we obtain the fused temperature fields and compare the accuracy of the temperature data before and after the fusion.Finally, the temporal and spatial variations in temperature in the Beijing-Tianjin-Hebei region from 1956 to 2005 are analyzed based on the HASM fusion results.

Study Region.
e location of the Beijing-Tianjin-Hebei region is shown in Figure 1(a) (between 36 °N-42.5 °N and 113 °E-120 °E). is region incorporates Beijing, Tianjin, and Hebei Provinces.Beijing and Tianjin have high levels of urbanization, and the central areas of these two cities are dominated by underlying urban areas.
e elevation of Hebei Province is high in the northwest and low in the southeast, with complex and diverse landforms.Population is mostly distributed in the low-elevation area, where the land cover is dominated by cities and farmland.e northwest region of Hebei Province is covered by mountain forests and grasslands.e Beijing-Tianjin-Hebei region is a key economic zone in China that is affected by climate change and, in turn, influences the trend of climate change due to its large population and the change in land cover.
erefore, studying the temporal and spatial changes in temperature in the Beijing-Tianjin-Hebei region can improve our understanding of the relationship between climate change and the local ecological environment.

Observation Data.
We obtained long-term  air temperature records from meteorological observation stations within and outside of the Beijing-Tianjin-Hebei region from the National Meteorological Information Center.Based on the HASM calculation, the additional data around the study area can significantly improve the HASM result [27].After screening, 64 stations with continuous observations met the research requirements during the study period (Figure 1(b)).Among them, 26 observation stations are located in the Beijing-Tianjin-Hebei region and 38 are distributed around the region.e air temperature records are processed into monthly means.

WRF.
As a next-generation mesoscale numerical weather prediction system, the WRF model has been widely applied in climate research [28,29]  information on high-resolution weather phenomena and are competitive options for use in WRF by reproducing precipitating convection and associated meteorological phenomena such as air temperature [31].We set the WRF model to output and process the hourly simulations into monthly values consistent with the processed observation data.

HASM.
e geometric properties of a surface, such as the lengths of curves, the angles of tangent vectors, the areas of regions, and the geodesics on surfaces, are called the first fundamental coefficients of the surface.
e local deformations of the surface, such as the deviations in the relevant points on the surface from the tangent plane, are called the second fundamental coefficients of the surface.
e first and second fundamental coefficients are also called intrinsic and extrinsic invariants, respectively, because the former is usually measured from outside of the surface and the latter is measured on the surface [32].According to the fundamental theorem of surface, the first and second fundamental coefficients uniquely identify a surface.HASM is developed based on this theory.
A surface, z, is defined as z � f(x, y). e first and second fundamental coefficients of the surface z are defined as follows: where equations ( 1)-( 3) and ( 4)-(6) determine the first and second fundamental coefficients, respectively.f x denotes the first-order partial derivative of the surface z with respect to the independent variable x. f xx denotes the second-order partial derivative of the surface z with respect to the independent variable x. f y , f yy , and f xy are defined similarly.Following the Gauss-Codazzi equations, the first fundamental coefficients (equations ( 1)-( 3)) and the second fundamental coefficients (equations ( 4)-( 6)) can be related as follows: where In the definition of HASM, equations ( 7)-( 9) are used to build the computational formula as follows: where A, B, and C represent equations ( 7)-( 9) in the coefficient matrix, respectively.d, q, and h represent equations ( 7)-( 9) in the right-hand vector, respectively, and n represents the number of iterations.S • z (n−1) � k is used to introduce the observation information into the surface, which is determined by the first equation of equation set (11). S denotes the sampling matrix, and k denotes the sampling vector.HASM has been applied in multiple research fields on global and regional scales [33][34][35][36][37][38][39][40][41].In the most recent research, the driving fields of HASM are created by statistical methods that usually lack the dynamic interpretation of the distribution of quantities and result in deviations in nonobservation points.erefore, the WRF model is applied in this paper to compensate for deficiency in the statistical method by dynamically downscaling the coarse-resolution CESM1 results.

Precision Test.
e observations from 26 stations located in the Beijing-Tianjin-Hebei region are chosen as the true values to test the accuracy of the temperature fields from the WRF simulation and the updated temperature fields after HASM fusion.Because the observations enter the HASM calculation as the accuracy control conditions, it is not appropriate to directly use the temperature observation data to verify the accuracy of the updated temperature field.

4
Advances in Meteorology erefore, we use the leave-one-out cross-validation method to test the accuracy of the HASM calculation results.at is, one point is extracted from 26 stations as the test point, and the other 25 points together with the 38 points outside the Beijing-Tianjin-Hebei region are input into the HASM calculation as the accuracy control conditions.e fusion accuracy of the HASM is verified by comparing the observations with the HASM output at the test point.e leave-one-out cross-validation is repeated 26 times until each of the 26 observation points serves as the test point once.
According to the meteorological definition, the 12 months of the year are divided into four seasons: spring (March to May), summer (June to August), autumn (September to November), and winter (December to February).Figures 2(a)-2(h) show the distribution of the bias in the seasonal mean temperature in the WRF simulations and HASM fusion results compared to the observations in 1956-2005.e WRF simulations in spring, autumn, and winter are significantly higher than the temperature observations, and the deviations along the west side of the study area are higher than those along the east side.In summer, the WRF model shows significantly better simulation results, with only a warm bias in the middle of the study area (higher than 0.5 °C) and a cool bias in the northeast (lower than −0.5 °C), which is accompanied by a deviation between 0.5 °C and −0.5 °C in other areas.e cross-validation shows that the accuracy of the HASM fusion results is significantly improved compared with that of the WRF model; not only is the number of high-deviation sites significantly reduced, but the degree of deviation is greatly decreased.In spring, autumn, and winter, the warm deviations in the WRF simulation in the western part of the study area are eliminated.Five stations in Tianjin and south of Tianjin improve in spring and autumn, whereas in winter, two of the five stations transform into cool biases.At the same time, during these three seasons, HASM has limited optimization ability in the northern part of the study area.
e point on the northernmost side of the study area changed from a weak deviation in the WRF model to a cool deviation in the HASM (lower than −1 °C).e reason may be that the estimation of nonobservation points in the HASM requires information from the surrounding observation points.However, in this case, the observation stations around this northernmost point are sparse, and there is no observation point on its north side, resulting in an increase in the deviation of the point.Due to the high simulation accuracy of the WRF model in summer, the HASM output does not change substantially.e distribution of biases in the annual mean temperature is shown in Figures 2(i e WRF simulation value is significantly warmer than the observation value, and the degree of deviation gradually increases from southeast to northwest.HASM significantly improves the WRF simulation, and the deviations at most of the observation points are reduced, especially the positive deviations in the western part of the study area.Similarly, the two points located on the northernmost side and along the coast of Tianjin have an increased bias due to the lack of surrounding observations, which illustrates that a boundary problem exists in the HASM method. To further quantify the accuracy, Tables 1-3 list the mean absolute error (MAE), root-mean-square error (RMSE), and spatial correlation coefficient (SCC) for the WRF-simulated temperature and HASM fusion results.e formulations of MAE, RMSE, and SCC are as follows: Tables 1 and 2 show that, after the observation data are fused, the accuracy of the HASM output is substantially higher than that of the original WRF simulation data.e deviations in the spring, autumn, winter, and annual means are reduced by at least half, and the RMSE is basically the same.At the same time, this result also reflects the high precision of the WRF simulation in summer in the Beijing-Tianjin-Hebei region, which has also been improved slightly by HASM.e calculation of MAE and RMSE is based on the bias (equations ( 12) and ( 13)); thus, the bias distribution in Figure 2 shows the MAE and RMSE between HASM and station data at different stations.e stations with greater biases have higher MAE and RMSE.e deviations of most stations in HASM are decreasing, which significantly reduces the MAE and RMSE of HASM.However, the observation station deviation on the northernmost side is continuously increasing in HASM; thus, there is a negative contribution to MAE and RMSE at this point.
In Table 3, the SCC indicates that the temperature fields before and after integration are the best in summer and autumn and slightly worse in spring and winter.e correlation between HASM fusion data and observation data is approximately 0.01-0.03higher than that with the WRF simulation data for each season, and the average annual temperature correlation increases by 0.02.e SCC indicates the degree of coincidence between the simulations and the observations with respect to the spatial distribution characteristics.e variation is only 0.01-0.03because the WRF is already highly accurate (at least 0.93 of SCC) and can simulate the spatial distribution of surface temperature very well.Although HASM has smaller MAE and RMSE, the spatial distribution characteristics of surface temperature simulated by HASM do not change greatly.From Figure 2, we can also see that, in a few stations, the biases increase, which also limits the improvement of the SCC of HASM compared with that of WRF.
Figure 3  was not signi cantly improved compared with that from 1956 to 1980.erefore, the temperature deviation in the WRF-simulated data may not be inherited from the initial eld but, rather, results from the misinterpretation of the atmospheric mechanism by the WRF model itself.Crossvalidation shows that the HASM greatly reduces the bias, with a slight cold deviation of no more than −0.2 °C.e seasonal temperature deviation characteristics are similar to those in Figure 3.
In summary, the spatial and temporal accuracies of the HASM fusion results are signi cantly higher than those of the WRF simulation.e fusion data contain information on the temperature distribution mechanism from the WRF and correction information from the observation data, which is highly suitable for long-term sequential climate analyses.In the next section, we use the HASM fusion data to analyze the temporal and spatial variations in temperature in the Beijing-Tianjin-Hebei region from 1956 to 2005.

Trend Analysis of Temperature.
e average temperatures of the Beijing-Tianjin-Hebei region are seasonally averaged, and their interannual variations are shown in Figure 4. e average regional temperature in the four seasons uctuates and increases over time, but the heating rates are signi cantly di erent.In spring, the heating rate   was approximately 0.039 °C/yr; the regional average minimum temperature appeared at the beginning of the study period (8.2  e temperature increases slowly in summer and autumn, but these rates are also higher than the global warming rate reported in IPCC5. To further analyze the change in average temperature, we divide the study period into ve periods: T1 (1956-1965), T2 (1966-1975), T3 (1976-1985), T4 (1986-1995), and T5 (1996-2005).Figure 5 shows the mean and variance in the decadal temperature over the ve periods for each season, which characterize the temperature change between the ve periods and the uctuations in temperature during each period.In spring from T1 to T3, the average temperature increased slightly but uctuated greatly, especially in the T3 period.From T4 to T5, the temperature increased rapidly, and the variance decreased, indicating that the interannual temperature changes were relatively stable.In summer, the temperature changes between the rst four periods were less than 0.2 °C and accompanied by slightly cooler temperatures in the T1-T3 period, and the interannual uctuation in temperature was small.e main warming occurred from T4 to T5, and the average temperature increased above 0.7 °C.However, during the T5 period, the interannual temperature was relatively unstable.e change in the average temperature in autumn was similar to that in spring.e rst three periods experienced little change, and the temperature increased rapidly in T4 and T5, with a heating rate of 0.4 °C/ 10 yr.Over time, the temperature in autumn gradually became unstable and the interannual changes increased.e temperature change in winter was much more complicated.After a brief cooling from T1 to T2, the T3, T4, and T5 periods heated up rapidly.
e interannual temperature changes in each period also experienced a trend of decreasing rst and then increasing dramatically in the T5 period.
Based on the above analysis, the four seasons reveal a rapid warming process from T3 to T5, which coincides with the period of rapid economic development after China's reform and opening policy.Except for spring, the regional average temperatures in summer, autumn, and winter have the largest variations during the T5 period, indicating that the temperature becomes unstable during this period, which may be related to the frequent occurrence of extreme weather in the context of global climate change.Zhai and Pan [42] studied the trend in temperature extremes in China and reported that, since the 1990s, the average numbers of warm days and warm nights exhibited increasing trends.Wang et al. [43] further con rmed that extreme temperature (highest and lowest) in North China, including the Beijing-  8 Advances in Meteorology Tianjin-Hebei region, showed a clear upward trend, and extreme warm events were much stronger than the extreme cold events.Figure 6 shows the interannual variation in the annual mean temperature and the average temperature change over ten years.
e overall warming rate of the annual mean during the study period reaches 0.035 °C/yr, and the main warming period occurs in the last 20 years, with a heating rate greater than 0.55 °C/10 yr.However, the decadal variation in temperature gradually decreases with time, which may be because the annual scale is too large, and the inuence of extreme weather on temperature is smoothed out.

Spatial Distribution of Decadal Changes.
In this section, the spatial distribution of temperature changes over time is analyzed.Following the de nition in the previous section, the study period is still divided equally into ve decades, the di erence between two adjacent decades is calculated, and the spatial evolution of the 10-year average temperature is obtained.
Figure 7 shows the spatial distribution of the decadal temperature changes in spring and summer.e overall temperature in the region increases signi cantly from T1 to T2 and from T3 to T5 in spring and from T4 to T5 in summer, which matches the ndings reported in Figure 5 well, but there are some di erences in the spatial distribution of temperature changes between di erent time periods and seasons.In the spring from T1 to T2, the northern part of the Beijing-Tianjin-Hebei region is mainly cooled, with a cooling range from −0.2 °C to 0 °C.As the latitude decreases, the air gradually becomes warmer, and the maximum temperature increase occurs in the cities of Xingtai and Handan in the southern part of the study area, exceeding 0.6 °C.e temperature distribution in the summer from T1 to T2 is similar to that in spring, but the temperature increase in the southern part of the study area is very weak.From T2 to T3, the spatial distribution and degree of temperature change in spring and summer are basically the same, showing the characteristics of slightly cooling in the south and warming in the north.From T3 to T4, the warming area covers the entire study area.In spring, the temperature rises to 0.4 °C north and east of the study area and in Beijing and Tianjin (speci cally, more than 0.6 °C in the eastern part of Tianjin and the northeastern edge of Hebei Province).
e temperature rise in this stage in summer is less than that in spring, with a small temperature increase of 0.4 °C in the western part of the study area and Beijing.During the study period, the maximum temperature increase occurred during the T4 to T5 stage.In the spring, most areas warmed above 0.8 °C, especially on the west side of the Beijing-Tianjin-Hebei region, where there was a large temperature increase of more than 1.2 °C.e temperature change from T4 to T5 in summer was slightly less than that in spring, but the distribution was similar, and a strong warming area also appeared in the west.1956-1965 1966-1975 1976-1985 1986-1995 1996-2005 Variance Average Temperature (°C) (a) Years 1956Years -1965Years 1966Years -1975Years 1976Years -1985Years 1986Years -1995Years 1996Years -2005 Variance Average  1956-1965 1966-1975 1976-1985 1986-1995 1996-2005 Variance Average  1956-1965 1966-1975 1976-1985 1986-1995 1996-2005 Variance Average

Advances in Meteorology
Figure 8 shows the spatial distribution of temperature changes over ten years in autumn and winter.From T1 to T2, the main warming areas in autumn and winter are distributed in the northwest, and most of the other areas mainly experience cooling temperatures, which are greater in winter than in autumn.In autumn from T2 to T3, the entire study area has a temperature increase between 0 °C and 0.4 °C.From T3 to T4, in some parts of the north and south areas and in Beijing, the temperature increases above 0.4 °C.In winter, from T2 to T4, the Beijing-Tianjin-Hebei region experiences   1956-1965 1966-1975 1976-1985 1986-1995 1996-2005 Temperature (°C) 10 Advances in Meteorology intense heating, especially east of 115 °E from T3 to T4, with the increase in temperature exceeding 1.2 °C and the range of the warming zone reaching its maximum of the four seasons.e intensity is also the strongest.e four seasons from T4 to T5 show a similar spatial distribution of temperature changes, but most warming areas in autumn and winter warm to a lesser extent than those in spring and summer, and the northeast regions of the study areas become slightly cool.A small range of intense warming areas above 1.2 °C also occurs in the western part of the study area in autumn and winter.e annual temperature changes over multiple decades encompass the characteristics of the four seasons (Figure 9).From T1 to T2, the central part of the Beijing-Tianjin-Hebei region is mainly cooled, with weak warming in the south and northwest regions of the study area.Most of the areas from T2 to T3 warm up between 0.2 °C and 0.4 °C.From T3 to T4, the temperature increase is further aggravated, especially in the northwest and eastern regions of the study area and in Beijing and Tianjin, where the warming exceeds 0.6 °C.From T4 to T5, the distribution characteristics of the increasing temperature from east to west during the four seasons reappear.West of the study area between 38 °N and 40 °N, the maximum temperature increase exceeds 1.2 °C.Existing research [44,45] has shown that, from the middle of the last century to the beginning of the 21st century, there has been strong warming of 0.5 °C/10 yr in Zhangjiakou and Baoding, which are in the western part of the Beijing-Tianjin-Hebei region.is strong warming is also shown in the present study, and the location where it occurred is consistent with previous studies.In addition, we further confirm that this warming phenomenon occurred only since the late 1990s and exhibited no seasonal differences.

Conclusion
In this paper, HASM is applied as a data fusion tool to obtain high-precision simulated temperature fields in the Beijing-Tianjin-Hebei region using simulated fields from the regional numerical weather model as the approximate driving fields and observed data from the meteorological stations as the accuracy control points.e main conclusions of this paper are as follows: (i) e WRF simulation of the temperature in the study area has relatively good accuracy, especially in summer and in the eastern part of the study area.However, due to the limitations of model performance, the overall deviation in the region is large.
e spatial and temporal accuracies of the HASM-fused data are significantly higher than those of the WRF simulations.e warm bias of the WRF simulation is largely eliminated, reducing the MAEs and RMSEs by half in spring, autumn, and winter, and the SCC is increased by approximately 0.02.Increased accuracy makes the fusion data more suitable for long-term climate analysis.e increase in individual biases reflects the boundary problem of HASM.(ii) e overall warming over 50 years is higher in spring and winter (0.039 °C/yr and 0.056 °C/yr, respectively), and the period from 1976 to 2005 is the period of rapid warming in the Beijing-Tianjin-Hebei region.Spring, summer, and autumn reached the highest heating rates in the last 10 years (from 1996 to 2005), while the winter warming period mainly occurred during the 1976-1995 period.e significant fluctuations in temperature between 1996 and 2005 in summer, autumn, and winter may be related to the frequent occurrence of extreme weather events.(iii) In the first 40 years (1956-1995), the spatial distribution of temperature changes is significantly different among seasons and decades and does not reflect any notable distribution patterns.One similarity appears from 1996 to 2005, when the temperature rises dramatically in the western part of the study area in all four seasons.e annual mean temperature also exhibits this warming phenomenon from T4 to T5.
Fusing simulations from the WRF model with observations from stations increases the applicable area for HASM considering that the driving fields used to be generated through statistical approaches in previous studies of climate change.However, due to the long-term study period, the available meteorological stations are limited, which may cause uncertainty in the analysis of temperature changes and makes it difficult to study historical climate change.e results of this article reveal significant warming over the last three decades, which deserves further study.In future research, the study period could be shortened to increase the available observation sites and further improve the accuracy of the data fusion.Advances in Meteorology

Figure 1 :
Figure 1: Study region: (a) location of the Beijing-Tianjin-Hebei region in China; (b) land cover of the study region derived from the European Space Agency (ESA) global land cover dataset.e red dots represent the locations of the available observation stations during the study period (1956 to 2005).

Figure 3 :
Figure 3: Interannual variability in the average temperature bias of the HASM and WRF data at 26 stations in the Beijing-Tianjin-Hebei region in 1956-2005.

Figure 5 :
Figure 5: Decadal mean and variance of the regional average temperature for each season: (a) spring; (b) autumn; (c) summer; (d) winter.

FIGURE 6 :
FIGURE 6: (a) Interannual variation in the annual mean temperature and (b) the average temperature change over ten years.

Table 1 :
Mean absolute error (MAE) values of the WRF model and HASM.

Table 3 :
Spatial correlation coe cient (SCC) values of the WRF model and HASM.
°C in 1956), and the highest temperature appeared in 2002 (12.0 °C).e warming rate in summer decreased to 0.02 °C/yr, with a maximum temperature of 20.9 °C in 1976 and a minimum temperature of 24.2 °C in 2000.e second lowest warming rate occurred in autumn at 0.026 °C/yr.e lowest temperature in autumn was 8.1 °C in 1956, and the highest temperature was 11.6 °C in 1998.Winter exhibited the highest increasing temperature rate of the four seasons (0.056 °C/yr).e lowest temperature occurred in 1969 (−8.6 °C), and the highest temperature was −2.7 °C in 2002.In the 50-year period, winter and spring have high rates of temperature changes.