Enhancing Water and Soil Resources Utilization via Wolfberry–Alfalfa Intercropping

The impact of the intercropping system on the soil–plant–atmosphere continuum (SPAC), encompassing soil evaporation, soil moisture dynamics, and crop transpiration, remains an area of uncertainty. Field experiments were conducted for two years in conjunction with the SIMDualKc (Simulation Dual Crop Coefficient) model to simulate two planting configurations: sole-cropped wolfberry (Lycium barbarum L.) (D) and wolfberry intercropped with alfalfa (Medicago sativa L.) (J). These configurations were subjected to different irrigation levels: full irrigation (W1, 75–85% θfc), mild deficit irrigation (W2, 65–75% θfc), moderate deficit irrigation (W3, 55–65% θfc), and severe deficit irrigation (W4, 45–55% θfc). The findings revealed that the JW1 treatment reduced the annual average soil evaporation by 32% compared with that of DW1. Additionally, mild, moderate, and severe deficit irrigation reduced soil evaporation by 17, 24, and 36%, respectively, compared with full irrigation. The intercropping system exhibited a more efficient canopy structure, resulting in reduced soil evaporation and alleviation of water stress to a certain extent. In terms of temporal dynamics, monocropping resulted in soil moisture levels from 1% to 15% higher than intercropping, with the most significant differences manifesting in the mid to late stages, whereas differences in the early stages were not statistically significant. Spatially, the intercropping system exhibited 7–19% lower soil water contents (SWCs) than sole cropping, primarily within the root water uptake zone within the 0–60 cm soil layer. The intercropping system showed an enhanced water absorption capacity for plant transpiration, resulting in a 29% increase in transpiration compared with sole cropping, thereby achieving water-saving benefits. These findings contribute to our understanding of the agronomic and environmental implications of intercropping wolfberry and alfalfa in arid regions and provide insights into optimizing water and soil resource management for sustainable agricultural practices.


Introduction
Wolfberry (Lycium barbarum L.) is a perennial deciduous shrub belonging to the Solanaceae family.It has a rich cultivation history spanning over a millennium in China.Renowned for its robust traits, such as salt and alkali resistance, adaptability to infertile soils, drought tolerance, and high economic returns [1], wolfberry is a significant agricultural commodity.The arid and semi-arid regions of Northwest China have emerged as primary production hubs for this remarkable crop owing to favorable climatic conditions characterized by ample sunlight and notable day-night temperature fluctuations [2].As the market demand for wolfberry continues to surge and local efforts to harness saline-alkali lands intensify, the cultivation footprint of wolfberry has consistently expanded.However, Plants 2024, 13, 2374 2 of 24 the current state of wolfberry cultivation in China remains in a developmental phase, marked by monocropping, meager vegetation cover, severe soil evaporation, and extensive field management practices.In tandem, these factors exacerbate secondary soil salinization.Intercropping systems featuring the harmonious coexistence of woody plants and grasses have emerged as a flagship approach to eco-agriculture [3].They stand out as a means to effectively curtail soil evaporation, elevate water utilization efficiency, and bolster economic return [4].Alfalfa (Medicago sativa L.), celebrated as the "forage king" and "feed queen", emerged as a natural choice for intercropping with wolfberry.Its robust regenerative abilities, palatability, and rich reservoir of high-quality crude protein render it an essential part of wolfberry-based agroforestry systems [5].Thus, investigation of the water-saving mechanisms within the wolfberry-alfalfa intercropping system is pivotal for the sustainable evolution of agroforestry practices.
Agroforestry systems that combine woody and herbaceous plants are widely implemented across various regions, including Europe, South America, the Saharan region of Africa, and South Asia.These systems exhibit diverse forms owing to variations in climatic conditions and management practices [6][7][8].Within agroforestry systems, the interplanting of herbaceous vegetation between woody plants in alleys optimizes land and solar resource utilization.This multifaceted approach not only enhances land coverage but also mitigates soil evaporation, fostering improved vegetation structure and soil enrichment strategies in line with the adage, "grasses to ameliorate the soil, and grasses to stimulate the trees" [8][9][10].The interaction between these two components in agroforestry systems introduces variations in soil moisture, soil evaporation, and crop transpiration.These changes are influenced by factors such as the choice of intercropped vegetation, prevailing climatic conditions, soil texture, and field management practices.Several studies suggest that, although herbaceous plants in agroforestry systems consume some soil moisture, they also contribute significantly to increased rainfall infiltration, reduced soil evaporation, and minimized soil erosion, ultimately resulting in stable or even augmented soil moisture levels [11,12].For instance, in the case of guinea grass (Panicum maximum Jacq.) intercropped with four tree species at varying row spacings, medium row spacing consistently demonstrated the highest soil moisture levels, accompanied by optimal growth and yields [13].Wen et al. [14] observed that intercropping of Alpinia oxyphylla Miq.(a medicinal herb) in oak plantations exhibited variable effects on water consumption rates and uptake across different soil layers, ultimately achieving an equilibrium in soil moisture distribution and overall soil moisture enhancement.In a semi-arid region of Liaoning, China, Bai et al. [15] reported that apricot (Prunus armeniaca L.) orchard intercropping systems positively influenced soil moisture conditions compared with sole cropping.Nonetheless, Oswald et al. [16] found that intercropping grass-legume mixtures with young peach (Prunus persica (L.) Batsch) trees in southern France led to more severe water stress than sole cropping, particularly during early growth when soil cover fabric was used.However, judicious intercrop cutting effectively mitigates water stress during mid-to late-growth stages [17].Similarly, Kurt et al. [18] documented that in a switchgrass (Panicum virgatum L.) and juvenile loblolly pine (Pinus taeda L.) intercropping system, competition exerted a greater negative impact than positive effects, resulting in reduced soil moisture and decreased yields for both switchgrass and juvenile loblolly pine when compared with sole cropping.In jujube (Ziziphus jujuba Mill.) orchards on the Chinese Loess Plateau, intercrop species and field management practices significantly influence the root distribution of agroforestry systems, thereby affecting the soil moisture conditions [19].These findings underscore the need to consider various factors when evaluating the impact of agroforestry systems on soil moisture dynamics, emphasizing the necessity of context-specific approaches to optimize water resource management within agroforestry practices.
Soil moisture plays a pivotal role in agroecosystems and is primarily influenced by two main consumption pathways: soil evaporation and crop transpiration.The introduction of intercropping into these systems reshapes the spatial distribution of plant canopies, thereby exerting a profound influence on the partitioning of soil evaporation and crop transpiration [20].Notably, intercropping often features taller crops, which, owing to their increased stature, can absorb more solar radiation, consequently leading to higher rates of transpiration when compared with shorter companion crops.For instance, a study conducted in the Hetao Irrigation District of Inner Mongolia revealed that maize (Zea mays L.), a tall crop, exhibited a substantial increase in transpiration, ranging from 29.3% to 42.1%, compared with shorter tomato plants [19].This difference in transpiration is predominantly attributed to belowground resource competition, with the root system of maize demonstrating a water uptake capacity 7.3 times greater than that of tomatoes (Solanum lycopersicum L.) [21].To better understand the nuances of soil moisture dynamics, Wang et al. [22] used an enhanced modeling approach to simulate transpiration allocation within agroforestry systems.Their findings indicated that apple trees (Malus pumila Mill.), with their greater stature, consumed 33.3-48.0%more water than cocksfoot grass (Dactylis glomerata L.).However, it is crucial to emphasize that the effect of intercropping on soil moisture is not universally consistent.For instance, in wheat (Triticum aestivum L.) and maize intercropping, adequate underground interaction between crop roots is conducive to reducing water consumption.Increasing maize planting density is beneficial to reduce the water consumption of the long-term subsurface complete interaction between wheat and maize [23].Within the agroforestry systems (crop + shrub) situated in central and southern Senegal, the introduction of shading mechanisms reduces direct solar exposure to the soil surface.After rainfall, there is an initial increase in the soil moisture content, followed by a period of pronounced soil evaporation, with a subsequent gradual decline in the rate of evaporation [24].In semi-arid tropical grasslands where taller trees predominate, the higher proportion of above-ground woody structures reduces evaporation due to "savings" from reduced solar exposure, significantly surpassing the "losses" from tree transpiration.In contrast, agroforestry systems, such as orchards, which are characterized by densely planted, shorter trees, experience a counterbalancing effect due to canopy interception and increased transpiration.Consequently, these dynamics can offset the favorable microclimatic effects.In pursuit of maximum productivity and economic returns, farmers tend to favor economically profitable fruit trees and similar choices [25].It is important to emphasize that soil surface water evaporation is a pivotal component of soil moisture balance, particularly in areas with sparse vegetation, where young trees may struggle to provide complete ground cover.This results in substantial soil moisture loss through ineffective evaporation mechanisms [26].Although intercropping has the potential to reduce soil evaporation within tree rows, it may also intensify soil transpiration [27].Radersma and Ong [28] conducted research in Kenya and reported a 35% reduction in soil evaporation within agroforestry systems compared with bare soil.However, Wang et al. [29] found that intercropping forage grass with jujube orchards in the hilly Loess Plateau partially offset reduced soil evaporation and increased infiltration, resulting in a certain degree of negative impact on soil moisture conditions.Effective field management practices such as timely grass-cutting have been identified as viable solutions for mitigating this impact.In conclusion, the perspectives of researchers, both nationally and internationally, regarding the influence of intercropping systems on soil water content (SWC), soil evaporation, and crop transpiration remain diverse.These disparities emphasize the necessity of regionspecific analyses that account for factors such as climate conditions, crop types, and field management practices.Importantly, the distinctive wolfberry-alfalfa intercropping system, particularly in the arid regions of Northwestern China, is not well-represented in the literature.Therefore, it is imperative to elucidate the spatial distribution of soil moisture within the wolfberry-alfalfa intercropping system, quantify soil evaporation and crop transpiration, explore their intricate interactions, and uncover the mechanisms governing soil-crop water regulation within the soil-plant-atmosphere continuum (SPAC).This undertaking will establish a robust theoretical foundation underpinning rational irrigation practices and sustainable advancement of the wolfberry-alfalfa intercropping system.
Therefore, the aims of this study were to (1) utilize the SIMDualKc (Simulation Dual Crop Coeffcient) model for simulation and validation to demonstrate its applicability in the study area and determine crop coefficients for each growth stage; (2) define the soil evaporation dynamics during the crop growth stages and quantitatively assess the influence of intercropping on soil evaporation; (3) compare the spatiotemporal variations in soil moisture between intercropping and sole cropping; and (4) determine the mechanisms governing soil-plant water dynamics and quantitatively assess the water-saving benefits of intercropping by evaluating crop transpiration across growth stages and their proportion to evaporation, in conjunction with the interactive impacts on soil evaporation and moisture.Through a comprehensive analysis of the water-saving mechanisms at play, we hope to provide valuable insights to the scientific community and promote the adoption of eco-friendly and economically viable agroforestry practices.

Calibration and Validation of the SIMDualKc (Simulation Dual Crop Coefficient)
The model parameters initially adopted the recommended values from FAO-56 [30].Adjustments were made to the parameters during the model calibration process, resulting in a close match between the SWCs (soil water contents) and observed values (Table 1).For alfalfa, a crop that can be harvested multiple times during the growing season, these parameters represent the crop coefficients for each subsequent growth period.K cb ini signifies the crop coefficient after the previous harvest, whereas K cb end represents the crop coefficient before the next harvest.Figure 1 illustrates the observed soil water contents (OSWCs) and simulated soil water contents (SSWCs) for both wolfberry monoculture and wolfberry-alfalfa intercropping under full irrigation treatments.The simulated values for all the treatments were evenly distributed around the 1:1 line.In 2019, for both intercropping and monoculture treatments, "b" values were approximately 0.99, with R-squared (R 2 ) values ranging from 0.88 to 0.89.The root mean square error (RMSE) varied from 0.0163 to 0.0267 m 3 m −3 , whereas the index of agreement (d IA ) ranged from 0.92 to 0.93.The Nash-Sutcliffe Efficiency (EF) remained consistently at 0.70.These results indicated a high level of calibration accuracy.When Plants 2024, 13, 2374 5 of 24 validated with 2020 data, "b" values ranged from 0.98 to 1.00, and R 2 values ranged from 0.86 to 0.87.The RMSE fell within the range of 0.0151-0.0193m 3 m −3 , whereas the d IA values ranged from 0.92 to 0.93.In summary, the soil moisture simulations generated by the SIMDualKc model met the required criteria for accuracy, as detailed in the Table 2 of model validation parameters.
under full irrigation treatments.The simulated values for all the treatments wer distributed around the 1:1 line.In 2019, for both intercropping and monocultu ments, "b" values were approximately 0.99, with R-squared (R 2 ) values ranging f to 0.89.The root mean square error (RMSE) varied from 0.0163 to 0.0267 m 3 m −3 , the index of agreement (dIA) ranged from 0.92 to 0.93.The Nash-Sutcliffe Efficie remained consistently at 0.70.These results indicated a high level of calibration a When validated with 2020 data, "b" values ranged from 0.98 to 1.00, and R 2 value from 0.86 to 0.87.The RMSE fell within the range of 0.0151-0.0193m 3 m −3 , wherea values ranged from 0.92 to 0.93.In summary, the soil moisture simulations gene the SIMDualKc model met the required criteria for accuracy, as detailed in the T model validation parameters.

Comparison of Crop Coefficients
In the wolfberry monoculture treatment (Figure 2a,c), plant growth was sluggish during the initial growth stage, characterized by relatively small leaves and weak transpiration, resulting in a base crop coefficient (K cb ini ) of 0.25.As the growth progressed to the rapid growth phase, leaf development accelerated, transpiration intensified, and the crop coefficient exhibited linear growth.By the mid-growth stage, the base crop coefficient (K cb mid ) had reached 0.9.During the early-to mid-growth stage, when the root system remains relatively shallow, soil moisture satisfactorily meets the plant's growth requirements, which leads to a close alignment between K cb and K cb act .As the late growth stage is set in, plant growth decelerates, and leaves turn yellow and drop, resulting in a reduction in the crop coefficient, with a base crop coefficient (K cb end ) of 0.6.In the late growth stage, after both plant and fruit development are complete, water consumption decreases and corresponding reductions in irrigation are made to prevent water wastage.In 2019 (Figure 2a), the late growth stage witnessed lower rainfall, causing water stress.Nevertheless, it has been concluded that plant and fruit development minimally affect the growth and yield of wolfberry.The soil evaporation coefficient (K e ) varied in response to irrigation and rainfall.During the early growth stage, when leaves are not fully developed and ground cover is limited, water consumption is predominantly driven by soil evaporation, resulting in relatively high K e values.Following irrigation and rainfall events, K e rapidly increased owing to strong evaporation, leading to elevated surface soil moisture content.Subsequently, soil evaporation entered an energy-limited phase, marked by evaporation occurring at its maximum rate, which induced significant K e fluctuations.In the mid-growth stage, as the leaves become nearly fully developed and ground cover increases, transpiration emerges as the primary water consumption process, causing K e to decrease and exhibit reduced fluctuations.The study findings indicated that, for the wolfberry monoculture, K e ranged from 0.43 to 1.02 during the initial growth stage and subsequently decreased to a range of from 0.12 to 0.41 from the mid-growth stage onward.
Plants 2024, 13, x FOR PEER REVIEW 7 of 25 that of the monoculture treatment during the growth period.However, owing to the higher ground cover during intercropping, soil evaporation was reduced, resulting in smaller fluctuations.For wolfberry-alfalfa intercropping, Ke ranged from 0.41 to 1.02 during the initial growth stage and decreased to a range of from 0.06 to 0.21 from the midgrowth stage onward.The interannual differences in the crop coefficients between 2019 and 2020 were relatively small.

Comparison of Soil Evaporation
The soil evaporation simulations for both wolfberry monoculture and intercropping under full irrigation showed a consistent distribution around the 1:1 line (Figure 3).In 2019, for both intercropping and monoculture treatments, the Nash-Sutcliffe efficiency coefficient (b) ranged from 0.89 to 0.93, whereas the coefficient of determination (R 2 ) ranged from 0.91 to 0.93.The RMSE varied from 0.51 to 0.54 m 3 m −3 , and the dIA ranged  For the wolfberry and alfalfa intercropping systems (Figure 2b,d), the dynamic interplay and competition between the two components introduced complexity to the variation in crop coefficients.Owing to the growth characteristics of alfalfa, which shorten the initial growth period of the intercropping system, it undergoes three harvests during its entire growth cycle.Consequently, the crop coefficient for the intercropping system exhibited a periodic three-phase variation curve, with each harvest representing the growth cycle.The overall trend of the crop coefficient within each growth cycle was similar to that of monoculture.During the initial growth stage of the first alfalfa harvest in the intercropping system, the base crop coefficient (K cb ini ) was relatively low, averaging 0.25.In the rapid growth phase, the crop coefficient increased linearly, reaching its maximum value, averaging 1.19, when the ground was completely covered by vegetation.After the alfalfa harvest, the base crop coefficient decreased sharply to 0.9.As alfalfa began to regrow, there was a brief decrease in the intercropping crop coefficient.This temporary decrease was due to the initial low height of the alfalfa plants, resulting in a significant difference in height between alfalfa and wolfberry, which negatively affected the intercropping crop coefficient.With the rapid growth of alfalfa, it entered the rapid growth phase.When the ground was completely covered by the second alfalfa harvest, the base crop coefficient reached a maximum value of 1.15.These results indicated that the time spans of the second and third harvest phases were longer than that of the first harvest phase.After the second harvest (early stage of the third harvest), K cb was 0.9; when the ground was fully covered during the third harvest, it reached 1.13.Before harvest (late growth stage), K cb was 1.08.In 2019 (b), the intercropping treatment also experienced mild water stress during the late growth stage, although the stress level was significantly lower than that of the wolfberry monoculture treatment.This reduction in stress can be attributed to intercropping, which increases ground cover and reduces soil evaporation, helping maintain adequate moisture levels in the root zone.Therefore, intercropping could alleviate a certain degree of water stress.The soil evaporation coefficient (K e ) exhibited a similar trend to that of the monoculture treatment during the growth period.However, owing to the higher ground cover during intercropping, soil evaporation was reduced, resulting in smaller fluctuations.For wolfberry-alfalfa intercropping, K e ranged from 0.41 to 1.02 during the initial growth stage and decreased to a range of from 0.06 to 0.21 from the mid-growth stage onward.The interannual differences in the crop coefficients between 2019 and 2020 were relatively small.

Comparison of Soil Evaporation
The soil evaporation simulations for both wolfberry monoculture and intercropping under full irrigation showed a consistent distribution around the 1:1 line (Figure 3).In 2019, for both intercropping and monoculture treatments, the Nash-Sutcliffe efficiency coefficient (b) ranged from 0.89 to 0.93, whereas the coefficient of determination (R 2 ) ranged from 0.91 to 0.93.The RMSE varied from 0.51 to 0.54 m 3 m −3 , and the dIA ranged from 0.95 to 0.96.The EF was in the range of from 0.82 to 0.84.These results highlighted the high calibration accuracy that was achieved.Validation using data from 2020 yielded similar outcomes, with b ranging from 0.83 to 0.89 and R 2 ranging from 0.87 to 0.93.The RMSE varied from 0.50 to 0.55 m 3 m −3 , and the dIA ranged from 0.92 to 0.94 (Table 3).The soil evaporation simulations based on the SIMDualKc model consistently met the required standards.

Soil Evaporation
In both 2019 and 2020, whether in the wolfberry monoculture or intercropped with alfalfa, the soil evaporation trends exhibited consistent patterns.Using 2020 as an example (Figure 4), all treatment groups demonstrated similar behavior in soil evaporation, characterized by higher values during the early growth stage and lower values in the later stages of crop development.Throughout the growing season, soil evaporation followed a pulsating curve, with fluctuations corresponding to rainfall and irrigation events.Among the different irrigation treatments, the curve for the fully irrigated conditions exhibited the smallest fluctuations.As the severity of the water deficit increased, soil evaporation exhibited greater variability.This phenomenon is attributed to increased water stress, which adversely affects the growth and development of crops.Severe water stress inhibits leaf growth, resulting in reduced soil cover, and, consequently, increased soil evaporation.In particular, after irrigation or rainfall events, when shallow soil moisture levels were high, soil evaporation rapidly increased but then gradually decreased as soil moisture levels declined, contributing to more pronounced fluctuations in soil evaporation.Intercropping treatments displayed a more rational spatial distribution of components with higher soil coverage, effectively reducing soil evaporation and mitigating fluctuations during the growing season.In the DW1 treatment (wolfberry monoculture), the average soil evaporation was 292.48 mm, whereas, in the JW1 treatment (wolfberry intercropped with alfalfa), it was 198.68 mm.Intercropping reduced ineffective soil evaporation by 32% compared with monoculture.In the monoculture treatments, soil evaporation decreased with increasing water deficits.Relative to fully irrigated conditions, soil evaporation decreased by 17% under mild deficit irrigation, 24% under moderate deficit irrigation, and 36% under severe deficit irrigation.However, there were no significant differences in soil evaporation among the intercropping treatments under different water deficit levels, indicating that intercropping can alleviate water stress through the interaction of different components.

Impact of Intercropping on Spatiotemporal Soil Moisture Dynamics
During the 2019-2020 growing season, a consistent pattern of soil moisture variation was observed as the cultivation of wolfberry progressed (Figure 5).SWCs initially decreased, followed by a gradual increase.In the early stages, particularly during germination, the SWCs were relatively high owing to ample spring irrigation, which ensured optimal wolfberry growth.However, as wolfberry plants matured, soil moisture levels declined, reaching their lowest levels during the initial fruiting stage.This decline corresponds to a crucial phase of fruit development marked by significant root water uptake, which facilitates fruit maturation.Among the DW treatments (DW1, DW2, DW3, and DW4), the average SWCs during the growth period ranged from 17.74% to 31.74%.For the JW treatments (JW1, JW2, JW3, and JW4), the average SWCs during the growth period varied from 16.25% to 31.44%.Increasing irrigation led to higher soil moisture levels, with DW1, DW2, and DW3 showing 23%, 16%, and 9% higher SWCs than DW4, respectively.Similarly, JW1, JW2, and JW3 exhibited 33%, 20%, and 7% higher SWCs than JW4 under increased irrigation.In general, single cropping (monoculture) resulted in higher soil moisture levels than intercropping.However, under sufficient irrigation conditions, no significant differences were observed between intercropping and single-cropping treatments.This lack of distinction can be attributed to the enhanced water absorption capacity of the intercropped plants through their root systems, rendering the advantage of intercropping less pronounced with ample irrigation.During the initial stages, when spring irrigation was sufficient, both single cropping and intercropping treatments exhibited relatively high soil moisture levels, with no significant differences.At this point, both wolfberry and intercropped alfalfa were in their early growth stages and characterized by relatively low water requirements.Nevertheless, as the growth period progressed, disparities between the single cropping and intercropping treatments gradually became more evident.Research findings indicate that soil moisture levels in single cropping treatments were consistently 1-15% higher than those in intercropping treatments.The temporal and vertical patterns of SWCs remained consistent across different growth stages from 2019 to 2020.Notably, the shallow soil layers exhibited more pronounced fluctuations in the SWCs, with the most significant changes occurring during fruit development.In contrast, the deeper soil layers demonstrated relatively stable moisture levels.Within the surface 0-20 cm soil layer, the SWCs initially ranged from 15% to 20%, primarily because of the combined effects of soil evaporation and water uptake by the crop roots.Early in the season, increased soil surface evaporation led to a decline in the SWCs.However, as the leaves reached full development and ground cover increased, evaporation rates diminished, resulting in soil moisture levels reaching approximately 20-25%.During this period, the depletion of shallow soil moisture was predominantly driven by soil evaporation, and there were no discernible differences in SWCs between the intercropped and monoculture treatments.In the 20-40 cm soil layer, a noteworthy reduction in soil moisture content was observed after the early flowering stage, reaching levels of between 7% and 12%.This phenomenon was attributed to the continuous growth of plant The temporal and vertical patterns of SWCs remained consistent across different growth stages from 2019 to 2020.Notably, the shallow soil layers exhibited more pronounced fluctuations in the SWCs, with the most significant changes occurring during fruit development.In contrast, the deeper soil layers demonstrated relatively stable moisture levels.Within the surface 0-20 cm soil layer, the SWCs initially ranged from 15% to 20%, primarily because of the combined effects of soil evaporation and water uptake by the crop roots.Early in the season, increased soil surface evaporation led to a decline in the SWCs.However, as the leaves reached full development and ground cover increased, evaporation rates diminished, resulting in soil moisture levels reaching approximately 20-25%.During this period, the depletion of shallow soil moisture was predominantly driven by soil evaporation, and there were no discernible differences in SWCs between the intercropped and monoculture treatments.In the 20-40 cm soil layer, a noteworthy reduction in soil moisture content was observed after the early flowering stage, reaching levels of between 7% and 12%.This phenomenon was attributed to the continuous growth of plant roots, which extended the water-absorbing zone deeper into the soil profile.During this stage, when water resources were limited, intercropped Lycium chinense had a more pronounced impact on SWCs due to the regrowth of alfalfa.Consequently, soil moisture levels in the intercropped areas were 12-22% lower than those in the monoculture treatments.After the fruit-setting stage, the SWCs in the 40-60 cm soil layer experienced a slight decrease, with reductions ranging from 2% to 5%.The 0-60 cm soil layer emerged as the primary zone for water uptake by plant roots, supporting fruit development.Conversely, the soil moisture content in the 80-100 cm soil layer remained relatively stable, exhibiting only minor declines during periods of robust plant growth, signifying reduced water extraction by the deeper roots.Intercropping resulted in a more balanced spatial distribution of soil moisture.The introduction of alfalfa between crop rows effectively optimized the utilization of soil moisture in the 0-60 cm soil layer.Notably, the SWCs within this layer in the intercropped plots was consistently lower, ranging from 7% to 19%, compared with the monoculture treatments (Figure 6).

Transpiration
During the budding stage, the plant is in its early growth phase, and water consumption is primarily dominated by soil evaporation with relatively low transpiration rates.In the monoculture treatment, the Transpiration-to-Evapotranspiration ratio (T/ET) stands at only 24.90%, whereas, in the intercropping treatment, it reached 36.87%.As plants enter the vigorous shoot growth phase, leaf development and increased photosynthetic activity led to a rapid rise in transpiration pull.In response, the root system actively absorbs the soil moisture for transpiration.At this stage, the soil moisture content experienced a rapid decline, shifting the water consumption from soil evaporation to plant transpiration.The T/ET values for the monoculture and intercropping treatments were 67.52% and 82.02%, respectively.Transpiration reached its peak and stabilized by the onset of the initial fruiting stage.With fully developed leaves and maximum ground coverage, evaporation rates decrease and plant water uptake is predominantly driven by transpiration.During the initial fruiting and subsequent stages, soil moisture remained consistently low with minimal fluctuations.In the monoculture treatment, the T/ET values were 73.99%, 74.20%, and 71.92% for the initial, peak, and late fruiting stages, respectively.Conversely, the intercropping treatment resulted in T/ET values of 83.58%, 83.22%, and 84.11% for the corresponding stages.Throughout the growth period, the intercropping treatment maintained a T/ET ratio of 76.90%, whereas the monoculture treatment exhibited a lower T/ET ratio of 64.74% (Table 4).The higher T/ET ratio in the intercropping treatment suggests that the spatial arrangement of multiple components was advantageous.It effectively reduced soil evaporation, as discussed earlier, leading to a 32% decrease in ineffective soil evaporation compared with monoculture.Additionally, the intercropping treatment demonstrated a 29% increase in transpiration compared with the monoculture, thereby achieving water-saving effects.Notably, there was no significant difference in the transpiration rates between the monoculture and intercropping treatments in 2019 and 2020.

Crop Coefficients
The variation pattern of crop coefficients during the growth stages is well documented and is characterized by a lower coefficient in the initial growth phase, followed by a linear increase during the rapid growth period.This coefficient reaches its maximum and stabilizes during the mid-growth phase, with a slight decline observed in the lategrowth phase [30].Wolfberry, which are multi-branched shrubs, undergo a relatively short period of rapid growth followed by an extended mid-growth phase [31].In a study by Yin et al. [32] using the single crop coefficient method, wolfberry exhibited a crop coefficient (K c ) of from 0.67 to 0.68 during the spring shoot growth period, from 1.14 to 1.16 during the fruit ripening period, and from 0.89 to 0.92 during the leaf-falling period.The dual crop coefficient method dissects K c into two components: the basal crop coefficient (K cb ) representing crop transpiration and the soil evaporation coefficient (K e ) representing soil evaporation.In this study, using the dual-crop coefficient method, the baseline crop coefficient (K cb ini ) for wolfberry in a monoculture was determined to be 0.25.During the mid-growth phase, K cb mid was 0.9, whereas, in the late growth phase, K cb end was 0.6.The soil evaporation coefficient (K e ) for wolfberry exhibited a range of from 0.43 to 1.02 in the early growth period and from 0.12 to 0.41 in the mid-growth period and beyond.Alfalfa is a perennial forage crop that displays cyclic variations in crop coefficients) [33].In the research of Liu et al. [34], the first and second crops of alfalfa were found to have K c ini values of 0.4, K c mid values of 0.95 and 1.2, and K c end values of 0.9 and 1.15, respectively.These coefficients followed an "S"-shaped trend, increasing initially and then decreasing within each crop cycle.In the wolfberry berry-alfalfa intercropping system, the crop coefficient dynamics align more closely with alfalfa, particularly during the mid-growth phase when alfalfa cutting exhibits periodic changes.In the initial growth phase, K cb was 0.25, whereas in the mid-growth phase, it reached 1.19.As alfalfa is typically harvested approximately one week after flowering, the crop coefficient of the intercropping system for wolfberry, K cb , remained at 0.9.With the second growth of alfalfa, the crop coefficient of the intercropping system exhibited linear growth, reaching 1.15 when alfalfa completely covered the soil surface.After cutting, in the third crop cycle, K cb was 0.9, and it reached 1.13 when the soil surface was entirely covered before harvesting (late growth phase), with a K cb of 1.08.

Soil Evaporation
After rainfall or irrigation, the SWCs initially increased and soil evaporation occurred at its maximum rate.During this phase, soil evaporation is primarily limited by the energy available at the surface.As cumulative evaporation increases, soil moisture decreases, leading to a reduction in soil evaporation rate, a phenomenon known as evaporation decline [35].In the early stages of crop development, when leaves are not fully formed and ground cover is minimal, soil evaporation rates are high.However, as the leaves reach full development, soil evaporation decreases and only increases again following irrigation or rainfall, exhibiting pulse-like fluctuations.These findings are consistent with those of Li [36], who studied evaporation patterns between corn rows in Tongliao.Li [36] found that, early in the growth period, soil evaporation peaked, especially during the seedling stage, accounting for 85% of total evaporation.As corn plants grow, evaporation rates decrease but increase after irrigation or rainfall events.Similar conclusions were drawn by Duan et al. [37] in their study of wolfberry.Using different models, they estimated soil evaporation during the wolfberry growth period to be between 119.56 mm and 139.72 mm.In contrast, our study found an average soil evaporation of 292.48 mm for monocultured wolfberry farming.This discrepancy may be attributed to the shorter planting duration and wider spacing of wolfberry plants in our study area, resulting in higher soil evaporation owing to the arid climate.Furthermore, we observed that, in wolfberry-alfalfa intercropping systems, the average soil evaporation was 198.68 mm, representing a 32% reduction in soil evaporation compared with monoculture wolfberry farming.These findings align with the research of Huo et al. [38] on apple-canola (Brassica napus L.) intercropping systems in the loess hilly region, where intercropping reduced soil evaporation by from 25% to 33%.Intercropping systems featuring tall and short crops can effectively optimize land resource utilization, achieve a more rational spatial layout, reduce soil evaporation, utilize soil moisture between crop rows, and enhance water use efficiency [39].Moderating irrigation can promote water compensation among intercropped components, encouraging crops to tap deeper into soil moisture [40].These results were consistent with our findings, indicating that intercropping can mitigate water stress through interactions between different components.

Soil-Crop Water Regulation
The soil-crop system primarily regulates soil moisture distribution through soil evaporation and crop transpiration [41].Soil evaporation leads to changes in soil moisture, affecting crop water uptake through root systems, which in turn influences crop transpiration.Transpiration is closely linked to crop growth, plant height, and leaf development.These interactions create a complex web of feedback mechanisms within the soil-crop system [42].Compared to monocultures, intercropping systems involve competition among different crop components in both underground and aboveground spaces.This competition inevitably affects the spatiotemporal dynamics of soil moisture.At the initial stages of growth, there was no significant difference in SWCs between the intercropping and monoculture treatments.During this period, water consumption by crops is primarily driven by soil evaporation, with transpiration accounting for a lower proportion of the total evapotranspiration [43].At this stage, when leaves were not fully developed and ground cover was minimal, the soil moisture distribution was similar between the intercropping and monoculture treatments.As crops continue to grow, especially in the mid-to-late stages, leaf development becomes complete and transpiration becomes the dominant mode of water consumption, particularly in intercropping treatments.In these cases, both wolfberry and alfalfa absorb significant amounts of soil moisture for transpiration, resulting in soil moisture levels 1-15% higher in monoculture than in intercropping, especially when irrigation amounts are substantial.Our study revealed that, in the early growth stages, the transpiration/evapotranspiration ratio (T/ET) was only 24.90% for monoculture but reached 36.87% for intercropping.In the mid-to-late stages, T/ET ranged from 64.05% to 76.65% for monoculture and from 76.85% to 85.09% for intercropping.On a spatial scale, the main differences in SWCs between intercropping and monoculture treatments occurred in the 0-60 cm soil layer, which is where active root water absorption occurs [44].There were no significant differences between intercropping and monoculture treatments in the 0-20 cm soil layer.This can be attributed to the fact that the 0-20 cm layer mostly consists of easily evaporable soil [30], in which evaporation is the dominant process and root water uptake plays a secondary role.During the early stages of intercropping, when alfalfa did not provide full ground cover, soil evaporation was high.However, during the mid-to-late stages, alfalfa cover reduced intercropping soil evaporation, whereas crop transpiration increased [45].During the germination phase, intercropping resulted in SWCs in the 20-40 cm layer, which was 12-22% lower than that in the monoculture.This can be attributed to the fact that, during this phase, alfalfa in intercropping is in its regrowth stage, absorbing a significant amount of soil moisture for its growth, causing lower soil moisture levels in intercropping than in monoculture.In the 40-60 cm layer, the soil moisture content decreased after the fruiting stage because of the transition of wolfberry from vegetative to reproductive growth.During this stage, increased transpiration led to a water potential gradient generated by transpiration pull [46], resulting in a from 2% to 5% decrease in SWCs.Quan et al. [47] found that covering apple orchards with sedges significantly converted soil evaporation into plant transpiration.In line with these findings, our study suggests that wolfberry-alfalfa intercropping, with its more rational spatial distribution, maintains soil moisture levels from 7% to 19% lower than monoculture in the primary root water uptake zone (0-60 cm).This allocation of more water for crop transpiration in intercropping leads to a 29% increase in transpiration compared with monoculture, effectively reducing soil evaporation and enhancing water-use efficiency.
Therefore, through a two-year experiment and an in-depth investigation into the soil-crop water regulation mechanisms within the wolfberry and alfalfa intercropping system, it was found that agroforestry practices could alleviate water stress to a certain extent, achieving water-saving goals.However, this study has yet to provide a clear understanding of root distribution within the intercropping system.Overlapping root systems and competition between components during their growth processes add complexity to the soil moisture distribution.To further elucidate the mechanisms of water regulation, future research should focus on clarifying root distribution within intercropping systems.

Site Description
This experiment was conducted at the High-Efficiency Water-Saving Irrigation Experimental Station of the Jingtaichuan Electric Power Pumping Management Bureau in Gansu Province, China (37 • 26 ′ N, 104 • 04 ′ E; average altitude of approximately 1562 m) from 2019 to 2020.The experimental area falls within a temperate arid continental climate zone, which is characterized by low rainfall, with an average annual precipitation of 185 mm.It also experiences high evaporation rates, with an average annual evaporation of 3028 mm.Additionally, the area has significant diurnal temperature variations, with an average annual temperature of 8.3 • C, ranging from a maximum of 35.6 • C to a minimum of −27 • C. The region experiences an annual average of 2652 h of sunlight and an annual average radiation of 6.18 × 10 5 J cm −2 .The frost-free period spans 191 days.The soil type in the experimental area was loam with a bulk density of 1.35 g cm −3 and a field capacity (θ fc ) of 24.1% (mass-based water content).The hydraulic characteristics of the soil are listed in Table 5.The soil nutrient content at the experimental site included organic matter at 1.32 g kg −1 , available nitrogen at 74.51 mg kg −1 , available phosphorus at 26.31 mg kg −1 , and a pH value of 8.11.In 2019, the total precipitation and average temperature during the entire growth period of wolfberry were 227 mm and 19.1 • C, respectively.In 2020, during the entire growth period of wolfberry, the total precipitation and average temperature were 101 mm and 19.7 • C (Figure 7).

Experimental Design
The experiment used a two-factor randomized block design.Factor one consisted of two planting modes: wolfberry monoculture (D) and wolfberry intercropped with alfalfa (J) (Table 6).Factor two involved water management, in which the irrigation plan targeted soil moisture within a designated range (based on the percentage of field capacity with a moistened layer depth of 60 cm): full irrigation (W1, 75-85%θfc), mild deficit irrigation (W2, 65-75%θfc), moderate deficit irrigation (W3, 55-65%θfc), and severe deficit irrigation (W4, 45-55%θfc).In late April 2017, two-year-old wolfberry seedlings were transplanted with north-south row planting, with a spacing of 150 cm between plants and 300 cm between rows.Simultaneously, alfalfa was strip-sown between the rows of wolfberry (in the

Experimental Design
The experiment used a two-factor randomized block design.Factor one consisted of two planting modes: wolfberry monoculture (D) and wolfberry intercropped with alfalfa (J) (Table 6).Factor two involved water management, in which the irrigation plan targeted soil moisture within a designated range (based on the percentage of field capacity with a moistened layer depth of 60 cm): full irrigation (W1, 75-85%θ fc ), mild deficit irrigation (W2, 65-75%θ fc ), moderate deficit irrigation (W3, 55-65%θ fc ), and severe deficit irrigation (W4, 45-55%θ fc ).In late April 2017, two-year-old wolfberry seedlings were transplanted with north-south row planting, with a spacing of 150 cm between plants and 300 cm between rows.Simultaneously, alfalfa was strip-sown between the rows of wolfberry (in the north-south direction) at a distance of 30 cm from the tree trunks, with a row spacing of 30 cm (Figure 8).Eight treatment groups were established with three replicates each, totaling 24 plots.Each plot had an area of 45 m 2 (7.5 m × 6.0 m), with two rows of wolfberry, each containing five plants, making a total of 10 plants per plot.After planting, sufficient water was applied to ensure the survival and growth of both the wolfberry and alfalfa.Water management practices commenced in April 2019, using a drip irrigation system.The field irrigation system included water storage tanks, water pump devices, filtration equipment, main pipes, branch pipes, and micro-tubes (embedded drip irrigation tubes with an inner diameter of 16.0 mm, drip emitter spacing of 0.3 m, and a flow rate of 2.0 L h −1 ).Drip irrigation tapes were placed on both sides of the wolfberry tree trunks at a spacing of 0.3 m.Each experimental plot was equipped with a water meter (precision of 0.0001 m 3 ) to control the irrigation volume, and fertilization was performed using a venturi fertilizer applicator.each containing five plants, making a total of 10 plants per plot.After planting, suffi water was applied to ensure the survival and growth of both the wolfberry and al Water management practices commenced in April 2019, using a drip irrigation sy The field irrigation system included water storage tanks, water pump devices, filtr equipment, main pipes, branch pipes, and micro-tubes (embedded drip irrigation with an inner diameter of 16.0 mm, drip emitter spacing of 0.3 m, and a flow rate of h −1 ).Drip irrigation tapes were placed on both sides of the wolfberry tree trunks at a ing of 0.3 m.Each experimental plot was equipped with a water meter (precision of 0 m 3 ) to control the irrigation volume, and fertilization was performed using a ventur tilizer applicator.

Meteorological Data
Meteorological data, including wind speed, atmospheric pressure, relative hum precipitation, temperature, solar radiation, and other relevant parameters, were auto ically collected at hourly intervals using a compact Davis meteorological station imp from the United States.

Soil Water Content
In each experimental plot, a time-domain reflectometer (TDR) probe extending length of 120 cm was positioned 30 cm from the base of the wolfberry plants.This arr ment was replicated three times.Soil water content (SWC), spanning the 0-120 cm layer, was systematically measured at 7-day intervals.These measurements were ducted using a PICO-BT TDR instrument, which is a German device manufacture

Meteorological Data
Meteorological data, including wind speed, atmospheric pressure, relative humidity, precipitation, temperature, solar radiation, and other relevant parameters, were automatically collected at hourly intervals using a compact Davis meteorological station imported from the United States.

Soil Water Content
In each experimental plot, a time-domain reflectometer (TDR) probe extending to a length of 120 cm was positioned 30 cm from the base of the wolfberry plants.This arrangement was replicated three times.Soil water content (SWC), spanning the 0-120 cm soil layer, was systematically measured at 7-day intervals.These measurements were conducted using a PICO-BT TDR instrument, which is a German device manufactured by IMKO.SWCs were calibrated at various growth stages of wolfberry plants using standard soil drilling and oven drying methods.A proactive irrigation regimen was promptly initiated in response to soil moisture levels falling below the predetermined lower threshold as per the specified moisture treatment requirements.

Soil Evaporation
In the plots designated for ample irrigation experiments, automatic inter-row soil evaporation gauges (manufactured by Beijing Time-Domain Technology Co., Ltd., Beijing, China) were placed at a distance of 300 cm from the base of the wolfberry plants.These devices comprised four key components: a solar-powered unit, weighing assembly, outer protective bucket for the soil column, and soil column itself.The soil columns were 25 cm in height and 20 cm in diameter, with data collection intervals set at 1 h intervals.Near the automatic inter-row soil evaporation gauges in the ample irrigation plots, micro-lysimeters (MLSs) were installed in triplicate to cross-verify the data generated by the automatic inter-row soil evaporation gauges.In addition, in other irrigation experiment plots and ample irrigation plots, micro-lysimeters were placed at the same locations for daily soil evaporation data collection.For the micro-lysimeters, Polyvinyl chloride (PVC) pipes were used as the inner and outer containers, each with a length of 10 cm.The outer cylinder has a diameter of 16 cm, whereas the inner cylinder has a diameter of 14 cm.The bottom was fitted with a breathable 500-mesh screen to prevent soil leakage [35].Weights were recorded daily at 17:00 using an electronic balance with a precision of 0.1 g.To ensure consistency with the in-field soil moisture, the soil was replaced every 5-7 days, and soil replacement was conducted following rainfall or irrigation events.

The SIMDualKc Model
The SIMDualKc (Simulation Dual Crop Coeffcient) model, rooted in the principles of water balance, uses the dual crop coefficient approach to effectively differentiate between two key components within its modeling framework: crop transpiration (K cb •ET 0 ) and soil evaporation (K e •ET 0 ).This model is capable of estimating both soil evaporation and crop transpiration and encompasses various modules, including soil, meteorological, crop, irrigation, and intercropping data [48].In the process of model application, crop transpiration was determined based on calibrated K cb values, with the model output providing the value of K e , signifying daily soil evaporation.
ET a = K c •ET 0 = (K s K cb + K e )ET 0 (1) where ET a is the actual crop evapotranspiration (mm d −1 ), K cb is the basal crop coefficient, K e is the soil evaporation coefficient, K s is the water stress reduction coefficient, and ET 0 is the reference evapotranspiration (mm d −1 ).The model divides the crop growth period into four stages: initial growth, rapid growth, mid-growth, and late growth.The initial values of crop coefficients for these stages (K cb ini , K cb mid , and K cb end ) can be obtained from the FAO-56 recommended table.However, when RH min is not equal to 45% or wind speed (u 2 ) is not equal to 2 m s −1 , the following correction formula is applied to adjust K cb mid and K cb end when they exceed 0.45, as recommended by [30]: The equation uses K cb (input) to denote the basic crop coefficient for K cb mid at a wind speed of u 2 = 2 m s −1 and a minimum relative humidity of RH min = 45%, as referenced in the recommended table.The values of u 2 and RH min correspond to the observed averages, whereas h represents the average crop height.
Soil evaporation predominantly occurs in regions with moist soil surfaces.It was assumed that soil evaporation beneath the plant canopy was within K cb •ET 0 , denoting soil evaporation as K e •ET 0 .The soil evaporation reached its maximum rate when the soil was adequately moist.However, the available energy for surface evaporation constrains the combined crop coefficient to K cb + K e ≤ K c max .As the soil surface dries, soil evaporation becomes directly proportional to the remaining soil moisture content [49].
In the equation, K e represents the soil evaporation coefficient, K cb is the base crop coefficient, K c max denotes the maximum K c value after rainfall or irrigation, K r is a dimensionless coefficient that depends on the cumulative depth of surface soil evaporation (or water consumption) and affects the reduction in evaporation, and few represents the ratio of bare and wet soil, indicating the percentage of the maximum soil evaporation surface area.
The intercropping crop coefficient K c was determined based on the surface area covered by the different intercropping crops and the height of the crop plants.
The equation involves f 1 and f 2 , which represent the surface areas covered by wolfberry and alfalfa, respectively, whereas h 1 and h 2 represent the heights of the wolfberry and alfalfa plants, respectively.
In this study, SIMDualKc requires that the input data consist of the following: (5) Intercropping data included the selected intercropping type, intercropping crop soil coverage ratio, and comprehensive data packages containing information on the intercropping crop, soil, meteorological conditions, and irrigation data.

Statistic Analysis
The SIMDualKc model based on the FAO dual crop coefficient requires input data: meteorological data, soil data, crop data, irrigation data, and intercropping data, and the model parameters are adjusted by giving the initial parameters K cb and soil water deficit coefficient.The regression coefficient (b), coefficient of determination (R 2 ), root mean square error (RMSE), discrete accuracy index of accuracy (d IA ), and model effectiveness coefficient (EF) were used to test the simulation effect of the model.

Conclusions
In conclusion, the SIMDualKc model provided accurate simulations of soil moisture and evaporation for both monoculture and intercropping systems of wolfberry and alfalfa.For monoculture wolfberry, the crop coefficients (K cb ) followed a pattern with values of 0.25 during the early growth stage (K cb ini ), 0.9 during the mid-growth stage (K cb mid ), and 0.6 in the late growth stage (K cb end ).During the rapid growth period, the soil evaporation coefficient (K e ) ranged from 0.85 in the early growth stage to 0.27.In the intercropping system, influenced by alfalfa growth characteristics, the crop coefficients exhibited a cyclic pattern.During the early growth stage, K cb ini was 0.25, which varied during harvests, reaching an average of 1.19 before each harvest, and decreased to 0.9.This pattern continued through the three cycles, resulting in a balanced distribution of soil moisture.The intercropping system effectively reduced soil evaporation by 34%, with soil evaporation decreasing from

Figure 1 .
Figure 1.Comparison of the measured and simulated soil moisture by calibration and val the model for Wolfberry monoculture treatment (a,c) in 2019 and 2020, respectively, inte between wolfberry and alfalfa (b,d) in 2019 and 2020.SSWC, simulated soil water conten observed soil water content.

Figure 1 .
Figure 1.Comparison of the measured and simulated soil moisture by calibration and validation of the model for Wolfberry monoculture treatment (a,c) in 2019 and 2020, respectively, intercropping between wolfberry and alfalfa (b,d) in 2019 and 2020.SSWC, simulated soil water content; OSWC, observed soil water content.

Figure 2 .
Figure 2. Standard and actual basal crop coefficient (Kcb, Kcb act), evaporation coefficient (Ke) for wolfberry monoculture in 2019 (a), 2020 (c) and intercropping between wolfberry and alfalfa in 2019 (b), 2020 (d), with depicting precipitation and irrigation events.The black polyline represents kcb, the black dashed line represents Kcb act, the red dashed line represents Ke, the black vertical line represents precipitation, and the blue vertical line represents irrigation

Figure 2 .
Figure 2. Standard and actual basal crop coefficient (K cb , K cb act ), evaporation coefficient (K e ) for wolfberry monoculture in 2019 (a), 2020 (c) and intercropping between wolfberry and alfalfa in 2019 (b), 2020 (d), with depicting precipitation and irrigation events.The black polyline represents k cb , the black dashed line represents K cb act , the red dashed line represents K e , the black vertical line represents precipitation, and the blue vertical line represents irrigation.

Figure 3 .− 1 )Figure 3 .
Figure 3.Comparison between simulated and observed soil evaporation: on the left, the chang progress of the observed and simulated soil evaporation for wolfberry monoculture (a,e) and in cropping between wolfberry and alfalfa (c,g) in 2019 and 2020; The figure on the right shows regression between the observed and calculated values of soil evaporation from wolfberry mo culture (b,f) and intercropping between wolfberry and alfalfa (d,h) in 2019 and 2020.The red dott line is the simulation of soil evaporation (SSE), and the black solid line is the observation of evaporation (OSE).The black vertical line represents precipitation, and the blue vertical line rep sents irrigation 2.4.Soil Evaporation

Figure 4 .
Figure 4. Simulated soil evaporation under different treatments: full irrigation (a), mild deficit (c), moderate deficit (e), and severe deficit (g), with wolfberry monoculture and full irrigation (b), mild deficit (d), moderate deficit (f), and severe deficit (h), with intercropping of wolfberry and alfalfa.The black vertical line represents precipitation, and the blue vertical line represents irrigation.

Figure 7 .
Figure 7. Daily precipitation and average temperature during field experiments in 2019 (a) and 2020 (b).

( 1 )
Crop data: basal crop coefficients (K cb ) of the initial, mid-season, and end-season stages under RH min = 45% and u 2 = 2 m•s −1 ; fraction of ground cover f c ; plant height h [m]; root depth Z r [m]; and soil water depletion fraction for non-stress (p).(2) Daily meteorological data: precipitation, u 2 [m•s −1 ], minimum and maximum air temperature, T min and T max [ • C], RH min [%] and P [mm]; reference evapotranspiration, ET 0 [mm] by calculation.(3) Soil data: the total available water, TAW [mm•m −1 ] refers to the capacity of a soil to retain water available to plants; effective depth of the evaporation layer, Z e [m]; total evaporable water, TEW [mm] and readily evaporable water REW [mm].(4) Irrigation data: fraction of soil wetted by irrigation, irrigation depth, and irrigation date.

Table 1 .
Initial and calibrated parameters used in SIMDualKc model.

Table 2 .
Parameters of the calibrated and verified soil moisture content data.

Table 3 .
Parameters of the calibrated and verified soil evaporation data.

Table 4 .
Transpiration (T) and transpiration ratio (T/ET) for each growth stage in the monoculture and intercropping treatments.

Table 5 .
Main soil hydraulic properties.

Table 5 .
Main soil hydraulic properties.