The Loess Plateau (LP) is located in the middle reaches of the Yellow River basin (33°43^′–41°16^′ N, 104°54^′–114°33^′ E) and represents a typical arid and semi-arid region in northern China (Fig. 1). The region covers an area of approximately 640 000 km² and is characterized by complex topography, with elevation generally decreasing from northwest to southeast and an overall elevation range exceeding 3000 m. The Loess Plateau is dominated by a temperate continental monsoon climate. The mean annual air temperature and precipitation are approximately 9.4 °C and 593.7 mm, respectively. Precipitation exhibits strong seasonality, with about 60 %–70 % of the annual total occurring between June and September, whereas winter precipitation is limited and contributes negligibly to the annual water budget (Wang et al., 2023). Under the influence of the East Asian monsoon system, precipitation is primarily driven by large-scale moisture transport, which is subsequently modified by regional land–atmosphere recycling processes. Precipitation over the Loess Plateau also shows pronounced interannual variability, with wet years receiving several times more precipitation than dry years and extreme rainfall events occurring frequently (Liu et al., 2025a).

Given the strong seasonality of both precipitation and vegetation activity, this study focuses on July and August as the primary analysis period. These two months correspond to the peak growing season on the Loess Plateau, during which leaf area index (LAI) and fractional vegetation cover generally reach their annual maxima, solar radiation is relatively strong, and soil moisture conditions remain favorable for active land–atmosphere exchange. As a result, evapotranspiration and vegetation–atmosphere interactions are most pronounced during this period (Wang et al., 2019; Li et al., 2024; Guo et al., 2025). Focusing on July and August therefore allows this study to better capture the period when vegetation-associated moisture recycling is expected to be strongest, providing a targeted basis for evaluating vegetation-associated hydrological effects during the peak growing season in this water-limited region. This seasonal focus is intended to isolate the period of strongest vegetation–atmosphere coupling rather than to represent the full annual water balance.

Monthly evapotranspiration (ET) and potential evapotranspiration (ET₀) data for July and August during 2000–2022 were obtained from the Global Land Evaporation Amsterdam Model (GLEAM, v3.6a) dataset, with a spatial resolution of 0.25° (Martens et al., 2017). GLEAM ET has been extensively evaluated against eddy covariance (EC) observations across China and shows good agreement at the monthly scale, with an average correlation coefficient of 0.86 (Yang et al., 2017).

Monthly precipitation data for July and August during 2000–2022 were derived from the ERA5 reanalysis dataset at a spatial resolution of 0.25° (Hersbach et al., 2020). ERA5 precipitation has been shown to capture the annual and seasonal variability of observed precipitation over China reasonably well (Jiao et al., 2021) and exhibits a strong correlation (r=0.92) with official national statistics at the 99 % confidence level (An et al., 2025). Before subsequent analyses, the ERA5 precipitation and GLEAM ET datasets were temporally aligned for July–August during 2000–2022 and resampled to a common 0.25° grid using nearest-neighbor interpolation to ensure spatial consistency.

For atmospheric moisture tracking, hourly meteorological variables from the ERA5 dataset were collected for June–August during 2000–2022. These variables included air temperature, precipitation, specific humidity, wind components, surface pressure, dew-point temperature, total column water vapor, and surface evaporation. Specific humidity and wind fields were extracted for 17 vertical pressure levels ranging from 1000 to 100 hPa, all at a spatial resolution of 0.25°. ERA5 data have been widely used in previous moisture-tracking studies over China and have shown robust performance (An et al., 2025; Xie et al., 2024; Zhang et al., 2024). Considering the atmospheric residence time of water vapor (van der Ent et al., 2013), meteorological data from June were included as auxiliary input to enable more complete tracking of the moisture contributing to precipitation in July and August.

Vegetation conditions were characterized using the High-Quality Leaf Area Index (HiQ-LAI) dataset, a globally consistent LAI product derived from MODIS LAI Collection 6.1 using a spatiotemporal information compositing algorithm (STICA) (Yan et al., 2024). Compared with the original MODIS LAI product, HiQ-LAI improves temporal stability and spatial consistency by integrating pixel-level quality information and spatiotemporal correlations while preserving the underlying radiative transfer framework. The dataset provides global coverage at an 8 d temporal resolution and a spatial resolution of 5 km for 2000–2022, and has been validated against ground-based observations, showing reduced noise and improved phenological consistency relative to the original MODIS LAI. In this study, HiQ-LAI data for July and August were extracted and aggregated to characterize vegetation conditions during the peak growing season over the Loess Plateau. The HiQ-LAI dataset is publicly available via Zenodo (Yan et al., 2024).

Additional variables used in the random forest modeling of the Budyko parameter n were derived from the ERA5-Land and ERA5 reanalysis datasets. Four-layer soil water content (swvl1–swvl4), near-surface wind components (u and v), and surface temperature were obtained from ERA5-Land at a spatial resolution of 0.1°. Surface runoff (src) and surface shortwave radiation (ssr) were obtained from the ERA5 hourly single-level dataset at a spatial resolution of 0.25°. Vegetation cover (VC) was characterized using the China regional fractional vegetation cover (FVC) dataset at 250 m resolution (Gao et al., 2022). All variables were spatially aggregated to a common 0.25° grid to maintain consistency with the Budyko framework.

This study employs an integrated framework combining atmospheric moisture tracking, vegetation-based diagnostics, and hydrological modeling to examine vegetation-associated precipitation recycling and its hydrological implications over the Loess Plateau during the peak growing season (Fig. 2a).

Atmospheric moisture tracking based on WAM-2layers for the peak growing season (July–August) during 2000–2022 shows that precipitation over the Loess Plateau (LP) is sustained by both external moisture advection and internal moisture recycling. The mean amount of precipitation originating from evaporation at individual grid cells and returning to the LP (Pv0) is 2.05 mm, indicating that terrestrial evaporation contributes measurably to regional precipitation. The dominant moisture transport pathways supplying precipitation to the LP exhibit a clear bimodal pattern (Fig. 3a). One pathway is associated with the mid-latitude westerlies, through which moisture from Central and Western Asia is transported eastward toward the LP. The other is related to southerly transport driven by the southwest and southeast monsoons, which carries moisture from the Indian and Pacific Oceans northward into the study region.

Quantitative source attribution further shows that land-sourced moisture dominates the precipitation supply to the LP. On average, terrestrial moisture contributes approximately 71.0 mm of precipitation, accounting for 86.5 % of total growing-season precipitation, whereas oceanic moisture contributes about 11.1 mm, corresponding to 13.5 % (Fig. 3b). This land-sourced moisture signal reflects moisture that has undergone substantial transformation through repeated evaporation–precipitation cycles during inland transport, rather than exclusively locally generated evaporation. Such a pattern suggests progressive modification of monsoon-derived moisture along continental transport pathways within a monsoon–inland transition zone, where ocean-origin vapor is continuously mixed with terrestrial evapotranspiration inputs before reaching the Loess Plateau.

In terms of recycling pathways, externally advected precipitation accounts for 70.3 mm (85.6 %), while internally recycled precipitation contributes 11.8 mm, representing 14.4 % of total precipitation over the LP. During the study period, internally recycled precipitation shows a trend of -0.073 mm yr⁻¹ (Fig. 3c). This result indicates a gradual weakening of the contribution of internal moisture recycling to growing-season precipitation over the LP.

To characterize the spatial patterns of vegetation conditions relevant to moisture recycling over the Loess Plateau, the vegetation-weighted index LAI_w was analyzed at the grid scale. The spatial distribution and temporal trends of LAI_w are shown in Fig. 4a–b. The multi-year mean LAI_w over the Loess Plateau is 0.45, and approximately 90.6 % of grid cells exhibit LAI_w values lower than 1, mainly in the northwestern part of the plateau where vegetation cover is relatively sparse. In contrast, grid cells with LAI_w values greater than 4 account for only 2.9 % of the region and are concentrated in localized areas of the southeastern Loess Plateau. At the regional scale, 48.7 % of grid cells show increasing trends in LAI_w, whereas 26.5 % show decreasing trends. Spatially, increases in LAI_w are mainly located in areas that already exhibit relatively high LAI_w values (Fig. 4b).

To further evaluate the statistical association between vegetation conditions and recycled precipitation, a grid-based multiple linear regression was performed with recycled precipitation (Prec) as the dependent variable and LAI_w, near-surface air temperature, and surface shortwave radiation as explanatory variables. Here, Prec represents the amount of precipitation over the Loess Plateau that originates from evaporation at a given grid cell and subsequently returns to the plateau through atmospheric moisture transport. The regional mean of the sensitivity coefficient of LAI_w (a) is 2.34, and 57.5 % of grid cells exhibit positive a values, indicating that higher LAI_w is statistically associated with greater recycled precipitation across much of the study area (Fig. 4c–d).

Figure 4e shows the interannual evolution of the vegetation-associated contribution to the evaporation–precipitation (E–P) cycle, decomposed into positive, negative, and net components. Over the study period, the mean positive contribution is 1.16 %, with a weak increasing trend of 0.007 % yr⁻¹. In contrast, the mean negative contribution is -0.97 %, with a trend of -0.026 % yr⁻¹. As a combined result of these opposing components, the net vegetation-associated contribution to the E–P cycle remains small, with a mean value of 0.12 % and a slight trend of -0.019 % yr⁻¹. Overall, these results show that positive and negative vegetation-associated influences on internal moisture recycling largely offset each other at the regional scale.

A pronounced drying tendency is observed over the Loess Plateau during the growing season from 2000 to 2022 (Fig. 5a). Interannual variations in precipitation (P), evapotranspiration (E), surface water availability (SWA), and the aridity index (AI) show that regional precipitation exhibits a non-significant trend of -0.0043 mm yr⁻¹, whereas evapotranspiration shows a significant increasing trend of 0.31 mm yr⁻¹. The contrasting changes in precipitation and evapotranspiration correspond to increasingly dry hydroclimatic conditions over the study period. Although the interannual variations in SWA and AI generally follow the temporal pattern of precipitation, both indices exhibit significant downward trends under increasing evapotranspiration, with SWA decreasing at a rate of 0.34 mm yr⁻¹ and AI decreasing at a rate of 0.008 yr⁻¹.

The spatial distribution of mean SWA and its temporal trends also show clear heterogeneity across the Loess Plateau (Fig. 5b–c). Higher SWA values are mainly distributed in the southeastern part of the region, whereas persistently low SWA values characterize the northwestern arid areas. Declining SWA trends occur in both relatively high-SWA and low-SWA regions, whereas increasing trends are mainly concentrated in the southwestern Loess Plateau. Although 51 % of grid cells show increasing trends in SWA and 49 % show decreasing trends, the pattern becomes asymmetric when statistical significance is considered. Among the 56.6 % of grid cells with significant SWA trends (P<0.05), 40 % exhibit significant declines, whereas only 16.6 % show significant increases. This indicates that the most robust regional signal is a widespread decline in SWA.

To quantify the hydrological implications of vegetation-associated recycled precipitation (P0), the P0 term estimated in Sect. 3.2 was subtracted from total precipitation and the resulting water-balance components were evaluated within the Budyko–RF framework. The results (Fig. 5d) show that the mean relative contributions of vegetation-associated recycled precipitation to P, E, and SWA are 0.059 %, 0.02 %, and 0.107 %, respectively. A consistent response pattern is observed across the study region: grid cells with positive P0 contributions to precipitation also show higher E and SWA, whereas negative P0 contributions are accompanied by concurrent reductions in both E and SWA. Overall, the positive hydrological contribution of P0 is relatively more evident in drier areas with lower SWA.

Figure 5e further shows a systematic dependence of the SWA contribution associated with vegetation-associated recycled precipitation on LAI_w. Under low-LAI_w conditions (LAIw<1), the contribution to SWA remains consistently positive, with a mean value of 0.2855 % and no negative values across all years, although it shows a trend of -0.0023 % yr⁻¹. As LAI_w increases, this positive contribution diminishes and gradually shifts toward negative values. For intermediate LAI_w conditions (1≤LAIw≤2.5), the mean contribution is close to zero (0.02565 %) and shows substantial interannual variability, with negative values occurring in 11 years; however, a modest positive trend (0.0035 % yr⁻¹) is still observed. In contrast, high-LAI_w conditions (LAIw≥2.5) are characterized by persistently negative contributions to SWA, with mean values of -0.0361 % for 2.5≤LAIw≤4 and -0.1495 % for LAIw>4, and the strongest negative trend (-0.0092 % yr⁻¹) occurs in the highest-LAI_w class. Overall, these results indicate a systematic shift from weakly positive to increasingly negative SWA contributions along the LAI_w gradient.

Recent global studies show that surface water availability (SWA) has increased over much of the terrestrial surface. This is mainly because precipitation has risen faster than evapotranspiration across more than 53 % of global land areas (Cui et al., 2022). In contrast, our results show a decline in SWA over the Loess Plateau during the study period (Fig. 5a). In this arid and semi-arid region, evapotranspiration increased faster than precipitation, reducing available surface water. Similar drying trends have been widely reported in dryland. Climate warming is expected to increase drought frequency and severity (Stevenson et al., 2022; Tabari et al., 2021). In addition, limited precipitation and enhanced evaporation often lead to declining regional water resources (Lai et al., 2022). Together, these findings indicate that water availability in in drylands is strongly constrained by low precipitation and high evaporative demand.

Under this hydroclimatic background, precipitation over the Loess Plateau is strongly controlled by large-scale monsoon moisture transport. It is continuously modified by land–atmosphere exchanges along its inland pathway. The WAM-2layers results show that Central China is the predominantly continental-influenced moisture source for precipitation over the Loess Plateau (Fig. 3a). This is consistent with previous Lagrangian studies. These studies indicate that summer moisture mainly comes from Central–Eastern China and northwestern continental Asia. In contrast, direct oceanic contributions from the East China Sea, South China Sea, Bay of Bengal, and Arabian Sea are minor, together accounting for only about 10 % of total precipitation (Hu et al., 2018).

The dominance of continental moisture sources reflects progressive changes in monsoon vapor during inland transport. Along this pathway, repeated evaporation–precipitation cycles and boundary layer recycling gradually reshape the original oceanic moisture signal. Recent studies also show that vegetation greening enhances evapotranspiration and strengthens regional moisture recycling (Cao et al., 2024). In this process, the land surface acts both as a moisture source and an intermediate contributor to atmospheric moisture through evapotranspiration and recycling processes. This reinforces the continental component in the monsoon–land coupled system. Therefore, the high fraction of land-sourced precipitation (86.5 %) should not be interpreted as a weak oceanic influence. Instead, it reflects the combined effects of monsoon moisture transport and terrestrial recycling. The East Asian monsoon provides the initial moisture input, while continental recycling determines the final precipitation composition over the Loess Plateau.

Within this externally controlled system, internal precipitation recycling contributes about 14.4 % of total precipitation during the study period (Fig. 3b–c). This value is higher than those reported in some previous studies, such as 7 % from the UTrack Lagrangian moisture-tracking model (Wang et al., 2023) and 9.16±0.35 % to 10.18±0.28 % from the Dynamic Precipitation Recycling Method (Tian et al., 2022). These differences mainly reflect methodological approaches and climatic context rather than physical inconsistency. Most previous studies use annual or long-term averages. In contrast, this study focuses on July–August, when vegetation activity, evapotranspiration, and land–atmosphere coupling are strongest over the Loess Plateau (Li et al., 2024; Guo et al., 2025).

This seasonal amplification agrees with recent findings that vegetation and land use changes regulate atmospheric moisture recycling by affecting evapotranspiration and atmospheric moisture residence time (Wang et al., 2023). In this framework, precipitation recycling is not a fixed ratio. It responds dynamically to land surface conditions and atmospheric moisture supply. Similarly, afforestation has been shown to strengthen land–atmosphere coupling and increase both external and internal moisture contributions to precipitation (Tian et al., 2022). Therefore, the relatively high recycling ratio in this study reflects enhanced seasonal coupling during the peak growing season, rather than overestimation. These results indicate that precipitation over the Loess Plateau is jointly controlled by monsoon moisture transport and terrestrial recycling processes.

Under water-limited conditions, vegetation affects regional water cycling through two competing pathways. First, enhanced vegetation activity increases evapotranspiration and strengthens land–atmosphere coupling. This can promote the return of locally evaporated moisture as precipitation (Sterling et al., 2013). Second, vegetation growth increases terrestrial water consumption. This reduces soil moisture and limits water availability for surface and subsurface hydrological processes (Piao et al., 2020). These two pathways coexist. As a result, vegetation-associated effects on internal moisture recycling do not translate into a uniformly positive impact on surface water availability.

Recent work by Tian et al. (2026) provides a complementary system-level view. They show that vegetation greening can enhance precipitation efficiency through albedo-driven energy feedbacks. These feedbacks increase surface radiation, atmospheric instability, and moisture-to-precipitation conversion efficiency. Their results also indicate that greening can simultaneously increase evapotranspiration and streamflow in semi-arid regions. This suggests a coupled energy–moisture mechanism that intensifies the atmospheric water cycle. The results are broadly consistent with the energy–moisture coupling perspective proposed by Zhang (2025), particularly in terms of enhanced atmospheric moisture recycling under greening conditions. We observe strengthened land–atmosphere coupling and enhanced internally recycled moisture. However, the effect on surface water availability depends on temporal scale and hydrological closure assumptions.At the regional scale, vegetation-associated effects on the evaporation–precipitation cycle largely cancel out over the long term. Vegetation can enhance precipitation during active growth periods. However, this is accompanied by increased evapotranspiration. As also suggested by Tian et al. (2026), this represents a redistribution of water within the system rather than a net gain in water resources. Consequently, the net effect on regional water availability is very small (Figs. 4e, 5d). This small value should be interpreted carefully. Here, we only quantify the vegetation contribution to precipitation generated from local evaporation and recycled within the Loess Plateau. Internally recycled precipitation accounts for only 14.4 % of total growing-season precipitation (Fig. 3c). Therefore, even a detectable vegetation effect on this component becomes strongly diluted when expressed relative to total precipitation and total surface water availability.

In addition, the regional mean signal reflects strong compensation between positive and negative vegetation effects. The strongest positive and negative contributions reach 12.2 % and -32.2 % of the internally recycled precipitation, respectively (Figs. 3c and 6a). Thus, the small net effect results not only from the limited share of internal recycling, but also from strong cancellation across vegetation conditions. This compensation is consistent with the nonlinear energy–moisture coupling described by Tian et al. (2026). In that framework, enhanced moisture conversion efficiency coexists with increased evaporative demand. This leads to spatially heterogeneous hydrological responses rather than uniform changes.

A small net effect does not imply that vegetation–atmosphere feedbacks are negligible. Instead, it indicates that under the dry conditions in the Loess Plateau, additional precipitation generated through recycling is generally insufficient to offset increased evapotranspiration caused by vegetation growth. As a result, vegetation does not substantially increase net surface water availability. Instead, it primarily reflects a redistribution of water within a tightly constrained land–atmosphere system.

In this context, the hydrological benefit of vegetation-associated recycling is strongly limited by background hydroclimatic conditions. This is different from humid or energy-limited regions, where vegetation-associated recycling more easily translates into net water gains (Cui et al., 2022; Meier et al., 2021). The comparison highlights the key role of climate setting in regulating vegetation–water coupling efficiency. It also suggests that small mean effects should not be interpreted as hydrological irrelevance. Rather, they reflect strong compensatory feedbacks under water-limited conditions.More importantly, vegetation effects on surface water availability show a clear nonlinear pattern along the LAI_w gradient. Two lines of evidence support this. First, the contribution of vegetation-associated recycled precipitation shifts from positive under low-LAI_w to negative under high-LAI_w. This indicates that the hydrological benefit of recycling is mainly limited to sparse vegetation. It weakens when vegetation density exceeds a moderate level (Figs. 5e and 6a). Second, the temporal evolution of positive and negative components shows asymmetric behavior. Positive contributions are small and intermittent, while negative contributions become more persistent with increasing LAI_w (Figs. 4e, 5e and 6b). This pattern is consistent with the limited water storage capacity and strong evaporative sensitivity in the Loess Plateau.

Within the broader framework of Tian et al. (2026), this nonlinear behavior reflects a trade-off between enhanced precipitation efficiency and increased evaporative demand under higher vegetation density. Under such conditions, further vegetation growth is more likely to intensify water consumption than to produce sustained gains in surface water availability. Therefore, the key scientific value of the small regional mean effect lies not in its magnitude, but in revealing the conditions under which vegetation-associated recycling shifts from weakly beneficial to increasingly negative under higher LAI_w conditions.

Because LAI_w is a derived diagnostic indicator rather than a directly observed physical variable, its formulation may raise concerns about whether the identified hydrological patterns depend on the weighting scheme. In particular, weighting LAI by normalized actual evapotranspiration could, in principle, influence the apparent strength of vegetation–hydrological coupling and thereby affect the inferred threshold-like behavior.

To test this, we repeated the full analytical workflow using the original, unweighted peak-growing-season LAI instead of LAI_w. This parallel analysis included the multiple linear regression for Prec and the Budyko-based attribution of hydrological effects. We then compared vegetation contributions to surface water availability (CSWA) between the LAI-based and LAI_w-based frameworks at both the regional mean scale and across vegetation classes defined in Sect. 3.3.

The results show that the main findings are highly robust to the choice of vegetation indicator. Interannual variations in P0 and CSWA are nearly identical between the two frameworks, with Pearson correlation coefficients of 0.99996 and 0.99995, respectively (Fig. 7). The multi-year mean values also show only minor differences, with an absolute difference of 0.013 %. This high consistency indicates that the observed hydrological response is not an artifact of the LAI_w weighting scheme. Instead, LAI_w mainly acts as a diagnostic rescaling of the vegetation gradient under different evaporative conditions. The dominant signal is still controlled by vegetation structure.

Importantly, the threshold-like behavior should be interpreted as a real regime transition driven by the interaction between vegetation density and evaporative demand. It is not produced by the LAI_w formulation. Rather, LAI_w helps to reveal how vegetation effects on moisture recycling and surface water availability change under different atmospheric conditions. This suggests that vegetation–water interactions in the Loess Plateau operate within a constrained hydrological regime. The system shifts from a weakly coupled state with slight positive effects to a strongly evaporative-dominated state with negative effects. This transition reflects limited water storage capacity and high evaporative sensitivity in the region.

From a management perspective, this shift implies an implicit water-limited carrying capacity for vegetation restoration. Beyond this level, further increases in vegetation density are more likely to increase evapotranspiration than to enhance surface water availability. However, this should not be interpreted as a fixed threshold for restoration. Instead, it represents a context-dependent response governed by hydroclimatic conditions in dryland environments.

This study provides a quantitative assessment of vegetation-associated precipitation recycling and its contribution to surface water availability during the peak growing season over the Loess Plateau. Several limitations should be acknowledged.

First, the seasonal focus on July–August. This period represents a peak-coupling stage within the monsoon-dominated growing season from June to September. The growing season accounts for 60 %–70 % of annual precipitation and coincides with the annual maxima of precipitation, evapotranspiration, and land–atmosphere coupling (Liang et al., 2015; Zhang et al., 2018, 2022b). During July–August, vegetation is fully developed and atmospheric moisture exchange is strong. In contrast, coupling is weaker in early summer because vegetation is not yet fully developed. It also weakens in late summer because of soil moisture depletion (Feng et al., 2016; Liang et al., 2015). Therefore, our analysis captures the period of strongest vegetation–atmosphere interaction rather than a full-season or annual mean state. The identified recycling signals should thus be interpreted as the upper bound of seasonal coupling intensity. Annual-scale water balance assessments, by contrast, reflect integrated hydrological responses across multiple seasonal phases (Wang et al., 2016).

Second, the analysis is limited to local hydrological responses within the Loess Plateau. It does not explicitly consider downwind effects caused by vegetation-induced moisture export. Vegetation-enhanced evapotranspiration may reduce local water availability. At the same time, it may contribute to precipitation in downstream regions through atmospheric transport (Greve et al., 2025). A more complete assessment would require a coupled local–downwind framework. Such a framework could quantify the broader hydroclimatic footprint of vegetation-induced moisture redistribution.

Third, uncertainty arises from both the input datasets and the attribution framework. Recycled precipitation and hydrological responses were derived from multiple reanalysis and remote-sensing products. These uncertainties may propagate through the moisture-tracking analysis and the Budyko–RF attribution framework. In addition, LAI_w is a diagnostic indicator rather than a directly observed physical variable. The P0 term also represents a statistical attribution of vegetation-associated recycled precipitation, rather than a strictly causal estimate from perturbation experiments. Therefore, the counterfactual analysis isolates a first-order hydrological effect. It does not resolve secondary feedbacks involving soil moisture, vegetation dynamics, or atmospheric circulation responses.

Future work should address these limitations in three ways. First, the analysis should be extended from the peak growing season to full seasonal and annual scales. Second, downwind moisture transport should be explicitly incorporated to capture the broader hydroclimatic footprint of vegetation. Third, uncertainty characterization should be improved for both input datasets and diagnostic attribution frameworks. Together, these efforts would help clarify how vegetation–atmosphere feedbacks evolve across spatial gradients and temporal scales in water-limited environments.