Paper describes a data-driven framework based on spatio-temporal ensemble machine learning to produce distribution maps for 16 forest tree species (Abies alba Mill., Castanea sativa Mill. , Corylus avellana L., Fagus sylvatica L., Olea europaea L., Picea abies L. H. Karst., Pinus halepensis Mill., Pinus nigra J. F. Arnold, Pinus pinea L., Pinus sylvestris L., Prunus avium L., Quercus cerris L., Quercus ilex L., Quercus robur L., Quercus suber L. and Salix caprea L.) at high spatial resolution (30 m). Tree occurrence data for a total of 3 million of points was used to train different Machine Learning (ML) algorithms: random forest, gradient-boosted trees, generalized linear models, k-nearest neighbors, CART and an artificial neural network. A stack of 585 coarse and high resolution covariates representing spectral reflectance (Landsat bands, spectral indices; time-series of seasonal composites), different biophysical conditions (i.e. temperature, precipitation, elevation, lithology) and biotic competition (other species distribution maps) was used as predictors for realized distributions, while potential distribution was modelled with environmental predictors only. Logloss and computing time were used to select the three best algorithms to train an ensemble model based on stacking with a logistic regressor as a meta-learner for each species. High resolution (30 m) probability and model uncertainty maps of realized distribution were produced for each species using a time window of 4 years for a total of 6 distribution maps per species for the studied period, while for potential distributions only one map per species was produced. Results of spatial cross validation show that Olea europaea and Quercus suber achieved the best performances in both potential and realized distribution, while Pinus sylvestris and Salix caprea achieved the worst. Further analysis shows that fine-resolution models consistently outperformed coarse resolution models (250 m) for realized distribution (average decrease in logloss: +53%). Realized distribution models achieved higher predictive performances than potential distribution ones. Importance of predictor variables differed across species and models, with the green band for summer and the NDWI and NDVI for fall for realized distribution and the diffuse irradiation and precipitation of the driest quarter being the most important and frequent for potential distribution. The ensemble model outperformed or performed as good as the best individual model in all potential species distributions, while for ten species it performed worse than the best individual model in modeling realized distributions. The framework shows how combining continuous and consistent EO time series data with state of the art ML can be used to derive dynamic distribution maps. The produced time-series occurrence predictions can be used to quantify temporal trends and detect potential forest degradation.