diff --git a/paper/paper.bib b/paper/paper.bib
index e41fa07a..57eb1801 100644
--- a/paper/paper.bib
+++ b/paper/paper.bib
@@ -299,3 +299,11 @@ @article{Tsujino2018
   volume  = {130},
   year    = {2018},
 }
+
+@article{Silvestri2024GPU,
+  author  = {Simone Silvestri and Gregory L. Wagner and Navid C Constantinou and Christopher Hill and Jean-Michel Campin and Andre Nogueira Souza and Siddhartha Bishnu and Valentin Churavy and John C Marshall and Raffaele Ferrari},
+  title   = {A {GPU}-based ocean dynamical core for routine mesoscale-resolving climate simulations},
+  journal = {ESS Open Archive},
+  year    = {2024},
+  doi     = {10.22541/essoar.171708158.82342448/v1},
+}
diff --git a/paper/paper.md b/paper/paper.md
index 66130b1d..3bba6eed 100644
--- a/paper/paper.md
+++ b/paper/paper.md
@@ -61,7 +61,7 @@ bibliography: paper.bib
 ## Regional ocean modeling
 
 In the ocean, fast and small-scale motions (from ~100m to ~100km varying at time scales of hours to days) play an important role in shaping the large-scale ocean circulation and climate (length scales ~10,000km varying at decadal time scales) [@Melet2022ch2; @deLavergne2022ch3; @Gula2022ch8].
-Despite the increase in computational power and the use of graphical processing units that can bring breakthrough performance and speedup [@silvestri2023oceananigansjl], there are always processes, boundary, or forcing features that are smaller than the model's grid spacing and, thus, remain unresolved in global ocean models.
+Despite the increase in computational power and the use of graphical processing units that can bring breakthrough performance and speedup [@Silvestri2024GPU], there are always processes, boundary, or forcing features that are smaller than the model's grid spacing and, thus, remain unresolved in global ocean models.
 Regional ocean models can be run at higher resolutions while limiting the required computational resources.
 
 A regional ocean model simulates the ocean only in a prescribed region, which is a subset of the global ocean.