This paper is of high quality, it exemplifies the complexity of detecting ocean eddies (eulerian versus lagrangian, importance of the code parameters of the regions considered), it acknowledges the fact that not all ocean eddies are the same (dynamically talking), and that their effect on plankton is complex. The method presented is well-described, powerful and would allow digging into these biophysical interactions in many different oceanographic settings. It shows that the common view following seminal papers by e.g. McGillicuddy and D. Chelton that consists of "cyclonic eddies -> isopycnals shaoling -> increase phytoplankton biomass" versus "anticyclonic eddies -> isopycnal deepening -> decreased phytoplankton biomass" is too simplistic and needed to be updated.