1 Introduction

Global surface temperature has increased 0.8 °C during the past century; four thirds of this increase occurred in the last three decades (Hansen et al. 2006). The acceleration in global warming and its associated changes in precipitation have already affected the ecosystems, biodiversity and human systems on Earth (Kotir 2011). Among the many impacts of climate change, the risk to agriculture is considered one of the most significant. Climate change affects agriculture and the food production system in many ways (Godfray et al. 2011). Crop production is affected by climatic variables such as rising temperatures, changing precipitation regimes and increased atmospheric CO2 levels; it is also affected by biological variables such as the lengths of the crop growth periods and the crop cycle. A longer life cycle was one of the most widely observed biological changes in the response of crops to climatic warming across the Northern Hemisphere during the twentieth century (Steltzer and Post 2009). Overall, the yields of wheat and maize have responded negatively to the recent warming since the 1980s at the global scale, although the yield response signals of other crops (e.g., rice) have still been unclear (Lobell and Field 2007). Global food production may need to increase by as much as 70 % in 2050 for a world of 9 billion people (Schmidhuber and Tubiello 2007). It will be a major challenge to fulfill this global food security goal under climate change. This is especially true for China as the world’s most populous country and the third largest in land area.

China is among the most affected countries by climate change (Turral et al. 2011). Although much research in China has focused on the impact of climate change on food production (Lin et al. 2005; Lobell et al. 2008; Piao et al. 2010), systematic studies that directly link climate change and food security are still lacking, especially at the national scale. Essentially, food security assessment needs to consider parameters involving not only food production but also food consumption, resulting in added complexity hence difficulty in data collection and model building (Schmidhuber and Tubiello 2007). The objectives of this paper are to: (1) simulate the yield of major food crops in China under the IPCC SRES A2 and B2 scenarios in 2030/2050; (2) assess the impact of climate change on food security in China, involving scenarios of population size, urbanization rate, cropland area, cropping intensity, technology development, etc.; and (3) formulate policy options to ensure China’s food security under climate change in 2030/2050.

2 Materials and methods

2.1 Research framework

A five-step approach has been adopted in this research. First, climatic, crop, soil, management, and socio-economic data were collected, manipulated, and used to simulate the yields of food crops. Second, the simulated yields were compared to the observed yields in order to validate the simulation process. Third, the food production capacities in 2030 and 2050 were estimated based on the scenarios of population growth, urbanization rate, cropland area, cropping intensity, new crop varieties, and management improvements, which were all designed in line with the socio-economic development pathways associated with climate change scenarios. Fourth, a food security index (FSI) was computed based on the analysis of the supply–demand relationship of food produced from agricultural crops. Finally, the effect of climate change on the FSI was quantitatively evaluated, and policy options toward long-term food security in China were formulated (Fig. 1).

Fig. 1
figure 1

The flow chart of the research framework. Food security assessment is conducted from a supply-demand point of view involving a series of scenarios

2.2 Data manipulation

To facilitate the simulation of the yields of major food crops including rice (Oryza sativa L.), wheat (Triticum aestivum L.), and maize (Zea mays L.), the study area was divided into grid cells of 404 rows by 482 columns, each of 100 km2. The cell number was used as an index to store input parameters in a relational database. Daily climatic variables of maximum and minimum temperature, precipitation, and sunshine duration from 680 stations were provided by the Chinese Meteorological Data Center over the period 1981–2000. These stations form a spatially representative climatic network over croplands in China (Fig. 2a). Daily solar radiation was derived from a linear relationship with n/N (actual to maximum ratio of daily hours of sunshine) as given by Pohlert (2004). Each grid cell was associated with the nearest meteorological station. The soil parameters (clay–silt–sand contents, bulk density, pH, and organic carbon content) were obtained from the National Soil Survey Office; soil moistures at field capacity and wilting point were derived using pedo-transfer functions (Pawls et al. 1982). A GIS procedure was followed to rasterize the soil mapping units into 10 × 10 km grid cells (Knox et al. 2000). Crop (cultivar, sowing date, length of crop cycle and phenological subdivisions, leaf area index at max. growth rate, harvest index, and rooting depth) and management parameters (sowing density, row width, transplanting date, fertilizer application date and quantity, irrigation date, and volume) were collected per experimental site (Xiong et al. 2008a) and assigned to the grid cells using the same approach as for the meteorological parameters. In total 522 yield trials from 119 experimental stations were used for, e.g., rice (Xiong et al. 2008a).

Fig. 2
figure 2

Annual yield growth rate of maize in China under the SRES A2 and B2 scenarios, calculated from the average simulated yields between the periods of 1961–1990 and 2011–2020 (a, d), between 2011–2020 and 2021–2030 (b, e), and between 2021–2030 and 2031–2040 (c, f). The maize cropping area is extracted from the MIRCA2000 dataset (Portmann et al. 2010) and upscaled to the spatial scale of 10 km. Black dots (a) represent climatic stations. NE Northeast, NCP North China Plain, SB Sichuan Basin

2.3 Yield simulation

The yields of rice, wheat, and maize were simulated per grid cell of cropland using the CERES-Rice, Wheat and Maize model, respectively, with the CO2 fertilization effect included (Lin et al. 2005; Lobell and Field 2008). The spatial distribution of cropland (Fig. 2) was assessed based on the MIRCA2000 dataset (Portmann et al. 2010). The CERES models are dynamic process-based mechanistic models that simulate daily phenological development and growth in response to environmental (climate and soil) and management factors (crop variety, planting conditions, N fertilization, irrigation, etc.). The CERES models have been validated in a wide variety of environment, ranging from temperate to tropical, across the world (Jones et al. 2003).

The models run for three times to simulate the irrigated and the rainfed yields at the grid scale. The first simulation run was performed for the calibration period of 1981–2000; the second run was performed for the reference period of 1961–1990, while the third run was performed for the projection period of 2011–2040. Climatic data of the reference and the projection periods was derived and downscaled from the PRECIS, or Providing Regional Climates for Impact Studies (Xu et al. 2006), model for China. The baseline scenario used the PRECIS simulated daily weather data for 1961–1990 with the CO2 concentration set at ~330 ppm. Two climate change scenarios, namely A2 and B2, for 2011–2040 were also generated. The IPCC SRES A2 scenario describes a very heterogeneous world of high population growth, slow economic development, and strong regional cultural identities, while the B2 scenario reflects a heterogeneous world with diverse technological change, low population growth and emphasis on local solutions to economic, social, and environmental sustainability problems. The SRES A2 and B2 scenarios were selected to represent two contrasting development pathways, pessimistic versus optimistic, for China, respectively.

The area-weighted average crop yield Y (t ha−1) is calculated at the grid scale from the simulated rainfed (Y rain) and irrigated (Y irri) yields (Ye et al. 2008):

$$ Y = p \cdot {Y_{{irri}}} + \left( {1 - p} \right) \cdot {Y_{{rain}}} $$
(1)

where p ∈ [0, 1] is the irrigation probability, \( p = {{{{S_{{irri}}}}} \left/ {{\left( {{S_{{irri}}} + {S_{{rain}}}} \right)}} \right.} \)with S irri and S rain being the irrigated and rainfed areas, respectively. The irrigated and rainfed crop areas were evaluated based on the MIRCA2000 dataset (Portmann et al. 2010).

2.4 Scenario building

The socio-economic and agronomic scenarios of population size, urbanization rate, cropland area, cropping intensity, and technology development (including new crop varieties and management improvement) were designed in line with the pessimistic versus optimistic development pathways represented by the SRES A2 and B2 emission scenarios, respectively. The settings of these scenarios are given in Table 1.

Table 1 Settings of the socio-economic and agronomic scenarios associated with the IPCC SRES A2 and B2 emission scenarios, including population size, urbanization rate, cropland area, multi-cropping index (MCI), and technology development in terms of yield enhancements resulting from crop variety and management improvements

The logistic growth model was employed on the historical Chinese population data (NBSC 2010) to quantify the future trends of population growth. The UNFPA (2010) projected size of 1.42 billion in 2050 and plateau time in 2025 set the key parameters of the low growth variant of the population scenario, corresponding to the SRES B2 emission scenario. The urbanization scenario was based on the historical trend of urban population growth of 1 % year−1 observed during 1978–2009 and on China’s national development strategy (Change and Brada 2006). Since the year 2000, cropland in China has been lost at a rate of 1.45 million hectares per year. Urbanization was known as the major driver of it. The future evolution of cropland area was evaluated by coupling the loss of cropland with urbanization rate, using a linear regression model (Ye and Van Ranst 2009). Cropping intensity, expressed in terms of the multi-cropping index (MCI) or the ratio between total sown area and the cropland area, increased steadily at a rate of 0.9 % year−1 during the past 20 years. Although double cropping or even triple cropping is a common practice in southern China, the national MCI average in 2009 was only 130 %. A study conducted by the National Office of Agricultural Regional Planning (NOARP 1997) suggested that the average MCI in China could reach 170 %, with a theoretical potential of 190 %.

Two forms of technology development were considered in scenario building. One was the breakthroughs in crop breeding technologies (Tester and Langridge 2010), which produces higher-yielding varieties, and the other was management improvements in terms of, e.g., raised efficiency of fertilizer use and associated field management techniques. Field trials from China indicated that use of inorganic fertilizer increased the yield of major food crops by 47.8 % (Jin 2011), although over fertilization also caused a number of eco-environmental problems such as GHG emission (Tester and Langridge 2010) and soil acidification (Guo et al. 2010). With bio-technology breakthroughs, gradual adoption of higher-yielding varieties is considered capable of raising crop yields at an average rate of 0.8 % year−1 before 2030 and 0.5 % year−1 thereafter (Zhao et al. 2008). Based on statistical modeling using data from 5 Chinese provinces during 1988–2007, Li et al. (2011) found that 0.7 % of the yield growth rate can be attributed to management improvements. These yield enhancement technologies are fully implementable under SRES B2, due to its intrinsic capacity to accommodate technological changes, and only partially implementable under A2 (Table 1).

2.5 Food security assessment

Food security was defined by the World Food Summit in 1996 as a “situation that exists when all people, at all times, have physical, social and economic access to sufficient, safe and nutritious food that meets their dietary needs and preferences for an active and healthy life.” Since then, numerous measures have been used to quantify the overall status and regional discrepancies in global food security (or insecurity), including the FAO indicators of undernourishment (Schmidhuber and Tubiello 2007). However, none of these have been widely accepted for the populous Asian countries. As in other Asian countries, food security in China almost exclusively means food self-sufficiency or even grain self-sufficiency (Jiang 2008). The grain self-sufficiency level of 95 % has recently been adopted as a strategic goal for maintaining food security in China.

It is important to note that: (a) China’s population is too large to rely on the world food market for its food supply, (b) the swift and fair distribution of large quantities of imported food is a major challenge given the scale of China, and (c) food is a special good and its self-sufficiency has domestic political and economic significance (Jiang 2008). Therefore, the relative food surplus is defined in this paper as the food security index, or FSI, in an attempt to measure the general status of food security in China:

$$ FSI = \frac{{{{s} \left/ {g} \right.} - d}}{d} \cdot 100 $$
(2)

where s is the per capita food supply, d is the per capita food demand, and g is the expected food self-sufficiency level (g = 0.95).

3 Results and discussion

3.1 Model validation

The performance of the CERES crop models in predicting crop yields in space and time has primary importance in long-term decision making. Thus, a lot of efforts have been made to adequately calibrate the CERES models for use in regional yield simulations in China and to thoroughly evaluate their performances in regional climate impact studies, as separately reported (Xiong et al. 2008a; 2008b). Year-by-year comparisons between the simulated yields (1981–2000) and the observed yields, which were collected and detrended (Easterling et al. 1996) from the agricultural meteorological experimental station network in China, showed that the simulated yields were in good agreement with the observed yields. In the case of maize (Xiong et al. 2007), for instance, yield records from four experimental stations representing major maize growing regions across China were used for comparison; they are Changjizhou (44°N, 87.43°E) in Xinjiang province, Zhengzhou (34.72°N, 113.65°E) in Henan province, Guangyuan (32.43°N, 105.85°E) in Sichuan province, and Harbin (45.75°N, 126.77°E) in Heilongjiang province. The results showed that 96 % of the variations in the observed yields was reproduced by the simulated yields, or R 2 = 0.96 (p < 0.01). The spatial pattern of the variability in observed yields was reasonably well captured by the model-simulated yields. A case study in Northeast China’s Jilin Province, in which the county-level census maize yields from 1981 to 2000 were first downscaled to the grid level and then compared to the simulated yields, showed that the simulated and the census yields were correlated across all grids (R 2 = 0.24, p < 0.05), and the simulated annual yields correlated to observed annual yields in 69 % of the grids at the confidence level of 95 %. Overall, these results are encouraging in regional yield simulation and provide reasonable confidence in the simulation results in space and time.

3.2 Yield growth rate under climate change

The rainfed and irrigated yields of major food crops were simulated annually at the grid scale during the reference period of 1961–1990 and during the projection period of 2011–2040. The simulated rainfed and irrigated yields were combined to derive the yield of a particular crop (e.g., maize) using Eq. 1. The multi-year mean of the simulated yield was calculated for the reference period and for each decade of the projection period. The compound annual yield growth rate (%), r, was computed from the mean yields using the following equations:

$$ {Y_2} = {Y_1} \cdot {\left( {1 + {{r} \left/ {{100}} \right.}} \right)^n} $$
(3)

thus

$$ r = 100 \cdot \left[ {\exp \left( {\frac{1}{n}\log \left( {{{{{Y_2}}} \left/ {{{Y_1}}} \right.}} \right)} \right) - 1} \right] $$
(4)

where Y 1 is the average crop yield during period 1, Y 2 is the average yield during period 2, and n is the distance between period 1 and period 2, counted in number of years from the middle year of period 1 to the middle year of period 2. For example, n = 40 between the periods 1961–1990 (centered at 1975) and 2011–2020 (center at 2015). As an example, compound annual yield growth rate of maize, evaluated either between the reference and the projection periods or at the decadal intervals within the projection period, under the SRES A2 and B2 scenarios are shown in Fig. 2. The crop-specific yield growth rates, aggregated at the national scale, are summarized in Table 2.

Table 2 Annual yield growth rates (%) of rice, wheat, and maize in China under SRES A2 and B2

The maize yield in Northeast China, for example, is projected to respond negatively to climate change during the 2015 (as the mid-point of the 2011–2020 decade) to 2025 (mid-point 2021–2030) period under the SRES A2 scenario, as shown by the dominating brownish colors in Fig. 2b. In another major region in food production in China, the Sichuan Basin, the impact is projected to be marginal to negligible. Unsurprisingly, similar spatial disparities are also observed for other periods of time (Fig. 2c) or under other scenarios (Fig. 2e, f). Overall, the maize yield is projected to increase at 0.2–0.3 % annually at the national scale during 2011–2040 under A2. This translates into a yield change of +10 % in 2040 over the 1961–1990 average. Under the B2 scenario, the maize yield is projected to first increase at 0.5 % annually during 2011–2030 and then decrease slightly at an annual rate of −0.05 % during 2031–2040. Temporal disparities in yield growth rates of wheat and rice under climate change are also observed (Table 2).

Our results show largely positive (or less negative) effects of climate change on crop yield in China, contrasting with earlier estimates in many cases. Using yield transfer functions, Parry et al. (2004), e.g., demonstrated that climate change would decrease cereal yield in China by 2.5–10 % and 2.5–5 % in 2020 under A2 and B2 scenarios, respectively, and by 5–10 % in 2050 under either A2 or B2. A recent assessment using the EPIC model (Wu et al. 2011) also revealed negative impacts of climate change on crop yield in China, especially in the most important food production region of North China Plain. One fundamental difference between this research and those mentioned above is that CO2 fertilization (Gosling et al. 2011; Lobell and Field 2008) was intentionally included in this simulation, while CO2 fertilization was excluded from others such as those of Parry et al. (2004) and Wu et al. (2011). The rationale to consider CO2 fertilization is that the rise of atmospheric CO2 concentration has already been an observed fact and it will continue to rise in the foreseeable future, despite uncertainties about the magnitude of this rise (Lin et al. 2005; Long et al. 2006). Elevated CO2 concentrations stimulate photosynthesis, by a margin of ~0.08 % for rice and wheat and ~0.05 % for maize per 1 ppm increase, leading to increased plant productivity and modified water and nutrient cycles. Past studies suggested that crop yield tends to increase under higher CO2 concentrations (Tubiello et al. 2007). Compared with the current atmospheric CO2 concentrations of ~380 ppm, crop yields increase at 550 ppm CO2 in the range of 10–20 % for C3 crops (e.g., rice and wheat) and 0–10 % for C4 crops (e.g., maize). The magnitude of the positive yield effects found here (Table 3, no technology) fell well in this range. Although the robustness of CO2 fertilization is being debated (Gosling et al. 2011; Long 2012), its yield effect has been confirmed by a variety of field experiments such as the free-air carbon dioxide enrichment (FACE) experiment.

Table 3 Crop yield and food production capacity under the considered socio-economic and agronomic scenarios in association with the SRES A2 and B2 emission scenarios to explicitly show the effect of technology development on crop yield and production capacity

The temporal variabilities of the simulated yield effects (Table 2, Fig. 2) at the decadal time scale can be explained as large-scale responses of crops to future temperature and precipitation trends as governed by climate scenarios used (Xiong et al. 2007). Historically, spatial variations of precipitation in major agricultural regions in China showed a 10-year north-south oscillation between the Yellow River Basin and the Yangtze River Basin regions. It has been observed that in the years when the Yellow River Basin receives more precipitation than the Yangtze River Basin, the Yangtze River Bain usually receives higher radiation due to lower chance of cloud cover. During these years, rice in the Yangtze River Basin could benefit from less precipitation plus higher radiation, while wheat and maize in the Yellow River Basin could benefit from more precipitation. On the other hand, in the years when the Yangtze River Basin receives more precipitation, the Yellow River Basin usually experiences higher frequency of draughts. The climatic conditions in these years are undesirable for both rice (higher precipitation, lower radiation) in the Yangtze River Basin and wheat and maize (lower precipitation, higher radiation) in the Yellow River Basin, as observed during late 1990s. This tempo-spatial pattern was reproduced by most regional climate models (RCMs) in China such as PRECIS (Xu et al. 2006) which was used to generate climate scenarios in this research. The decadal temporal pattern of yield variability described here is usually more obvious under scenarios accommodating more extreme events (e.g., A2) than under more mild scenarios (e.g., B2), as confirmed by Table 1.

3.3 Yield

In order to isolate the effect of climate change and to maintain comparability, the yields of rice, wheat and maize during the projection period of 2011–2040 were evaluated by applying the annual yield growth rate on the census yield in the baseline year of 2009 (NBSC 2010). The yields during 2041–2050 were extrapolated from the projected yields in 2040, assuming the same yield growth rate during 2041–2050 as during 2031–2040. The projected yields were then aggregated to the national scale and averaged at the decadal time scale and given in Table 3.

Climate change is simulated to have moderate positive effects on the yields of major food crops in China at the decadal time scale from 2020 through 2050 (Table 3, no technology). The maize yield, for instance, is projected to increase by ~10 %, from 5.3 t ha−1 in 2009 to 5.8 t ha−1 in 2050, under the A2 scenario. Under the B2 scenario, a 4 % increase in maize yield is projected. Overall, crop yields are projected to increase in 21 out of 24 cases (3 crops by 4 decadal intervals by 2 scenarios, Table 3, no technology). The wheat yield is projected to be mostly stagnated or decrease by a small margin under either A2 or B2. The average yield of all food crops taken together (i.e., staple grains plus tubers and beans; the proportion of tubers and beans in total food production, 10 % in 2009, was assumed constant throughout the projection period) is projected to increase by 11 % from 4.9 t ha−1 in 2009 to 5.4 t ha−1 in 2050 under A2, assuming that the sown area proportions of individual crops in 2009 are kept unchanged during the entire projection period, while this yield increase is projected to be 4 % under B2, meaning that short-to-medium-term yield growth is more likely to be achieved under the A2 scenario which assumes higher emission levels (320 % more CO2 in 2100 than in 2000), as also observed by others (e.g., Parry et al. 2005). Agricultural production will likely benefit from a more balanced development pathway as assumed under B2, but this benefit may probably only be achieved over longer terms.

3.4 Food production capacities

China’s food production capacity (Table 3, with technology development) tends to increase over the projection period under both A2 and B2 scenarios. China will be able to achieve a production of 572 and 635 million tons (MT) from food crops under the A2, while this production is 615 and 646 MT under the B2 scenarios, in 2030 and 2050, respectively, compared to the production of 531 MT in 2009. It is clear that this increase in food production capacity is overwhelmingly attributed to production intensification as indicated by a 30 % increase in the multi-cropping index values during 2009–2050, especially under the B2 scenario (Table 3). The moderate positive effect of climate change on food yield (Tables 2 and 3) will not be able to offset the negative effect of the loss of cropland on food production in China over the projection period. China would only be able to achieve a production of 497 MT in 2050 under the B2 scenario, i.e., a 7 % drop from the 2009 level, should the multi-cropping index value be kept at the same level as in 2009, while this production would be 550 MT in 2050 under A2, meaning a 4 % increase over the baseline level. This mild increase in total food supply is insufficient to feed 17 % more people (Table 1). As a result, per capita supply would drop from ~400 kg to ~350 kg under either A2 or B2 scenario. This shows that agricultural intensification will be the inevitable choice for countries like China to ensure food security over the long run (Godfray et al. 2011; Tilman et al. 2002).

3.5 Food security analysis

The food security index (FSI) values in China, evaluated from a food supply–demand point of view (Ye and Van Ranst 2009) using census and estimated data for the pre-2009 and model simulated crop yields under socio-economic and agronomic scenarios in association with the SRES A2 and B2 climate change scenarios, are shown in Fig. 3. Historical variations in food security are well captured by the FSI curve. China’s food security status was significantly improved soon after the long-lasting wars that ended in the late 1940s. At the end of the first 5-year plan, the FSI increased from −2.4 in 1949 to 31.4 in 1957, showing that the supply–demand relationship turned from a 2.4 % deficit to a 31.4 % surplus. The peak FSI value of 38.5 appeared in 1984, coinciding with the record harvest of 390 MT in the same year. Although higher productions (~500 MT) were achieved consecutively during 1996–1999, the FSI values in the same period were not higher than that of 1984, reflecting the combined effects of a larger population and a higher standard of living. Extreme climatic events and natural hazards, which caused notable production losses during 2000–2003, were responsible for the second largest drop in the FSI values after the period of the Great Leap Forward (1957–1961). China has achieved record harvests for six consecutive years since 2004, reaching the level of 530 million tons in 2009. However, the average FSI level of 18.8 during the first decade of the twenty-first century is considerably lower than the average level of 31.9 in 1990s or 26.3 in 1980s, showing the dragging effect of steadily increasing consumption levels on FSI. This suggests that food utilization (safe, balanced and nutritious food, etc.; see Schmidhuber and Tubiello 2007) is gaining momentum and attention is needed on how to integrate it into existing food security assessment frameworks (e.g., Fig. 1) which took food availability as the primary indicator.

Fig. 3
figure 3

Relative food surplus (%) in per capita terms, or the food security index (FSI), as evaluated using census data during 1949-2009 and as projected in 2030 and 2050. Horizontal bars represent the average FSI levels at the decadal time scale

Statistical analysis was conducted to reveal the cause–effect relationship between crop yield and the FSI (Fig. 4). Results show that the correlation between crop yield and FSI is much weaker than the correlation between the year-to-year changes of the crop yield and the FSI. Only 22 % of the variations in the absolute values of the FSI can be explained by the variations in the average yield of all food crops taken together, or R 2 = 0.22 (p < 0.001), as shown in Fig. 4c. However, as much as 82 % of the variations in the year-to-year changes of the FSI values can be explained by the variations in the annual growth rate of the crop yield (Fig. 4f). It is important to note that the year-to-year changes of FSI, given by the first difference of the FSI series as inspired by the approach of Lobell and Field (2008) in relating CO2 growth rate to crop yield annual changes, is measured in terms of percent change of relative food surplus by definition, which is essentially the growth rate of the relative food surplus. The statistical results suggest that the annual growth rate of crop yield is a much better indicator of food security than crop yield per se, meaning that yield improvement on the yearly basis has great significance in ensuring food security for countries with a growing population, as in the case of China during the pre-2025 or pre-2040 era, depending on the choice population growth pathways associated with the B2 or A2 climate change scenarios, respectively (Table 1).

Fig. 4
figure 4

Cause-effect relationships between crop yield and the food security index (FSI) values show that yield growth rate is a better indicator of food security than absolute yield. a Time plot of FSI; b time plot of average yield of food crops; c scatter plot of FSI versus yield, overlaid with the trend line (R 2 = 0.22); d time plot of the yield-to-year FSI changes (ΔFSI), evaluated as the first difference (Lobell and Field 2008) of the FSI series; e time plot of the annual growth rate of the average yield of food crops (ΔYield), evaluated as the ratio (%) between the first difference of yield and absolute yield; f scatter plot of ΔFSI versus ΔYield (R 2 = 0.82)

The Chinese population is projected to plateau and decrease within the course of the projection period, despite the difference on exact timing (Table 1). The FSI is predicted to drop sharply from 24.2 in 2009 to 10.2 and −4.5 in 2030 under the B2 and A2 scenarios, respectively (Fig. 3). This drop can be explained by the relative importance of population versus other, yield enhancement factors. Under the A2 scenario, population grows quickly by a margin of 240 million people during 2009–2030, which is equivalent to the population size of Indonesia as the fourth most populous country of the world. Alongside the high population growth, the socio-economic and agronomic development takes a sub-optimal pathway. The multi-cropping index, for example, increases by 9 units from 130 % in 2009 to 139 % in 2030. Due to the input limitations, only half of the available higher-yielding crop varieties can be adopted and crop management can only be exercised half efficiently as compared to the B2 scenario. Aggregate, technology development will result in a 0.75 % yield enhancement effect annually. Consequently, average crop yield will increase by 7 % in 2030 as compared to 2009 with the yield effect of climate change included. Despite all these minor positive effects, per capita supply decreases all the way down from 398 kg in 2009 to 363 kg in 2030. In contrast, per capita consumption increases steadily from 337 kg in 2009 to 400 kg in 2030 (not shown). This shows that population growth has enormous impact on food security (Hopfenberg and Pimentel 2001). During the period of 2030–2050, the FSI is predicted to increase from −4.5 in 2030 to 7.1 in 2050 under the A2 and from 10.2 to 20.0 under the B2 scenario, respectively, coinciding with the projected decrease in population size during the period (Table 1), reaffirming the controlling effect of population growth on food security in populous countries like China. Therefore, as a countermeasure for food security, the yield growth rate should be maintained at a level higher than that of population growth. Hopfenberg and Pimentel (2001) argued with their analysis that crop yield needs to grow 1–2 % faster than population growth rate in order to secure food supply. Alternatively, population control should be prioritized where proper yield growth rate cannot be sustained (Ehrlich and Ehrlich 2009).

A comparison between the simulated FSI curves under the A2 and B2 associated scenarios shows that the socio-economic and agronomic pathways in association with the A2 and B2 emission scenarios have significant impact on FSI. The average distance between these two curves during 2011–2050 is evaluated to be 13 units, meaning that food surplus linked with a more balanced development pathway (B2, low population growth, stimulation of technological change and emphasis on environmental sustainability) is 13 % higher than a more stagnated pathway (A2, high population growth, low economic development, low regional coordination). This difference can translate into 76 million tons of additional grain harvests in 2030 under B2. The average distance between the simulated FSI curve under B2, for example, and the FSI curve under B2 but excluding technology development during 2011–2050 is evaluated to be 7 units (Fig. 3). This reveals that yield improvements realized by technology development alone may have raised the food surplus level during 2011–2050 under B2 by 7 %, suggesting that technology development is one critical means to raise food security level through yield growth rate maintenance (Alston et al. 2009).

3.6 Policy options

We advise the following research priorities and policy reforms in order to ensure food security under climate change: (1) Modeling. The food system is complex, and interventions often have unintended and detrimental effects on food security. There is a clear need to bring together economists and scientists in modeling the food system (Godfray et al. 2011). Biophysical-economic models using food price (Alston et al. 2009) as the coupling parameter are one of the priorities that have the potential to meet the needs of agro-environmental to socio-economic decision makings. Multi-model ensemble simulation of crop yield under climate change (Challinor et al. 2009; Gosling et al. 2011) is much needed to control and to quantify the uncertainties linked with the obtained results. Coordinated use of these models will hopefully provide a promising platform for food security assessment (Rötter et al. 2011); (2) Technology. Past experiences show that breeding and agronomic improvements have on average achieved a linear increase in global food production at an average rate of 32 million tons per year for the past few decades. Contemporary, large-scale environmental change poses great threats to the continuity of this linear trend, but it is not prohibiting. Internationally coordinated efforts are in urgent need to speed up the breeding of new crop varieties which are more efficient in utilizing water (Molden and Sakthivadivel 1999) and nitrogen fertilizers (Jin 2011), and drought and heat tolerant (Challinor et al. 2010; Long 2012); (3) Extension. In most circumstances in China, extension services are the key to swift and effective implementation of new technologies. However, the extension service network needs to be reinvented to build the relevant skills base among food producers, and better economic incentives need to be introduced to make it worthwhile for farmers to adopt new technologies and/or practices. The extension services need to be revitalized to guarantee the successful adoption of such technologies into the smallholder systems in China’s food production (Piao et al. 2010).

4 Conclusions

An integrated assessment of food security in China based on the simulated effects of climate change on crop yield in 2050 was presented in this paper. The assessment involved a range of socio-economic and agronomic scenarios which were developed following two contrasting development pathways associated with the IPCC SRES A2 and B2 climate change scenarios. The food security index (FSI, or relative food surplus per capita) was found to be a simple yet powerful and robust tool in quantitative cause–effect food security analysis. The obtained results showed that climate change only has a moderate positive effect on food security as compared to other factors like cropland area and population growth. Overall, the socio-economic development pathways exert great impact on the future trends of food security. Particularly, the more balanced development pathway with emphasis on environmental sustainability in association with the B2 emission scenario was found much superior in ensuring food security than the other A2-associated pathway. Using regression analysis that correlates FSI and crop yield, we found that yield growth rate is a much better indicator of food security than crop yield per se, meaning that yield improvement on the yearly basis has great significance in ensuring food security for countries with a growing population like China. The relatively lower than expected FSI values during the first decade of the twenty-first century suggested that per capita food consumption in China is undergone a steady growth. Here we argue that the gravity center of food security in China starts shifting from the availability aspects to the utilization aspects. More attention is therefore needed on food safety and balanced nutrition issues in order to ensure food security in China over the long run.