diff --git a/docs/sphinx/source/whatsnew/v0.16.2.rst b/docs/sphinx/source/whatsnew/v0.16.2.rst
index 537f8db20e..9c0c22a172 100644
--- a/docs/sphinx/source/whatsnew/v0.16.2.rst
+++ b/docs/sphinx/source/whatsnew/v0.16.2.rst
@@ -25,6 +25,9 @@ Documentation
 * Clarified in :py:func:`pvlib.tracking.singleaxis` that the tracker
   rotation minimizes the angle of incidence (AOI) to maximize the use of
   direct normal irradiance (DNI). (:issue:`2314`, :pull:`2869`)
+* Fix the :py:func:`pvlib.spa.solar_position` docstring, which listed
+  outputs ``elevation`` and ``apparent_elevation`` in the wrong order.
+  (:issue:`2886`, :pull:`2887`)
 
 
 Testing
diff --git a/pvlib/spa.py b/pvlib/spa.py
index dba1e2ae84..f7b2736b1f 100644
--- a/pvlib/spa.py
+++ b/pvlib/spa.py
@@ -1075,8 +1075,8 @@ def solar_position(unixtime, lat, lon, elev, pressure, temp, delta_t,
     Numpy Array with elements:
         apparent zenith,
         zenith,
-        elevation,
         apparent_elevation,
+        elevation,
         azimuth,
         equation_of_time