Land surface temperature (LST) is a key variable in surface-atmosphere energy and water exchanges. The main goals of this study are to (i) evaluate the LST of the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim and ERA5 reanalyses over Iberian Peninsula using the Satellite Application Facility on Land Surface Analysis (LSA-SAF) product and to (ii) understand the main drivers of the LST errors in the reanalysis. Simulations with the ECMWF land-surface model in offline mode (uncoupled) were carried out over the Iberian Peninsula and compared with the reanalysis data. Several sensitivity simulations were performed in a confined domain centered in Southern Portugal to investigate potential sources of the LST errors. The Copernicus Global Land Service (CGLS) fraction of green vegetation cover (FCover) and the European Space Agency’s Climate Change Initiative (ESA-CCI) Land Cover dataset were explored. We found a general underestimation of daytime LST and slightly overestimation at night-time. The results indicate that there is still room for improvement in the simulation of LST in ECMWF products. Still, ERA5 presents an overall higher quality product in relation to ERA-Interim. Our analysis suggested a relation between the large daytime cold bias and vegetation cover differences between (ERA5 and CGLS FCocver) with a correlation of −0.45. The replacement of the low and high vegetation cover by those of ESA-CCI provided an overall reduction of the large Tmax biases during summer. The increased vertical resolution of the soil at the surface, has a positive impact, but much smaller when compared with the vegetation changes. The sensitivity of the vegetation density parameter, that currently depends on the vegetation type, provided further proof for a needed revision of the vegetation in the model, as there is a reasonable correlation between this parameter and the Tmax mean errors when using the ESA-CCI vegetation cover (while the same correlation cannot be reproduced with the original model vegetation). Our results support the hypothesis that vegetation cover is one of the main drivers of the LST summertime cold bias in ERA5 over Iberian Peninsula.

Deep learning (DL) methods have recently garnered attention from the climate change community for being an innovative approach to downscaling climate variables from Earth system and global climate models (ESGCMs) with horizontal resolutions still too coarse to represent regional- to local-scale phenomena. In the context of the Coupled Model Intercomparison Project phase 6 (CMIP6), ESGCM simulations were conducted for the Sixth Assessment Report (AR6) of the Intergovernmental Panel on Climate Change (IPCC) at resolutions ranging from 0.70 to 3.75^∘. Here, four convolutional neural network (CNN) architectures were evaluated for their ability to downscale, to a resolution of 0.1^∘, seven CMIP6 ESGCMs over the Iberian Peninsula – a known climate change hotspot, due to its increased vulnerability to projected future warming and drying conditions. The study is divided into three stages: (1) evaluating the performance of the four CNN architectures in predicting mean, minimum, and maximum temperatures, as well as daily precipitation, trained using ERA5 data and compared with the Iberia01 observational dataset; (2) downscaling the CMIP6 ESGCMs using the trained CNN architectures and further evaluating the ensemble against Iberia01; and (3) constructing a multi-model ensemble of CNN-based downscaled projections for temperature and precipitation over the Iberian Peninsula at 0.1^∘ resolution throughout the 21st century under four Shared Socioeconomic Pathway (SSP) scenarios. Upon validation and satisfactory performance evaluation, the DL downscaled projections demonstrate overall agreement with the CMIP6 ESGCM ensemble in magnitude for temperature projections and sign for the projected temperature and precipitation changes. Moreover, the advantages of using a high-resolution DL downscaled ensemble of ESGCM climate projections are evident, offering substantial added value in representing regional climate change over Iberia. Notably, a clear warming trend is observed in Iberia, consistent with previous studies in this area, with projected temperature increases ranging from 2 to 6 ^∘C, depending on the climate scenario. Regarding precipitation, robust projected decreases are observed in western and southwestern Iberia, particularly after 2040. These results may offer a new tool for providing regional climate change information for adaptation strategies based on CMIP6 ESGCMs prior to the next phase of the European branch of the Coordinated Regional Climate Downscaling Experiment (EURO-CORDEX) experiments.

Earth observations were used to evaluate the representation of land surface temperature (LST) and vegetation coverage over Iberia in two state-of-the-art land surface models (LSMs) – the European Centre for Medium-Range Weather Forecasts (ECMWF) Carbon-Hydrology Tiled ECMWF Scheme for Surface Exchanges over Land (CHTESSEL) and the Météo-France Interaction between Soil Biosphere and Atmosphere model (ISBA) within the SURface EXternalisée modeling platform (SURFEX-ISBA) for the 2004–2015 period. The results showed that the daily maximum LST simulated by CHTESSEL over Iberia was affected by a large cold bias during summer months when compared against the Satellite Application Facility on Land Surface Analysis (LSA-SAF), reaching magnitudes larger than 10 ^∘C over wide portions of central and southwestern Iberia. This error was shown to be tightly linked to a misrepresentation of the vegetation cover. In contrast, SURFEX simulations did not display such a cold bias. We show that this was due to the better representation of vegetation cover in SURFEX, which uses an updated land cover dataset (ECOCLIMAP-II) and an interactive vegetation evolution, representing seasonality. The representation of vegetation over Iberia in CHTESSEL was improved by combining information from the European Space Agency Climate Change Initiative (ESA-CCI) land cover dataset with the Copernicus Global Land Service (CGLS) leaf area index (LAI) and fraction of vegetation coverage (FCOVER). The proposed improvement in vegetation also included a clumping approach that introduces seasonality to the vegetation cover. The results showed significant added value, removing the daily maximum LST summer cold bias completely, without reducing the accuracy of the simulated LST, regardless of season or time of the day. The striking performance differences between SURFEX and CHTESSEL were fundamental to guiding the developments in CHTESSEL highlighting the importance of using different models. This work has important implications: first, it takes advantage of LST, a key variable in surface–atmosphere energy and water exchanges, which is closely related to satellite top-of-atmosphere observations, to improve the model's representation of land surface processes. Second, CHTESSEL is the land surface model employed by ECMWF in the production of their weather forecasts and reanalysis; hence systematic errors in land surface variables and fluxes are then propagated into those products. Indeed, we showed that the summer daily maximum LST cold bias over Iberia in CHTESSEL is present in the widely used ECMWF fifth-generation reanalysis (ERA5). Finally, our results provided hints about the interaction between vegetation land–atmosphere exchanges, highlighting the relevance of the vegetation cover and respective seasonality in representing land surface temperature in both CHTESSEL and SURFEX. As a whole, this work demonstrated the added value of using multiple earth observation products for constraining and improving weather and climate simulations.

Cities concentrate people, wealth, emissions, and infrastructure, thus representing a challenge and an opportunity for climate change mitigation and adaptation. This urgently demands for accurate urban climate projections to help organizations and individuals to make climate-smart decisions. However, most of the large ensembles of global and regional climate model simulations do not include sophisticated urban parameterizations (e.g., EURO-CORDEX; CMIP5/6). Here, we explore this shortcoming in ERA5 (the latest generation reanalysis from the European Centre for Medium-Range Weather Forecasts) and in a simulation with the SURFEX (Surface Externalisée) land surface model employing the widely used bulk bare rock approach. The city of Paris is considered as a case study. Subsequently, we apply a more complex urban scheme – SURFEX coupled to the Town Energy Balance (TEB) urban canopy model to assess its benefits on characterizing the Paris urban climate. Both simulations and ERA5 were compared to the LSA SAF (Satellite Application Facility on Land Surface Analysis) land surface temperature product to evaluate the simulation of Parisian surface urban heat island (SUHI). Our results show a significant added value of SURFEX-TEB in reproducing the SUHI during the daytime and the UHI during both the daytime and nighttime (with overall reductions in the bias and root mean square error and improvements in the representation of the statistics of the SUHI/UHI displayed by the Perkins skill score or S score). The improvement in the simulated SUHI is lower during the nighttime due to the lack of land–atmosphere feedbacks in the proposed offline framework. Nonetheless, the offline SURFEX-TEB framework applied here clearly demonstrates the added value of using more comprehensive parameterization schemes to simulate the urban climate and, therefore, allowing the improvement of urban climate projections.

Earth observations were used to evaluate the representation of Land Surface Temperature (LST) and vegetation coverage over Iberia in two state-of-the-art land surface models (LSMs) – the European Centre for Medium Range Weather Forecasting (ECMWF) Carbon-Hydrology Tiled ECMWF Scheme for Surface Exchanges over Land (CHTESSEL) and the Météo-France Interaction between Soil Biosphere and Atmosphere model (ISBA) within the SURface EXternalisée modelling platform (SURFEX-ISBA) for the 2004–2015 period. The results show that the daily maximum LST simulated by CHTESSEL over Iberia is affected by a large cold bias during summer months when compared against the Satellite Application Facility on Land Surface Analysis (LSA-SAF), reaching magnitudes larger than 10 °C over wide portions of central and southwestern Iberia. This error is shown to be tightly linked to a misrepresentation of the vegetation cover. In contrast, SURFEX simulations did not display such a cold bias. We show that this was due to the better representation of vegetation cover in SURFEX, which uses an updated land cover dataset (ECOCLIMAP-II) and an interactive vegetation evolution, representing seasonality. The representation of vegetation over Iberia in CHTESSEL was improved by combining information from the European Space Agency Climate Change Initiative (ESA-CCI) land cover dataset with the Copernicus Global Land Service (CGLS) Leaf Area Index (LAI) and fraction of vegetation coverage (FCOVER). The proposed improvement in vegetation also includes a clumping approach that introduces seasonality to the vegetation cover. The results show significant added value, removing the daily maximum LST summer cold bias completely, without reducing the accuracy of the simulated LST, regardless of season or time of the day. The striking performance differences between SURFEX and CHTESSEL were fundamental to guide the developments in CHTESSEL highlighting the importance of using different models. This work has important implications: first, it takes advantage of LST, a key variable in surface-atmosphere energy and water exchanges, which is closely related to satellite top-of-atmosphere observations, to improve model’s representation of land surface processes. Second, CHTESSEL is the land surface model employed by ECMWF in the production of their weather forecasts and reanalysis, hence systematic errors in land surface variables and fluxes are then propagated into those products. Indeed, we show that the summer daily maximum LST cold bias over Iberia in CHTESSEL is present in the widely used ECMWF fifth generation reanalysis (ERA5). Finally, our results provide hints into the interaction between vegetation land-atmosphere exchanges, highlighting the relevance of the vegetation cover and respective seasonality in representing land surface temperature in both CHTESSEL and SURFEX. As a whole, this work demonstrates the added value in using multiple earth observation products for constraining and improving weather and climate simulations.

Understanding and simulating urban climate processes, as well as how climate change affects cities is crucial for designing effective mitigation and adaptation strategies and policies. However, producing climate projections at the city scale requires very high-resolution physically-based models, which are computationally demanding and time-consuming to run. Deep Learning (DL) downscaling and offline simulations of land surface models offer cost- and time-effective alternatives.Here, we present a comparison between DL-based downscaling and offline simulations performed using the SURFEX land surface model for the city of Paris, France. Two lightweight 3-layer Convolutional Neural Network (CNN) architectures are trained to downscale ECMWF ERA5 reanalysis for the 2004-2012 period. The CNNs generate hourly predictions of 2-meter temperature (T2m) at point-level (using data from 24 in-situ observational stations) and Land Surface Temperature (LST) at a spatial resolution of ~5 km, respectively. The SURFEX land surface model (versions 8.1 and 9.0) is run at two different spatial resolutions (5 km and 1 km) for the 2013-2022 period. DL and SURFEX output are compared in terms of their representation of T2m, LST, and their respective urban heat island (UHI), surface urban heat island (SUHI) and extremes, in present climate (2013-2022). DL-based downscaling presents improved performance metrics in relation to SURFEX. DL also presents a diurnal cycle closer to the observations. Both DL-downscaling and SURFEX replicate the Parisian UHI effect, described in previous studies and in the observational data used to train the CNNs. This work supports the use of DL-based downscaling for urban climate studies as a viable alternative to more computationally and time-heavy approaches.Acknowledgements: This work is supported by FCT, I.P./MCTES through national funds (PIDDAC): LA/P/0068/2020 - https://doi.org/10.54499/LA/P/0068/2020, UID/50019/2025, https://doi.org/10.54499/UID/PRR/50019/2025, UID/PRR2/50019/2025. The authors would like also to acknowledge the project “Elaboração do Plano Municipal de Ação Climática de Barcelos" (PMACB).Frederico Johannsen was supported by FCT, I.P with the doctoral grant with the reference UI/BD/151498/2021 and DOI identifier 10.54499/UI/BD/151498/2021.