Empirical vegetation indices derived from spectral reflectance data are widely used in remote sensing of the biosphere, as they represent robust proxies for canopy structure, leaf pigment content, and, subsequently, plant photosynthetic potential. Here, we generalize the broad family of commonly used vegetation indices by exploiting all higher-order relations between the spectral channels involved. This results in a higher sensitivity to vegetation biophysical and physiological parameters. The presented nonlinear generalization of the celebrated normalized difference vegetation index (NDVI) consistently improves accuracy in monitoring key parameters, such as leaf area index, gross primary productivity, and sun-induced chlorophyll fluorescence. Results suggest that the statistical approach maximally exploits the spectral information and addresses long-standing problems in satellite Earth Observation of the terrestrial biosphere. The nonlinear NDVI will allow more accurate measures of terrestrial carbon source/sink dynamics and potentials for stabilizing atmospheric CO₂ and mitigating global climate change.

Changing vegetation cover alters the radiative and non-radiative properties of the surface. The result of competing biophysical processes on Earth’s surface energy balance varies spatially and seasonally, and can lead to warming or cooling depending on the specific vegetation change and background climate. Here we provide the first data-driven assessment of the potential effect on the full surface energy balance of multiple vegetation transitions at global scale. For this purpose we developed a novel methodology that is optimized to disentangle the effect of mixed vegetation cover on the surface climate. We show that perturbations in the surface energy balance generated by vegetation change from 2000 to 2015 have led to an average increase of 0.23 ± 0.03 °C in local surface temperature where those vegetation changes occurred. Vegetation transitions behind this warming effect mainly relate to agricultural expansion in the tropics, where surface brightening and consequent reduction of net radiation does not counter-balance the increase in temperature associated with reduction in transpiration. This assessment will help the evaluation of land-based climate change mitigation plans.

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.

We estimate global terrestrial gross primary production (GPP) based on models that use satellite data within a simplified light-use efficiency framework that does not rely upon other meteorological inputs. Satellite-based geometry-adjusted reflectances are from the MODerate-resolution Imaging Spectroradiometer (MODIS) and provide information about vegetation structure and chlorophyll content at both high temporal (daily to monthly) and spatial (∼1 km) resolution. We use satellite-derived solar-induced fluorescence (SIF) to identify regions of high productivity crops and also evaluate the use of downscaled SIF to estimate GPP. We calibrate a set of our satellite-based models with GPP estimates from a subset of distributed eddy covariance flux towers (FLUXNET 2015). The results of the trained models are evaluated using an independent subset of FLUXNET 2015 GPP data. We show that variations in light-use efficiency (LUE) with incident PAR are important and can be easily incorporated into the models. Unlike many LUE-based models, our satellite-based GPP estimates do not use an explicit parameterization of LUE that reduces its value from the potential maximum under limiting conditions such as temperature and water stress. Even without the parameterized downward regulation, our simplified models are shown to perform as well as or better than state-of-the-art satellite data-driven products that incorporate such parameterizations. A significant fraction of both spatial and temporal variability in GPP across plant functional types can be accounted for using our satellite-based models. Our results provide an annual GPP value of ∼140 Pg C year ${}^{-1}$ for 2007 that is within the range of a compilation of observation-based, model, and hybrid results, but is higher than some previous satellite observation-based estimates.