# Atmospheric-ocean coupling drives prevailing and synchronic dispersal patterns of marine species with long pelagic durations

Eduardo Ramirez-Romero, Angel Amores, David Diaz, Anabel Muñoz, Ignacio Catalan, Juan Carlos Molinero, Andres Ospina-Alvarez

<figure><img src="image_1.png" /></figure>

## To cite this version:

Eduardo Ramirez-Romero, Angel Amores, David Diaz, Anabel Muñoz, Ignacio Catalan, et al. "Atmospheric-ocean coupling drives prevailing and synchronic dispersal patterns of marine species with long pelagic durations." Scientific Reports, 2023, 13 (1), pp.2366. 10.1038/s41598-023-29543-7. hal-03994928

HAL Id: hal-03994928

https://hal.umontpellier.fr/hal-03994928

Submitted on 5 Jun 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Atmospheric-ocean coupling drives prevailing and synchronic dispersal patterns of marine species with long pelagic durations

Eduardo Ramirez-Romero<fnref n="1" />, Angel Amores², David Diaz⁴, Anabel Muñoz⁴, Ignacio A. Catalan³, Juan Carlos Molinero⁵ & Andres Ospina-Alvarez³✉

Dispersal shapes population connectivity and plays a critical role in marine metacommunities. Prominent species for coastal socioecological systems, such as jellyfish and spiny lobsters, feature long pelagic dispersal phases (LPDPs), which have long been overlooked. Here, we use a cross-scale approach combining field surveys of these species with a high-resolution hydrodynamic model to decipher the underlying mechanisms of LPDP patterns in northwestern Mediterranean shores. We identified basin-scale prevailing dispersal routes and synchronic year-to-year patterns tightly linked to prominent circulation features typical of marginal seas and semiclosed basins, with an outstanding role of a retentive source area replenishing shores and potentially acting as a pelagic nursery area. We show how the atmospheric forcing of the ocean, a marked hydrological driver of the Mediterranean Sea, modulates dispersal routes and sources of LPDP at interannual scales. These findings represent a crucial advance in our understanding of the functioning of metapopulations of species with LPDP in marginal seas and may contribute to the effective management of coastal ecosystem services in the face of climate change.

Marine populations inhabit a fragmented seascape where dispersal, i.e., the movement of organisms or their off-spring away from their source to a settlement or nursery site, shapes gene flow and structures the connectivity and persistence of populations, thereby sculpting metacommunities<fnref n="1" />, dynamics and seashore biogeography². Given the pivotal role of dispersal, a mechanistic understanding is essential for both theoretical (i.e., evolution) and management (i.e., conservation) reasons, as rapid anthropogenic impact is altering marine seascapes³, affecting species distribution and consequently nature's contributions to people and human wellbeing⁴.

Nearshore ecosystems host a wide diversity of invertebrates featuring pelagic larval phases prone to dispersal by ocean currents. This process connects spawning locations with suitable nursery areas or settlement habitats. Dispersal takes place at multiple scales and is driven by a complex array of processes and interactions, namely, the lifespan of oceanographic features, temperature⁵, species life-history traits, including planktonic larval duration (PLD), and behavior (e.g., vertical migration)⁶,⁷. Efforts to unveil population connectivity have thus far mainly focused on short-term scales, ≈ 2–4 weeks⁶, using proxies combining PLD and genetic metrics⁸, and high-resolution biophysical coupled models⁹. Overall, at these scales, dispersion has been described as a highly stochastic and uncertain process⁹, while at long-term scales, the process remains less understood. Indeed, the dispersal and connectivity of species with LPDP (> 4 months) have long been overlooked, partly due to traditional research mainly focused on the late phase of dispersal, i.e., arrival at nursery and stranding areas, settlement and temporal correlations¹⁰⁻¹², and partly due to the difficulty of making synoptic observations over broad spatial and temporal scales. Long-lasting plankton forms are shaped by interlinked biological, i.e., growth and mortality, and physical processes ranging from transient features, namely, fronts and mesoscale eddies, to basin-scale persistent physical phenomena, i.e., geostrophic currents¹³. Hence, dispersal and connectivity over broad scales may be nurtured by heterogeneous oceanic processes presenting different predictability levels⁹,¹⁴⁻¹⁶. In this regard,

<footnote n="1">¹Instituto de Ciencias Marinas de Andalucía, ICMAN, CSIC, República Saharaui, 4, Puerto Real, 11519 Cádiz, Spain. ²Department of Physics, University of the Balearic Islands, Crta. Valldemossa, km 7.5, 07122 Palma, Illes Balears, Spain. ³Mediterranean Institute for Advanced Studies, IMEDEA, CSIC-UIB, Miquel Marquès, 21, 07190 Esporles, Illes Balears, Spain. ⁴Instituto Español de Oceanografía, IEO, CSIC, Centre Oceanogràfic de les Balears, Moll de Ponent sn, 07015 Palma, Illes Balears, Spain. ⁵Marine Biodiversity Exploitation and Conservation, MARBEC, IRD, CNRS, Ifremer, Université de Montpellier, 610101 Sète, France. ✉email: aospina.co@me.com</footnote>
genetic studies can integrate to a certain extent the size and dynamics of metapopulations, including the effect of adults migrations (e.g. for fish)¹⁷,¹⁸ but may fail to capture underlying mechanisms linking hydrodynamic processes and species life cycles and derived interannual changes in population size. To fill this gap, we propose a complementary cross-scale study combining field surveys and state-of-the-art hydrodynamic modeling to reconstruct potential dispersal routes and entire dispersal kernels and to gain a mechanistic understanding of processes and scales shaping metapopulations.

In densely populated basins such as the Mediterranean Sea, some iconic species providing worthy natural contributions to people are featured by the LPDP. An example is the European spiny lobster (*Palinurus elephas*), a major target of artisanal fisheries since ancient times¹⁹, which faces overexploitation and is included in the red list of the International Union for Conservation of Nature's Red List of Threatened Species²⁰. Another prominent species also characterized by LPDP is the mauve stinger jellyfish (*Pelagia noctiluca*), which plays a key role in the pelagic food web and severely affects economic activities, including fisheries, aquaculture and tourism²¹,²². These species challenge conservation and management actions in the Mediterranean Sea, and therefore, resolving their dispersal patterns and connectivity has widespread implications for coastal communities and pivotal sectors. Tourism, particularly based on "beach and sun", is an outstanding example in the Balearic Islands, contributing approximately 45% to the gross domestic product and 32% of jobs in the prepandemic scenario²³.

We hypothesize that (i) due to the long pelagic lifespan of these species, large-scale persistent circulation features, typical of semiclosed basins, may play a prominent role in dispersal and connectivity in combination with stochastic processes. If this is true, (ii) dispersal patterns might be shared among different species with LPDP yielding synchronous interannual patterns. In addition, (iii) this physical-biological coupling may be intricately linked with the atmospheric forcing in the North Atlantic, acknowledged as a major driver of thermohaline circulation in the western Mediterranean²⁴,²⁵. The confirmation of these hypotheses might open opportunities for basin-scale conservation and management policies in this basin and might encourage analogous analyses in other marginal seas and semiclosed basins. To test these hypotheses, we used two model organisms with similar long pelagic lifespans. We combined the most comprehensive dataset of interannual records of both species in the Western Mediterranean shores and the state-of-the-art Lagrangian drifter scheme coupled with a high-resolution hydrodynamic model. We simulated back-in-time trajectories of particles to identify dispersal route pathways and source areas over a decade.

**Biological models.** The European spiny lobster *P. elephas* is a valued, charismatic species that inhabits shallow rocky shores typical in Mediterranean ecosystems. By virtue of its long pelagic phase or phyllosomal larvae, ca. 8 months, compared to other meroplanktonic invertebrates, this species offers an excellent framework to assess mechanisms shaping long-term dispersal. The settlement of early benthic juveniles occurs in shallow rocky areas around late spring-early summer in this region, from May to August²⁶,²⁷.

The mauve stinger jellyfish is the most abundant holoplanktonic scyphozoan in the Mediterranean Sea, the main stinging species and one of the most worrying jellyfish due to its impact on human leisure and economic activities²¹. This species lives ca. 9 months²⁸, being mainly associated with offshore nutrient-rich waters²⁹, i.e., marine fronts and stranding events are closely linked with sea current variations³⁰,³¹. Outbreaks of this species are associated with warm and dry late spring-early summers in the western Mediterranean³², although recent investigations have shown an enhanced bloom frequency and a permanent presence of *P. noctiluca* in the northernmost basin areas³³,³⁴.

**Physical oceanography.** The Mediterranean Sea is a semiclosed basin and is considered a small-scale ocean. The western Mediterranean sub-basin is subject to well-identified circulation and production patterns, largely differing from the eastern sub-basin, mainly due to topographic constraints³⁵. The excess evaporation in the basin is compensated with a continuous inflow of fresher Atlantic Waters (AW) through the Strait of Gibraltar (Fig. 1a). AW flows along the marginal southern and eastern continental slopes following an anticlockwise circuit, and AW is progressively transformed by atmospheric interaction into saltier and colder modified Atlantic water (MAW), particularly in the northernmost basin. From the Ligurian Sea toward the Balearic Sea, the Northern Current (NC) is the dominant geostrophic feature that closes the circulation cell at surface waters. The Gulf of Lions (GoL) is one of the few ocean regions where open-sea deep convection leads to deep or intermediate water formation, mainly occurring every winter. These phenomena are forced by ocean-atmosphere exchanges and persistent and cold N-NW winds²⁵. The Balearic Sea is considered a particular water mass, sharing similar dynamics and features with other sub-basins, such as South China or the Caribbean Seas³⁶. Through the Balearic Channels, there is a meridional exchange between new AW and colder, saltier, and productive MAW transported by the NC, forming an oceanographic density front, the Balearic front. Frontal emplacement is mainly driven by the blocking or exchange of water masses through Balearic channels. The meridional position of this front and channel exchanges are crucial to the marine ecosystem's derived response around the islands³⁶,³⁷.

## Methods

**Simulation specifications and Lagrangian scheme.** We used a state-of-the-art 3D regional climate simulation as the background hydrodynamic model, where a Lagrangian model is coupled to simulate drifting patterns and sources. Specifically, we used the Nucleus for European Modeling of the Ocean NEMO-MED36v75 (version 3.2) covering the whole Mediterranean basin for the period 2003-2013, using a horizontal resolution of 1/36° (≈2–3 km) and 75 z-vertical levels. For the simulations set carried out here, only the first 300 m were considered, and the domain was restricted to the western Mediterranean (Fig. 1). NEMO-MED36v75 successfully reproduces the main circulation features in the Mediterranean basin (surface circulation, coastal current shape and velocities, and deep convection) and mesoscale dynamics in the NW Mediterranean³⁸⁻⁴⁰.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Average sea surface circulation in the western Mediterranean and biological time series locations. The lobster settlement index (LSI, red circles) belongs to different marine protected areas (MPAs) at two sites, one in the northernmost basin (NCat) and three in the Balearic Islands (North, NW and South Mallorca, green area and figure inset) (a). The jellyfish stinging index (JSI, blue circles) covered three tourist hot-spot locations in the Balearic Islands with different orientations and exposed to different water masses (E Ibiza, N Mallorca and S Menorca). Squares denote FAO Geographical Statistical SubAreas (GSAs) used to interpret the connectivity of LSI lobster series (see "Methods"). The Gulf of Lions (GoL, light orange area) and Balearic Sea (BS, light green area) indicate the areas where ocean-atmosphere fluxes were extracted for comparison against field records. A 10-year climatology of drifter sources (%) reaching either NCat (b) or Mallorca (c) is displayed as the particle location eight months before the observations (origin). Source geographic locations were aggregated in monthly density maps of 0.25 × 0.25°. NC, Northern Current in red. LS, Ligurian Sea. SoG, Strait of Gibraltar. The figure and superimposed maps were created using "M Map: A mapping package for MATLAB", v.1.4m (www.eoas.ubc.ca/~rich/map.html) in MATLAB v.R2010b (www.mathworks.com).</figcaption></figure>

A spatially explicit individual-based model was used to simulate the dispersal and arrival of individuals along the western Mediterranean coasts. The model was forced with the 3D hydrodynamic model daily velocity fields using a customized version of the open source modeling tool ICHTHYOP⁴¹. As detailed biological information on the movement of these species through all their planktonic phase is lacking, passive Lagrangian particles (without buoyancy) were subjected to advective and diffusive processes that conditioned their horizontal and vertical movement in the water column.

To establish the appropriate number of virtual individuals to be released in a single simulation, multiple trials were conducted in which the number of released individuals were progressively increased (5000, 10,000, 15,000, 20,000, 50,000, 100,000, 125,000, and 150,000). The average and standard deviation of the results were calculated for each trial, and the point at which these statistics stabilized was determined⁴². Through this process, it was determined that a release of at least 100,000 individuals per simulation was appropriate (p < 0.001). Therefore,
the simulations were set up with 134,000 individuals starting within the two closest cells to the shoreline (≈ 5 km) with a randomly distributed vertical distribution throughout the Western Mediterranean basin. It was assumed that further repetition was not necessary to obtain unbiased results for each simulation. Backward-in-time trajectories of the drifters were simulated to assess the dispersal route and source areas. Each simulation covered 240 days (8 months) following the pelagic duration of the chosen species here²⁸,⁴³. The simulation set included 26 simulations per year, starting every year from the same odd weeks and covering 10 years, from 2004 to 2013 (260 simulations). Furthermore, the initial positions of the individuals were always the same in each simulation; thus, we ran 260 trajectories of the same individual.

Individuals were filtered from initial positions following available field records of LPDP, with ca. 1300 at North Catalonia (NCat) and 3000 (Mallorca) trajectories per location and simulation (Fig. 1a). Subsequently, simulation outputs were analyzed in two different ways:

1. Density maps of each site were computed summing all source positions of individuals using a grid of 1/4° of horizontal resolution. Different simulations were gathered to reconstruct source monthly maps over the 10 years (120 maps) to increase the number of simulated individuals. Using a reference date at the 15th day of each month, we used moving 3-month windows (reference date ± 1.5 months) gathering 6 simulations for each density map, meeting a balance to include enough individuals/simulations but not reaching an excessive smoothing of oceanic processes. Different time windows were tested, obtaining similar results. Then, an empirical orthogonal function (EOF) analysis⁴⁴ was applied to the set of density maps. We used the first EOF and principal components (PC) that explained the highest percentage of variability to extract the main spatiotemporal pattern of the source areas.

2. From each simulation, regarding connectivity for lobster settlement, we filtered the drifters at sources falling within lobster adult habitats (shallowest 200 m)⁴³ and aggregated them using the geographic subareas (GSAs) of the FAO General Fisheries Commission for the Mediterranean. We computed and averaged the percentage of individuals from each GSA (Fig. 1a) during the lobster settlement season (May to August)²⁶. In the absence of more complex biological data, we hypothesized that key atmosphere-ocean fluxes should explain a significant fraction of the simulation results, including their interannual variability. Subsequently, we tested these proposed driving mechanisms against the best available field records of LPDP.

The simulation was originally forced by ARPERA atmospheric fields⁴⁵ but in this case, we used the fluxes from the longer ERA5 atmospheric reanalysis to extend the analysis out of the simulation period, firstly checking the similar behaviour of both atmospheric reanalysis⁴⁶, considering positive fluxes from the ocean toward the atmosphere. We extracted the monthly averaged fluxes (2001–2018) of (i) the ocean heat loss, using here the average during winter months (JFM) from the GoL area²⁵ (Fig. 1a, light orange) and (ii) evaporation-precipitation budget in spring (April) from the Balearic Sea (Fig. 1a, light green).

The strength of the link between the atmosphere-ocean fluxes and main dispersal patterns via the first PC (PC1) was also assessed through robust correlation analysis. We quantified the probability density distribution of Pearson correlation coefficients obtained by bootstrap resampling. This analysis involved random pairwise sampling with replacement, where each time series was resampled 10,000 times. The number of elements in each bootstrap sample equals the number of elements in the original dataset. Detailed results of this analysis can be found in the SI (Supplementary Table S2).

**Lobster Settlement Index.** The Lobster Settlement Index (LSI) was assessed at two main sites, including one MPA at NCat and three MPAs at Balearic Island (north, west and south of Mallorca) (Fig. 1a). Early benthic juveniles were recorded once a year during summer from 2001 to 2018, and more specific details can be found in reference²⁷.

No live animals were used in this study; accordingly, ethical approval was not required.

**Jellyfish stinging index.** We identified summers with a high abundance of jellyfish arriving at the Balearic shores. Here, we proposed an index based on the number of first aids given to bathers due to jellyfish stings as a proxy of arriving events of *P. noctiluca*, the main stinging species in this basin²¹,⁴⁷. To this aim, we used a database originally developed by the Regional Government Emergency Area of the Balearic Islands. Lifeguard daily reports of first aid and causes in beach locations were compiled in the database. We identified and extracted three coherent time series without gaps corresponding to three beach sites with high potential jellyfish impact on the bathers, one on each major island of the Balearic archipelago: Ibiza, Mallorca and Menorca (Fig. 1a). Each time series comprised the number of first aids, including jellyfish stings aggregated over the tourism peak season (July August), covering the 2009–2018 period. To avoid biases, the number of stings was normalized over the number of first aids excluding the stings as a proxy of the number of bathers at each beach. Furthermore, we computed each site stinging series annual anomaly to avoid bias due to local factors.

**Consent to participate.** The authors voluntarily agree to participate in this research study.

## Results and discussion

**Sources and dispersal routes in the NW Mediterranean.** We show that NCat and Mallorca shores shared source areas in the northernmost basin, following the NC path from the Ligurian Sea (Fig. 1b,c). However, Mallorca shore sources also extended along the entire southwestern basin from the Strait of Gibraltar
(Fig. 1c). The northeastern GSAs (GSAs 8, 9, 10 and 11.2) showed a large match and tight correlations with the observed LSI at both sites, although this area contributed to the replenishment of NCat and Mallorca only at low percentages, with peaks ca. 10% and 1.5%, respectively. Other GSAs showed a higher potential contribution as source areas (Fig. 1b,c; Supplementary Fig. S1), although they showed low or even negative correlations with LSI (Supplementary Fig. S1).

Our simulations unveiled that the Ligurian Sea acts as a source (origin) of drifters reaching NW Mediterranean shores via NC, also shaping drifter routes whose pelagic lifespan is shorter than 8 months (Fig. 2b,d). Therefore, the Ligurian Sea may act as a “pelagic nursery area” sensu Kough et al.⁴⁸ due to physical processes known to trap plankters long enough to foster productivity for high trophic levels, including marine mammals, seabirds and fish⁴⁹. This partly explains why this area has been identified as a biodiversity hot spot in the Mediterranean and declared the Pelagos Sanctuary for Mediterranean marine mammals⁵⁰.

**Ocean-atmosphere coupling mold dispersal and connectivity.** Source areas feeding NCat and Mallorca shores are closely linked with the ocean-atmosphere system playing out in the North Atlantic, the influence of which shapes interannual hydrographic variations²⁴ and dispersal patterns in the western Mediterranean. For NCat, a marked seasonal modulation in the source area was captured. During spring and summer (from May to August), variability in source areas was mainly ascribed to the Ligurian Sea and to both current branches encircling Corsica (GSA8, Fig. 1a) feeding the NC (Fig. 3a,c in red). At interannual scales, principal component 1 (PC1, temporal pattern) for spring and summer mainly portrayed the ocean heat loss winter variability in GoL (Fig. 4a). Harsh and windy winters yield intense ocean heat loss around the GoL, promoting deep convection and intermediate or deep-water formation. These processes enhance the strength of the NC²⁵,⁵¹, ultimately increasing the contribution of the northeastern shores as source areas (Figs. 2a,b and 3; Supplementary Fig. S2). Surprisingly, we found that these conditions, in particular the intense heat loss, were concurrent with major lobster settlement peaks in NCat, as shown by the close correlation between winter heat loss and LSI (R² = 0.33; p-value = 0.020; *n* = 16), whereas lower settlement events occurred after warm winters, when heat loss declined (Fig. 4b).

Our simulations uncovered dipole sources and dispersal routes feeding Mallorca shores, involving (i) the circulation cell of MAW encompassing NC, the Balearic Sea and the Ligurian Sea (red) and (ii) the inflowing fresher southwestern AW (blue) (Fig. 3d). A marked interannual signal (Fig. 3a) suggests a seesaw-like effect promoting an alternance in main source areas, ca. 2-3 years. This pattern mirrors the water balance (evaporation-precipitation, E-P) in spring in the Balearic Sea (Fig. 4c). An excess of evaporation (E > P) in this sub-basin drives the northward progress of fresher AW surpassing the Balearic channels. Thus, AW acts as the main source of arriving drifters in the Balearic Islands (negative PC1 and spatial values, blue in Fig. 3c,d) while blocking the arrival of MAW and constraining the contribution of northern sources (Supplementary Fig. S2). In contrast, rainy springs (P > E) increased the influence of productive MAW and northern source areas (positive PC1 and spatial values, red in Fig. 3c,d; Supplementary Fig. S2), ultimately resulting in marked peaks of both LPDS arrivals, as suggested

<figure><img src="image_5.png" /></figure>

**Figure 2.** Spatiotemporal connectivity patterns and dispersal kernels. Comparison of observed LSI anomalies for NCat and Mallorca (red line, Panels a and c, errors are 1 SD) vs. percentage of individuals arriving from adult lobster habitat GSAs (0–200 m) (blue bars) (Fig. 1). The percentage of individuals was averaged over the settlement season (May to August) to compare with settlement field records. Correlations LSI vs. total % individuals from GSAs (8, 9, 10, 11.2): NCat (*r* = 0.89; p-value = 0.0003; *n* = 10); Mallorca (*r* = 0.81; p-value = 0.004; *n* = 9). Dispersal kernels were calculated for drifting times spanning from 1 to 9 months (panels b, d).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Prevailing sources and interannual patterns of LPDP arriving at the target areas after eight months. Temporal (**a**,**c**) and spatial (**b**,**d**) representation of the first principal component (PC1) of empirical orthogonal function (EOF) analysis conducted on source areas of individuals based on 120 monthly origin maps (2003–2014) (more details in “Methods”). The results are for N. catalonia (**a**,**b**) and Mallorca (**c**,**d**). The total explained variance was 29% for N. Catalonia and 44% for Mallorca. Positive PC and EOF values are in red, and negative values are in blue. Original source maps can be reconstructed using the climatology (Fig. 1b,c) and adding the product of PC1 by EOF values. Therefore, intensified influence as source areas emerges when blue/red areas (from EOF value) and PC value match (blue/red). The figure and superimposed maps were created using “M Map: A mapping package for MATLAB”, v.1.4m (www.eoas.ubc.ca/~rich/map.html) in MATLAB v.R2010b (www.mathworks.com).</figcaption></figure>

by LSI and JSI field records (Fig. 4d). We showed evidence of a synchronous interannual pattern involving different LPDS (LSI vs. JSI, $R^2 = 0.65$; p-value = 0.005, *n* = 10). This synchrony in coastal arrival, explained largely through physical forces, supports the observed association between stinging jellyfish and lobster phyllosoma, which may use the umbrella of the jellyfish as a shelter, food source and transport vector⁵².

**A mechanistic explanation of the interannual variability in western Mediterranean metapopulations.** Our approach allowed us to mechanistically test our hypotheses, first unveiling the pervasive role of Ligurian Sea-NC and AW as prevailing dispersal routes of the LPDP in the W Mediterranean Sea (Figs. 1 and 3). Second, the dispersal patterns of different species shared with the LPDP in the W Mediterranean were further modulated at the interannual scale by atmospheric-ocean coupling, but the latter responded to a site-specific atmospheric driver (Fig. 4). This approach opens opportunities for understanding the processes and scales shaping metapopulation structure and dynamics and hints to test similar hypotheses in other marginal seas and semienclosed basins to boost conservation-management policies for marine resources.

The bloom frequency of *P. noctiluca* in the western Mediterranean apparently has increased along with environmental changes in the region experienced in recent decades³³. However, little synoptic information exists on its population structure and dynamics in this basin²¹,³¹, although this information is essential for jellyfish management actions, particularly where economic activities are threatened. Current knowledge suggests that a single mauve stinger jellyfish population exists in the whole basin⁵³ and that interannual variability in stranding is partly coupled with climate³²,³³ and stochastic variations³⁰,⁵⁴. Short-term forecasting can effectively simulate local jellyfish stranding events⁵⁴,⁵⁵; however, we show that a comprehensive description and potential forecast of jellyfish blooms at larger scales cannot be achieved without fully resolving source areas and dispersal routes during the whole lifespan of the species⁵⁶. The explanatory mechanisms proposed here open the possibility for potential forecasting of swarm dispersal and stranding events, stressing the potential use of North Atlantic climate variability for the W Mediterranean (Fig. 4). Furthermore, these findings enable a future design of research efforts, e.g., description of processes at retention areas, to disentangle the relative role of physical and biological components on the interannual variability of these blooms.
<figure><img src="image_7.png" /></figure>

**Figure 4.** Atmosphere-ocean coupling driving dispersal patterns and proposed mechanisms explaining LPDP interannual variability. Density distribution of the Pearson correlation coefficients between ocean-atmosphere fluxes and PC1 (Fig. 3a,c) for summers/settlement season of jellyfish and lobster, which was obtained by bootstrap resampling (10,000 times) for the two sites: (a) Heat loss (HL) vs. PC1 for North Catalonia (NCat) (Fig. 3a). (c) (Evaporation-precipitation, E–P) vs. PC1 for Mallorca (Fig. 3c). Linear relationship of the atmosphere-ocean fluxes and LPDP field records for both sites, Lobster Settlement Index (LSI) and Jellyfish Stinging Index (JSI). (b) For NCat, HL vs. LSI ($R^2$ = 0.33; p-value = 0.021, *n* = 16). (d) In Mallorca, (E–P) vs. LSI ($R^2$ = 0.33; p-value = 0.003, *n* = 14) and (E-P) vs. JSI ($R^2$ = 0.66; p-value = 0.005, *n* = 10). More details in Supplementary Table S2, Supplementary Information). We found synchronic patterns for both LPDP: LSI vs. JSI ($R^2$ = 0.65; p-value = 0.005, *n* = 10).

Our findings on spiny lobster settlement offer a mechanistic explanation of its interannual variability and emphasize that fisheries management areas and political borders do not fit the metapopulation dynamics of species with the LPDP, particularly in marginal seas with several riverine countries⁴⁸,⁵⁷. Hence, these results are essential for tailoring the managing actions of spiny lobster fisheries in the Western Mediterranean. Genetic studies suggest a common lobster larval pool for the northernmost basin in the NW Mediterranean, which is supported by our simulations that further unveiled a potential influx of individuals from southern or Atlantic sources feeding the Balearic Islands population⁵⁸,⁵⁹. Our results bear vast implications for conservation and resource management policies, as the connectivity patterns uncovered by our simulations showed that on the NW Spanish coast, MPAs might benefit from spiny lobster population recovery in the Sardinia area (GSA 11.2 in Figs. 1a and 2a, c; Supplementary Figs. S1, S2)⁶⁰, thus stressing the need for wide-scale multinational efforts to develop effective resource management and conservation policies. Ignoring source areas will jeopardize the sustainable management of the already threatened spiny lobster population of the Mediterranean Sea.

Our results on the physical influence on LPDP dispersal patterns set a baseline to connect physical forcing with future climate trends. We demonstrated that winter conditions in the GoL area are crucial for LPDP dispersal and connectivity, thus explaining interannual variations. Projections for 2030⁶¹ suggest oceanographic shifts leading to an increase in (i) the strength and duration of deep convection and (ii) an enhancement of mesoscale activity and NC circulation due to meridional wind component intensification. These changes may affect the metapopulation features of species with LPDP in the northwestern Mediterranean and therefore should be considered in management scenarios.

## Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Received: 30 September 2022; Accepted: 6 February 2023

Published online: 09 February 2023

## References

1. Guichard, F., Levin, S. A., Hastings, A. & Siegel, D. Toward a dynamic metacommunity approach to marine reserve theory. *BioScience* 54(11), 1003. https://doi.org/10.1641/0006-3568(2004)054[1003:tadmat]2.0.co;2 (2004).
2. Wieters, E. A., Gaines, S. D., Navarrete, S. A., Blanchette, C. A. & Menge, B. A. Scales of dispersal and the biography of marine predator-prey interactions. *Am. Nat.* **171**(3), 405–417. https://doi.org/10.1086/527492 (2008).

3. Martinez-Moreno, J. et al. Global changes in oceanic mesoscale currents over the satellite altimetry record. *Nat. Clim. Change* https://doi.org/10.1038/s41558-021-01006-9 (2021).

4. van Gennip, S. J. et al. Going with the flow: The role of ocean circulation in global marine ecosystems under a changing climate. *Glob. Change Biol.* **23**(7), 2602–2617. https://doi.org/10.1111/gcb.13586 (2017).

5. O'Connor, M. I. et al. Temperature control of larval dispersal and the implications for marine ecology, evolution, and conservation. *Proc. Natl. Acad. Sci. U.S.A.* **104**(4), 1266–1271. https://doi.org/10.1073/pnas.0603422104 (2007).

6. Cowen, R. K. & Sponaugle, S. Larval dispersal and marine population connectivity. *Ann. Rev. Mar. Sci.* **1**(1), 443–466. https://doi.org/10.1146/annurev.marine.010908.163757 (2009).

7. Ospina-Alvarez, A., Parada, C. & Palomera, I. Vertical migration effects on the dispersion and recruitment of European anchovy larvae: From spawning to nursery areas. *Ecol. Model.* **231**, 65–79. https://doi.org/10.1016/j.ecolmodel.2012.02.001 (2012).

8. Selkoe, K. A. & Toonen, R. J. Marine connectivity: A new look at pelagic larval duration and genetic metrics of dispersal. *Mar. Ecol. Prog. Ser.* **436**, 291–305. https://doi.org/10.3354/meps09238 (2011).

9. Siegel, D. A. et al. The stochastic nature of larval connectivity among nearshore marine populations. *Proc. Natl. Acad. Sci. U.S.A.* **105**(26), 8974–8979. https://doi.org/10.1073/pnas.0802544105 (2008).

10. De Lestang, S. et al. What caused seven consecutive years of low puerulus settlement in the western rock lobster fishery of Western Australia?. *ICES J. Mar. Sci.* **72**, i49–i58. https://doi.org/10.1093/icesjms/fsu177 (2015).

11. Linnane, A. et al. Evidence of large-scale spatial declines in recruitment patterns of southern rock lobster *Jasus edwardsii*, across south-eastern Australia. *Fish. Res.* **105**(3), 163–171. https://doi.org/10.1016/j.fishres.2010.04.001 (2010).

12. Briones-Fourzán, P., Candela, J. & Lozano-Alvarez, E. Postlarval settlement of the spiny lobster Panulirus argus along the Caribbean coast of Mexico: Patterns, influence of physical factors, and possible sources of origin. *Limnol. Oceanogr.* **53**(3), 970–985. https://doi.org/10.4319/lo.2008.53.3.0970 (2008).

13. Haury, L. R., McGowan, J. A. & Wiebe, P. H. Patterns and processes in the time-space scales of plankton distributions. In *Spatial Pattern in Plankton Communities* (ed. Steele, J. H.) 277–327 (Springer US, 1978). https://doi.org/10.1007/978-1-4899-2195-6_12.

14. Cowen, R. K., Paris, C. B. & Srinivasan, A. Scaling of connectivity in marine populations. *Science* **311**(5760), 522–527. https://doi.org/10.1126/science.1122039 (2006).

15. Kavanaugh, M. T. et al. Seascapes as a new vernacular for pelagic ocean monitoring, management and conservation. *ICES J. Mar. Sci.* **73**(7), 1839–1850. https://doi.org/10.1093/icesjms/fsw086 (2016).

16. Ospina-Alvarez, A., Weidberg, N., Aiken, C. M. & Navarrete, S. A. Larval transport in the upwelling ecosystem of central Chile: The effects of vertical migration, developmental time and coastal topography on recruitment. *Prog. Oceanogr.* **168**, 82–99. https://doi.org/10.1016/j.pocean.2018.09.016 (2018) http://www.sciencedirect.com/science/article/pii/S0079661117300800.

17. Palumbi, S. Population genetics, demographic connectivity, and the design of marine reserves. *Ecol. Appl.* **13**(1 Supplement), S146–S158 (2003).

18. Barahona, M. et al. Environmental and demographic factors influence the spatial genetic structure of an intertidal barnacle in central-northern Chile. *Mar. Ecol. Prog. Ser.* **612**, 151–165. https://doi.org/10.3354/meps12855 (2019) http://www.int-res.com/abstracts/meps/v612/p151-165/.

19. Spanier, E. et al. A concise review of lobster utilization by worldwide human populations from prehistory to the modern era. *ICES J. Mar. Sci.* **72**(May), i7–i21. https://doi.org/10.1093/icesjms/fsv066 (2015).

20. IUCN. *Palinurus elephas*: Goñi, R.: The IUCN Red List of Threatened Species 2014: e.T169975A1281221. Tech. Rep., International Union for Conservation of Nature (2013). http://www.iucnredlist.org/details/169975/0. Type: dataset.

21. Canepa, A. et al. *Pelagia noctiluca* in the Mediterranean sea (eds Pitt, K. A. & Lucas, C. H.) In *Jellyfish Blooms*, Vol. 9789400770 237–266 (Springer Netherlands, 2014). https://doi.org/10.1007/978-94-007-7015-7_11.

22. Bosch-Belmar, M. et al. Jellyfish blooms perception in Mediterranean finfish aquaculture. *Mar. Policy* **76**, 1–7. https://doi.org/10.1016/j.marpol.2016.11.005 (2017).

23. Exceltur. Impactur baleares 2014. Tech. Rep., EXCELTUR - Govern de les Illes Balears, Madrid (2014).

24. Vignudelli, S., Gasparini, G. P., Astraldi, M. & Schiano, M. E. A possible influence of the North Atlantic Oscillation on the circulation of the Western Mediterranean Sea. *Geophys. Res. Lett.* **26**(5), 623–626. https://doi.org/10.1029/1999GL900038 (1999).

25. Somot, S. et al. Characterizing, modelling and understanding the climate variability of the deep water formation in the North-Western Mediterranean Sea. *Clim. Dyn.* **51**(3), 1179–1210. https://doi.org/10.1007/s00382-016-3295-0 (2018).

26. Díaz, D., Marí, M., Abelló, P. & Demestre, M. Settlement and juvenile habitat of the European spiny lobster *Palinurus elephas* (Crustacea: Decapoda: Palinuridae) in the Western Mediterranean Sea. *Sci. Mar.* **65**(4), 347–356. https://doi.org/10.3989/scimar.2001.65n4347 (2001).

27. Muñoz, A. et al. Exploration of the inter-annual variability and multi-scale environmental drivers of European spiny lobster, *Palinurus elephas* (Decapoda: Palinuridae) settlement in the NW Mediterranean. *Mar. Ecol*. https://doi.org/10.1111/maec.12654 (2021).

28. Malej, A. & Malej, M. Population dynamics of the jellyfish *Pelagia noctiluca* (Forsskal, 1775) In *Marine Eutrophication and Population Dynamics* (eds Colombo, G., Ferrari, I., V., C. & R., R.) 215–219 (Olsen and Olsen, 1992).

29. Ottmann, D. et al. Abundance of *Pelagia noctiluca* early life stages in the western Mediterranean Sea scales with surface chlorophyll. *Mar. Ecol. Prog. Ser.* **658**, 75–88. https://doi.org/10.3354/meps13423 (2021).

30. Benedetti-Cecchi, L. et al. Deterministic factors overwhelm stochastic environmental fluctuations as drivers of jellyfish outbreaks. *PLoS One* **10**(10), 1–16. https://doi.org/10.1371/journal.pone.0141060 (2015).

31. Licandro, P. et al. A blooming jellyfish in the northeast Atlantic and Mediterranean. *Biol. Lett.* **6**(5), 688–691. https://doi.org/10.1098/rsbl.2010.0150 (2010).

32. Goy, J., Morand, P. & Etienne, M. Long-term fluctuations of *Pelagia noctiluca* (Cnidaria, Scyphomedusa) in the western Mediterranean Sea. Prediction by climatic variables. Deep Sea Res. Part A *Oceanogr. Res. Pap.* **36**(2), 269–279 (1989). https://doi.org/10.1016/0198-0149(89)90138-6.

33. Yahia, M. N. D. et al. Are the outbreaks timing of *Pelagia noctiluca* (Forsskal, 1775) getting more frequent in the Mediterranean basin?. *ICES Cooper. Res. Rep.* **300**, 8–14 (2010).

34. Ferraris, M. et al. Distribution of *Pelagia noctiluca* (Cnidaria, Scyphozoa) in the Ligurian Sea (NW Mediterranean Sea). *J. Plankton Res.* **34**(10), 874–885. https://doi.org/10.1093/plankt/fbs049 (2012).

35. Millot, C. Circulation in the Western Mediterranean Sea. *J. Mar. Syst.* **20**(1–4), 423–442. https://doi.org/10.1016/S0924-7963(98)00078-5 (1999).

36. Galarza, J. A. et al. The influence of oceanographic fronts and early-life-history traits on connectivity among littoral fish species. *Proc. Natl. Acad. Sci.* **106**(5), 1473–1478. https://doi.org/10.1073/pnas.0806804106 (2009).

37. Fernández de Puelles, M. L. & Molinero, J. C. Decadal changes in hydrographic and ecological time-series in the Balearic Sea (western Mediterranean), identifying links between climate and zooplankton. *ICES J. Mar. Sci.* **65**(3), 311–317. https://doi.org/10.1093/icesjms/fsn017 (2008).

38. Arsouze, T. et al. CIESM (ed.) Sensibility analysis of the Western Mediterranean Transition inferred by four companion simulations. (ed. CIESM) EGU General Assembly Conference Abstracts, Vol. 1 of EGU General Assembly Conference Abstracts, 13073 (2013).
39. Amores, A., Jordà, G., Arsouze, T. & Le Sommer, J. Up to what extent can we characterize ocean eddies using present-day gridded altimetric products?. *J. Geophys. Res. Oceans* **123**(10), 7220–7236. https://doi.org/10.1029/2018JC014140 (2018).

40. Waldman, R. et al. Impact of the mesoscale dynamics on ocean deep convection: The 2012–2013 case study in the northwestern mediterranean sea. *J. Geophys. Res. Oceans* **122**(11), 8813–8840. https://doi.org/10.1002/2016JC012587 (2017).

41. Lett, C. et al. A Lagrangian tool for modelling ichthyoplankton dynamics. *Environ. Model. Softw.* **23**(9), 1210–1214. https://doi.org/10.1016/j.envsoft.2008.02.005 (2008).

42. Brickman, D. & Smith, P. C. Lagrangian stochastic modeling in coastal oceanography. *J. Atmos. Ocean. Technol.* **19**(1), 83–99. https://doi.org/10.1175/1520-0426(2002)019<0083:LSMICO>2.0.CO;2 (2002).

43. Goñi, R. & Latrouite, D. Review of the biology, ecology and fisheries of Palinurus spp. species of European waters: *Palinurus elephas* (Fabricius, 1787) and *Palinurus mauritanicus* (Gruvel, 1911). *Cahiers de Biol. Mar.* **46**(2), 127–142 (2005).

44. Bjornsson, H. & Venegas, S. A manual for EOF and SVD analyses of climatic data. Tech. Rep. CCGCR Report No. 97-1, McGill's Centre for Climate and Global Change Research (C2GCR) (1997).

45. Herrmann, M., Somot, S., Sevault, F., Estournel, C. & Déqué, M. Modeling the deep convection in the northwestern mediterranean sea using an eddy-permitting and an eddy-resolving model: Case study of winter 1986–1987. *J. Geophys. Res. Oceans* **113**(C4) (2008). https://doi.org/10.1029/2006JC003991.

46. Hersbach, H. et al. ERA5 monthly averaged data on single levels from 1979 to present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). **10**, 252–266 (2019). https://doi.org/10.24381/cds.f17050d7.

47. Bernard, P., Berline, L. & Gorsky, G. Long term (1981–2008) monitoring of the jellyfish *Pelagia noctiluca* (Cnidaria, Scyphozoa) on Mediterranean Coasts (Principality of Monaco and French Riviera). *J. Oceanogr. Res. Data* **4**(1), 1–10 (2011).

48. Kough, A. S., Paris, C. B. & Butler, M. J. IV. Larval connectivity and the international management of fisheries. *PLoS One* **8**(6), 1–12. https://doi.org/10.1371/journal.pone.0064970 (2013).

49. Sandvik, H. et al. Modelled drift patterns of fish larvae link coastal morphology to seabird colony distribution. *Nat. Commun.* **7**(May), 1–8. https://doi.org/10.1038/ncomms11599 (2016).

50. Notarbartolo-Di-Sciara, G., Agardy, T., Hyrenbach, D., Scovazzi, T. & Van Klaveren, P. The Pelagos Sanctuary for Mediterranean marine mammals. *Aquat. Conserv. Mar. Freshw. Ecosyst.* **18**(4), 367–391. https://doi.org/10.1002/aqc.855 (2008).

51. Astraldi, M., Gasparini, G. P., Vetrano, a. & Vignudelli, S. Hydrographic characteristics and interannual variability of water masses in the central Mediterranean: A sensitivity test for long-term changes in the Mediterranean Sea. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **49**(4), 661–680 (2002). https://doi.org/10.1016/S0967-0637(01)00059-0.

52. Muffett, K. & Miglietta, M. P. Planktonic associations between medusae (classes Scyphozoa and Hydrozoa) and epifaunal crustaceans. *PeerJ* **9**, e11281. https://doi.org/10.7717/peerj.11281 (2021) https://peerj.com/articles/11281.

53. Stopar, K., Ramšak, A., Trontelj, P. & Malej, A. Lack of genetic structure in the jellyfish *Pelagia noctiluca* (Cnidaria: Scyphozoa: Semaeostomeae) across European seas. *Mol. Phylogenet. Evol.* **57**(1), 417–428. https://doi.org/10.1016/j.ympev.2010.07.004 (2010).

54. Berline, L., Zakardjian, B., Molcard, A., Ourmières, Y. & Guihou, K. Modeling jellyfish *Pelagia noctiluca* transport and stranding in the Ligurian Sea. *Mar. Pollut. Bull.* **70**(1–2), 90–99. https://doi.org/10.1016/j.marpolbul.2013.02.016 (2013).

55. Prieto, L., Macías, D., Peliz, A. & Ruiz, J. Portuguese Man-of-War (*Physalia physalis*) in the Mediterranean: A permanent invasion or a casual appearance? *Sci. Rep.* **5** (2015). https://doi.org/10.1038/srep11545.

56. Houghton, J. D. R. et al. Identification of genetically and oceanographically distinct blooms of jellyfish. *J. R. Soc. Interface* **10**(80), 20120920–20120920. https://doi.org/10.1098/rsif.2012.0920 (2013).

57. Segura-García, I. et al. Reconstruction of larval origins based on genetic relatedness and biophysical modeling. *Sci. Rep.* **9**(1), 1–9. https://doi.org/10.1038/s41598-019-43435-9 (2019).

58. Elphie, H., Raquel, G., David, D. & Serge, P. Detecting immigrants in a highly genetically homogeneous spiny lobster population (*Palinurus elephas*) in the northwest Mediterranean Sea. *Ecol. Evol.* **2**(10), 2387–2396. https://doi.org/10.1002/ece3.349 (2012).

59. Babbucci, M. et al. Population structure, demographic history, and selective processes: Contrasting evidences from mitochondrial and nuclear markers in the European spiny lobster *Palinurus elephas* (Fabricius, 1787). *Mol. Phylogenet. Evol.* **56**(3), 1040–1050. https://doi.org/10.1016/j.ympev.2010.05.014 (2010).

60. Cau, A. et al. European spiny lobster recovery from overfishing enhanced through active restocking in Fully Protected Areas. *Sci. Rep.* **9**(1) (2019). https://doi.org/10.1038/s41598-019-49553-8.

61. Macías, D., Garcia-Gorriz, E. & Stips, A. Deep winter convection and phytoplankton dynamics in the NW Mediterranean Sea under present climate and future (Horizon 2030) scenarios. *Sci. Rep.* **8**(1), 1–15. https://doi.org/10.1038/s41598-018-24965-0 (2018).

## Acknowledgements

ER-R is grateful for the funding from Govern de les Illes Balears-Conselleria d'Innovació, Recerca i Turisme, Programa Vicenç Mut. DD acknowledges MAREMATING project (CTM2016-77027-R). AM acknowledges pre-doctoral FPI Fellowship from Conselleria d'Innovació, Recerca i Turisme of the Government of the Balearic Islands co-financed by the European Social Fund, as part of the FSE 2014-2020 Operational Program. AO-A acknowledges "Juan de la Cierva Incorporación 2020" (IJC2020-044266-I) fellowship funded by MCIN/AEI/10.13039/501100011033 and the European Union "NextGeneration EU/PRTR".

## Author contributions

E.R.-R., A.A., A.O.-A. and I.A.C. designed the research. D.D. and A.M. carried out the field samplings. E.R.-R., A.A., A.O.-A. performed the simulations. E.R.-R., A.A. and J.C.M. analysed the data. All the authors contributed writing the text.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-023-29543-7.

**Correspondence** and requests for materials should be addressed to A.O.-A.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023