diff --git a/autotest/test_mf6.py b/autotest/test_mf6.py
index 9a1a43b6d..7ea5ab6b0 100644
--- a/autotest/test_mf6.py
+++ b/autotest/test_mf6.py
@@ -2332,6 +2332,8 @@ def test_remove_model(function_tmpdir, example_data_path):
             break
 
 
+@requires_pkg("shapely")
+@requires_exe("triangle")
 def test_flopy_2283(function_tmpdir):
     # create triangular grid
     triangle_ws = function_tmpdir / "triangle"
diff --git a/autotest/test_model_splitter.py b/autotest/test_model_splitter.py
index a05708ec4..09de92752 100644
--- a/autotest/test_model_splitter.py
+++ b/autotest/test_model_splitter.py
@@ -158,8 +158,8 @@ def test_model_with_lak_sfr_mvr(function_tmpdir):
     np.testing.assert_allclose(new_heads, original_heads, err_msg=err_msg)
 
 
-@requires_pkg("pymetis")
 @requires_exe("mf6")
+@requires_pkg("pymetis")
 @pytest.mark.slow
 def test_metis_splitting_with_lak_sfr(function_tmpdir):
     sim_path = get_example_data_path() / "mf6" / "test045_lake2tr"
@@ -848,8 +848,7 @@ def test_unstructured_complex_disu(function_tmpdir):
 
 
 @requires_exe("mf6")
-@requires_pkg("pymetis")
-@requires_pkg("scipy")
+@requires_pkg("pymetis", "scipy")
 def test_multi_model(function_tmpdir):
     from scipy.spatial import KDTree
 
@@ -1301,6 +1300,7 @@ def build_gwt_model(sim, gwtname, rch_package):
 
 
 @requires_exe("mf6")
+@requires_pkg("pymetis")
 def test_timeseries(function_tmpdir):
     sim = MFSimulation(
         sim_name="np001",
diff --git a/autotest/test_plot_cross_section.py b/autotest/test_plot_cross_section.py
index a0f058e66..c478b6374 100644
--- a/autotest/test_plot_cross_section.py
+++ b/autotest/test_plot_cross_section.py
@@ -217,6 +217,7 @@ def test_plot_limits():
     plt.close(fig)
 
 
+@requires_pkg("shapely")
 def test_plot_centers():
     from matplotlib.collections import PathCollection
 
diff --git a/autotest/test_plot_map_view.py b/autotest/test_plot_map_view.py
index 386242d1d..bf6befb88 100644
--- a/autotest/test_plot_map_view.py
+++ b/autotest/test_plot_map_view.py
@@ -276,6 +276,7 @@ def test_plot_limits():
     plt.close(fig)
 
 
+@requires_pkg("shapely")
 def test_plot_centers():
     nlay = 1
     nrow = 10