Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Mesoscale circulation determines broad spatio-temporal settlement patterns of lobster

  • Paulina Cetina-Heredia ,

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Visualization, Writing – original draft

    p.cetinaheredia@unsw.edu.au

    Affiliation Regional and Coastal Oceanography Laboratory, School of Mathematics and Statistics, UNSW Australia, Sydney, Australia

  • Moninya Roughan,

    Roles Conceptualization, Funding acquisition, Project administration, Writing – review & editing

    Affiliation Regional and Coastal Oceanography Laboratory, School of Mathematics and Statistics, UNSW Australia, Sydney, Australia

  • Geoffrey Liggins,

    Roles Conceptualization, Investigation, Methodology, Writing – review & editing

    Affiliation Department of Primary Industries, NSW Fisheries, Sydney, New South Wales, Australia

  • Melinda A. Coleman,

    Roles Conceptualization, Funding acquisition, Writing – review & editing

    Affiliation Department of Primary Industries, NSW Fisheries and National Marine Science Centre, Coffs Harbour, New South Wales, Australia

  • Andrew Jeffs

    Roles Conceptualization, Funding acquisition, Methodology, Writing – review & editing

    Affiliation Institute of Marine Science, and School of Biological Sciences, University of Auckland, Auckland, New Zealand

Correction

2 Apr 2019: Cetina-Heredia P, Roughan M, Liggins G, Coleman MA, Jeffs A (2019) Correction: Mesoscale circulation determines broad spatio-temporal settlement patterns of lobster. PLOS ONE 14(4): e0214996. https://doi.org/10.1371/journal.pone.0214996 View correction

Abstract

The influence of physical oceanographic processes on the dispersal of larvae is critical for understanding the ecology of species and for anticipating settlement into fisheries to aid long-term sustainable harvest. This study examines the mechanisms by which ocean currents shape larval dispersal and supply to the continental shelf-break, and the extent to which circulation determines settlement patterns using Sagmariasus verreauxi (Eastern Rock Lobster, ERL) as a model species. Despite the large range of factors that can impact larval dispersal, we show that within a Western Boundary Current system, mesoscale circulation explains broad spatio-temporal patterns of observed settlement including inter-annual and decadal variability along 500 km of coastline. To discern links between ocean circulation and settlement, we correlate a unique 21- year dataset of observed lobster settlement (i.e., early juvenile & pueruli abundance), with simulated larval settlement. Simulations use outputs of an eddy-resolving, data-assimilated, hydrodynamic model, incorporating ERL spawning strategy and larval duration. The latitude where the East Australian Current (EAC) deflects east and separates from the continent determines the limit between regions of low and high ERL settlement. We found that years with a persistent EAC flow have low settlement while years when mesoscale eddies prevail have high settlement; in fact, mesoscale eddies facilitate the transport of larvae to the continental shelf-break from offshore. Proxies for settlement based on circulation features observed with satellites could therefore be useful in predicting broadscale patterns of settlement orders of magnitudes to guide harvest limits.

Introduction

The influence of physical oceanographic processes on the dispersal of larvae is critical for understanding the ecology of species, evolution of marine communities, and for predicting settlement into fisheries to aid long-term sustainable harvest [1]. For most species with dispersive larval stages, however, we know little about the specific factors influencing dispersal during their time in the ocean. Despite an increase in studies using biophysical models to explore larval dispersal, the extent to which different factors such as ocean climate, or species larval traits, can explain patterns of settlement remains largely hypothetical. This is mostly due to the paucity of long term observed settlement data required to test the accuracy of biophysical models that seek to represent settlement. Comparison of biophysical model estimates against long-term observations can be a powerful tool for discerning the key factors involved in structuring patterns of settlement which often drives subsequent recruitment to important fisheries [2]. Indeed, the ability to predict the arrival of larvae onshore and subsequent recruitment into a fishery can aid in setting and adapting harvest limits effective for sustainable management [35].

Palinurid lobsters, commonly known as spiny or rock lobsters, are key components of coastal ecosystems [6, 7], and form the basis of highly lucrative fisheries in many parts of the world with a combined global value of harvests in excess of US$1 billion a year [8]. As a result of their high market value, dedicated fisheries, often with intense harvesting pressure and with specific management focus, are common features for most spiny lobster species [9]. Spiny lobster populations worldwide display dramatic inter-annual and spatial variability in settlement and recruitment [3, 10]. The relationship between settlement abundance and stock populations for some lobster species, such as Panulirus argus, distributed across extensive coastlines of the Atlantic, is hindered by post-settlement processes [1113]. However, for other species, direct observations of recruitment to the fishable stock have been related to settlement of post-larvae (pueruli), suggesting that settlement provides a leading indicator of the abundance of lobsters several years later [1417]. For instance, in Western Australia, pueruli settlement of Panulirus cygnus, western rock lobster, is correlated with adult abundance 4 years later, and annual measures of pueruli settlement are used as a fundamental variable in managing subsequent landings [14]. The strength of pueruli settlement has therefore utility in predicting recruitment and calculating harvest limits [15]; consequently, long-term monitoring of lobster settlement is commonplace [14, 18, 19]. An understanding of drivers of pueruli settlement will aid improved prediction of settlement and its variability over large spatial scales, thereby guiding future monitoring programs and management of lobster fisheries.

Sagmariasus verreauxi (Eastern Rock Lobster, ERL) inhabit a stretch of uninterrupted coastline along southeast Australia. In this region, the East Australian Current (EAC), a swift Western Boundary Current (WBC), meanders poleward aligned with the coast year-round, feeding eddies and driving the regional main patterns of circulation [20]. The EAC separates from the continent typically between 30.7–32.4°S but can extend further south to 38°S [21]; the formation and shedding of mesoscale eddies throughout the year displace the EAC separation latitude [21, 22], and influence the regional circulation with a ~120 days periodicity [23, 24]. A previous modelling study of ERL dispersal under contemporary and future climate conditions in the EAC proposed that a future poleward displacement of the mean EAC separation latitude will induce a similar poleward shift of maximum settlement [25]. However, a quantitative validation against observations of the extent to which the EAC determines settlement and its variability along southeast Australia remains unknown.

Here, we quantitatively test whether observed settlement of ERL along southeastern Australia can be estimated from larval dispersal simulations that rely on mesoscale circulation and pelagic larval duration. ERL settlement has been measured for 21 years along this coast, providing an invaluable data set for validation of simulated settlement and for discerning the relative importance of physical processes on pueruli settlement. We use outputs of an eddy-resolving hydrodynamic model with data assimilation to simulate larval dispersal over a 21 year period. The larval trajectories are used to: 1) examine the role of ocean currents in explaining observed spatial and temporal variability of ERL settlement, and 2) diagnose circulation features facilitating larval supply to the continental shelf-break in this dynamic region.

Methods

Study species

Sagmariasus verreauxi (Eastern Rock Lobster, ERL), has a life cycle with benthic adults living in coastal waters and a pelagic larval phase spent in open ocean. Hatching of eggs, carried by adult female ERL occurs simultaneously across spawning grounds (i.e., in southeast Australia between 28°-33°S), from November to January peaking in December (G. Liggins pers. comm.). The hatched naupliosoma larvae moult into phyllosoma larvae and drift in the ocean for 8–12 months before metamorphosing into pueruli that can actively migrate across the shelf seeking suitable habitat in which to settle [26, 27]. The subsequent recruitment of juveniles re-stocks the adult population, which shows no genetic differentiation along southeast Australia [28].

ERL puerulus monitoring program

In 1995 the Department of Primary Industry Fisheries of New South Wales (NSW, 28–37°S) established an annual survey to monitor the abundance of ERL pueruli recruiting to the coast. One of the general objectives was to understand the spatial and temporal variability in recruitment and the influence of environmental factors on settlement [26], [29]. As part of the monitoring program ERL pueruli collectors have been deployed to survey settlement at four locations (three replicate sites within each location with three collectors at each site) extending along the coast of south-east Australia: Coffs Harbour (30.3°S), Tuncurry (32.1°S), Sydney (33.8°S), and Ulladulla (35.4°S) [26], (Fig 1).

thumbnail
Fig 1.

a) Map of the study area showing the settlement monitoring locations Coffs Harbour (CH), Tuncurry (Tu), Sydney (Sy) and Ulladulla (Ul) indicated by the empty circles; the 150 m isobath is shown in grey and the area within which settlement was consider to occur is delimited by the blue line (i.e. 150 m isobath along the southeast coast of Australia and the 50 m isobath along the south coast of Australia). b) Top view map encompassing spawning grounds; the grey dots show a top view of all possible release locations and sections shown in panels c-e are indicated with black thick lines. c-e) vertical sections across 28, 31.1 and 33°S, where the continental shelf is ~55 km, the narrowest (~11 km) and the widest (77 km) respectively, showing the depth distribution of release locations (grey dots) for a day (1 Dec 2010) as an example. Black dots indicate all possible release locations.

https://doi.org/10.1371/journal.pone.0211722.g001

The collectors and their use in the long term pueruli settlement surveys are thoroughly described in [30]; briefly, they are sea-weed-type collectors sampled every 4 weeks during the first quarter of the lunar month from August to January. Monthly settlement at each location is estimated using the mean number of pueruli captured in the collectors.

The ERL monitoring program has produced a 21-year time series of ERL settlement along a latitudinal range extending over 500 km. The number of pueruli settling at each collector has been recorded and used to compute a mean number of pueruli settled at each location and survey month (from August to January within each year). This study examines inter-annual variability in settlement along ~500 km of coastline along southeast Australia using yearly settlement. Yearly settlement is derived from the summation of monthly settlement measured from all pueruli collectors for each of the four locations from August to January. For example, pueruli settlement surveyed between August 1995 and January 1996 is considered to contribute to settlement in 1996.

Hydrodynamic model

ERL larval transport by ocean currents is simulated by advecting larvae with velocity outputs from the Ocean Forecasting Australian Model Bluelink ReANalysis (OFAM BRAN 3p5). This is a near-global hydrodynamic model with data assimilation developed to hindcast and forecast realistic upper ocean conditions. BRAN3p5 has a spatial horizontal resolution of 0.1° for all longitudes between 75°S and 75°N that allows it to resolve eddies, and 51 z* vertical levels with a 5 m vertical resolution in the first 40 m of water depth, 10 m vertical resolution to 200 m deep, and a decreasing vertical resolution spanning 120–1000 m between vertical layers from 200 m to the sea floor [31]. To better represent bottom topography a ‘partial cell’ scheme (z*) is implemented [32].

OFAM BRAN 3p5 assimilates sea surface temperature (SST), sea level anomaly (SLA) and in-situ temperature and salinity observations to regularly adjust the model state to match the observations. Data assimilation is particularly useful in eddy-resolving models because the generation and evolution of eddies is modulated by unpredictable instabilities; therefore data assimilation is useful in order to reproduce eddies in the correct place and time or with accurate intensity and characteristics [31]. The design of BRAN 3p5 allows it to reliably reproduce mesoscale eddies which are known to influence larval dispersal [25, 33], and therefore this model is an ideal tool to study larval dispersal.

ERL larval transport simulations

Larval transport simulations are conducted with the Connectivity Modelling System, CMS, [34], a Lagrangian tracking model. Because the study focuses on the impact of ocean circulation on larval transport and consequent settlement, the Lagrangian simulations represent passive larvae and do not account for spatio-temporal changes in breeding stock abundance and productivity. Nevertheless, particle release, particle advection time, and rules for settlement are based on the ERL habitat, life cycle, and larval phase duration. Specifically, particles are seeded along the natural spawning grounds from 28 to 33°S and between the coastline and the 150 m isobath daily from 15 November until the 15 January; data supplied by NSW commercial lobster fishers indicate that, during the period 2009–2017, >97% of the reported catch of ovigerous females occurred within these latitudes and depths; and data from a NSW fishery-independent survey indicates that 90% of ovigerous females shed their eggs during December (G. Liggins pers. comm.). Larvae are advected with simulated ocean currents throughout their pelagic stage (8–12 months, [26] and all larval trajectories are recorded over 12 months. As a proxy for settlement we model the delivery of larvae onto the continental shelf-break due to the narrow shelf width (as low as 16 km within the study domain) and the model spatial resolution (0.1°). After hatching, ERL larvae take 8–12 months to develop into a pueruli competent to settle. Thus, larvae are considered to settle if at any time 8–12 months after release they are anywhere between the coastline and the 150 m isobath along the east coast, or between the coastline and the 50 m isobath within Bass Strait (Fig 1A); the trajectory of larvae after settlement is ignored. Due to paucity of data on larval survival in the ocean for larvae of any lobster species [35], we assume mortality is spatially and temporally uniform.

The vertical distribution of lobster larvae in the ocean has been examined for some species; the western lobster, Panulirus cygnus, migrates vertically to increasing depths with age [36]; similar changes in the vertical distribution of the southern rock lobster, Jasus edwardsii, have been observed with mid-stage larvae concentrated within the top 20 m and late stage larvae moving down to 100 m [37]. It is expected that larvae of this study species, Sagmariasus verreauxi, undergo vertical migration; however, due to lack of empirical data describing their vertical distribution, the larval transport simulations do not implement specific vertical migration. Instead, they allow for passive vertical movement using 3D velocity fields for advection as well as an initial distribution encompassing all depths. Dispersal studies have shown larval behavior can impact settlement [38]; for instance, [39] finds that ontogenetic vertical migration of the Caribbean lobster, Panulirus argus, larvae diminishes dispersal distances and increases settlement abundance. This is partly due to differences in current direction or speed with depth [40]. However, unlike the ocean region inhabited by P. argus larvae, the dominant circulation features in the study region, such as the EAC and mesoscale eddies, are all coherent with depth; with the EAC prevailing over thousands of meters [41], and mesoscale eddies extending over 100s of meters from the surface which is the primary habitat of these larvae [42, 43]. Thus, main circulation features are expected to induce similar particle trajectories irrespective of their depth distribution; particularly, over the top 200 m of the water column, a range within which the larvae of all rock lobster species have been found to reside regardless of developmental stage ([35]. In the region, [44] argue that surface dispersal patterns of king prawn along southeast Australia would not change significantly from those at depth. By simulating the dispersal of larvae without behavior we test the extent to which advection by ocean currents can explain observed broad patterns of settlement. Implementation of larval behavior however, could improve settlement predictions from simulations if empirical data were available to inform models.

We released the same number of particles at each source-latitude (i.e., every 0.1° between 28–33°S) and at the centre of the grid cells of the hydrodynamic model; the number of particles seeded at each source-latitude is determined by the number of grid cells between the coastline and the 150 m isobath at the latitude where the continental platform is the narrowest (i.e., 13 grid cells at 31.1°S). The release longitudes and depths at each 0.1° of latitude are randomly chosen daily amongst all possibilities (i.e., throughout the water column and at any longitude between the coastline and the 150 m isobath). Because field observations of the vertical distribution of larvae are lacking, we released particles at random depths but setting a condition that warranted particles being released throughout the water column rather than mostly at shallow depths; i.e. particles releases must occur below and above 50m. Across any given latitudinal section from the coast to the 150 m isobath, there are more possible release locations at shallow than at deeper waters. Thus, we chose a 50 m threshold depth to ensure releases occurred at most depths throughout the water column (Fig 1C–1E). We seeded 13 particles at each 0.1° of latitude from 28–33°S; this implies that 663 particles were released daily for a total of 40,443 particles over a spawning season. Larval transport predictions can be sensitive to the number of particles released (e.g. [45] [46]; a previous study of dispersal in the region corroborates that settlement patterns obtained from dispersal simulations of ERL larvae do not change beyond 455 virtual larvae released per day [25]. The number of releases in the present study is above such threshold. Particles are advected offline with a time-step of 1 hour and with the 3D velocity fields generated by the hydrodynamic model plus a random velocity vector to account for circulation structure at spatial scales smaller than that captured by the hydrodynamic model (i.e., sub-grid scale diffusion). The velocity vector associated with diffusion is scaled by a horizontal diffusion coefficient set to 8.8 m2 s-1 following [47] which has been previously used in dispersal studies with similar hydrodynamic model outputs [25].

Spatial and temporal variability analysis

Comparison of observed and simulated settlement.

Total simulated larval settlement along the coast of southeastern Australia was obtained for each year between 1995–2015 and each 0.1° of latitude; settlement at each 0.1° of latitude was measured as the total amount of larvae that reached any grid cell between the coast and the 150 m isobath at any time 8–12 months after release. Although monthly variability in settlement has been observed, this study focuses on a larger temporal scale, i.e., annual and longer.

The ratio of simulated settlement north and south of 32°S is compared against the ratio of observed settlement in northern (Coffs Harbour and Tuncurry) and southern (Sydney and Ulladulla) locations. North and south simulated settlement is obtained adding settlement occurring at all grid cells north and south of 32°S. Observed north settlement is obtained adding monitored settlement in Coffs Harbour and Tuncurry and observed south settlement is obtained adding monitored settlement in Sydney and Ulladula.

In addition, inter-annual variability of observed settlement summed across all monitored locations is compared against the addition of simulated settlement occurring within all four 0.1° latitude degree closest to the monitored locations (i.e., Coffs Harbour 30.3°S, Tuncurry 32.3°S, Sydney 33.9°S, and Ulladulla 35.5°S). Observed and simulated annual settlement are also compared separately at individual locations in order to discern the spatial scale over which settlement estimates agree with observations (e.g., along southeast Australia or at specific locations). To correlate observed and simulated inter-annual settlement, yearly settlement is first normalized by total settlement across years of observations and simulations respectively. That is, yearly settlement is divided by total settlement over the 21 years and multiplied by 100. Non-oceanographic processes such as a gradual change in size of the breeding stock over time can cause long-term variation in observed settlement. To minimise the effect of such long-term factors on our analysis, we compute yearly changes in settlement and correlate the observed and simulated time series of year-to-year differences in settlement. If inter-annual differences in settlement due to ocean currents are much greater in magnitude than inter-annual differences, due to changes in the abundance of breeding stock and egg production then this approach effectively removes the confounding influence of long-term changes in egg production from the analysis. This approach therefore examines the importance of dispersal by ocean currents in determining year-to-year variations in settlement. Changes in yearly settlement are computed from the normalized yearly settlement time series of observed and simulated settlement respectively.

To examine dominant periodicities or temporal variability over the 21-years, a wavelet analysis is conducted on both the observed and simulated time series of yearly normalized settlement. The wavelet analysis is computed using code from [48], and utilising the bias correction of [49]. Finally, to elucidate if there is a single or more ERL population(s), a mean connectivity matrix was composed across the 21 years; each year, the number of larvae released at each source latitude and arriving to each sink latitude was quantified. Source latitudes span 28–33° S every 0.1° and sink latitudes 26–38° S every 0.1°.

Modelled ocean circulation and simulated settlement.

To interpret settlement in the context of ocean processes, the spatial variability of yearly settlement along the coast obtained from simulations is compared to the variability of the modelled circulation along the coast. To account for circulation variability we focus on two factors; 1) the proximity to the shelf break of a coherent unidirectional flow (i.e., EAC), and 2) the spatial structure of persistent flow identified with empirical orthogonal functions (EOFs).

Here, we explore if regions where the EAC is far from the coast, separates from the coast, or has a dominant component perpendicular to the coastline, coincide with regions of simulated and observed higher settlement. With this purpose we detect daily locations where the EAC departs from an alongshore path as it travels polewards over a time period encompassing the direct observations from the long-term pueruli survey program (1995–2015, Fig 2A).

thumbnail
Fig 2.

a) Map depicting the EAC meandering as it flows polewards and separating from the continent. Snapshot of modelled Sea Surface Height (SSH, colorbar), and the EAC SSH isoline (black contour); the white dots indicate latitudes where the EAC deviates from the 150 m isobath for more than 70 km; the black dot indicates where the EAC reattaches (comes within 70 km from the 150 m isobath after having deviated upstream), the arrows represent the geostrophic currents derived from modelled SSH. b) Across years mean (thick line) and standard deviation (blue shade) of the percentage of settlement at each 0.1°of latitude relative to total settlement across latitudes. The black line shows cumulative settlement from north to south. c) Frequency histogram of the EAC separation latitude over 21 years (1995–2015, bars) and cumulative frequency from north to south (black line).

https://doi.org/10.1371/journal.pone.0211722.g002

Specifically, we record the latitudes where the core of the current is more than 70 km eastward of the 150 m isobath or farther, as well as locations where the stream re-attaches or comes within 70 km of the 150 m isobath (Fig 2A). To identify the EAC core we find the Sea Surface Height (SSH) isoline corresponding to the maximum poleward geostrophic velocity at 25°S between the coast and 160°E. This approach is based on the method used [21] to identify the core of the jet and the latitude where the EAC separates from the continent.

A qualitative comparison between circulation structure and annual settlement is used to discern flow patterns associated with settlement for individual years; particularly, to recognize if there is a common flow pattern for years that show either high or low settlement. The Empirical Orthogonal Functions (EOF’s) capture the dominant modes of variability of the circulation; for instance, the prevalence over time of the EAC, or an eddy field. EOF’s are therefore used to identify the dominant circulation structure of each individual year. To compute the EOF’s we use surface velocity fields over the settling time period of each year; this analysis is designed to reveal the dominant circulation mode during the delivery of larvae to the continental shelf-break rather than that dictating trajectories throughout the pelagic stage, or across the 21-years of observed settlement. For instance, settlement in 1996, which corresponds to larvae hatched in November 1994-January 1995, is qualitatively compared against EOF’s calculated from velocity fields 8–12 months after being hatched (from July 1995 to January 1996). Finally, to further test the role of eddies as a mechanism facilitating settlement, we computed the proportion of larvae transported into the continental shelf-break by an eddy relative to total settlement for a specific year (2003).

Results

Spatial variability

Latitudinal variability in observed and simulated settlement.

Observed settlement over the last 21 years (1995–2015) consistently shows low settlement in northern locations (Coffs Harbour at 30.3°S and Tuncurry at 32.3°S) and high settlement in southern locations (Sydney at 33.9°S, and Ulladulla at 35.5°S); the mean across years of the south to north observed settlement ratio is 11.9°8.4 S.E. Simulated settlement also reveals higher settlement in south compared to north locations every year; nevertheless the south to north ratio, based on simulated settlement north and south of 32°S, (1.85±0.83 S.E.) is almost an order of magnitude smaller than the observed.

Despite evident inter-annual variability, there is a distinct pattern of low simulated settlement (<20% relative to total settlement) between 28–32.3°S (encompassing the northern locations where settlement is surveyed), high settlement (>20% relative to total settlement) between 32.4–36°S, (encompassing the southern locations where settlement is surveyed) and maximum settlement between 32.4–33.8°S (>30% relative to total settlement, Fig 2B). This broad pattern is in agreement with observed settlement patterns. At individual monitoring locations however, 21-year long time series of observed and simulated settlement are not significantly correlated (Fig 3).

thumbnail
Fig 3.

Twenty-one year (1995–2015) time series of observed (black) and simulated (blue) settlement at each monitoring location, correlation coefficients (r) and corresponding p values are shown for correlations between 1995–2015 (top left of each panel), 1995–2005 (mid-left), and 2010–2015 (mid-right). Open bars are used for 2006–2009 to facilitate distinction of decadal periods before and after, for which correlation coefficients and p values are included.

https://doi.org/10.1371/journal.pone.0211722.g003

Latitudinal variability in simulated settlement, EAC meandering and separation.

The latitudinal range of simulated maximum settlement lies within a region where the core of the EAC often deviates from an alongshore path and is not in close proximity of the coast. In fact, a cumulative frequency histogram shows that within the region where the simulated settlement is highest (~32.4–33.8°S), the EAC core is at least 70 km east of the 150 m isobath 40–88% of the time over 21 years (1995–2015, Fig 2C); thus, maximum simulated settlement occurs within typical EAC separation latitudes.

Temporal variability in observed and simulated settlement

Over the 21 years (1995–2015), observed and simulated settlement summed across all four locations (i.e. Tuncurry, Coffs Harbour, Sydney, Ulladula) and corresponding latitudes in the model respectively fluctuates over an order of magnitude; specifically, between 1.9–10.9% and 0.7–12.5% of the total settlement across the 21 years for the observations and simulations respectively (Fig 4A). Thus, the degree of variability along southeast Australia in observed settlement is well matched by the simulations.

thumbnail
Fig 4.

a) Twenty-one year (1995–2015) time series of observed (black) and simulated (blue) settlement added across monitoring locations, correlation coefficients (r) and corresponding p values are shown for correlations between 1995–2015 (r = 0.1, p>>0.05), 1995–2005 (r = 0.78, p < 0.01 filled bars), and 2010–2015 (r = 0.78, p = 0.065 filled bars). Open bars are used for 2006–2009 to facilitate distinction of decadal periods before and after, for which correlation coefficients and p values are included. b) Time series of observed and simulated yearly change in settlement, correlation coefficient (r = 0.63, p < 0.01). c) Observed vs. simulated yearly differences in settlement; filled symbols indicate years when observed and simulated differences agree (i.e., either both increase or decrease), while empty symbols indicate when they show opposite trends (i.e. one decrease and the other increases and vice versa); years with largest discrepancy are emphasized with an asterisk. Settlement and differences in settlement are presented in percentage relative to total settlement across all years for the observed and simulated settlement respectively.

https://doi.org/10.1371/journal.pone.0211722.g004

The correlation between observed and simulated annual settlement (i.e., total settlement each year) over the entire 21 years is not significant at individual monitoring locations or latitudes (Fig 3) or for all for monitoring locations or latitudes summed together, with the simulation only explaining 1% of the variation in annual settlement for data summed across all four monitoring locations or latitudes (Pearson’s test, r = 0.11, p = 0.61, Fig 4A). It is noteworthy that the time series of observed and simulated settlement before and after 2007–2008 oscillate together; in fact, correlations between observed and simulated annual settlement over these two time periods (1995–2005, and 2010–2015) are high, explaining 61% of the variation (both with r = 0.78, p < 0.01 and p = 0.07 respectively, Fig 4A). In contrast, year-to-year changes in simulated settlement (20 data points), a measure of annual settlement partially corrected for long-term changes in breeding stock size, explain 40% of the year-to-year observed variation (i.e., increase or decrease in settlement from one year to the next), and are significantly correlated (r = 0.63, p < 0.01, Fig 4B).

Importantly, years of observed and simulated highest and lowest settlement coincide. For instance, 1999 and 2011 are years characterised by high settlement while 2005 and 2015 are characterized by low settlement, and this is well captured in the simulations (Fig 4). Additionally, a wavelet analysis identifies peaks of settlement occurring every four years in both the observed and simulated 21-year time series of total settlement (Fig 5).

thumbnail
Fig 5. Global wavelet spectrum of simulated and observed settlement showing peaks of energy between 2–3 years and at 4 years.

https://doi.org/10.1371/journal.pone.0211722.g005

Discussion

This study uses an ocean hydrodynamic model coupled with a particle tracking model to simulate the dispersal of ERL larvae. Comparison against 21 years of larval settlement observations allows diagnosis of physical mechanisms driving larval dispersal and supply to the shelf-break, and quantifying the extent to which hydrodynamic processes determine patterns of settlement variability at different temporal and spatial scales. Despite the large range of factors that can impact larval dispersal, we show that along southeast Australia, within a Western Boundary Current (WBC) system, mesoscale circulation explains broad (inter-annual to decadal) spatio-temporal patterns of observed settlement of our model species Sagmariasus verreauxi (Eastern Rock Lobster, ERL) over 500 km of coastline.

Previous studies have proposed a number of mechanisms driving larval dispersal in the region [25, 33]; using observations of lobster settlement over a 21 year period, this study corroborates for the first time the impact of such mechanisms on larval delivery to the coast by comparing simulated with observed settlement. In agreement with the hypothesis proposed by [25], we demonstrate that patterns of settlement abundance are influenced by the presence and coherence of the poleward flowing EAC and the region where it deviates from the coast. More specifically, low and high settlement locations coincide with regions upstream and downstream of the latitude where the EAC separates (Fig 2); similarly, years with low settlement occur when a coherent EAC persists, while high settlement relates to a flow dominated by eddies, or flow perpendicular to the shelf break caused by a discontinuous EAC, EAC meandering, or separation from the continent (Fig 6). In addition, meso-scale eddies which often form where the EAC veers east, but are also present intermittently all along the coast as far north as 25°S, are a mechanism facilitating the transport of larvae onshore into the continental shelf-break.

thumbnail
Fig 6.

a-b) Spatial structure of the first mode of the Empirical Orthogonal Functions (EOFs) showing speed (colorbar) and velocity fields (arrows) for settlement time period in 1999 (a) and 2005 (b). Superimposed is settlement at each 0.1 latitude degree (grey line). c-d) temporal structure of the first mode of the Empirical Orthogonal Functions (EOFs) for the settlement time period in 1999 (c) and 2005 (d); the left hand axis shows settlement over time across all latitudes (grey dots).

https://doi.org/10.1371/journal.pone.0211722.g006

Nevertheless, simulations based solely on the advection of passive larvae by mesoscale circulation do not capture observed settlement at smaller spatial scales (i.e., individual monitoring locations). Intricacies including exact differences in settlement abundance at individual monitoring locations might only be resolved when implementing the impact of local forcing, sub-mesoscale scale circulation, swimming behaviour of pueruli, and the interaction and effects of the environment on larval survival. This is further discussed in the model limitations section below.

Spatial variability in settlement

Effect of the EAC jet, meandering, and separation on settlement.

The EAC facilitates advection of larvae along the coast and connectivity amongst locations separated over large distances, while its meandering and pinching off of eddies induces larval re-circulation and self-seeding. Although connectivity from upstream into downstream locations is dominant, there is also some larval exchange occurring from south to north (Fig 7A and 7B); see also [33, 50, 51]. These results support recent genetic studies suggesting there is a single ERL meta-population along southeast Australia with widespread gene flow [28]. The absence of spatial genetic structure for a variety of species with a planktonic stage in this region [52, 53] indicates that patterns of connectivity are widely influenced by the prevailing circulation in this system irrespective of larval life traits [25, 33].

thumbnail
Fig 7.

a) Trajectory of a particle that settles upstream, the black and grey dots show the release and settlement locations respectively, the particle trajectory is shown for 12 months rather than being truncated at settlement, i.e., at ~9 months after release). b) Mean connectivity matrix (1995–2015), the white indicates values that correspond to self-recruitment (the source and sink latitudes are the same), the colorbar indicates the number of larvae in logarithmic scale. c) Particle re-circulation within eddies and transport into the continental shelf break. Snapshot on the 13/October/2003 showing a simulated settlement pulse and preceding particle re-circulation within eddies with Sea Surface Height (SSH) on the background, dots indicate the position of particles at settlement and lines show their trajectories in the previous 10 days; only particles that reach settlement are shown.

https://doi.org/10.1371/journal.pone.0211722.g007

As a well-formed jet with few instabilities, the EAC has the potential to retain larvae within its core, transport them poleward while inhibiting across-shore transport; and hence, prevent the supply of larvae to the continental shelf-break for subsequent settlement. A modelling study used particle trajectories to show that a jet has a retaining effect and can only leak particles when it experiences instabilities [54]. The EAC can act as a jet, and its potential to inhibit exchange between offshore and coastal waters has been inferred from particle trajectories that rarely move across the EAC [33]. Supporting the hypotheses of 1) a jet inhibiting across-jet transport, and 2) the EAC acting as a jet, we find that the observed low settlement in the north coincides with a latitudinal range (28–31.9°S) where the EAC is more often a continuous poleward stream. Conversely, in the southern section, encompassing the region where the EAC typically separates from the continent, there is high settlement (Fig 2).

From the observed and simulated patterns of low and high settlement up and downstream of the EAC separation respectively; we infer that settlement of ERL in the region where the EAC separates (most frequent at 32 and 33°S, based on 21 years of OFAM BRAN 3p5 data, Fig 2C) benefits from two processes: 1) the advection of larvae by the EAC from upstream spawning locations, and 2) the occurrence of flow meandering that facilitates transport into the continental shelf-break. Similarly, along the coast of Western Australia, the poleward flowing Leeuwin current increases the downstream transport of larvae [55], and settlement in southern regions when the current is strong [56].

Although low, the occurrence of settlement in the north of the region indicates that the EAC is not always coherent within this section of the coast but it also meanders and is interrupted by the presence of eddies. In fact, within this section (28–32°S), the EAC sometimes deviates from an alongshore pathway, having its core located more than 70 km away from the coastline 20–40% of the time over the 21 year period (Fig 2). The EAC meandering and formation of eddies as far north as 25°S has the potential to induce transport from offshore into the continental shelf-break, and flow reversals which would in turn enable settlement in the northern section of the region. Monitoring locations are distributed up and downstream of the EAC separation and are therefore able to capture settlement associated with different circulation regimes.

Effect of mesoscale eddies on settlement.

The impact of eddies retaining larvae and shaping dispersal patterns has been well documented in other regions [5761]. For instance, in New Zealand, high abundances of rock lobster larvae observed in plankton samples coincided with an area of strong eddy circulation [62]. Similarly, the presence of eddies offshore of the Middle Florida Keys induces counter-current anomalies and influx of pueruli [63]. Moreover, [64] suggest that frontal eddies off southeastern Australia have the potential to entrain and retain fish larvae over the pelagic phase allowing recruitment back to the coast. Also, within our study region, modelling studies have shown that eddies retain particles [33], and that areas with more eddy activity relate to higher settlement [25].

This study corroborates the role of eddies in the distribution of larvae and in their role enhancing delivery to the continental shelf-break. Trajectories of larvae arriving to the continental shelf-break in pulses (i.e. in considerable quantities over a few days) consistently show a curved pathway that indicates their transport within an eddy from offshore towards the continental shelf (Fig 7). In addition, settlement peaks along the coast are associated with regions where eddies prevail (Fig 6).

Temporal variability in settlement

Our results show that our model has useful capacity for the prediction of the inter-annual scale of settlement (Fig 4) and of dominant periodicities in settlement variability (Fig 5) over broad spatial scales (i.e., total settlement along southeast Australia). The model captures annual settlement over time periods of up to 10 years as well as year-to-year observed fluctuations in settlement for the 21-year monitoring period (1995–2015). Thus, larval dispersal simulations aid in identifying the physical processes associated with settlement. Years with low settlement (e.g. 2005) tend to have a first EOF that shows a coherent EAC, whereas years with high settlement (e.g. 1999) have first EOF’s that show a velocity field with large spatial variability and eddies (Fig 6). Thus, the dominance of a coherent EAC jet during the period when larvae settle (8–12 months after spawning) consistently translates into years of low settlement; while years when flow is eddy dominated have high settlement (Fig 6). In agreement, we find that the delivery of larvae to the continental shelf-break by eddies alone contributes ~ 30% of simulated settlement in a year.

Although circulation variability patterns relate well to temporal settlement patterns of the ERL along the southeast coast of Australia, there are some years when modelled settlement departs from the observed, suggesting variability in settlement is not determined by ocean circulation alone. We find different inter-decadal trends in observed and simulated settlement, which may be related to the changing size of the breeding stock. An increase in recorded catch per unit effort (CPUE) suggests growth of the breeding population, which could in turn induce the positive trend in observed settlement. The CPUE between 1995–2005 remains relatively constant with highest values of ~2.5 kg per trap-month; in following years CPUE increased continuously, at least until 2010, reaching ~4.5 kg per trap-month in 2007–2008 (Montgomery and Liggins 2013). Interestingly, observed and simulated annual settlements are significantly correlated when partitioned into two periods (1995–2005, and 2010–2015) that disregard the years when CPUE increased monotonically (2006–2010). Hence, the assumption of a constant spawning stock in our simulations could be provoking the dissimilarities we find between observed and simulated long-term settlement trends. Indeed, observed and simulated year-to-year settlement differences, in which long-term fluctuations are minimized, are in better agreement. Thus, incorporating spawning population size on larval dispersal simulations could yield more accurate estimates of annual settlement abundance.

In addition, our simulations consider that mortality is uniform in space and time, and do not include larval responses to the environment (e.g., behavioural or physiological responses). Little is known about changes in the interaction of larval traits with a changing oceanic environment over long time periods (inter-decadal), or the extent to which biological factors affect dispersal. [25] showed that under the A1B future climate change scenario, temperature changes will be favourable for lobster larval survival and consequent settlement particularly at 30°S and further south along southeast Australia. Conversely, for the northern region of the ERL range the seawater temperatures for larvae are increasingly less favourable [8]. Similar climate driven changes (e.g., temperature increases) have occurred in the recent decades along other Western Boundary Currents [65]. Therefore, the impact on settlement of factors that are not considered in our simulations (e.g., effect of temperature on larval survival), may cause the observed (positive) long-term trends in settlement in the south-eastern regions of this study.

Model limitations

This study aims to model the delivery of larvae to the continental shelf-break, and shows that the arrival of pueruli to the coast is influenced by meso-scale circulation; nevertheless, local hydrodynamics, not accounted for in our simulations, may also influence the delivery of larvae and subsequent larval settlement. For instance, [66] found that fluctuations of spiny lobster settlement along the Mexican-Caribbean reflect the influence of mesoscale oceanographic processes, but the impact of waves must be considered over smaller time-scales to estimate settlement within reef lagoons. Similarly, winds that induce onshore transport have been correlated with high settlement [18, 67]. The role of other physical processes on cross shelf recruitment have also been identified in a number of recent studies correlating the settlement and recruitment of spiny lobsters with oceanographic and climatological variables including stokes drift and local wind forcing [10, 66, 68, 69].

In addition, horizontal swimming behaviour of late-stage larvae and pueruli may also play a role in determining arrival on the coast. For instance, [70] find that mid-stage larval distributions of mid-stage larvae of J. edwardsii are well described by passive drift alone, but other mechanisms inducing shoreward transport are needed to explain the shoreward movement of the subsequent late-stage larvae. Furthermore, spiny lobster pueruli are known to be capable of extended periods of cross-shelf swimming of tens to hundreds of kilometres required to reach shallow coastal reefs where they settle [35, 71, 72]. In this study we have not implemented any of these more complex scenarios and yet the predictive capability of settlement inter-annual variation across all sites is significant over decadal periods; however, it prevents the simulations from predicting variation in settlement at individual locations.

Implications for management

Simulated larval trajectories show connectivity throughout southeast Australia (Fig 4) in agreement with genetic studies, which identify a single ERL population in the region. This implies that estimates of settlement over broad spatial scales (total across southeast Australia) are useful to set catch quotas. The relationships we find between the EAC coherence, presence of eddies, and the observed spatio-temporal patterns of settlement suggest that proxies for broad spatio-temporal patterns settlement can be developed from products such as satellite altimetry that capture the EAC variability and the meso-scale eddy field. Such proxies have potential to provide fisheries managers with a useful tool to anticipate broad settlement patterns, and future recruitment into the fishery. This would allow greater capability to set fisheries quotas and sustainably manage this highly valued fishery into the future.

Acknowledgments

The NSW Department of Primary Industries and shareholders in the NSW Rock Lobster Fishery. The research was supported by the Australian Research Council Linkage Project No. LP150100064 to MR, MAC, AJ and GL, and an OECD Co-operative Research Programme Fellowship to AJ. This research is a contribution to the Worldwide Universities Network project; Ocean Eddies in a Changing Climate: Understanding the Impact of Coastal Climates and Fisheries Production.

References

  1. 1. Kough AS, Paris CB, Butler MJIV. Larval Connectivity and the International Management of Fisheries. PLOS ONE. 2013;8(6):e64970. pmid:23762273
  2. 2. Metaxas A, Saunders M. Quantifying the “Bio-” Components in Biophysical Models of Larval Transport in Marine Benthic Invertebrates: Advances and Pitfalls. The Biological Bulletin. 2009;216(3):257–72. pmid:19556593.
  3. 3. Linnane A, James C, Middleton J, Hawthorne P, Hoare M. Impact of wind stress anomalies on the seasonal pattern of southern rock lobster (Jasus edwardsii) settlement in South Australia. Fisheries Oceanography. 2010;19(4):290–300.
  4. 4. Penn JW, Caputi N, de Lestang S. A review of lobster fishery management: the Western Australian fishery forPanulirus cygnus, a case study in the development and implementation of input and output-based management systems. ICES Journal of Marine Science: Journal du Conseil. 2015;72(suppl 1):i22–i34.
  5. 5. Reiss H, Hoarau G, Dickey-Collas M, Wolff WJ. Genetic population structure of marine fish: mismatch between biological and fisheries management units. Fish and Fisheries. 2009;10(4):361–95.
  6. 6. Edgar G J. Predator-prey interactions in seagrass beds. III. Impacts of the western rock lobster Panulirus cygnus George on epifaunal gastropod populations. Journal of Experimental Marine Biology and Ecology. 1990;139(1):33–42. https://doi.org/10.1016/0022-0981(90)90036-C
  7. 7. Provost EJ, Kelaher BP, Dworjanyn SA, Russell BD, Connell SD, Ghedini G, et al. Climate-driven disparities among ecological interactions threaten kelp forest persistence. Global Change Biology. 2017;23(1):353–61. pmid:27392308
  8. 8. Fitzgibbon Q, Ruff N, Tracey S. C. Battaglene S. Thermal tolerance of the nektonic puerulus stage of spiny lobsters and implications of ocean warming2014. 173–86 p.
  9. 9. Liggins GW. Eastern Rock Lobster Sagmarisus verreauxi. Book. Canberra: Fisheries Research and Development Corporation, 2016.
  10. 10. Hinojosa IA, Gardner C, Green BS, Jeffs A, Leon R, Linnane A. Differing environmental drivers of settlement across the range of southern rock lobster (Jasus edwardsii) suggest resilience of the fishery to climate change. Fisheries Oceanography. 2017;26(1):49–64.
  11. 11. Wahle R, Steneck R. Recruitment Habitats and Nursery Grounds of the American Lobster Homarus Americanus: A Demographic Bottleneck?1991. 231–43 p.
  12. 12. Wahle RA, Incze LS. Pre- and post-settlement processes in recruitment of the American lobster. Journal of Experimental Marine Biology and Ecology. 1997;217(2):179–207. https://doi.org/10.1016/S0022-0981(97)00055-5.
  13. 13. Butler MJ, Steneck RS, Herrnkind WF. Juvenile and Adult Ecology. Lobsters: Biology, Management, Aquaculture and Fisheries: Blackwell Publishing Ltd; 2007. p. 263–309.
  14. 14. Caputi N, Chubb CF, Brown RS. Relationships between Spawning Stock, Environment, Recruitment and Fishing Effort for the Western Rock Lobster, Panulirus cygnus, Fishery in Western Australia. Crustaceana. 1995;68(2):213–26.
  15. 15. Gardner C, Frusher SD, Kennedy RB, Cawthorn A. Relationship between settlement of southern rock lobster pueruli, Jasus edwardsii, and recruitment to the fishery in Tasmania, Australia. Marine and Freshwater Research. 2002;52(8):1271–5. https://doi.org/10.1071/MF01032.
  16. 16. Linnane A, McGarvey R, Gardner C, Walker TI, Matthews J, Green B, et al. Large-scale patterns in puerulus settlement and links to fishery recruitment in the southern rock lobster (Jasus edwardsii), across south-eastern Australia. ICES Journal of Marine Science. 2014;71(3):528–36.
  17. 17. Breen PA, Booth JD. Puerulus and juvenile abundance in the rock lobsterJasus edwardsiiat Stewart Island, New Zealand (Note). New Zealand Journal of Marine and Freshwater Research. 1989;23(4):519–23.
  18. 18. Incze LS, Naimie CE. Modelling the transport of lobster (Homarus americanus) larvae and postlarvae in the Gulf of Maine. Fisheries Oceanography. 2000;9(1):99–113.
  19. 19. Ehrhardt NM, Fitchett MD. Dependence of recruitment on parent stock of the spiny lobster, Panulirus argus, in Florida. Fisheries Oceanography. 2010;19(6):434–47.
  20. 20. Ridgway KR, Dunn JR. Mesoscale structure of the mean East Australian Current System and its relationship with topography. Progress in Oceanography. 2003;56(2):189–222. https://doi.org/10.1016/S0079-6611(03)00004-1.
  21. 21. Cetina-Heredia P, Roughan M, van Sebille E, Coleman MA. Long-term trends in the East Australian Current separation latitude and eddy driven transport. Journal of Geophysical Research: Oceans. 2014;119(7):4351–66.
  22. 22. Bull CYS, Kiss AE, Jourdain NC, England MH, van Sebille E. Wind Forced Variability in Eddy Formation, Eddy Shedding, and the Separation of the East Australian Current. Journal of Geophysical Research: Oceans. 2017;122(12):9980–98.
  23. 23. Wilkin John L, Zhang Weifeng G. Modes of mesoscale sea surface height and temperature variability in the East Australian Current. Journal of Geophysical Research: Oceans. 2007;112(C1).
  24. 24. Schaeffer A, Roughan M, Wood JE. Observed bottom boundary layer transport and uplift on the continental shelf adjacent to a western boundary current. Journal of Geophysical Research: Oceans. 2014;119(8):4922–39.
  25. 25. Cetina-Heredia P, Roughan M, van Sebille E, Feng M, Coleman MA. Strengthened currents override the effect of warming on lobster larval dispersal and survival. Glob Chang Biol. 2015;21(12):4377–86. pmid:26268457.
  26. 26. Montgomery SS, Craig JR. Distribution and abundance of recruits of the eastern rock lobster (Jasus verreauxi) along the coast of New South Wales, Australia. New Zealand Journal of Marine and Freshwater Research. 2005;39(3):619–28.
  27. 27. Jeffs AG, Gardner C, Cockcroft A. Jasus and Sagmariasus Species. Lobsters: Biology, Management, Aquaculture and Fisheries: John Wiley & Sons, Ltd; 2013. p. 259–88.
  28. 28. Woodings LN, Murphy NP, Doyle SR, Hall NE, Robinson AJ, Liggins GW, et al. Outlier SNPs detect weak regional structure against a background of genetic homogeneity in the Eastern Rock Lobster, Sagmariasus verreauxi. Marine Biology. 2018;165(12):185.
  29. 29. Montgomery SS, Liggins GW. Recovery of the eastern rock lobster Sagmariasus verreauxi off New South Wales, Australia. Marine Biology Research. 2013;9(1):104–15.
  30. 30. Montgomery S S., Craig J R. Effects of moon phase and soak time on catches of Jasus (Sagmariasus) verreauxi recruits on collectors2003. 847–51 p.
  31. 31. Oke PR, Sakov P, Cahill ML, Dunn JR, Fiedler R, Griffin DA, et al. Towards a dynamically balanced eddy-resolving ocean reanalysis: BRAN3. Ocean Modelling. 2013;67:52–70.
  32. 32. Adcroft A, Hill C, Marshall J. Representation of Topography by Shaved Cells in a Height Coordinate Ocean Model. Monthly Weather Review. 1997;125(9):2293–315.
  33. 33. Roughan M, Macdonald HS, Baird ME, Glasby TM. Modelling coastal connectivity in a Western Boundary Current: Seasonal and inter-annual variability. Deep Sea Research Part II: Topical Studies in Oceanography. 2011;58(5):628–44.
  34. 34. Paris CB, Helgers J, van Sebille E, Srinivasan A. Connectivity Modeling System: A probabilistic modeling tool for the multi-scale tracking of biotic and abiotic variability in the ocean. Environmental Modelling & Software. 2013;42(Supplement C):47–54. https://doi.org/10.1016/j.envsoft.2012.12.006.
  35. 35. Phillips BF, Booth JD, Cobb JS, Jeffs AG, McWilliam P. Larval and Postlarval Ecology. Lobsters: Biology, Management, Aquaculture and Fisheries: Blackwell Publishing Ltd; 2007. p. 231–62.
  36. 36. Rimmer DW, Phillips BF. Diurnal migration and vertical distribution of phyllosoma larvae of the western rock lobster Panulirus cygnus. Marine Biology. 1979;54(2):109–24.
  37. 37. Bradford RW, Bruce BD, Chiswell SM, Booth JD, Jeffs A, Wotherspoon S. Vertical distribution and diurnal migration patterns of Jasus edwardsii phyllosomas off the east coast of the North Island, New Zealand. New Zealand Journal of Marine and Freshwater Research. 2005;39(3):593–604.
  38. 38. Goldsteina Jason S, Butler Mark J. Behavioral enhancement of onshore transport by postlarval Caribbean spiny lobster (Panulirus argus). Limnology and Oceanography. 2009;54(5):1669–78.
  39. 39. Butler MI, Paris CB, Goldstein JS, Matsuda H, Cowen RK. Behavior constrains the dispersal of long-lived spiny lobster larvae. Marine Ecology Progress Series. 2011;422:223–37.
  40. 40. Paris CB, Cowen RK, Lwiza KMM, Wang D-P, Olson DB. Multivariate objective analysis of the coastal circulation of Barbados, West Indies: implication for larval transport. Deep Sea Research Part I: Oceanographic Research Papers. 2002;49(8):1363–86. https://doi.org/10.1016/S0967-0637(02)00033-X.
  41. 41. Sloyan BM, Ridgway KR, Cowley R. The East Australian Current and Property Transport at 27°S from 2012 to 2013. Journal of Physical Oceanography. 2016;46(3):993–1008.
  42. 42. Rykova T, Oke PR, Griffin DA. A comparison of the structure, properties, and water mass composition of quasi-isotropic eddies in western boundary currents in an eddy-resolving ocean model. Ocean Modelling. 2017;114:1–13. https://doi.org/10.1016/j.ocemod.2017.03.013.
  43. 43. Roughan M, Keating SR, Schaeffer A, Cetina Heredia P, Rocha C, Griffin D, et al. A tale of two eddies: The biophysical characteristics of two contrasting cyclonic eddies in the East Australian Current System. Journal of Geophysical Research: Oceans. 2017;122(3):2494–518.
  44. 44. Everett Jason D, Sebille E, Taylor Matthew D, Suthers Iain M, Setio C, Cetina‐Heredia P, et al. Dispersal of Eastern King Prawn larvae in a western boundary current: New insights from particle tracking. Fisheries Oceanography. 2017;26(5):513–25.
  45. 45. Simons RD, Siegel DA, Brown KS. Model sensitivity and robustness in the estimation of larval transport: A study of particle tracking parameters. Journal of Marine Systems. 2013;119–120:19–29. https://doi.org/10.1016/j.jmarsys.2013.03.004.
  46. 46. Monroy P, Rossi V, Ser-Giacomi E, López C, Hernández-García E. Sensitivity and robustness of larval connectivity diagnostics obtained from Lagrangian Flow Networks. ICES Journal of Marine Science. 2017;74(6):1763–79.
  47. 47. Okubo A. Oceanic diffusion diagrams. Deep Sea Research and Oceanographic Abstracts. 1971;18(8):789–802. https://doi.org/10.1016/0011-7471(71)90046-5.
  48. 48. Torrence C, Compo GP. A Practical Guide to Wavelet Analysis. Bulletin of the American Meteorological Society. 1998;79(1):61–78.
  49. 49. Liu Y, Liang XS, Weisberg RH. Rectification of the Bias in the Wavelet Power Spectrum. Journal of Atmospheric and Oceanic Technology. 2007;24(12):2093–102.
  50. 50. Coleman MA, Roughan M, Macdonald HS, Connell SD, Gillanders BM, Kelaher BP, et al. Variation in the strength of continental boundary currents determines continent-wide connectivity in kelp. Journal of Ecology. 2011;99(4):1026–32.
  51. 51. Coleman MA, Cetina-Heredia P, Roughan M, Feng M, van Sebille E, Kelaher BP. Anticipating changes to future connectivity within a network of marine protected areas. Global Change Biology. 2017;23(9):3533–42. pmid:28122402
  52. 52. Banks SC, Piggott MP, Williamson JE, Bové U, Holbrook NJ, Beheregaray LB. OCEANIC VARIABILITY AND COASTAL TOPOGRAPHY SHAPE GENETIC STRUCTURE IN A LONG-DISPERSING SEA URCHIN. Ecology. 2007;88(12):3055–64. pmid:18229840
  53. 53. Coleman MA, Kelaher BP. Connectivity among fragmented populations of a habitat-forming alga, Phyllospora comosa (Phaeophyceae, Fucales) on an urbanised coast. Marine Ecology Progress Series. 2009;381:63–70.
  54. 54. Harrison CS, Siegel DA. The tattered curtain hypothesis revised: Coastal jets limit cross-shelf larval transport. Limnology and Oceanography: Fluids and Environments. 2014;4(1):50–66.
  55. 55. Griffin DA, Wilkin JL, Chubb CF, Pearce AF, Caputi N. Ocean currents and the larval phase of Australian western rock lobster, Panulirus cygnus. Marine and Freshwater Research. 2001;52(8):1187–99. https://doi.org/10.1071/MF01181.
  56. 56. Caputi N. Impact of the Leeuwin Current on the spatial distribution of the puerulus settlement of the western rock lobster (Panulirus cygnus) and implications for the fishery of Western Australia. Fisheries Oceanography. 2008;17(2):147–52.
  57. 57. Lee TN, Clarke ME, Williams E, Szmant AF, Berger T. Evolution of the Tortugas Gyre and its Influence on Recruitment in the Florida Keys. Bulletin of Marine Science. 1994;54(3):621–46.
  58. 58. Nakata H, Kimura S, Okazaki Y, Kasai A. Implications of meso-scale eddies caused by frontal disturbances of the Kuroshio Current for anchovy recruitment. ICES Journal of Marine Science. 2000;57(1):143–52.
  59. 59. Lindo-Atichati D, Bringas F, Goni G, Muhling B, Muller-Karger FE, Habtes S. Varying mesoscale structures influence larval fish distribution in the northern Gulf of Mexico. Marine Ecology Progress Series. 2012;463:245–57.
  60. 60. Simons RD, Nishimoto MM, Washburn L, Brown KS, Siegel DA. Linking kinematic characteristics and high concentrations of small pelagic fish in a coastal mesoscale eddy. Deep Sea Research Part I: Oceanographic Research Papers. 2015;100:34–47.
  61. 61. Funes-Rodríguez R, Ruíz-Chavarría JA, González-Armas R, Durazo R, Guzmán-del Proó SA. Influence of Hydrographic Conditions on the Distribution of Spiny Lobster Larvae off the West Coast of Baja California. Transactions of the American Fisheries Society. 2015;144(6):1192–205.
  62. 62. Booth J D., Stewart R. Distribution of phyllosoma larvae of the red rock lobster Jasus edwardsii off the east coast of New Zealand in relation to the oceanography2017.
  63. 63. Yeung C, Jones D L., Criales MM, Jackson T, Richards W J. Influence of coastal eddies and counter-currents on the influx of spiny lobster, Panulirus argus, postlarvae into Florida Bay2001. 1217–32 p.
  64. 64. Mullaney TJ, Suthers IM. Entrainment and retention of the coastal larval fish assemblage by a short-lived, submesoscale, frontal eddy of the East Australian Current. Limnology and Oceanography. 2013;58(5):1546–56.
  65. 65. Wu L, Cai W, Zhang L, Nakamura H, Timmermann A, Joyce T, et al. Enhanced warming over the global subtropical western boundary currents. Nature Clim Change. 2012;2(3):161–6. http://www.nature.com/nclimate/journal/v2/n3/abs/nclimate1353.html-supplementary-information.
  66. 66. Briones-FourzÁn P, Candela J, Lozano-Álvarez 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. Limnology and Oceanography. 2008;53(3):970–85.
  67. 67. Linnane A, Gardner C, Hobday D, Punt A, McGarvey R, Feenstra J, et al. Evidence of large-scale spatial declines in recruitment patterns of southern rock lobster Jasus edwardsii, across south-eastern Australia. Fisheries Research. 2010;105(3):163–71. https://doi.org/10.1016/j.fishres.2010.04.001.
  68. 68. Feng M, Caputi N, Penn J, Slawinski D, de Lestang S, Weller E, et al. Ocean circulation, Stokes drift, and connectivity of western rock lobster (Panulirus cygnus) population. Canadian Journal of Fisheries and Aquatic Sciences. 2011;68(7):1182–96.
  69. 69. de Lestang S, Caputi N, Feng M, Denham A, Penn J, Slawinski D, et al. What caused seven consecutive years of low puerulus settlement in the western rock lobster fishery of Western Australia? ICES Journal of Marine Science: Journal du Conseil. 2015;72(suppl 1):i49–i58.
  70. 70. Chiswell S, Booth J D. Chiswell S, Booth JD. Rock lobster Jasus edwardsii larval retention by the Wairarapa Eddy off New Zealand. Mar Ecol Prog Ser 183: 227–24401999.
  71. 71. Jeffs AG, Montgomery JC, Tindle CT. How do spiny lobster post‐larvae find the coast? New Zealand Journal of Marine and Freshwater Research. 2005;39(3):605–17.
  72. 72. Wilkin JL, Jeffs AG. Energetics of swimming to shore in the puerulus stage of a spiny lobster: Can a postlarval lobster afford the cost of crossing the continental shelf? Limnology and Oceanography: Fluids and Environments. 2011;1(1):163–75.