Unveiling the Dynamics and Influence of Water Footprints in Arid Areas: A Case Study of Xinjiang, China

: A prerequisite for the rational development and utilization of regional water resources is the measurement of water stress. In this study, from the perspective of water footprints, we took the proportion of the agricultural water footprint within the total water resource usage of Xinjiang (hereafter referred to as XJ) as an example to measure its water stress index and explore the state of water stress in the region and its corresponding driving factors. The ESDA method was applied to characterize the spatial patterns of and correlations with water stress. The effects of different factors on the spatial differentiation between the water footprint and water stress were quantified using the LMDI and geoprobes, respectively. The results showed that (1) both the agricultural water footprint and the water stress index in XJ showed an upward trend, the spatial distribution of water stress was uneven, and the regional pressure difference between the east and the west was greater than that between the north and the south; (2) the water stress index has an obvious negative spatial correlation, fluctuations in its discrete nature have been enhanced, and the number of spatially correlated prefectures is decreasing; (3) water consumption of CNY 10,000 GDP, GDP per capita, and total CO 2 emissions have the most significant impact on the evolution of the agricultural water footprint in XJ. Meanwhile, spatial variations in water stress are mainly determined by the area of cultivation, the area of natural oasis, and the proportion of water used in agriculture. Analysis of the characteristics of and factors influencing water stress in XJ from the perspective of its agricultural water footprint provides a new perspective for further analyzing the actual state of the water footprint and water stress in XJ and supplies a reference basis for the decision-makers of the XJ government.


Introduction
Water is the source of life.Rapid population growth and increasing human production and consumption activities have made long-valuable freshwater resources even more scarce, and conflict between the supply of and demand for water resources is intensifying, with increasing pressure on the water resource system, especially in arid and semi-arid areas [1].In China, agriculture accounts for 74.2% of its total water use, and China's water stress will face unprecedented challenges in the next 20 years as the demand for water resources are further developed and water is increasingly utilized [2].XJ is a typical, inland arid region in northwestern China, and also a typical region where socio-economic development is dominated by agricultural production; therefore, water resource exploitation in agricultural production in XJ affects the ecosystem [3].Highly water-intensive and low-income agriculture and animal husbandry industries bring very significant demand for water resources and ecological pressure [4,5].How much water is used in agriculture and whether it is used rationally or not directly affect the degree of water stress, and the degree of water stress also constrains the exploitation of agricultural water resources.Therefore, adopting the agricultural water footprint as a measure of agricultural water consumption, evaluating its pressure in terms of agricultural water consumption versus regional broad water resources [6], and studying its spatial autocorrelation characteristics and the strength of the role of influencing factors are of great practical significance for exploring the evolutionary driving force of the pressure on water resources and the optimal allocation of water resources.
The concept of water stress was first proposed in 1992 and usually refers to the pressure on the use of water resources in a certain region [7].With society's concern for sustainable development and the measurement of water stress by the World Resources Institute (WRI), the study of water stress has gradually attracted the attention of scholars both at home and abroad.There has therefore been a concomitant gradual increase in the measurement and analysis of water stress on different scales [8,9].Pfister, S. et al. [10] developed a method for assessing the environmental impact of freshwater consumption.This method can be used within most of the existing life-cycle impact assessment (LCIA) methods.Their research found that water consumption may dominate the aggregated life-cycle impacts of cotton textile production in arid regions.Therefore, consideration of water consumption is crucial in LCA studies that include water-intensive products, such as agricultural goods.Gai, L. [11] estimated the water footprint of production and water stress in China during the years from 1985 to 2009.The results showed that high to severe water stress existed mainly in mega-cities and agricultural areas located in the downstream areas of the Yellow River and the Yangtze River in North and Central China.Zeng, Z. [12] developed a simple approach to assessing water scarcity considering both water quantity and quality.The results showed that Beijing has made very significant progress in mitigating water scarcity, and these achievements were made possible by the intensive efforts of water-saving measures and wastewater treatment.With the continuous development of and improvements in water stress assessment methodologies, the supplydemand ratio method has become more widely used among the various measurement methods [13,14].Currently, there has been a paradigm shift in the field of global water consumption, the core element of which is the water footprint theory [15].The introduction of the water footprint theory has provided new methods and ideas for analyzing the real demand, occupancy, and pressure measurement of water resources in agriculture from a socio-hydrological perspective [16].It can accurately reflect the use of water resources in agricultural production, and thus the water footprint has been widely incorporated into water stress measurement as an important indicator of water consumption [17].
Currently, most of the water stress studies undertaken in XJ have focused on the construction of an evaluation index system using physical water and related data to evaluate and analyze the vulnerability [18], security [19], and carrying capacity [20] of its water resources; however, there have been fewer studies on the spatial distribution of water stress and the factors affecting it based on accounting for the agricultural water footprint.To compensate for the lack of such studies in XJ, this paper classifies crop water footprints (including blue, green, and gray water) and animal water footprints from the perspective of social water cycle fluxes and takes the total amount as the actual consumption of water resources.In order to explore the real pressures faced by the water resource system, this study used the ratio of the real consumption of regional water resources to the amount of water resources in the region to measure the water stress index.Meanwhile, the ESDA method was used to characterize the spatial patterns of and correlations with water stress in XJ's cities and towns from 2000 to 2020, geoprobes were used to explore the strength of the role of the spatial distribution characteristics of water stress, and the LMDI was used to analyze the driving force of the effects of various factors on the evolution of the water footprint.The driving mechanism of the water footprint and water stress evolution and a sustainable development strategy to mitigate water stress in XJ are proposed.The results of this study can provide references and ideas for scientific planning and the sustainable management of water resources in XJ, as well as for the formulation of future agricultural development policies.

Study Site Description
The XJ region is located in northwestern China.It is located in the hinterland of the Eurasian continent, with a geographic position of 34 The XJ region is characterized by an extensive distribution of mountains, deserts, oases, glaciers, and snow [21].The study area covers an area of about 1.66 × 10 6 km 2 .There are 14 regions, namely Urumqi, Altay, Bozhou, Changji, Karamay, Tacheng, Yili, Aksu, Bazhou, Kashgar, Hotan, Kexu, Hami, and Turfan.The geographical factors of the Tianshan Mountains and the human factors of XJ have also divided these regions into the northern, southern, and eastern border areas (Figure 1) [22].XJ is far from the ocean and has a temperate continental arid climate, with average annual hours of sunshine of about 2800 h.The average annual evaporation potential ranges from 1600 mm to 2300 mm, which is much more than the average annual precipitation (157.4 mm) [23].The total water resources of XJ in 2020 were 83.16 × 10 9 m 3 , of which 78.77 × 10 9 m 3 were surface water resources, 50.35 × 10 9 m 3 were underground water resources, and 45.96 × 10 9 m 3 were double-counted between surface water and underground water.Its local primary industry (agriculture, forestry, animal husbandry, and fishery) production water consumption was 49.61 × 10 9 m 3 , which dominated (93.4%) in its total water consumption (53.11 × 10 9 m 3 ) [24].
method was used to characterize the spatial patterns of and correlations with water stress in XJ's cities and towns from 2000 to 2020, geoprobes were used to explore the strength of the role of the spatial distribution characteristics of water stress, and the LMDI was used to analyze the driving force of the effects of various factors on the evolution of the water footprint.The driving mechanism of the water footprint and water stress evolution and a sustainable development strategy to mitigate water stress in XJ are proposed.The results of this study can provide references and ideas for scientific planning and the sustainable management of water resources in XJ, as well as for the formulation of future agricultural development policies.

Study Site Description
The XJ region is located in northwestern China.It is located in the hinterland of the Eurasian continent, with a geographic position of 34°09′ N~49°08′ N, 73°25′ E~96°24′ E. The XJ region is characterized by an extensive distribution of mountains, deserts, oases, glaciers, and snow [21].The study area covers an area of about 1.66 × 10 6 km 2 .There are 14 regions, namely Urumqi, Altay, Bozhou, Changji, Karamay, Tacheng, Yili, Aksu, Bazhou, Kashgar, Hotan, Kexu, Hami, and Turfan.The geographical factors of the Tianshan Mountains and the human factors of XJ have also divided these regions into the northern, southern, and eastern border areas (Figure 1) [22].XJ is far from the ocean and has a temperate continental arid climate, with average annual hours of sunshine of about 2800 h.The average annual evaporation potential ranges from 1600 mm to 2300 mm, which is much more than the average annual precipitation (157.4 mm) [23].The total water resources of XJ in 2020 were 83.16 × 10 9 m 3 , of which 78.77 × 10 9 m 3 were surface water resources, 50.35 × 10 9 m 3 were underground water resources, and 45.96 × 10 9 m 3 were doublecounted between surface water and underground water.Its local primary industry (agriculture, forestry, animal husbandry, and fishery) production water consumption was 49.61 × 10 9 m 3 , which dominated (93.4%) in its total water consumption (53.11 × 10 9 m 3 ) [24].

Meteorological, Water Resource, and Statistical Data Sources
The basic data for this study included meteorological data, water resource data, and statistical data, of which the meteorological data (precipitation, air temperature, average wind speed, sunshine hours, and relative humidity) were obtained from the data of 88 meteorological stations provided by the China Meteorological Data Sharing Service Network (http://data.cma.cn/,accessed on 1 August 2021).The hydrological data (water resources, primary industry water consumption, domestic water consumption, and total water consumption) were obtained from the XJ Water Resources Bulletin (2001-2020) provided by the Water Resources Bureau of XJ Uygur Autonomous Region.The statistical data (resident population number, GDP, value added of CNY 10,000 of industry, cultivated area, meat production of livestock, fertilizer application, etc.

Methodology 2.3.1. Agricultural Water Footprint
In this study, based on the water footprint theory of Hoekstra [25] and others, we calculated the agricultural water footprint (AWF) in XJ, including the crop water footprint (WF crop ) and the animal water footprint (WF animal ).The calculation formula is as follows: where AWF is the agricultural water footprint (m 3 ); WF crop is the crop water footprint (m 3 ); and WF animal is the animal water footprint (m 3 ).
(1) Crop Water Footprint The water footprint of crops is the amount of water resources consumed by crops in the course of their growth.Depending on the type of water consumed and its impact on the environment, it can be further classified into a blue water footprint, green water footprint, or gray water footprint, which are calculated using the following formulas: where WF crop is the water footprint of the crops (m 3 ); WF crop−blue is the blue water footprint of the crops (m 3 ); WF crop−green is the green water footprint of the crops (m 3 ); and WF crop−grey is the gray water footprint of the crops (m 3 ).In this paper, major crops were selected for the water footprint calculation based on the XJ Uygur Autonomous Region Statistical Yearbook and the XJ Production and Construction Corps Statistical Yearbook.The crops included 15 crops such as rice, wheat, coarse cereals, soybeans, cotton, oil plants, sugar beets, vegetables, melons, potatoes, alfalfa, grapes, apples, fragrant pears, and red jujube.According to the statistics, it can be seen that the selected crops accounted for 88.7% of the total cultivated area in XJ.Therefore, in this paper, the sum of the water footprints of the above 15 crops was taken as the blue-green water footprint of the crops in XJ.Firstly, this paper calculated the evapotranspiration (ET c ) and effective rainfall (P e f f ) during the reproductive period of each major crop in the study area.The Penman-Monteith method, recommended by the Food and Water Agriculture Organization of the United Nations (FAO), was used to calculate ET c (mm), which can accurately calculate ET c under different regions and climatic conditions [26].P e f f (mm) was calculated using the CROPWAT model recommended by the FAO.Secondly, the blue water footprint and green water footprint of the different crops in each region were calculated separately based on the source of water use using the formulas shown below [27]: Water 2024, 16, 1164 5 of 19 where P i is the total yield (t) of crop i. Y refers to the crop yield per unit area (t/ha).Coefficient 10 is a coefficient for converting water depth into water volume per unit area of land area.∑ is the cumulative amount of blue (or green) water from the planting period to the harvesting period.lgp is the length of the growing season (d).
The gray water footprint of crops mainly refers to the amount of water required to dilute underutilized fertilizer to its standard concentration after its application through leaching it into rivers via irrigation water or rainfall, which produces polluting ions.In this paper, we analyzed the gray water footprint of the crops using nitrogen and phosphorus fertilizers, which are the most applied in XJ, as the main sources of polluting ions.The same water body can dilute multiple pollutants at the same time, so the gray water footprint is determined by the amount of water required to dilute the largest pollutant ion.The calculation formula is shown below [28]: where The animal water footprint is defined as the water footprint of each animal throughout its life cycle from birth to slaughter, including water consumed in the form of virtual and direct water for growth, processed feed, drinking water, and washing services, calculated using the formula shown below [29]: where WF animal is the animal water footprint (m 3 ); WF ani− f eed is the animal feed water footprint (m 3 ); WF ani−drink is the animal drinking water footprint (m 3 ); and WF ani−serve is the animal service water footprint (m 3 ).In this paper, the main animal products were selected according to the XJ Uygur Autonomous Region Statistical Yearbook and the XJ Production and Construction Corps Statistical Yearbook for the purposes of the water footprint calculation.The animal products included cattle, horses, camels, pigs, goats, sheep, chickens, and rabbits; these comprised eight categories of livestock, using live animal data to remove the growth and feed processing water used and the overlap with the water footprint of the feed crops.The calculation of the formula is shown in the following equation [30]: where Q mix denotes the amount of water required for mixed feed (m 3 •d −1 ); C ij is the weight of feed j consumed by the ith animal (t•d −1 ); and SWC j denotes the amount of virtual water contained in feed i (m 3 •t −1 ).The virtual water volume of the feed crops was calculated using the crop water footprint method, weighted according to the ration (weight) of different feed crops; Q i d and Q i s are the amount of water required as drinking water and for washing services for animal i per day (m 3 •d −1 ), respectively.

Water Resource Pressure Assessment Index
The water resource pressure index (WRPI) reflects the degree of pressure on regional water resources, characterizing the load on water resources from human production and living processes and the relative scarcity of regional water resources.In order to explore the actual pressure on the water resources in each region of XJ, the actual consumption of water resources was measured based on the results of regional water footprint calculations, and the ratio of the actual consumption of regional water resources (AWF) to the total water resources in a broad sense (WA all ) was used to express the pressure on the water resources.The calculation formula is shown below: where WRPI is the water resource stress index; AWF is the total water footprint of regional agricultural production (m 3 ); and WA all is the total water resources of the region in a broad sense (m 3 ).When WRPI > 1, that means that the water consumption of agricultural production has exceeded the amount of available water resources, and the larger the value, the more serious the local water resource pressure; similarly, when WRPI ≤ 1, the smaller the value, the smaller the water resource pressure.

Exploratory Spatial Data Analysis
Exploratory spatial data analysis (ESDA) is a method used to analyze spatial data and express them visually, which can intuitively reveal characteristics of spatial data such as their correlation, aggregation, and hot/cold spots [31].Global spatial autocorrelation is usually chosen to characterize the spatial correlation of the geographic phenomena across an entire region, which reflects the overall trend in the spatial correlation of the observed variables across the whole study area.The commonly used measure is global Moran's I statistic, and the I statistic index takes a value between −1 and 1.The formula is shown below [32]: where i ̸ = j, N is the number of study objects; X i is the observation; X is the mean value of W ij ; and W ij is the spatial weighting matrix between the study objects i and j, with a spatial proximity of 1 and non-proximity of 0.
Although global spatial autocorrelation analysis can reveal the degree of dependence on objects as a whole, it ignores possible local smoothing; therefore, local spatial autocorrelation must be chosen to reveal the autocorrelation of the local regional units in the neighboring space and to reflect the degree of spatial agglomeration of a certain element in detail, which is represented by a LISA agglomeration map.In this paper, local Moran's I index was used to measure the heterogeneity of the spatial elements between the regional units i and j, and the calculation formula is shown below [33]: where Z i and Z j are the normalization of the observations in spatial cells i and j, respectively; and W ij is the spatial weight.

Logarithmic Mean Divisia Index
Logarithmic Mean Divisa Index (LMDI) decomposition is a complete, residual-free decomposition analysis method part of the Divisia index method [34].Its core idea is to decompose a target variable into a combination of several influencing factor changes so that the magnitude of the influence of each factor can be identified, and then the factors that contribute more can be determined [35].Based on the research purpose of exploring the drivers of agricultural water footprint changes in XJ, this paper decomposed them into a population effect (pop-eff ), life effect (life-eff ), structural effect (struc-eff ), scale effect (scale-eff ), technology effect (tech-eff ), economic effect (econ-eff ), ecological effect (ecol-eff ), and environmental effect (env-eff ).Among them, pop-eff shows the change in p ij , life-eff shows the change in l ij , struc-eff shows the change in st ij , scale-eff shows the change in sc ij , tech-eff shows the change in T ij , econ-eff shows the change in econ ij , ecol-eff shows the change in ecol ij , and env-eff exhibits the change in env ij .According to the LMDI full decomposition model, the formulas for each factor effect are shown below [36,37]: Water 2024, 16, 1164 where ω f ij denotes the weight; p j denotes the population of j area (10 4 ); l j denotes the domestic water consumption of j area (10 8 m 3 ); sc j denotes the area of crop cultivation in j area (10 3 ha); env j denotes the natural oasis area in j area (km 2 ); env j denotes the total carbon dioxide emissions of j area (10 4 t); AW j denotes the agricultural water use in j region (m 3 ); TW j denotes the total water use in j region (m 3 ); Y j denotes the gross domestic product in j region (CNY 10 4 ); and "∆WF" denotes the change in the water footprint of agriculture in XJ.The superscripts "0" and "t" denote the base year and target year of this study, respectively.

Geographical Detector
A geoprobe is a spatial statistical method used to explore the degree of influence of geographic factors on the spatial differentiation of a geographic phenomenon and to analyze the causes of such differentiation.In this paper, the factor detection part is used to explore the influence of different factors on the spatial distribution of water stress in XJ's prefectures and cities.The formula for calculating the degree of influence of geographic factors is shown below [38]: where q is the value of the degree of influence of influencing factors on water stress in XJ; h = 1, 2, . ..; L is the classification number of the influencing factors; N is the number of prefectural and municipal administrative divisions in XJ; N h denotes the number of prefectural and municipal administrative regions classified as h; σ 2 is the variance in the interstate water stress index in XJ's prefectures; and σ 2 h is the variance in the regional water stress index classified as h.The value of q is [0, 1], and the larger the value of q is, the greater the influence of this factor on water stress in XJ. q = 0 means that it is not affected by this factor.The results show that the total AWF increased by 124.0% over the past 20 years (WF crop and WF animal increased by 147.9% and 48.9%, respectively).The regional differences in the annual average value of AWF were obvious; in order from highest to lowest are the southern border (27.74 × 10 9 m 3 ), the northern border (21.67 × 10 9 m 3 ), and the eastern border (2.43 × 10 9 m 3 ).Among them, the AWF at the southern border was 11.42 times higher than that at the eastern border.

Temporal Evolution Analysis of the Agricultural Water Footprint
In the agricultural production process, the WF animal multi-year average was 22.4 per cent of the AWF multi-year average, so it can be seen that WF crop is a major component of AWF (~77.6%).Of the three components that make up WF crop in each region, blue water dominates (74.6%), with relatively small proportions of green water (14.0%) and gray water (11.4%).Similar to WF crop , WF crop−blue showed a gradual increase in all zones from 2000 to 2020, with the southern border contributing the most to the components of WF crop−blue (60.0%), followed by the northern border (34.2%) and the eastern border (5.8%).Compared with WF crop−blue , the interannual changes in WF crop−green and WF crop−grey showed an overall relatively small increase in all regions.
2020.The results show that the total  increased by 124.0% over the past 20 years ( and  increased by 147.9% and 48.9%, respectively).The regional differences in the annual average value of  were obvious; in order from highest to lowest are the southern border (27.74 × 10 9 m 3 ), the northern border (21.67 × 10 9 m 3 ), and the eastern border (2.43 × 10 9 m 3 ).Among them, the  at the southern border was 11.42 times higher than that at the eastern border.

Spatio-Temporal Evolution Analysis of Water Resource Pressure
We calculated the water resource pressure index (WRPI) for each prefecture in XJ from 2000 to 2020.In this paper, due to space limitations, we only selected data for 2000, 2005, 2010, 2015, and 2020 to map the changes in the spatial patterns of the WRPI according to five levels, as follows: lower (0.08 ≤ WRPI < 0.5), low (0.51 ≤ WRPI < 1.0), medium (1.01 ≤ WRPI < 1.5), high (1.51 ≤ WRPI < 2.0), and higher (2.01 ≤ WRPI < 6.0) (Figure 3).In the agricultural production process, the  multi-year average was 22.4 per cent of the  multi-year average, so it can be seen that  is a major component of  (~77.6%).Of the three components that make up  in each region, blue water dominates (74.6%), with relatively small proportions of green water (14.0%) and gray water (11.4%).Similar to  ,  showed a gradual increase in all zones from 2000 to 2020, with the southern border contributing the most to the components of  (60.0%), followed by the northern border (34.2%) and the eastern border (5.8%).Compared with  ,the interannual changes in  and  showed an overall relatively small increase in all regions.

Spatio-Temporal Evolution Analysis of Water Resource Pressure
We calculated the water resource pressure index (WRPI) for each prefecture in XJ from 2000 to 2020.In this paper, due to space limitations, we only selected data for 2000, 2005, 2010, 2015, and 2020 to map the changes in the spatial patterns of the WRPI according to five levels, as follows: lower (0.08 ≤ WRPI < 0.5), low (0.51 ≤ WRPI < 1.0), medium (1.01 ≤ WRPI < 1.5), high (1.51 ≤ WRPI < 2.0), and higher (2.01 ≤ WRPI < 6.0) (Figure 3).With the development of the social economy and the intensification of the frequency of human water extraction activities, the overall water stress in XJ has shown an increasing trend year by year (Figure 3).The results show that from 2000 to 2020, the water resource pressure index increased by 187.1%.Among them, water resources in Turpan (1.83) were under high pressure in 2000, followed by Changji (1.19), Kashgar (0.68), and Aksu (0.53), while the rest of the region had a low water resource pressure index, and the pressure index at the junction of eastern and northern XJ was relatively high.In 2005, the water stress in Karamay (increased by 1.75) and Tacheng (increased by 0.23) increased more than in previous years, resulting in a relatively high stress index in the northern border area; meanwhile, in the southern border area, the situation remained the same as in 2000.In With the development of the social economy and the intensification of the frequency of human water extraction activities, the overall water stress in XJ has shown an increasing trend year by year (Figure 3).The results show that from 2000 to 2020, the water resource pressure index increased by 187.1%.Among them, water resources in Turpan (1.83) were under high pressure in 2000, followed by Changji (1.19), Kashgar (0.68), and Aksu (0.53), while the rest of the region had a low water resource pressure index, and the pressure index at the junction of eastern and northern XJ was relatively high.In 2005, the water stress in Karamay (increased by 1.75) and Tacheng (increased by 0.23) increased more than in previous years, resulting in a relatively high stress index in the northern border area; meanwhile, in the southern border area, the situation remained the same as in 2000.In 2010, the stress index increased in Bortala, Urumqi, Aksu, and Kashgar, and the water stress eased in Karamay and Tacheng.In this period, the water stress was higher in both the eastern (1.04) and southern (0.51) regions than in the northern region (0.37).In 2015, the water resource pressure in southern XJ, northern XJ, and eastern XJ had increased to varying degrees compared to the previous period, and this increase in pressure was also mainly characterized in the different prefectures mentioned above.By 2020, water stress in XJ had increased rather than decreased, reaching the strongest stress level (0.89) in the past 20 years.All regions had indices higher than 1.0 except for the six prefectures of Altay, Urumqi, Ili, Kizilsu, Hotan, and Bayingol, which had relatively low stress levels, with the highest WSI in Karamay (5.82) and the second and third highest being Changji (3.45) and Aksu (2.19), respectively.Overall, over the last 20 years, the WSI in the eastern XJ region has generally been higher than that in the southern and northern XJ regions, while the Karamay, Changji, and Aksu regions have been the regions with higher WSIs among the many prefectures of XJ.These results indicate that the water stress index of XJ's prefecture and municipalities passed the 95% confidence interval assumption within the period from 2000 to 2005, but the Moran index was relatively small, and the spatial positive autocorrelation effect was relatively insignificant.Within the period from 2005 to 2020, the Moran index was less than 0, indicating its negative spatial autocorrelation, i.e., water stress had obvious discrete spatial characteristics.In 2000-2005 and 2010-2015, the degree of aggregation was weakened, and the discrete nature was increasing; meanwhile, in 2005-2010 and 2015-2020, the degree of aggregation was enhanced, and the discrete nature was slightly weakened.There were clear fluctuations in dispersion across the study periods, with the sharpest fluctuations occurring between 2000 and 2005.

Local Spatial Autocorrelation Analysis of Water Resource Pressure in XJ
In order to reveal the discrete spatial phenomenon of water resource pressure in each province, LISA clustering maps of the water resource pressure in XJ in 2000, 2005, 2010, 2015, and 2020 were plotted based on the local Moran's I index and a Moran scatterplot (Figure 4).The results of LISA plots are often categorized into High-High Cluster (H-H), Low-Low Cluster (L-L), High-Low Cluster (H-L), Low-High Cluster (L-H), and not significant, where the HH and LL types both exhibit a spatially positive correlation type.The HL and LH types both exhibit spatially negative correlation types, and not significant indicates no obvious clustering or discrete features.In the blocks of five-year periods, 6, 5, 4, 3, and 3 prefecture cities and municipalities were in a state of spatial positive/negative correlation (HH, HL, or LH clustering), respectively.Overall, they followed a decreasing trend, and the change trend was the same as the G-Moran's I trend, indicating that the discrete spatial status of water resource pressure in XJ's cities and towns showed an intensifying trend over time.In 2005 and 2010, the HH agglomeration was consistent, and in the remaining years, only one prefecture experienced HH agglomeration.In 2000, two prefectures had an HL concentration.For 2005-2020, there was only one geographic area with HL agglomeration.The number of prefectures with LH agglomeration was the largest in those five years, decreasing gradually from three prefectures in 2000 to one prefecture in 2020.
throughout XJ, with Altay, Tacheng, and Bayingol in this state in 2000; then, Altay, Tacheng, and Bayingol transformed from this state to a discrete state in sequence in the next 5 to 15 years until 2020, when only Hami, which is adjacent to Turpan, was transformed from a discrete state to an LH agglomeration state.Over the past 20 years, with the transformation of agglomeration-disagglomeration states among XJ's prefectures, water stress agglomeration has weakened, and fluctuations in disagglomeration have increased.Moreover, the number of spatially correlated prefectures has decreased, the resource mismatch has increased, and the differences in the water resource background conditions and water use between prefectures have gradually become larger.

Analysis of Factors Influencing the Evolution of the Production Water Footprint and Water Resource Pressure
Considering the comprehensive and accessible nature of the data, four influencing factors, comprising people's lives, agricultural development, economy and technology, and the ecological environment, were selected, including the following relationships between effects: a total population-population effect (X1); a domestic water consumptionlife effect (X2); a proportion of agricultural water consumption-structure effect (X3); a cultivation area-scale effect (X4); a water consumption of CNY 10,000 GDP-technology effect (X5); a GDP per capita-economic effect (X6); a natural oasis area-ecological effect (X7); and a total carbon dioxide emissions-environmental effect (X8).The LMDI method was used to decompose the contribution rate of each factor to the change in the water footprint, and the geoprobe method was used to explore the factors influencing the regional differences in XJ's water stress index and analyze the intensity of their roles.

Contribution Value of the Agricultural Water Footprint's Driving Factors in XJ
In order to further explore the factors influencing the agricultural water footprint in XJ, the contribution value of each driver of the water footprint in XJ was calculated with 2000 as the base year and 5 years as the time period (Figure 5).The total effect of each factor in XJ from 2005 to 2020 has always been positive, but the overall trend has declined, The water stress HH catchment area is mainly located at the junction of the eastern and northern borders, principally including the cities of Turpan and Urumqi, forming an area with high values for the water stress index.The region began to expand in 2005 and then shrunk to only one place in 2015, mainly due to a shortage of precipitation, which resulted in high water stress.This has become a resource shortcoming in the process of socio-economic development.Water stress HL agglomeration is mainly distributed in the western part of the southern border, and there are only a few prefectures in this state of water stress, including Kashgar and Aksu in 2000; meanwhile, only Kashgar remained in the HL agglomeration state after 2005.Water stress LH agglomeration is scattered throughout XJ, with Altay, Tacheng, and Bayingol in this state in 2000; then, Altay, Tacheng, and Bayingol transformed from this state to a discrete state in sequence in the next 5 to 15 years until 2020, when only Hami, which is adjacent to Turpan, was transformed from a discrete state to an LH agglomeration state.Over the past 20 years, with the transformation of agglomeration-disagglomeration states among XJ's prefectures, water stress agglomeration has weakened, and fluctuations in disagglomeration have increased.Moreover, the number of spatially correlated prefectures has decreased, the resource mismatch has increased, and the differences in the water resource background conditions and water use between prefectures have gradually become larger.

Analysis of Factors Influencing the Evolution of the Production Water Footprint and Water Resource Pressure
Considering the comprehensive and accessible nature of the data, four influencing factors, comprising people's lives, agricultural development, economy and technology, and the ecological environment, were selected, including the following relationships between effects: a total population-population effect (X 1 ); a domestic water consumption-life effect (X 2 ); a proportion of agricultural water consumption-structure effect (X 3 ); a cultivation area-scale effect (X 4 ); a water consumption of CNY 10,000 GDP-technology effect (X 5 ); a GDP per capita-economic effect (X 6 ); a natural oasis area-ecological effect (X 7 ); and a total carbon dioxide emissions-environmental effect (X 8 ).The LMDI method was used to decompose the contribution rate of each factor to the change in the water footprint, and the geoprobe method was used to explore the factors influencing the regional differences in XJ's water stress index and analyze the intensity of their roles.

Contribution Value of the Agricultural Water Footprint's Driving Factors in XJ
In order to further explore the factors influencing the agricultural water footprint in XJ, the contribution value of each driver of the water footprint in XJ was calculated with 2000 as the base year and 5 years as the time period (Figure 5).The total effect of each factor in XJ from 2005 to 2020 has always been positive, but the overall trend has declined, and 2015 was the turning point for the trend in the total value of the total effect.From a structural point of view, the absolute value of the share of each factor in the amount of change in the agricultural water footprint is as follows, in descending order: technical effect (−3.35) > economic effect (2.09) > environmental effect (0.95) > scale effect (0.54) > life effect (0.43) > population effect (0.37) > ecological effect (−0.14) > structural effect (0.09).Among them, the technology effect and economic effect accounted for a much larger proportion than the contributions of the other factors, with mean values of −33.15 × 10 9 m 3 and 20.70 × 10 9 m 3 , respectively.Among the many factors, the economic effect, the environmental effect, and the population effect have always promoted the growth of the agricultural water footprint; on the contrary, the technology effect and the ecological effect always inhibit it, while the life effect, the structure effect, and the scale effect play different individual roles in the growth of the water footprint at different stages.Overall, the influence of the technology and economic effects dominates the process of growth in the agricultural water footprint in XJ.

Intensity of Factors Influencing Water Resource Pressure in XJ
Firstly, the natural breakpoint method was applied to divide the eight indicators into 10 levels, and then geodetic factor detection was used to calculate the intensity of the role of each indicator (q-value) in order to obtain the calculation results (Table 1).The main indicators affecting the spatial differentiation of the water stress index did not change significantly in different periods; however, the intensity of the role of each indicator factor changed significantly.Among them, the indicator with the highest-intensity role in the spatial differentiation of the water resource pressure in XJ was the cultivated area (0.665) in 2000, and per capita GDP (0.223) had the lowest-intensity role in the same year.Over time, each indicator showed an increasing or decreasing trend in the overall intensity of its effect on the spatial differentiation of water resource pressure during the 20-year period.The indicators showing an increasing trend in descending order were domestic water use (135.3%)> GDP per capita (93.7%) > total CO2 emissions (29.2%) > area of natural oases (9.4%).Meanwhile, the indicators showing a decreasing trend in descending order were water use per CNY 10,000 of GDP (−18.8%) > total population (−10.3%)> share of water used in agriculture (−2.4%) > area under cultivation (−0.2%).
In 2000, specifically, the most important factor contributing to the spatial differentiation of water stress was the area under cultivation, followed by the total population and then the proportion of water used in agriculture, with GDP per capita playing the weakest role.In 2005, the main factor shifted from the area under cultivation to the total population, with the intensity of the area under cultivation decreasing over the five-year period after 2005, while the intensity of the natural oasis areas increased significantly, and the intensity of total CO2 emissions was the lowest in the same year.In 2010, identical to the The population effect of the change in the agricultural water footprint in XJ generally ranges from 3.01× 10 9 to 5.38 × 10 9 m 3 , showing a fluctuating upward trend; this indicates that although the population effect has had a relatively small impact on the growth of the water footprint in XJ, it has always played a facilitating role.The livelihood effect was a reverse driver of the water footprint's growth before 2005 but became a facilitator after 2005, with an average annual contribution of 4.21 × 10 9 m 3 to the agricultural water footprint.As with the population effect, the livelihood effect has had a relatively small impact on the water footprint changes in XJ.Among the many factors, the structure effect has had the smallest impact on the change in the agricultural water footprint in XJ, with an average annual contribution value of only 0.9 × 10 9 m 3 .This indicates that although there have been changes in the water use structure in XJ, their overall effect is not obvious.The scale effect is similar to the livelihood effect, which positively drove the water footprint before 2015 but inhibited it after 2015.The average contribution of the scale effect to the changes in XJ's agricultural water footprint declined most significantly overall, by −2.2 percent.The technology effect has driven the increase in XJ's agricultural water footprint in the reverse direction.The average annual contribution of the technology effect to the water footprint from 2005 to 2020 was as high as −33.15 × 10 9 m 3 , with a mean contribution value of −335%.This indicates that the efficiency of water use in XJ improved across the period, and water resources were developed and utilized in a more intensive and efficient way than in the base period.The economic effect is the main factor that has caused the growth in the agricultural water footprint in XJ.The economic effect was a positive driver from 2005 to 2020, with an average annual contribution value of 20.7 × 10 9 m 3 and a mean contribution value of 209%.In addition, it is the higher share and rapid growth of the economic effect that has caused the total water footprint of XJ to continue to show an upward trend despite the continuous improvements in water use efficiency.The ecological effect is generally in the range of −0.33 × 10 9 ~−2.02 × 10 9 m 3 , which means it has not played an obvious role; however, in the same way as the technological effect, the ecological effect has always inhibited the growth of the agricultural water footprint in XJ.The environmental effect shows a fluctuating upward trend in the period 2005-2020, with an average annual contribution of 9.44 × 10 9 m 3 to the growth in the agricultural footprint, contributing 95% of the mean value.This makes it the third most influential factor after the technological and economic effects.

Intensity of Factors Influencing Water Resource Pressure in XJ
Firstly, the natural breakpoint method was applied to divide the eight indicators into 10 levels, and then geodetic factor detection was used to calculate the intensity of the role of each indicator (q-value) in order to obtain the calculation results (Table 1).The main indicators affecting the spatial differentiation of the water stress index did not change significantly in different periods; however, the intensity of the role of each indicator factor changed significantly.Among them, the indicator with the highest-intensity role in the spatial differentiation of the water resource pressure in XJ was the cultivated area (0.665) in 2000, and per capita GDP (0.223) had the lowest-intensity role in the same year.Over time, each indicator showed an increasing or decreasing trend in the overall intensity of its effect on the spatial differentiation of water resource pressure during the 20-year period.The indicators showing an increasing trend in descending order were domestic water use (135.3%)> GDP per capita (93.7%) > total CO 2 emissions (29.2%) > area of natural oases (9.4%).Meanwhile, the indicators showing a decreasing trend in descending order were water use per CNY 10,000 of GDP (−18.8%) > total population (−10.3%)> share of water used in agriculture (−2.4%) > area under cultivation (−0.2%).In 2000, specifically, the most important factor contributing to the spatial differentiation of water stress was the area under cultivation, followed by the total population and then the proportion of water used in agriculture, with GDP per capita playing the weakest role.In 2005, the main factor shifted from the area under cultivation to the total population, with the intensity of the area under cultivation decreasing over the five-year period after 2005, while the intensity of the natural oasis areas increased significantly, and the intensity of total CO 2 emissions was the lowest in the same year.In 2010, identical to the position in 2000, the cultivated area was still the most important factor, while the intensity of domestic water consumption was significantly higher than the other factors.The intensity of the population, agricultural water use, and natural oasis area had decreased significantly, and the weakest-intensity role was played by GDP per capita.In 2015, the cultivated area was slightly more important than it had been five years previously, and the factors with an intensity greater than 0.6 were oasis area, domestic water use, and agricultural water use, with CO 2 emissions having the lowest intensity.In 2020, although the intensity index had decreased slightly, the main influencing factor was still the planted area, while the intensity of total CO 2 emissions had risen significantly and was one of the top three factors, together with natural oasis areas and planted areas, and the intensity of water consumption per CNY 10,000 GDP reached its lowest value in 20 years and was also the weakest factor in that year.
In conclusion, although the main factors affecting the spatial differentiation of water stress have shifted between cultivated areas and population in recent years, it is undeniable that the intensity of the effect of cultivated on the spatial differentiation of water stress has remained high over the past 20 years.In addition to the number of people, the proportion of agricultural water use and the area of natural oases continue to have a greater impact on the spatial differentiation of water stress.In general, for XJ, the expansion of its irrigation scale, its population growth, the imbalance of the industrial (water use) structure, and the evolution of natural-artificial oases had the greatest influence on the spatial differentiation of water stress, while GDP per capita had the weakest influence, followed by water use per CNY 10,000 GDP.

Key Factors Affecting Changes in the Agricultural Water Footprint and Water Stress and Suggestions for Addressing Them
As shown by the historical development trend, AWF shows a continuous growth trend in terms of spatial and temporal changes (WF crop and WF animal are also in these states).If corresponding regulatory policies are not implemented, the growth of WF crop−blue , driven through the expansion of the cultivation area and irrigation scale, will lead to a reduction in natural oasis areas, which will further lead to the obvious decline in biohabitat quality, water productivity, and the carbon sequestration capacity of vegetation [39].On the other hand, an increase in agricultural fertilizer's application, due to WF crop−grey 's growth, will not only place enormous pressure on water use for the dilution of hazardous chemicals, affecting the surface and groundwater quality to some extent, but will also cause a deterioration in soil quality and a decline in crop yields due to the continued accumulation of salts in the soil [40].If there is a decline in the marginal effects of the objective input costs (e.g., agricultural mechanization and cropping systems, etc.), the scale of agricultural (plantation and livestock) production in XJ will expand further.Therefore, from the perspective of sustainable development in the future, the government has proposed policies to address this issue in order to maintain an appropriate scale of agricultural development.Under this policy, WF crop−blue , WF crop−grey , and WF animal will decrease with the gradual implementation of the policy.In addition, in future sustainable development scenarios, the government should pay special attention to the phenomenon of the erosion of natural oases through the growth of artificial oases, and avoiding the severity of this phenomenon as much as possible will effectively improve the problem of deteriorating habitat quality.Existing studies have shown that the implementation of such management measures can increase the overall ecosystem service value of the region [41,42].
The industrial development of the XJ region is dominated by agriculture, and thus agricultural water use occupies a pivotal position in the process of water use for socioeconomic development in XJ [43].The impact of the share of agricultural water use on water stress in XJ is significant, and thus an accurate evaluation of agricultural water use in XJ's states is fundamental to the evaluation of water stress.The use of water footprinting methods to assess agricultural water use efficiency and water stress has been widely adopted worldwide [44].However, it has been demonstrated during the development and utilization of existing water resources that the implementation of water-saving technologies will improve the water use efficiency of decision-makers; however, the development of water-saving irrigation systems, the optimization of water allocation in irrigation districts, and the development of irrigation water demand forecasting technologies may lead to a continued high level of socio-ecological water demand and use, a further crowding-out of natural vegetation, and increased water consumption by water bodies, thus illustrating the "Jevons paradox" of water saving applied in the agricultural context [45].
The degree of the match between irrigation water use and the distribution pattern of effective irrigated areas is one of the most important indicators affecting the water stress index.The amount of water resources consumed throughout the economic production process in XJ is more closely related to the scale of crop cultivation in terms of time and space.Thus, we believe that the most effective way to rationally and optimally allocate water and alleviate water stress would be to focus on the spatial matching of water and land resources and their sensitivity [46].In short, water resource management policy-makers can alleviate the structural imbalance between water resources and land resources through measures such as promoting the transfer of agricultural land and changing the scale of agricultural land management, thereby alleviating the existing water resource pressure across XJ.On the other hand, in terms of cropping structure, farmers in all regions of XJ tend to grow cash crops (e.g., cotton and jujube) rather than food crops in order to maintain their smooth income and standard of living.As a result, the large-scale production of high water-consuming cash crops further increases the pressure on local water resources, so it is important to encourage farmers to increase the proportion of low-water-consuming crops and to grant or increase subsidies to farmers.
Therefore, in the implementation of future development policies, the government should focus on China's "One Belt, One Road" initiative to optimize and develop the product structure of the processing and service industries, and through adjustments to the industrial structure, it can effectively reduce the ineffective consumption of water resources and ensure the high-quality development of the socio-economics of the whole region, so as to achieve a win-win situation [47].From a macro perspective, governors should pay more attention to improving agricultural production while maintaining a benign ecological environment, developing better resource carrying capacity assessment tools to accurately assess the appropriate carrying capacity of each region, and improving the integrated management of water and land resources [48].From an engineering perspective, the government can alleviate the current imbalance between agricultural development and water use through the implementation of inter-regional water transfer projects and integrated planning for planting scale control, thereby reducing the water stress index.

Discussion on the Driving Mechanism of Agricultural Water Footprint and Water Resource Pressure
Through this study and analyses of its drivers, we have obtained preliminary results that explore the spatial evolution trend in the agricultural water footprint, water stress, and its main drivers in XJ from 2000 to 2020.Firstly, the three major factors that make the highest contribution to the change in the agricultural water footprint in XJ are the technological effect, the economic effect, and the environmental effect.Among them, the technological effect plays an inhibiting role, while the economic and environmental effects have instead positively driven the growth in the water footprint.Secondly, the three most important factors contributing to the spatial differentiation in water stress in XJ are the cultivated areas, the natural oasis areas, and the agricultural water use ratio.Although the intensity of these drivers has been quantified, the mechanisms driving the spatial evolution of water stress are not yet fully understood.
The emergence and evolution of anything are closely related to its original developmental basis.New water stress patterns generally represent the inheritance and development of existing patterns.Therefore, the basis of historical development is one of the main driving mechanisms for the differences in water stress among XJ's states.From 2000 to 2020, although there was a certain degree of randomness and instability in the spatial distribution of the water resource pressure across XJ, the overall trend showed that it was concentrated in the western part of the southern border and the border between the eastern and northern parts of the border.One of the reasons for this is that the prefecture and municipalities in this region have always been constrained by its water resource conditions in various ways, such as the degree of development of its irrigation districts, the characteristics of its industrial structure, and obvious disadvantages in terms of water vapor's abundance and scarcity.
Regional geography is one of the important factors that influences the evolution of regional water stress patterns.Through its influences on the sequence of and opportunities for regional water resource development, it is related to the development speed of regional water resource development and utilization and the status of the local spatial autocorrelation of water pressure.From 2000 to 2015, in terms of the spatial autocorrelation of the water resource pressure differences between the cities and towns in XJ, HL agglomeration and LH agglomeration were always mainly concentrated in the southern XJ region, and regional water resources development and use perform relatively poorly in terms of overall correlation.The ecological environment of the southern border region is relatively fragile, the degree of exploitation and utilization has gradually increased but in a more relaxed manner, and the main water resource background conditions do not have any significant disadvantages compared with those of the northern border.However, the sufficient and positive conditions for crop cultivation and production have led to the rapid development of the scale of agricultural irrigation in the southern XJ region, and its industrial structure is far less reasonable than that in other regions of XJ.This has further led to a massive outflow of local physical water resources to the rest of the country in the form of virtual water, resulting in an inefficient economic use of resources and ultimately creating a significant gap between the pressure on water resources, the competitive environment of economic development, and the dynamism of economic development here in comparison with other regions.
Regional development policy is an important driving force in the evolution of regional water stress differentials and plays a very important role in shaping regional water stress patterns.In a sense, regional development policy actually reflects a specific priority, which is the government's intention and goals for regional development, i.e., balanced or unbalanced regional development, formed after comprehensive research and analysis.From 2001 to 2010, the state proposed the Western Development Strategy, and XJ was one of the key regions in the implementation of this strategy because it was rich in resources but its people's income level was generally low [49].After the 16th Party Congress in 2002, the Party Central Committee and the State Council agreed that XJ's development required national input and thus implemented a series of poverty alleviation programs and deployments focusing on farmland water conservancy in XJ.However, this directly contributed to the growth in crop acreage, and the significant increase in the agricultural water footprint has further increased the proportion of water used in agriculture, exacerbating the difficult water resource carrying situation.After 2010, the implementation of the XJ aid policy further enabled XJ in its growth and development from a previously weak industrial base.With the support of other provinces across the country, its development has been further focused on arable land reclamation and comprehensive agricultural development, especially in the southern XJ region, in the form of constructing water-efficient zones and developing eco-agricultural demonstration parks [50].Up to the present day, the effect of such policy factors on the formation of spatial patterns in water resource use continues to be amplified.

Conclusions
This paper took the XJ region as a case study.Working from the perspective of its agricultural water footprint, the regional water stress index was accounted for in the results of an agricultural water footprint assessment based on the past 20 years.Trend analysis and ESDA were applied to analyze the characteristics of the spatial pattern evolution and spatial correlation of the water stress index in XJ's prefectures and cities from 2000 to 2020.The LMDI and geodetectors were used to detect the factors influencing the evolution in the agricultural water footprint and the spatial differentiation in water stress, respectively.The conclusions are as follows: (1) From 2000 to 2020, XJ's agricultural water footprint was in a state of continuous growth, with the total amount more than doubling.Among the components of the agricultural water footprint, the crop water footprint increased more than the animal water footprint, and significant growth has occurred mainly in the southern and northern regions.During the study period, changes in the water stress in XJ's cities and towns were unstable, with super-high water stress indices concentrated in the eastern XJ region, and the overall trend was upward.The spatial distribution of the water stress index was uneven, concentrated in the west and east of XJ, and the regional differences between the north and the south were smaller than those between the east and the west.Regions such as Karamay, Changji, Turpan, Aksu, and Kashgar had high multi-year average water stresses.(2) The spatial correlation of the water resource pressure index is obvious, its negative correlation is significant, and the intensity of the discrete state shows fluctuating changes.The water resource pressure HH agglomeration area is mainly concentrated in the junction zone of east XJ and north XJ, and the HL agglomeration area is mainly distributed in the western part of south XJ; however, the spatial changes in these two agglomerations were relatively small, while the changes in the LH agglomeration area were larger.During the period of 2000-2020, the state of water resource pressure agglomeration in XJ weakened, the fluctuation of the discrete state was enhanced, the number of spatially correlated prefectures decreased, the mismatch continued to increase, and the differences were gradually highlighted.Therefore, the XJ government should pay more attention to the variability in the water stress between regions to avoid increasing polarization.(3) In terms of its contribution to driving the change in the agricultural water footprint in XJ, the total effect of each factor was always positive during the period of 2005-2020, but the effect value generally showed a downward trend.Among them, 2015 was the turning point for the change in the trend in the total effect value, in which the technological effect (water consumption of CNY 10,000 GDP), the economic effect (GDP per capita), and the environmental effect (total carbon emission of CO 2 ) had the most significant impact on the evolution of the agricultural water footprint in XJ.Therefore, the government should formulate policies to encourage the innovation and development of agriculture-related technologies.In this way, it can promote the suppression of technological effects on the agricultural water footprint, thereby reducing water stress in XJ.The intensity of the effect of different factors on the spatial differentiation of the water stress index in XJ was not the same in different time periods.The spatial differentiation of water resources in XJ has mainly been determined by the expansion of its irrigation scale (cultivated area), the growth of its population (total population), the imbalance of its industrial structure (share of agricultural water use), and the evolution of the area of natural oases.
This study can provide new perspectives for further analyzing the actual state of the water footprint and water stress in XJ and provide a reference basis for decision-makers in the XJ government.Future research in this field should pay more attention to the following aspects: firstly, the spatial matching of the agricultural water footprint and land resources and their sensitivity; secondly, accurate assessment of a suitable regional carrying capacity and development of the assessment tools; thirdly, feasibility studies on reducing water stress through integrated planning of transregional water transfer projects and cultivation scale control.

Figure 1 .
Figure 1.The topographic map of XJ, including north XJ, south XJ, and east XJ.Figure 1.The topographic map of XJ, including north XJ, south XJ, and east XJ.

Figure 1 .
Figure 1.The topographic map of XJ, including north XJ, south XJ, and east XJ.Figure 1.The topographic map of XJ, including north XJ, south XJ, and east XJ.
) were obtained from the XJ Statistical Yearbook (2001-2021) and the XJ Production and Construction Corps Statistical Yearbook (2001-2021) provided by the XJ Bureau of Statistics and XJ Production and Construction Corps Bureau of Statistics.

Figure 2
Figure 2 depicts the trend in XJ's agricultural water footprint (AWF) from 2000 to 2020.The results show that the total AWF increased by 124.0% over the past 20 years (WF crop and WF animal increased by 147.9% and 48.9%, respectively).The regional differences in the annual average value of AWF were obvious; in order from highest to lowest are the southern border (27.74 × 10 9 m 3 ), the northern border (21.67 × 10 9 m 3 ), and the eastern border (2.43 × 10 9 m 3 ).Among them, the AWF at the southern border was 11.42 times higher than that at the eastern border.In the agricultural production process, the WF animal multi-year average was 22.4 per cent of the AWF multi-year average, so it can be seen that WF crop is a major component of AWF (~77.6%).Of the three components that make up WF crop in each region, blue water dominates (74.6%), with relatively small proportions of green water (14.0%) and gray water (11.4%).Similar to WF crop , WF crop−blue showed a gradual increase in all zones from 2000 to 2020, with the southern border contributing the most to the components of WF crop−blue (60.0%), followed by the northern border (34.2%) and the eastern border (5.8%).Compared with WF crop−blue , the interannual changes in WF crop−green and WF crop−grey showed an overall relatively small increase in all regions.

Figure 2 .
Figure 2. Recent trends in agricultural water footprint in each district of the XJ.North XJ (abbreviated as N); south XJ (S); east XJ (E).

Figure 2 .
Figure 2. Recent trends in agricultural water footprint in each district of the XJ.North XJ (abbreviated as N); south XJ (S); east XJ (E).
Water 2024, 16, x FOR PEER REVIEW 9 of 20
Water Resource Pressure 3.3.1.Global Spatial Autocorrelation Analysis of Water Resource Pressure in XJ The global Moran's I index and significance level p-value were calculated based on a neighborhood queen spatial weight matrix for the 2000, 2005, 2010, 2015, and 2020 time cross-section data.This revealed the spatial clustering/dispersion of the water stress indices in XJ's prefectures and cities in the last 20 years.The calculated Moran's I and p-values were 0.1089 (p = 0.0428), −0.2851 (p = 0.0042), −0.2109 (p = 0.0032), −0.3301 (p = 0.0025), and −0.3258 (p = 0.0029).

Figure 4 .
Figure 4.A LISA cluster map of water resource pressure in XJ from 2000 to 2020 (a-e).

Figure 4 .
Figure 4.A LISA cluster map of water resource pressure in XJ from 2000 to 2020 (a-e).

20 Figure 5 .
Figure 5.The contribution value of water footprint's driving factors in XJ.

Figure
Figure The contribution value of water footprint's driving factors in XJ.
α TN , α TP are the fertilizer loss coefficients of nitrogen and phosphate fertilizers, respectively, taken from the First National Pollution Source Census-Manual of Fertilier Loss Coefficients of Agricultural Pollution Sources; L nit− f er , L pro− f er are the discounted pure application amounts (kg) of the nitrogen fertilizer and phosphate fertilizer, respectively;

Table 1 .
The effect intensity of each influencing factor in different years.