Widespread decline in critical soil moisture threshold of plant water stress across global drylands

Python version 3.11, R version 4.4.1, and Matlab2024a were used.
1. drydown_identification.m was executed to identify soil dry-down events and to save the corresponding datasets. Then, SM_threshold_estimation.R was run to estimate soil moisture thresholds.
2. Figure1.py was used to perform trend analysis of soil moisture thresholds and to generate Figure 1. Figure2.R was executed to generate Figure 2. Figure3.R was used to conduct machine-learning–based and multivariate linear regression attribution analyses and to generate Figure 3. Finally, Figure4.py and Figure5.py were executed to generate Figures 4 and 5, respectively.