Spatial Variation of Surface Energy Fluxes Due to Land Use Changes across China

We estimate the heat flux changes caused by the projected land transformation over the next 40 years across China to improve the understanding of the impacts of land dynamics on regional climate. We use the Weather Research and Forecasting (WRF) model to investigate these impacts in four representative land transformation zones, where reclamation, overgrazing, afforestation, and urbanization dominates the land use and land cover changes in each zone respectively. As indicated by the significant variance of albedo due to different land use and cover changes, different surface properties cause great spatial variance of the surface flux. From the simulation results, latent heat flux increases by 2 and 21 W/m 2 in the reclamation and afforestation regions respectively. On the contrary, overgrazing and urban expansion results in decrease of latent heat flux by 5 and 36 W/m 2 correspondingly. Urban expansion leads to an average increase of 40 W/m 2 of sensible heat flux in the future 40 years, while reclamation, afforestation, as well as overgrazing result in the decrease of sensible heat flux. Results also show that reclamation and overgrazing lead to net radiation decrease by approximately 4 and 7 W/m 2 respectively, however, afforestation and urbanization lead to net radiation increase by 6 and 3 W/m 2 respectively. The simulated impacts of projected HLCCs on surface energy fluxes will inform sustainable land management and climate change mitigation.


Introduction
Land use/cover change (LUCC) affects climate system at local, regional and global scales through various biogeochemical and biogeophysical processes [1][2][3].LUCC would alter the physical properties of land surface, such as albedo, fractional vegetation coverage, moisture content and surface roughness.These physical properties determine the absorption, emission, and exchange of energy at the Earth's surface [1,4], represented by the fluxes of energy, radiation and moisture into atmosphere, etc. [5][6][7][8][9] and consequently affect climate system [10][11][12].Furthermore, the interaction between land surface and atmosphere mainly works through the exchanges of momentum, energy and water, among which the exchange of energy plays a key role in affecting climate system.In this context, it is of great importance to investigate the effects of LUCC on surface energy fluxes.
The effects of LUCC on climatic change have raised growing concern and have been widely studied [1,[13][14][15].There have been many studies that attempted to reveal the impacts of human-induced land cover changes (HLCCs) on the regional climate in China [3,14,16,17].However, most of these studies only focused on the effects of HLCCs on temperature and precipitation [18][19][20][21]; little is known about the effects on surface energy fluxes so far, such as latent heat flux and sensible heat flux.Moreover, previous studies on China just concentrated on one typical case study areas, for example, Northeast China, Northwest China, Jiangxi Province and the Pearl River Delta, which are characterized by different types of LUCC.To date, an integrated analysis at the national level in China remained underexplored, which is needed to make comprehensive policies regarding land planning and management.
In order to address this knowledge gap, we estimate the heat flux changes caused by land transformation over the next 40 years across China to improve the understanding of the impacts of land dynamics on regional climate.In this study, we use the Weather Research and Forecasting (WRF) model to investigate these impacts in four representative land transformation zones where reclamation, overgrazing, afforestation, and urbanization dominants the land use and land cover changes separately.The remainder of the paper is structured as follows.First, we describe the study area and introduce the main datasets used in this study (Section 2).Next, we present the methodology in Section 3. We then analyze and discuss the simulation results in detail (Section 4).Finally, we highlight the major findings in Section 5.

Selection of the Study Area
We measured the regional variance of LUCC by Land Use Dynamic Degree [22].Specifically, we overlap the vector data of land use dynamics from 2000 to 2008 on the 10 km × 10 km grid data, and then we calculate and summarize the percentages of each changing type of land use within each grid.To present the dominant changing types of land use in the study areas, four zones are selected as the representatives for four types of land transformation: cultivated land reclamation, overgrazing, afforestation and urbanization.Land use changes in China have obvious regional differentiation characteristics during 2000 to 2008 (Figure 1).Specifically, the grassland was reclaimed to cultivated land in northeast plain, eastern Inner Mongolia, northwest arid areas as well as part of traditional oasis agricultural region (Zone I).Western Inner Mongolia, loess plateau area were mainly characterized by grassland degradation (Zone II).Land use changes in mountainous Southern China, especially Jiangxi province, were dominated by reforestation or expansion of grasslands due to farmland abandonment (Zone III).In the Huang-Huai-Hai Plain and the southeastern coastal areas (especially the Bohai Economic Rim and the Yangtze River Delta), there is an apparent expansion tendency of residential, industrial, and mining land uses (Zone IV).The four zones are of heterogeneous structure and several types of land use changes occurred in each zone simultaneously.However, we only presented the impacts of dominant type of land use changes on surface heat flux in each zone.In other words, we merely take account of reclamation, overgrazing, afforestation, and urbanization in zones I-IV correspondingly.The reason why we selected these four representative types of land use changes is twofold.On the one hand, these four typical types of land use changes happened widely in China.On the other hand, these land transformations have significant climatic impacts according to the previous studies.In detail, the large scale cultivated land reclamation in China started since the last century due to the continual population growth and economic development.Such rapid expansion of cultivated land suggests an intense clearance of natural vegetation, which exerts significant impacts on the climate.Grassland accounts for 40% of the total land area of China and the grassland degradation and desertification in China have been greatly intensified due to overgrazing, overexploitation of the natural resources, rapid population growth and expansion of road network in the past decades [23].The serious overgrazing led to grassland degradation, and caused severe grassland productivity degradation and increased the frequency of extreme climate events, such as droughts and fierce freeze-up, and consequently led to acute contradiction between the human and nature in this region [24,25].A number of afforestation projects have been carried out in China, such as "Grain for Green Project" and "One Big Four Small Project", which has made great contribution to the its forestry area.Afforestation is one of the most noticeable human activities and has affected the climate not only as a carbon sink but also by changing the land surface thermal properties [26].Rapid urbanization dramatically alters the land surface properties, and thus affects the surface energy fluxes, and consequently exerts great influence on the regional climate [27][28][29][30].To sum up, impacts of typical HLCCs on regional climate and environment should be explored in detail, which can shed light on sustainable land management and climate change mitigation [31].

Data Source
The underlying surface data of the four case study areas are presented in this section.The dataset of future land use and land cover from 2010 to 2100 in China was simulated with the Dynamics of Land System (DLS) model [32] and derived from the study by Yuan et al. [33].The land use/land cover dataset of 2010 and 2050 was classified according to the land use and land cover classification system of the United States Geological Survey (USGS).In this study, the 1 km × 1 km grid data was integrated into the 10 km × 10 km grid data by resampling.
In  For Zone I, cultivated land in 2010 is approximate 450,000 km 2 , which accounts for 36.7% of the total land area of this zone, with an increasing rate of 0.3% from 2010 to 2050.In Zone II, the percentage of grassland decreases from the 37% to 36%.While the forest increases by about 0.22% in Zone III, and the proportion of urban area soars from 8% to 9.5% in Zone IV.
The dataset of RCP6.0 simulated by the fifth phase of the Coupled Model Intercomparison Project (CMIP5), was used as the atmospheric forcing data for the WRF-ARW model.This dataset includes the data such as air temperature, specific humidity, sea level pressure, eastward wind, northward wind, and geopotential height of years 2010 to 2050.

General Procedure
The WRF Model is a state-of-the-art atmospheric simulation system based on the Fifth-Generation Penn State/NCAR Mesoscale Model (MM5).The Advanced Research WRF (WRF-ARW) model is applied to simulate the impacts of LUCC on the surface energy balance in this study (Figure 3).The future spatial variation of surface energy fluxes due to land use changes can be simulated based on different underlying surface data, and the impacts of land use changes on the climate can be quantified through the difference of the two simulation results between the baseline test and sensibility test as follows: where

Solar Radiation Algorithm
According to the surface energy budget equation, there is a close relationship between the surface net radiation, land surface albedo, downward shortwave radiation, downward long-wave radiation and land surface emissivity [34].Energy comes in to the system when sunlight penetrates the top of the atmosphere, and it goes out in two ways: reflection by clouds, aerosols, or the surface; and thermal radiation-heat emitted by the surface and the atmosphere, including clouds (Figure 4).The net radiation can be expressed as follows:

R n = (1-albedo) SWDOWN + GLW + LW out
(2) where R n is the net surface radiation as a function of the surface albedo, incoming short-wave flux (SWDOWN) and long wave flux (GLW) radiation, and outgoing long wave radiation (LW out ); σ is the Stefan-Boltzmann constant (σ = 5.67 × 10 −8 J s −1 m −2 K −4 ) and T is the land surface temperature.The net radiation is the balance between incoming and outgoing energy at the top of the atmosphere.It is the total energy that is available to influence the climate.This study focuses on how the underlying surface change affects the land surface albedo, downward shortwave radiation, downward long-wave radiation, and land surface emissivity to clarify the key influencing mechanism of the future land use changes on the regional surface radiation across China.Surface radiation and energy balance are coupled through R n .The principal use of this energy is in the phase of changing of water (latent heat flux), changing the temperature of the air (sensible heat flux) and subsurface (ground heat flux) [35].The surface energy budget can be expressed as follows: where R n is the net radiation (downwelling minus upwelling solar and infrared radiation); HFX is sensible heat flux; LE is latent heat flux; and GFX is heat flux into the soil.

Spatial Variation of Flux Changes
During 2010-2050, there are various changes for heat flux among these four representative zones.Figure 5 illustrates the spatial changes of latent heat flux, sensible heat flux and surface net radiation (red represents increase and blue represents decrease).Over all three maps, in the overgrazing area (Zone II), each type of heat flux shows a downward trend while in other three zones, they present upward trends.As expected, changes in surface albedo caused by urban expansion may actually be greater than that caused by cultivated land reclamation, overgrazing and afforestation.The change of energy balance is mainly about decrease of net shortwave radiation and surface sensible heat flux due to increase of surface albedo.As a consequence of widespread afforestation, plant transpiration and surface evapotranspiration will increase and further latent heat flux will increase correspondingly.
It is found that latent heat flux increases about 2 W/m 2 in the Northeast China (Zone I) due to conversion from other land to cultivated land, while in the overgrazing area it mostly decreases by about 5 W/m 2 (Zone II) (Figure 5).In the afforestation area (Zone III), the Southeast China, the increase in latent heat flux will be mostly less than about 21 W/m 2 .Meanwhile, latent heat flux decreases by mostly less than about 36 W/m 2 in the surrounding urban expansion area (Zone IV) due to conversion from other land to built-up area.During 2010-2050, sensible heat flux will decrease by approximate 4, 2 and 13 W/m 2 in the reclamation area (Zone I), the overgrazing area (Zone II) and afforestation area (Zone III).However, in urban expansion area (Zone IV), sensible heat flux will increase by about 40 W/m 2 and the increase of albedo is associated with sensible heat flux increasing.For Zone I and II, the net radiation will decrease by 4 and 7 W/m 2 respectively.On the contrast, in Zone III and IV, it will increase by 6 and 3 W/m 2 respectively during 2010-2050.

Annual Average Heat Flux
Since each case zone has different types of land use/cover changes, we simulate the changes of heat flux under different types of LUCC in the four case zones.We compare and analyze the differences of heat flux simulations in specified grids with and without changes of land use types.Specifically, in the cultivated land reclamation area, we figure out the grids with the conversions from grassland and forest area to cultivated land, and in the overgrazing area, the specified grids are mainly with the changes to cultivated land, forest area and built-up area at the expense of grassland.In afforestation area, we focus on the land transformation to forest area (such as from cultivated land, grassland and built-up area), while in the urbanization area, the statistics are mainly about the area with the cultivated land, grassland and forest converted to built-up area.Figure 6 shows the effect of LUCC on the heat flux from 2010 to 2050, including latent heat flux and sensible heat flux and net radiation.From Figure 6, we find that the net radiation is the highest heat flux while the sensible one is the lowest.In Zone IV, the urban expansion area, latent heat flux with urban expansion and sensible heat without urbanization stay at the bottom.In the former two areas, the differences between the two types simulation are less than 10 W/m 2 for latent heat flux, sensible heat flux and even net radiation.Different from latent heat flux in Zone I, with the situation of land reclamation, the net radiation and sensible heat flux will both be lower than that without it.On the contrary, for the Zone III, latent heat flux has less discrepancy, which is less than 10 W/m 2 , but difference of sensible heat flux between the area with and without afforestation is more than 30 W/m 2 , and as for the net radiation, the difference is about 20 W/m 2 .Under the urban expansion, sensible heat flux increases by about 40 W/m 2 while latent heat flux decreases by about 40 W/m 2 .

Monthly Average Heat Flux
In order to figure out the monthly changes between 2010 and 2050, we count and analyze the three kinds of heat flux in the four case zones on a monthly basis (Table 1).In the cultivated land reclamation area and overgrazing area, the changes are small with less than 10 W/m 2 .Latent heat flux in Zone I will increase by about 1.8 W/m 2 on average and in June it will increase by 3.2 W/m 2 , which is maximum through the whole year.However, sensible heat flux and net radiation will decrease by 3.8 and 4.0 W/m 2 on average, respectively, and the extreme maximum values will appear in August and October with the minimum ones appearing in June and July, respectively.In the Zone II, both types of heat flux will decrease and the smallest change is the change in sensible heat flux.Because of vegetation changes in Zone III and IV, changes in the heat flux present distinctly.Since the forest area increases, latent heat flux and net radiation increase by about 20 and 4.7 W/m 2 respectively, while sensible one will reduce by about 13.6 W/m 2 .In contrast, in the Zone IV due to the green reducing, latent heat flux will decrease by 36 W/m 2 with the maximum decrease appearing in September and minimum appearing in May while sensible one increase by 39 W/m 2 , with minimum and maximum increase of 38 and 41 W/m 2 appearing in July and May respectively), and the net radiation only changes no more than 1.5 W/m 2 during this period.

Conclusions
In this study, we investigate the effects of HLCCs on regional climate in typical reclamation, overgrazing, afforestation, and urbanization regions in China with the aid of WRF model.Simulation results show that future HLCCs patterns have significant impacts on the surface energy fluxes.From the simulation results, there is significant variance of albedo under different HLCCs, which suggests that albedo of different underlying surfaces results in spatial variance of the radiation flux.Results show that latent heat flux increase by 2 and 21 W/m 2 in the case of cultivated land reclamation and afforestation respectively.On the contrary, overgrazing and urban expansion leads to latent heat flux decrease by 5 and 36 W/m 2 correspondingly.Simulation results also indicate that urban expansion has significant effect on sensible heat flux, and it leads sensible heat flux to have an average increase of 40 W/m 2 during 2010-2050.Reclamation, afforestation and overgrazing have opposite trends, i.e., decrease sensible heat flux.In addition, net radiation would decrease by approximate 4 and 7 W/m 2 in reclamation and overgrazing dominated region respectively; and increase by 6 and 3 W/m 2 in afforestation and urbanization regions separately.
This study also reveals the seasonal pattern of the impacts of HLCCs on surface heat flux.The considerable changes of latent and sensible heat flux and net radiation appear in winter and spring, especially in the afforestation and urbanization areas.Furthermore, different HLCCs have various effects on monthly variation of latent and sensible heat flux and net radiation.The change tendency of sensible heat flux is consistent with net radiation under various underlying surface types.Urban expansion, afforestation, cultivated land reclamation and overgrazing have different impacts on latent heat flux and the effect of urban expansion is by far the most significant of all.The simulated impacts of projected HLCCs on surface energy fluxes in a relatively long run will shed light on screening climate change mitigation strategies through regulating land use activities.

Figure 1 .
Figure 1.Main land use and land cover change in China, 2000-2008.

Figure 2 ,
four bars show the dynamic area of the dominant land use and cover changes in these four study areas from 2010 to 2050, and the percentages attached to the bars indicate the proportions of the representative land use types to the total area of this region in 2010 and 2050 correspondingly.

Figure 2 .
Figure 2. Area and percentage of dominant land use and land cover in each zone.
refers to the effects of future land use change on the surface energy fluxes; i indexes the annual or monthly average surface energy fluxes; and P and B denotes to the simulation results obtained with the WRF-ARW model based on underlying surface in 2010 and 2050 respectively in this study.

Figure 3 .
Figure 3. General procedure followed in the study.

Figure 4 .
Figure 4. Diagram of the energy balance.

Figure 5 .
Figure 5. Changes in latent heat flux (a), sensible heat flux (b) and surface net radiation (c) in the four case zones during 2010-2050 (Unit: W/m 2 ).