Abstract. In this work, an ensemble of machine learning algorithms was trained using stratified sampling from an existing European-scale biomass map from 2018 to predict an updated version for 2020. The objective of stratification is to make sure that the full range of biomass values is represented. The sampled biomass values from 2018 were filtered to remove areas that did were subject to forest disturbances between 2018 and 2020. This information was available from forest cover/loss/gain maps derived from satellite imagery. We train using a total of 49 features derived from the following sources: bioclimatic data, maps of land-cover, tree cover, tree height, annual composites of vegetation indices per pixel (EVI and NDVI) obtained from Sentinel-2, radar backscatter median annual values from Sentinel-1 and ALOS-2, and the ALOS DSM (3D) elevation grid. A model was created dividing Europe into 19 tiles to limit variability due to very different bioclimatic zones. The result is a raster with 100 m × 100 m resolution and an estimated value of biomass (Mg ha−1) at each node. Overall results on validation data over Europe report a root mean square error (RMSE) of 32.4 Mg ha−1 and a mean absolute error (MAE) of 21.5 Mg ha−1; when considering single tiles, the largest RMSE was 54.7 Mg ha−1 in tile D2, which can be explained by the very high variance of climate, environment, terrain topography and biomass values as the tile enclosed the Alpine region and the western part of Eastern Europe.