Regional model ICON-D2 The DWD's ICON-D2 model is a forecast model which is operated for the very-short range up to +27 hours (+45 hours for the 03 UTC run). Due to its fine mesh size, the ICON-D2 especially provides for improved forecasts of hazardous weather conditions, e.g. weather situations with high-level moisture convection (super and multi-cell thunderstorms, squall lines, mesoscale convective complexes) and weather events that are influenced by fine-scale topographic effects (ground fog, Föhn winds, intense downslope winds, flash floods). The model area of ICON-D2 covers the whole German territory, Benelux, Switzerland, Austria and parts of the other neighbouring countries at a horizontal resolution of 2.2 km. In the vertical, the model defines 65 atmosphere levels. The fairly short forecast periods make perfect sense because of the purpose of ICON-D2 (and its small model area). Based on model runs at 00, 06, 09, 12, 15, 18 and 21 UTC, ICON-D2 provides new 27-hour forecasts every 3 hours. The model run at 03 UTC even covers a forecast period of 45 hours. The ICON-D2 forecast data for each weather element are made available in standard packages at our free DWD Open Data Server, both on a rotated grid and on a regular grid. Regional ensemble forecast model ICON-D2 EPS The ensemble forecasting system ICON-D2 EPS is based on the DWD's numerical weather forecast model ICON-D2 and currently includes 20 ensemble members. All ensemble members are calculated at the same horizontal grid spacing as the operational configuration of ICON-D2 (2.2 km). Like ICON-D2, the ICON-D2 EPS ensemble system provides forecasts up to +27 hours for the same model area (up to +45 hours based on the 03 UTC run). For generating the ensemble members, some of the features of the forecasting system are changed. The method currently used to generate the ensemble members involves varying the - lateral boundary conditions - initial state - soil moisture - and model physics. For varying the lateral boundary conditions and the initial state, forecasts from various global models are used. The ICON-D2 EPS is provided on the DWD Open Data Server in the native triangular grid. Note: All previously COSMO-D2 based aviation weather products have been migrated to ICON-D2 on 10.02.2021. However, the familiar design of these products remains unchanged.

The Copernicus DEM is a Digital Surface Model (DSM) which represents the surface of the Earth including buildings, infrastructure and vegetation. The original GLO-30 provides worldwide coverage at 30 meters (refers to 10 arc seconds). Note that ocean areas do not have tiles, there one can assume height values equal to zero. Data is provided as Cloud Optimized GeoTIFFs. Note that the vertical unit for measurement of elevation height is meters. The Copernicus DEM for Europe at 1000 meter resolution (EU-LAEA projection) in COG format has been derived from the Copernicus DEM GLO-30, mirrored on Open Data on AWS, dataset managed by Sinergise (https://registry.opendata.aws/copernicus-dem/). Processing steps: The original Copernicus GLO-30 DEM contains a relevant percentage of tiles with non-square pixels. We created a mosaic map in https://gdal.org/drivers/raster/vrt.html format and defined within the VRT file the rule to apply cubic resampling while reading the data, i.e. importing them into GRASS GIS for further processing. We chose cubic instead of bilinear resampling since the height-width ratio of non-square pixels is up to 1:5. Hence, artefacts between adjacent tiles in rugged terrain could be minimized: gdalbuildvrt -input_file_list list_geotiffs_MOOD.csv -r cubic -tr 0.000277777777777778 0.000277777777777778 Copernicus_DSM_30m_MOOD.vrt In order to reproject the data to EU-LAEA projection while reducing the spatial resolution to 1000 m, bilinear resampling was performed in GRASS GIS (using r.proj) and the pixel values were scaled with 1000 (storing the pixels as Integer values) for data volume reduction. In addition, a hillshade raster map was derived from the resampled elevation map (using r.relief, GRASS GIS). Eventually, we exported the elevation and hillshade raster maps in Cloud Optimized GeoTIFF (COG) format, along with SLD and QML style files.

osm: Buildings according to OSM polygons and the Copernicus impervious build-up layer (2018), aggregated and rasterized first to 10m spatial resolution and after downsampled to 30m by spatial average.

Overview: Actual Natural Vegetation (ANV): probability of occurrence for the Silver fir in its realized environment for the period 2000 - 2020 Traceability (lineage): This is an original dataset produced with a machine learning framework which used a combination of point datasets and raster datasets as inputs. Point dataset is a harmonized collection of tree occurrence data, comprising observations from National Forest Inventories (EU-Forest), GBIF and LUCAS. The complete dataset is available on Zenodo. Raster datasets used as input are: harmonized and gapfilled time series of seasonal aggregates of the Landsat GLAD ARD dataset (bands and spectral indices); monthly time series air and surface temperature and precipitation from a reprocessed version of the Copernicus ERA5 dataset; long term averages of bioclimatic variables from CHELSA, tree species distribution maps from the European Atlas of Forest Tree Species; elevation, slope and other elevation-derived metrics; long term monthly averages snow probability and long term monthly averages of cloud fraction from MODIS. For a more comprehensive list refer to Bonannella et al. (2022) (in review, preprint available at: https://doi.org/10.21203/rs.3.rs-1252972/v1). Scientific methodology: Probability and uncertainty maps were the output of a spatiotemporal ensemble machine learning framework based on stacked regularization. Three base models (random forest, gradient boosted trees and generalized linear models) were first trained on the input dataset and their predictions were used to train an additional model (logistic regression) which provided the final predictions. More details on the whole workflow are available in the listed publication. Usability: Probability maps can be used to detect potential forest degradation and compositional change across the time period analyzed. Some possible applications for these topics are explained in the listed publication. Uncertainty quantification: Uncertainty is quantified by taking the standard deviation of the probabilities predicted by the three components of the spatiotemporal ensemble model. Data validation approaches: Distribution maps were validated using a spatial 5-fold cross validation following the workflow detailed in the listed publication. Completeness: The raster files perfectly cover the entire Geo-harmonizer region as defined by the landmask raster dataset available here. Consistency: Areas which are outside of the calibration area of the point dataset (Iceland, Norway) usually have high uncertainty values. This is not only a problem of extrapolation but also of poor representation in the feature space available to the model of the conditions that are present in this countries. Positional accuracy: The rasters have a spatial resolution of 30m. Temporal accuracy: The maps cover the period 2000 - 2020, each map covers a certain number of years according to the following scheme: (1) 2000--2002, (2) 2002--2006, (3) 2006--2010, (4) 2010--2014, (5) 2014--2018 and (6) 2018--2020 Thematic accuracy: Both probability and uncertainty maps contain values from 0 to 100: in the case of probability maps, they indicate the probability of occurrence of a single individual of the target species, while uncertainty maps indicate the standard deviation of the ensemble model.

The Land Cover Map of Europe 2017 is a product resulting from the Phase 2 of the S2GLC project. The final map has been produced on the CREODIAS platform with algorithms and software developed by CBK PAN. Classification of over 15 000 Sentinel-2 images required high level of automation that was assured by the developed software. The legend of the resulting Land Cover Map of Europe 2017 consists of 13 land cover classes. The pixel size of the map equals 10 m, which corresponds to the highest spatial resolution of Sentinel-2 imagery. Its overall accuracy was estimated to be at the level of 86% using approximately 52 000 validation samples distributed across Europe. Related publication: https://doi.org/10.3390/rs12213523

222: Cultivated parcels planted with fruit trees and shrubs, intended for fruit production, including nuts. The planting pattern can be by single or mixed fruit species, both in association with permanently grassy surfaces.

Overview: osm: Harbours rasterized from OSM landuse polygons, first to 10m spatial resolution and after downsampled to 30m by spatial average. Traceability (lineage): The class-wise layers of this dataset were extracted from OpenStreetMap data downloaded from geofabrik.de and aggregated based on labels assigned to the volunteered geographical information objects. Scientific methodology: nan Usability: The extracted classes can be used to preprocess training data (as detailed in Witjes et al., 2022 (in review, preprint available at https://doi.org/10.21203/rs.3.rs-561383/v3 ). Users are advised to remember the potential inconsistencies in volunteered geographical information, however: Some regions of Europe have been less consistently mapped in OpenStreetMap. This may introduce bias in any subsequent modelling. Uncertainty quantification: nan Data validation approaches: This dataset has not been validated Completeness: Volunteered geographical information often more complete in regions with more active contributors. It is likely that this dataset contains many omission errors in regions of Europe where OpenStreetMap is used less intensively. Consistency: Volunteered geographical information often more complete in regions with more active contributors. It is likely that this dataset contains many omission errors in regions of Europe where OpenStreetMap is used less intensively. Positional accuracy: The rasters have a spatial resolution of 30m Temporal accuracy: The maps are based on an extract from 2020. Thematic accuracy: The 30m pixels of each OSM extract map have values ranging from 0-100, indicating the density aggregated from 10m pixels where rasterized objects burned the value 100 in a 0-value raster.

Overview: Potential Natural Vegetation (PNV): potential probability of occurrence for the Olive tree from 2018 to 2020 Traceability (lineage): This is an original dataset produced with a machine learning framework which used a combination of point datasets and raster datasets as inputs. Point dataset is a harmonized collection of tree occurrence data, comprising observations from National Forest Inventories (EU-Forest), GBIF and LUCAS. The complete dataset is available on Zenodo. Raster datasets used as input are: monthly time series air and surface temperature and precipitation from a reprocessed version of the Copernicus ERA5 dataset; long term averages of bioclimatic variables from CHELSA; elevation, slope and other elevation-derived metrics and long term monthly averages snow probability. For a more comprehensive list refer to Bonannella et al. (2022) (in review, preprint available at: https://doi.org/10.21203/rs.3.rs-1252972/v1). Scientific methodology: Probability and uncertainty maps were the output of a spatiotemporal ensemble machine learning framework based on stacked regularization. Three base models (random forest, gradient boosted trees and generalized linear models) were first trained on the input dataset and their predictions were used to train an additional model (logistic regression) which provided the final predictions. More details on the whole workflow are available in the listed publication. Usability: Probability maps are particularly useful when compared with existing products of potential distribution of species or when combined with maps of realized distribution: gaps in potential and realized distribution can be identified and used as information for future programs of tree planting or forest restoration. Uncertainty quantification: Uncertainty is quantified by taking the standard deviation of the probabilities predicted by the three components of the spatiotemporal ensemble model. Data validation approaches: Distribution maps were validated using a spatial 5-fold cross validation following the workflow detailed in the listed publication. Completeness: The raster files perfectly cover the entire Geo-harmonizer region as defined by the landmask raster dataset available here. Consistency: Areas which are outside of the calibration area of the point dataset (Iceland, Norway) usually have high uncertainty values. This is not only a problem of extrapolation but also of poor representation in the feature space available to the model of the conditions that are present in this countries. Positional accuracy: The rasters have a spatial resolution of 30m. Temporal accuracy: The maps cover the period 2018 - 2020 Thematic accuracy: Both probability and uncertainty maps contain values from 0 to 100: in the case of probability maps, they indicate the probability of occurrence of a single individual of the target species, while uncertainty maps indicate the standard deviation of the ensemble model.

ERA5-Land total precipitation monthly time series for Mauritania at 30 arc seconds (ca. 1000 meter) resolution (2019 - 2023) Source data: ERA5-Land is a reanalysis dataset providing a consistent view of the evolution of land variables over several decades at an enhanced resolution compared to ERA5. ERA5-Land has been produced by replaying the land component of the ECMWF ERA5 climate reanalysis. Reanalysis combines model data with observations from across the world into a globally complete and consistent dataset using the laws of physics. Reanalysis produces data that goes several decades back in time, providing an accurate description of the climate of the past. Total precipitation: Accumulated liquid and frozen water, including rain and snow, that falls to the Earth's surface. It is the sum of large-scale precipitation (that precipitation which is generated by large-scale weather patterns, such as troughs and cold fronts) and convective precipitation (generated by convection which occurs when air at lower levels in the atmosphere is warmer and less dense than the air above, so it rises). Precipitation variables do not include fog, dew or the precipitation that evaporates in the atmosphere before it lands at the surface of the Earth. This variable is accumulated from the beginning of the forecast time to the end of the forecast step. The units of precipitation are depth in metres. It is the depth the water would have if it were spread evenly over the grid box. Care should be taken when comparing model variables with observations, because observations are often local to a particular point in space and time, rather than representing averages over a model grid box and model time step. Processing steps: The original hourly ERA5-Land data has been spatially enhanced from 0.1 degree to 30 arc seconds (approx. 1000 m) spatial resolution by image fusion with CHELSA data (V1.2) (https://chelsa-climate.org/). For each day we used the corresponding monthly long-term average of CHELSA. The aim was to use the fine spatial detail of CHELSA and at the same time preserve the general regional pattern and fine temporal detail of ERA5-Land. The steps included aggregation and enhancement, specifically: 1. spatially aggregate CHELSA to the resolution of ERA5-Land 2. calculate proportion of ERA5-Land / aggregated CHELSA 3. interpolate proportion with a Gaussian filter to 30 arc seconds 4. multiply the interpolated proportions with CHELSA Using proportions ensures that areas without precipitation remain areas without precipitation. Only if there was actual precipitation in a given area, precipitation was redistributed according to the spatial detail of CHELSA. The spatially enhanced daily ERA5-Land data has been aggregated to monthly resolution, by calculating the sum of the precipitation per pixel over each month. File naming: ERA5_land_monthly_prectot_sum_30sec_YYYY_MM_01T00_00_00_int.tif e.g.:ERA5_land_monthly_prectot_sum_30sec_2023_12_01T00_00_00_int.tif The date within the filename is year and month of aggregated timestamp. Pixel values: mm * 10 Scaled to Integer, example: value 218 = 21.8 mm Projection + EPSG code: Latitude-Longitude/WGS84 (EPSG: 4326) Spatial extent: north: 28:18N south: 14:42N west: 17:05W east: 4:49W Temporal extent: January 2019 - December 2023 Spatial resolution: 30 arc seconds (approx. 1000 m) Temporal resolution: monthly Lineage: Dataset has been processed from original Copernicus Climate Data Store (ERA5-Land) data sources. As auxiliary data CHELSA climate data has been used. Software used: GRASS GIS 8.3.2 Format: GeoTIFF Original ERA5-Land dataset license: https://cds.climate.copernicus.eu/api/v2/terms/static/licence-to-use-copernicus-products.pdf CHELSA climatologies (V1.2): Data used: Karger D.N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R.W., Zimmermann, N.E, Linder, H.P., Kessler, M. (2018): Data from: Climatologies at high resolution for the earth's land surface areas. Dryad digital repository. http://dx.doi.org/doi:10.5061/dryad.kd1d4 Original peer-reviewed publication: Karger, D.N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R.W., Zimmermann, N.E., Linder, P., Kessler, M. (2017): Climatologies at high resolution for the Earth land surface areas. Scientific Data. 4 170122. https://doi.org/10.1038/sdata.2017.122 Representation type: Grid Processed by: mundialis GmbH & Co. KG, Germany (https://www.mundialis.de/) Contact: mundialis GmbH & Co. KG, info@mundialis.de Acknowledgements: This study was partially funded by EU grant 874850 MOOD. The contents of this publication are the sole responsibility of the authors and don't necessarily reflect the views of the European Commission.

Overview: Aerosol Optical Depth (AOD) at 550 nm Traceability (lineage): The data source used to generate this dataset are daily Aerosol Optical Depth (AOD) at 550 nm with 1 km spatial resolution from the MCD19A2 product of MODIS collection 6 for the years 2018-2020 (https://lpdaac.usgs.gov/products/mcd19a2v006/), daily modeled AOD at 550 nm with 80 km spatial resolution from Copernicus Atmosphere Monitoring Service (CAMS) for the years 2018-2020, and elevations from the worldwide digital surface model of the Japan Aerospace Exploration Agency (JAXA) (https://www.eorc.jaxa.jp/ALOS/). NASA’s Aerosols Robotic Network (AERONET) ground measurements were used to validate this dataset. Scientific methodology: The outcome of the scheme was a geo-harmonized atmospheric dataset for aerosol optical depth (GHADA) at 550 nm with 1 km spatial resolution and full coverage over Europe. Daily AOD maps were created by training an optimized space-time extra trees model for each year 2018-2020. The results have shown that our trained models reach a prediction accuracy up to 95% when predicting the missing values in the MODIS MCD19A2 product. Usability: AOD maps can be used for future air quality studies concerning Europe. Uncertainty quantification: nan Data validation approaches: GHADA was validated using AOD measurements from AERONET stations across Europe. Completeness: The raster files cover the entire Geo-harmonizer region. Consistency: nan Positional accuracy: 1 km spatial resolution Temporal accuracy: Daily maps for 3 years, 2018-2020. Thematic accuracy: nan