diff --git a/CHANGELOG.rst b/CHANGELOG.rst
index c287bae6..ff9b03ed 100644
--- a/CHANGELOG.rst
+++ b/CHANGELOG.rst
@@ -13,6 +13,7 @@ Bug fixes
 ^^^^^^^^^
 * Fixed bug affecting `GriddedForcing`, where `station_idx` in the call to ``nc_specs`` was set to ``1`` instead of ``None``. (PR #501)
 * Fixed bug in the `run` method, where `overwrite=False` was not being respected. (PR #503)
+* Pin `scipy` below v1.16.0 due to a breaking change that affects `statsmodels` below v0.14.4. (PR #521)
 
 Internal changes
 ^^^^^^^^^^^^^^^^
diff --git a/environment-dev.yml b/environment-dev.yml
index c9d2605a..d32ff9e9 100644
--- a/environment-dev.yml
+++ b/environment-dev.yml
@@ -30,7 +30,7 @@ dependencies:
   - pyproj >=3.3.0
   - rasterio
   - rioxarray
-  - scipy >=1.11.0
+  - scipy >=1.11.0,<1.16.0 # See: https://github.com/statsmodels/statsmodels/issues/9584
   - shapely >=2.0.0
   - spotpy >=1.6.1
   - statsmodels >=0.14.2
diff --git a/pyproject.toml b/pyproject.toml
index b8348821..07635ced 100644
--- a/pyproject.toml
+++ b/pyproject.toml
@@ -52,7 +52,7 @@ dependencies = [
   "pydantic >=2.11",
   "pydap >=3.4.0,<3.5.5", # pydap 3.5.5 is not currently supported by `xarray` (v2025.3.1)
   "pymbolic >=2024.2",
-  "scipy >=1.11.0",
+  "scipy >=1.11.0,<1.16.0", # See: https://github.com/statsmodels/statsmodels/issues/9584
   "spotpy >=1.6.1",
   "statsmodels >=0.14.2",
   "typing-extensions",