diff --git a/docs/sphinx/source/whatsnew/v0.15.3.rst b/docs/sphinx/source/whatsnew/v0.15.3.rst
index 23d89b06e5..a6a66c24da 100644
--- a/docs/sphinx/source/whatsnew/v0.15.3.rst
+++ b/docs/sphinx/source/whatsnew/v0.15.3.rst
@@ -12,6 +12,7 @@ Deprecations
 ~~~~~~~~~~~~
 
 
+
 Bug fixes
 ~~~~~~~~~
 
diff --git a/docs/sphinx/source/whatsnew/v0.16.0.rst b/docs/sphinx/source/whatsnew/v0.16.0.rst
index 6d62b9e1c6..879cf25641 100644
--- a/docs/sphinx/source/whatsnew/v0.16.0.rst
+++ b/docs/sphinx/source/whatsnew/v0.16.0.rst
@@ -15,7 +15,9 @@ Breaking Changes
 
 Deprecations
 ~~~~~~~~~~~~
-
+* Deprecate :py:func:`pvlib.irradiance.king`.
+  Use other diffuse transposition models in :py:mod:`pvlib.irradiance` instead.
+  (:issue:`2636`, :pull:`2783`)
 
 Bug fixes
 ~~~~~~~~~
diff --git a/pvlib/irradiance.py b/pvlib/irradiance.py
index 98bd948286..67e860b249 100644
--- a/pvlib/irradiance.py
+++ b/pvlib/irradiance.py
@@ -987,6 +987,12 @@ def reindl(surface_tilt, surface_azimuth, dhi, dni, ghi, dni_extra,
     return sky_diffuse
 
 
+@deprecated(
+    since="0.15.3",
+    removal="0.17.0",
+    name="pvlib.irradiance.king",
+    alternative="other diffuse transposition models in pvlib.irradiance",
+)
 def king(surface_tilt, dhi, ghi, solar_zenith):
     '''
     Determine diffuse irradiance from the sky on a tilted surface using
diff --git a/tests/test_irradiance.py b/tests/test_irradiance.py
index fbe9906358..8607074d59 100644
--- a/tests/test_irradiance.py
+++ b/tests/test_irradiance.py
@@ -250,8 +250,9 @@ def test_reindl(irrad_data, ephem_data, dni_et):
 
 
 def test_king(irrad_data, ephem_data):
-    result = irradiance.king(40, irrad_data['dhi'], irrad_data['ghi'],
-                             ephem_data['apparent_zenith'])
+    with pytest.warns(pvlibDeprecationWarning, match='king'):
+        result = irradiance.king(40, irrad_data['dhi'], irrad_data['ghi'],
+                                 ephem_data['apparent_zenith'])
     assert_allclose(result, [0, 44.629352, 115.182626, 79.719855], atol=1e-4)