diff --git a/ci/environment.yml b/ci/environment.yml
index 979b5117ff..b7c87939cf 100644
--- a/ci/environment.yml
+++ b/ci/environment.yml
@@ -47,6 +47,7 @@ dependencies:
   - bokeh ==3.5.1
   - gilknocker ==0.4.1
   - openssl >1.1.0g
+  - rasterio >=1.4.0
   - rioxarray ==0.17.0
   - h5netcdf ==1.3.0
   - xesmf ==0.8.7
diff --git a/tests/geospatial/test_zonal_average.py b/tests/geospatial/test_zonal_average.py
index 3eec34d2d7..b74c93b574 100644
--- a/tests/geospatial/test_zonal_average.py
+++ b/tests/geospatial/test_zonal_average.py
@@ -37,8 +37,9 @@ def test_nwm(
         subset = ds.zwattablrt.sel(time=time_range)
 
         counties = rioxarray.open_rasterio(
-            s3.open("s3://nwm-250m-us-counties/Counties_on_250m_grid.tif"),
+            "s3://nwm-250m-us-counties/Counties_on_250m_grid.tif",
             chunks="auto",
+            opener=s3.open,
         ).squeeze()
 
         # Remove any small floating point error in coordinate locations