In this interesting and original study, the authors present an ensemble Machine Learning (ML) model for the prediction of the habitats' suitability, which is affected by the complex interactions between living conditions and survival-spreading climate factors. The research focuses in two of the most dangerous invasive mosquito species in Europe with the requirements' identification in temperature and rainfall conditions. Though it is an interesting approach, the ensemble ML model is not presented and discussed in sufficient detail and thus its performance and value as a tool for modeling the distribution of invasive species cannot be adequately evaluated.

The authors use an Ensemble Approach (ENAP) based on 10 timely ML algorithms, aiming to draw up the habitats' maps for both species of mosquitoes. Ensemble methods are meta-algorithms that combine several techniques into a unique predictive model to decrease variance. For example, in Bagging different training data subsets are randomly drawn—with replacement—from the entire training dataset, to train a different classifier. In Boosting, resampling is strategically geared to provide the most informative training data for each consecutive classifier, or to improve predictions. Stacking, involves training to combine the predictions of several other learning algorithms (Zhou, 2012).

Unlike a statistical ensemble in statistical mechanics which is usually infinite, a ML ensemble consists of only a concrete finite set of alternative models, but typically allows for much more flexible structures to exist among those alternatives. Perhaps one of the earliest works on ensemble systems is the paper by Dasarathy and Sheela (1979). They first introduced an ENAP for partitioning the feature space, using two or more classifiers, in a divide-and-conquer fashion. Over a decade later, Hansen and Salamon (1990) showed the variance reduction property of an ENAP. They managed to improve the generalization performance of an ANN by using an ensemble of similarly configured ANN. But it was Schapire's work that has put the ENAP at the center of ML research, as he has proven that a strong classifier can be generated by combining weak classifiers (Schapire, 1990). Finally, Buisson et al. suggested that attention should be paid to the use of predictions ensembles resulting from the application of several statistical methods. Forecasted impacts should always be provided with an assessment of their uncertainty (Buisson et al., 2010).

Unfortunately, the authors of this interesting paper, do not offer a deep description of the proposed ENAP and it is not clear if their approach can cover the main points of the ensemble techniques. For example, the proposed ENAP convert species' probability of occurrence into binary presence-absence data using a predefined threshold. Assessing models based on presence only data, it is difficult to learn the overall species occurrence probability, based on false or misleading information or unjustified simplifying assumptions, because there is typically no validation data with true presences and absences (Hastie and Fithian, 2013). The ENAP that was proposed cannot surmount this problem, it only makes it more hidden.

Current practices in Alien Species Distribution Modeling (ASDM) algorithms (Lorena et al., 2011; Duan et al., 2014; Shabani et al., 2016), include Profile Methods (BIOCLIM, ENFA) (Lorena et al., 2011; Duan et al., 2014; Shabani et al., 2016), Regression-based techniques (GLM, MARS) (Lorena et al., 2011; Duan et al., 2014; Shabani et al., 2016), ML techniques (MAXENT, ANN, SVM) (Lorena et al., 2011; Duan et al., 2014; Shabani et al., 2016).

A widely used and effective method in ASDM involves creating ML ensembles' models (Duan et al., 2014). The two most important advantages of ENAP focus on the fact that they offer better prediction and more stable and robust models, as the overall behavior of a multiple model is less noisy than a corresponding single one (Kuncheva, 2004; Zhou, 2012). For example, in Zhang and Zhang (2012) the authors propose an effective ENAP to assess the impacts of predictor variables and ASDM. In Daliakopoulos et al. (2017) the Random Forest EANP has proven that it can provide a better understanding of facilitating and limiting factors of alien species presence, both for research and management purposes. Finally, Lauzeral et al. (2012) proposes an iterative ENAP to ensure noise absence and hence to improve the predictive reliability of ensemble modeling of species distributions.

Some of the most important points related to the operation and use of the ENAP that should be included and discussed thoroughly by the authors are presented below:

It is worth noting that in general an ENAP can lead to much better prediction results, while offering generalization. This is one of the key issues in the field of ML, as it can reduce bias and variance and it has the potential to eliminate overfitting (Kuncheva, 2004; Zhou, 2012). Moreover, it implements robust predictive models capable of responding to high complexity problems such as those of spreading invasive species (Demertzis and Iliadis, 2015, 2017). However, the development of these models should not be done in a black box mode research and it should be accompanied by a set of in-depth analysis regarding key training and operation decision points, thus allowing critical readers to fully and thoroughly evaluate the proposed methodology and to promote research in the broader scope. Finally, there are cases where wide variety of comparatively model-free forecasting methods outperforms the correct mechanistic data-driven model. However, according to Moustakas (2017) “if one simply relies on data-driven science, several components of scientific methods will be made poorer.”

KD and LI conceived of the presented idea. KD and LI developed the theoretical formalism, performed the analytic calculations and performed the numerical simulations. KD, LI, and V-DA verified the analytical methods. All the authors contributed to the final version of the manuscript. LI supervised the project.