The ratio of plant transpiration to total terrestrial evapotranspiration (T/ET) captures the role of vegetation in surface-atmosphere interactions. However, several studies have documented a large variability in T/ET. In this paper, we present a new T/ET dataset (also including transpiration, evapotranspiration data) for China from 1981 to 2015 with spatial and temporal resolutions of 0.05° and 8 days, respectively. The T/ET dataset is based on a model-data fusion method that integrates the Priestley-Taylor Jet Propulsion Laboratory (PT-JPL) model with multivariate observational datasets (transpiration and evapotranspiration). The dataset is driven by satellite-based leaf area index (LAI) data from GLASS and GLOBMAP, and climate data from the Chinese Ecosystem Research Network (CERN). Observational annual T/ET were used to validate the model, withR²and RMSE values were 0.73 and 0.07 (12.41%), respectively. The dataset provides significant insight into T/ET and its changes over the Chinese terrestrial ecosystem and will be beneficial for understanding the hydrological cycle and energy budgets between the land and the atmosphere.

Terrestrial ecosystems provide multiple services interacting in complex ways. However, most ecosystem services (ESs) models (e.g., InVEST and ARIES) ignored the relationships among ESs. Process‐based models can overcome this limitation, and the integration of ecological models with remote sensing data could greatly facilitate the investigation of the complex ecological processes. Therefore, based on the Carbon and Exchange between Vegetation, Soil, and Atmosphere (CEVSA) models, we developed a process‐based ES model (CEVSA‐ES) integrating remotely sensed leaf area index to evaluate four important ESs (i.e., productivity provision, carbon sequestration, water retention, and soil retention) at annual timescale in China. Compared to the traditional terrestrial biosphere models, the main innovation of CEVSA‐ES model was the consideration of soil erosion processes and its impact on carbon cycling. The new version also improved the carbon‐water cycle algorithms. Then, the Sobol and DEMC methods that integrated the CEVSA‐ES model with nine flux sites comprising 39 site‐years were used to identify and optimize parameters. Finally, the model using the optimized parameters was validated at 26 field sites comprising 135 site‐years. Simulation results showed good fits with ecosystem processes, explaining 95%, 92%, 76%, and 65% interannual variabilities of gross primary productivity, ecosystem respiration, net ecosystem productivity, and evapotranspiration, respectively. The CEVSA‐ES model performed well for productivity provision and carbon sequestration, which explained 96% and 81% of the spatial‐temporal variations of the observed annual productivity provision and carbon sequestration, respectively. The model also captured the interannual trends of water retention and soil erosion for most sites or basins.