diff --git a/imod/prepare/reproject.py b/imod/prepare/reproject.py
index 621be048d..3124ddc1f 100644
--- a/imod/prepare/reproject.py
+++ b/imod/prepare/reproject.py
@@ -1,3 +1,5 @@
+import warnings
+
 import affine
 import numpy as np
 import xarray as xr
@@ -46,17 +48,26 @@ def reproject(
     **reproject_kwargs,
 ):
     """
+
+    .. attention::
+
+        This function is deprecated and will be removed in a future version. Use
+        `rioxarray.reproject
+        <https://corteva.github.io/rioxarray/stable/rioxarray.html#rioxarray.raster_array.RasterArray.reproject>`_
+        instead.
+
+    .. warning::
+
+        The `like` paramater does not seem to resample anymore as expected. For
+        regridding, consider using `xugrid to regrid instead.
+        <https://deltares.github.io/xugrid/examples/regridder_overview.html>`_
+
     Reprojects and/or resamples a 2D xarray DataArray to a
     different cellsize or coordinate system.
 
-    * To resample to a new cellsize in the same projection: provide only ``like``.
-    * To only reproject: provide only ``src_crs`` and ``src_crs``.
+    * To reproject: provide only ``src_crs`` and ``dst_crs``.
     * To reproject and resample to a specific domain: provide ``src_crs``, ``src_crs``, and ``like``.
 
-    Note: when only ``like`` is provided, Cartesian (projected) coordinates are a
-    ssumed for resampling. In case of non-Cartesian coordinates, specify
-    ``src_crs`` and ``dst_crs`` for correct resampling.
-
     Parameters
     ----------
     source: xarray DataArray
@@ -65,6 +76,13 @@ def reproject(
     like: xarray DataArray
         Example DataArray that shows what the resampled result should look like
         in terms of coordinates. Must contain dimensions ``y`` and ``x``.
+
+        .. warning::
+
+            The `like` parameter does not seem to resample anymore due to
+            changes in rasterio. For regridding, consider using `xugrid to regrid
+            instead.
+            <https://deltares.github.io/xugrid/examples/regridder_overview.html>`_
     src_crs: string, dict, rasterio.crs.CRS
         Coordinate system of ``source``. Options:
 
@@ -105,38 +123,32 @@ def reproject(
 
     Examples
     --------
-    Resample a DataArray ``a`` to a new cellsize, using an existing DataArray ``b``:
-
-    >>> c = imod.rasterio.reproject(source=a, like=b)
-
-    Resample a DataArray to a new cellsize of 100.0, by creating a ``like`` DataArray first:
-    (Note that dy must be negative, as is usual for geospatial grids.)
-
-    >>> dims = ("y", "x")
-    >>> coords = {"y": np.arange(200_000.0, 100_000.0, -100.0), "x": np.arange(0.0, 100_000.0, 100.0)}
-    >>> b = xr.DataArray(data=np.empty((200, 100)), coords=coords, dims=dims)
-    >>> c = imod.rasterio.reproject(source=a, like=b)
-
     Reproject a DataArray from one coordinate system (WGS84, EPSG:4326) to another (UTM30N, EPSG:32630):
 
-    >>> c = imod.rasterio.reproject(source=a, src_crs="EPSG:4326", dst_crs="EPSG:32630")
+    >>> c = imod.prepare.reproject(source=a, src_crs="EPSG:4326", dst_crs="EPSG:32630")
 
     Get the reprojected DataArray in the desired shape and coordinates by providing ``like``:
 
-    >>> c = imod.rasterio.reproject(source=a, like=b, src_crs="EPSG:4326", dst_crs="EPSG:32630")
+    >>> c = imod.prepare.reproject(source=a, like=b, src_crs="EPSG:4326", dst_crs="EPSG:32630")
 
     Open a single band raster, and reproject to RD new coordinate system (EPSG:28992), without explicitly specifying ``src_crs``.
     ``src_crs`` is taken from ``a.attrs``, so the raster file has to include coordinate system metadata for this to work.
 
     >>> a = rioxarray.open_rasterio("example.tif").squeeze("band")
-    >>> c = imod.rasterio.reproject(source=a, use_src_attrs=True, dst_crs="EPSG:28992")
+    >>> c = imod.prepare.reproject(source=a, use_src_attrs=True, dst_crs="EPSG:28992")
 
     In case of a rotated ``source``, provide ``src_transform`` directly or ``use_src_attrs=True`` to rely on generated attributes:
 
     >>> rotated = rioxarray.open_rasterio("rotated_example.tif").squeeze("band")
-    >>> c = imod.rasterio.reproject(source=rotated, dst_crs="EPSG:28992", reproject_kwargs={"src_transform":affine.Affine(...)})
-    >>> c = imod.rasterio.reproject(source=rotated, dst_crs="EPSG:28992", use_src_attrs=True)
+    >>> c = imod.prepare.reproject(source=rotated, dst_crs="EPSG:28992", reproject_kwargs={"src_transform":affine.Affine(...)})
+    >>> c = imod.prepare.reproject(source=rotated, dst_crs="EPSG:28992", use_src_attrs=True)
     """
+
+    warnings.warn(
+        "imod.prepare.reproject is deprecated and will be removed in a future version.",
+        DeprecationWarning,
+    )
+
     # Make sure the rio accessor is avaible.
     import rioxarray  # noqa pylint: F401