1. Introduction
With global climate change, continuous population growth, and rapid urbanization, food security issues and policies remain a subject of concern to the international community. China feeds about 18% of the global population using 8% of its farmland [
1]. The fundamental reality of more people and less farmland in China demonstrates that food security is crucial to the lasting political stability in China. With the impact of the epidemic and rising uncertainty in the international trade environment, the issues of ensuring baseline of food security and grasping the initiative in food security have become more prominent. As an important factor affecting food security, the change of farmland area has received more attention. Over the past 40 years of the reform and opening up in China, irreversible non-agricultural changes in a large amount of farmland have taken place according to the progress of fast industrialization and urbanization, which has led to the decrease of farmland area, and a threat to food security [
2]. Strict observation of the red line of 1.2 million km
2 of farmland has become a national political task to ensure food security. However, with the continuous increase of total population and urbanization, it is very difficult to increase the food supply by increasing farmland area. On the contrary, it is more feasible to ensure national food security by improving the level of intensive use of existing farmland [
3]. Multiple cropping is an important aspect of the intensive use of farmland [
4]. From 1986 to 1995, the increased grain yields attributed to the increasing multiple crop index (MCI) of farmland, accounted for one third of the average annual total grain yields (429 billion kg) in China [
5]. About 12% of global farmland applied multiple cropping in 2000. In addition, 34%, 13%, and 10% of rice, wheat, and maize crops, respectively, utilize multiple cropping, demonstrating the importance of such cropping systems for cereal production [
6]. Moreover, compared with reclaimed farmland, existing farmland possesses better production conditions. Multiple cropping is, therefore, an effective way to increase the grain yields and ensure food security [
3,
7,
8].
With the expansion of the connotation of food security, the objectives of researches on multiple cropping of farmland in different countries have also shown differences. In European and North American countries with a high level of economic and social development, the researches mainly focus on the impact of multiple cropping on pest control and soil improvement. The conclusions based on field experiments prove that increasing the level of MCI of farmland can increase diversity, thereby contributing to pest control and reducing herbicide intensity [
9]. Exploring different multiple cropping modes can effectively enhance organic matter and microbial activities in the soil, thereby developing organic agriculture and obtaining a higher income [
10]. In South America, Africa, and Asia, where the level of economic and social development is relatively low, the goal of related researches is mainly to increase grain yields [
11,
12,
13]. In recent years, however, it has begun to shift to the direction of balanced nutrition [
14]. It is worth noting that in Asia, human-land contradiction is very serious. The production system characterized by smallholder determines the necessity of increasing MCI to increase the grain yields and incomes [
15]. Therefore, Asia has become a region of focus for researches on MCI [
16,
17]. In India, the zoning map of rice MCI was drawn and used to estimate the irrigation demand of different zones to provide a scientific basis for policy evaluation [
18]. From north to south in China, great differences in terms of crop types and cropping systems are exhibited among eight temperature zones [
19,
20]. Influenced by agricultural production conditions and socio-economic development, the MCI in the major grain producing areas [
21,
22] and the rice-growing areas, where “double cropping to single cropping” is common [
23,
24], has noticeably declined in recent years.
Reasons for this are summarized into the following four aspects. Firstly, marginal incomes earned via a multi-cropping system decrease significantly as a result of the increasing production cost. MCI was changed from multiple cropping to single cropping to maximize the economic benefits [
25]. Secondly, the labor marginal incomes from non-agricultural employment are much higher than those of agricultural production for Chinese farmers. Farmers are more inclined to transfer more labor time and production resources to part-time or non-agricultural production activities [
26], thus resulting in seasonal or year-round abandonment of farmland especially in labor-intensive cash crops and regions closer to urban areas [
22]. Thirdly, more farmers may face a poor harvest after using the “double cropping to single cropping” method, since those that plant double-cropping rice may be exposed to the intensive damages of insects, birds, and animals [
27]. Fourthly, the adjustment of food policies will also cause changes of MCI through incentives and constraints on the planting behaviors of farmers. In a word, decreasing multi-cropping level and even abandoning cultivation, is a rational choice of farmers under low planting efficiency [
28].
Existing researches mainly concern the influence of MCI on food security, and use MCI as the input variable to calculate variations in farmland area and grain yields. These prove that multi-cropping system can indeed increase the outputs of corn and rice [
29], while decreasing the multi-cropping level can inhibit, and even decrease the growth of food output. This makes maintaining the self-sufficiency of cereal a challenge. However, some researches pointed out that the improvement of the multi-cropping system might influence the resource ecological environment. The practices in Pakistan prove that around 51% and 13% of water inefficiency are present under multiple and sole cropping systems, respectively [
30]. The expansion of the multi-cropping system increased agricultural greenhouse gas emissions in the North Plain and neighboring regions in China [
31], and the growth of the annual mean temperature, in return, can influence the growth of crops [
32,
33]. Evidently, pursuing high MCI blindly, and ignoring the water and temperature conditions would work against the increasing of grain yields and the sustainable development of the ecological environment [
34].
The cropping system of China is not only experiencing a decline of MCI, but also facing the risk of spatial mismatch between cropping system and natural production conditions (including water, soil, gas, etc.). Firstly, there are abundant water and heat resources in South China. Historically, the food supply pattern entailed “transport from south regions to north regions”, but now has changed to “transport from north regions to south regions”, thus increasing consumption of farmland resources [
35]. Secondly, the location of large and medium cities often highly overlaps with that of high-quality farmland [
36]. A considerable amount of farmland with high-quality water and heat resources is occupied by urban construction sprawl, while the reclaimed farmland with poor production conditions is used to compensate for the loss of high-quality farmland with fertile soil and high MCI. The imbalance of the quality and production capacity of farmland has threatened China’s food security [
37]. In this regard, some studies have measured the potential multiple cropping index (PMCI) of farmland in China based on water and temperature conditions, which is the theoretical highest MCI of farmland based on the natural environment conditions [
3]. Based on PMCI, some researches inferred the most sown area [
38] and grain yields [
1] under the optimal cropping system.
Scientific analysis of the relationship between multi-cropping system and potential multi-cropping system is conducive to deepening our understanding of farmland use and the scientific exploration of the potential of farmland, as well as providing references and supports for the implementation of a strategy that “stores foods in farmland”. Based on the normalized difference vegetation index (NDVI), this paper will analyze the spatio-temporal characteristics of the MCI of farmland in China from 2000 to 2018. The distortion of water-land resources will be judged by the gap of MCI and PMCI. Finally, suggestions will be put forward to give full play to the production potential of high-quality farmland so as to achieve a win-win situation for food security and ecological security. Compared with existing researches on MCI of farmland [
16,
17,
18,
19,
20,
21,
22,
23,
24], one of the innovations of this study is the problem of increasingly serious farmland abandonment introduced into the study of multi-cropping system. We will further divide the decline of MCI into “seasonal” abandonment and year-round abandonment, so as to respond to attentions on abandonment of farmland in China’s farmland protection system. Another innovation was the delineation of four types of grain producing zones, namely key development zone, potential growth zone, appropriate development zone, and restricted development zone, and the provision of references for optimization of food production layout and benefit compensation mechanism design.
2. Materials and Methods
2.1. Materials
NDVI, also called the standardized vegetation index, is a comprehensive reflection of vegetation type, coverage form, and growth conditions in unit pixel. The value of NDVI is determined by the vegetation coverage and leaf area index (LAI). The physical growth processes of crop sowing, seedling, heading, maturing, and harvest in a year reflect fluctuations of NDVI with time, and peaks correspond to the time phases when the biomass of crop populations is the largest. According to this principle, the MCI of farmland is gained by extracting the peaks number of NDVI in one year. NDVI ranges between minus 1 and 1. Specifically, a negative value represents that a surface is covered by cloud, water, or snow; 0 represents rocks or naked soils; a positive value indicates vegetation coverage, which increases with the increase coverage [
39]. In this study, the monthly (January to December) NDVI sequence from 2000 to 2018 is generated by the maximum value combination based on continuous time series of SPOT/VEG satellite remote sensing data. The spatial resolution of NDVI was 1 km × 1 km.
The spatial distribution data of potential multi-cropping system in China is estimated by the Global Agro-Ecological Zones (GAEZ) model developed by the FAO and IIASA together based on data of DEM, soil, farmland, and meteorological. On this basis, the ideal cropping system can be realized for the farmland. The potential multi-cropping system data includes single cropping, double cropping, and triple cropping in a year, with a spatial resolution of 10 km × 10 km.
In this study, the farmland grid data in five phases (2000, 2005, 2010, 2015, and 2018) were used to restrict the identification range of cropping system in farmland and eliminate interferences of other land use types. The spatial resolution of it is 1 km × 1 km. The number of farmland grids has been decreasing continuously since 2000, and it experienced a sharp reduction from 2005 to 2010 and since 2015 (
Figure 1).
The above three types of data are provided by the Data Registration and Publishing System of Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (
https://www.resdc.cn/data.aspx?DATAID=254, accessed on 12 October 2019). They cover 31 provinces in China, except Taiwan, Hong Kong, and Macau, as shown in
Figure 2. In addition, it is necessary to introduce the locations of several important agricultural areas used in the paper. Huang-Huai-Hai Plain is composed of Hebei, Beijing, Tianjin, most of Henan, and northern Anhui and Jiangsu. It is the most important grain producing area in China due to its balanced rain and heat and flat terrain. The Loess Plateau includes Shanxi, Ningxia, north of the Qinling Mountains in Shaanxi, and southeastern Gansu. The terrain of this area is complex and diverse, and the ecological environment is fragile.
2.2. Methods
Step 1: Extract monthly NDVI sequence of farmland
This study concerns the MCI of farmland. Hence, the NDVI dataset of farmland was extracted only to eliminate interferences of other land use types. Firstly, the grids with farmland attribute in the dataset of national land use were clipped to build up a mask of farmland. Secondly, the mask of farmland was overlapped and spatial registration with NDVI data from January to December using ArcGIS software, thus extracting the NDVI sequences of farmland.
Step 2: Reconstruct NDVI sequence through the Savitzky-Golay (S-G) filter
The NDVI sequence of farmland had noise caused by atmospheric interferences or other reasons, thus making the NDVI value lower than the true value. Hence, the S-G filter was used for further smoothing and denoising of the NDVI sequence. As a result, high-quality NDVI sequence, which represented the growth trend of crops, was gained. The S-G filter is a convolutional smoothing approach based on the least square method [
40], and it performs the polynomial least square fitting to the adjacent values in a local window. The S-G filter needs two parameters, which are the width of the smoothing window (m) and the degree of the polynomial (d). It requires that the m is shorter than the length of the NDVI sequence and is an odd number, and d is less than m. The larger the m and the smaller the d, the smoother the filtering result, but it also possibly eliminate more real details. Attributes of each farmland grid were the NDVI values from January to December and the sequence length was 12. According to the principle of parameter determination, three filtering windows of (m = 3, d = 2), (m = 5, d = 3) and (m = 5, d = 4) were chosen.
Step 3: Combine the original curve and fitting curve of NDVI
This step was used to maintain high values and decrease abnormal low values. When the original value of the NDVI was higher than the fitting value obtained in Step 2, the original NDVI value was retained. Otherwise, when the original NDVI value was smaller than the fitting value, it was replaced by the fitting value, and the NDVI curve was rebuilt.
Step 4: Calculate the fitting effect coefficient
This step was used to judge the fitting effect between the newly NDVI sequence obtained in Step 3 and the original NDVI sequences. The smaller coefficient indicates the better fitting effect. The calculation formula of the fitting effect coefficient is as follows:
where
is the fitting effect coefficient of S-G filter; n is the length of NDVI sequence (which is 12);
is the serial number of elements in the NDVI sequence;
and
are the fitting value and original value of the NDVI of the element
, respectively;
is the weight of element
;
is the absolute residual error between the original value and fitting value of the NDVI of element
; and
is the maximum value in
. The
of the three filtering windows in Step 2 was compared and that of the (m = 5, d = 4) was the smallest. The coefficient of (m = 5, d = 4) was 0.054, which was 0.014 and 0.013 lower than that of (m = 3, d = 2) and (m = 5, d = 3), respectively. The fitting accuracy of filtering windows of (m = 5, d = 4) increased by about 20% (
Table 1). Therefore, the fitting values of the original NDVI data were performed by filtering windows (m = 5, d = 4).
Step 5: Peak recognition
The first-order differential method was used to recognize the peaks and valleys of the NDVI sequence. In the NDVI sequence of three successive months that first rises and then falls, the middle value was recognized as peak. On the contrary, in the NDVI sequence of three successive months that first falls and then rises, the middle value was recognized as valley. Meanwhile, a statistical analysis on the peak number of each grid was carried out to represent the MCI in each farmland grid.
Step 6: Eliminate interference peaks
In order to reduce the error of the MCI, this paper set up some criteria to remove interference peaks : ① The NDVI value of each peak was higher than 0.4, which was the empirical value of relevant researches [
41]; ② to remove interference peaks, the occurrence time of real peaks was limited from April to October; and ③ the NDVI difference between peak and its adjacent two valleys cannot be smaller than 20% of the difference between the maximum and minimum in the NDVI sequences of 12 months. When one peak could not meet the above three criteria at the same time, it was regarded as an interference peak and deleted.
4. Discussion and Conclusions
Based on the NDVI data, this paper calculated the MCI of farmland in China from 2000 to 2018, and explored the spatio-temporal characteristics of MCI. In addition, the spatial optimization scheme of farmland was put forward according to the gap between MCI and PMCI. The conclusions are drawn as below: from 2000 to 2018, China’s MCI of farmland underwent the fluctuation of rising first, then falling, before rising continuously. These fluctuations were closely associated with the agricultural support policies enforced in different stages and farmers’ reduction in the intensity and utilization of farmland, owing to the low income earned from growing grain. The areas with high MCIs in China were situated in the major grain producing areas, such as Huang-Huai-Hai plain and the southern areas including Guangdong and Guangxi provinces. The proportion of counties with declining MCIs from 2009 to 2018 was lower than that from 2000 to 2009, but the counties with declining MCIs in the later stage were chiefly situated in the major grain producing areas, which exert adverse influence on food security. Compared with the PMCI, the utilization intensity of farmland in northern China was high, while most areas in southern China boasted great potential to increase the MCI. Four different regions and relative optimizing countermeasures were proposed.
In consideration of China’s basic national conditions, i.e., more people and less farmland, small-scale production restricts the increase of the income from agricultural production. Despite that, the overall MCI in China has shown a rising trend since 2000, with still more than 30% of the counties having experienced a downturn. Firstly, China’s floating population was 376 million in 2020, and most of them were rural migrants. Along with the acceleration of China’s industrialization and urbanization, a large number of laborers will still be transferred to cities and towns, and the input of agricultural labor will continue to be decreased. Farmland in Chinese agricultural areas takes the form of collective property rights, and some migrant farmers choose to make use of farmland in an extensive manner for fear of the loss of their rights and interests arising from land transfer. Secondly, government departments have not yet controlled the abandonment of farmland or the change of two seasons to one season according to legal instruments. Following the present development trend, the MCI of some counties will take a downturn trend in the future. It is necessary to attach enormous importance to the phenomenon of decreasing the MCI in principal grain producing areas to avoid food security problems as a result of a large-scale occurrence. Thirdly, non-agricultural farmland is inevitable due to rapid urbanization. In the pursuit of the balance of arable land, low quality “new” arable land reduces the multiple cropping of arable land, which leads to ecological problems.
There is an interactive relationship between the MCI of farmland and the natural environment. In the long term, the high level of climatic variability affects the MCI by affecting the PMCI. With the improvement of farmland irrigation facilities, the water supply limitations noticeably decreased under the irrigated scenario. The growth of the annual mean temperature was identified as the main reason underlying the increase of the PMCI. Furthermore, the area found to be suitable only for single cropping farming decreased, while the area suitable for triple cropping farming increased significantly from the 1960s to the 2000s. The magnitude of change of the PMCI showed a pattern of increase both from northern China to southern China and from western China to eastern China. However, the fluctuations of the MCI calculated in this paper are mainly related to agricultural policies and farmers’ decision-making, and were not always increasing. Furthermore, the spatial pattern of the MCI does not change with longitude and latitude, and the proportion of single cropping farmland increased significantly. Clearly, the spatio-temporal characteristics of the MCI and PMCI variability are not consistent. It can be seen that from a short-term perspective, the MCI variability is mainly affected by farmers’ planting behavior, whose goal is to maximize economic benefits, even at the expense of ecological resources and the environment. Among them, the impact on water has received widespread attention.
Thus far, there is no conclusive conclusion on whether increasing MCI can improve the efficiency of water utilization. Studies in Pakistan [
30] and Brazil [
44] proved that multi-cropping system improved water use efficiency. This situation is likely due to the fact that crop species selection is determined, in large part, based on farmers’ financial interests, but not necessarily on which crop is the most suitable. However, the spatial mismatch between the MCI and PMCI in China is bound to increase the overall water consumption. Southern China possesses favorable water and heat conditions, but low efficiency in terms of the utilization of farmland, while northern China possesses high intensity utilization of farmland. In particular, groundwater has been over-exploited in the Huang-Huai-Hai plain for the development of irrigation agriculture. As a result, serious underground funnels appeared in this area.
Some suggestions were or will be put forward to solve the problem of increased water consumption. Firstly, the Chinese government has begun to enforce measures to close pumping wells in some areas of Hebei, and decrease over-exploitation of groundwater by way of cropping rotation. Secondly, following the concept of the green development of agriculture, China’s agricultural production should optimize the pattern of utilization of farmland in the future, and reinforce the determination of yield by water in areas with high utilization intensity. Meanwhile, the government should increase the utilization intensity of farmland in key development zones. In southern China, it is necessary to issue agricultural support policies, such as comprehensive production subsidies, agricultural socialization services, and financial loans, etc., to encourage and support the adjustment of planting structure or land transfer, and promote the shift from single cropping to multiple cropping, or the rotation of food crops and cash crops. Thirdly, the distortion in the implementation of policy of “pothook of city construction land increase and rural residential land decrease” has caused a decline in quality and implicit decrease in quantity of farmland, respectively. To solve this problem, policies can be adjusted to establish a supplementary mechanism based on farmland productivity. The increase in the MCI of farmland can be used to offset the amount of farmland balance index.
In this paper, by comparing the MCI of farmland calculated based on NDVI data with statistics, we proved that the MCI of farmland, gained by ways of S-G filtering, peak extraction, and peak elimination, is credible and endowed with a reference and popularization value. Grain production in China is mainly carried out by smallholders, and the farmland vegetation is complex, diverse, and irregularly distributes on the surface, particularly in mountainous and hilly areas where there are few large-scale agricultural crops of the same type. Since remote sensing images are made up of regular grids of equal size and influenced by mixed pixels, there will inevitably be certain errors in our results. In future research, perfecting the selection of extraction methods and making use of higher resolution data, in order to raise the accuracy of the MCI of farmland, should be attempted.