diff --git a/doc/references.rst b/doc/references.rst
index 013ca3182..9b1b9badf 100644
--- a/doc/references.rst
+++ b/doc/references.rst
@@ -7,10 +7,10 @@ References
 .. [Cooper2000] Cooper, G.R.J. (2000), Gridding gravity data using an equivalent layer, Computers & Geosciences, Computers & Geosciences, doi:`10.1016/S0098-3004(99)00089-8 <https://doi.org/10.1016/S0098-3004(99)00089-8>`__
 .. [Dampney1969] Dampney, C. N. G. (1969). The equivalent source technique. Geophysics, 34(1), 39–53. doi:`10.1190/1.1439996 <https://doi.org/10.1190/1.1439996>`__
 .. [Forste_etal2014] Förste, Christoph; Bruinsma, Sean.L.; Abrikosov, Oleg; Lemoine, Jean-Michel; Marty, Jean Charles; Flechtner, Frank; Balmino, G.; Barthelmes, F.; Biancale, R. (2014): EIGEN-6C4 The latest combined global gravity field model including GOCE data up to degree and order 2190 of GFZ Potsdam and GRGS Toulouse. GFZ Data Services. doi:`10.5880/icgem.2015.1 <http://doi.org/10.5880/icgem.2015.1>`__
-.. [Fukushima2019] Fukushima, T. (2019). Fast computation of prismatic gravitational field. doi:`10.13140/RG.2.2.30598.93766 <https://doi.org/10.13140/RG.2.2.30598.93766>`__
+.. [Fukushima2020] Fukushima, T. (2020). Speed and accuracy improvements in standard algorithm for prismatic gravitational field. Geophysical Journal International. doi:`10.1093/gji/ggaa240 <https://doi.org/10.1093/gji/ggaa240>`__
 .. [Grombein2013] Grombein, T., Seitz, K., Heck, B. (2013), Optimized formulas for the gravitational field of a tesseroid, Journal of Geodesy. doi:`10.1007/s00190-013-0636-1 <https://doi.org/10.1007/s00190-013-0636-1>`__
 .. [Hofmann-WellenhofMoritz2006] Hofmann-Wellenhof, B., & Moritz, H. (2006). Physical Geodesy (2nd, corr. ed. 2006 edition ed.). Wien ; New York: Springer.
 .. [Nagy2000] Nagy, D., Papp, G. & Benedek, J.(2000). The gravitational potential and its derivatives for the prism. Journal of Geodesy 74: 552. doi:`10.1007/s001900000116 <https://doi.org/10.1007/s001900000116>`__
-.. [Nagy2002] Nagy, D., Papp, G. & Benedek, J.(2000). Corrections to "The gravitational potential and its derivatives for the prism". Journal of Geodesy (2002) 76: 475. doi:`10.1007/s00190-002-0264-7 <https://doi.org/10.1007/s00190-002-0264-7>`__
+.. [Nagy2002] Nagy, D., Papp, G. & Benedek, J.(2002). Corrections to "The gravitational potential and its derivatives for the prism". Journal of Geodesy. doi:`10.1007/s00190-002-0264-7 <https://doi.org/10.1007/s00190-002-0264-7>`__
 .. [Vajda2004] Vajda, P., Vaníček, P., Novák, P. and Meurers, B. (2004). On evaluation of Newton integrals in geodetic coordinates: Exact formulation and spherical approximation. Contributions to Geophysics and Geodesy, 34(4), 289-314.
 .. [TurcotteSchubert2014] Turcotte, D. L., & Schubert, G. (2014). Geodynamics (3 edition). Cambridge, United Kingdom: Cambridge University Press.
diff --git a/harmonica/forward/prism.py b/harmonica/forward/prism.py
index 7b4e82b2d..dc035fabb 100644
--- a/harmonica/forward/prism.py
+++ b/harmonica/forward/prism.py
@@ -31,7 +31,7 @@ def prism_gravity(
     inside the prism.
 
     This implementation makes use of the modified arctangent function proposed
-    by [Fukushima2019]_ (eq. 12) so that the potential field to satisfies
+    by [Fukushima2020]_ (eq. 12) so that the potential field to satisfies
     Poisson's equation in the entire domain. Moreover, the logarithm function
     was also modified in order to solve the singularities that the analytical
     solution has on some points (see [Nagy2000]_).
@@ -280,7 +280,7 @@ def safe_atan2(y, x):
     gravitational field of the prism satisfies the Poisson's equation.
     Therefore, it guarantees that the fields satisfies the symmetry properties
     of the prism. This modified function has been defined according to
-    [Fukushima2019]_.
+    [Fukushima2020]_.
     """
     if x != 0:
         result = np.arctan(y / x)