diff --git a/src/evalhyd-python.cpp b/src/evalhyd-python.cpp
index 9a7bb358e5750432c5b58de01ee8c2eaada4feaf..8cf2c41198db7ea4ccb4b3c57b50e7e364e6c404 100644
--- a/src/evalhyd-python.cpp
+++ b/src/evalhyd-python.cpp
@@ -98,9 +98,7 @@ PYBIND11_MODULE(evalhyd, m)
                    function with a negative exponent (since none are defined
                    for 0). If not provided (or set to default value -9),
                    one hundredth of the mean of the streamflow observations
-                   is used as value for epsilon, as recommended by
-                   `Pushpalatha et al. (2012)
-                   <https://doi.org/10.1016/j.jhydrol.2011.11.055>`_.
+                   is used as value for epsilon.
 
                 t_msk: `numpy.ndarray`, optional
                    1D array of mask(s) used to generate temporal subsets of
@@ -204,9 +202,7 @@ PYBIND11_MODULE(evalhyd, m)
                    function with a negative exponent (since none are defined
                    for 0). If not provided (or set to default value -9),
                    one hundredth of the mean of the streamflow observations
-                   is used as value for epsilon, as recommended by
-                   `Pushpalatha et al. (2012)
-                   <https://doi.org/10.1016/j.jhydrol.2011.11.055>`_.
+                   is used as value for epsilon.
 
                 t_msk: `numpy.ndarray`, optional
                    2D array of mask(s) used to generate temporal subsets of