Spatiotemporal origin of soil water taken up by vegetation

Dynamic process of ecosystem water use efficiency and response to drought in the Yellow River Basin, China

Ecosystem water use efficiency (WUE) is a crucial indicator of the impact of climate change on terrestrial ecosystems, reflecting the balance between biological processes (photosynthesis and transpiration) and physical processes (evapotranspiration). However, the response mechanisms and driving processes of WUE to drought remain to be further understood. In this study, we analyzed the spatial and temporal dynamics and response mechanisms of WUE in the Yellow River Basin (YRB) using data on Gross Primary Productivity (GPP), Evapotranspiration (ET) and Standardized Precipitation Evapotranspiration Index (SPEI), which revealed the cumulative effect of drought on WUE and assessed the ecosystem's resilience. The study results showed that (1) GPP, ET and WUE in the YRB exhibited a significant increasing trend, with 63.04 % of the area showing a marked increase in WUE. (2) GPP was the dominant factor influencing WUE in 65.36 % of the area, particularly in cropland and grassland, while ET was more influential in forested areas. Vapor pressure deficit (VPD) was identified as the principal driver affecting vegetation GPP in semi-arid and semi-humid regions of the YRB. In contrast, soil moisture (SM) was the limiting factor in arid areas. (3) 71.00 % of the WUE in the basin was affected by drought cumulative effects, with an average cumulative duration of 4.5 months. Arid regions experienced the most extended duration of 7.29 months, compared to 3.05 months in semi-humid regions. (4) 74.85 % of the regional ecosystems exhibited ecological resilience to drought, particularly in the source areas of the western basin of the YRB. Shrublands have the highest drought resilience among vegetation types, while grasslands have the lowest. The resilience of each climatic zone was in the order of semi-humid, semi-arid, and arid order. This study comprehensively analyzed of the spatial and temporal dynamics and response mechanisms of WUE in the YRB, offering a new perspective and scientific basis for understanding and predicting the ecosystem response to climate change.

Dataset of soil hydraulic parameters in the Yellow River Basin based on in situ deep sampling

Soil hydraulic parameters are vital for precisely characterizing soil hydrological processes, which are critical indicators for regulating climate change effects on terrestrial ecosystems and governing feedbacks between water, energy, and carbon-nitrogen cycles. Although many studies have integrated comprehensive soil datasets, data quality and cost challenges result in data completeness deficiencies, especially for deep soil information. These gaps not only impede methodological endeavours but also constrain soil parameter-based ecosystem process studies spanning from local profiles to global earth system models. We established a soil dataset across the entire Yellow River Basin (YRB) (795,000 km2) using high-density in situ field sampling. This observation-based dataset contains records of soil texture (2924), bulk density (2798), saturated hydraulic conductivity (2782), and water retention curve parameters (1035) down to a maximum depth of 5 m. This dataset, which extends the recorded data range for deep soil hydraulic parameters, is valuable as a direct data resource for environmental, agronomical and hydrological studies in the YRB and regions with similar pedological and geological backgrounds around the world.

Long-term species-level observations indicate the critical role of soil moisture in regulating China's grassland productivity relative to phenological and climatic factors

As a sensitive indicator of climate change and a key variable in ecosystem surface-atmosphere interaction, vegetation phenology, and the growing season length, as well as climatic factors (i.e., temperature, precipitation, and sunshine duration) are widely recognized as key factors influencing vegetation productivity. Recent studies have highlighted the importance of soil moisture in regulating grassland productivity. However, the relative importance of phenology, climatic factors, and soil moisture to plant species-level productivity across China's grasslands remains poorly understood. Here, we use nearly four decades (1981 to 2018) of in situ species-level observations from 17 stations distributed across grasslands in China to examine the key mechanisms that control grassland productivity. The results reveal that soil moisture is the strongest determinant of the interannual variability in grassland productivity. In contrast, the spring/autumn phenology, the length of vegetation growing season, and climate factors have relatively minor impacts. Generally, annual aboveground biomass increases by 3.9 to 25.3 g∙m2 (dry weight) with a 1 % increase in growing season mean soil moisture across the stations. Specifically, the sensitivity of productivity to moisture in wetter and colder environments (e.g., alpine meadows) is significantly higher than that in drier and warmer environments (e.g., temperate desert steppes). In contrast, the sensitivity to the precipitation of the latter is greater than the former. The effect of soil moisture is the most pronounced during summer. Dominant herb productivity is more sensitive to soil moisture than the others. Moreover, multivariate regression analyses show that the primary climatic factors and their attributions to variations in soil moisture differ among the stations, indicating the interaction between climate and soil moisture is very complex. Our study highlights the interspecific difference in the soil moisture dependence of grassland productivity and provides guidance to climate change impact assessments in grassland ecosystems.

A unifying modelling of multiple land degradation pathways in Europe

Land degradation is a complex socio-environmental threat, which generally occurs as multiple concurrent pathways that remain largely unexplored in Europe. Here we present an unprecedented analysis of land multi-degradation in 40 continental countries, using twelve dataset-based processes that were modelled as land degradation convergence and combination pathways in Europe's agricultural (and arable) environments. Using a Land Multi-degradation Index, we find that up to 27%, 35% and 22% of continental agricultural (~2 million km2) and arable (~1.1 million km2) lands are currently threatened by one, two, and three drivers of degradation, while 10-11% of pan-European agricultural/arable landscapes are cumulatively affected by four and at least five concurrent processes. We also explore the complex pattern of spatially interacting processes, emphasizing the major combinations of land degradation pathways across continental and national boundaries. Our results will enable policymakers to develop knowledge-based strategies for land degradation mitigation and other critical European sustainable development goals.

Tree water uptake patterns across the globe

Plant water uptake from the soil is a crucial element of the global hydrological cycle and essential for vegetation drought resilience. Yet, knowledge of how the distribution of water uptake depth (WUD) varies across species, climates, and seasons is scarce relative to our knowledge of aboveground plant functions. With a global literature review, we found that average WUD varied more among biomes than plant functional types (i.e. deciduous/evergreen broadleaves and conifers), illustrating the importance of the hydroclimate, especially precipitation seasonality, on WUD. By combining records of rooting depth with WUD, we observed a consistently deeper maximum rooting depth than WUD with the largest differences in arid regions - indicating that deep taproots act as lifelines while not contributing to the majority of water uptake. The most ubiquitous observation across the literature was that woody plants switch water sources to soil layers with the highest water availability within short timescales. Hence, seasonal shifts to deep soil layers occur across the globe when shallow soils are drying out, allowing continued transpiration and hydraulic safety. While there are still significant gaps in our understanding of WUD, the consistency across global ecosystems allows integration of existing knowledge into the next generation of vegetation process models.

Temperate woody species across the angiosperm phylogeny acquire tolerance to water deficit stress during the growing season

Understanding the capacity of temperate trees to acclimate to limited soil water has become essential in the face of increasing drought risk due to climate change. We documented seasonal - or phenological - patterns in acclimation to water deficit stress in stems and leaves of tree species spanning the angiosperm phylogeny. Over 3 yr of field observations carried out in two US arboreta, we measured stem vulnerability to embolism (36 individuals of 7 Species) and turgor loss point (119 individuals of 27 species) over the growing season. We also conducted a growth chamber experiment on 20 individuals of one species to assess the mechanistic relationship between soil water restriction and acclimation. In three-quarters of species measured, plants became less vulnerable to embolism and/or loss of turgor over the growing season. We were able to stimulate this acclimatory effect by withholding water in the growth chamber experiment. Temperate angiosperms are capable of acclimation to soil water deficit stress, showing maximum vulnerability to soil water deficits following budbreak and becoming more resilient to damage over the course of the growing season or in response to simulated drought. The species-specific tempo and extent of this acclimatory potential constitutes preadaptive climate change resilience.

Analysis of soil water use by exotic and native vegetation in a semi-arid area and their associated interspecific competition

The exotic vegetation used in dryland vegetation restoration projects is characterized by its fast-growing and deep-rooted system, which enables it to expedite the restoration of ecosystem functions and enhance biodiversity. However, the interspecific relationship between exotic and native vegetation and soil water uptake in these restored ecosystems remains unclear, limiting our ability to evaluate the succession process and sustainability of restored ecosystems. In this study, stable isotope techniques and a proportional similarity index were used to investigate soil water use strategies and interspecific relationships between exotic and native vegetation. The results showed significant differences between the soil water use strategies of both exotic and native vegetation between seasons and species, where the proportions of deep soil water (30-100 cm) used by exotic shrubs (Caragana korshinskii) and exotic grass (Medicago sativa) were significantly higher than those used by the co-occurring native grass (Stipa bungeana) (p < 0.05). As soil water storage declined, exotic vegetation increased its utilization of deep soil water, whereas native grasses relied more on surface water (0-10 cm). This suggests that deep-rooted exotic vegetation has greater adaptability and access to water resources than shallow-rooted native vegetation. However, a prolonged decline in soil water storage led to increased competition for surface soil water (0-30 cm) between the exotic and native vegetation. This may increase the risk of degradation of exotic vegetation, particularly in situations with lower soil water content in the deep layers. Overall, this study highlights the variation in water-use strategies and interspecies relationships between exotic and native vegetation and their implications for ecosystem succession, which provides valuable insights for developing future vegetation restoration strategies and managing restored ecosystems.

Rooting depth and xylem vulnerability are independent woody plant traits jointly selected by aridity, seasonality, and water table depth

Evolutionary radiations of woody taxa within arid environments were made possible by multiple trait innovations including deep roots and embolism-resistant xylem, but little is known about how these traits have coevolved across the phylogeny of woody plants or how they jointly influence the distribution of species. We synthesized global trait and vegetation plot datasets to examine how rooting depth and xylem vulnerability across 188 woody plant species interact with aridity, precipitation seasonality, and water table depth to influence species occurrence probabilities across all biomes. Xylem resistance to embolism and rooting depth are independent woody plant traits that do not exhibit an interspecific trade-off. Resistant xylem and deep roots increase occurrence probabilities in arid, seasonal climates over deep water tables. Resistant xylem and shallow roots increase occurrence probabilities in arid, nonseasonal climates over deep water tables. Vulnerable xylem and deep roots increase occurrence probabilities in arid, nonseasonal climates over shallow water tables. Lastly, vulnerable xylem and shallow roots increase occurrence probabilities in humid climates. Each combination of trait values optimizes occurrence probabilities in unique environmental conditions. Responses of deeply rooted vegetation may be buffered if evaporative demand changes faster than water table depth under climate change.

Southern Hemisphere dominates recent decline in global water availability

Global land water underpins livelihoods, socioeconomic development, and ecosystems. It remains unclear how water availability has changed in recent decades. Using an ensemble of observations, we quantified global land water availability over the past two decades. We show that the Southern Hemisphere has dominated the declining trend in global water availability from 2001 to 2020. The significant decrease occurs mainly in South America, southwestern Africa, and northwestern Australia. In the Northern Hemisphere, the complex regional increasing and decreasing trends cancel each other, resulting in a negligible hemispheric trend. The variability and trend in water availability in the Southern Hemisphere are largely driven by precipitation associated with climate modes, particularly the El Niño-Southern Oscillation. This study highlights their dominant role in controlling global water availability.

Diagnosing evapotranspiration responses to water deficit across biomes using deep learning

Accounting for water limitation is key to determining vegetation sensitivity to drought. Quantifying water limitation effects on evapotranspiration (ET) is challenged by the heterogeneity of vegetation types, climate zones and vertically along the rooting zone. Here, we train deep neural networks using flux measurements to study ET responses to progressing drought conditions. We determine a water stress factor (fET) that isolates ET reductions from effects of atmospheric aridity and other covarying drivers. We regress fET against the cumulative water deficit, which reveals the control of whole-column moisture availability. We find a variety of ET responses to water stress. Responses range from rapid declines of fET to 10% of its water-unlimited rate at several savannah and grassland sites, to mild fET reductions in most forests, despite substantial water deficits. Most sensitive responses are found at the most arid and warm sites. A combination of regulation of stomatal and hydraulic conductance and access to belowground water reservoirs, whether in groundwater or deep soil moisture, could explain the different behaviors observed across sites. This variety of responses is not captured by a standard land surface model, likely reflecting simplifications in its representation of belowground water storage.

Dietary sources apportionment and health risk assessment for trace elements among residents of the Tethys-Himalayan tectonic domain in Tibet, China

Dietary intake of toxic elements (TEs) and essential trace elements (ETEs) can significantly impact human health. This study collected 302 samples, including 78 food, 104 drinking water, 73 cultivated topsoil, and 47 sedimentary rock from a typical area of Tethys-Himalaya tectonic domain. These samples were used to calculate the average daily dose of oral intake (ADDoral) and assess the health risks of five TEs and five ETEs. The results indicate that grain and meat are the primary dietary sources of TEs and ETEs for local residents. The intake of manganese (Mn) and copper (Cu) is mainly from local highland barley (66.90% and 60.32%, respectively), iron (Fe) is primarily from local grains (75.51%), and zinc (Zn) is mainly from local yak meat (60.03%). The ADDoral of arsenic (As), Mn, Fe and Zn were found to be higher than the maximum oral reference dose in all townships of study area, indicating non-carcinogenic health risks for local residents. Additionally, lead (Pb) and nickel (Ni) in 36.36% townships, and Cu in 81.82% townships were above the maximum oral reference dose, while As posed a carcinogenic risk throughout the study area. The concentrations of As, mercury (Hg), Pb, Mn, Cu Fe and selenium (Se) in grains were significantly correlated with those in soils. Moreover, the average concentrations of As in Proterozoic, Triassic, Jurassic and Cretaceous was 43.09, 12.41, 15.86 and 6.22 times higher than those in the South Tibet shell, respectively. The high concentrations of TEs and ETEs in the stratum can lead to their enrichment in soils, which, in turn, can result in excessive intake by local residents through the food chain and biogeochemical cycles . To avoid the occurrence of some diseases caused by dietary intake, it is necessary to consume a variety of exotic foods, such as high-selenium foods, foreign rice and flour in order to improve the dietary structure.

Deep roots mitigate drought impacts on tropical trees despite limited quantitative contribution to transpiration

Deep rooting is considered a central drought-mitigation trait with vast impact on ecosystem water cycling. Despite its importance, little is known about the overall quantitative water use via deep roots and dynamic shifts of water uptake depths with changing ambient conditions. Knowledge is especially sparse for tropical trees. Therefore, we conducted a drought, deep soil water labeling and re-wetting experiment at Biosphere 2 Tropical Rainforest. We used in situ methods to determine water stable isotope values in soil and tree water in high temporal resolution. Complemented by soil and stem water content and sap flow measurements we determined percentages and quantities of deep-water in total root water uptake dynamics of different tree species. All canopy trees had access to deep-water (max. uptake depth 3.3 m), with contributions to transpiration ranging between 21 % and 90 % during drought, when surface soil water availability was limited. Our results suggest that deep soil is an essential water source for tropical trees that delays potentially detrimental drops in plant water potentials and stem water content when surface soil water is limited and could hence mitigate the impacts of increasing drought occurrence and intensity as a consequence of climate change. Quantitatively, however, the amount of deep-water uptake was low due to the trees' reduction of sap flow during drought. Total water uptake largely followed surface soil water availability and trees switched back their uptake depth dynamically, from deep to shallow soils, following rainfall. Total transpiration fluxes were hence largely driven by precipitation input.

Programmable Microfluidics Enabled by 3D Printed Bionic Janus Porous Matrics for Microfluidic Logic Chips

Numerous structures have been functionally optimized for directional liquid transport in nature. Inspired by lush trees' xylem that enable liquid directional transportation from rhizomes to the tip of trees, a new kind of programmable microfluidic porous matrices using projection micro-stereolithography (PµSL) based 3D printing technique is fabricated. Structural matrices with internal superhydrophilicity and external hydrophobicity are assembled for ultra-fast liquid rising enabled by capillary force. Moreover, the unidirectional microfluidic performance of the bionic porous matrices can be theoretically optimized by adjusting its geometric parameters. Most significantly, the successive programmable flow of liquid in a preferred direction inside the bionic porous matrices with tailored wettability is achieved, validating by a precisely printed liquid displayer and a microfluidic logic chip. The programmable and functional microfluidic matrices promise applications of patterned liquid flow, displayer, logic chip, cell screening, gas-liquid separation, and so on.

A multi-satellite framework to rapidly evaluate extreme biosphere cascades: The Western US 2021 drought and heatwave

The increasing frequency and intensity of climate extremes and complex ecosystem responses motivate the need for integrated observational studies at low latency to determine biosphere responses and carbon-climate feedbacks. Here, we develop a satellite-based rapid attribution workflow and demonstrate its use at a 1-2-month latency to attribute drivers of the carbon cycle feedbacks during the 2020-2021 Western US drought and heatwave. In the first half of 2021, concurrent negative photosynthesis anomalies and large positive column CO2 anomalies were detected with satellites. Using a simple atmospheric mass balance approach, we estimate a surface carbon efflux anomaly of 132 TgC in June 2021, a magnitude corroborated independently with a dynamic global vegetation model. Integrated satellite observations of hydrologic processes, representing the soil-plant-atmosphere continuum (SPAC), show that these surface carbon flux anomalies are largely due to substantial reductions in photosynthesis because of a spatially widespread moisture-deficit propagation through the SPAC between 2020 and 2021. A causal model indicates deep soil moisture stores partially drove photosynthesis, maintaining its values in 2020 and driving its declines throughout 2021. The causal model also suggests legacy effects may have amplified photosynthesis deficits in 2021 beyond the direct effects of environmental forcing. The integrated, observation framework presented here provides a valuable first assessment of a biosphere extreme response and an independent testbed for improving drought propagation and mechanisms in models. The rapid identification of extreme carbon anomalies and hotspots can also aid mitigation and adaptation decisions.

Plant growth forms dictate adaptations to the local climate

Adaptive radiation is a significant driver of biodiversity. Primarily studied in animal systems, mechanisms that trigger adaptive radiations remain poorly understood in plants. A frequently claimed indicator of adaptive radiation in plants is growth form diversity when tied to the occupation of different habitats. However, it remains obscure whether morphological adaptations manifest as growth form diversity per se or as its constituent traits. We use the classic Aeonium radiation from the Canary Islands to ask whether adaptation across climatic space is structured by growth form evolution. Using morphological sampling with site-associated climate in a phylogenetic context, we find that growth forms dictate adaptations to the local environment. Furthermore, we demonstrate that the response of specific traits to analogous environments is antagonistic when growth forms are different. This finding suggests for the first time that growth forms represent particular ecological functions, allowing the co-occurrence of closely related species, being a product of divergent selection during evolution in sympatry.

Coupling water cycle processes with water demand routes of vegetation using a cascade causal modeling approach in arid inland basins

Vegetation degradation is the key cause of land desertification in arid areas. Water stress is one of the most critical factors leading to vegetation degradation. The water needed for vegetation growth is inseparable from the water cycle processes. It is a new scope to reveal the vegetation water demand mechanisms from the water cycle processes. Water cycle processes in arid inland basins can be conceptually separated as RFA (runoff formation area) and RCA (runoff consumption area). In this study, both the water demand mechanisms of natural vegetation and farmland were discovered by creatively constructing the vegetation water demand route model. The TRB (Tarim River Basin), a typical arid inland basin system that RFA is separated from RCA, is considered as the study area. The tendency and relevance of water cycle factors and NDVI were detected. The dominant factors of vegetation growth were identified. According to the interaction causality of water cycle factors and vegetation, the PLS-SEM (partial least squares structural equation models) were constructed in RFA and RCA. Results displayed that SMroot (root-zone soil moisture), groundwater and precipitation were the dominant water sources for natural vegetation in RFA. The water demand for natural vegetation in RCA mainly came from SMroot and that for farmland mainly came from SMsurf (surface soil moisture). New findings showed that blue and green water circulations were more active in RFA than in RCA. Natural vegetation had better adaptability and resilience to water shortages compared with farmland. The higher effect of vegetation on AET (actual evapotranspiration) denoted the better growth status. It is contributed to the rational utilization of water resources in arid basins.

Observation-based assessment of secondary water effects on seasonal vegetation decay across Africa

Local studies and modeling experiments suggest that shallow groundwater and lateral redistribution of soil moisture, together with soil properties, can be highly important secondary water sources for vegetation in water-limited ecosystems. However, there is a lack of observation-based studies of these terrain-associated secondary water effects on vegetation over large spatial domains. Here, we quantify the role of terrain properties on the spatial variations of dry season vegetation decay rate across Africa obtained from geostationary satellite acquisitions to assess the large-scale relevance of secondary water effects. We use machine learning based attribution to identify where and under which conditions terrain properties related to topography, water table depth, and soil hydraulic properties influence the rate of vegetation decay. Over the study domain, the machine learning model attributes about one-third of the spatial variations of vegetation decay rates to terrain properties, which is roughly equally split between direct terrain effects and interaction effects with climate and vegetation variables. The importance of secondary water effects increases with increasing topographic variability, shallower groundwater levels, and the propensity to capillary rise given by soil properties. In regions with favorable terrain properties, more than 60% of the variations in the decay rate of vegetation are attributed to terrain properties, highlighting the importance of secondary water effects on vegetation in Africa. Our findings provide an empirical assessment of the importance of local-scale secondary water effects on vegetation over Africa and help to improve hydrological and vegetation models for the challenge of bridging processes across spatial scales.

Feedback in tropical forests of the Anthropocene

Tropical forests are complex systems containing myriad interactions and feedbacks with their biotic and abiotic environments, but as the world changes fast, the future of these ecosystems becomes increasingly uncertain. In particular, global stressors may unbalance the feedbacks that stabilize tropical forests, allowing other feedbacks to propel undesired changes in the whole ecosystem. Here, we review the scientific literature across various fields, compiling known interactions of tropical forests with their environment, including the global climate, rainfall, aerosols, fire, soils, fauna, and human activities. We identify 170 individual interactions among 32 elements that we present as a global tropical forest network, including countless feedback loops that may emerge from different combinations of interactions. We illustrate our findings with three cases involving urgent sustainability issues: (1) wildfires in wetlands of South America; (2) forest encroachment in African savanna landscapes; and (3) synergistic threats to the peatland forests of Borneo. Our findings reveal an unexplored world of feedbacks that shape the dynamics of tropical forests. The interactions and feedbacks identified here can guide future qualitative and quantitative research on the complexities of tropical forests, allowing societies to manage the nonlinear responses of these ecosystems in the Anthropocene.

Impact of Snowpack on the Land Surface Phenology in the Tianshan Mountains, Central Asia

The accumulation and ablation processes of seasonal snow significantly affect the land surface phenology in a mountainous ecosystem. However, the ability of snow to regulate the alpine land surface phenology in the arid regions is not well described in the context of climate change. The impact of snowpack changes on land surface phenology and its driving factors were investigated in the Tianshan Mountains using the land surface phenology metrics derived from satellited products and a snow dataset from downscaled regional climate model simulations covering the period from 1983 to 2015. The results demonstrated that the annual mean start of growing season (SOS) and length of growing season (LOS) experienced a significant (p < 0.05) decrease and increase with a rate of −2.45 days/decade and 2.98 days/decade, respectively. The significantly advanced SOS and increased LOS were mainly seen in the Western Tianshan Mountains and Ili Valley regions with elevations from 2500 to 3500 m a.s.l and below 3000 m a.s.l, respectively. During the early spring, the significant decline in snow cover fraction (SCF) could advance the SOS. In contrast, snowmelt amount and annual maximum snow water equivalent (SWE) have an almost equally substantial positive correlation with annual maximum vegetation greenness. In particular, the SOS of grassland was the most sensitive to variations of snow cover fraction during early spring than that of other vegetation types, and their strong relationship was mainly located at elevations from 1500 to 2500 m a.s.l. Its greenness was significantly controlled by the annual maximum snow water equivalent in all elevation bands. Both decreased SCF and increased temperature in the early spring caused a significant advance of the SOS, consequently prolonging the LOS. Meanwhile, more SWE and snowmelt amount could significantly promote vegetation greenness by regulating the soil moisture. The results can improve the understanding of the snow ecosystem services in the alpine regions under climate change.

Widespread increasing vegetation sensitivity to soil moisture

Global vegetation and associated ecosystem services critically depend on soil moisture availability which has decreased in many regions during the last three decades. While spatial patterns of vegetation sensitivity to global soil water have been recently investigated, long-term changes in vegetation sensitivity to soil water availability are still unclear. Here we assess global vegetation sensitivity to soil moisture during 1982-2017 by applying explainable machine learning with observation-based leaf area index (LAI) and hydro-climate anomaly data. We show that LAI sensitivity to soil moisture significantly increases in many semi-arid and arid regions. LAI sensitivity trends are associated with multiple hydro-climate and ecological variables, and strongest increasing trends occur in the most water-sensitive regions which additionally experience declining precipitation. State-of-the-art land surface models do not reproduce this increasing sensitivity as they misrepresent water-sensitive regions and sensitivity strength. Our sensitivity results imply an increasing ecosystem vulnerability to water availability which can lead to exacerbated reductions in vegetation carbon uptake under future intensified drought, consequently amplifying climate change.