Codes for "Latitudinal divergence in runoff responses to global forestation due to forest-atmosphere feedbacks"

Attached please find the analysis codes of "Latitudinal divergence in runoff responses to global forestation due to forest-atmosphere feedbacks".Above all analysis codes are based on Matlab R2023a, including:(1) PET_Budyko_n_Estimation.m: Code for calculating PET via multiple methods and inverting Budyko parameter n.(2) dR_Budyko_decomposition.m: Code for decomposing runoff changes into precipitation, PET and n change contributions.(3) MoistureBudgetDecomposition.m: Code for atmospheric moisture budget decomposition of precipitation changes.(4) Readme.txt: Brief description of the three analysis codes.If you use this code in your research, please cite our paper:Kan et al. Latitudinal divergence in runoff responses to global forestation due to forest-atmosphere feedbacks Nature Communications (2026). https://doi.org/10.1038/s41467-026-68945-9