Projected Impacts of Bioenergy-Demand-Induced Land Use and Cover Changes on Regional Climate in Central Europe

1 Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, No. 11A, Datun Road, Anwai, Beijing 100101, China 2 Center for Chinese Agricultural Policy, Chinese Academy of Sciences, Beijing 100101, China 3 Faculty of Information Engineering, China University of Geosciences, No. 388, Lumo Road, Wuhan 430074, China 4 School of Economics, Renmin University of China, Beijing 100872, China 5 State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, Beijing 100048, China 6China Institute of Water Resources and Hydropower Research, Beijing 100048, China


Introduction
Energy shortfalls have been increasingly serious all over the world; meanwhile, the renewable energy is rapidly developing but has not become a significant source of energy yet [1].Since the 1990s, the United States, high-energy-consumption countries in Europe, and the agricultural products trading countries such as Brazil have begun to develop the renewable bioenergy [2].As the main type of the renewables, the worldwide governments are also promoting its development in order to mitigate the climatic effects of the consumption of massive fossil fuel.At the same time, it is becoming a scientific research hotspot in recent decades.Some previous researches on the effectiveness of bioenergy have estimated the potential impacts of the development of renewable energy according to the change of greenhouse gas (GHG) emission due to land use/land cover change (LUCC) and by means of life cycle analysis (LCA) [3].LUCC has great impacts on GHGs, and it also modifies the surface energy and water balance [4,5] through influencing the near-surface temperature and precipitation.Besides, it serves as an additional driving force of the climate change at the global [6][7][8][9] and regional [10][11][12] scales.This study mainly focused on the effects of bioenergy development on LUCC, which influences the biogeophysical processes of the land surface and subsequently impacts the climate change.
The bioenergy has showed high potential in coping with the worldwide energy crisis and increasingly serious environmental problems since it is famous for being renewable, biodegradable, nontoxic, and environmentally friendly [13].Most researchers have studied the impacts of the growth Advances in Meteorology and production of the bioenergy with the GBEP common methodological framework for GHG lifecycle analysis of bioenergy [14].An additional 15% of greenhouse emission results from the land use change, especially the change linked to deforestation caused by the expansion of agricultural land.Croplands have higher albedo than that of forests, and there is lower absorption of solar radiation over deforested areas [15].Due to bioenergy production needing cultivated land or land [16,17], it changes land use structure [18], affecting the carbon sinks at the same time and then influences climate change, including the change of temperature, precipitation, and evaporation.What is more is that the Global Bioenergy Partnership [19] has reported that specific significant methodological uncertainties and methodological choices that can significantly affect the indicator values are related to indirect land use change and base year to measure the LUCC.
Nevertheless, most of studies indicate there may be some (potential) surplus of agricultural land in the EU in the near future, while the cropland in Central Europe has shown a downward trend since 1970 [20].The size of this surplus, however, depends on a large number of factors.The general conclusion is that the future availability of land for bioenergy production will be largely determined by the agricultural policy and market development.Another general conclusion is that the above-mentioned studies focused on the aggregated amount, which only provide limited insight into developments at the regional scale.A number of studies provide important information on the environmental effects of bioenergy production at the plot scale.However, these studies failed to capture the effects of bioenergy production at the landscape scale [21].
The climatic impacts of LUCC through influencing the terrestrial biogeophysical processes have attracted more and more attention in the current scientific domain [22][23][24].LUCC can alter the convective rainfall and other regional weather patterns such as the surface flux of moisture and sensible and latent heat fluxes [25].It is very difficult to summarize the effects of land cover change on the climate since the climatic impacts of different biogeophysical effects may offset each other [26].Besides, in comparison to the impacts of land cover change at the global or annual scales, the impacts of land cover change at the regional scale have not been well represented in the annual or global average statistics since it may often show opposite signs [27].In fact, the coarse resolution of the Global Circulation models (GCMs) seriously limits their capability to capture the mesoscale features that play a pivotal role in regional dynamics [12].By comparison, the Regional Circulation models (RCMs) are more suitable for the studies at the regional scale.
Currently, an increasing number of researchers have studied the effects of LUCC on the regional climate with the dynamic downscaling [20,28,29], with different original data, conversion technologies, and applications.The expansion of agricultural land at the middle latitudes leads to the decrease of the range of daily average temperature in many areas [30].Many scholars studied effects of LUCC on the regional climate with the RCMs inversion, while the parameterization of the process is not completely targeted on the bioenergy in some researches.This study mainly aims to clarify how the production of bioenergy influences LUCC and its impacts on the climate change.This study firstly reviewed how the development of bioenergy in Central Europe influences land use such as conversion from grasslands and forest lands to cultivated lands.Then the impacts of LUCC on the temperature, precipitation, and heat balance were detected and analyzed on the basis of simulation with the Weather Research and Forecasting (WRF) model.

Methods
This study explored the LUCC under the influence of demand of bioenergy production for land resources with a multiscale and multimodel approach and simulated its effects on the climate change.Firstly, the land use structure from 2010 to 2050 was simulated with the Agriculture and Land Use model (AgLU) in GCAM.The AgLU model is a top-down economic model with enough structures to simulate the global land use change and the consequent carbon emission over one century [31].These simulations were implemented with and without a carbon policy which was represented with a positive carbon price.Changes in the carbon price create incentives for commercial production of biomass, which affects the distribution of other land types and subsequently influences the carbon emission from land use change.Commercial production of biomass provides a link between the agricultural and energy systems.In this study, the main driving forces of land use change were analyzed, including the economic growth, income growth, demand for food growth, population growth, climate change, and especially demand for bioenergy production.The core of the LUCC lies in the distribution mechanism between the cropland, grassland, forest land, and unused land and the maximization of regional economic profit.
Demand equation (1) for each of the crops, processed crops, and animal products is expressed with the unit of total calories demanded per year as follows: CDR ,, =  , ⋅ (  ,,  ,,=baseyear ) where CDR ,, is the demand for calories (total calories demanded per year in a region at time ),  can be crops, processed crops, or animal products,  is a region index,  is the point in time of the calculation,  , is a scaling coefficient (unitless) to calibrate the price and income feedback terms in the base year,  ,, is the market price of the commodity ($/calories consumed per person per day),   is the price elasticity,  ,, is per capita income index {$ GDP/population/($ GDP/population in base year)},   is the income elasticity for the calories demanded which implies it may differ from income elasticity in the demand equation for transportation, industry, or buildings (  ),  , is the total population by region (which differs from the population index, the ratio of the regional population at time  divided by the regional population in the base year, in the demand equation for energy services in the transportation and building sector), and  ,, is calories consumed per person per day (multiplied by 365 to get yearly consumption).
Equation ( 2) is to compute demand for two types of forest products, industrial wood and fuel wood: where WDR ,, is industrial wood and fuel wood in cubic meters demanded.Selection of land use is based on maximizing economic return at each location.Profit per hectare is equal to revenue (yield per hectare times price received) minus production cost (yield per hectare times nonland cost per unit of output).This relationship is shown in (3), where  is an index for land use type, l is the region index, and  is an index for geographical location within a region: where  ,,, is the economic return of the land as a profit rate ($/ha-yr),  ,, is the market price for the product produced by land use  (units $/yield units: calories or m 3 ),  ,, is the nonland cost per unit of output in land use (units are $/yield units: calories or m 3 ).P,, is the forward price for forest products which is the time lag between planting and harvest.Consider where  is the interest rate ($/$, that is, unitless).With specific assumptions on the functional form of the yield distribution, the share of land allocated to use  is given by a logit share equation: where  is a positive parameter that determines the rate that land shares change in response to a change in profit rate.Land use for a specific purpose is calculated based on this logit-based share of total land: Then the Dynamics of Land System (DLS) model was used to spatially simulate the LUCC through identifying conversions between cropland and other land cover types during the period of 2010-2050.DLS, a good model for allocation of land use pattern depending on the specific demand, is a collection of programs that simulates the changes in the land use pattern through conducting scenario analysis of land use change in the study area [32,33].Firstly, the effects of geophysical conditions and socio-economic environment on the spatial patterns of cropland area were explored, and the predominant driving factors on the spatial pattern of cropland were identified.Secondly, the changing trend of the predominant factors influencing spatial patterns of cropland was predicted through considering the historical characteristics and current status at the regional scale.Thirdly, a proper scenario was selected to predict the balance between the supply and demand of land resources.The net changes of cropland area, which are constrained by the supply and demand of land among various land-consuming sectors, were allocated into each pixel according to the land use structure.Finally, the cropland area was spatially allocated at 30 km × 30 km grid scale, and the spatial pattern was finally generated.
Finally, the climate change was simulated with the Weather Research and Forecasting (WRF) model in order to assess the impacts of LUCC on the regional climate.The WRF model, a regional climate model for higher regional projection, is a next-generation mesoscale numerical weather prediction system designed to serve both atmospheric research and operational forecasting needs.It is featured by two dynamical cores, a data assimilation system and a software architecture that allows the parallel computation and system extensibility.The model serves a wide range of meteorological applications across scales ranging from meters to thousands of kilometers.The model version 3.5 which was developed by the National Center for Atmospheric Research (NCAR) includes three land surface schemes in Table 1: the simple soil thermal diffusion (STD) scheme, the Noah scheme, and the Rapid Update Cycle (RUC) scheme [34].The lateral boundary forcing data came from the NCEP/FN dataset and was updated every 6 hours.In the parameterization scheme of the physical process in the model, the cumulus parameterization scheme adopted the Grell-Devenyi ensemble scheme, the boundary layer process scheme was YSU, and the shortwave radiation scheme was the CAM scheme, while the land surface process scheme was Noah land surface model.The boundary buffer was set to be 4 layers of grid points, and the boundary conditions adopted the relaxation scheme.The time interval of the model integration was set to be 5 minutes, and that of the radiation process and cumulus convection was half an hour and 5 minutes, respectively.There were 27 layers in the vertical direction, and the atmospheric pressure at the top layer was 50 hPa.
The WRF model was used to simulate the climate change from 2040-2050 with the first year as spin up.There was a set of tests which was used to simulate the climate of 2010 to validate the method.A 10-year simulation ensured that we would see how surface climate changes with respect to land cover change, while we found that a 1-year spin up was enough to initialize the soil moisture.The simulation area was centered at 50 ∘ N and 10 ∘ E and was projected on a Lambert conformal grid covering almost the whole of Germany; the area covers 32 and 21 grid points in the longitudinal and latitudinal directions, respectively.The spatial resolution of the model is 30 km × 30 km; the main topographic features of the simulation area are captured at this fine resolution.

Data Processing
3.1.Biofuel in Europe.The Organization for Economic Cooperation and Development (OECD) gives out the statistics of the contribution of renewables to energy supply (Figure 1), which showed an obviously increasing trend during 1990-2010.
There has been substantial growth of bioenergy in the EU over the past decades, especially the bioenergy for the modern energy carriers' electricity and the transport sector.The proportions of renewable energy produced in the EU in 2005 were as follows: 66.1% from biomass, 22.2% from hydropower, 5.5% from wind power, 5.5% from geothermal energy, and 0.7% from solar power.In general, bioenergy plays an essential role in the European ambition to increase the share of renewable and domestic energy.During 1990-2010, the thermal power increased by approximately 2% per year, bioelectricity increased by approximately 9% per year, while the bioenergy production increased about eight times (with an increasing rate over 20% per year).
The contribution of bioenergy to the energy supply in the EU in 1990 was a bit less than 2000 PJ, accounting for two-thirds of the total renewable energy production or 4% of the total energy supply in the EU.There will be a significant increase of the market share of bioenergy in energy supply in the EU, with the aim of reaching 1% on average (it has doubled in two years).In addition, there is significant difference in the market share of bioenergy in energy supply among different countries in the EU; for example, that of Germany (3.8%) and Sweden (2.2%) has achieved the reference value.
In order to reach the ambitious target of a 20% share of energy from renewable sources in the overall energy consumption amount, the EU plans to make efforts on the electricity, heating and cooling sectors, and biofuels.In the transportation industry, which is almost exclusively dependent on oil, the EU commission hopes that the share of biofuels in the overall fuel consumption will have reached 10% by 2020.So it was assumed that the share of biofuels may achieve 20% by 2050 according to this trend in this study.

History of LUCC and Climate Change.
With the rapid development of bioenergy, the land use pattern in Europe has also gradually changed.The land use change caused by the development of bioenergy is mainly represented by the conversion from grassland and forest to cropland.According to OECD statistics, about 27.5% of land in Europe is arable and permanent cropland, and permanent grassland accounts for only 16.2%, while the forest accounts for 32.4%.We also found in the History Today (2001) that the cropland in Central Europe increased from 41.9 Mha in 1900 to 56.9 Mha in 1960, and since then it fell to 46.7 Mha in 2005 with the growth rates per year of 0.41% from 1900 to 1950 and −0.03% between 1950 and 2000.The annual average temperature has been 10 ∘ C, and the annual precipitation has ranged from 700 mm to 1000 mm in Central Europe during 1900 to 2010.According to the statistics of temperature and precipitation data in the past [35], the temperature showed an obvious increasing trend with the increment of 1.3 ∘ C (Figure 2(a)), and there is also an increase in the coefficient of variation, while there is no obvious changing trend of the precipitation with little decrease of 15 mm (Figure 2(b)).The climate in the study area is characterized by the peaks and valleys of temperature and precipitation appearing during the same period.The minimum values of temperature and precipitation both appear in January, reaching about 2 ∘ C and 80 mm, respectively, while their maximum values generally appear in July, reaching about 18 ∘ C and 100 mm, respectively.

Land Use Change.
In Central Europe, the main land use type is cropland.At the beginning of the 21st century, cropland accounted for about 67% of the total land area and forest accounted for 17%, while grassland only accounted for 10%.The cropland is mainly distributed in the northern part of the study area, while the forest is generally located in the southern part where there are mountainous areas.The current land use situation and the dominant conversion types are shown in Figure 3. Based on the GCAM model, the trend of land use change in the next 40 years was simulated.The result indicated that it is mainly characterized by the conversion from cropland to other land use types in the northwest part, while the conversion from forest into cropland dominates the southern part.Nearly 30% of the grassland and 2.7% of the forest land will be converted into cropland in the study area in the future.However, in the process of the development of biomass energy, not only the other land use types will be converted into cropland but also some cropland will be converted into grassland, while on the whole, the cropland area in Central Europe will increase to some extent.Part of grasslands, forests, and some other land will be converted into cropland, which will cause the change of land use and eventually influence the climate change in Europe (Table 2).The trend of land use change in the next 40 years was obtained according to the simulation result with DLS, which was presented with the red and yellow boxes as main conversion from cropland and conversion into cropland, respectively (Figure 3).The result indicated that there will be some obvious land use change in the study area.In the northwestern part, a lot of cropland will be converted into other land use types, while in the southern part, some forests and grasslands will be converted into cropland.The projected land use in 2050 will be presented in Figure 4, and obviously cropland has increased and forest has decreased.

Climate Change in 2010-2050.
To validate the precision of simulation result with the WRF model, the simulation data of year 2010 were compared with the observation.The comparison result indicated that the accuracy of the projection of temperature (Figure 5(a)) was higher than that of precipitation (Figure 5(b)).In spring and winter from January to April and November to December, simulated temperature is a little far away from the observed one, and the temperature difference of winter is higher than other seasons, and the simulation values were closer to the observations in the summer and autumn.The simulation result indicated that the simulated temperature from January to April and from November to December was a little far away from the observed one.Besides, the difference between the simulated and observed temperature in the winter is greater than that in other seasons, while the simulated and observed ones in the summer and autumn are more consistent.On the whole, the trends of the simulated and observed temperature were consistent, both of which show an increasing trend from January and February, reaching the peak in August and then show a decreasing trend.The comparison between the simulated and observed precipitation indicated that there is significant difference between simulated and observed values.According to the observation data, the precipitation is the lowest in March, and the peak occurs in June, reaching about 1600 mm.Then it declines to less than 500 mm in November.By comparison, the simulated precipitation fluctuates over the year, with higher values occurring in the summer and autumn.
As is shown in Table 3, the minimum temperature in 2050 appears in December, reaching 2.56 ∘ C, while the maximum temperature appears in July, reaching 21.73 ∘ C. The average temperature is 12.38 ∘ C, while the counterpart in 2010 is 11.19 ∘ C. The minimum and maximum of precipitation occur in November and July, reaching 80 mm and 463 mm, respectively.In terms of cumulative precipitation, there is 3346 mm in 2050 and 5749 mm in 2010.By comparison, the heat flux is much more stable all year round and so is the latent heat.During 2040-2050, the temperature shows an increasing trend, with an increment of 0.2 ∘ C, while the precipitation will decrease by 65 mm per month.The latent heat shows an increasing trend (Figure 6), with an increasing rate of 0.4% per year, while the heat flux stays at a stable status.

Conclusions.
The result of this study indicated that the increase of production of bioenergy in Central Europe would lead to more demand for cropland in the southern part of the study area.The increasing demand of bioenergy for land resources would lead to LUCC which will change the land surface and subsequently influence the future climate.The LUCC in the future was simulated with DLS, and the impacts of LUCC on the climate were analyzed on the basis of projection with the WRF model.The simulation result indicated that the region land use change is mainly characterized by the conversion from forest and grassland to farmland by 2050, which makes the latent heat flux and temperature show an increasing trend, while precipitation will show a decreasing trend.

Discussions.
This study explored the impacts of LUCC on climate change at the landscape scale under the influence of demand of bioenergy production for land resources.The results indicated that expansion of cropland is mainly due to the human demand for bioenergy, which is consistent with the current energy strategies and policies [36,37].Once the Common Agricultural Policy (CAP) is implemented, there may be some reform and rationalization in agriculture due to the change of current land use, productivity of agricultural land, and the expected major transitions of the agricultural sector in Central Europe [38,39], which may provide some opportunities for the alternative crops.Besides, the production of bioenergy may also increase the intensity of land use through improving the technology and education [40].
The spatial pattern of production of bioenergy represented with the ecological indicators of farmland quality [41] will be essential to the assessment of the climate effects of the development of bioenergy as well.Conversion from forest to cropland or pasture reduces the aerodynamic roughness of the landscape and decreases both the capture of precipitation on the canopy and the root extraction of soil moisture; these changes tend to decrease evaporation and hence reduce the fluxes of moisture and latent heat from the surface to the atmosphere, which acts to increase the temperature near the surface.In addition, the land use change is related to not only human activities but also the natural processes, which have not been taken into consideration in this study.What is more is that the earth is a complex system that involves various factors that directly influence the LUCC and subsequently influence the climate change as well as various factors that directly influence the climate change.Therefore, there are still some uncertainties in the result of this study, and it is necessary to carry out further researches on the sensitivity and uncertainty of the impacts of LUCC on the climate change.

Table 1 :
Physical process and parameterization scheme.

Table 2 :
Transition matrix of LUCC with grid count and percent in Central Europe during 2010 to 2050.

Table 3 :
Monthly projected values of climatic factors in Central Europe in 2050.