Close Home
iForest - Biogeosciences and Forestry
vol. 7, pp. 372-384
Copyright © 2014 by the Italian Society of Silviculture and Forest Ecology
doi: 10.3832/ifor1170-007

Collection: RegioResources21
“Spatial information and participation of socio-ecological systems: experiences, tools and lessons learned for land-use planning”
Guest Editors: Daniele La Rosa, Carsten Lorz, Hannes Jochen König, Christine Fürst

Research Articles

Future land use and food security scenarios for the Guyuan district of remote western China

Lin Zhen (1)Corresponding author, Xiangzheng Deng (2), Yunjie Wei (1), Qunou Jiang (3), Yingzhi Lin (2), Katharina Helming (4), Chao Wang (1), Hannes Jochen König (4), Jie Hu (1)

Introduction 

“Land use” is understood to mean human activities that exhibit a spatial dimension and that change the bio-geophysical conditions of land and the associated environment ([16]). Land use is a key factor that determines a landscape’s performance with respect to socioeconomic functions such as land-based production, infrastructure, and housing ([26]). Significant changes in land use and land management have taken place in China since 1978, and one of the most important drivers of these changes has been human intervention through the implementation of land-use policies. One particularly dramatic example is the Sloping Land Conversion Program (SLCP), the initial phase of China’s Grain for Green program, which was initiated in 1999 by the State Council ([33]) and which was fully implemented in 2000. Excessive reclamation of land for agriculture in environmentally fragile frontier regions had led to a series of natural disasters ([23]). For instance, severe flooding in areas surrounding the Yangtze River in 1998 caused an economic loss of respectively 13.5 billion RMB (renminbi) in Shandong Province and 156.8 billion RMB (US$ 19.6 billion) overall, and at least 1526 people died ([50], [22], [18]). In response, the Chinese government realized that it had become necessary to take immediate action to alleviate these water-related problems and the related soil erosion that occurred in agricultural land on steep slopes ([44]). Based on their experience with the pilot program, the government recognized that far more marginal agricultural land should be retired from agricultural use and revegetated using either trees or grassland vegetation ([11]). Therefore, the original program evolved into a more widespread program to combat soil erosion, alleviate the deterioration of natural ecosystems, safeguard water resources, and promote sustainable development in rural China ([35]).

The goal of the SLCP is to convert agricultural land on steep slopes or heavily degraded land into forest or grassland ([31], [14]), and to restore degraded forest and grassland; the ultimate goal is to prevent soil erosion by reducing the area of exposed soil. The steepness of the slope (15° or above) is one of the main criteria for the selection of farmlands to be converted into forest or grassland. The SLCP stipulates that farmers who convert degraded farmland on steep slopes into either grassland or “ecological forests” (orchards or plantations of trees with medical value, and timber-producing forests that will not be harvested because of their role in protecting slopes) will be compensated with: (1) an annual in-kind subsidy of grain to compensate them for their lost crops; (2) a cash subsidy; and (3) free seedlings and grass seed for use in revegetation at the beginning of the planting period. To account for differences in regional average yields, the grain subsidy was set at 2250 kg ha-1 in the Yangtze River Basin and 1500 kg ha-1 in the Yellow River Basin. The annual cash subsidy is 300 RMB (about US$ 44) per hectare of eligible land to pay for tending and miscellaneous expenses on the newly converted land. The grain and cash subsidies last for 8 years if ecological forests are planted and for 5 years or 2 years if economic forests or grasses are planted, respectively ([45]). In addition, the SLCP provides farmers with a one-time cash subsidy of 750 RMB ha-1 (about US$ 110 ha-1) to purchase seeds or tree seedlings for planting.

The policy has been welcomed by farmers because the cash and food compensation generally exceeded their net production, and in some cases even exceeded their gross yields on the targeted land ([4]). It is possible, however, that as arable land is converted into forest and grassland and built-up area, the loss of agricultural land will conflict with China’s national food security goals. Although little was heard about this conflict in the late 1990s through 2003, rising grain prices in early 2004 led policymakers and scholars to debate the role of these land conversions in rising food prices, and to question whether the land conversion might threaten China’s food supply ([27], [14]). Among other responses in early 2004, the State Council addressed the importance of slowing down the conversion of natural vegetation to arable land. Some researchers believed that this policy reversal was unnecessary and could slow regional economic growth ([17]). Another group of scholars claimed that the change was critical to maintain national food security ([2], [39], [46]).

In China, food security is looked at from both quantitative and qualitative aspects, from quantity aspect it includes three perspectives: per capita food or grains; national total quantity; and the concept of regional quantity. From quality aspect, food security implies the nutrition level that can be met from protein foods consumption. A multitude of studies have addressed the concerns raised by changes in land use. Most were based on historical observations (e.g., [30], [23], [36], [48], [6]), on investigations of specific land-use change processes and their driving forces in local case studies (e.g., [15], [47], [55], [53]), or on the loss of agricultural land and its causes ([30], [3], [31]). A study of land use dynamics based on land-use inventories in Italy revealed that land conversion into forest over time was similar in magnitude to land conversion into urban areas (e.g., [7]). In some recent papers, König et al. ([19], [20], [21]) assessed the impact of the SLCP on land-use functions in the Guyuan District of the Ningxia Hui Autonomous Region of north-central China and Inner Mongolia, and found that the program had improved the environment by mitigating soil erosion, but that the region had little economic potential in non-agricultural sectors, and the bottom-up stakeholder participation is essential for ensuring successful implementation of ecosystem conservation program. In particular, the forestry sector provided only limited job opportunities to local people. Zhen et al. ([52]) investigated the changing food consumption patterns and impact on land requirement as a result of SLCP. The SLCP requires integrated forest management to minimize the risk of unsustainable water consumption by agriculture while contributing economic and social benefits ([19]). There has been little or no systematic investigation of future land-use patterns and the associated impacts on food security if the SLCP continues.

Land conversion under the SLCP is determined by two primary variables: the slope of the land (i.e., its eligibility for the program) and the subsidy provided to the farmers (i.e., the incentive to participate). In the present study, our goal was to analyze land-use changes and the effects on grain production under the baseline (current) scenario and under a conservation policy scenario, each with three growth-rate scenarios for each driving factor, and to analyze the impact of the resulting land-use change on China’s food security.

The paper starts with a brief introduction to the study site, the Guyuan District of north-central China. We then describe our research methodology, including the parameters and variables used for scenario development, and provide a detailed description of the baseline and policy scenarios. We then explain the dynamics of land systems (DLS) modeling method used to predict future conditions in each scenario. We conclude by presenting the major results of our analysis and discussing their implications.

Study area (Guyuan District) 

Guyuan District is located in a remote, environmentally fragile mountainous frontier region of north-central China (Fig. 1). The area has a mean annual rainfall of 483 mm and an annual potential evapotranspiration of 1250 to 2000 mm, with only rain-fed agriculture possible in most of the district; however, irrigated agriculture is possible in small areas. Farmers in the area have limited resources and mostly practice subsistence agriculture. Their per capita holdings of arable land, forest, and grassland are 0.27, 0.28, and 0.31 ha, respectively, which are higher than the corresponding national averages for farmers of 0.11, 0.11, and 0.24 ha ([29]). Major crops grown in the area are cereal crops such as winter wheat (Triticum aestivum), summer maize (Zea mays), and potato (Solanum tuberosum). The remaining land is used mainly for the cultivation of millet (Panicum miliaceum), buckwheat (Fagopyrum esculentum), oilseed plants (Brassica napus), oats (Avena sativa), beans (including field beans, Vicia faba; common beans, Phaseolus vulgaris; and peas, Pisum sativum), sweet potato (Ipomoea batatas), leafy green vegetables (e.g., lettuce, celery), fruit trees, and fodder crops. As of 2009, the population density was 142 persons km-2, and out of the total of 1.5 million people, some 89% of the population was employed as farmers.

Fig. 1 - Location of the case study area, Guyuan District of the Ningxia Hui Autonomous Region.

Methods 

Dependent variables and factors influencing land use in Guyuan District

The original land-use map of Guyuan district was derived from Landsat TM and ETM images of the middle 2000s including the images for 2004/2005 with the spatial resolution of 30 × 30 m. The images were processed by geometric correction, image registration, etc. using the software of ENVI ([8]). Land use in the study area were extracted by means of a supervised classification using the software ArcGIS® (ESRI, Redlands, CA, USA) based on differences in the spectral, textural, and color characteristics of the image, which were expressed by the land use code (1 for arable land, 2 for forest, 3 for grassland, 4 for water area, 5 for built-up area and 6 for unused land). The results of the image interpretation were then verified through field surveys; more than 200 plots in the study area were chosen to test the interpretation accuracy, proved through the kappa value. The overall precision of all kinds of land uses was more than 92%. After the land use information was extracted, it was converted to the grid data with the spatial resolution of 300 × 300 m. When the ground-truth observation revealed that a cell contained multiple land uses, we used the land use that covered the majority of the cell area as the land-use type in that cell. Land uses were then reclassified and expressed using a special representation scheme in which each 300 × 300 m cell of the grid was expressed as a binary value (1 or 0) to indicate the existence or absence (respectively) of each type of land use in that cell ([8]).

The dynamics of a land-use system are potentially influenced by many factors. In this case study, we divided these factors into four categories: geophysical, climatic, proximity, and socioeconomic. These factors are generally accepted as determinants of the land-use distribution and as the driving factors responsible for land-use change ([37]). The land-use types and the influencing factors used for our analysis of the land-use change scenarios are listed in Tab. 1. All of these factors were used for simulation of our scenarios. The relationships between the binary variable for each land use and the factors influencing these land uses were identified, measured, and then considered in our simulation of land-use changes. We explored the relationships between the land uses and the factors influencing them by means of spatial regression analysis.

Tab. 1 - List of land-use types and influencing factors considered in the Guyuan District case study.

First of all, all of the factors were collected and processed at the pixel level with a spatial resolution of 300 × 300 m. Among them, the climatic data acquired from the China Meteorological Bureau were interpolated by means of kriging ([8]). Geophysical data were provided by the Data Center of the Chinese Academy of Sciences, which were all pre-made in the 300 × 300 m grid level. Proximity data was estimated using national highway and provincial capital data with the software ArcGIS®. The spatial distribution of population density was simulated at the grid level by combining remotely sensed data, meteorological data, soil data, and statistical data by means of a grid-generation model ([24]). Finally, all of the factors and land use data were inputted into the software GeoDa ([8]) to carry out the spatial regression analysis.

Scenario development

In the present study, we developed scenarios for land-use changes based on a literature review and expert judgments that relied on broad consultation with experts from the Chinese research community (e.g., the Chinese Academy of Sciences, Chinese Academy of Social Sciences, China Agricultural University, Peking University, Bureau of Survey and Mapping of China), with policymakers (e.g., the Administration for China’s Agenda 21, Ministry of Agriculture, Ministry of Environmental Protection, State Bureau of Forestry), and with local decision makers (e.g., from the Commission of Development and Reform, Bureau of Agriculture, Bureau of Forestry, Bureau of Environmental Protection, Bureau of Statistics).

A list of drivers for the baseline scenarios and variables for environmental conservation policy scenario was prepared based on our knowledge and experience about the local situations. Relevant officials at the aforementioned organizations were then requested to rank the most important drivers and variables. Altogether 20 experts were involved in the interview, and a final list with main drivers for the baseline scenario and variables for the policy scenario was obtained (refer to Tab. 2 and Tab. 3). Based on these inputs, we developed three levels of growth in the driving factors, to simulate land-use changes in Guyuan District from 2005 to 2020.

Tab. 2 - Baseline scenario and variables.
Tab. 3 - Key variables used to develop the alternative environmental conservation policy scenarios and resulting land-use changes in Guyuan District from 2010 to 2020. (†): Net present value derived from the subsidy in 2005, considering the estimated inflation rate ([8]), multiple-cropping index, and governmental investment in land-use conservation in Guyuan District.

Baseline scenario and growth scenarios

In the baseline scenario, we assumed that the current subsidy for participants in the SLCP will not change, and that no new areas of sloping land will be targeted for conversion under the SLCP. We identified four main drivers for the baseline scenario: population density (persons km-2), urbanization trend (% of urban population), government’s investment in agricultural research and development (R&D, in RMB), and gross domestic product (GDP) per unit land area (10 000 RMB km-2). The high-, moderate-, and low-growth rate values for each of those four drivers were estimated based on their past trends and expert judgment from the interview (Tab. 2).

Environmental conservation policy scenarios and growth scenarios

Policy interventions to change land use are determined through adjustments of the values of various policy variables. In the present study, the environmental conservation policy variables were the subsidy paid to farmers for each unit of converted land (RMB ha-1) and the proportion of the total land area (%) with a slope > 15°. For the scenario analysis, we set the values of these two policy variables as high, moderate, and low in order to examine the resulting land-use changes (Tab. 3). The variables that will affect the land-use changes were defined subjectively according to the projected values of these variables and the plans of the local government. In the conservation policy scenarios, land-use changes and driving factors from the baseline scenarios were used as starting point, and then the effects of the subsidy and area of sloping land were added. The specific values of each policy variable (Tab. 3) used for the simulation were determined as follows:

  • Subsidy: Under the moderate growth scenario, we accounted for the estimated inflation rate ([8]), the multiple-cropping index (i.e., the area-weighted proportion of each crop in the study area), and the potential increase in land conservation ([29]), and defined total subsidies of 2589 and 3273 RMB ha-1 in 2010 and 2020, respectively. The subsidies under the high- and low-growth scenarios were assumed to be 20% higher and lower, respectively, than the value in the moderate-growth scenario. This percentage was chosen based on local land conversion plans ([5]).
  • Proportion of the land area with a slope > 15°: We used the value for the base year (2005) for the moderate growth scenario in 2010, and for 2020, assumed a total increase of 0.24% based on the current annual plan for conversion arable land on slopes into grassland and forest ([5]). The values for the high- and low-growth scenarios were set to 5% higher and lower, respectively, than the level in the moderate-growth scenario based on local land conversion plans ([5]).

Modeling framework

We adopted the DLS modeling framework developed by Deng et al. ([8]) to simulate land-use changes in our scenarios. Based on the conversion of land-use and its effects (CLUE) model developed by Verburg et al. ([40]), DLS embeds two special features: scenario analysis for land demand under open system (i.e., an initial set of parameter values) and a spatial econometric regression between the explained variables and the explanatory variables. DLS regards the land-use changes occurring under land demand driven by economic and social developments to be an open system with a regional extent. Pixel-specific land-use changes are influenced by various demographic, societal, economic, location, and geophysical factors. We conducted equilibrium analyses of the land supply and demand by sectors and by sub-regions under an open system ([8]), and used the results to simulate the land-use changes at a regional extent.

Given the complexity of land-use systems, the DLS model incorporates parameters accounting for technological changes, lifestyle changes (e.g., food consumption patterns), economic growth (e.g., GDP), population growth, and urbanization. Scenarios for land-use change were developed based on these factors, and the changes in land use were spatially disaggregated into individual cells of the grid based on the simulated relationships between land uses and their influencing factors. There are three main modules in DLS: a spatial regression module that identifies the relationships between land uses and their influencing factors; a scenario analysis module that identifies the land-use changes required by the regional demand for each land use; and a spatial disaggregation module that allocates the overall regional land-use changes to the disaggregated grid cells.

Population projections and the grain-balance analysis

We predicted the future population of Guyuan District to let us compare the grain supply (based on the predicted land-use changes) with the grain demand (based on the predicted population changes) in 2005, 2010, and 2020. We used the following equation for the population projection (eqn. 1):

\begin{equation} P = P_{0} - A e^{-y} \end{equation}

where P is the total population in the projection year, P0 is the population in the initial year of the simulation (2005), y represents the exact year number for simulation (i.e., the number of years between the initial year and the projection year), and A is the annual population growth rate. In this study, we defined population growth rates of 0.58, 0.40, and 0.17 % for the high-, moderate-, and low-growth scenarios.

We calculated grain equivalents based on the total grain production and the total grain demand of a given population level. Grains included winter wheat, spring-sown maize, millet, buckwheat, and potatoes. Other crops which did not account for a significant proportion of the total harvest were omitted.

Results 

Analysis of driving factors for the dynamics of land-use change

Tab. 4 summarizes the regression coefficients for the factors influencing the land uses of arable land, forest, grassland, water area, built-up areas, and unused land in Guyuan District in 2005. The standard errors in parentheses indicate that the influencing factors vary in their effect on the land use dynamics. The smaller the standard error, the stronger the parameter’s effect on the land-use dynamics. For instance, the cumulative temperatures above 0 and 10 °C, elevation, population density, and per capita income strongly affect the change in the areas of arable land, forest, and grassland, and the sunshine duration and phosphorus concentration in the soil strongly affect the area of arable land.

Tab. 4 - Estimation of the pixel-specific regression model with spatial lag terms for explained variables. Standard errors are included in parentheses; the total number of observations was 124 927. Constant terms were estimated but not reported here. (†): W_y represents the spatial lag terms for the explained variables for cultivated land, forest, grassland, water area, built-up area, and unused land.

Dynamic simulation of land-use changes

Among the six land-use categories defined by the Chinese Academy of Sciences ([24]), the dominant land uses of Guyuan in 2005 were grassland (46.8%) and arable land (43.2%), followed by forest (7.4%), built-up land (1.7%), and bodies of water and other unused land (<1%). Unused land is land which cannot be used directly, including saline-alkali soil land, marshland, sandy soil land, and bare exposed rock or gravel (Fig. 2).

Fig. 2 - Spatial distribution of the main land uses in 2005, based on supervised classification of Landsat TM data.

Results from the scenario analysis show that the balance among the major land-use types in Guyuan District did not change greatly under any of the scenarios, with arable land and grassland remaining the dominant land uses (accounting for more than 80% of the total land area), followed by forest, built-up land, bodies of water, and unused land (Tab. 5). In comparison with the baseline scenario, clear differences in the trends for the high-, moderate-, and low-growth scenarios were obtained (Tab. 5, Fig. 3 and Fig. 4).

Tab. 5 - Land-use changes (in terms of area × 103 ha) under the three baseline scenarios and the three environmental conservation scenarios from 2005 to 2020. Unused land is not listed due to its small and relatively stable land area (ca. 100 ha) in all scenarios. (†): Numbers in parentheses represent the increase (+) and decrease (-) in the area or percentage of each of land-use category compared with the value in the base year. (‡): Percentages of the total land area.
Fig. 3 - Simulated land-use pattern in Guyuan District under the baseline scenarios.
Fig. 4 - Simulated land-use pattern in Guyuan District under the environmental conservation scenarios.

Changes in land use under the influence of the conservation policy scenarios were characterized by a steady decrease in the area of arable land, with the proportion of the total area accounted for by arable land in 2010 decreasing (as compared with the value in the base year) by 3.1, 3.6, and 4.2 % in the low-, moderate-, and high-growth scenarios, respectively; the corresponding decreases by 2020 were 6.0, 6.9, and 8.2%, respectively. About half of the decrease in the area of arable land is expected to originate from an increase in built-up areas, whereas other losses will occur as a result of conversion to horticultural uses (e.g., parks and gardens), and rehabilitation of land through reforestation or grassland restoration (Fig. 5 - [14], [39], [48]). The decrease in the area of arable land is clearly greatest in the high-growth scenario. The relatively high subsidy for the conversion of arable land into grassland and forest, as well as restrictions on crop cultivation on slopes > 15°, will decrease the area of arable land compared with the baseline scenario. Most of the converted arable land is distributed near built-up land and sloping land (Fig. 2-4). Compared with the baseline scenario, the decrease in arable land is greater under the conservation policy scenarios. For example, in the high-growth scenario, 4.2 % (2010) and 8.2 % (2020) of the area of arable land will be lost under the environmental conservation policy, versus 1.6 and 3.0 %, respectively, under the baseline scenario.

Fig. 5 - Land use transitions in Guyuan under baseline and environmental conservation scenarios.

The area of grassland decreased at the same rates in the low- and moderate-growth conservation scenarios in 2010 (0.3%) and 2020 (0.5%), and increased in the high-growth scenario, with a slightly higher rate in 2020 (0.2%) than in 2010 (0.1%). This suggests that the increase in grassland strongly depends on government incentives (high subsidy) and regulations (which mandate the conversion of farmland on steep slopes into grassland and forest). Existing large areas of grassland also slow the rate of further expansion. Decreases in the area of grassland resulted mainly from its conversion into arable land and forest, whereas increases resulted from conversion of forest and arable land into grassland ([48], [54]). However, in the baseline scenario, the area of grassland continues to decrease under all three growth rates.

Forest and built-up areas showed a steady expansion. The area of forest increased by 2.5, 3.1, and 3.6 % by 2010 under the low-, moderate-, and high-growth scenarios, respectively, with corresponding values of 4.9, 6.1, and 6.9 % in 2020. Built-up areas increased by 0.7, 0.6, and 0.5 % in 2010 in the low-, moderate-, and high-growth scenarios, respectively, higher than national average of 0.25 % from 2002 -2008, due to the fact that the initial value of built-up areas was relatively smaller in less developed western China like Guyuan than rest of the country; therefore, the net area increase would be higher compared to the country average ([42]). The built-up areas will increase by 1.5, 1.2, and 0.9 % under respective scenarios by 2020. The area of forest expanded mostly in the central mountains of Guyuan District, and resulted primarily from the conversion of nearby grassland (Fig. 3 and Fig. 4 - [54]). The continuous increase in the area of forest resulted from the governmental policy, in that the benefits to be obtained from forested land (e.g., forest products) created from arable land will be owned by the farmers ([33]). This expansion of forest is higher in the conservation policy scenarios than in the baseline scenarios (Tab. 5).

Urban and rural settlements and other built-up areas would expand dramatically under the baseline scenario, and the newly expanded urban land would mainly appear near existing urban areas or residential centers. This increase would be the fastest in the low-growth scenario, followed by the moderate- and high-growth scenarios, in both 2010 and 2020. This would result from the strong enforcement of the policy to increase forest and grassland areas. Even so, no large-scale, interconnected urban areas will develop in Guyuan District by the end of 2020, and the increase in the built-up area under the conservation policy will be less than that under the baseline scenario (Tab. 5).

The change in the area of bodies of water in Guyuan District, which is affected by annual variations in precipitation, would be relatively small compared with the changes in other land uses (i.e., 0.8 % of total land use). Therefore, the area and spatial distribution of bodies of water would remain largely unchanged during the simulation period in both the baseline and conservation policy scenarios. Unused land accounts for a very small proportion of the total land area (0.04%), and its area will increase slightly (e.g., to about 0.10% in all scenarios) or will remain unchanged. Biophysical constraints such as low precipitation will limit changes in water bodies, and the poor quality of the soils in unused land (which has undergone severe soil erosion, or which has a high abundance of rocks and stones) limit the potential other uses of this land.

Impacts of land-use changes on grain production and food security

We analyzed changes in food security based on projections of regional agricultural production and population growth. Based on the land-use changes under the baseline and conservation policy scenarios, we estimated grain balances (surpluses or deficits). Tab. 6 summarizes the scenario for the area cultivated, the total supply (production) and demand for grains, and the resulting net grain balance. Compared with the base year (2005), the arable area in the baseline scenario increased by 3.2, 3.4, and 3.5 % in 2010 under the low-, moderate-, and high-growth scenarios, respectively, versus corresponding increases of 3.6, 4.3, and 5.0 % in 2020. Therefore, if current land-use drivers continue to act at the predicted levels in the future, the arable area will continue to increase. However, the arable area in the conservation policy scenario in 2010 will decrease by 0.5, 1.0, and 1.5 % under the low-, moderate-, and high-growth scenarios, respectively, versus corresponding decreases of 1.5, 3.0, and 4.4 % by 2020. The decrease in the arable area under the conservation policy scenario will result from continued implementation of the SLCP.

Tab. 6 - Grain production and demand in Guyuan District from 2005 to 2020 under the baseline (BL) and conservation policy (CP) scenarios. Numbers in parentheses represent the percentage increase (+) or decrease (-) in comparison with the base year (2005).

Grain demand will continue to increase between 2005 and 2020, primarily because of population growth, but also because of our assumption that consumption will increase as prosperity increases due to economic growth. The grain balance, which represents the gap between total production and total demand, will decrease continuously until 2020 in all scenarios and for all rates of growth. There will be a surplus of grain in all scenarios in 2010, but the surplus will be lower under the conservation policy scenarios. By 2020, the grain balance will become negative in all scenarios (i.e., a deficit), indicating that regional grain production can no longer meet the increasing demand from the region’s population. However, the deficit will increase with increasing growth rate, and will be larger under the conservation policy scenarios than in the baseline scenarios. Therefore, implementation of the environmental conservation policy, coupled with a decrease in the area of arable land, population growth, and increasing prosperity, will lead to food shortages that must be met by imports from other regions. The impacts of improved cultivation technology on the grain yields has been identified and represented in scenario analyses based on historical trends and predicted future technological development ([5]).

Total grain production decreased in the conservation policy scenarios compared to the baseline scenarios. For example, in 2010, the decrease was 3.7, 4.4, and 5.1 %, respectively, in the low-, moderate-, and high-growth scenarios, versus corresponding values of 5.5, 7.9 and 10.1 %, respectively, in 2020, suggesting a greater decrease in production in the high-growth scenario. Therefore, the decrease in the area of arable land under the conservation policy scenario will reduce the area available for grain production and increase the risk of grain shortages compared with the baseline scenario. For example, in the simulation year of 2010, the total grain loss in the high-growth conservation scenario, in comparison with the low-growth scenario, is 5400 t (with a loss of arable land of 1.5 and 0.5 % respectively in comparison with the base year); in 2020, the total grain loss increases to 16 800 t in the high-growth conservation scenario compared to the low-growth scenario, with the percentage decrease in the area of arable land of 4.4 and 1.5 %, respectively (Tab. 6). We calculated the land area required to meet the demand for each crop based on the predicted productivity per unit area for that crop (Tab. 7). Minor shifts in land use can be seen for each major crop. The land area of the major crops decreased by 0.8, 1.6, and 3.3 %, respectively, in the low-, moderate-, and high-growth conservation scenarios, which suggests that the decreased area of arable land will affect each of the major crops under cultivation.

Tab. 7 - Land demands for the production of major crops by local farmers. (†): Percentage decrease of land area from 2005 to 2020.

Discussion 

Land-use dynamics

Between 40 and 90% of the decrease in the area of arable land in Guyuan District is estimated to have resulted from the SLCP ([6]). However, the history of land-use change in Guyuan District shows that the area of arable land increased by 8.3% (39 111 ha) from 1986 to 2000, and 97.9% of the new arable land was converted from grassland and the rest was converted from former bodies of water (mostly where the water had dried up or was drained) or forest ([6], [51]). This increase mainly resulted from the high proportion of total GDP accounted for by agriculture (more than 50%), the high proportion of the rural population (more than 85%), and low agricultural productivity during this period ([6]). The high reliance of the local economy on agricultural production motivated farmers to reallocate land from grassland and forest to farmland to improve the grain supply. Our study revealed that, although the arable land area is decreasing, only minor decreases in the area of agricultural land are expected based on simulations (Tab. 7). There are two main reasons for this: the regulations for basic farmland protection place tight restrictions on land conversions from protected agricultural land to non-agricultural uses; and, because the study area is located in a dry and remote mountainous region, all of the farmer households rely heavily on farming for their livelihood, so they protect their agricultural land carefully. The production of this land is maintained by means of high inputs (e.g., fertilizer) rather than by expansion of the arable area.

Though differences between the analyzed scenarios seem minor, they are significant enough to provide guidance for policy development. We found that the difference in the area of arable land between the high- and low-growth environmental protection scenarios would be nearly 1.1×104 ha in 2020, which accounts for more than 3.0% of the total arable area in Guyuan District. The Chinese government has defined 1.2×108 ha as the minimum area (the “red warning line”) for arable land in China to ensure food self-sufficiency. The implication of this standard for each of China’s districts and counties is that the area of arable land reserved for grain cultivation must be strictly protected.

Grain balance

Grain yield in Guyuan District has increased at an annual rate of 9% from 1971 to 2005, mainly based on improvements in production technology and increased production inputs ([6]). As a result, the district has maintained self-sufficiency in food production for the past 30 years. Since the current land-use conversion process is expected to last for an additional 10 to 30 years ([14]), the impacts on the total grain balance should generally be small, but they will nonetheless be significant at the local level. For example, our simulations revealed a severe grain deficiency by the end of the simulation period (2020) under both the baseline and the conservation policy scenarios. Some authors have concluded that the ongoing conversion of arable areas into other uses under the SLCP will significantly reduce total production ([12], [34]), but that agricultural productivity will increase in north-central China as a result of more intensive agricultural practices and the development of improved cultivation technology for the remaining land ([14]). However, our study considers the predicted average grain yield in the scenario analysis, and the result shows that even though grain yield (including winter wheat, spring-sown maize, millet, buckwheat, and potatoes) is increasing, the grain demand is increasing as a result of population growth and grain production is simultaneously being undermined by the loss of farmland for ecological conservation. Accordingly, the grain losses in our simulations lead to serious concerns about regional food security and may lead to a need to import food from outside the Guyuan District.

Land-use change and its impact on food production have been studied previously. For instance, Belsky & Siebert ([1]) observed results similar to ours in Indonesia, and concluded that the conversion of swidden fields into cacao plantations, which produce a cash crop that is vulnerable to global market trends, constrains future food production opportunities and raises significant concerns over local food security and agricultural sustainability. Case studies on land conversion and food security in western China have produced a range of conclusions. For example, Feng et al. ([13]) found that if land conversion continues as planned, a total of 400×104 t of grains will be lost, but that this will not have negative impacts on food security because production by the remaining farmland can be enhanced by increased agricultural inputs. In addition, the improvement of ecological conditions brought about by the increased forest and grassland areas will decrease the frequency of natural disasters and help to secure grain yields, despite reductions in the area of arable land and an increasing population. Lu ([25]) and Zhang & Wu ([49]) argued that under current land conversion plans, food shortages may occur only in areas where food demand and supply are currently closely balanced or where food security risks already exist. Dong et al. ([9]) compared food security with and without land conversion, and found that in Guyuan District the grain supply was more secure for farmers who participated in the SLCP than those who did not, because farmers were compensated for their converted land in the form of grains and cash. However, they also found that local production should be improved to ensure long-term food security. Xi & Cao ([43]) studied land use in China’s loess plateau and found that at a lower per capita grain consumption level (e.g., 190 kg year-1), the current rate of arable land conversion would not jeopardize food security by 2010, but if per capita food demand increases to 230 kg year-1, food shortages would occur. Therefore, both food demand and production should be considered when tackling food security issues.

Our scenario analysis quantified the consequences of policies that directly or indirectly affect land-use change decisions and agricultural production by accounting for land-use change and changes in grain demand from the increasing population. It therefore provides valuable information to facilitate both systematic management of the existing land resources and planning and policy development for future arable land conversion directions and target quantities of each land use. Moreover, it reveals the trade-offs between conservation policies and food security in a sensitive region of China.

Methodology

The DLS approach provides an effective framework for simulating the dynamics of a land-use system. However, it should be noticed that the simulation results are not meant to predict the actual future land use; instead, they are projections designed to explore the impacts of assumed scenario conditions by allowing researchers to quantitatively visualize the results of scenario descriptions ([41]). Analysis of the driving factors responsible for land-use change is a significant part of the DLS approach, in which the influencing factors are aggregated and processed at a pixel level, and then estimated by means of spatial regression analysis. In the present study, all variables were resampled to a pixel size of 300 × 300 m before performing the regression analysis. This may be problematic because some variables (e.g., climatic and socioeconomic variables) have a very different original resolution. Climatic and socioeconomic variables are derived from the station or counties, so there are some errors when allocated into the grid pixel level. However, such error cannot be avoided if the regression analysis is carried out at the pixel level.

The validity of a model is mainly determined by the specific characteristics of the case study and by the quality of the input data ([41]). Using 2010 data from Chinese statistical yearbooks ([28]), we compared actual grain production in Guyuan District with the values predicted in our baseline scenarios. We found that our estimate in the moderate-growth baseline scenario was within 5% of the actual values. This suggests that our modeling approach is reasonable and should be refined in future research to provide even better estimates. However, the uncertainty resulting from reducing the complexity inherent in land-use systems could be reduced by supplying a more flexible interface allowing the users to input spatial regression results by means of more robust model specification and by developing more than one kind of scenario to simulate the dynamics of land-use systems under various conditions.

Visualization and decision support

The visualization of land use trajectories for different scenarios can support policy discussions about the future land-use patterns in Guyuan District specifically, and for all of China more generally, by supporting the identification of priority areas for intervention and by testing the potential consequences of certain policy options. Our comparison of the results for the baseline and conservation policy scenarios in the present study revealed that the area of arable land will decrease and the areas of forest and grassland will increase. This is likely to become an important issue in terms of the region’s food security. Many studies in western China indicated that such land-use changes are common ([23], [14], [38], [48]). From the results of this study, we can provide the following advices to the local decision makers to ensure food security while conserving the natural environment:

  1. Population control: population natural growth rate in the study area is 1.18%, which is higher than the national average (0.479% - [32]), it is thus important to continue population control policy.
  2. Strict restrictions on occupation of arable land by other uses, in order to protect basic farmland areas.
  3. Increase government investment in agricultural R&D for increasing total production through improved crop varieties and cultivation techniques.
  4. Strengthening the monitoring and management of grassland and forestland converted from arable land, so that the natural environment could be conserved and farmers income from forest products can be increased.

Conclusions 

Between 2010 and 2020, the land-use dynamics in Guyuan District will be characterized by a decrease in the area devoted to agriculture and an increase in forest and grassland cover, accompanied by urban expansion. The increased areas of forest and grassland are distributed mostly in the central mountains of the district, and most of the new urban development is concentrated in downtown and residential areas, where most of the infrastructure and economic opportunities are located. The area of agricultural land has fluctuated over the past 20 years, but has only decreased regularly since implementation of the SLCP in 2000 in Guyuan district, after which a 9.3% reduction in the area of agricultural land was observed from 2000 to 2005. For the other major land use types from 1986 to 2000, the areas of forest, urban, and unused land increased by 8.0, 8.7, and 6.0 %, respectively, but the area of grassland decreased by 8.1%. As the urban population increased from 6.6% of the total in 1989 to 12.3% in 2007 ([10], [29]), the area of urban land also increased. The implementation of the SLCP since 2000 has been the major process that has permitted the expansion of the areas of forest, grassland, and urban land.

Our study considered the predicted grain yield (including winter wheat, spring-sown maize, millet, buckwheat, and potatoes) in the scenario analysis, revealing a severe grain deficiency by the end of the simulation period (2020) under both the baseline scenarios and the conservation policy scenarios; such effect will be more severe under the high-growth scenarios and in the conservation scenarios. Therefore, food security is at risk in Guyuan District. However, decision makers could take some strategies to protect the region’s food security, the main recommendations include control of population growth, increase investment in agricultural R&D, improvement in monitoring and management of grass and forest, and strict protection of basic farmland.

Our scenario analysis quantified the consequences of policies that directly or indirectly affect land-use change decisions and agricultural production. Thus, they provide valuable information to support systematic management of the existing land resources and planning and policy development to guide future arable land conversion directions and the target quantities of each primary land use, and reveal the trade-offs between conservation policies and food security in a sensitive region of China.

Acknowledgments 

This research was funded by the “Key Project for the Strategic Science Plan in IGSNRR of CAS (No. 2012ZD007), Technical Support Program of the Ministry of Science and Technology of China (No. 2013BAC 03B04), and the Ten Year’s Eco-environment Assessment Project of MOST (STSN-14-00, National Key Project for Basic Research no. 2014CB95430201). We extend our heartfelt thanks to Geoff Hart for his assistance with the English language and for other valuable comments.

References

(1)
Belsky JM, Siebert SF (2003). Cultivating cacao: Implications of sun-grown cacao on local food security and environmental sustainability. Agriculture and Human Values 20: 277-285.
::CrossRef::Google Scholar::
(2)
Brown LR (1995). Who will feed China? Wake-up call for a small planet. WW Norton& Company, New York, USA, pp. 163.
::Google Scholar::
(3)
Caradec Y, Lucas S, Vidal C (1999). Agricultural landscapes: over half of Europe’s territory maintained by farmers. European Commission’s Directorate General for Agriculture and Rural Development, web site.
::Online::Google Scholar::
(4)
CCICED (2004). China’s agricultural and rural development in the new era: challenges, opportunities and policy options’ policy recommendations and executive summary. China Council for International Cooperation on Environment and Development, Center for Chinese Agricultural Policy, Chinese Academy of Sciences, Beijing, China, pp. 23. [in Chinese with English summary]
::Online::Google Scholar::
(5)
CDRG (2007). 11-12th five-year development plan of Guyuan district. Commission of Development and Reform of Guyuan, Ningxia Hui Autonomous Region, Guyuan, China, pp. 35. [in Chinese with English summary]
::Google Scholar::
(6)
Chen CC (2008). Study on the ecosystem carrying capacity in Jinghe River Basin. PhD Dissertation, Graduate University of Chinese Academy of Sciences, Beijing, China, pp. 184. [in Chinese with English summary]
::Google Scholar::
(7)
Corona P, Barbati A, Tomao A, Bertani R, Valentini R, Marchetti M, Fattorini L, Perugini L (2012). Land use inventory as framework for environmental accounting: an application in Italy. iForest 5: 204-209.
::CrossRef::Google Scholar::
(8)
Deng XZ, Su HB, Zhan JY (2008). Integration of multiple data sources to simulate the dynamics of land systems. Sensors 8: 620-634.
::CrossRef::Google Scholar::
(9)
Dong M, Zhong FN, Wang GJ (2005). An empirical study on grain for grass and food security in less developed regions: a case study of Ningxia Hui Autonomous region. China Population, Resources and Environment 15 (1): 104-108. [in Chinese with English summary]
::Google Scholar::
(10)
DPS (1991). China population statistics yearbook 1990. Department of Population Statistics, National Bureau of Statistics of China, Scientific and Technical Documents Publishing House, Beijing, China, pp. 771. [in Chinese with English summary]
::Google Scholar::
(11)
Du SF (2001). Environmental economics. Encyclopedia Press, Beijing, China, pp. 102. [in Chinese with English summary]
::Google Scholar::
(12)
Feng ZM, Li XL (2000). The stratagem of cultivated land and food supplies security: storing food in land - raising the comprehensive productivity of land resource of China. Geography and Territorial Research 16 (3): 1-5. [in Chinese with English summary]
::Google Scholar::
(13)
Feng ZM, Zhang PT, Song Y (2002). Food security: the impact of land conversion from farmland to forest or grassland on grain production in Northwest China. Journal of Natural Resources, 17 (3): 299-306. [in Chinese with English summary]
::Google Scholar::
(14)
Feng ZM, Yang YZ, Zhang YQ, Zhang PT, Li YQ (2005). Grain-for-green policy and its impacts on grain supply in West China. Land Use Policy 22 (4): 301-312.
::CrossRef::Google Scholar::
(15)
Hao SL, Cheng BC, Zhao XM, Yu Q (2005). Dynamic analysis on the spatial pattern of landscape of small watershed in loess hilly-gully region: a case study of Shanghuang test area in Guyuan. Science of Soil and Water Conservation 3 (2): 107-112. [in Chinese with English summary]
::Google Scholar::
(16)
Helming K, Tscherning K, Koenig B, Sieber S, Wiggering H, Kuhlman T, Wascher D, Perez-Soba M, Smeets P, Tabbush P, Dilly O, Huttl R, Bach H (2008). Ex ante impact assessment of land use changes in European regions - the SENSOR approach. In: “Sustainability Impact Assessment of Land Use Changes” (Helming K, Perez-Soba M, Tabbush P eds). Springer-Verlag, Berlin, Germany, pp. 77-105.
::Google Scholar::
(17)
Huang JK, Deng XZ, Rozelle S (2004). Cultivated land conversion and bioproductivity in China. In: “Remote sensing and modeling of ecosystems for sustainability” (Gao W, Shaw R eds). Proceedings of the SPIE, Bellingham, Washington, DC, USA, pp. 135-148.
::Google Scholar::
(18)
Jiang T, Shi YF (2003). Global climatic warming, the Yangtze floods and potential loss. Advances in Earth Sciences 18 (2): 277-284. [in Chinese with English summary]
::Google Scholar::
(19)
König HJ, Zhen L, Helming K (2012). Assessing the impact of the sloping land conversion program on rural sustainability in Guyuan, western China. Land Degradation & Development (early view).
::CrossRef::Google Scholar::
(20)
König HJ, Uthes S, Schuler J, Zhen L, Purushothaman S, Suarma U, Sghaier M, Makokha S, Helming K, Sieber S, Chen L, Brouwer F, Morris J, Wiggering H (2013). Regional impact assessment of land use scenarios in developing countries using the FoPIA approach: findings from five case studies. Journal of Environmental Management 127: 556-564.
::CrossRef::Google Scholar::
(21)
König HJ, Podhora A, Helming K, Zhen L, Wang Z, Wübbeke J, Baumeister T, Du B, Yan H (2014). Confronting international research topics with local stakeholder perspectives on multifunctional land use: the case of Inner Mongolia, China. iForest (early view)
::CrossRef::Google Scholar::
(22)
Li AT (1999). 1998 Yangtze flood and fight against it. Yangtze River 30 (1): 1-7. [in Chinese with English summary]
::Google Scholar::
(23)
Lin GCS, Samuel PSH (2003). China’s land resources and land-use change: insights from the 1996 land survey. Land Use Policy 20: 87-107.
::CrossRef::Google Scholar::
(24)
Liu JY, Liu ML, Deng XZ, Zhuang DF, Zhang ZX, Luo D (2002). The land-use and land-cover change database and its relative studies in China. Journal of Geographical Sciences 12 (3): 275-282.
::CrossRef::Google Scholar::
(25)
Lu DD (1997). Report on regional development of China. Commercial Press, Beijing, China, pp. 333. [in Chinese with English summary]
::Google Scholar::
(26)
Mander U, Helming K, Wiggering H (2007). Multifunctional land use: meeting future demands for landscape goods and services. In: “Multifunctional land use: meeting future demands for landscape goods and services” (Mander U, Wiggering H, Helming K eds). Springer, Berlin, Germany, pp. 1-13.
::Google Scholar::
(27)
Ministry of Land and Resources (2004). Regulation on implementing protection of cultivated land nationwide. Beijing, China, pp. 38. [in Chinese with English summary]
::Google Scholar::
(28)
National Bureau of Statistics (2010). Statistic Year Book of China. China Statistic Press, Beijing, China, pp. 1032. [in Chinese]
::Google Scholar::
(29)
OGC (2007). Guyuan yearbook 2007. Office of Guyuan Chorography, Chinese Cultural and Historical Press, Beijing, China, pp. 443. [in Chinese with English summary]
::Google Scholar::
(30)
Parés-Ramos I, Gould W, Aide TM (2008). Agricultural abandonment, suburban growth, and forest expansion in Puerto Rico between 1991 and 2000. Ecology and Society 13 (2): 1.
::Online::Google Scholar::
(31)
Ramankutty N, Foley JA, Olejniczak NJ (2002). People on the land: changes in global population and croplands during the 20th century. Ambio 31 (3): 251-257.
::Online::Google Scholar::
(32)
SBGC (2011). Report of the 6th population census of Guyuan. Statistical Bureau of Guyuan City, Guyuan city of Ningxia Hui Autonomous Region, China, pp. 38. [in Chinese with English summary]
::Google Scholar::
(33)
State Council (1998). Circular for conservation of forest resources and prohibition of deforestation and occupation of forest land. State Council, Beijing, China, pp. 2. [in Chinese with English summary]
::Google Scholar::
(34)
Tan S, Peng B (2003). Examination of cultivated land utilization to guarantee our grain security and its recent adjustment thoughtfulness. Economic Geography 23 (3): 371-374. [in Chinese with English summary]
::Google Scholar::
(35)
Tao R, Xu ZG, Xu JT (2004). Grain for green project, grain policy and sustainable development. Social Sciences in China 6: 25-38. [in Chinese with English summary]
::Google Scholar::
(36)
Thapa GB, Rasul G (2006). Implications of changing national policies on land use in the Chittagong Hill Tracts of Bangladesh. Journal of Environmental Management 81 (4): 441-453.
::CrossRef::Google Scholar::
(37)
Turner BL, Moss RH, Skole D (1993). Relating land use and global land-cover change: a proposal for an IGBP-HDP core project. IGBP Report No. 24/HDP Report No. 5, International Geosphere-Biosphere Programme, Stockholm, Sweden, pp. 65.
::Online::Google Scholar::
(38)
Uchida E, Rozelle S, Deng XZ, Huang JK (2004). Changes in cultivated land due to urbanization: China in an international perspective. CCAP Working Paper, Paper Presented in CCICED, Agricultural and Rural Development Working Meeting, Center for Chinese Agricultural Policy, Chinese Academy Sciences, Beijing, China, pp. 56.
::Google Scholar::
(39)
Verburg PH, Chen YQ, Veldkamp TA (2000). Spatial explorations of land use change and grain production in China. Agriculture, Ecosystems & Environment 82: 333-354.
::CrossRef::Google Scholar::
(40)
Verburg PH, Soepboer W, Veldkamp A, Limpiada R, Espaldon V (2002). Modeling the spatial dynamics of regional land use: the CLUE-S model. Environment Management 30: 391-405.
::CrossRef::Google Scholar::
(41)
Verburg PH, Schulp CJE, Witte N, Veldkamp TA (2006). Downscaling of land use change scenarios to assess the dynamics of European landscapes. Agriculture, Ecosystems & Environment 114: 39-56.
::CrossRef::Google Scholar::
(42)
Wang WG, Pang XX, Song YX, Li RZ (2012). The spatial differential features of construction land changes in China. Areal Research and Development 31 (1): 110-115. [in Chinese with English abstract]
::Online::Google Scholar::
(43)
Xi J, Cao MM (2008). Grain-for-Green program and its impacts on food security in hilly and gully area of Loess plateau in northern Shaanxi province: take Mizhi county in Yulin city as a case. Journal of Arid Land Resources and Environment 22 (8): 165-169. [in Chinese with English summary]
::Google Scholar::
(44)
Xu JT, Cao YY (2001). Converting steep cropland to forest and grassland: efficiency and prospects of sustainability. International Economic Review 2: 56-60.
::Google Scholar::
(45)
Xu ZG, Bennett MT, Tao R, Xu JT (2004). China’s sloping land conversion programme four year on: current situation and pending issues. International Forestry Review 6: 317-327.
::CrossRef::Google Scholar::
(46)
Yang H, Li XB (2000). Cultivated land and food supply in China. Land Use Policy 17: 73-88.
::CrossRef::Google Scholar::
(47)
Yang L, Xie GD, Zhen L, Leng YF, Guo GM (2005). Spatio-temporal changes of land use in Jinghe watershed. Resources Science 27(4): 26-32. [in Chinese with English summary]
::Online::Google Scholar::
(48)
Yang L (2006). Landscape pattern changes and driving force analysis in Jinghe Watershed. PhD Dissertation, Graduate University of Chinese Academy of Sciences, Beijing, China, pp. 108. [in Chinese with English summary]
::Google Scholar::
(49)
Zhang LC, Wu CC (1997). The forecast of the regional pattern changes in China’s grain production in the early 21st century. National Territory and Land Economics 1: 16-19. [in Chinese with English summary]
::Google Scholar::
(50)
Zhang RY, Yang XL (1998). History and brief analysis of drying up in Yellow River. Yellow River 20 (10): 38-40. [in Chinese with English summary]
::Google Scholar::
(51)
Zhang XB, Jiang SL, Shangguan ZP, Wang LX, Liu B (2006). Study on the eco-agriculture construction and sustainable development in Jinghe watershed. Agricultural Research in the Arid Areas 24 (3): 138-142. [in Chinese with English summary]
::Google Scholar::
(52)
Zhen L, Cao S, Cheng S, Xie G, Wei Y, Liu X, Li F (2010). Arable land requirements based on food consumption patterns: Case study in rural Guyuan District, Western China. Ecological Economics 69 (7): 1443-1453.
::CrossRef::Google Scholar::
(53)
Zhen L, Cao S, Wei Y, Dilly O, Liu X, Li F, Koenig H, Tscherning K, Helming K (2009). Comparison of sustainability issues in two sensitive areas of China. Environmental Science & Policy 12 (8): 1153-1167.
::CrossRef::Google Scholar::
(54)
Zhen L, Xie GD, Yang L, Cheng SK (2007). Land use dynamics, farmers’ preference and policy implications in the Jinghe watershed of remote Northwestern China. Outlook on Agriculture 32 (2): 127-135.
::CrossRef::Google Scholar::
(55)
Zhen L, Xie GD, Yang L, Cheng SK, Guo GM (2005). Land-use change dynamics, driving forces and policy implications in Jinghe watershed of Western China. Resources Science 27 (4): 33-37. [in Chinese with English summary]
::Google Scholar::

Download

Paper Contents

Paper Sections

Paper Figures

Paper Tables

 
 
Close

 

Zhen L, Deng X, Wei Y, Jiang Q, Lin Y, Helming K, Wang C, König HJ, Hu J (2014).
Future land use and food security scenarios for the Guyuan district of remote western China
iForest - Biogeosciences and Forestry 7: 372-384. - doi: 10.3832/ifor1170-007
Close
First Previous Next Last
 
Close
© iForest

Download Reference

Paper ID# ifor1170-007
Title Future land use and food security scenarios for the Guyuan district of remote western China
Authors Zhen L, Deng X, Wei Y, Jiang Q, Lin Y, Helming K, Wang C, König HJ, Hu J
Format
Close Download