diff --git a/envs/environment.fixed.yaml b/envs/environment.fixed.yaml
index 1c7aeaad..94f2c9ce 100644
--- a/envs/environment.fixed.yaml
+++ b/envs/environment.fixed.yaml
@@ -335,6 +335,7 @@ dependencies:
 - requests=2.28.1
 - retry=0.9.2
 - rich=12.5.1
+- rioxarray=0.13.3
 - rtree=1.0.0
 - s2n=1.0.10
 - scikit-learn=1.1.1
diff --git a/envs/environment.yaml b/envs/environment.yaml
index 490268e8..74a63eb8 100644
--- a/envs/environment.yaml
+++ b/envs/environment.yaml
@@ -29,6 +29,7 @@ dependencies:
 - pandas
 - geopandas>=0.11.0
 - xarray
+- rioxarray
 - netcdf4
 - networkx
 - scipy
diff --git a/scripts/build_ship_raster.py b/scripts/build_ship_raster.py
index 0ae134fd..fcfa632d 100644
--- a/scripts/build_ship_raster.py
+++ b/scripts/build_ship_raster.py
@@ -41,7 +41,7 @@ import logging
 import os
 import zipfile
 
-import xarray as xr
+import rioxarray as rio
 from _helpers import configure_logging
 from build_natura_raster import determine_cutout_xXyY
 
@@ -59,7 +59,7 @@ if __name__ == "__main__":
 
     with zipfile.ZipFile(snakemake.input.ship_density) as zip_f:
         zip_f.extract("shipdensity_global.tif")
-        with xr.open_rasterio("shipdensity_global.tif") as ship_density:
+        with rio.open_rasterio("shipdensity_global.tif") as ship_density:
             ship_density = ship_density.drop(["band"]).sel(
                 x=slice(min(xs), max(Xs)), y=slice(max(Ys), min(ys))
             )