Getting the leaves right matters for estimating temperature extremes

Abstract. Atmospheric reanalyses combine observations and models through data assimilation techniques to provide spatio-temporally continuous fields of key surface variables. They can do so for extended historical periods whilst ensuring a coherent representation of the main Earth system cycles. ERA5 and its enhanced land surface component, ERA5-Land, are widely used in Earth system science and form the flagship products of the Copernicus Climate Change Service (C3S) of the European Commission. Such land surface modelling frameworks generally rely on a state variable called leaf area index (LAI), representing the number of leaves in a grid cell at a given time, to quantify the fluxes of carbon, water and energy between the vegetation and the atmosphere. However, the LAI within the modelling framework behind ERA5 and ERA5-Land is prescribed as a climatological seasonal cycle, neglecting any interannual variability and the potential consequences that this uncoupling between vegetation and atmosphere may have on the surface energy balance and the climate. To evaluate the impact of this mismatch in LAI, we analyse the corresponding effect it has on land surface temperature (LST) by comparing what is simulated to satellite observations. We characterise a hysteretic behaviour between LST biases and LAI biases that evolves differently along the year depending on the background climate. We further analyse the repercussions for the reconstructed climate during more extreme conditions in terms of LAI deviations, with a specific focus on the 2003, 2010 and 2018 heat waves in Europe for which LST mismatches are exacerbated. We anticipate that our results will assist users of ERA5 and ERA5-Land data in understanding where and when the larger discrepancies can be expected, but also guide developers towards improving the modelling framework. Finally, this study could provide a blueprint for a wider benchmarking framework for land surface model evaluation that exploits the capacity of LST to integrate the effects of both radiative and non-radiative processes affecting the surface energy.