Next Article in Journal
Combined Effect of Ultrasound and Microwave Power in Tangerine Juice Processing: Bioactive Compounds, Amino Acids, Minerals, and Pathogens
Next Article in Special Issue
The Impact of Off-Farm Employment on Farmland Production Efficiency: An Empirical Study Based in Jiangsu Province, China
Previous Article in Journal
A Multimethod Study on Kitchen Hygiene, Consumer Knowledge and Food Handling Practices at Home
Previous Article in Special Issue
On-Farm Spatial Composition, Management Practices and Estimated Productivity of Urban Farms in the San Francisco Bay Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study on Spatial-Temporal Differentiation and Influencing Factors of Agricultural Water Footprint in the Main Grain-Producing Areas in China

1
School of Management and Economics, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
2
School of Public Administration, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
*
Author to whom correspondence should be addressed.
Processes 2022, 10(10), 2105; https://doi.org/10.3390/pr10102105
Submission received: 27 September 2022 / Revised: 13 October 2022 / Accepted: 13 October 2022 / Published: 17 October 2022
(This article belongs to the Special Issue Processes in Urban Farming and Food Security)

Abstract

:
It is an urgent scientific issue to explore the spatial and temporal differentiation and impact indicators of the agricultural water footprint in major grain-producing areas. Therefore, this study tries to use the water footprint theory to implement top-down calculation of the agricultural water footprint in major grain-producing regions from 2000 to 2019 and investigate the various impacts on the agricultural water footprint under the influence of spatial-temporal effects using spatial autocorrelation and the spatial Dubin model. The results indicate that from 2000 to 2019, the overall agricultural water footprint of China showed a fluctuating downward trend in an inverted N shape and demonstrated high–high and low–low spatial aggregation characteristics. There are notable characteristics, including high spatial dependence, spatial barriers, and path locking of the agricultural water footprint in most provinces and regions of the main grain-producing areas. Policy factors, water-saving technologies, social development, economic development, and industrial structure adjustment are all significantly and negatively correlated with the increase in the agricultural water footprint, while agricultural production and natural factors have a significant positive relationship with the agricultural water footprint. The spatial spillover effect of water-saving technologies, industrial restructuring, agricultural production, and natural factors is powerful. Therefore, a rationally optimized industrial structure, strengthened regional linkage of water resources management and control, and the promotion of efficient water infrastructure technology are important ways to inhibit the agricultural water footprint.

1. Introduction

Water and food security are two of the vital global challenges of the future [1]. As a large, populous, and agricultural country, China is confronted with more severe pressures on water and food security. Future population growth, intensified competition for industry-used water, and climate change will all challenge water use in agriculture [2]. The 13 major grain-producing provinces, which are mainly responsible for China’s agricultural production, supplied up to 78.89% of the country’s grain production in 2019 but used less than 40% of water resources [3]. Under the current Chinese grain supply and demand structure, the main grain-producing areas have significant advantages in grain production in terms of water and soil resources, and they are bearing increasing pressure on agricultural water usage while exporting massive internal water resources to the external areas. However, considering the current accelerating frequency of extreme weather caused by climate change and the long-term ambiguity caused by COVID-19 since early 2020, agricultural ecological environment problems in China such as the degradation of arable land quality, agricultural non-point source pollution, soil erosion, and groundwater over-extraction have become increasingly prominent [4]. Agricultural water resources flow in regional agricultural product trade in the form of embedded water in the production, consumption, and trade process [5]. The water footprint theory can effectively quantify the real demand for water resources in regional agriculture, and it transforms the emphasis on traditional water resources from water resource supply management to product demand management. It broadens the traditional water resource evaluation system and provides a decision basis for the rational and efficient utilization of water resources [6]. In view of the great threat posed by a series of realistic dilemmas of agricultural water use to China’s food security, it is necessary to use the water footprint theory to explore the spatial-temporal differences of agricultural water footprint and its influencing factors in China’s major grain-producing areas.
In relevant studies on agricultural water footprint assessment methods, foreign assessment methods are mainly based on virtual water theory and life cycle theory, and the spatial scale of assessment contains global, national, provincial, and municipal levels. Abdullah et al. [7] argue that spatial variations of the water footprint of agricultural products can be employed to mitigate regional and global problems concerning water scarcity. After confirming the virtual water content of agricultural products, for the first time, Chapagain and Hoekstra [8] systematically calculated the agricultural water footprint from 1997 to 2001 at the global level, and further explored virtual water flows through commodity trade between countries. Then, Boulay et al. [9] proposed an agricultural product water footprint calculation method based on the life cycle methodology and discussed the impact of land use form changes on blue water resources. Rodriguez et al. [10] found that the gray water footprint accounted for a higher proportion when they analyzed the potato water footprint in the pam-pean region of Argentina, so they suggested improving the efficiency of fertilizer application and regulating the regional agricultural production water footprint with more professional technical support and stricter groundwater control measures. Abdullah et al. [11] demonstrated the value of the water footprint accounting framework for such river basins with the case study of the Upper Euphrates River basin. Zhou et al. [12] believed that reforming rice cropping systems could mitigate the water footprint and enhance grain yield in central China. Tamea et al. [13] established a synthesis database for country-oriented water footprints based on 370 agricultural products to support the national and global agricultural water assessment. Deihimfard [14] quantified the water footprint of agroecological systems in 17 Iranian provinces located in arid or semi-arid regions. Hossain [15] calculated the water footprints of local fruits and vegetables in Victoria according to the Australian Water Footprint Assessment Manual framework.
In research on the influencing factors of the agricultural water footprint, many scholars tend to consider the influence of natural environmental factors and production technology factors on the agricultural water footprint. Bocchiola [16] explored the effects of natural factors such as precipitation amount, soil conditions, humidity, and wind speed on the water footprint of crops. Cao et al. [17] found that irrigation patterns and social factors are the main factors affecting the agricultural water footprint. In terms of production technology, Ababaei and Etedali [18] found that the total power of farm machinery and the application of agricultural materials such as pesticides and fertilizers have an outstanding influence on the blue water footprint of wheat, and technology has a notable positive influence on the green water footprint of wheat. The study of Rao et al. [19] indicated that improper irrigation, low fertilization, and improper farm water management resulted in a high water footprint of crops. Munro et al. [20] reduced the water footprint in crop production by applying soil mulching and drip irrigation and analyzed the influencing factors of spatial differentiation in crops, planting patterns, and soil and field management factors.
At present, studies on the assessment methods and influencing factors of the agricultural water footprint are of great value, but there are certain shortcomings, as follows: Most of the previous research was conducted at national or provincial levels, but research setting the main grain-producing areas as the research object is rare; few studies have concentrated on the spatial and temporal differentiation of the agricultural water footprint, especially on the spatial distribution and significant congregation of the agricultural water footprint; study on the influencing factors of the agricultural water footprint failed to take the influence of spatial and temporal differences into consideration, and it is even rarer to adopt the spatial and temporal effects into the influencing factors analysis model. Therefore, based on the assessment of the agricultural water footprint in the main grain-producing areas, this paper utilizes the spatial autocorrelation method to explore the spatial and temporal differentiation characteristics, adopts the spatial and temporal effects in the influencing factor analysis model, and applies the spatial Durbin model to explore the driving impacts of each influencing factor on the agricultural water footprint.
The other parts of this paper are structured as follows: Section 2 contains the data methodology and data sources. Section 3 illustrates the empirical results and analysis. Section 4 demonstrates the conclusions and implications.

2. Research Methods and Data

2.1. Theoretical Framework

The measurement of the agricultural water footprint needs to consider the real consumption of water resources of regional agriculture and the cross-regional trade of virtual water, including the internal water footprint and the external water footprint. This paper explores the spatial-temporal differentiation characteristics of the agricultural water footprint by applying the SDM model and discusses the spatial effect as a critical influencing factor. The agricultural water footprint, as an indicator used to evaluate the effect of regional water resource utilization, is also affected by regional policies, water-saving technologies, production levels, social development, economic development, natural factors, etc. Among the above factors, the adjustment of policy factors, water-saving technologies, social development, economic development, and other factors could increase the efficiency with which water is used, thereby restraining the increase in the agricultural water footprint. This research analyzes the spatial-temporal differentiation of the agricultural water footprint in the main grain-producing areas in China based on the assessment of the agricultural water footprint, adopts the spatial factors in the influencing factor analysis model, and explores the influencing indicators of the agricultural water footprint in order to provide effective suggestions and solutions for breaking the spatial path dependence of the regional agricultural water footprint and exploring approaches for the regionally linked water resources management in the future. The theoretical framework of this paper is shown in Figure 1.

2.2. Study Region

The main grain-producing area refers to the exclusive economic zone where natural conditions such as climate, geography, and soil, and economic conditions such as human resources, technology, and materials are more conducive to the cultivation of grain crops with high proportion, large production, and output of food crops. In 2003, China’s Ministry of Finance designated 13 provinces as major grain-producing areas covering Northeast China, Huang-Huai-Hai, and the Yangtze river, which consist of Inner Mongolia, Jilin Province, Liaoning Province, Heilongjiang Province, Hebei Province, Jiangsu Province, Anhui Province, Shandong Province, Henan Province, Jiangxi Province, Hubei Province, Hunan Province, and Sichuan Province, the distribution of which is shown in Figure 2.

2.3. Agricultural Water Footprint

The agricultural water footprint is the total of the internal and external water footprint, in which the internal water footprint is the local agricultural production water minus the virtual water of local agricultural products, while the external water footprint is the entire virtual water of agricultural products from other places consumed by local residents [21]. The calculation formulas are listed as follows:
A W F j t = I A W F j t + E A W F j t = A W F P j t A W F O j t + A W F I j t = A W F P j t A W F N j t
In the above formula, t stands for the year, j stands for different regions, A W F j t stands for the agricultural water footprint (100 million m3), I A W F j t stands for the internal agricultural water footprint (100 million m3), and E A W F j t stands for the external agricultural water footprint (100 million m3). A W F P j t stands for the total amount of water required to produce agricultural products (100 million m3), namely the total crop production water footprint. A W F O j t stands for the virtual water output of agricultural products (100 million m3), A W F I j t stands for the agricultural virtual water input (100 million m3), and A W F N j t stands for the agricultural virtual water net output.

2.3.1. Agricultural Production Water Footprint

The agricultural production water footprint refers to the water consumption for producing a certain crop in a region within a specific time. The crop water footprint is the actual utilization sum of water resources during the growing process of a certain crop in a specific area, which can be distributed into the gray water footprint, blue water footprint and green water footprint [22]. The blue water and green water footprints are used for crop growth consumption of irrigation water and effective precipitation via the growth of crops, while the gray water footprint is the amount of freshwater consumed by diluting the pollutants produced in the agricultural production to meet the environmental water quality standards [23].
A W F P ij t = A W F P i j t g r e e n + A W F P i j t b l u e + A W F P i j t g r e y
A W F P i j t g r e e n = P j e t × S i j t γ i j t
P m 250 mm , P e m = P m ( 125 0.2 P m ) / 125 P m > 250 mm , P e m = 125 + 0.1 P m
A W F P i j t b l u e = I R j t × S i j t G = I R j t × S j I R t S i j t S j t S
A W F P i j t g r e y = α × A R i j t C max C nat
In the above formula, S i j t stands for the sown area (1000 hectares) of crop i in period i and area j , γ i j t is the multi-cropping index, P j e t stands for the effective precipitation amount in period t and area j , and P m and P e m stand for the monthly precipitation amount and monthly effective precipitation amount (mm). I R j t stands for the water consumption for unit farmland irrigation in period t and area j (m3/ha), S i j G t stands for the irrigated areas of crops (1000 hectares), S j I R t stands for effective irrigated areas (1000 hectares), S j S t stands for the total crops sown areas (1000 hectares), and S i j t stands for the sown area of crop i (1000 hectares). α stands for the nitrogen fertilizer leaching rate, namely the proportion of water pollution caused by nitrogen fertilizer, valued as 10%, A R i j t stands for the pure nitrogen element amount in nitrogen fertilizer consumption (kg) of crop i in each region, C max is valued as 0.01 kg/m3, and C n a t is valued as 0.

2.3.2. Agricultural Virtual Water Net Exports

The agricultural virtual water net exports represent the net output of virtual water generated when regional agricultural products are exported and imported at a certain time [24]. The calculation formula is listed as follows:
A W F N i j t = A W F O i j t A W F I i j t = W F P i j t ( O i j t I i j t ) = A W F P i j t Y i j t ( O i j t I i j t )
In the above formula, A W F N i j t stands for the net output of agricultural virtual water, A W F O i j t is the agricultural virtual water output (100 million m3), and A W F I i j t is the agricultural virtual water input (100 million m3). W F P i j t is the unit crop production water footprint (m3/kg), Y i j t is the crop yield (kg), and O i j t and I i j t are the output and input (kg) of regional agricultural products, respectively.

2.4. Spatial Analysis Methods

2.4.1. Spatial Autocorrelation Method

The spatial autocorrelation method involves both global and local aspects, which are measured by the global Moran’s I index and local Moran’s I index separately [25]. The global Moran’s I index generally reflects the spatial correlation characteristics of the whole region, and the local Moran’s I index is used for evaluating the spatial aggregation characteristics in local areas.
The global spatial autocorrelation analysis can measure whether there is a spatial correlation at the overall level of the spatial unit. The calculation formula is computed as:
I = i = 1 n j = 1 n w i j ( x i x ¯ ) ( x j x ¯ ) S 2 i = 1 n j = 1 n w i j
In the above formula, I is the global Moran’s I index, W i j is the spatial weight matrix, X i and X j  are attribute values of each spatial unit, x ¯ is the mean value, S 2 is the sample variance, and n stands for the number of spatial units.
The local spatial autocorrelation method analyzes the aggregation characteristics of local space and specific attribute space. The calculation formulas are as follows:
I = x i x _ S 2 j i n w i j ( x j x _ )
S 2 = 1 n ( x i x ¯ ) 2
When the local Moran’s I index is greater than 0, it is believed that the high (low) valued regions and the high (low) valued adjacent regions are clustered. When the local Moran’s I index is less than 0, it is considered that the high (low) valued regions and the low (high) valued adjacent regions are clustered.

2.4.2. Spatial Dubin Model

The spatial econometric model considering both temporal and spatial effects can effectively analyze the spatial influence of each factor on the explained variable [26]. This paper selects the spatial Dubin model (SDM) with fixed space and time to conduct spatial econometric analysis on the factors of the agricultural water footprint in the main grain-producing regions. The model is constructed as follows:
ln Y i t = ρ j = 1 n W i j ln Y j t + β ln X i t + φ j = 1 n W i j ln X j t + μ i + v t + ε i t
In the above formula, Y i t stands for the explained variable, namely the agricultural water footprint value, X i t stands for the explanatory variable, namely the influencing factors that affect the agricultural water footprint, W i t stands for the spatial weight matrix, and μ i and v t stand for spatially fixed effects and the time-fixed effects.

2.5. Variable Selection and Data Sources

Referring to existing research, this study incorporates seven variables, including policy factors, agricultural water-saving technologies, agricultural production, economic and social development, and natural factors, into the influencing factor system, which is shown in Table 1. Data in this paper mainly come from the 2000–2019 related statistical yearbook, and statistical yearbooks of various provinces and regions. The precipitation data come from the related platform. The irrigation water data come from the related bulletins of various provinces and regions. The nitrogen fertilizer application amount for various crops per unit area comes from the National Agricultural Product Cost and Benefits Data Compilation. This paper adopts the linear prediction method to make up for the missing data during the sample period. The price-related indicators are deflated according to indicators in the year 2000. Similarly, missing data are also filled in using the linear prediction method.

3. Empirical Results and Analysis

3.1. Spatial-Temporal Differentiation Analysis of Agricultural Water Footprint in the Main Grain-Producing Areas

3.1.1. Calculation of Unit Production Water Footprint of Main Grain Crops in Main Grain-Producing Areas

The unit water footprint and its itemized water footprint values of major grain crops in the major grain-producing areas in the past 20 years were calculated, as shown in Table 2. By comparing the average unit water footprint of each crop, the water consumption types of each crop can be accurately classified. According to the unit agricultural production water footprint, the water consumption types of each crop are further divided into the following three categories: When WFP > 2, the unit crop production water footprint of crops is categorized as high water consumption; when 1 < WFP < 2, the unit crop production water footprint of crops is categorized as medium water consumption; when WFP < 1, the unit crop production water footprint of crops is categorized as low water consumption [27]. According to the average value of the per unit crop production water footprint, rice and corn are medium-water-consumption crops, while beans and wheat are high-water-consumption crops. It can be found that in the composition of the water footprint per unit production, the per unit green water footprint of all food crops was higher than the blue water footprint, indicating that precipitation restricted food production.
It can be seen from Figure 3 that both the grain yield and production value in major grain-producing areas increased year by year from 2000 to 2019, which was mainly due to the increase in the yield of rice and corn, especially the high economic value of corn, which led to an increase in the agricultural production value in the main grain-producing areas.

3.1.2. Time Series Variation Characteristics of the Agricultural Water Footprint in the Main Grain-Producing Areas

When calculating the agricultural water footprint in 13 major grain-producing regions from 2000 to 2019, this paper compares the production and consumption data of grain crops in each main grain-producing region. It is assumed that when the consumption of grain crops is less than the output of grain crops, the excess is exported to outside areas. Based on this, the agricultural water footprint values calculated in this research are shown in Figure 4 and Figure 5.
Figure 6 reflects the time-series variation of the agricultural water footprint in the main grain-producing regions from 2000 to 2019. Taking the main grain-producing regions as a whole, the agricultural water footprint showed a fluctuating downward trend in an inverted N shape. From 2000 to 2005, the overall agricultural water footprints in the main grain-producing regions were comparatively high, and the regional differences were obvious. Agriculture water consumption was high. From 2006 to 2010, the high-water-consumption areas gradually moved closer to the low-water-consumption areas, and the agricultural water footprint differences between regions were reduced. The agricultural water footprint rebounded slightly from 2011 to 2016, and it gradually stabilized below 15 billion cubic meters after 2016. Compared with the beginning of the period, agricultural water consumption decreased to varying degrees.

3.1.3. Spatial Distribution Differences of Agricultural Water Footprint in the Main Grain-Producing Regions

Based on the measurement results of the global spatial autocorrelation Moran’s I index, Table 3 shows the global Moran’s I index significance level of the agricultural water footprint in the main grain-producing regions from 2000 to 2019 is between 0.000 and 0.016 and all below 1%, excluding the year 2019. It means that there is a significant positive spatial autocorrelation for the agricultural water footprint in the main grain-producing regions in most years, and there is an obvious spatial agglomeration effect from the overall perspective. From 2000 to 2019, the global Moran’s I index of the main grain-producing areas showed a slow wavelike decrease, indicating that the spatial autocorrelation of the agricultural water footprint in the main grain-producing areas has weakened since 2000, but the overall spatial aggregation effect has become stronger.
Regarding local spatial autocorrelation, the four quadrants of the local Moran’s I scatter diagram stand for the four modes of spatial aggregation, namely high–high aggregation, low–high aggregation, low–low aggregation, and high–low aggregation (Figure 7).
On the whole, many regions are scattered in the first and third quadrants, displaying a strong spatial agglomeration effect as a whole. With the change in time, the spatial aggregation characteristics of the main grain-producing areas showed high–high aggregation and low–low aggregation, and the trend of polarization was obvious. Provinces and regions with a low agricultural water footprint, on the one hand, were similar in agricultural production and consumption. On the other hand, they transformed into low-value aggregation areas through regional cooperation, water-saving technologies, and resource management policies, and generated strong radiation and driving effects on neighboring areas. The main grain-producing provinces and regions with a high agricultural water footprint were usually high–high aggregation areas because of their high population density, high pressure on agricultural production and water resources development, and the impact of the high agricultural water footprint in adjacent areas.
In terms of various periods, during the 10th Five-Year Plan period, high–high agglomeration areas include Shandong, Jiangsu, Henan, Anhui, and Sichuan, while low–low agglomeration areas include Inner Mongolia, Jilin, Liaoning, and Heilongjiang, low–high agglomeration areas include Jiangxi and Hubei, and high–low agglomeration areas include Hebei and Hunan. During the 11th Five-Year Plan period, the spatial aggregation characteristics of the agricultural water footprint in the main grain-producing areas were strengthened, revealing distinct polarization. Except for Inner Mongolia and Northeast China, which are low–low aggregation areas, all the other provinces are high–high aggregation areas. During the 12th Five-Year Plan period, Jiangxi and Hubei moved to the second and fourth quadrants, respectively. During the 13th Five-Year Plan period, Hubei and Hunan moved to the second and fourth quadrants, respectively, compared with the previous period. The Northeast region and Huang-Huai-Hai region demonstrate more agglomeration effects, while the Yangtze river basin shows more discrete effects. Judging from the distribution of the main grain-producing provinces in various agglomeration areas, most provinces are quite stable, which means that there are notable characteristics including high spatial dependence, spatial barriers, and path locking of the agricultural water footprint in most provinces and regions of the main grain-producing areas.

3.2. Analysis of Influencing Factors of Agricultural Water Footprint in the Main Grain-Producing Areas

3.2.1. Model Test

To ensure the accuracy of the influencing factors model for the agricultural water footprint in the main grain-producing regions and increase the credibility of the estimated results, this study carried out the spatial LM test, Hausman random effect, fixed-effect test, Wald test, and likelihood ratio LR test. Results show that the Spatial Dubin Model with fixed space and time should be adopted. The test results are listed in Table 4.

3.2.2. Estimation Results and Analysis

Spatial Dubin model estimation results of the whole main grain-producing areas show that policy factors, water-saving technologies, social development, economic development, and industrial structure are significantly negatively correlated with the agricultural water footprint, while agricultural production and natural factors are markedly positively correlated with agricultural water footprint.
The policy factor that regards the proportion of irrigation and water conservancy expenditure in the total water conservancy investment as a proxy variable passes the 5% significance level test, and the coefficient is significantly negative. It means that the construction of agricultural irrigation and water resource facilities led by the central and local governments can restrain the agricultural water footprint increases.
The coefficient of water-saving technologies is −0.043, and it passes the 10% significance level test, confirming that water-saving irrigation technology can effectively restrain the agricultural water footprint of the main grain-producing regions and play a positive role in water sustainable utilization.
The coefficient of agricultural production, although under 1%, is markedly positive, meaning that the development of the agricultural production level with the total power of agricultural machinery as a proxy variable will lead to increases in the agricultural water footprint. The modernization of agricultural machinery and their increasing utilization rate will improve crop yields, thereby consuming more water resources.
The coefficient of social development is 0.215, and it passes the 1% significance level test. It means that with the acceleration of the urbanization process and the growth of the city’s population, on one hand, agricultural water will be cut down. On the other hand, the change in the consumption structure of residents in the city and rural areas results in the need for higher-quality agricultural products as well as changes in the agricultural water footprint accordingly.
The coefficient of natural factors with the amount of precipitation as a proxy variable is distinctly positive, and it passes the 1% significance level test. According to the regression results, the annual precipitation amount has a significantly positive effect on the change in the agricultural water footprint. It means that areas with more abundant precipitation consume more green water in agricultural production and have a higher agricultural water footprint.
The coefficient of economic development is −0.424, and it passes the 1% significance level test. It shows that because of the economic growth, the financial, technical, and resource support brought by regional economic development to agricultural production, especially the application of water-saving technologies and the improvement of irrigation facilities, could all effectively restrain the increase in the agricultural water footprint.
The coefficient of industrial structure is −0.211 and it passes the 1% significance level test. Industrial restructuring is the main factor that restrains an increase in the agricultural water footprint. The increase in agricultural value-added production, the competitiveness enhancement of the agricultural industry in both domestic and foreign markets, and the sustainable development of agriculture are the main manifestations of agricultural industrial restructuring. The high-quality development of the agricultural economy in the main grain-producing regions will help to restrain the increase in the agricultural water footprint.
In order to further explore the regional differences of various influencing factors on the agricultural water footprint in the main grain-producing regions, this study divides the main grain-producing areas into the northern region and the southern region (Table 5). According to the spatial Dubin model estimation results at the subregional level, the coefficients of policy factors in the northern and southern regions are both significantly negative, which are −0.038 and −0.030, respectively. It means that agricultural irrigation and water conservancy expenditures pose deeper influences on restraining the agricultural water footprint in the northern region than in the southern region. The coefficients of water-saving technologies in the northern and southern regions both pass the 5% significance level test, and the absolute coefficient value of the southern part of the region is significantly larger than that in the northern region. It means that the driving effect of water-saving technologies in the northern region is significantly lower than that in the southern region. The coefficients of social development and natural factors in both northern and southern regions are consistent overall, and both of them pass the significance test. Different from the whole situation, the coefficients of agricultural production and industrial restructuring in the northern region and the coefficient of economic development in the southern region failed to pass the significance level test.

3.2.3. Decomposition Effect Analysis

Judging from the decomposition results for the whole main grain-producing area, the spatial spillover effects of water-saving technologies, agricultural production, natural factors, and industrial structure are strong. The details are listed as follows.
The direct effect coefficient of policy factors is negative, and it passes the 5% significance level test, but the indirect effect coefficient fails the test. This means that the increase in agricultural irrigation and water conservancy expenditure mainly restrains the increase in the local agricultural water footprint. The self-interested nature of local fiscal policies to support agriculture determines that all kinds of agriculture-supporting policies could only affect local agricultural production, and their spatial spillover effects are not strong. Although the capital expenditure for the construction of irrigation facilities can effectively reduce the local agricultural water footprint to some extent, there is a large tendency to form path dependence and locking features in the region in the long run.
The direct and indirect effect coefficients of water-saving technologies are negative, and they pass the 10% and 5% significance level tests. It means that the development of water-saving technologies is not only beneficial to suppress the increase in the local agricultural water footprint, but also to inhibit the increase in the agricultural water footprint in adjacent areas through the negative spillover. Due to the positive externality of the technology, the local advanced irrigation and water-saving technologies spread to the outside and the advanced water resources management knowledge is diffused outward. Therefore, advanced water-saving technologies and knowledge will help to control the global agricultural water footprint.
The direct and indirect effect coefficients of agricultural production are both positive, and they pass the 1% and 5% significance level tests, respectively. This means that the extensive application of agricultural machinery will not only increase the local agricultural water footprint but also increase the agricultural water footprint in neighboring sites through positive spillover. The development of agricultural machinery technology in global regions has promoted the widespread application of agricultural machinery and equipment. In particular, the cross-regional services of agricultural machinery and equipment promote the mutual development of agricultural machinery in local areas, thereby increasing crop yields in local areas.
The direct effect coefficient of social development is −0.211 and it passes the 1% significance level test, while the indirect effect coefficient is not significant. It means that urbanization has the greatest influence on local agricultural water usage. Due to the similarity of household consumption between the adjacent regions in major grain-producing regions, changes in the consumption patterns and structures of residents have the most significant impacts on local agricultural water footprint.
The direct and indirect effect coefficients of natural factors are significantly positive under the 1% and 5% levels, respectively. It means that due to the similar natural climatic characteristics between the adjacent areas in the main grain-producing areas, the amount of precipitation has consistently positive effects on the agricultural water footprint of local and adjacent areas. The increase in local precipitation produces abundant surface water and groundwater resources, which cause positive spillovers through the water circulation system, and relieve the agricultural water pressure in adjacent regions.
The direct effect coefficient of economic development is significantly negative, but the indirect effect coefficient fails the significance level test. It means that economic development can significantly inhibit the local agricultural water footprint, but its spatial spillover effect is not strong enough to generate a negative effect on the agricultural water footprint in adjacent areas. The regional linkage of economic development does not work effectively, which works against the decrease in the agricultural water footprint in the main grain-producing regions.
The direct and indirect effect coefficients of the structure of the industry are both negative, and they pass the 1% and 5% significance level tests, respectively. It means that industrial restructuring can significantly inhibit the agricultural water footprint of local and adjacent areas.
Based on the regional decomposition results of the main grain-producing regions (Table 6), the direct and indirect effect coefficients of water-saving technologies and economic development in the northern region are both positive, and they all pass the 5% significance level test. It means that the growth of water-saving irrigation areas and economic development in the northern region are beneficial to restraining the increase in the local agricultural water footprint, and in the meantime, suppress the increase in the agricultural water footprint in adjacent areas through the negative spatial spillover effect. The direct effect coefficients of policy factors, social development, and natural factors all pass the 5% significance level test, which are −0.040, −0.329, and 0.303, respectively, while the indirect effect coefficients failed the test. It means that the increase in irrigation water conservancy expenditures and the improvement of the urbanization level in the northern region have outstanding effects on restraining the increase in the local agricultural water footprint. The direct and indirect effect coefficients of agricultural production in the southern region are both positive, and both of them pass the 1% significance level test. It means that the improvement of local agricultural mechanization in the southern region not only elevates the local agricultural water consumption but also promotes the increase in the agricultural water footprint in the adjacent areas through the positive spillover effect. The direct effect coefficients of policy factors, water-saving technologies, social development, natural factors, and industrial structure pass the 10% significance level test, and their indirect effect coefficients failed the test. It means that the spillover effect of these influencing indicators in the southern region is not strong enough.

4. Conclusions and Implications

4.1. Main Conclusions

This study refers to agricultural data of the main grain-producing areas from 2000 to 2019, uses the spatial autocorrelation analysis to explore the spatial and temporal variation features of the agricultural water footprint based on regional agricultural water footprint assessment, and finally, builds the spatial Dubin model to analyze the influencing factors of the agricultural water footprint in the major grain-producing regions. The main conclusions are listed as follows:
(1)
Judging from the time series features of the agricultural water footprint in the major grain-producing regions, it shows a fluctuating downward trend in an inverted N shape. From 2000 to 2019, the agricultural water footprint in the main grain-producing regions was polarized. The internal differences narrowed from an overall aspect, and the low-value provinces and regions were developing rapidly.
(2)
Based on the spatial distribution differences of the agricultural water footprint in the major grain-producing regions, the agricultural water footprint is significantly positively spatially autocorrelated from the whole, and it shows high–high and low–low spatial aggregation characteristics in local areas, with a distinct polarization trend. For example, the northeast region presents low–low aggregation while the Huang-Huai-Hai region presents high–high aggregation. From 2000 to 2019, there was strong spatial dependence, spatial barriers, and path-locking characteristics for the agricultural water footprint in most areas of the main grain-producing regions.
(3)
On basis of the influencing factors of the agricultural water footprint in the main grain-producing areas, policy factors, water-saving technologies, social development, economic development, and industrial restructuring could dramatically restrain the increase in the agricultural water footprint. The negative spatial spillover effect of water-saving technologies and industrial structure is strong, and the positive spatial spillover effect of agricultural production and natural factors is powerful, both of which could significantly affect the agricultural water footprint in the adjacent areas.

4.2. Policy implications

(1)
Promoting the development of the agricultural economy with a rational and optimized industrial structure. Industrial restructuring with the proportion of agricultural output value to regional GDP as a proxy variable is helpful to inhibit the agricultural water footprint. It is necessary to control water consumption from the source of agricultural production to realize the sustainable use of agricultural water resources. In the aspect of industrial restructuring, the main grain-producing areas must adjust their internal agricultural cropping structure, decreasing the production of crops with higher water consumption and lower economic value. They also need to adjust the external agricultural industrial structure and vigorously promote other agricultural industries with strong value-added capabilities.
(2)
Strengthening the management and control of water resources via a regional-linked management mechanism. The agricultural water footprint in the main grain-producing areas presents high–high and low–low spatial aggregation characteristics, and there is stable spatial dependence for the agricultural water footprint in most provinces. All major grain-producing provinces and regions should improve the agricultural water footprint assessment system, establish regionally linked agricultural water resources management mechanisms, and reform significant methods for agricultural water resources management in order to break the spatial barriers and path locking of the agricultural water footprint, drive the reduction in the agricultural water footprint from low-value areas to adjacent areas, and eventually achieve low-value aggregation of the agricultural water footprint in global areas.
(3)
Popularizing water-saving irrigation with complete and efficient water conservancy facilities and technologies. Increasing the proportion of water-saving irrigation areas in arable land, popularizing the application of water-saving technologies, and increasing expenditure on irrigation and water conservancy facilities are all beneficial to curb the agricultural water footprint in major grain-producing regions. The major grain production provinces and regions should provide financial support for the construction of farmland water-saving irrigation facilities and the promotion of high-efficiency water-saving technologies by making full use of the agricultural water-saving project funds provided by the government, and actively guiding the non-governmental funds at the same time. They must enhance the construction of inter-regional and cross-regional farmland water conservancy facilities to establish a “powerful framework” for the development of water-saving agriculture. Meanwhile, they need to positively promote new and efficient water-saving irrigation technologies to “soften the veins” for the development of water-saving agriculture.
In the future, major grain-producing areas need to actively explore regional linkage of water resources management policies with the agricultural water footprint theory, scientifically regulate agricultural water footprint, break the spatial dependence of regional agricultural water footprint, and achieve “low-low” aggregation in the whole region. At the same time, major grain-producing areas should pay greater attention to the R&D, promotion, and utilization of efficient water-saving technology in the future and accelerate the development of water-saving agriculture to solve the water resource dilemma.

Author Contributions

Conceptualization, Y.W. and Z.S.; methodology, Y.W.; software, Y.W. and Q.Z.; validation, Z.S. and Q.Z.; formal analysis, Y.W.; investigation, Y.W. and Q.Z.; resources, Z.S.; data curation, Z.S. and Q.Z.; writing–original draft preparation, Y.W.; writing–review and editing, Y.W.; visualization, Q.Z.; supervision, Y.W.; project administration, Z.S; funding acquisition, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Soft Science Research Project of Henan Province, “Study on generation mechanism and realization path of water-carbon-production nexus of Henan Province” (222400410160); Philosophy and Social Science Planning Annual Program of Henan Province, “Research on value chain optimization of energy-intensive manufacturing industries in Henan Province under the dual goals of carbon reduction and development” (2022BJJ065).

Data Availability Statement

Data available in a publicly accessible repository that does not issue DOIs. Publicly available datasets were analyzed in this study. This data can be found here: https://www.epsnet.com.cn/index.html#/Index.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Orimoloye, I.R.; Belle, J.A.; Olusola, A.O.; Busayo, E.T.; Ololade, O.O. Spatial assessment of drought disasters, vulnerability, severity and water shortages: A potential drought disaster mitigation strategy. Nat. Hazards 2020, 105, 2735–2754. [Google Scholar] [CrossRef]
  2. Tian, J.; Liu, J.; Wang, J.; Li, C.; Nie, H.; Yu, F. Trend analysis of temperature and precipitation extremes in major grain producing area of China: Trends of climate extremes in major grain production area of China. Int. J. Climatol. 2016, 37, 672–687. [Google Scholar] [CrossRef]
  3. Song, C.-X.; Liu, R.-F.; Oxley, L.; Ma, H.-Y. Do farmers care about climate change? Evidence from five major grain producing areas of China. J. Integr. Agric. 2019, 18, 1402–1414. [Google Scholar] [CrossRef]
  4. Li, M.; Xu, W.; Zhu, T. Agricultural water allocation under uncertainty: Redistribution of water shortage risk. Am. J. Agric. Econ. 2018, 101, 134–153. [Google Scholar] [CrossRef]
  5. Murphy, E.; De Boer, I.J.; Van Middelaar, C.E.; Holden, N.M.; Shalloo, L.; Curran, T.P.; Upton, J. Water footprinting of dairy farming in Ireland. J. Clean. Prod. 2017, 140, 547–555. [Google Scholar] [CrossRef] [Green Version]
  6. Wichelns, D. Advising Morocco: Adopting recommendations of a water footprint assessment would increase risk and impair food security for the country and its farmers. Water Int. 2018, 43, 762–784. [Google Scholar] [CrossRef]
  7. Abdullah, M.; Zehra, A. Spatial analysis of blue and green water footprints of agricultural crop patterns: Turkey. Water Manag. 2021, 174, 291–308. [Google Scholar]
  8. Chapagain, A.K.; Hoekstra, A.Y. Water footprints of nations. J. Bank. Financ. 2004, 27, 1427–1453. [Google Scholar]
  9. Boulay, A.-M.; Bouchard, C.; Bulle, C.; Deschênes, L.; Margni, M. Categorizing water for LCA inventory. Int. J. Life Cycle Assess. 2011, 16, 639–651. [Google Scholar] [CrossRef]
  10. Rodriguez, C.I.; Galarreta de Ruiz, V.A.; Kruse, E.E. Analysis of water footprint of potato production in the pam-pean region of Argentina. J. Clean. Prod. 2015, 90, 91–96. [Google Scholar] [CrossRef]
  11. Abdullah, M.; Emre, I.; Ertug, E. Water resources management of large hydrological basins in semi-arid regions: Spatial and temporal variability of water footprint of the Upper Euphrates River basin. Sci. Total Environ. 2022, 846, 157396. [Google Scholar]
  12. Zhou, Y.; Liu, K.; Harrison, M.T.; Fahad, S.; Gong, S.; Zhu, B.; Liu, Z. Shifting rice cropping systems mitigates ecological footprints and enhances grain yield in central China. Front. Plant Sci. 2022, 13, 895402. [Google Scholar] [CrossRef] [PubMed]
  13. Tamea, S.; Tuninetti, M.; Soligno, I.; Laio, F. Virtual water trade and water footprint of agricultural goods: The 1961–2016 CWASI database. Earth Syst. Sci. Data 2021, 13, 2025–2051. [Google Scholar] [CrossRef]
  14. Deihimfard, R.; Rahimi-Moghaddam, S.; Collins, B.; Azizi, K. Future climate change could reduce irrigated and rainfed wheat water footprint in arid environments. Sci. Total Environ. 2021, 807, 150991. [Google Scholar] [CrossRef] [PubMed]
  15. Hossain, I.; Imteaz, M.A.; Khastagir, A. Water footprint: Applying the water footprint assessment method to Australian agriculture. J. Sci. Food Agric. 2020, 101, 4090–4098. [Google Scholar] [CrossRef]
  16. Bocchiola, D. Impact of potential climate change on crop yield and water footprint of rice in the Po valley of Italy. Agric. Syst. 2015, 139, 223–237. [Google Scholar] [CrossRef]
  17. Cao, X.; Huang, X.; Huang, H.; Liu, J.; Guo, X.; Wang, W.; She, D. Changes and driving mechanism of water footprint scarcity in crop production: A study of Jiangsu Province, China. Ecol. Indic. 2018, 95, 444–454. [Google Scholar] [CrossRef]
  18. Ababaei, B.; Etedali, H.R. Estimation of water footprint components of Iran’s wheat production: Comparison of global and national scale estimates. Environ. Processes 2014, 1, 193–205. [Google Scholar] [CrossRef] [Green Version]
  19. Rao, J.H.; Hardaha, M.K.; Vora, H.M. The water footprint assessment of agriculture in Banjar river watershed. Curr. World Environ. 2019, 14, 476–488. [Google Scholar] [CrossRef]
  20. Munro, S.A.; Fraser, G.C.; Snowball, J.D.; Pahlow, M. Water footprint assessment of citrus production in South Africa: A case study of the Lower Sundays River Valley. J. Clean. Prod. 2016, 135, 668–678. [Google Scholar] [CrossRef]
  21. Shtull-Trauring, E.; Bernstein, N. Virtual water flows and water-footprint of agricultural crop production, import and export: A case study for Israel. Sci. Total Environ. 2018, 622–623, 1438–1447. [Google Scholar] [CrossRef] [PubMed]
  22. Mekonnen, M.M.; Hoekstra, A.Y. The green, blue and grey water footprint of crops and derived crop products. Hydrol. Earth Syst. Sci. Discuss 2011, 8, 763–809. [Google Scholar] [CrossRef] [Green Version]
  23. Le Roux, B.; Van der Laan, M.; Gush, M.B.; Bristow, K.L. Comparing the usefulness and applicability of different water footprint methodologies for sustainable water management in agriculture. Irrig. Drain. 2018, 67, 790–799. [Google Scholar] [CrossRef]
  24. Barbosa, E.A.; Matsura, E.E.; dos Santos, L.N.; Gonçalves, I.Z.; Nazário, A.A.; Feitosa, D.R. Water footprint of sugarcane irrigated with treated sewage and freshwater under subsurface drip irrigation, in Southeast Brazil. J. Clean. Prod. 2017, 153, 448–456. [Google Scholar] [CrossRef]
  25. Hu, L.; Chun, Y.; Griffith, D.A. Uncovering a positive and negative spatial autocorrelation mixture pattern: A spatial analysis of breast cancer incidences in Broward County, Florida, 2000–2010. J. Geogr. Syst. 2020, 22, 291–308. [Google Scholar] [CrossRef]
  26. Song, J.; Geng, L.; Fahad, S.; Liu, L. Fiscal decentralization and economic growth revisited: An empirical analysis of poverty governance. Environ. Sci. Pollut. Res. 2022, 29, 28020–28030. [Google Scholar] [CrossRef]
  27. Qi, R.; Geng, Y.; Zhu, Q.H. Evaluation of regional water resources utilization based on water footprint method. J. Nat. Resour. 2011, 26, 486–495. [Google Scholar]
Figure 1. Theoretical framework of this paper.
Figure 1. Theoretical framework of this paper.
Processes 10 02105 g001
Figure 2. Study area considered to estimate the agricultural water footprint.
Figure 2. Study area considered to estimate the agricultural water footprint.
Processes 10 02105 g002
Figure 3. Grain yield and agricultural production value in major grain-producing areas from 2000 to 2019.
Figure 3. Grain yield and agricultural production value in major grain-producing areas from 2000 to 2019.
Processes 10 02105 g003
Figure 4. Agricultural water footprint of the main grain-producing provinces in typical year (100 million m3).
Figure 4. Agricultural water footprint of the main grain-producing provinces in typical year (100 million m3).
Processes 10 02105 g004
Figure 5. Agricultural water footprint in the main grain-producing provinces in the typical year.
Figure 5. Agricultural water footprint in the main grain-producing provinces in the typical year.
Processes 10 02105 g005
Figure 6. Time-series variation of the agricultural water footprint in the main grain-producing regions from 2000 to 2019.
Figure 6. Time-series variation of the agricultural water footprint in the main grain-producing regions from 2000 to 2019.
Processes 10 02105 g006
Figure 7. Moran I scatter diagram of agricultural water footprint in the main grain-producing areas from 2005 to 2019.
Figure 7. Moran I scatter diagram of agricultural water footprint in the main grain-producing areas from 2005 to 2019.
Processes 10 02105 g007
Table 1. Selection and statement of influencing factor indicators.
Table 1. Selection and statement of influencing factor indicators.
TypeIndexStatementSymbolUnitData Source
dependent variableagricultural water footprint Y 100 million m3calculated in former part of this paper
Independe-nt variablepolicy factorsproportion of irrigation and water conservancy expenditure in the total amount of water conservancy investmentX1100 million CNYrelated Yearbook
water-saving technologiesproportion of water-saving irrigation areas in the whole arable landX21000 hectaresrelated Yearbook
agricultural productiontotal power of agricultural machineryX310,000 KWrelated Yearbook
social developmenturbanization rateX4%related Yearbook
natural factorsannual precipitation amountX5mmrelated Platform
economic developmentper Capita GDPX6CNY/personChina Statistical Yearbook
industry structureproportion of agricultural output value in total regional output valueX7%China Statistical Yearbook
Table 2. Water footprint per unit production of main grain crops in main grain-producing areas (m3/kg).
Table 2. Water footprint per unit production of main grain crops in main grain-producing areas (m3/kg).
CropsGreen Water FootprintBlue Water FootprintGrey Water FootprintTotal Water Footprint
Grain Crops0.920.410.261.59
Rices0.690.310.231.23
Wheats1.530.670.242.44
Corns0.930.410.241.58
Beans2.611.140.093.84
Table 3. Global Moran’s I index of agricultural water footprint in the main grain-producing areas.
Table 3. Global Moran’s I index of agricultural water footprint in the main grain-producing areas.
YearIE (I)sd (I)zp-Value
20000.563−0.0830.1973.2880.001
20010.576−0.0830.1983.3330.001
20020.622−0.0830.2003.5200.000
20030.576−0.0830.2013.2760.001
20040.517−0.0830.2032.9580.003
20050.599−0.0830.1993.4310.001
20060.582−0.0830.2003.3280.001
20070.569−0.0830.1993.2730.001
20080.520−0.0830.1983.0500.002
20090.480−0.0830.1962.8800.004
20100.742−0.0830.1954.2290.000
20110.609−0.0830.1973.5140.000
20120.423−0.0830.1962.5760.010
20130.578−0.0830.1963.3820.001
20140.589−0.0830.1953.4510.001
20150.575−0.0830.1943.4030.001
20160.518−0.0830.1963.0670.002
20170.489−0.0830.1912.9980.003
20180.513−0.0830.1963.0360.002
20190.386−0.0830.1962.3990.016
Table 4. Spatial econometric model test.
Table 4. Spatial econometric model test.
Test TypeTest ContentStatistical Valuep-ValueConclusion
LM TestLM (SEM) test7.3190.007choose the SDM model
Robust LM (SEM) test12.3430.000
LM (SAR) test21.3260.000
Robust LM (SAR) test26.3490.000
Hausman TestHausman test36.7900.000choose the SDM model and reject random effect
Fixed Effect TestLR (ind or both) test58.6600.000choose the SDM model, reject region fixed effect and time fixed effect, it is better with both fixed time and space
LR (time or both) test220.0800.000
Wald TestWald test SAR26.0800.000choose the SDM model
Wald test SEM28.3300.000
LR TestLR test SAR24.9000.000SDM model cannot be degraded to SAR model
LR test SEM26.7600.000SDM model cannot be degraded to SEM model
Table 5. Regression results of spatial Durbin model on the impact factors of agricultural water footprint in major grain-producing regions.
Table 5. Regression results of spatial Durbin model on the impact factors of agricultural water footprint in major grain-producing regions.
VariableMain Grain-Producing AreasNorthern RegionSouthern Region
Coefficientp-ValueCoefficientp-ValueCoefficientp-Value
lnX1−0.0270.031−0.0380.020−0.0300.054
lnX2−0.0430.053−0.0760.012−0.1010.031
lnX30.1730.000−0.0790.3830.2220.000
lnX4−0.2150.000−0.3490.000−0.2150.000
lnX50.2600.0010.3120.0010.3070.000
lnX6−0.4240.000−0.6540.0010.1690.692
lnX7−0.2110.000−0.1270.226−0.5860.000
W × lnX1−0.0010.9690.0400.287−0.0220.451
W × lnX2−0.0860.063−0.1790.003−0.1140.185
W × lnX30.2100.0360.2170.4050.5150.000
W × lnX4−0.0920.079−0.3300.046−0.0590.104
W × lnX50.2960.0170.1980.263−0.1460.312
W × lnX6−0.2000.361−1.0500.007−0.5070.001
W × lnX7−0.2890.010−0.3380.0800.0080.972
Table 6. Decomposition effect results of influencing factors.
Table 6. Decomposition effect results of influencing factors.
RegionVariableDirect EffectIndirect EffectTotal Effect
Coefficientp-ValueCoefficientp-ValueCoefficientp-Value
Main Grain-
producing Regions
lnX1−0.0270.0380.0010.971−0.0260.400
lnX2−0.0420.075−0.0820.048−0.1230.014
lnX30.1680.0000.1980.0290.3660.000
lnX4−0.2110.000−0.0730.113−0.2840.000
lnX50.2550.0010.2680.0240.5230.000
lnX6−0.4190.000−0.1550.483−0.5740.033
lnX7−0.2020.000−0.2640.019−0.4660.002
Northern Region lnX1−0.0400.0170.0450.1690.0050.889
lnX2−0.0650.033−0.1530.002−0.2180.000
lnX3−0.0990.2790.2290.3320.1300.570
lnX4−0.3290.000−0.2320.118−0.5610.002
lnX50.3030.0020.1240.4670.4270.002
lnX6−0.5950.001−0.8400.013−1.4350.003
lnX7−0.1030.293−0.2830.108−0.3850.138
South Region lnX1−0.0290.067−0.0150.495−0.0430.139
lnX2−0.0960.038−0.0830.211−0.1790.088
lnX30.2050.0000.3960.0000.6010.000
lnX4−0.2120.000−0.02870.205−0.2400.000
lnX50.3120.000−0.1410.2110.1700.070
lnX60.1920.646−0.3990.004−0.2060.612
lnX7−0.5810.0000.0640.738−0.5170.037
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Su, Z.; Zhang, Q. A Study on Spatial-Temporal Differentiation and Influencing Factors of Agricultural Water Footprint in the Main Grain-Producing Areas in China. Processes 2022, 10, 2105. https://doi.org/10.3390/pr10102105

AMA Style

Wang Y, Su Z, Zhang Q. A Study on Spatial-Temporal Differentiation and Influencing Factors of Agricultural Water Footprint in the Main Grain-Producing Areas in China. Processes. 2022; 10(10):2105. https://doi.org/10.3390/pr10102105

Chicago/Turabian Style

Wang, Yun, Zhaoxian Su, and Qingqing Zhang. 2022. "A Study on Spatial-Temporal Differentiation and Influencing Factors of Agricultural Water Footprint in the Main Grain-Producing Areas in China" Processes 10, no. 10: 2105. https://doi.org/10.3390/pr10102105

APA Style

Wang, Y., Su, Z., & Zhang, Q. (2022). A Study on Spatial-Temporal Differentiation and Influencing Factors of Agricultural Water Footprint in the Main Grain-Producing Areas in China. Processes, 10(10), 2105. https://doi.org/10.3390/pr10102105

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop