Vegetation plays a fundamental role in modulating the exchange of water, energy, and carbon fluxes between the land and the atmosphere. These exchanges are modelled by Land Surface Models (LSMs), which are an essential part of numerical weather prediction and data assimilation. However, most current LSMs implemented specifically in weather forecasting systems use climatological vegetation indices, and land use/land cover datasets in these models are often outdated. In this study, we update land surface data in the ECMWF land surface modelling system ECLand using Earth observation-based time varying leaf area index and land use/land cover data, and evaluate the impact of vegetation dynamics on model performance. The performance of the simulated latent heat flux and soil moisture is then evaluated against global gridded observation-based datasets. Updating the vegetation information does not always yield better model performances because the model’s parameters are adapted to the previously employed land surface information. Therefore we recalibrate key soil and vegetation-related parameters at individual grid cells to adjust the model parameterizations to the new land surface information. This substantially improves model performance and demonstrates the benefits of updated vegetation information. Interestingly, we find that a regional parameter calibration outperforms a globally uniform adjustment of parameters, indicating that parameters should sufficiently reflect spatial variability in the land surface. Our results highlight that newly available Earth-observation products of vegetation dynamics and land cover changes can improve land surface model performances, which in turn can contribute to more accurate weather forecasts.

A long-standing challenge in studying the global carbon cycle has been understanding the factors controlling inter–annual variation (IAV) of carbon fluxes related to vegetation photosynthesis and respiration, and improving their representations in existing biogeochemical models. Here, we compared an optimality-based mechanistic model and a semi-empirical light use efficiency model to understand how current models can be improved to simulate IAV of gross primary production (GPP). Both models simulated hourly GPP and were parameterized for (1) each site–year, (2) each site with an additional constraint on IAV (CostIAV), (3) each site, (4) each plant–functional type, and (5) globally. This was followed by forward runs using calibrated parameters, and model evaluations at different temporal scales across 198 eddy covariance sites. Both models performed better on hourly scale than annual scale for most sites. Specifically, the mechanistic model substantially improved when drought stress was explicitly included. Most of the variability in model performances was due to model types and parameterization strategies. The semi-empirical model produced statistically better hourly simulations than the mechanistic model, and site–year parameterization yielded better annual performance for both models. Annual model performance did not improve even when parameterized using CostIAV. Furthermore, both models underestimated the peaks of diurnal GPP in each site–year, suggesting that improving predictions of peaks could produce a comparatively better annual model performance. GPP of forests were better simulated than grassland or savanna sites by both models. Our findings reveal current model deficiencies in representing IAV of carbon fluxes and guide improvements in further model development.

The spatial contribution to the global land-atmosphere carbon dioxide (CO\textsubscript{2}) exchange is crucial in understanding and projecting the global carbon cycle, yet different studies diverge on the dominant regions. Informing land models with observational data is a promising way to reduce the parameter and structural uncertainties and advance our understanding. Here, we develop a parsimonious diagnostic process-based model of land carbon cycles, constraining parameters with observation-based products. We compare CO\textsubscript{2} flux estimates from our model with observational constraints and Trends in Net Land-Atmosphere Carbon Exchange (TRENDY) model ensemble to show that our model reasonably reproduces the seasonality of net ecosystem exchange (NEE) and GPP and interannual variability (IAV) of NEE. Finally, we use the developed model, TRENDY models, and observational constraints to attribute variability in global NEE and gross primary productivity (GPP) to regional variability. The attribution analysis confirms the dominance of Northern temperate and boreal regions in the seasonality of CO\textsubscript{2} fluxes. Regarding NEE IAV, we identify a significant contribution from tropical savanna regions as previously perceived. Furthermore, we highlight that tropical humid regions are also identified as at least equally relevant contributors as semi-arid regions. At the same time, the largest uncertainty among ensemble members of NEE constraint and TRENDY models in the tropical humid regions underscore the necessity of better process understanding and more observations in these regions. Overall, our study identifies tropical humid regions as key regions for global land-atmosphere CO\textsubscript{2} exchanges and the inter-model spread of its modeling.

Natural and anthropogenic disturbances act as important drivers of tree mortality, shaping the structure, composition and biomass distribution of forests. Disturbance regimes may emerge from different characteristics of disturbance events over time and space. We design a model- based experiment to investigate the links between disturbance regimes at the landscape scale and spatial features of biomass patterns. The effects on biomass of a wide range of disturbance regimes are simulated by varying three different parameters, i.e. μ (probability scale), α (clustering degree), and β (intensity slope) that shape the extent, frequency, and intensity of disturbance events, respectively. A simple dynamic carbon cycle model is used to simulate 200 years of plant biomass dynamics in response to circa +2000 different disturbance regimes, depending on the different combinations of μ, α, and β. Each parameter combination yields a spatially explicit estimate of plant biomass for which sixteen synthesis statistics are estimated on the spatial distributions of biomass, including information-based and texture features. Based on a multi-output regression approach we link these synthesis statistics with additional gross primary production (GPP) constraints to retrieve the three disturbance parameters. In doing so we evaluate the confidence in inferring disturbance regimes from spatial distributions of biomass. Our results show that all three parameters can be confidently retrieved. The Nash-Sutcliffe efficiency for the prediction of the μ, α, and β is 97.3%, 96.6%, and 97.9%, respectively. A feature importance analysis reveals that the distribution statistics dominate the prediction of μ and β, while features quantifying texture have a stronger connection with α. Overall, this study clarifies the association between biomass patterns emerging from different underlying disturbance regimes, while overcoming the previously found equifinality between mortality rates and total biomass. Given the links between decadal vegetation dynamics and the uncertainties in the role of terrestrial ecosystems in the global biogeochemical cycles, a better understanding and the quantification of disturbance regimes would improve our current understanding of controls and feedback at the biosphere-atmosphere interface in the current Earth system models.

Vegetation turnover time (τ) is a central ecosystem property to quantify the global vegetation carbon dynamics. However, our understanding of vegetation dynamics is hampered by the lack of long-term observations of the changes in vegetation biomass. Here we challenge the steady state assumption of τ by using annual changes in vegetation biomass that derived from remote-sensing observations. We evaluate the changes in magnitude, spatial patterns, and uncertainties in vegetation carbon turnover times from 1992 to 2016. We found that the forest ecosystem is close to a steady state at global scale, contrasting with the larger differences between τ under steady state and τ under non-steady state at the grid cell level. The observation that terrestrial ecosystems are not in a steady state locally is deemed crucial when studying vegetation dynamics and the potential response of biomass to disturbance and climatic changes.

Vegetation turnover time (τ) is a central ecosystem property to quantify the global vegetation carbon dynamics. However, our understanding of vegetation dynamics is hampered by the lack of long-term observations of the changes in vegetation biomass. Here we challenge the steady state assumption of τ by using annual changes in vegetation biomass that derived from remote-sensing observations. We evaluate the changes in magnitude, spatial patterns, and uncertainties in vegetation carbon turnover times from 1992 to 2016. We found that the forest ecosystem is close to a steady state at global scale, contrasting with the larger differences between τ under steady state and τ under non-steady state at the grid cell level. The observation that terrestrial ecosystems are not in a steady state locally is deemed crucial when studying vegetation dynamics and the potential response of biomass to disturbance and climatic changes.

Hydrological interactions between vegetation, soil, and topography are complex, and heterogeneous in semi-arid landscapes. This along with data scarcity poses challenges for large-scale modelling of vegetation-water interactions. Here, we exploit metrics derived from daily Meteosat data over Africa at ca. 5 km spatial resolution for ecohydrological analysis. Their spatial patterns are based on Fractional Vegetation Cover (FVC) time series and emphasise limiting conditions of the seasonal wet to dry transition: the minimum and maximum FVC of temporal record, the FVC decay rate and the FVC integral over the decay period. We investigate the relevance of these metrics for large scale ecohydrological studies by assessing their co-variation with soil moisture, and with topographic, soil, and vegetation factors. Consistent with our initial hypothesis, FVC minimum and maximum increase with soil moisture, while the FVC integral and decay rate peak at intermediate soil moisture. We find evidence for the relevance of topographic moisture variations in arid regions, which, counter-intuitively, is detectable in the maximum but not in the minimum FVC. We find no clear evidence for wide-spread occurrence of the “inverse texture effect”’ on FVC. The FVC integral over the decay period correlates with independent data sets of plant water storage capacity or rooting depth while correlations increase with aridity. In arid regions, the FVC decay rate decreases with canopy height and tree cover fraction as expected for ecosystems with a more conservative water-use strategy. Thus, our observation-based products have large potential for better understanding complex vegetation–water interactions from regional to continental scales.

Next to precipitation, secondary water sources emerging from shallow groundwater and lateral redistribution of soil moisture, together with soil properties modulating their accessibility are highly important in water-limited ecosystems. However, effects of these land-associated secondary inputs are not well known over large domains given the mismatch of spatial scales of processes. Here, we quantify the role of land properties on the spatial variations of seasonal decay rate of vegetation cover over water-limited regions of Africa, using machine learning. Over the study domain, 17 % of these variations are directly attributed to land properties, and 16 % are attributed to interaction effects of land properties with climate and vegetation. Locally, total land attributed variations account for more than 60 % in hotspots with different land properties like shallow groundwater, complex topography, and favourable soil properties. Our findings lend empirical evidence for the importance of local-scale secondary water inputs over large domains.