l Evapotranspiration (ET) was also estimated based on multiple datasets, including (1) the Global Land Evaporation (GLEAM v3.6a), which was a combination of gauge-based, reanalysis, and satellite-based data, spanning the 42 years 1980-2021 at a spatial resolution of 0.25°. (2) PML_V2 dataset, coupling with the GPP process, at 8-day intervals from 2002 to 2019 and a spatial resolution of 0.05°. (3) Satellite-estimated ET from Song et al., which gave monthly ET from 2001 to 2020 at 0.01° spatial resolution.