diff --git a/R/biomod2_classes_0.R b/R/biomod2_classes_0.R
index e4713ca7..9087fcfc 100644
--- a/R/biomod2_classes_0.R
+++ b/R/biomod2_classes_0.R
@@ -509,6 +509,10 @@ setMethod('BIOMOD.options.dataset', signature(strategy = 'character'),
               argsval <- bm_Tuning(model = mod, tuning.fun = tuning.fun, do.formula = TRUE, do.stepAIC = TRUE
                                    , bm.options = BOM, bm.format = bm.format, calib.lines = calib.lines
                                    , metric.eval = ifelse(mod == "MAXENT", "auc.val.avg", "TSS"))
+              
+              # To avoid problem with ENMeval and "bad closing" of foreach loop
+              env <- foreach:::.foreachGlobals
+              rm(list = ls(name = env), pos = env)
             }
             BOD <- new('BIOMOD.options.dataset', BOM, args.values = argsval)