Next Article in Journal
Influence of Groundwater Discharge on Temporal Evolution in Two Wetlands of an Intensely Anthropized Area: Analysis Using an Integrated Approach
Next Article in Special Issue
Water Footprint of Rangeland Beef Production in New Mexico
Previous Article in Journal
Evaluation and Scenario Prediction of the Water-Energy-Food System Security in the Yangtze River Economic Belt Based on the RF-Haken Model
Previous Article in Special Issue
Past, Present, and Future of Virtual Water and Water Footprint
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Understanding the Spatial-Temporal Changes of Oasis Farmland in the Tarim River Basin from the Perspective of Agricultural Water Footprint

1
College of Water and Architectural Engineering, Shihezi University, Shihezi 832000, China
2
State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, Department of Water Resources, China Institute of Water Resources and Hydropower Research, Beijing 100038, China
3
Ministry of Ecology and Environment Center for Satellite Application on Ecology and Environment, Beijing 100094, China
4
Power China Chengdu Engineering Co., Ltd., Chengdu 610072, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this study and share first authorship.
Water 2021, 13(5), 696; https://doi.org/10.3390/w13050696
Submission received: 26 November 2020 / Revised: 24 February 2021 / Accepted: 1 March 2021 / Published: 4 March 2021
(This article belongs to the Special Issue Water Footprint Assessment Research)

Abstract

:
The Tarim River Basin in China has predominantly assumed the task of commodity cotton and other high water-intensive crop production in recent years. The spatial matching status of agricultural water and land resources is a prerequisite for local economic development. This paper provides an insight into the spatiotemporal variation trends of agricultural production water footprint and oasis farmland in the Tarim River Basin. The degree of spatial mismatching between oasis farmland and crop production water footprints studied in this paper found how the crop water footprint affected the change in oasis farmland area by sensitivity analysis. Time series data covering the period of 1990–2015 were used for the study. The results showed that the annual variation of crop production water footprint and oasis farmland area have experienced upward trends in Tarim River Basin. The blue water makes the largest contribution to the components of the crop production water footprint in each district (all exceeded 77%). The crop production water footprint and oasis farmland area tend to aggregate towards the eastern region. The level of spatial mismatch between the blue water footprint and farmland area fluctuated during the study period, but it was gradually remedied after 2000, while the spatial mismatch between green water footprint and farmland area gradually worsened. The number of districts with mid and high sensitivity to changes in blue water footprint continuously increased during 1990–2005, which revealed that the change in blue water footprint has an increasing influence on oasis farmland. The results can provide operable recommendations for efficient use of water resources, maintaining oasis suitable farmland scale and agricultural sustainable development in the Tarim River Basin.

1. Introduction

The increasingly global scale of environmental degradation in the Anthropocene has engendered broad discussion about the conflicts between water and land resources in inland arid oases [1]. Shortages and inconsistencies between available water and farmland resources are long-term and radically restrictive factors in the process of agriculture modernization [2]. The Tarim River Basin (TRB) is the largest inland basin in China, and within the background of the strategy Western Development Program, the TRB government has undertaken the task of commodity cotton production [3]. This task involves a total cotton yield accounting for more than 48.4% of the production in China. Especially in the past 25 years, because of the impact of high-intensity agricultural exploitation, the newly added farmland area was 12.94 × 103 km2 in TRB, and about 88.7% of the newly added area has been converted from natural grassland [4]. In addition, irrigation water has long-term accounted for more than 90% of the total gross requirements. Soil erosion and pollution caused by the large-scale occupation of agricultural water and land resources have caused the quality of the wetland environment to degrade in the oasis [5]. The volume and reasonable utilization or not of water both directly impact the production efficiency and utilization pattern of oasis farmland, while the degree of reclamation of oasis farmland restricts the development and utilization of water resources [6]. Therefore, making spatial-temporal changes to farmland and agricultural water, and the matching characteristics research between them, is of great practical significance for the optimal allocation of oasis agricultural water and land resources to sustain agricultural production capacity.
Thus far, multiple studies on arid regions have focused on the spatial-temporal distribution characteristics and the matching of water and land resources. Li et al. (2016) evaluated this situation concerning agricultural water and land resources in Yan’an city of China by constructing Lorenz curve and calculating agricultural water and land resource mismatch model [7]. Their research found that the mismatch degree between agricultural water and land resources was higher than that at the national average level for the same period. Obvious spatial differences were observed in the degree of mismatching, which generally improved from south to north. Du et al. (2019) analyzed the matching degree between irrigation water and farmland resources in Ningxia Province of China from the perspective of a Gravity Center model. Results showed that the gravity centers of irrigation water and effective irrigated area both moved in the southeast direction, and their matching degree was superior to that between irrigation water and farmland area [8]. Zhang et al. (2020) studied the matching pattern between renewable water resources and farmland area in Central Asia by constructing a Gini Coefficient model. The authors concluded that the large spatial differences in the matching degree in water distribution and utilization among Central Asian countries, along with the overexploitation of land resources, have ultimately led to the serious water crisis [9]. As shown above, the spatial-temporal distribution and matching situation between water and land resources has been mainly calculated using the following three methods: (i) studying the equilibrium status of regional water and land resources by drawing a Lorenz curve and building Gini Coefficient model, (ii) investigating the spatial–temporal change trend of water and land resources based on the Gravity Center model, and (iii) based on the statistical yearbook’s water consumption and land area data, analyzing the matching index of water and land resources using the water consumption per unit land area as the measurement indicator. Although the aforementioned methods are somewhat successful for measuring the spatial-temporal changes and matching of water and land resources, certain shortcomings have been observed: (i) the Lorenz curve and Gini Coefficient can represent only the overall matching degree between water and land resources, and it reflects the relative relationship between each sub-region in the study area but cannot reveal the absolute matching situation; (ii) research on the spatial-temporal changes of water and land resources include the distance of the center of gravity shift as well as the direction distribution, spatial expansion and aggregation, and thus, the center of gravity model provides unilateral and incomplete results when the center of gravity shift of water and land resources are used only as the explanatory variables; (iii) the matching index can display the relative spatiotemporal ratio between water and land resources, and the land area parameters obtained from the statistical yearbook cannot directly express the specific geographical spatial images. The water resource parameters mainly included water resource volume [10], available water resource volume [11] and irrigation water volume [12]. Currently, a paradigm shift has been observed in the field of global water consumption, with the core element represented by the production water footprint theory [13]. The introduction of production water footprint theory provides new methods and ideas to analyze the real demand and occupation of oases to agricultural water resources from the perspective of socio-hydrology. It can accurately reflect the use of water resources in agricultural production [14]. Veettil and Mishra (2016) considered that, broadly defined, water consumption in agricultural production (crop production water footprint) can be divided into two parts: blue water footprint and green water footprint. The blue water footprint refers to the consumption of surface water and groundwater. The green water footprint refers to the consumption of rainwater insofar as it does not become runoff [15]. Hoekstra and Mekonnen (2012) estimated that the water footprint used for global agricultural production is approximately 7106 Gm3/a, including blue water footprint at 920 Gm3/a and green water footprint at 6186 Gm3/a [16]. The research results of Zhang et al. (2018) indicated that more than 90% of the agricultural production water footprint was dependent on the blue water footprint in Xinjiang Province [17]. Thus, the water resource volume, available water resource volume and irrigation water volume are all unilaterally used as characteristic parameters of the water resource consumption of a region, and the matching situation should be considered from the perspective of the blue–green production water footprint.
This paper investigates the spatial-temporal changes and matching degree of agricultural water footprint and oasis farmland resources. The present study has the following objectives: (a) to reveal the recent spatial-temporal change trends of oasis farmland resources and crop production water footprint, (b) to construct the measuring model of the spatial mismatching index and Gini coefficient between the oasis farmland resources and crop production water footprint, (c) to analyze the impacts of crop production water footprint change on oasis farmland area, and (d) to establish sustainable strategies to increase the efficiency of blue-green water use and ensure suitable oasis farmland scale. TRB is used as a case study, and the analysis is based on a 25-year period from 1990 to 2015. The results can offer novel thoughts for the sustainable management of oasis farmland and for drafting future agricultural development policies.

2. Study Area, Data and Methodology

2.1. Study Site Description

The Tarim River Basin (TRB) in northwest China comprises the catchment areas of 9 stream sources and the main stream of the Tarim River. The basin geographical boundary lies between 71°39′ E and 93°45′ E longitude and 34°20′ N to 43°39′ N latitudes. The TRB is approximately 1.02 million km2, covering five districts, namely Aksu, Bazhou, Kashgar, Hotan and Kezhou (Figure 1). Among all the districts, Kashgar makes the largest contribution (~35.8%) to the water footprint of crop production, and Kezhou makes the smallest contribution (~2.4%). The studied catchment experiences an arid-continental climate with an annual average temperature of 10.6~11.5 °C and 2400~3250 h of sunshine. Mean annual evaporation in the watershed ranges from around 1500 to 2500 mm, far exceeding the mean annual precipitation (<80 mm). TRB’s artificial oasis area accounts for approximately 4.2% of the basin area. About 96.8% of the TRB’s freshwater diversions were for agricultural irrigation. Of the water diverted for irrigation, most (~87.5%) was taken from surface water, with the balance coming from groundwater [18,19].

2.2. Data Sources

2.2.1. Land Use Data Acquisition and Processing

We selected the Landsat scenes acquired for TRB in 1990, 1995, 2000, 2005, 2010 and 2015 epochs to reveal the changing characteristics of artificial oasis farmland area (http://www.gscloud.cn/, accessed on 10 April 2019). The remote sensing images were collected between August and September of the corresponding years to allow a more accurate vegetation interpretation due to cloud cover being at its lowest. The digital elevation model (DEM) data were from the ASTER GDEM V2 data in the geospatial data cloud with a nominal spatial resolution of 30 m. We have carried out large series of calibrations, such as geometric rough correction, radiometric calibration, geometric accurate correction, atmospheric correction and removal of cloud (shadow) pixels, to avoid geometric or radiometric calibration errors. A remote-sensing classification software eCognition 8.7 was performed for the multiresolution segmentation and ground information extraction of the remote-sensing image. These temporally aggregated images were then combined with the high-resolution satellite data in Google Earth. Visual interpretation was performed with reference to stereoscopic aerial photographs and field surveys. It is estimated to have a high classification accuracy of 90–93%.

2.2.2. Meteorological, Water Resources and Statistical Data

Meteorological data (sunshine hours, wind speed, humidity, effective precipitation and air temperature) were obtained from 37 meteorological stations supplied by China Meteorological Data Sharing Service Network (http://data.cma.cn/, accessed on 10 May 2019). The Xinjiang Division of Water Resources is the provider of most of the data used for this study, such as municipal and industrial water resources development and utilization, agricultural and ecological water demand and water consumption. Additionally, the other statistical data were sourced directly from the Xinjiang Statistical Service, such as demographic information, crop yield, planting area and effective irrigation area.

2.3. Methodology

2.3.1. Water Footprint of Crop Production

We calculated the production water footprint of 15 crops: rice, wheat, coarse cereals, soybeans, cotton, oil plants, sugar beets, vegetables, melons, potatoes, alfalfa, grapes, apples, fragrant pears and red jujube. These 15 crop types were estimated to occupy 97.4% of the total cultivated area in the TRB. Therefore, we took the sum of all these 15 crops as the TRB’s crop production water footprint. First, we calculated the ETc (evapotranspiration) and Peff (effective rainfall) of these crops in each district during the growth period, respectively. The theoretical method commonly used to calculate ETc (mm) uses the Penman–Monteith method, which can accurately calculate ETc under different regions and climatic conditions [20]. Peff (mm) uses the CROPWAT model developed by the Food and Agriculture Organization of the United Nations (FAO). Second, we calculated the WFpro (crop production water footprint), which can be divided into WFblue (blue water footprint) and WFgreen (green water footprint) from the water consumption perspective. The WFpro (m3), WFblue (m3) and WFgreen (m3) are expressed as follows [21]:
W F p r o = n = 1 n W F i ( b l u e , g r e e n )
W F b l u e = ( 10 × d = 1 lg p max ( 0 , E T c P e f f ) / Y ) × P i
W F g r e e n = ( 10 × d = 1 lg p min ( E T c , P e f f ) / Y ) × P i
In Equation (1) to Equation (3), Pi is the total output of crop i (t). Y refers to the crop yield per unit area (t/ha). The factor 10 is the conversion coefficient of water depth into water per unit area of land area. ∑ is the cumulative amount of blue water (or green water) from planting period to harvest period. lgp is the length of the growth period (d).

2.3.2. Standard Deviation Ellipse Analysis

The Standard Deviation Ellipse (SDE), or directional distribution analysis, has unique advantages in exploring the spatial-temporal changes in resource factors. It can visually and accurately reveal these resource factors’ evolution characteristics in two-dimensional space [22,23,24,25]. Therefore, by introducing the concept of the SDE, the spatial-temporal changes model of oasis agricultural water and land resources in TRB was constructed. The model judges the main trend direction and dispersion degree of the distribution of water and land resources in the oasis through the elliptical coverage area. We also calculated the moving direction and deviation distance of the center of gravity of WFpro and oasis farmland area (OFa) and analyzed the change trajectories of their centers of gravity. If the study area consists of n units and (xi, yi)t is the geometric coordinate of i-th unit (i = 1, 2, 3,…, n) at time t, the WFpro’s gravity center coordinate G ( x ¯ , y ¯ ) at time t can be written as:
G ( x ¯ , y ¯ ) = ( i = 1 n x i × m i t / i = 1 n m i t , i = 1 n y i × m i t / i = 1 n m i t )
where m i t is the attribute value of the i-th unit at time t.
The moving distance of the gravity center (D, km) can be calculated by:
D = C × ( x ¯ t 2 x ¯ t 1 ) 2 + ( y ¯ t 2 y ¯ t 1 ) 2
where ( x ¯ t 1 , y ¯ t 1 ) and ( x ¯ t 2 , y ¯ t 2 ) are the gravity center coordinate at the time t1 and t2 of the study and C is a constant with the value of 111.111 (km), denoting the conversion coefficient from longitude and latitude coordinates on earth to plane distance.
The deviation angle (θ) in the SDE model is the angle of clockwise rotation. The deviation angle (θ), principal axis and auxiliary axis of the SDE can be calculated by:
tan θ = [ ( i = 1 n x i 2 i = 1 n y i 2 ) + ( i = 1 n x i 2 i = 1 n y i 2 ) 2 + 4 ( i = 1 n x i y i ) 2 ] / ( 2 i = 1 n x i y i )
σ x = i = 1 n ( m i t x i cos θ m i t y i sin θ ) 2 / i = 1 n ( m i t ) 2
σ y = i = 1 n ( m i t x i sin θ m i t y i cos θ ) 2 / i = 1 n ( m i t ) 2
where x i and y i are the deviations of gravity center and σ x and σ y are the standard deviations of principal axis and auxiliary axis, respectively.

2.3.3. Theory of Water Footprint and Farmland Matching

The Gini Coefficient (GC) proposed by economist Gini Corrado was first applied to the study of income disequilibrium based on the Lorentz curve [26]. Because the spatial distribution of natural resources shows diversity, the GC is also applicable in the study of the patterns of oasis water and land resource matching. To investigate the balanced state between WFpro and oasis farmland area (OFa), the GC can be calculated as:
G C = 1 [ l × ( 2 a ω ) ]
where l and ω refer to, respectively, the proportions of OFa and WFpro in each district, and a reference to the cumulative percentage of the OFa of the TRB. The GC ranges from 0 to 1. A greater degree of equilibrium between WFpro and OFa will lead to a smaller GC and vice versa.
Following the study of Gobillon et al. (2007) [27], spatial mismatch is used to describe the imbalance between water and land resources in space, which has a profound impact on regional water management and land-use policies. The spatial mismatch index was introduced to measure the spatial relationship between WFpro and OFa. The formula is as follows:
S M I ( W F p r o _ O F a ) i = ( W F i i = 1 n W F i O F i i = 1 n O F i ) × 100
S M I ( W F p r o _ O F a ) = i = 1 n | S M I ( W F p r o _ O F a ) i |
where SMI(WFpro_OFa)i refers to the spatial mismatch index of region i between WFpro and OFa, WFi and OFi refer to, respectively, the region i’s crop production water footprint (km3) and oasis farmland area (km2), and S M I ( W F p r o _ O F a ) refers to the total level of spatial mismatch in TRB. A higher value of SMI(WFpro_OFa)i means there is more WFpro in each of unit oasis farmland of region i, and a reduced SMI(WFpro_OFa)i means the oasis farmland resource is utilized extensively and the WFpro is weakened in region i. In order to judge whether the spatial relationship between WFpro and OFa is matching in each region, we chose the Jenks Natural Breaks Classification to help us set standard values of the results [28]. If | S M I ( W F p r o _ O F a ) i | is less than the standard value, it shows that the relationship between WFpro and OFa is a spatial match in region i and vice versa.

2.3.4. Sensitivity Index (SI)

A sensitivity assessment index (SI) was introduced to quantitatively analyze oasis farmland’s sensitivity to WFpro. It aims to explain the variations in OFa caused by the change in WFpro for investigating the potential impacts between them. The formula is as follows:
S I i = | ( O F t 2 O F t 1 ) / O F t 1 ( W F t 2 W F t 1 ) / W F t 1 |
where SIi represents the value of sensitivity index in region i, OFt1 and OFt2 are the OFa at the start and end of study period (km2), and WFt1 and WFt2 refer to, respectively, WFpro at the beginning and end of study period (m3). A higher value of SIi reflects a higher sensitivity of OFa to WFpro change, which means that a great change in OFa could be altered by the subtle change in WFpro [29].

3. Results

3.1. Crop Production Water Footprint and Farmland Resources

Figure 2 depicts the changing trends of the crop production water footprint (WFpro) and oasis farmland area (OFa) in the TRB. Results showed that the total WFpro increased by 184.2% (WFblue and WFgreen increased by 185.4% and 175.0%, respectively) in the past 25 years. The annual average value of WFpro showed obvious differences between districts; in descending order of this value, the disctricts were Kashgar, Aksu, Bazhou, Hotan and Kezhou. Among them, the WFpro of Kashgar was 14.96 times larger than that of Kezhou. Similarly to WFpro, the WFblue of each district showed a progressively increasing trend from 1990 to 2015. Kashgar made the largest contribution (36.1%) to the components of the WFblue among all the districts in TRB. Aksu and Bazhou, ranked in second and third place, were 32.9% and 15.9%. Compared to the WFblue, the inter-annual variability of the WFgreen showed a relatively slightly increasing trend in each district. The WFpro in each region was dominated by blue water (~87.3%), and the proportion of green water was relatively small (~12.7%).
OFa increased by 51.9% from 1990 to 2015. The inter-annual variability of each district had similar changing trends of moderately increasing during the entire study period. The annual average OFa ranking of the districts was in descending order in Kashgar, Aksu, Bazhou, Hotan and Kezhou. The sum of OFa in Kashgar and Aksu accounts for more than 70% of the OFa in TRB. Overall, it indicated that Kashgar and Aksu constitute the two main districts of the WFpro and OFa development in the TRB.

3.2. Spatial Distribution Characteristics of Crop Production Water Footprint and Oasis Farmland

As shown in Figure 3, the spatial distribution of WFpro and OFa in TRB is not uniform and presents the characteristics of large changes of WFpro and OFa in Aksu and Kashgar, and a smaller change of WFpro and OFa in other districts. The change in WFpro and OFa of some districts show a degree of mismatch. For example, the changes in WFpro in Aksu and Kashgar were roughly consistent (both greater than 5 km3), but the change in OFa in Kashgar was only 63.9% of that of Aksu.
In the past 25 years, the movement pattern of WFpro’s center of gravity was roughly inconsistent with that of OFa. The gravity center of WFpro moved more obviously in longitude. In terms of moving direction, the gravity centers of WFpro and OFa both moved to the southeast between 2000 and 2005 but were different in other periods. In terms of moving distance, the location characteristics of OFa gradually expanded 36.01 km to the southeast, while the center of gravity of WFpro moved 18.02 km to the northeast during the entire period. The centers of gravity of WFpro and OFa both moved to the east, showing that WFpro and OFa in eastern regions accounted for the increasing proportion in the whole study area.

3.3. Spatial Matching Analysis of Crop Production Water Footprint and Oasis Farmland

The results of spatial mismatch index (SMI) and Gini coefficient (GC) in TRB are shown in Figure 4. From 1990 to 2005, the SMI (WFgreen_OFa) presents an increasing trend, while the SMI (WFblue_OFa) presented a roughly decreasing trend. In the same period, the distance between the centers of gravity of WFpro and OFa was closer, as seen in Figure 3. These showed that the spatial mismatch between WFpro and OFa was remedied because agricultural water was dominated by WFblue in TRB. From 2005 to 2010, both the SMI and GC of (WFblue_OFa) presented a decreasing trend, and the spatial match between WFblue and OFa was improved again. The GC (WFblue_OFa) was on a decline after 2010, but the SMI sharply increased from 8.34 to 13.09, which reflected that the degree of equilibrium between WFblue and OFa decreased. On the whole, the variation tendency of the SMI and GC reveal that the spatial mismatch between WFblue and OFa in TRB was gradually remedied after 2000. However, the spatial mismatch between WFgreen and OFa was gradually exacerbated.
Using ArcGIS 10.2, we set −2.0 and 2.0 as the standard values to visualize the matching patterns between WFpro and OFa at each district based on Natural Breaks Classification (Jenks), as shown in Figure 5. From 1990 to 2015, the OFa with serious spatial mismatch and inefficient use of blue water was mainly distributed in Hotan, and the OFa with serious spatial mismatch and inefficient use of green water was mainly distributed in Aksu. From 2010 to 2015, the OFa with inefficient blue water use extended to Kashgar. The OFa with highly efficient blue water utilization was mainly distributed in Aksu and Bazhou. From 2000 to 2010, The SMI (WFblue_OFa) of Aksu and Bazhou started to decrease, coupled with the increase in SMI (WFgreen_OFa). As can be seen in Figure 3, the distance between the center of gravity of WFpro and that of OFa was constantly getting closer, indicating that agricultural water use efficiency of Aksu and Bazhou improved. According to the results shown in Figure 3 and Figure 4, we also found that the improvement of Hotan with agricultural water use efficiency was the reason why the SMI (WFpro_OFa) and GC (WFpro_OFa) both decreased and the spatial mismatch in Hotan decreased from 2000 to 2010. Until 2015, the poor spatial mismatch in Hotan was remedied due to the relative stability of increasing trends in the WFpro and agricultural water use efficiency.

3.4. Sensitivity Analysis of Crop Production Water Footprint and Oasis Farmland

We calculated the sensitivity index (SI) of OFa to WFpro of each district in TRB, as shown in Figure 6. The sensitivity was divided into four levels: 0< |SI|≤ 0.5 indicates that OFa is non-sensitivity to changes in WFpro, 0.5 < |SI| ≤ 1.0 indicates low sensitivity, 1.0 < |SI| ≤ 1.5 indicates moderate sensitivity, and |SI| > 1.5 indicates high sensitivity.
The number of districts with moderate and high sensitivity of WFblue_OFa first increased and then decreased, while the quantity of districts with non- and low sensitivity first decreased and then increased during the study period. During 2000–2005, there were four districts with moderate and high sensitivity. In comparison, all districts in other periods had non- and low sensitivity. Especially in 1990–1995 and 2005–2010, the OFa in all districts was non-sensitive to WFblue change in TRB, while only two districts in TRB showed non-sensitivity during 2010–2015 (Kashgar and Aksu). In addition, during the study period, there were no districts with the high sensitivity of OFa to green water footprint change in TRB. During 1995–2000 and 2000–2005, there was one moderate-sensitive area in each. During 2010–2015, the number of mid-sensitive areas rose to two. The OFa in all districts were mainly non-sensitive or low-sensitive to green water footprint change during other periods. From the perspective of spatial distribution, the districts with moderate and high sensitivity of OFa to the WFblue in TRB were gradually expanded westward. Especially during 2000–2005, the OFa in Aksu and Kashgar, located in the northwest of the TRB, was significantly more sensitive to WFblue change than the other three districts. However, the sensitivity of OFa to WFgreen change was relatively insignificant in spatial distribution.

4. Discussion

4.1. Main Factors Influencing Crop Water Footprint and Farmland Resource Changes

Both the WFblue and OFa continue to increase and aggregate in spatial-temporal changes based on the historical trends. Without the implementation of the corresponding regulatory policies, the growth of WFblue driven by the expansion of OFa will lead to a reduction in ecological land, which will result in obvious declines in habitat quality, water economic production efficiency and carbon sequestration [30]. As the marginal efficiency of external inputs factors (urbanization and agricultural mechanization) has declined, the agricultural production scale has continued to expand in TRB. In future sustainable development scenarios, it is recommended that the government proposes a plan to avoid long-term expansion of oasis farmland, which is the most effective in maintaining the suitable oasis farmland scale. Under this policy scenario, both WFblue and WFgreen will decline. In addition, future sustainable development policies should also effectively maintain and improve carbon sequestration and habitat quality by reducing the occupation of ecological land by farmland. Existing studies have shown that the implementation of such management measures can enhance the overall ecosystem service value [31,32]. The development of water-saving technology has proven to be highly effective in improving water efficiency for policymakers. The water-saving irrigated area in the TRB has increased by 7033 km2 and is irrigated by drip irrigation systems covering the period between 2005 and 2015. With the widespread adoption of water-saving engineering technologies such as prevention of channel leakage, water delivery with low-pressure pipe and drip irrigation systems, the WFpro per unit yield of some crops in the TRB has declined [33]. However, water-saving irrigation scheduling, optimized water allocation in irrigation areas and irrigation water demand forecasting technology have not been improved. These can be the key drivers of increasing local WFblue. The global green plus blue global WF benchmark methodology has been widely used when assessing agricultural water efficiency [34]. Nevertheless, the TRB has a typical temperate continental arid climate, characterized by low rainfall, high evapotranspiration and extremely scarce green water resources. Therefore, the TRB dominated by irrigated agriculture has lower green water production efficiency as compared to the humid regions dominated by rain-fed agriculture [35]. More attention should be paid to the relationship between agricultural production and ecological restoration during the implementation of future development policies, developing better environmental flows assessment tools to accurately estimate the extent of water scarcity per catchment and to improve integrated water and land resources management level [36,37].
In terms of improving water efficiency, we suggest that more attention should be paid to the relationship between agricultural production and ecological restoration, and to improving integrated water and land resources management levels during the implementation of future development policies. This can ensure reasonable land utilization and development under the constraints of the regional water resource carrying capacity and can reduce the pressure caused by excessive WFpro while satisfying agricultural development. OFa was shown to be significantly affected by the climatic characteristics and water quality conditions in Hotan and Aksu. Moreover, the change in land use transition and crop plant structure can also change the OFa in the two districts by affecting the WFpro [38]. Due to the accelerated process of urbanization, a large amount of oasis farmland was inefficiently utilized, abandoned or transferred into other land types [39]. Farmers in Aksu and Bazhou interested only in profit potential would therefore be more inclined to plant economic crops (e.g., cotton and red jujube) rather than food crops to keep their income and living standard more stable. The large-scale production of these high water-intensive economic crops have further increased the local WFpro per unit yield.

4.2. Implications and Suggestions from Spatial Matching Theory and Sensitivity Analysis

Although the relationship between irrigation water and effective irrigated area distribution patterns has been noted by a number of research studies, the correlation between the actual water use of crops and farmland landscapes has not yet been fully appreciated [40]. It was believed that the amount of water resource consumed by humans in the entire economic production process (including physical water and virtual water) and the scale of farmland were more closely related in time and space. In the present study, the WFpro and OFa of TRB were moving from the northern regions with a scarcity of water resource conditions to eastern regions, while the economic water productivity per unit area of crops in the eastern regions did not increase significantly. In the long run, this inconsistent distribution pattern will pose risks in terms of local agricultural production. In TRB, the phenomenon of spatial mismatch between WFpro and OFa is more obvious, which can cause problems such as an imbalance in regional water and land structures and ecological risks. Meanwhile, the increase in the number of moderate- and high-sensitivity districts also revealed that the WFpro has a greater impact on OFa changes. However, with the scarcity of high-quality farmland and the rapid development of the primary industry, it has become increasingly difficult to implement Returning Farmland to Forest Program (RFFP) by controlling the WFblue [41].
The Chinese government has implemented a series of policy transmission mechanisms to rationally allocate water and farmland resources. These decision-making sessions have proven effective and feasible in some pilot regions with water shortages [42]. Measures such as changing the farmland operation scale and promoting farmland transfer are based on the successful experience of the above-mentioned pilot regions. Briefly, the TRB’s policymakers can try to take these measures to alleviate the structural imbalance between WFblue and OFa. Policymakers can use subsidy policies to realize such a transformation of TRB from small-scale to moderate-scale agricultural production. Promoting the farmland transfer may enhance the contribution of WFpro per unit to the improvement of crop productivity and benefit to some extent, including subsidizing moderate-scale farmland management, improving the registration system of farmland contractual management rights, standardizing control of irrigation water use and encouraging innovation in forms of farmland transfer. The cross-regional water transfer and farmland scale control integrated planning should be carried out to alleviate the current disparities in agricultural development and water use. In oasis areas with vulnerable ecologies, crop rotation and fallow systems are gaining greater attention for sustainable intensification of agro-ecosystem. It is essential to encourage farmers to increase the proportion of low-water-intensive crops planted and render or increase subsidies for farmers. The role of markets should be given full play in the optimization of farmland resource allocation to conduct a pattern of Transferrable Development Rights transaction (TDR) [43]. Moreover, modern oasis land-use management practices may achieve a win–win effect between the improvement in blue water utilization efficiency and quality assurance of oasis farmland. WFgreen is a kind of water consumption restricted by the natural conditions of the region; unlike WFblue, its impact on the change in OFa was limited. At present, the effective management of green water can only be achieved via virtual water trade [44]. The government should focus on China’s “One Belt One Road” initiative so as to optimize and develop the product structure of processing and service industries. Compared with the agriculture sector, the direct blue–green physical water consumption and virtual water export volume per unit output in secondary and tertiary industries were relatively low [45].

5. Conclusions

Using TRB as a case study, this paper assessed the spatial-temporal variability of WFpro and OFa over the period 1990–2015. Furthermore, the potential impact of WFpro changes in OFa was examined, which facilitates an improved understanding regarding the blue–green water use efficiency.
On the whole, the OFa was more strongly associated with WFblue change than WFgreen. In 1990–2015, the increasing trend of WFblue was more remarkable in various districts compared to that of WFgreen. The average annual WFblue apparently grew more quickly than WFgreen in TRB. The remarkable, significant increase in the TRB’s WFblue in the study period mainly occurred in Kashgar and Aksu districts, where the oasis farmland area displays a continuous trend of expansion. The analysis changes in WFpro and OFa in the TRB demonstrated that the SDE model gives insights into spatial distribution characteristics analysis. Our findings show that the spatial aggregation of WFpro and OFa in the eastern districts of TRB continuously increased.
Assessing the spatial mismatch between WFblue and OFa in TRB, we found that the overall level of spatial mismatch fluctuated during the study period. After 2000, the spatial mismatch between WFblue and OFa in TRB was gradually remedied, but the spatial mismatch between WFgreen and OFa gradually deteriorated. The districts with the most inefficient use of blue water were mainly distributed in the oasis farmland of Hotan. With the popularization of water-saving technologies in oasis farmland areas of Aksu and Bazhou, these two districts became more dominant over time in the regional pattern of crop production and agricultural security in the TRB.
The spatial-temporal pattern of sensitivity showed that the number of districts with moderate and high sensitivity to changes in the WFblue increased continuously and the spatial distribution of districts shifted from the east to the west during 1990–2005. The change in WFblue has a great increasing influence on OFa in TRB with time, especially in the western districts. The continuous acceleration of urbanization and the excessive reclamation of farmland after 2005 led to the inefficient use of a large amount of oasis farmland, which lowered the SI (WFblue_OFa) in the TRB. For the sustainability of the oasis agricultural economy, TRB’s government should improve the local agricultural blue–green water use efficiency by changing the oasis farmland operation pattern and promoting farmland transfer. These ways can not only increase the industrial diversity of agriculture but also have great importance for controlling oasis farmland scale.

Author Contributions

J.Y. and X.D. wrote the manuscript and prepared figures. A.L., X.H., H.G. and J.D. advised the study design and data analyses. J.Z. and C.R. collected the data. A.L. and X.D. revised the paper. All authors discussed the results and commented on the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The authors recognize the research support of the Ministry of Science and Technology of the People’s Republic of China through grants 2016FYA0601600 and 2017YFC0404300. This work was also supported by the National Natural Science Foundation of China (U1803244).

Acknowledgments

Conversations with Xu Zhao were helpful in this research. Comments from Pieter van Oel were critical in bringing important focus to this study. The authors thank Yawen Shao from the University of Melbourne for polishing the language of the manuscript. We are grateful for the time and valuable comments of the anonymous reviewers in improving this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xue, X.; Liao, J.; Hsing, Y.; Huang, C.; Liu, F. Policies, land use, and water resource management in an arid oasis ecosystem. Environ. Manag. 2015, 55, 1036–1051. [Google Scholar] [CrossRef]
  2. Fitton, N.; Alexander, P.; Arnell, N.; Bajzelj, B.; Calvin, K.; Doelman, J.; Gerber, J.S.; Havlik, P.; Hasegawa, T.; Herrero, M.; et al. The vulnerabilities of agricultural land and food production to future water scarcity. Glob. Environ. Chang. 2019, 58, 101944. [Google Scholar] [CrossRef]
  3. Zhao, X.; Othmanli, H.; Schiller, T.; Zhao, C.; Sheng, Y.; Zia, S.; Muller, J.; Stahr, K. Water use efficiency in saline soils under cotton cultivation in the Tarim River Basin. Water 2015, 7, 3103–3122. [Google Scholar] [CrossRef]
  4. Huang, S.; Wortmann, M.; Duethmann, D.; Menz, C.; Shi, F.; Zhao, C.; Su, B.; Krysanova, V. Adaptation strategies of agriculture and water management to climate change in the Upper Tarim River basin, NW China. Agric. Water Manag. 2018, 203, 207–224. [Google Scholar] [CrossRef]
  5. Xue, L.; Wang, J.; Zhang, L.; Wei, G.; Zhu, B. Spatiotemporal analysis of ecological vulnerability and management in the Tarim River Basin, China. Sci. Total Environ. 2019, 649, 876–888. [Google Scholar] [CrossRef]
  6. Waldron, B.; Gui, D.; Liu, Y.; Feng, L.; Dai, H. Assessing water distribution and agricultural expansion in the Cele Oasis, China. Environ. Monit. Assess. 2020, 192, 1–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Li, H.; Zhou, W.; Zhuang, Y.; Cong, M.; Liu, B. Agricultural water and soil resources matching patterns and carrying capacity in Yan’an City. Trans. Chin. Soc. Agric. Eng. 2016, 32, 156–162. (In Chinese) [Google Scholar]
  8. Du, J.; Yang, Z.; Wang, H.; Yang, G.; Li, S. Spatial–temporal matching characteristics between agricultural water and land resources in Ningxia, Northwest China. Water 2019, 11, 1460. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, Y.; Yan, Z.; Song, J.; Wei, A.; Sun, H.; Cheng, D. Analysis for spatial-temporal matching pattern between water and land resources in Central Asia. Hydrol. Res. 2020, 10, 1–20. [Google Scholar]
  10. Yang, Q.; Yang, X.L.; Qiao, C.; Zuo, Y.; Cui, Y.; Jin, J. Matching analysis of water and land resources in Anhui Province based on Gini coefficient. Chin. Agric. Sci. Bull. 2016, 32, 72–76. (In Chinese) [Google Scholar]
  11. Jiang, Q.X.; Fu, Q.; Wang, Z.L.; Jiang, N. Spatial matching patterns of land and water resources in Sanjiang Plain. J. Nat. Resour. 2011, 26, 270–277. (In Chinese) [Google Scholar]
  12. Zheng, J.; Zhao, X.; Cao, X.; Sun, S.; Zhang, L. Study on spatiotemporal matching pattern of agricultural water and land resources in Hetao irrigation district. Res. Soil Water Conserv. 2015, 22, 132–136. [Google Scholar]
  13. Sun, S.; Wu, P.; Wang, Y.; Zhao, X.; Liu, J.; Zhang, X. The impacts of interannual climate variability and agricultural inputs on water footprint of crop production in an irrigation district of China. Sci. Total Environ. 2013, 444, 498–507. [Google Scholar] [CrossRef]
  14. Chapagain, A.K.; Hoekstra, A.Y. The blue, green and grey water footprint of rice from production and consumption perspectives. Ecol. Econ. 2011, 70, 749–758. [Google Scholar] [CrossRef]
  15. Veettil, A.V.; Mishra, A.K. Water security assessment using blue and green water footprint concepts. J. Hydrol. 2016, 542, 589–602. [Google Scholar] [CrossRef]
  16. Hoekstra, A.Y.; Mekonnen, M.M. The water footprint of humanity. Proc. Natl. Acad. Sci. USA 2012, 109, 3232–3237. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, P.; Deng, X.; Long, A.; Hai, Y.; Li, Y.; Wang, H.; Xu, H. Impact of social factors in agricultural production on the crop water footprint in Xinjiang, China. Water 2018, 10, 1145. [Google Scholar] [CrossRef] [Green Version]
  18. Xinjiang Uygur Autonomous Region Statistics Bureau. Xinjiang Statistical Yearbook; China Statistics Press: Beijing, China, 2019. (In Chinese) [Google Scholar]
  19. Ling, H.; Guo, B.; Zhang, G.; Xu, H.; Deng, X. Evaluation of the ecological protective effect of the “large basin” comprehensive management system in the Tarim River basin, China. Sci. Total Environ. 2019, 650, 1696–1706. [Google Scholar] [CrossRef]
  20. Yu, J.; Long, A.; Deng, X.; He, X.; Zhang, P.; Wang, J.; Hai, Y. Incorporating the red jujube water footprint and economic water productivity into sustainable integrated management policy. J. Environ. Manag. 2020, 269, 110828. [Google Scholar] [CrossRef]
  21. Hoekstra, A.Y.; Chapagain, A.K.; Aldaya, M.M.; Mekonnen, M.M. The Water Footprint Assessment Manual: Setting the Global Standard; Earthscan: London, UK, 2011. [Google Scholar]
  22. Liu, H.; Fang, C.; Miao, Y.; Ma, H.; Zhang, Q.; Zhou, Q. Spatio-temporal evolution of population and urbanization in the countries along the Belt and Road 1950–2050. J. Geogr. Sci. 2018, 28, 919–936. [Google Scholar] [CrossRef] [Green Version]
  23. Yang, Z.; Xiu, Q.; Chen, D. Historical changes in the port and shipping industry in Hong Kong and the underlying policies. Transp. Policy. 2019, 82, 138–147. [Google Scholar] [CrossRef]
  24. Ma, J.; Chen, H.; Gao, X.; Xiao, J.; Wang, H. African swine fever emerging in China: Distribution characteristics and high-risk areas. Prev. Vet. Med. 2020, 175, 104861. [Google Scholar] [CrossRef] [PubMed]
  25. Ayana, E.K.; Ceccato, P.; Fisher, J.R.B.; DeFries, R. Examining the relationship between environmental factors and conflict in pastoralist areas of East Africa. Sci. Total Environ. 2016, 557, 601–611. [Google Scholar] [CrossRef]
  26. Groves-Kirkby, C.J.; Denman, A.R.; Phillips, P.S. Lorenz curve and Gini coefficient: Novel tools for analysing seasonal variation of environmental radon gas. J. Environ. Manag. 2009, 90, 2480–2487. [Google Scholar] [CrossRef]
  27. Gobillon, L.; Selod, H.; Zenou, Y. The mechanisms of spatial mismatch. Urban Stud. 2007, 44, 2401–2427. [Google Scholar] [CrossRef] [Green Version]
  28. Naghibi, S.A.; Moghaddam, D.D.; Kalantar, B.; Pradhan, B.; Kisi, O. A comparative assessment of GIS-based data mining models and a novel ensemble model in groundwater well potential mapping. J. Hydrol. 2017, 548, 471–483. [Google Scholar] [CrossRef]
  29. Han, Z.; Song, W.; Deng, X. Responses of ecosystem service to land use change in Qinghai Province. Energies 2016, 9, 303. [Google Scholar] [CrossRef] [Green Version]
  30. Gong, B.; Liu, Z. Assessing impacts of land use policies on environmental sustainability of oasis landscapes with scenario analysis: The case of northern China. Landsc. Ecol. 2020, 1–20. [Google Scholar] [CrossRef]
  31. Zhang, Z.; Xia, F.; Yang, D.; Huo, J.; Wang, G.; Chen, H. Spatiotemporal characteristics in ecosystem service value and its interaction with human activities in Xinjiang, China. Ecol. Indic. 2020, 110, 105826. [Google Scholar] [CrossRef]
  32. Wang, Y.; Zhang, S.; Zhen, H.; Chang, X.; Shataer, R.; Li, Z. Spatiotemporal evolution characteristics in ecosystem service values based on land use/cover change in the Tarim River Basin, China. Sustainability 2020, 12, 7759. [Google Scholar] [CrossRef]
  33. Long, A.; Zhang, P.; Hai, Y.; Deng, X.; Li, J.; Wang, J. Spatio-temporal variations of crop water footprint and its influencing factors in Xinjiang, China during 1988–2017. Sustainability 2020, 12, 9678. [Google Scholar] [CrossRef]
  34. Mekonnen, M.M.; Hoekstra, A.Y. Water footprint benchmarks for crop production: A first global assessment. Ecol. Indic. 2014, 46, 214–223. [Google Scholar] [CrossRef] [Green Version]
  35. Vanham, D.; Mekonnen, M.M. The scarcity-weighted water footprint provides unreliable water sustainability scoring. Sci. Total Environ. 2021, 756, 143992. [Google Scholar] [CrossRef]
  36. Vanham, D.; Hoekstra, A.Y.; Wada, Y.; Bouraoui, F.; de Roo, A.; Mekonnen, M.M.; van de Bund, W.J.; Batelaan, O.; Pavelic, P.; Bastiaanssen, W.G.M.; et al. Physical water scarcity metrics for monitoring progress towards SDG target 6.4: An evaluation of indicator 6.4. 2 “Level of water stress”. Sci. Total Environ. 2018, 613, 218–232. [Google Scholar] [CrossRef] [PubMed]
  37. Mekonnen, M.M.; Hoekstra, A.Y. Four billion people facing severe water scarcity. Sci. Adv. 2016, 2, e1500323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Fang, G.; Chen, Y.; Li, Z. Variation in agricultural water demand and its attributions in the arid Tarim River Basin. J. Agric. Sci. 2018, 156, 301–311. [Google Scholar] [CrossRef]
  39. Yang, Y.; Liu, Y. Spatio-temporal analysis of urbanization and land and water resources efficiency of oasis cities in Tarim River Basin. J. Geogr. Sci. 2014, 24, 509–525. [Google Scholar] [CrossRef]
  40. Chai, J.; Wang, Z.; Yang, J.; Zhang, L. Analysis for spatial-temporal changes of grain production and farmland resource: Evidence from Hubei Province, central China. J. Clean. Prod. 2019, 207, 474–482. [Google Scholar] [CrossRef]
  41. Liu, Y.; Chen, Y. Impact of population growth and land-use change on water resources and ecosystems of the arid Tarim River Basin in Western China. Int. J. Sustain. Dev. World Ecol. 2006, 13, 295–305. [Google Scholar] [CrossRef]
  42. Chinese Ministry of Agriculture; Chinese Ministry of Finance. Guiding Opinions on Adjusting the Three Agricultural Subsidies Policies. Available online: http://nys.mof.gov.cn/zhengfuxinxi/czpjZhengCeFaBu_2_2/201505/t20150522_1237664.html (accessed on 28 December 2020).
  43. Yildiz, O. Object-based modeling of restrictions on the sale of agricultural land. Land Use Policy 2019, 82, 538–549. [Google Scholar] [CrossRef]
  44. Zhang, X.; Liu, J.; Zhao, X.; Yang, H.; Deng, X.; Jiang, X.; Li, Y. Linking physical water consumption with virtual water consumption: Methodology, application and implications. J. Clean. Prod. 2019, 228, 1206–1217. [Google Scholar] [CrossRef]
  45. Gao, D.; Long, A.; Yu, J.; Xu, H.; Su, S.; Zhao, X. Assessment of inter-sectoral virtual water reallocation and linkages in the Northern Tianshan Mountains, China. Water 2020, 12, 2363. [Google Scholar] [CrossRef]
Figure 1. Location, crop production water footprint and oases of the study basin.
Figure 1. Location, crop production water footprint and oases of the study basin.
Water 13 00696 g001
Figure 2. Recent trends in crop production water footprint and farmland area in each district of the Tarim River Basin (TRB).
Figure 2. Recent trends in crop production water footprint and farmland area in each district of the Tarim River Basin (TRB).
Water 13 00696 g002
Figure 3. Change and migration trajectories of crop production water footprint (WFpro) and oasis farmland area (OFa) in TRB during 1990 to 2015.
Figure 3. Change and migration trajectories of crop production water footprint (WFpro) and oasis farmland area (OFa) in TRB during 1990 to 2015.
Water 13 00696 g003
Figure 4. The Gini Coefficient and spatial mismatch index of crop water footprint and oasis farmland.
Figure 4. The Gini Coefficient and spatial mismatch index of crop water footprint and oasis farmland.
Water 13 00696 g004
Figure 5. The spatial-temporal pattern of spatial mismatch index from 1990 to 2015. (a) The spatial mismatch index between blue water footprint and oasis farmland area (SMI (WFblue_OFa)); (b) The spatial mismatch index between green water footprint and oasis farmland area (SMI (WFgreen_OFa)); (c) The spatial mismatch index between crop production water footprint and oasis farmland area (SMI (WFpro_OFa)).
Figure 5. The spatial-temporal pattern of spatial mismatch index from 1990 to 2015. (a) The spatial mismatch index between blue water footprint and oasis farmland area (SMI (WFblue_OFa)); (b) The spatial mismatch index between green water footprint and oasis farmland area (SMI (WFgreen_OFa)); (c) The spatial mismatch index between crop production water footprint and oasis farmland area (SMI (WFpro_OFa)).
Water 13 00696 g005
Figure 6. The spatial-temporal pattern of sensitivity from 1990 to 2015.
Figure 6. The spatial-temporal pattern of sensitivity from 1990 to 2015.
Water 13 00696 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Long, A.; Yu, J.; Deng, X.; He, X.; Gao, H.; Zhang, J.; Ren, C.; Du, J. Understanding the Spatial-Temporal Changes of Oasis Farmland in the Tarim River Basin from the Perspective of Agricultural Water Footprint. Water 2021, 13, 696. https://doi.org/10.3390/w13050696

AMA Style

Long A, Yu J, Deng X, He X, Gao H, Zhang J, Ren C, Du J. Understanding the Spatial-Temporal Changes of Oasis Farmland in the Tarim River Basin from the Perspective of Agricultural Water Footprint. Water. 2021; 13(5):696. https://doi.org/10.3390/w13050696

Chicago/Turabian Style

Long, Aihua, Jiawen Yu, Xiaoya Deng, Xinlin He, Haifeng Gao, Ji Zhang, Cai Ren, and Jie Du. 2021. "Understanding the Spatial-Temporal Changes of Oasis Farmland in the Tarim River Basin from the Perspective of Agricultural Water Footprint" Water 13, no. 5: 696. https://doi.org/10.3390/w13050696

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