diff --git a/mnp/species_models/habitat_suitability.py b/mnp/species_models/habitat_suitability.py
index c6c614b71b5befe948f7b149baf488769c7eeac6..8bf126af819367bc76b014f2b1bed70b16384a63 100644
--- a/mnp/species_models/habitat_suitability.py
+++ b/mnp/species_models/habitat_suitability.py
@@ -60,9 +60,10 @@ class Nature(HSIProcedure):
                 if land_type in land_types.keys():
                     # Just add land type suitabilities if not running with environmentals
                     if environmental_factor_array is None:
-                        hsi += land_types[land_type].multiply(
-                            self._suitabilities[land_type]
-                        )
+                        if self._suitabilities[land_type] >= self._hsi_threshold:
+                            hsi += land_types[land_type].multiply(
+                                self._suitabilities[land_type]
+                            )
 
                     # Otherwise, multiply hsi with environmental array
                     else: