Overestimating Impacts of Urbanization on Regional Temperatures in Developing Megacity: Beijing as an Example

Land-use and land cover changes may have important local, regional, and global climatic impacts by modifying the underlying land surface conditions, which in turn inﬂuence the exchange of energy and moisture between the land surface and atmosphere. Many studies have shown that urbanization has contributed to climate warming, and the amount of warming has varied. As the capital of China and one of the world’s megacities, Beijing has experienced rapid urbanization over the past 30years. In this study, we quantitatively investigated the impacts of urbanization on regional temperatures based on observations from meteorological stations and National Centers for Environmental Prediction (NCEP) reanalysis data and overestimating of the impacts were found. Comparing the temperature trends of land-use types, forest showed stronger inhibitory eﬀects on temperature increase ( − 0.085 ° C/10a). Cropland also had a negative eﬀect on climate warming yearly and seasonally, especially in winter ( − 1.133 ° C/10a) and spring ( − 0.299 ° C/10a). Conversely, the urban area showed strong warming eﬀects (0.438 ° C/10a). The conversion of cropland to urban land appeared to show the highest warming trend (0.548 ° C/10a). However, the cooling eﬀect of forest and grassland with high vegetation coverage inhibited climatic warming attributed to rapid urbanization. In addition, planting trees or grass along roadsides and increasing green parks and green roofs can also suppress surface warming. Therefore, the actual warming eﬀects of urbanization on temperatures were overestimated in megacities or urban agglomeration regions. The results showed that the green space and landscape conﬁguration should be considered in urban planning to increase green space and reduce the inﬂuence of urban heat island eﬀect.


Introduction
Increased greenhouse gases represent one of the primary factors underlying global climate change. However, increasing evidence indicates that land-use and land cover (LULC) changes, which fundamentally alter the fluxes of solar and thermal infrared radiation, sensible and latent heat, the movement of water between the subsurface and atmosphere, and the exchange of momentum between the land surface and atmosphere, may also have important local, regional, and global climatic implications [1][2][3]. In certain cases, the climatic response to LULC changes may exceed the contribution from increasing greenhouse gases [4,5].
However, fully separating the climatic impacts of LULC changes and global warming is difficult. Climate models such as the community land model (CLM), the community climate system model (CCSM), the dynamic global vegetation model (DGVM), the Lund-Potsdam-Jena model for managed land (LPJmL), the general circulation model (GCM), and the simplified parameterizations primitive equation dynamics model (SPEEDY) were used to simulate, couple, or compare for estimating the effects of LULC on climate change [6][7][8]. e simulations were based on land cover or vegetation change experiments and comparisons of climatic parameters, such as temperature and precipitation, simulated from actual and potential land-use scenarios.
However, the results varied for the simulation schemes from different models, which may further increase the uncertainties of the study. Many other studies investigated the effects of land-use change on climate primarily via observations from meteorological stations combined with statistical data or remote sensing images [9][10][11][12]. However, field data have limited coverage, and there are uncertainties associated with observations from satellites.
Urbanization has increased the surface areas available for absorbing solar energy [13], which may induce a large warming bias in regional or local temperature [14][15][16][17][18][19][20][21]. is can affect the formation of the urban heat island (UHI) [22,23], which impacts local, regional, and even global climate [24][25][26]. Gallo et al. evaluated the diurnal temperature range (DTR) at weather observation stations with respect to the predominant LULC and found that stations associated with rural-related LULC typically displayed the greatest observed DTR, whereas stations associated with urban-related LULC displayed the lowest observed DTR. In addition, the warming trend caused by urbanization was 0.06°C per century [9]. Hansen et al. compared surface air temperature changes that occurred in the United States and globally from 1900 to 2000, and they reclassified rural, small town, and urban stations based on satellite measurements of night light intensity instead of population.
eir results indicated that the warming trend caused by urbanization was 0.15°C per century [10]. us, the results differed significantly depending on whether population data or satellite measurements of night lights were used to classify urban and rural areas. Additionally, the definition of urban and rural stations based on population differed considerably among countries, and light intensity varied significantly among regions in various stages of economic development. Kalnay et al. proposed the observation minus reanalysis (OMR) method, which used the difference between trends in the observed surface temperatures and the corresponding trends in a reconstruction of the surface temperatures based on a reanalysis of the global weather to estimate the impact of land-use changes on surface warming [24]. e OMR was developed because of the lack of surface observations of temperature, moisture, and wind over land in the National Centers for Environmental Prediction (NCEP)-National Center for Atmospheric Research (NCAR) 50-year reanalysis (NNR) dataset. However, atmospheric vertical soundings of wind and temperature strongly influence the NNR, and surface temperatures were estimated from the atmospheric values. As a result, OMR has the ability to remove the signals of large-scale temperature variations induced by greenhouse gases while retaining the local signals [24]. Based on the OMR method, Kalnay [27][28][29]. Lim et al. investigated the sensitivity of surface climate change to land types over the Northern Hemisphere and found that surface warming was higher in areas that are barren, anthropogenically developed, or covered with needle-leaf forests, whereas surface warming appeared to be suppressed in highly vegetated areas [30]. ese aforementioned studies confirmed the robustness of the OMR method for detecting nonclimatic changes at the station level and providing a quantitative estimate of additional warming trends associated with LULC changes. In China, a number of studies related to the warming effects of urbanization have been performed [31][32][33][34]. However, most of these studies were on a global, country, or regional scale, and city-scale studies were infrequent. In this study, we apply the OMR method to quantitatively investigate the impacts of urbanization on regional temperatures based on observations from meteorological stations and NCEP reanalysis data, taking Beijing as an example. e study aimed to investigate the overestimation or underestimation of the effects of urbanization on regional temperatures.

Study Area.
Beijing, both China's capital and one of the world's megacities, has experienced rapid urbanization over the past 30 years. Because of the city's expansion, the land use has changed considerably, with many croplands converted to urban areas. In 2015, Beijing had an area of 16,410 km 2 and a population of 21.70 million. Beijing has a dry, monsoon-influenced humid continental climate. e monthly daily average temperature is −3.7°C in January and 26.2°C in July. e precipitation averages approximately 570 mm annually, with approximately three-fourths of the total occurring from June to August. e primary vegetation zone in Beijing is warm temperate deciduous broad-leaved forest.

Meteorological Observations.
e meteorological observation data in China are primarily sourced from the China Meteorological Data Service Center (http://data.cma. cn/), where daily data are collected for temperature, precipitation, pressure, evaporation, wind direction, wind speed, and sunshine duration from 756 national basic meteorological stations throughout mainland China. However, the station density is sparse, and only 2 national basic stations are located in Beijing. In addition to the national basic stations, we also collected daily temperature data from 17 other local meteorological stations from rural counties around Beijing provided by the Beijing Meteorological Service. All 19 meteorological stations have formed a relatively dense weather observation network, and periodic continuous observations have been conducted over a long period ( Figure 1, Table 1). e daily temperature dataset covered the period from 1979 to 2010. Repeated verifications and quality controls such as data consistency validations, abnormal value elimination, and empty value filling were conducted by the data provider. We aggregated the daily observations into monthly, seasonal, and annual mean values, and the seasonal and annual anomalies were calculated by removing the 32-year (1979-2010) mean annual cycle.

Reanalysis Data.
e monthly reanalysis temperature data were obtained from the National Centers for 2 Advances in Meteorology

Environmental
Prediction/Department of Energy (NCEP_DOE) Atmospheric Model Intercomparison Project (AMIP)-II Reanalysis (R-2) (1.9°× 1.9°) covering a period from January 1979 to December 2010, which represented an updated and human error-corrected version of the NCEP-NCAR reanalysis (R-1) dataset [35,36]. For each meteorological station, the monthly mean R-2 temperatures were interpolated to their location (longitude and latitude) using a bilinear interpolation method. en, the seasonal and annual mean values were aggregated from the monthly value. Similar to the observational data, seasonal and annual anomalies were calculated by removing the 32-year (1979-2010) mean annual cycle. To ensure the reliability of the R-2 data, we assessed the performance of the R-2 temperatures relative to the observational data. e annual anomalies for the observations and R-2 data were compared, and a good correlation was observed for the interannual variability ( Figure 2). e maximum correlation coe cient was 0.86 at the Foyeding station, and the minimum correlation coe cient was 0.52 at the Tongzhou station ( Table 2). Data in 1990Data in , 2000Data in , and 2010. Quantitative analysis of the e ects of LULC change on temperature at a city scale has rarely been performed, and raster data with 1 km × 1 km or coarser spatial resolution were typically used.

LULC
us, capturing LULC changes at the city scale is di cult, and additional detailed vector data are required to support the research. In this study, LULC data at a 1 : 100,000 scale in 1990, 2000, and 2010 were collected from China's Land-Use/ Cover Datasets (CLUDs), which were provided by the National Resources and Environmental Scienti c Data Center (RESDC) of the Chinese Academy of Sciences (CAS).  ese data were classified as farmland, forestland, grassland, built-up land, water body, and unutilized land according to the resources and utilization property. ese data were produced by human-computer interactive interpretations of Landsat ematic Mapper (TM) digital images. For areas not covered by Landsat TM data or covered by poor quality data, supplemental data were used from the CCD multispectral data from the Huanjing-1 satellite (HJ-1). After interpretations were performed, nationwide field surveys were conducted predominantly in fall for northern China and spring for southern China. e accuracy of the 6 landuse classes was above 94.3%, which met the required user mapping accuracy at the 1 : 100,000 scale [37][38][39]. In this study, the LULC dataset in Beijing was extracted from CLUDs using ArcMap 10.3, and the first four land-use types (farmland, forestland, grassland, and urban land) would be discussed.

Station Classification Scheme.
To quantify the effects of LULC change on climate, we classified the stations associated with different LULCs according to the workflow shown in Figure 3. First, we extracted LULC data in 1990, 2000, and 2010 for buffer zones with a radius of 3 km [32] centered at each meteorological station ( Figure 4). Second, we calculated the areal ratios of LULC to the buffer zone of each station. Finally, we calculated the LULC ratio changes for 1990-2000 and 2000-2010 ( Figure 5). When the areal fraction change of any land cover for a given station was greater than 20%, it was defined as a station with dramatic LULC changes; otherwise, it was designated as a station with stable LULC types. For the stable stations, if the areal fraction of any land cover in the 3 km buffer zone exceeded 70%, that land cover type was considered to be the dominant land cover type with the greatest influence on temperature changes adjacent to the station. e positions of eight stations (Table 1) were changed, which would change the underlying surface of the meteorological stations and affect observation records. erefore, only stable stations without relocation were applied for station classification.
According to the classification scheme, six stable stations (SS) and five changed stations (CS) were selected from all the stations (Tables 3 and 4). e SS included forestland, urban land, and mixed type stations, and the LULC changes of CS were primarily converted from farmland to urban land. However, there were no stations that could represent cropland stations according to the criterion. Cropland was the second largest land-use type in Beijing; therefore, it was important to investigate its impact on regional climate change. In 1990, the dominant land-use type surrounding  the stations of Daxing and Pinggu was cropland, and the areal proportion of cropland in the 3 km buffer zones was 73% and 77%, respectively. Although much of the cropland around the two stations was changed to urban land after 1990, in order to make the best use of the available information, we selected Daxing and Pinggu as the representative stations for cropland during 1979-1990 (Table 3).

Observation Minus Reanalysis Method.
Many studies have shown that the OMR analysis method proposed by Kalnay and Cai can better reflect the impact of LULC change on climate; thus, this method was used for our study. Based on temperature data from 1979 to 2010, firstly, we calculated the annual and seasonal average temperature anomalies for the observed and reanalyzed temperature. en, we calculated the trend of the temperature anomalies for the observed and reanalyzed temperature, respectively, using the linear trend analysis method, which was written as follows: where t slope represents the temperature anomaly trends of reanalysis (t rea ) and observation (t obs ), n is the number of years, and t j is each reanalysis or observed year temperature anomaly.
Finally, we calculated the OMR temperature trends as follows: where T OMR is the temperature trends of OMR temperature; T rea is the temperature trends of reanalyzed temperature; and T obs is the temperature trends of observed temperature. e temperature trend is not only related to the type but also the area of the land cover types or changes. In this study, we did not use the traditional averaged temperature trend from all meteorological stations applied in previous studies; instead, we used the area-weighted temperature trend of the entire study area because it could better reflect regional climate changes. Formula is calculated as where T is the area-weighted temperature trend; T i is the temperature trend of a typical LULC or LULC change, such as the forest temperature trend or the temperature trend of the LULC change from cropland to urban land; A i is the area of the typical LULC or LULC change; and A is the total area of the entire region.

Temperature Changes over the Last 30 Years. From 1979
to 2010, the observed annual mean temperature fluctuated with an upward trend. From the 1980s to 1990s, the observed temperature had a larger increase of 0.69°C, whereas from the 1990s to 2000s, the observed temperature had a smaller increase of 0.12°C (Figure 6(a)). Conversely, the reanalyzed annual mean temperature had a different fluctuation compared with the observed temperature ( Figure 6(b)). e reanalyzed temperature change fluctuated upwards, with an increase of 0.62°C from the 1980s to 1990s but declined by 0.38°C from the 1990s to 2000s. Figure 7 showed the observed, reanalyzed, and OMR temperature anomaly trends for each station. For the observed temperature (Figure 7(a)), all stations presented remarkable warming trends with an average of 0.463°C/10a. Spatially, stations in the suburban area showed a larger increase in temperature, and the nearest suburban stations to the city center, such as Daxing, Tongzhou, and Shunyi exhibited the most significant warming trends, with values of 0.95°C/10a, 0.93°C/10a, and 0.85°C/10a, respectively. ese trends were higher than those of stations in the urban area, such as Shijingshan, Haidian, and Beijing. e warming trends for the Shangdianzi, Zhaitang, Xiayunling, and Tanghekou stations, which were located in rural areas, were not remarkable; Tanghekou station showed the lowest trend of 0.004°C/10a.
For the reanalyzed temperature (Figure 7(b)), all stations showed warming trends, with an average of 0.249°C/10a. irteen stations had less remarkable warming trends, with values between 0.14°C and 0.20°C/10a. e higher trends occurred at the Shunyi, Daxing, and Huairou stations, and the Shunyi station trend reached 0.59°C/10a.
For the OMR temperature (Figure 7(c)), the trends were determined by subtracting the data from Figure 7(a) and Figure 7(b). Most stations had higher observed temperature values than the reanalyzed temperatures; thus, the corresponding OMR temperatures primarily exhibited positive values. Only the observed temperature trends in the Xiayunling, Zhaitang, and Tanghekou stations were lower than the reanalyzed temperatures. erefore, the OMR temperature trends were negative, which indicated a decreasing trend.  e seasonal OMR temperature trends of forestland varied considerably. In summer and autumn, the OMR temperatures showed decreasing trends. However, in spring and winter, the OMR temperatures showed increasing trends. e decreasing trends were higher than the increasing trends; thus, the annual OMR trends of the forestland also showed decreasing trends, with a value of −0.085°C/10a. Additionally, compared with the BVs, both the seasonal and annual forestland trends were lower, which indicated that forestland inhibited a temperature increase.
For the cropland, the annual and seasonal OMR trend values were negative, which indicated that cropland had a   strong inhibiting e ect on climate warming. e annual OMR trend was −0.114°C/10a, which was the lowest among the 4 dominant land-use types shown in Table 5. In winter and spring, the inhibiting impacts on climate warming were the most remarkable, and the OMR trend values were −1.133°C/10a and −0.299°C/10a, respectively.
For the urban area, the annual OMR temperature trend was 0.438°C/10a, which was the highest among the 4 types. All OMR temperatures in the four seasons showed increasing trends, and the trend values in summer and autumn were lower than those in spring and winter. In addition, all the trend values were much higher than the BVs. As a result, the urban area showed strong warming e ects on temperature changes.
For the mixed land cover combining forestland, grassland, and cropland, the annual OMR temperature showed an increasing trend; however, the trend was low (0.025°C/10a). e seasonal OMR temperature changes were similar to that of forestland in summer and autumn, which showed decreasing OMR temperatures trends. In spring and winter,    e annual and seasonal trend values were much lower than the BVs. Hence, the mixed land cover also inhibited temperature increases.
e dominant LULC change in Beijing from 1990 to 2010 was from farmland to urban land, with 742.3 km 2 converted from 1990 to 2000, which accounted for 65.5% of the total conversion area, and 489.4 km 2 converted from 2000 to 2010, which accounted for 80.6% of the total conversion area for this period (Figure 8). erefore, five changed stations (CSs) were selected to analyze the effects of urbanization on local temperatures. Table 6 shows the OMR temperature change trends of the five CSs. e results showed that as urbanization dramatically increased over the past 20 years (1990-2010), the OMR annual temperatures for all stations showed increasing trends, and the trend values were much higher than the entire study period .
ese findings indicated that the warming effect was enhanced as the farmland was converted to urban land. From 1990 to 2010 when intense LULC changes occurred, the temperature trend during the first decade (1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000) was much higher than that during the second decade (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010), which indicated that the warming trend was slowing down. is decrease may have been due to decreases in the conversion area from cropland to urban area. We analyzed the relationship between the proportion of changed area and the annual OMR temperature trends from 1990 to 2000 and 2000 to 2010 (Figure 9). e results showed a positive correlation, with a higher conversion rate from farmland to urban area corresponding to a higher temperature trend.

Combined Effects of Land Cover Types and Urbanization on Regional Temperature Change.
e unchanged (Figure 10(a)) and changed (Figure 10(b)) land cover areas for 1990-2010 were integrated into one map (Figure 10(c)) based on the following 6 categories: unchanged forest (UF), unchanged cropland (UC), unchanged urban area (UU), unchanged water area (UW), unchanged mixed type (UM) (combination of forestland, grassland, and cropland), and urbanization from cropland (CCU). e area of UW was small and was not been considered.
e Fishnet tool in ArcMap 10.3 was used to create a net of rectangular cells with predefined size. In this study, the cell size was 6 km, and the area of each fishnet cell was close to the aforementioned buffer zone (centered at the station location with a 3 km radius) in Section 2.5. An intersecting analysis was performed between the integrated map and the fishnet, and the area percentage of UF, UC, UU, UM, and CCU in each cell was calculated to identify the predominant type for each cell for the classification.
Based on the results in Tables 5 and 6, the temperature changing trend was set for each cell in Figure 10(c), and the area-weighted OMR trend (AW-OMR) ( Table 7) of the entire Beijing area was calculated using the area-weighted method presented in equation (3). e UF covered the largest area in Beijing, but its AW-OMR was only −0.034°C/10a, which was the lowest among the 5 categories in Table 7. e UC covered 22.3% of the area in Beijing, and its AW-OMR was also negative (−0.025°C/ 10a). e AW-OMRs of UU, UM, and CCU were all positive, but the AW-OMRs of UU and UM were only 0.007°C/10a and 0.006°C/10a, respectively, which seemed to have a considerably lower impact on climate warming in Beijing. For UU, the areal ratio was the lowest; thus, the warming impact of UU was also very low, although the OMR trend of UU was the second largest. For UM, the areal ratio was the second largest, but the temperature increase was the lowest, so the UM also had little impact on regional climate warming in Beijing. e AW-OMR of CCU was the largest, nearly twice that of the UF. However, the areal ratio of CCU was only 12.2%, which was nearly one-third that of the UF and one-half that of the UC. ese results indicated that UF and UC had considerable inhibiting effects on temperature increase in Beijing, and the CCU contributed greatly to climate warming. In contrast, the warming impact of UU and UM can be ignored.

Climate Effects of Land Cover and Vegetation Coverage.
Different land cover types showed varying trends of temperature change in Beijing, which were consistent with the results of Fall et al. and Lim et al., who found that surface warming was greater for areas that are barren or anthropogenically developed, whereas surface warming appeared to be suppressed by highly vegetated areas [3,30]. Satellite observations of the GIMMS NDVI dataset at 16-day intervals from 1982 to 2006 showed that the annual and summer mean NDVI of all stations, which were extracted using the bilinear sampling method at each station, were in the following sequence: forest > mixed land cover > farmland > urban land ( Figure 11). e annual mean NDVI and NDVI trend were negatively correlated with the OMR temperature trend, and the correlation ratio was −0.42 and −0.78, respectively ( Figure 12). Higher levels of vegetation coverage correspond to stronger transpiration, and the additional radiation energy would be partitioned into latent heat, and the sensible heat would be low. erefore, as the NDVI increased, the OMR temperature showed a decreasing trend.
Based on the above findings, we further revealed that the declining trend of the average annual temperature of broadleaved deciduous forest was primarily due to the strong cooling trend in summer and autumn. e differences in the temperature trends among the four seasons were also partially related to the variations in vegetation coverage and precipitation. In summer, precipitation was higher than that in other seasons, and soil moisture would also be higher; therefore, most of the radiative energy would be absorbed at the ground surface and used for physical evaporation and transpiration of water.
us, the partitioned latent heat would be higher while the sensible heat would be lower, which could explain why the temperature trend was lower in the wet season.
Forestland and grassland showed stronger inhibitory effects on temperature increases relative to other land cover types, and the continuous implementation of ecological Advances in Meteorology restorations in Beijing may play an important role in their e ects. To improve and optimize the ecological and environmental conditions of Beijing, Tianjin, and surrounding areas; protect forest and grass vegetation; and restrain deserti cation expansion, the Chinese government launched the Beijing-Tianjin sandstorm source control project in 2001. In Beijing, 6 rural counties, including Mentougou, Huairou, Miyun, Yanqing, Pinggu, and Changping, were involved in this project. From 2001 to 2010, reforestation by planting and aerial seeding covered 1.24 × 10 5 hm 2 of mountainous area, and the sparsely forested land closed for natural regeneration totaled 1.12 × 10 5 hm 2 . Nearly 7900 people had emigrated from the mountainous region to suburban areas, thereby reducing anthropogenic activities and disturbances on these fragile ecosystems [40]. For forest and grassland, although the land cover was not changed, the vegetation coverage of these areas had been greatly increased because of ecological restoration programs. Such changes in vegetation likely led to the inhibiting e ects on increase in temperature for these areas.
Cropland had a strong inhibiting e ect on temperature increase, especially in spring and winter. In spring, precipitation was less abundant and was insu cient for crop growth, so most of the crops required irrigation, especially for winter wheat. Irrigation allowed more energy to be partitioned into latent heat than into sensible heat, which caused the cooling e ect. In winter, from the observations of the two representative stations (Pinggu and Daxing), the maximum snow depth increased during 1979-1990, which would cause an increase in surface albedo and lead to decreased net radiation and temperature.
To inhibit the warming rate, continuous implementation of ecological restoration projects had been shown to be e ective by increasing the regional vegetation coverage and forest canopy density (Figures 11 and 12). Higher vegetation coverage corresponds to greater transpiration, and more radiation energy will be partitioned into latent heat, which will suppress the surface warming [3,6,30,41,42]. Additionally, greening projects in urban areas, such as planting vegetation along roadsides, increasing green spaces, can also limit surface warming. Meanwhile, the coverage of green vegetation can protect impervious surface from solar radiation, thereby reducing the net solar radiation absorbed by the surface, and less surface heat will be transported to the atmosphere [43][44][45].

Overestimating Impacts of Urbanization on Regional
Temperatures. Land cover that changed from cropland to urban area showed the highest warming trend, which may be       Advances in Meteorology the result of a decline in vegetation coverage, with the conversion of cropland to urban area and a change in the underlying surface from soil and vegetation to an impervious surface. Impervious surfaces are generally composed of bricks, tiles, concrete, and asphalt. Although these materials have high solar radiation re ectance, they also have more rapid heat transfer and a higher heat-trapping capacity relative to soil [43,44]. Solar radiation is directly absorbed by the surface, and much of the surface heat is transported to the atmosphere in the form of sensible heat, which leads to a high rate of increasing temperature. e combined climate e ects of land cover and changes measured 0.02°C/10a in Beijing, which indicated that the temperature showed increase trend under rapid urbanization but was far less than that estimated in previous studies [15-17, 31, 32, 45, 46]. Yang [32]. ere may be two reasons for the di erences between our results and those of the previous studies. China has experienced rapid urbanization and economic growth since the Chinese economic reform. From 1978 to 2000, the number of small towns soared from 2,176 to 20,312, nearly double that of the world average during this period. e number of cities increased from 190 to 663. e environment surrounding the meteorological stations in China had changed greatly due to urban expansions. Sun et al. found that after the 1980s, 109 national meteorological stations formerly located in rural areas were encroached upon by cities, and an additional 61 stations were near the cities [47]. Most previous studies [17,31,32,[47][48][49] only selected the national meteorological stations for analysis, so the results were greatly in uenced by urbanization, and the e ects on increased temperature may had been overestimated. In this study, more nonurban stations such as Zhaitang, Foyeding, Xiayuling, Tanghekou, and Shangdianzi were included in this analysis, which can more objectively explain the regional impacts of di erent LULC types on temperature change. Additionally, previous studies [17,31,32,50,51] used the arithmetic mean trend of all stations to represent the regional impacts of LULC. In this study, the combined climate e ects were estimated with an area-weighted method based on the OMR trend and area proportion of individual LULC and LULC change, which more accurately reflects the regional climate change.
rough the method, the cooling effects of large areas of forestland with high vegetation cover were highlighted. e area-weighted OMR trend was only 0.02°C/10a, which indicated that the climate warming trend attributed to rapid urbanization was greatly inhibited.

Conclusion
We investigated the combined effects of land cover and urbanization on climate warming using the OMR method in Beijing. Our calculations were based on the meteorological observations, NCEP reanalysis data, and high spatial resolution LULC data. e OMR temperature trend demonstrated that forest and cropland can both inhibit temperature increases. e cooling effect of forest was more remarkable in summer, whereas the cooling effect of cropland was more remarkable in winter. Urban areas showed increasing trends in temperature change.
e conversion of cropland to urban areas resulted in the highest warming trend, and higher conversion ratios correspond to more remarkable warming. e combined regional impacts of land cover and urbanization characterized by the AW-OMR trend was only 0.020°C/10a. Our results indicate that urbanization was still the main reason for temperature increase, whereas the forest, grassland, and cropland had a remarkable cooling effect in Beijing. Planting vegetation along roadsides and increasing green parks and green roofs can also suppress surface warming. erefore, in megacities or urban agglomeration regions, the actual warming effects of urbanization on temperatures appeared to be overestimated. Our study showed that green space and landscape configuration should be incorporated into urban planning to increase green space and reduce the influence of urban heat island.

Data Availability
In this paper, we mainly used the meteorological observation data, reanalysis data, and land-use/cover datasets. e meteorological observation data used to support the findings of this study were supplied by the National Meteorological Information Center (http://data.cma.cn/) under license and so cannot be made freely available. Requests for access to these data should be made to data@cma.gov.cn. e reanalysis data used to support the findings of this study were supplied by the National Centers for Environmental Prediction/Department of Energy (NCEP_DOE) Atmospheric Model Intercomparison Project (AMIP)-II Reanalysis (R-2) (http://www.emc.ncep.noaa.gov/) under license and so cannot be made freely available. Requests for access to these data should be made to cfs@noaa.gov. e land-use/cover datasets used to support the findings of this study were supplied by the National Resources and Environmental Scientific Data Center of the Chinese Academy of Sciences (http://resdc.cn/) under license and so cannot be made freely available. Requests for access to these data should be made to xuxl@lreis.ac.cn.

Conflicts of Interest
e authors declare that there are no conflicts of interest.