Spatial and Temporal Variations in Soil Moisture for a Tamarisk Stand under Groundwater Control in a Hyper-Arid Region

: In hyper-arid regions, soil moisture’s role in ecohydrological processes can differ significantly from that in arid or semi-arid ecosystems. We investigated the spatial–temporal dynamics of soil moisture and its relationship with groundwater depths in a 200 m × 300 m phreatophytic tamarisk stand in the lower basin of the Tarim River, a hyper-arid zone in China. Soil moisture proﬁles, from the surface to the water table, were derived using drilling and oven-drying techniques. Over a three-year period, the soil moisture at multiple depths was continuously monitored in a speciﬁc plot using nine frequency domain reﬂectometry (FDR) sensors. Our results indicate a correlation between horizontal variations in soil moisture and groundwater depths (GWDs). Nevertheless, anomalies in this correlation were observed. Variations in horizontal soil moisture were strongly inﬂuenced by the clay content in the soil, with ﬁner soils retaining more moisture. Despite varying GWDs, soil moisture proﬁles remained consistent, with no distinct correlation between them. Soil moisture exhibited stability across layers, with noticeable changes only adjacent to the water table. These results imply that in hyper-arid environments, soil texture primarily governs soil moisture distribution. However, the limited spatial and temporal scopes in our dataset, constrained by the region’s inhospitable conditions, necessitate further investigation. Future work should prioritize amalgamating diverse data sources to devise a region-speciﬁc soil moisture model for in-depth analysis of hyper-arid regions.


Introduction
Soil moisture is a key determinant of the health of terrestrial ecosystems.It plays a fundamental role in the feedback between the Earth and the atmosphere, as well as in all aspects of vegetation growth and composition [1].At the same time, soil moisture is a key element in the hydrological process at the land surface.It redistributes precipitation by affecting infiltration to runoff and directly influences the evapotranspiration process [2][3][4].The space-time characterization of the soil moisture dynamics may provide significant insight into the understanding of the hydrological, meteorological, and ecological processes acting over a variety of scales.Moreover, knowledge of the temporal and spatial variability of soil moisture is critical for understanding and predicting processes such as the partitioning of received energy at the land surface, vegetation water stress in a heterogeneous landscape, soil carbon, and nitrogen cycles, etc. [5].
The temporal dynamics of soil moisture are typically considered as a stochastic process, which is mainly driven by stochastic precipitation events [6,7].

Lots of investigations into
Water 2023, 15, 3403 2 of 17 temporal variations and the time stability of soil moisture with different ecosystems or under different climate conditions were carried out through ground or remote sensing-based data [8][9][10][11][12][13].Stochastic soil moisture dynamic models were developed for different conditions of soil, topography, and vegetation [14][15][16].The spatial dynamics of soil moisture are mainly influenced by soil texture, vegetation, and topographic patterns [17][18][19][20][21].The impact of spatial heterogeneity of soil properties on the soil moisture budget has been extensively investigated using the scaling theory of soil hydraulic properties [22][23][24][25].These studies show that the equivalent soil parameters, which were defined to capture the effects of spatial heterogeneity on the soil water budget [26], depend not only on the spatial characteristics of soils but also on the climate, particularly the precipitation processes [6,27,28].Most of this research shows that rain plays a driving role in the temporal and spatial variability of soil moisture.
In hyper-arid regions, where annual precipitation is usually less than 60 mm [29], the impacts of rain on soil water dynamics can be neglected.Widely spread sandy dunes appear in these regions due to extreme drought, suggesting that soil water is scarce in these regions.However, there are some inland rivers, such as the Tarim River and the Heihe River in arid northwestern China.These inland rivers foster the desert riparian ecosystems along the rivers, indicating that soil water is enough and maintains the growth and development of the desert vegetation in these riparian zones under hyper-arid climates.Obviously, the soil water dynamics in these desert riparian ecosystems are controlled by groundwater or flooding water rather than rain.However, the impacts of groundwater on soil moisture dynamics in these desert riparian ecosystems under hyper-arid climates have been rarely studied.
As another water source to soil water in a vadose zone, where groundwater depth is shallow, groundwater is particularly important for revealing the spatial and temporal dynamics of soil moisture [30][31][32][33].Some stochastic models have been developed for studying the interactions between rainfall, water table, and vegetation [34][35][36][37][38]. Based on these models, the stationary probability distributions of soil moisture at different depths influenced by capillary flux can be investigated for different soils, climates, and vegetation.However, few field measurements have been carried out to investigate the soil moisture dynamics controlled by groundwater.The field measurement results can improve our understanding of soil moisture dynamics under the control of groundwater, provide valuable information and datasets to validate the different kinds of simulation models of soil moisture and support the ecohydrological study in hyper-arid regions.
In this study, we conducted a field investigation of spatial and temporal variability in soil moisture in a tamarisk (Tamarix spp.) stand in the lower Tarim River basin of China.This area is hyper-arid with a mean annual precipitation of less than 50 mm.The influence of rain on soil moisture can be ignored in this area due to the rare rain.The water source of the studied tamarisk stands is mainly from groundwater rather than flood watering for at least fifty years because this area had been dried out for more than 30 years before artificial water conveyance began in the year 2000 [39].Moreover, the artificial water conveyance project did not cause any flooding to occur over the riverbanks.Tamarisk is a typical xeric phreatophytic shrub, groundwater-dependent, and not sensitive to precipitation [40].The groundwater-controlled and phreatophyte-covered stand provides an ideal field investigation site to explore the soil moisture dynamics under the control of the groundwater and analyze their controlling mechanism in hyper-arid regions.
Our objectives are (1) to reveal the temporal and spatial characteristics of soil moisture distribution in hyper-arid desert riparian vegetation, and (2) to explore the effects of environmental factors, such as climate, water table, and soil properties on spatial and temporal variations of soil moisture under the control of groundwater.Our study will obtain an insight into soil moisture dynamics under the specific environment without the influence of rain, which can deepen our understanding of the interactions between soil, vegetation, and atmosphere in hyper-arid regions, and provide data for ecological research in extreme arid regions.

Site Descriptions and Stand Choice
(1) Site descriptions The Tarim River is located in the Tarim Basin and is the longest inland river in China.Its main catchment has a length of 1321 km and covers an area of 17,600 km 2 .The lower basin of the Tarim River is approximately 250 km long, from the Daxihaizi reservoir to the terminal Taitmar Lake, and lies between the Taklamagan Desert and Kuluk Desert.The terrain is relatively flat.The narrow but long desert riparian forest in the lower basin stretches along the two sides of the Tarim River with a width of about 5 km or less.It is called the "green corridor" and protected the ancient "Silk Road" and the current Nation Road.The vegetation types in the lower basin can mainly be classified into two categories: Populus euphratica woodland and tamarisk shrubland.The Populus euphratica is an arbor species, while the tamarisk is a shrub species.Both are accompanied by small amounts of other shrubs and herbs, such as Lycium ruthenicum, Hal-imodendron halodendron, Phragmites australis, Apocynum venetum, Alhagi sparsifolia, Karelinia caspica, and Glycyrrhiza inflata in the two communities.The coverage of the plant community is low, less than 0.3 in most of this area, and the spatial change is obvious.
One of the distinctive features of the study area is the hyper-arid climate.According to the meteorological records of the Tieganlik Weather Station in the lower Tarim River basin, the multiyear average annual precipitation was approximately 33 mm.However, the annual potential evapotranspiration can reach more than 2000 mm.It belongs to the continental warm temperate hyper-arid climate specifically, one of the driest areas in China and Eurasia.Precipitation generally appears in June, July, and August; in other months, it is scarcely observed.
Due to scarce precipitation, plant growth and soil moisture dynamics are controlled by groundwater in this region.Its basic process of water cycling can be simply summarized as a process of the water movement from the river water to groundwater, and then, to the surface evapotranspiration, in which the water movement in the groundwater-soilplant-air system is one of the critical processes controlling the water cycling in the lower reaches [41].However, the proportion and role of soil moisture in the process of water cycling is unclear.Downstream of the Tarim River had been set off for more than 30 years in the last century on account of high strength utilization of water resources in the upper and middle reach of the Tarim River.The water table fell, and vegetation declined completely.The water table increased and vegetation recovered gradually after uninterrupted water transport since 2000.Soil types can be divided into four categories, sand, sandy loam, loam, and silt loam according to the international textural grade in the lower basin, although the texture of the soil layers is heterogeneous in soil profiles.
(2) Stand choice A tamarisk stand (87 • 54 N, 40 • 27 E), which covered a rectangle area of 200 m × 300 m and had high spatial variability of vegetation coverage (Figure 1), was chosen to investigate the temporal and spatial variations in soil moisture.The vertical distance from the Tarim River channel to the southwest corner of the stand was 175 m, and to the northwest corner of the stand was 300 m.The stand was generally flat, with an altitude of 846 m.The tamarisk community was mainly composed of Tamarix ramosissima, Tamarix hispida, and Tamarix elongat, with small amounts of Halimodendron halodendron (Pall.)Voss, Glycyrrh izainfata, A. venetum L., Kareliniacaspia(Pall.)Less, Alhagi sparsifolia, and Lycium ruthenicum Murr.under the tamarisk canopy.The mean vegetation coverage of the stand was 0.45, with a maximum value of 0.57 and a minimum value of 0.03.The main soil type was silt loam but the soil texture was heterogeneous in different soil layers.The average soil bulk density of the top layer is approximately 1.3 g cm −3 .The stand was a deep groundwater site, with a mean groundwater depth of approximately 5 m in the stand.loam but the soil texture was heterogeneous in different soil layers.The average soil bulk density of the top layer is approximately 1.3 g cm −3 .The stand was a deep groundwater site, with a mean groundwater depth of approximately 5 m in the stand.To reveal the spatial variations in soil moisture within the stand, we divided the stand into 100 sampling plots based on an approximate 20 m × 30 m rectangular plot (Figure 1).Each plot position is represented by rows (Row) and columns (Column), for example, plot (R2, C3) refers to the plot in the 2nd row and 3rd column.In each plot, avoiding plants, we selected one soil sampling point on open ground in the middle or near the middle of the plot.In total, 100 soil moisture samples were obtained.Latitude and longitude information for each sampling point was recorded using a handheld GPS.Samples were taken in a sequence, row by row from south (R1) to north (R10), and within each row, the sequence was column by column from west (C1) to east (C10).The sampling period was from October 6 to October 30.
A manual drilling method was used to obtain the soil samples.For each plot, vertical profile sampling began from the 10 cm surface layer, followed by sampling at 50 cm, and subsequently every 50 cm interval, until the depth reached the groundwater level.The criterion for identifying the groundwater level was when free water overflows at the drill bit, and the drill cannot proceed to sample at deeper levels, thereby leading to soil samples being washed out by the saturated water flow.In this study, the depth of the saturated soil layer was regarded as the groundwater depth (GWD), neglecting the height of the capillary fringe.
Soil samples were packed into aluminum containers.The oven-drying method was used to obtain the gravimetric soil water content (g g −1 ) of each soil sample.In addition, the soil particle diameters of all soil samples were obtained by using an MS-2000 laser  To reveal the spatial variations in soil moisture within the stand, we divided the stand into 100 sampling plots based on an approximate 20 m × 30 m rectangular plot (Figure 1).Each plot position is represented by rows (Row) and columns (Column), for example, plot (R2, C3) refers to the plot in the 2nd row and 3rd column.In each plot, avoiding plants, we selected one soil sampling point on open ground in the middle or near the middle of the plot.In total, 100 soil moisture samples were obtained.Latitude and longitude information for each sampling point was recorded using a handheld GPS.Samples were taken in a sequence, row by row from south (R1) to north (R10), and within each row, the sequence was column by column from west (C1) to east (C10).The sampling period was from October 6 to October 30.
A manual drilling method was used to obtain the soil samples.For each plot, vertical profile sampling began from the 10 cm surface layer, followed by sampling at 50 cm, and subsequently every 50 cm interval, until the depth reached the groundwater level.The criterion for identifying the groundwater level was when free water overflows at the drill bit, and the drill cannot proceed to sample at deeper levels, thereby leading to soil samples being washed out by the saturated water flow.In this study, the depth of the saturated soil layer was regarded as the groundwater depth (GWD), neglecting the height of the capillary fringe.
Soil samples were packed into aluminum containers.The oven-drying method was used to obtain the gravimetric soil water content (g g −1 ) of each soil sample.In addition, the soil particle diameters of all soil samples were obtained by using an MS-2000 laser particle analyzer (Malvern, UK).The soil textures of all soil samples were determined based on the information of the soil particle diameters. (

2) Data processing
To analyze the spatial variations in soil moisture in the horizontal direction within the tamarisk stand, the mean gravimetric soil water content (θ, g g −1 ) for each plot was generated by calculating the arithmetic mean of soil moisture content at various depths on the vertical profile (Equation ( 1)).
where θ i was the gravimetric soil water content at various depths on the vertical profile and n was the sample number of soil moisture on the vertical profile.
The basic statistics of the mean gravimetric soil water content (g g −1 ), corresponding groundwater depths, and soil physical clay content of all sampling points were shown in Table 1.The spatial patterns of soil moisture and corresponding groundwater depths in the stand were generated by ordinary kriging interpolation with a spherical model based on 100 θ data using the software ArcGIS 10.6.In order to explore the variations in soil moisture in the vertical profile within the tamarisk stand, 100 sampling points were classified into five groups, according to different groundwater depths: 3.4-4 m, 4-4.5 m, 4.5-5 m, 5-5.5 m, and 5.5-6.1 m.For each groundwater depth level, the vertical profile characteristics of soil moisture were expressed through arithmetically averaging the gravimetric soil water contents of different sampling points with identical sampling depths.
where SD represented the standard deviation of the gravimetric soil water contents of sampling points, and Mean was the mean value of the gravimetric soil water contents of sampling points.

Analysis of Temporal Variations of Soil Moisture
To reveal the temporal dynamics of soil moisture within the tamarisk stand, a series of instruments were installed at the plot (R6, C5) in the stand (Figure 1), where the tamarisk grew well, to continuously monitor multiple parameters, such as the soil volumetric water content (SWC, m 3 m −3 ), groundwater levels, precipitation, and evapotranspiration.The soil volumetric water contents were measured by nine sensors (Model: FDS100, Unism, Beijing, China) based on the frequency domain reflectometry (FDR), which were buried at various depths of 0.3, 0.6, 1.0, 1.5, 2.0, 2.5, 3.0, 4.0, and 5.0 m.To ensure the measurement accuracy, all nine FDR sensors were calibrated on-site using the oven-drying method.To continuously monitor variations in the groundwater table, an observation well was drilled at this location with a depth of 20 m.An automatic groundwater level logger (CTD-Diver, Eijkelcamp, Netherlands) was used to record the groundwater depth.An automatic rain gauge (RM Young Inc., Washington, DC, United States) was used to record precipitations.Surface evapotranspiration was obtained from the flux datasets.Detailed data processing was introduced in the study of Yuan et al. [44].The observation data above Water 2023, 15, 3403 6 of 17 had a synchronous sampling frequency of half an hour.The basic statistics of soil moisture and groundwater depth from the continuous monitoring site were shown in Table 2.The horizontal spatial pattern of soil moisture in the stand is shown in Figure 2a.The gravimetric soil water content (θ) varied from 0.05 to 0.25 g g −1 , with a coefficient of variation (CV) of 33%, thereby suggesting moderate variability.There was a general decrease in soil moisture as the distance from the river channel increased (Figure 2a).The primary determinant was that the tamarisk stand in our study was under the control of the groundwater.As one moves further from the river channel, there is a concomitant decline in groundwater levels (Figure 2b) and, consequently, a decrease in soil moisture content in the direction perpendicular to the river channels (Figure 2a).However, abnormal relationships between soil moisture and GWD were also detected.In the plots around plot (R4, C1) and plot (R2, C2), which were close to the river channel, the GWD was shallower (Figure 2b) but the gravimetric soil water content (θ) was lower (Figure 2b).Additionally, in the plots around plot (R6, C6) and plot (R1, C6), which were relatively far from the river channel, the GWD was deeper (Figure 2b) but the θ was higher (Figure 2a).These abnormal relationships indicated that there were other factors influencing horizontal spatial variations of soil moisture, and their impacts might even exceed that of the GWD.The spatial distribution of the soil moisture was similar to that of the soil physical clay content (Figure 2a,c), suggesting that soil texture could determine the soil water content in this area.We analyzed the relationship between the physical clay content and moisture in the soil.The result showed that there was a significant positive correlation between soil physical clay contents and the θ values (R 2 = 0.636, p < 0.01).The more soil physical clay content it had, the higher the soil moisture was (Figure 3).
The spatial distribution of the soil moisture was similar to that of the soil physical clay content (Figure 2a,c), suggesting that soil texture could determine the soil water content in this area.We analyzed the relationship between the physical clay content and moisture in the soil.The result showed that there was a significant positive correlation between soil physical clay contents and the  values (R 2 = 0.636, p < 0.01).The more soil physical clay content it had, the higher the soil moisture was (Figure 3).

The Variations of Soil Moisture in Vertical Profiles under Different Groundwater Depths
The GWD ranged from 3.4 m to 6.1 m in the tamarisk stand and generally increased with the increasing distance to the river channel located at the west side of the stand (Figure 2b).Five levels of GWDs were classified according to the range of GWDs in the stands: 3.4-4.0m, 4.0-4.5 m, 4.5-5.0m, 5.0-5.5 m, and 5.5-6.1 m.The general soil moisture profile characteristics were similar between different GWD levels (Figure 4).The mean soil water content (SWC) increased from the surface to the water table.However, there was a soil layer at an approximate depth of 1.0-2.0m, in which the SWC was relatively higher than that in the adjacent deeper layer.Based on the observed soil moisture profiles, four soil layers that expressed the soil moisture profile distribution characteristics could be identified in the stand: the surface layer, the shallow layer, the middle layer, and the deep layer.The surface layer is located at a depth of 0-0.3 m, and the depths of the other three layers depend on the GWD levels.The shallow layer was approximately at a depth of 0.5-2.0m, the middle layer was approximately at a depth of 1.5-2.5 m, and the deep layer was approximately below 2.0 m.

The Variations of Soil Moisture in Vertical Profiles under Different Groundwater Depths
The GWD ranged from 3.4 m to 6.1 m in the tamarisk stand and generally increased with the increasing distance to the river channel located at the west side of the stand (Figure 2b).Five levels of GWDs were classified according to the range of GWDs in the stands: 3.4-4.0m, 4.0-4.5 m, 4.5-5.0m, 5.0-5.5 m, and 5.5-6.1 m.The general soil moisture profile characteristics were similar between different GWD levels (Figure 4).The mean soil water content (SWC) increased from the surface to the water table.However, there was a soil layer at an approximate depth of 1.0-2.0m, in which the SWC was relatively higher than that in the adjacent deeper layer.Based on the observed soil moisture profiles, four soil layers that expressed the soil moisture profile distribution characteristics could be identified in the stand: the surface layer, the shallow layer, the middle layer, and the deep layer.The surface layer is located at a depth of 0-0.3 m, and the depths of the other three layers depend on the GWD levels.The shallow layer was approximately at a depth of 0.5-2.0m, the middle layer was approximately at a depth of 1.5-2.5 m, and the deep layer was approximately below 2.0 m.In the surface layer, the SWC was extremely low, only averaging 0.7% g g −1 , and exhibited little spatial variability.In the other three layers, the mean SWC was 13.3%, 11.5%, and 22.5% g g −1 for the shallow layer, the middle layer, and the deep layer, respectively.According to the SWC values, the surface, shallow, middle, and deep layers could be called the extremely dry layer, relatively wet layer, relatively dry layer, and wet layer, respectively.This soil moisture profile distribution was consistent between different GWD In the surface layer, the SWC was extremely low, only averaging 0.7% g g −1 , and exhibited little spatial variability.In the other three layers, the mean SWC was 13.3%, 11.5%, Water 2023, 15, 3403 9 of 17 and 22.5% g g −1 for the shallow layer, the middle layer, and the deep layer, respectively.According to the SWC values, the surface, shallow, middle, and deep layers could be called the extremely dry layer, relatively wet layer, relatively dry layer, and wet layer, respectively.This soil moisture profile distribution was consistent between different GWD levels (Figure 4), suggesting that the groundwater had no effect on the soil moisture profile shapes when the GWD was deeper than 3.0 m.However, the depths of the relative wet layers and relative dry layers were different between different GWD levels (Figure 4), indicating the effect of groundwater on the soil moisture in the profile.On the other hand, apart from the surface layer, the SWC at the same depth in the stand ranged largely with the volatility of 20 ± 5% g g −1 (Figure 4), indicating the high spatial variability of soil moisture in the stand.
The soil moisture profile could be simply divided into three layers, the surface layer (depth of 0-0.2 m), the middle unsaturated layer, and the saturated layer.The variabilities in soil moisture content were various with different layers (Table 3).The soil moisture content exhibited strong variations in the surface layer, which was much greater than those both in the middle unsaturated layer (moderate variations) and in the saturated layer (weak variations).The depth of the soil profile with extreme values of CV did not display special patterns in the middle unsaturated layer.In addition, the depth with the minimal value of CV was much deeper than that with the maximum value of CV in the middle unsaturated layer.The SWC in the profile was affected by the groundwater table overall.The soil volumetric water content (SWC) hourly values of three consecutive days in the three phenological stages (greening period (GP), maturity period (MP), and senescence period (SP) and non-growing season of tamarisk were selected to show the diurnal variations in soil moisture and compared with corresponding water table variations (Figure 5).Because the temporal variations in the soil moisture in most soil layers were similar, we only selected the SWC values at the depths of 0.3, 4.0, and 5.0 m, to exhibit the temporal variations.The SWC values at the depths of 0.3, 4.0, and 5.0 m represented the variations in the surface layer, the middle layer, and the layer near the water table, respectively.
In the non-growing season, the soil moisture in all the soil layers did not exhibit any diurnal variations (Figure 5a,e,i).The soil moisture at the depths of 0.3 m and 4.0 m did not exhibit any diurnal variation in the whole growing season (Figure 5b-d,f-h).However, the soil moisture at the depth of 5.0 m showed obvious diurnal fluctuation in MP, in which the SWC at the depth of 5.0 m decreased in the daytime and got recharged at night (Figure 5k).This diurnal fluctuation occurred in MP but did not in other phenological stages in the growing season (Figure 5j,l).The GWD showed diurnal fluctuation in the whole growing season (Figure 5m-p).The GWD diurnal fluctuations were weak in GP (Figure 5n) and SP (Figure 5p) but obvious in MP (Figure 5o). the three phenological stages (greening period (GP), maturity period (MP), and sene cence period (SP) and non-growing season of tamarisk were selected to show the diurn variations in soil moisture and compared with corresponding water table variations (F ure 5).Because the temporal variations in the soil moisture in most soil layers were simil we only selected the SWC values at the depths of 0.3, 4.0, and 5.0 m, to exhibit the tempo variations.The SWC values at the depths of 0.3, 4.0, and 5.0 m represented the variatio in the surface layer, the middle layer, and the layer near the water table, respectively.In the non-growing season, the soil moisture in all the soil layers did not exhibit a diurnal variations (Figure 5a,e,i).The soil moisture at the depths of 0.3 m and 4.0 m d In the study stand, surface evapotranspiration (ET) occurred in growing seasons but was very weak in the non-growing seasons [41,44].The ET processes and the relevant root uptake would determine the daily processes of the soil water and groundwater in the absence of the influence of precipitation under a hyper-arid climate.The diurnal variations in the soil water in different layers and the water table implied that the groundwater was the only water source for plants and that the plant root uptake should mainly occur near the water table.

Seasonal Variations
Three years of data on the soil moisture at different depths (0.3 m, 4.0 m, and 5.0 m), the corresponding GWD, and daily ET data are shown in Figure 6.The soil moisture in most of the soil layers hardly changed in the first two years, except in the layer near the water table, where the SWC in the non-growing season remained at 51%, although after the onset of the growing season, it would fall gradually until approximately 39%, and then remain stable (Figure 6a).In the late growing season, the SWC would rise gradually to 51%.This course was typical in the first year of the observation.However, in the second year, the date of the onset of the SWC decrease occurred later and finished earlier than in the first year (Figure 6a).
In the three years, the water table began to fall with the emergence of ET, when the plant growing season began, yet turned to rise in the middle of August (Figure 6b,c).The ET processes and the sudden increase in the lateral recharge to the groundwater from the Tarim River because of the water conveyance in August controlled the seasonal dynamics of the groundwater table in the stand [44].
water table, where the SWC in the non-growing season remained at 51%, although after the onset of the growing season, it would fall gradually until approximately 39%, and then remain stable (Figure 6a).In the late growing season, the SWC would rise gradually to 51%.This course was typical in the first year of the observation.However, in the second year, the date of the onset of the SWC decrease occurred later and finished earlier than in the first year (Figure 6a).The GWD in the second year of observation was shallower than in the first year of observation (Figure 6b).The capillary and water potential gradient in the shallower groundwater in the second year of observation could sustain longer-term saturated soil water at the depth of 5.0 m.Therefore, a shorter period of SWC decrease in the second year of observation could be expected.The seasonal course suggested that the SWC near the water table was controlled by the plant water uptake and the groundwater table.
It should be worth noting that some of the SWC observations in the third year were missing due to the damaged instrument.The seasonal variations in SWC were only analyzed for the first two years.

Impacts of Precipitation and Soil Evaporation
Precipitation is an important process controlling the soil water dynamics.However, precipitation is very scarce in our study area, thus, its impact on soil water dynamics can usually be ignored.We analyzed the variations in the SWC at the depth of 0.3 m during the precipitation periods.On 7 June in the second year of SWC observations, heavy rain of 22.3 mm occurred during the period from 1:30 am to 10:30 am, which was the largest precipitation process in the three years.However, the SWC at the depth of 0.3 m remained stable in the following two days (Figure 7) and did not respond to the heavy precipitation event.This indicated that the rainfall infiltration did not reach the layer of 0.3 m in depth.The strong evaporative demand in the study area could cause the rainwater in surface soil to evaporate quickly before the rainwater could reach the soil layers deeper than 0.3 m.During these heavy precipitation events, the stable SWC in the layer at the depth of 0.3 m indicated that most of the precipitation events could not influence the soil water dynamics in the whole soil profile, except the surface soil.
stable in the following two days (Figure 7) and did not respond to the heavy precipitation event.This indicated that the rainfall infiltration did not reach the layer of 0.3 m in depth.The strong evaporative demand in the study area could cause the rainwater in surface soil to evaporate quickly before the rainwater could reach the soil layers deeper than 0.3 m.During these heavy precipitation events, the stable SWC in the layer at the depth of 0.3 m indicated that most of the precipitation events could not influence the soil water dynamics in the whole soil profile, except the surface soil.Soil evaporation is an important factor influencing the dynamics of the soil moisture in the surface and shallow layers.Our water and heat flux observation in the stand indicated that the soil evaporation was weak in the hyper-arid regions [41,44].The surface ET in the non-growing seasons was less than 0.1 mm d −1 , especially in April before the growing seasons (Figure 8).In April, the air temperature and wind speed were close to them in the growing season and the potential evapotranspiration (ET 0 ) was almost the same as that in the growing season, although ET was still less than 0.1 mm d −1 .Until the onset of the growing season, the surface ET started to increase (Figure 8).This phenomenon indicated that the surface ET was mainly from plant transpiration, and that the soil evaporation was very weak.The occurrence of the extremely dry surface layer in this area can prevent soil evaporation, thereby reducing the loss of soil moisture.Soil evaporation is an important factor influencing the dynamics of the soil moisture in the surface and shallow layers.Our water and heat flux observation in the stand indicated that the soil evaporation was weak in the hyper-arid regions [41,44].The surface ET in the non-growing seasons was less than 0.1 mm d −1 , especially in April before the growing seasons (Figure 8).In April, the air temperature and wind speed were close to them in the growing season and the potential evapotranspiration (ET0) was almost the same as that in the growing season, although ET was still less than 0.1 mm d −1 .Until the onset of the growing season, the surface ET started to increase (Figure 8).This phenomenon indicated that the surface ET was mainly from plant transpiration, and that the soil evaporation was very weak.The occurrence of the extremely dry surface layer in this area can prevent soil evaporation, thereby reducing the loss of soil moisture.Although the climate is extremely arid in the study area, the soil layers below a depth of 0.5 m still maintain abundant soil water (Figure 4).The occurrence of high SWC in the  Although the climate is extremely arid in the study area, the soil layers below a depth of 0.5 m still maintain abundant soil water (Figure 4).The occurrence of high SWC in the soil layer under the hyper-arid climate was mainly due to the weak soil evaporation, while the water consumption by the tamarisk in this area was mainly from groundwater [44].The high SWC in the soil layer probably exists widely in hyper-arid regions because, in those areas, scarce precipitation and high evaporation demand can create an extremely dry land surface, which hinders the vapor diffusion from the soil layers below the land surface; therefore, retaining much of the water storage in the soil.
Owing to the little precipitation, the groundwater was the only water source of soil most of the time.Therefore, the water table fluctuation will be the main factor in controlling the soil water dynamics in the study area.However, the slow drop or rise in the water table did not influence the soil water dynamic changes in the soil profile (Figure 6), suggesting that only rapid and large fluctuations in the water table can probably be able to clearly influence the soil water dynamics.Because the groundwater table seldom fluctuated rapidly and largely in the natural environment, our observation results indicated that the soil water temporal variation, only controlled by the groundwater, would be insignificant in most circumstances.
Our observation results showed temporal variations in soil water contents, which were mainly found near the water table, suggesting that the water exchange below ground predominantly occurred between the groundwater and the soil water near the water table.These results indicated that the effects of the tamarisk on the soil water dynamics in the study area were mainly exhibited in the deep soil layer near the water table.Most of the vegetation in the hyper-arid regions is usually the phreatophyte species, which use the soil water in deep layers or groundwater.We believe that the soil water dynamics in the shallow layers will hardly be influenced by vegetation.

Soil Moisture Profile under Control of Groundwater and Its Ecological Implications
There was no infiltration process in the soil profile in our study area, which is an important process in determining the soil water dynamics in most circumstances [45].The upward flow is most often the only form of water movement in the soil layer in this area, unless extremely heavy rain events occur, which can occasionally appear at a hundred years in scale.In this upward flow in the groundwater-soil-plant-air (GSPA) system, groundwater was the input water, and surface ET was the output water.
Our three-year observation results of the temporal processes of soil water contents in different soil layers showed that the soil moisture nearly did not change with time, except for in the layer near the water table.Another study on groundwater ET in this stand indicated that the seasonal courses of groundwater ET and surface ET were consistent and their amounts were nearly equal [46].These results indicated that the upward flow of water in the GSPA system was in either a steady or quasi-steady state.In this steady or quasi-steady state, the groundwater evapotranspiration rate was approximately equal to the surface evapotranspiration rate, and the soil water content remained stable in this system.The soil layer (including the plant) seemed to act as a water channel to the land surface, and in this channel, the water flow rate did not change.
Although the steady or quasi-steady upward flow occurred in the specific stand that we studied, we argue that the steady or quasi-steady upward flow of soil water in the GSPA system is probably universal in hyper-arid regions.Three reasons can support this argument.First, the groundwater table seldom changes rapidly and stochastically in natural environments in these areas.A stable rather than stochastic water table variation process provides a steady background of water environments in hyper-arid regions.Second, vegetation can adapt to the stable water environment in hyper-arid regions, and consequently, evolve and form a steady water use strategy; thus, will not disturb the steady-state flow.Third, soil evaporation is an unsteady-state process because it is mainly controlled by variable weather, but soil evaporation is weak in hyper-arid regions and hardly influences soil water dynamics.We believe that steady-state soil water dynamics is a common process in hyper-arid regions, and quite different from those in most other areas in the world, where soil moisture dynamics are a stochastic process induced by rain events [47], far away from the steady-state.

Soil Moisture Profile under Hyper-Arid Climate and Its Ecological Implications
The soil moisture profiles under different GWD levels exhibited similar curve shapes in the stand and could be divided into four layers: extremely dry surface layer, relatively wet shallow layer, relatively dry middle layer, and wet deep layer.This curve shape in the soil moisture profile is probably typical for a homogeneous soil profile under a hyper-arid climate when the GWD is deep enough, for example, more than 3.0 m.In this soil moisture profile (Figure 9), the surface layer is extremely dry due to the hyper-arid climate.In the shallow layer, which is located below the surface layer, the soil water content depends on the ability and magnitude of the water supply from the land surface, such as occasional heavy rain or flooding in a riparian zone.The soil moisture in the middle and deep layers is controlled by groundwater and its capillary rise, thereby exhibiting a decreasing trend from the water table to up.
Water 2023, 15, x FOR PEER REVIEW 15 of 18 climate when the GWD is deep enough, for example, more than 3.0 m.In this soil moisture profile (Figure 9), the surface layer is extremely dry due to the hyper-arid climate.In the shallow layer, which is located below the surface layer, the soil water content depends on the ability and magnitude of the water supply from the land surface, such as occasional heavy rain or flooding in a riparian zone.The soil moisture in the middle and deep layers is controlled by groundwater and its capillary rise, thereby exhibiting a decreasing trend from the water table to up.Terrestrial ecosystem in hyper-arid regions is usually a groundwater-dependent ecosystem [48].Groundwater depth influences the structure and functions of desert ecosystems in those areas [49].Our observational results also indicated that the groundwater is the only water source of plant water uptake.However, we think the roles of soil moisture in ecosystem development in hyper-arid regions have not gotten enough recognition.We believe that the soil moisture in the shallow layer is a critical factor in influencing the spatial structure of desert ecosystems because it probably determines the conditions for seed germination in hyper-arid regions.
We suppose that the soil moisture in the shallow layer determines the seed sprout and the GWD determines the growth of the desert vegetation.Since the soil moisture exhibits a positive linear relationship with the physical clay content in the soil, it means that the finer the soil texture is, the higher the desert vegetation productivity probably is in hyper-arid regions.A "positive texture effect" of biogeography in hyper-arid regions probably exists, which is quite different from that in arid and semi-arid regions, where an "inverse texture effect" can usually be observed [50].This hypothesis of the "positive texture effect" needs to be verified through widespread field investigations in hyper-arid re- Terrestrial ecosystem in hyper-arid regions is usually a groundwater-dependent ecosystem [48].Groundwater depth influences the structure and functions of desert ecosystems in those areas [49].Our observational results also indicated that the groundwater is the only water source of plant water uptake.However, we think the roles of soil moisture in ecosystem development in hyper-arid regions have not gotten enough recognition.We believe that the soil moisture in the shallow layer is a critical factor in influencing the spatial structure of desert ecosystems because it probably determines the conditions for seed germination in hyper-arid regions.
We suppose that the soil moisture in the shallow layer determines the seed sprout and the GWD determines the growth of the desert vegetation.Since the soil moisture exhibits a positive linear relationship with the physical clay content in the soil, it means that the finer the soil texture is, the higher the desert vegetation productivity probably is in hyper-arid regions.A "positive texture effect" of biogeography in hyper-arid regions probably exists, which is quite different from that in arid and semi-arid regions, where an "inverse texture effect" can usually be observed [50].This hypothesis of the "positive texture effect" needs to be verified through widespread field investigations in hyper-arid regions.

Conclusions
The spatial and temporal variations in soil moisture in the studied stand exhibited distinct and unique characteristics.Horizontal spatial variations in soil water contents were highly related to the soil texture, and finer soils had higher soil water contents.The soil moisture vertical profiles were similar at different groundwater depths, and the relationship between the soil water content and groundwater depth was not obvious.Except in the layer nearest the water table, the soil moisture in the other layers varied minimally over time.Nominal soil evaporation and special water use strategies by phreatophytic species in the hyper-arid climate caused unique, temporal soil moisture processes.We suggest that water flows via groundwater-soil-plant-air systems in hyper-arid regions that are in steady or quasi-steady states.As a result, soil moisture generally remains constant over time, and the amount of water in the soil is mainly determined by the soil texture in hyper-arid regions.
In the presented study, several limitations warrant acknowledgment: (1) The challenging natural environment in the study area posed difficulties for data collection, leading to a dataset characterized by a brief time series and a confined spatial range.This restricts the generalizability of our conclusions.For more comprehensive research in the future, we advocate for the amalgamation of diverse data sources-including meteorological, topographical, and soil property data-to develop a model elucidating the spatiotemporal distribution of soil moisture in hyper-arid regions.
(2) This investigation did not account for the groundwater level fluctuations in response to incoming river water-a factor that could influence the study's findings.The groundwater in the Tarim River's lower reaches primarily originates from incoming river water replenishments.However, the dynamics and magnitude of these replenishments are not yet well-defined and will be the focus of subsequent research endeavors.

Figure 1 .
Figure 1.The location of the study area.

Figure 1 .
Figure 1.The location of the study area.

Figure 3 .
Figure 3. Linear relationship between the gravimetric soil water content (SWC) and soil physical clay content.

Figure 3 .
Figure 3. Linear relationship between the gravimetric soil water content (SWC) and soil physical clay content.

Figure 5 .
Figure 5.Diurnal variations of soil moisture (SWC) in different layers.(a-d) the depth of 0.3 m; h) the depth of 4.0 m; (i-l) the depth of 5.0 m; (m-p) the corresponding water tables in differe phenological stages of tamarisk.

Figure 5 .
Figure 5.Diurnal variations of soil moisture (SWC) in different layers.(a-d) the depth of 0.3 m; (e-h) the depth of 4.0 m; (i-l) the depth of 5.0 m; (m-p) the corresponding water tables in different phenological stages of tamarisk.

Figure 6 .
Figure 6.Seasonal variations in soil moisture (SWC) in different layers and the corresponding water table (GWD) and surface evapotranspiration (ET) values for the three years of observations.The gray regions represent the growing seasons.(a) SWC.(b) GWD.(c) ET.

Figure 6 .
Figure 6.Seasonal variations in soil moisture (SWC) in different layers and the corresponding water table (GWD) and surface evapotranspiration (ET) values for the three years of observations.The gray regions represent the growing seasons.(a) SWC.(b) GWD.(c) ET.

Figure 7 .
Figure 7. Temporal variation in soil water content (SWC) at a depth of 0.3 m after a heavy precipitation event.
Water 2023, 15, x FOR PEER REVIEW 13 of 18

Figure 7 .
Figure 7. Temporal variation in soil water content (SWC) at a depth of 0.3 m after a heavy precipitation event.

Figure 8 .
Figure 8. Variations in surface evapotranspiration (ET) and reference evapotranspiration (ET0) before and after the beginning of the tamarisk growing season.
Quasi-Steady-State Upward Flow of Soil Water under Hyper-Arid Climate

Figure 8 .
Figure 8. Variations in surface evapotranspiration (ET) and reference evapotranspiration (ET 0 ) before and after the beginning of the tamarisk growing season.

1 .
Quasi-Steady-State Upward Flow of Soil Water under Hyper-Arid Climate

Figure 9 .
Figure 9. Schematic drawing of the general soil moisture profile in homogeneous soil from a hyperarid climate.

Figure 9 .
Figure 9. Schematic drawing of the general soil moisture profile in homogeneous soil from a hyperarid climate.

Table 1 .
The basic statistics of soil moisture (g g −1 ), corresponding groundwater depths, and soil physical clay content of all sampling points.

Table 2 .
The basic statistics of soil moisture (m 3 m −3 ) and groundwater depth (m) from the continuous monitoring site.

Table 3 .
Vertical variabilities in soil moisture content at different groundwater depths.GWD denoted the groundwater depth.N denoted the number of samples.