Finding lobsters: investigating a period of unusually low settlement of Panulirus cygnus by using larval dispersal modelling
Jessica Kolbusz![https://orcid.org/0000-0003-2779-451X](/media/client/orcid_16x16.png)
A Oceans Graduate School and the UWA Oceans Institute, The University of Western Australia, Crawley, WA 6009, Australia.
B School of Biological Sciences and the UWA Oceans Institute, The University of Western Australia, Crawley, WA 6009, Australia.
C Western Australian Fisheries and Marine Research Laboratories, Department of Primary Industries and Regional Development, Government of Western Australia, Hillarys, WA 6025, Australia.
Marine and Freshwater Research 74(1) 20-38 https://doi.org/10.1071/MF22117
Submitted: 9 June 2022 Accepted: 29 September 2022 Published: 15 November 2022
© 2023 The Author(s) (or their employer(s)). Published by CSIRO Publishing. This is an open access article distributed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND)
Abstract
Context: Spiny lobsters are critically affected by ocean current movements over their long journey from hatching to swimming post-larvae (puerulus).
Methods: Isolating three different larval stages of the Western rock lobster (Panulirus cygnus), we have applied three-dimensional larval-dispersal modelling techniques
Aims: To better understand a period of unusually low puerulus settlement in the late 2000s.
Key results: Model results showed that transport mechanisms assisting pueruli in reaching nearshore habitats may have been inhibited by a strengthening of the southward-flowing Leeuwin Current over the austral summer at the settlement stage. As a result, modelled pueruli were transported significantly less northward, in some cases southward, during and after the period of unusually low settlement. Hatching-stage trajectories also suggest increased dispersal off the coast during 2007 and 2008, transporting modelled larvae further offshore and possibly contributing to decreased recruits into the fishery.
Conclusions: Our study distinguishes the influence that ocean circulation patterns off Western Australia have on three different larval stages.
Implications: Larval dispersal modelling of pueruli, therefore, provides a possible explanation for the changes in puerulus settlement over the past two decades and future considerations.
Keywords: Capes Current, fisheries, larval dispersal, Leeuwin Current, Panulirus cygnus, particle tracking, puerulus, settlement patterns, western rock lobster.
Introduction
Lagrangian analysis of ocean circulation and pathways is an essential tool for understanding the dynamic motions of seawater and the nutrients, heat and biological matter within it (van Sebille et al. 2018; Hawkins et al. 2019; Jones et al. 2020). From a practical and theoretical perspective, it forms the backbone of understanding dispersal throughout the ocean at local or vast spatial and temporal scales. For example, the distribution of marine species during their early larval stages can have varying consequences through high mortality because of current movements away from recruitment locations (Chiswell and Booth 2008) or survival through retention in eddy systems (Hill et al. 1996; Cetina-Heredia et al. 2019) and connecting populations (Chiswell et al. 2003; Cowen and Sponaugle 2009; Nolasco et al. 2018; Pires et al. 2020). Furthermore, understanding these mechanisms for a single species influences their spatial management scales and their likely responses to environmental disturbances or shifts (Kvile et al. 2018; Clavel-Henry et al. 2020).
In the life cycle of spiny lobster species, the pelagic larval stages are of utmost importance to their survival (Fig. 1). The length of this offshore journey varies from several months (Palinurus elephas) to 2 years (Jasus edwardsii; Booth and Phillips 1994). Those with long pelagic larval durations can be transported up to thousands of kilometres offshore. However, behavioural changes at different development stages can control the dispersal, assisting their return to preferential juvenile habitats (Butler et al. 2011). These behaviour patterns include vertical migration regulated by light changes. During the day, larvae (phyllosoma) sink to darker depths and, by night, return to prey-rich surface waters, defined as diurnal vertical migration (DVM). These depths during day and night also can change with age, known as ontogenetic vertical migration (OVM), and vary with the moon phases (Rimmer and Phillips 1979). Phyllosomas are within offshore current systems throughout these migrations (Fig. 1). They build up their lipid reserves to undergo their final metamorphosis into a puerulus (Wilkin and Jeffs 2011; Stanley et al. 2015). At this stage, they have developed swimming abilities. In addition, they are nektonic (non-feeding; Jeffs et al. 2005; Espinosa-Magaña et al. 2018), where their biochemical energy reserves assist their transport across continental shelves and towards reef or seagrass habitats (Phillips et al. 1978).
For the spiny lobster species endemic to the western coast of Australia, the western rock lobster, Panulirus cygnus (P. cygnus), the pelagic larval duration begins at hatching between September and February every year (Fig. 1). Phyllosoma are first immersed in the ocean between 40 and 100 m (Melville-Smith et al. 2009) and, subsequently, driven off the continental shelf by alongshore winds, through Ekman transport, and the region’s dominating eastern boundary current, the Leeuwin Current (LC; Phillips et al. 1981; Feng et al. 2011; Wang et al. 2015). The LC eddies retain the larvae as zooplankton while they feed on soft-bodied prey with a preference for chaetognaths (Saunders et al. 2012; Säwström et al. 2014). They undergo a series of moults with nine distinct life stages, with OVM over 8–10 months in the south-eastern Indian Ocean (Braine et al. 1979). Then, metamorphosis to puerulus occurs before crossing the continental shelf over approximately a month (Fig. 1; Chittleborough and Thomas 1969; Feng et al. 2010, 2011). The puerulus settlement season, whereby these recruits enter the fishery, is from May to April each year (de Lestang et al. 2012). Pueruli numbers are monitored at eight sites within the fishery every new moon (approximately monthly), with the lowest light and highest potential for collection (Phillips et al. 1978). These numbers form a puerulus index (PI) for fisheries management because, at this stage, they have ‘recruited’ into the western rock lobster fishery (de Lestang et al. 2012).
The role of physical oceanographic processes on the survival of P. cygnus phyllosoma during these vulnerable stages has been previously investigated through larval dispersal methods (Griffin et al. 2001; Feng et al. 2011; Caputi et al. 2018). Griffin et al. (2001) used Lagrangian methods in an individual-based model (IBM) in 2001. Modelling two sample years of low, high and average PI, they found a weak LC during early stages dispersed larvae west and north because of the wind. Conversely, as predicted pathways of puerulus crossed the continental shelf, a strong LC caused them to move south. Despite these differences and expected links to environmental factors, no changes were observed in the model PI (Caputi and Brown 1993; Caputi 2008). Therefore, the modelling by Griffin et al. (2001) suggested that biological responses from the hydrodynamic processes were the likely reason for the recorded natural variation in puerulus settlement. Subsequently, Feng et al. (2011) used six categories of larval stage P. cygnus to define these within an IBM. This model utilised Bluelink Reanalysis (BRAN) forecasting (Schiller et al. 2008) to assess varying oceanographic conditions as a form of larval dispersal for the larval-puerulus stage of P. cygnus. Caputi et al. (2018) built on this information by completing 24 model simulations to investigate the decrease in PI since the late 2000s. The model was forced with global and Australia-wide predicted atmospheric data with a resolution of 10-km grids (Caputi et al. 2018).
In 2008 and 2009, contrary to environmental conditions (Brown 2009), an unexpected decrease in the PI was observed, prompting a change in fishery management (de Lestang et al. 2015). After the PI recovered, different patterns in the puerulus settlement were observed. Differences include increased settlement in the north of the fishery (>31°S, Lancelin) and increased settlement in the second half of the season (November–March), as opposed to the previous peak being in October (Kolbusz et al. 2021a). A single cause for the unusually low puerulus settlement and its subsequent changes is unknown, and it is likely that multiple factors are involved. A period of sustained cooler water and low energetic conditions in the south-eastern Indian Ocean overlapped with the low settlement seasons (Kolbusz et al. 2021b). During early austral-summer months, there has been a consistent increase in the strength of the LC, coinciding with a decrease in northward wind-driven Cape Current (CC) c. 2008–09, which may explain the shifts in early settlement after this period (Pattiaratchi and Hetzel 2020; Pattiaratchi and Siji 2020; Kolbusz et al. 2021a). Whether these changes in settlement patterns have occurred primarily as a result of shifting oceanographic conditions or whether a biological mechanism, such as mature females spawning later, is yet to be fully understood.
Understanding why these changes have occurred for future climate-change scenarios and fishery management is essential. By using hindcast (in the past) oceanographic data, we can model the movement of particles and modify their movements to replicate those of a P. cygnus larvae at various developmental stages. As referred to previously, this has been completed using individual-based biophysical larval models (Feng et al. 2011; Caputi et al. 2018). Results showed several factors likely at play, one being onshore current systems mismatching with larvae leaving and returning to the coast (Caputi et al. 2018). Taking into account the finding of low energetic conditions in the south-eastern Indian Ocean (Kolbusz et al. 2021b), we have drawn our attention to the movement of the larvae with some biological considerations. We draw attention to the differences in larval movement before (2000–01), during (2007–09) and after (2010–13) the lower-than-expected puerulus settlement.
Our study aims to explain temporal variability (inter-annual and within a season) in settlement across latitudinal ranges for the early and late settlement periods of P. cygnus. We used three-dimensional (3-D) Lagrangian larval dispersal methods forced with 3-D advection fields from a numerical model. In particular, we simulated the dispersal of phyllosoma on hatching, transport in the open ocean, and pueruli movements towards recruitment areas. We also evaluated differences in transport temporally (early and late) and spatially (north and south) within each settlement season.
Materials and methods
Study region
Off the western coast shelf of Australia, there is the southward-flowing Leeuwin Current (LC, present in the upper 0–300 m) and on the shelf is the northward-flowing wind-driven Capes Current (CC, surface) during the summer months (Fig. 2). The LC system includes the Leeuwin Undercurrent (LUC, >300 m) and the Ningaloo Current (NC) north of 24°S (Fig. 2; Woo and Pattiaratchi 2008; Pattiaratchi and Woo 2009). The LC is a poleward flowing eastern boundary current generated by an oceanic pressure gradient between north-western Australia and the Southern Ocean (Fig. 2c; Hamon 1965; Pearce 1991). This pressure gradient increases during a La Niña period, resulting in a stronger LC, warmer sea-surface temperature (SST) and increased rainfall over the Australian continent (Pattiaratchi and Buchan 1991; Feng et al. 2003; Telcik and Pattiaratchi 2014; Wijeratne et al. 2018; Pattiaratchi and Hetzel 2020; Pattiaratchi and Siji 2020). The opposite is true for El Niño. The mean southward volume transport of the LC peaks around 32.8°S, owing to inflow from the South Indian Counter Current (Wijeratne et al. 2018).
The LC, and associated currents, develop instabilities with the variations in topography over the extent of the current, forming eddies, meanders and offshoots (Andrews 1977; Pearce and Griffiths 1991; Cosoli et al. 2020). The region is associated with the most energetic eddy field along an eastern boundary (Pattiaratchi and Siji 2020). The Abrolhos Islands, the narrowing of the continental shelf slope south around 30°S, and the Perth Canyon off 32°S are major topographic features for the preferential generation of mesoscale eddies between 28 and 33°S (Fig. 2a; Feng et al. 2005; Rennie et al. 2007; Meuleners et al. 2008; Huang and Feng 2015; Cosoli et al. 2020). These eddies can persist for almost 8 months, with a mean radius of ∼100 km and generally retain their original formation (Fang and Morrow 2003; Cosoli et al. 2020). Contrary to other locations, the anticyclonic warm-core eddies of the south-eastern Indian Ocean have higher productivity than do cyclonic cold-core eddies because of entrainment of water from the continental shelf during their formation. As a result, eddies play a fundamental role in controlling species distribution and nutrient fluxes (Greenwood et al. 2007; Moore et al. 2007). The increasing strength of the LC is transferred to increases in meanders and eddies along the continental shelf (Feng et al. 2005).
Surface current variability on the continental shelf within the study region is predominantly wind-driven (Fig. 2a; Smith et al. 1991; Pattiaratchi et al. 1997; Rafiq et al. 2020). Southerly wind stresses overcome the alongshore pressure gradient during the summer months, moving upper layers offshore and favouring upwelling onto the continental shelf (Pearce and Pattiaratchi 1999). The CC is identified through these cooler waters, always inshore of the LC and usually inshore of the 50-m contour. It is initiated at approximately 34°S, flowing northward towards 27°S (Fig. 2a; Gersbach et al. 1999; Woo and Pattiaratchi 2008). Variations in the southerly wind stress and the Australasian monsoon season (January–March) over the austral summer cause offshore migration and weaker LC (Pearce and Pattiaratchi 1999; Woo and Pattiaratchi 2008; Feng et al. 2010; Pattiaratchi and Siji 2020). However, a weaker secondary peak in the LC can occur over December–January (Wijeratne et al. 2018). Thus, the LC migrates offshore over these predominantly sea-breeze-dominated summer months, whereas it floods the shelf during autumn–winter and dominates the distribution of water masses (Cresswell et al. 1989; Pattiaratchi et al. 1997; Woo and Pattiaratchi 2008).
Three-dimensional (3-D) hydrodynamic model
Our Lagrangian larval dispersal model was forced with hourly advection fields from a 3-D hydrodynamic model on the basis of the Regional Ocean Modelling System (ROMS). ROMS (ozROMS, the ROMS configuration for Australia) was run in hindcast mode for Australia and the Indo-Pacific region between 2000 and 2018 (Wijeratne et al. 2018). OzROMS is a fully 3-D, high spatial-resolution (2–3 km) circulation model configured to include the entire Australian continental shelf, slope and the adjacent deep ocean by using ROMS. Further validation of the hindcast dataset is available in the extended validation section of the Supplementary material (Fig. S3–S11). Forcing for ozROMS includes tides, winds, air-sea exchanges and external forcing through the boundaries. These features are a distinct advantage for our application because prior larval dispersal studies of P. cygnus used a coarser resolution along the continental shelf. The details specific to the model set-up, hydrodynamic forcing and validation are found in Wijeratne et al. (2018).
Larval dispersal model
Larval dispersal was completed using the open-source larval transport model Ichthyop (Lett et al. 2008). Ichthyop can track the trajectories of modelled larvae, or particles, in three dimensions through the ozROMS hydrodynamic model forcing at hourly intervals, including inputs for swimming particles, diurnal vertical migration and ontogenetic vertical migration (Lett et al. 2008; Robson et al. 2017). ‘Particles’, or in our case, simulated P. cygnus larvae, are released in the model domain, and the effects of the physical environment, particularly the currents, on the particles can be tracked.
Given the longevity of the P. cygnus larval phase and the uncertainty surrounding their ontogenetic vertical migration (OVM) patterns offshore, we broke down the larval dispersal model into the following three scenarios representative of the three life stages (Fig. 3): (1) hatching (Run H: early stage phyllosoma leaving shelf waters); (2) offshore (Run O: mid-stage phyllosoma); and, (3) settlement (Run S: early pueruli moving shoreward). We deliberately did not replicate the full larval phase because this has been addressed before (Griffin et al. 2001; Feng et al. 2011; Caputi et al. 2018) and is far more subjective to larval behavioural assumptions. Run H simulated newly hatched phyllosoma in their first 90 days, released on the sea bed at 3-km spacing between the 40 and 80-m depth contours, between 26 and 34°S, giving 58 520 particles released each season (Feng et al. 2011). Run O simulated phyllosoma already offshore for 60 days, released at 40-m depth at 120-km spacing in depths greater than 300 m over the whole model domain (22–36°S, out to 108°E), giving 20 070 particles released each season. Run S comprised the following two simulations: (1) puerulus once they have crossed the continental shelf (released on the 80-m contour) as passive particles; and, (2) puerulus at the continental shelf (released on the 200-m contour) as particles with directional swimming capabilities, giving 69 967 particles released each season. In both Run S simulations, particles were released at 40-m depth at a 3-km spacing between 26 and 34°S. The 3-km spacing was chosen on the basis of the size of the ozROMS horizontal resolution (Wijeratne et al. 2018).
![]() |
Existing research has suggested that shifts in the oceanographic patterns over the past 20 years have contributed to the changes in puerulus settlement patterns and perhaps the unusually low settlement in 2008 and 2009 (Kolbusz et al. 2021a). Therefore, representative years for before (2000, 2001), during (2007–09) and after (2010–12) the settlement decrease, which intersected with ozROMS hourly data availability, were simulated for computation efficiency (Robson et al. 2017). The P. cygnus larval-cycle timing denotes that September to February of any given year is the hatching scenario data for the following settlement season (Fig. 3). Therefore, the settlement scenario is for the settlement season in question (Fig. 3). Season throughout refers to the puerulus settlement season (May–April).
Larval behaviour
All particles were forced by hourly ozROMS velocity fields in u and v directions and simulated using the Runge Kutta numerical advection scheme with a time step of 15 s. The small internal time step was chosen so that the Courant–Friedrichs–Lewy (CFL) criterion was not violated, causing instabilities when calculating w velocities. Hourly velocity fields were used, allowing for tides to be resolved, although the tidal range is small (∼0.5 m). The horizontal dispersion rate was set to 1 × 10−9 m2 s−1, as recommended within Ichthyop (Peliz et al. 2007; Lett et al. 2008). Particles were released every 10 days for the hatching (Run H) and offshore scenarios (Run O) and every 5 days for the settlement scenario (Run S; Fig. 3). The larval dispersal model was set to output every 12 and 24 h for Runs O and H respectively.
Within the settlement scenarios, particles were successfully ‘settled’ if they reached depths of less than 40 m, moved for at least 1 day, and stayed in depths of less than 40 m for at least 2 days. By contrast, if they had been moving for over 21 days and not crossed the 40-m contour, they were assumed to be perished because of their nektonic state and inability to survive (Feng et al. 2011).
Particles that moved outside the open ocean boundary (22–36°S, 108–117°E, Fig. 2) were considered out of bounds and removed. Particles were also removed if temperatures were outside survival bounds of 10–26°C. Particles that hit the coast were considered ‘beached’ and removed, unless it was within an appropriate settlement timeframe in Run S. If particles did not leave the continental shelf (>200 m) within 15 days of Run H, they were removed to account for high predation on the shelf (Caputi et al. 2018).
Sensitivity testing and additional releases
Ichthyop did not allow for both ontogenetic (OVM) and diurnal (DVM) vertical migration to co-occur. However, for the first 90 days of P. cygnus life, DVM occurs to a lesser extent, allowing only OVM methods to represent reality and testing against the passive Run H (Griffin et al. 2001; Feng et al. 2011). Ichthyop uses a method similar to that employed by the Connectivity Modelling System (CMS) to model OVM (Paris et al. 2017). The user inputs a probability matrix where the particles should be at time intervals after release. Following an approach similar to that of Feng et al. (2011), the particles stay between 0 and 20 m for the first 30 days. Between 30 and 90 days, the particles move between 0 and 80 m. Additional simulations were completed using OVM for Abrolhos and Big Bank (Fig. 2), key spawning areas (Feng et al. 2011). Particles were released in ozROMS at grid cells between 40- and 30-m depths at the Abrolhos Islands and between 100- and 80-m depth at Big Bank (Table 1).
![]() |
In Run S, swimming scenarios used the swimming action built within Ichthyop. This java code was modified to include a swimming direction instead of a random direction, which is the default. The direction remained at 64° from north towards the coast (Jeffs et al. 2005; Feng et al. 2011; Hinojosa et al. 2016). Three swimming speeds were tested, namely, 0.05, 0.1 and 0.15 m s−1. Given that puerulus not likely to continuously swim at one speed (Phillips and Olsen 1975), 0.05 m s−1 was considered the base Run S for subsequent follow-up results. As an additional sensitivity test, DVM was modelled for the Run S 0.05 m s−1 scenario, with model puerulus near the surface at night and 20 m deep during the day (Feng et al. 2011).
Analysis and interpretation
Run H: hatching scenario
The differences in the dispersal of particles after their first 3 months of life was the key response variable for Run H. Therefore, after 90 days, the locations of the particles for each iteration of Run H were collated into 1° bins. This was completed for early (September–December) and late (January–March) releases to correspond to the early and late spawning events, resulting in early and late settlements in the following year. The number of simulated phyllosoma within each degree bin was transformed into a percentage of the total present after 90 days to enable comparison among years.
Run O: offshore scenario
For Run O, particle dispersal was mapped at 10, 30 and 60 days after release. These were transformed into a percentage of particles still present within the number of days to make all years comparable.
Run S: settlement scenario
Each Run S output was processed by extracting only successfully recruited simulated puerulus that settled between 27.5 and 32.5°S to align with the locations of observed puerulus settlement (de Lestang et al. 2012). If particles settled after October, they were treated as settling in the latter portion of the season because of the changes in settlement patterns characterised by Kolbusz et al. (2021a). The latitudinal distribution of settled particles was explored through histograms for passive particles and those with swimming capabilities. These histograms were completed for the early and late portions of the season. Grouping the release latitude and settlement latitudes of model puerulus into connectivity matrices gave a spatially comparable indication of how the particles moved before successful settlement.
Given the results of the connectivity matrices, the shift between the release and settlement latitudes was explored further. Further exploration was completed only for Run S 0.05 m s−1 scenario because the number of successful particles settling was higher than passive, but less different when the speed was increased. The following three fixed factors were used to test changes in the mean latitude shift from release to settlement location before, during and after the low puerulus settlement: in the early or late season (two levels; early and late), the location of settlement (two levels; north, >31°S; and south, <31°S) and before, after or during the period of low settlement (three levels; before, during and after). For example, if particles settled north of 31°S, they were said to have settled north of the fishery (Kolbusz et al. 2021a). The mean latitude shift for each successful model puerulus recruit was averaged for both the early and late settlement season and separated into the northern and southern settlement. The years were treated as random and used as replicates nested within period (eight levels; 2000–01, 2007–12). An ANOVA was conducted to test for changes in the mean latitude shift of particles using the PERMANOVA+ add-on package for PRIMER (ver. 6, see https://www.primer-e.com/, using 9999 permutations; Anderson et al. 2008) and a Euclidean distance measure. Significant terms and interactions were investigated using a posteriori pairwise comparison (see Anderson 2001).
A comparison of DVM particles was completed with the Run S 0.05 m s−1 as the baseline. In addition, the contrast between the DVM and non-DVM particles was explored by dividing them into early or late successful settlements, as was completed for all Run S simulations.
Results
Only ‘alive’ particles, defined as those currently active, were considered in the analysis for all scenarios. There was an average of 88% mortality within Run H in passive-transport simulations and 85% in the OVM transport simulations (Supplementary Table S1). There was an average of 44% mortality within Run O (Supplementary Table S2) and 70% mortality or not-recruited within Run S (Supplementary Table S3).
Hatching scenario (Run H)
Run H found that during and after the period of unusually low settlement, particles were dispersing to a greater extent in both early and late releases (Fig. 4). However, in 2001 and 2007, the higher-density area was further north (Fig. 4). After 2007, the late release particles shifted southward. Before the low puerulus settlement, there were specific areas of higher-density particle retention, which were likely to be within eddy features (Fig. 4). Both passive and ontogenetic verticle migration (OVM) particles had high mortality for beaching in 2007 and 2008 (38%; Table S1).
Ontogenetic vertical migration releases
After 90 days, the locations of OVM particles were reasonably similar to those transported passively, with both displaying less concentration of particles after 2007 (Fig. 5g–l). After this time, dispersal to the north and south is more extensive for early and late releases throughout the unusually low settlement. There was a considerable dispersal of particles to the north in 2007, dissimilar to the corresponding passively transported simulation. Interestingly, in 2012 after the period of unusually low settlement, the higher-density locations were close to the shelf. There is a higher-density offshore at 33−34°S, but it is still more dissipated than earlier.
Key spawning location releases
Run H particles released from the Abrolhos and Big Bank with OVM show similar dispersal patterns in the early releases (Fig. 6), aside from 2012, where wider dispersal occurred. In the early releases of 2007, a large portion (∼10%) of particles ended up at 32°S offshore, which is inconsistent with the other years (Fig. 6e). Late releases (Fig. 6) showed higher spatial variation than did early releases, with the highest-density area moving further north and inshore with each simulation. Still, by 2012, the particles were widely dissipated.
Offshore scenario (Run O)
Offshore releases of particles, representing the mid-stage phyllosoma life-stage over the winter months, showed only minor differences among years (Fig. S1). Despite the same number of particles, a large portion (∼20%) left the boundary in 2000 after only 10 days. Approximately 5% of particles died as a result of being outside of the temperature bounds (Table S2). In 2011, a higher particle density was along the shelf region at 34°S after 30 days. Although by 60 days, the density of particles increased in areas off the shelf up to 112°E (Supplementary Fig. S1). The same patterns were seen in 2000. There was a reduction in relative density in the northern-most sections during 2011 and 2012 compared with the earlier years (Fig. S1).
Settlement scenario (Run S)
Passive particles
The maximum number of particles to settle in Run S (passive) was in 2009 (4362, Fig. 7a), with the majority of particles settling in the earlier portion of the season. However, there was an increase in settlement in the latter portion of the season after 2001 (Fig. 7a). An average of 74% of particles did not reach the nearshore, with 20% outside of temperature bounds (Table S3).
Before and after the low puerulus settlement, the peak latitude of settlement distribution peaked at 28.5°S. During the low settlement, this peak was further south at 30°S (Fig. 8a). The latitudinal distribution of the late settlement was similar to that of the early settlement (Supplementary Fig. S2). In both early and late portions of the seasons, there was also a drop in the settlement between 29.5 and 29°S each year (Fig. 9a, Fig. S2).
Swimming particles
The number of simulated puerulus settling increased when swimming capabilities were added (Fig. 7b–d), with the conditions in the early portion of the season (September–November) being more favourable than those later (December–March) in the season (Fig. 8c–d). However, an increase in swimming speed did not considerably increase the survival of particles (Fig. 7, Table S3). Furthermore, settlement decreased over the early and late portions of the season in 2007 and 2008 and was most prominent for the 0.05 m s−1 simulations (Fig. 7b).
The latitudinal distribution of modelled settlement occurred in two prominent peaks, similar to the non-swimming scenario (Fig. 8a). However, the inclusion of swimming caused the distributions to widen into a more uniform spread across the latitudes. This change was apparent as the swimming speed increased. In most scenarios and seasons (including passive), very little settlement occurred north of 28°S, and there was also less settlement just south of 29°S. The latter portion of the season had similar distribution patterns for all swimming particles (Fig. S2).
Source–sink relationships
For the passive model puerulus, during the early portion of the season, the most important release region was ∼27 to 28°S (Fig. 9a). However, this was further south for the latter portion of the season, ∼30 to 31°S for 2000–01 and 2007–09. It shifted to be more similar to the early portion with a northern release latitude (Fig. 9a). The spread of the release-settlement groupings moves further below the black line (where release and settlement latitudes are the same) over time in both the early and later portions of the season (Fig. 9a). As a result, this shows that the particles are being transported further South before eventually settling.
When swimming speed was added (0.05 m s−1), as particles moved directly toward the coast, their release and settlement latitudes were more similar (black line, Fig. 9b). The most important release latitude was 30–32°S for early settlement. This was less distinguishable for the latter settlement, where the latter settlement grouping was more dispersed than the earlier (Fig. 9b). However, both early and late portions of the season showed that the shift in latitude between release and settlement moved southward during and after the period of observed reduced settlement (Fig. 9b).
A comparison of the latitudinal shift between the release and settlement latitude was calculated for all 0.05 m s−1 swimming particles (Supplementary Table S4). Before the low settlement, particles from the north and south were typically transported northward (Fig. 10a). However, during and after the unusually low puerulus settlement, particles were typically transported southward if they originated in the north (during t = 5.87, Pmc < 0.01, after t = 3.999, Pmc < 0.05), and more weakly northward if they originated in the south (during t = 8.153, Pmc < 0.005, after t = 7.024, Pmc < 0.01; Fig. 10a). Over the whole simulation period, transport in the north was not significantly different over the early and late portions of the season. However, in the south, the source–sink pattern of the late portion of the season changed to a significantly northward transport (t = 4.512, P < 0.005; Fig. 10b).
Overall, particles were transported more southwards in the north during and after the period of low settlement (Fig. 10). Trajectory mapping further supports this pattern, where settled particles in a predominant location in the north are more likely to originate from further northwards during and after the period of low settlement (Fig. 11). In comparison to the south, where before the period of low settlement, particles typically originate from a range of latitudes along the continental shelf in both early and late releases (Fig. 12); however, over time there has been a tendency for particles to originate from further north, particularly in the late releases (Fig. 12).
Diurnal vertical migration
The pattern of settlement, allowing for diurnal vertical migration (DVM) behaviour of particles (Fig. 13), was relatively comparable to simulations without DVM (Fig. 7b). However, on average, there was 15% less settlement in later periods than earlier in the settlement season (Fig. 13). In addition, the late settlement showed a clear declining pattern until 2007 and then incrementally recovered to past the previous maximum by 2012 (Fig. 13). Again, this varies from the simulations without DVM (Fig. 7b), where the late settlement exhibited as much temporal variation. An additional difference was an abrupt decline in the early settlement in the 2010 season (Fig. 13), which was not as pronounced without DVM (Fig. 7b).
![]() |
Discussion
Our results suggest that the variations in observed puerulus settlement between 2000 and 2012 are a product of variation in the localised driving forces of the Leeuwin (LC) and Capes (CC) currents over the summer months. The years of low settlement coincided with a low energy period in the south-eastern Indian Ocean, where there were neutral ENSO conditions over an extended period (Pattiaratchi and Siji 2020; Kolbusz et al. 2021b). In 2011, when there was a strong La Niña, recovery in puerulus numbers began, but not to the expected extent based on the strength of the LC (de Lestang et al. 2015). This is surmised in Fig. 14. We have demonstrated that the cause of settlement failure could be identified by isolating the different early life-cycle phases of P. cygnus and the forcing mechanisms at each dispersal phase.
Our settlement results (Run S) have shown that particles are being transported from locations substantially and significantly further north during the period of unusually low settlement compared with prior. This trend also continued after the period of unusually low settlement (2010–12). The latitudinal shift for north-settling particles was similar in the early and late portions of the season. However, particles moved significantly further north to settle during the latter part of the season in the south. This matches observed changes in puerulus settlement, with a noticeable drop in the PI in the southern sites during and after the recruitment failure (Kolbusz et al. 2021a). In 2010 (the late portion being in the early 2011 summer) and 2011, the LC was strong over the later portions of those seasons, thereby weakening the CC by flooding the shelf (Benthuysen et al. 2014); this is likely to have decreased the number of observed puerulus compared with what was typically expected during a strong La Niña (Caputi 2008). The difference in the early and late latitudinal shifts also matched the differences in recovery of the puerulus index after the period of unusually low settlement experienced between the north and south of the fishery (Kolbusz et al. 2021a).
Whether north or south, an increase in the degree difference between the release and settlement latitude may imply that the distance travelled by the model pueruli has increased. The distance travelled for a puerulus is highly important for its survival because of its nektonic state (Lemmens 1994; Fitzgibbon et al. 2014). The further they travel, the higher the utilisation of energy reserves. With the increased dispersal of particles shown in our hatching runs, phyllosomas reaching the continental shelf to undergo metamorphosis may already have depleted energy reserves if they were not within eddy systems (Saunders et al. 2012; Waite et al. 2019). Model pueruli that successfully settled in the south of the fishery, particularly during the late portion of the season, travelled further to reach settlement locations. In the north, despite there being a southward shift in the direction they travelled to get to the settlement location (because of the LC summer increase), the distance travelled was more uniform over the periods. The differences in distance travelled may explain the more significant recovery in PI observed in the northern parts of the fishery (Kolbusz et al. 2021a).
The division of the larval cycle into three simulations (Run H, Run O and Run S), drew attention to the dispersal dynamics of P. cygnus larvae, bypassing biological considerations such as temperature-dependent growth (Liddy et al. 2004; Steele and Neuheimer 2022), primary productivity (Caputi et al. 2001) and food availability (Waite et al. 2019). On the eastern coast of Australia, decreases in the pelagic larval duration for Melicertus plebejus (prawn; Everett et al. 2017) and Pomatomus saltatrix (bluefish; Schilling et al. 2020) have been due to degree-day growth changes, in combination with mesoscale processes. Complete biophysical models considering degree-day growth have also been completed for P. cgynus (Feng et al. 2011; Caputi et al. 2018); however, the culmination of this research and the above considerations highlight the complexities surrounding the years of low settlement of P. cgynus. The difference in the movement of particles during and after the period of low settlement, combined with an increase in temperature from the LC, would have reduced the intermoult period and possibly their size (Liddy et al. 2004). Therefore, a decrease in their pelagic larval duration during early settlement could have led to a decreased CC (pre-November), possibly discouraging the larvae off the shelf.
After 3 months, the wider dispersal of particles by 2012 highlights that there has been less entrainment of particles within eddy structures, likely reducing their return pathway to the coast. The entrainment of phyllosoma and other long-lived larvae within eddies is regarded as positive for prey availability and transport (Mullaney and Suthers 2013; Cetina-Heredia et al. 2019; Woodings et al. 2019; Clavel-Henry et al. 2020). It has previously been suggested that the period of unusually low settlement may have occurred because of a ‘rogue’ eddy system entraining newly hatched larvae from the key spawning areas and transporting them offshore (Brown 2009; Säwström et al. 2014; Waite et al. 2019). However, Run H has shown that this is unlikely to have occurred. Instead, there is a tendency for particles to become increasingly dispersed leading up to, during and after the period of low settlement, with limited evidence of entrainment in a single eddy structure.
Particles released offshore during austral winter (Run O) were also influenced by the strength of the LC and an associated increase in eddy kinetic energy (EKE). The number of particles lost from the simulations was highest in 2000 when EKE in the system peaked (Pattiaratchi and Siji 2020; Kolbusz et al. 2021b). There was a similarly increasingly southern distribution of particles in 2011 to 2000, which coincided with another peak in EKE (Pattiaratchi and Siji 2020; Kolbusz et al. 2021b). The only simulation where we observed an increased density of particles overlapping with the flow of the LC was after 30 days in 2011 when the LC was undoubtedly the strongest (Boening et al. 2012; Benthuysen et al. 2014; Wijeratne et al. 2018). Alongside the increased dispersal in Run H, the strong LC in 2011 may have been too strong for phyllosoma to survive. Fisheries managers expected a significant increase in settlement given the strong La Niña conditions in 2011 (de Lestang et al. 2015). However, only a small increase in PI was observed (Kolbusz et al. 2021a). The LC also had a more considerable spatial extent (Huang and Feng 2015) throughout the unusually low settlement, which may have resulted in inadequate energy in the system to increase the PI substantially.
Despite minor inter-annual differences, the number of particles that successfully settled during Run S remained relatively consistent. This was expected because we did not model transport over the complete life-cycle of the larvae (but see Feng et al. 2011; Caputi et al. 2018). A slight decrease was recorded over 2007 and 2008, in line with the period of low settlement when swimming was added (Fig. 7b); however, not in 2009, where a low PI was also recorded (de Lestang et al. 2015; Kolbusz et al. 2021a). The lower number of successfully settled particles in the late portion of the season in Run S, including the diurnal vertical migration (DVM) sensitivity runs, suggests that the weak secondary LC peak over January and February inhibits puerulus from reaching the nearshore. The same is not likely to occur over the early portion of the season where the CC is flowing northward, and LC would usually be weaker (Gersbach et al. 1999; Wijeratne et al. 2018).
Separating the successfully settled particles from Run S into their locations along the coast showed areas where shelf and nearshore processes assist puerulus transport to likely settlement locations. More temporal variability in settlement latitudes was seen within the early settlement than the late, although, overall, they showed similar regions of increased settlement. Both passive- and swimming-particle simulations indicated a region in the north (28.2–28.8°S), where settlement decreases throughout the low settlement period, causing the distribution to skew southward. This north region includes Dongara, where the PI has historically been the highest and decreased significantly during the period of unusually low settlement (de Lestang et al. 2015). This suggests that the transport conditions during the period of the unusually low settlement were not conducive to transporting puerulus towards these key latitudes within the fishery.
Consistent with the literature on other spiny lobster species, our research showed the influences of nearshore oceanographic processes on the settlement of pueruli. For Jasus edwardsii, south-westerly wind-driven surface currents also assist transport to the South Australian coastline (Hinojosa et al. 2017). In addition, hurricanes and their associated sea-level changes have been correlated with Caribbean spiny lobster (Panurlius argus) settlement patterns, suggesting that the conditions assist the movement of recruits onshore (Briones-Fourzán et al. 2008). Throughout the unusually low settlement for P. cygnus, the boundary current, the LC, was stronger over the summer, flooding the shelf, and reducing the CC strength. At the same time, it functions to assist puerulus settlement. The same occurs within the East Australia Current, where the jet-like stream can act as a barrier to the settlement of Sagmariasus verreauxi (Roughan et al. 2011; Cetina-Heredia et al. 2015). In essence, it is a combination of offshore abundance and condition, followed by active swimming and oceanographic shelf processes that drive the successful settlement of P. cygnus and other spiny lobster species globally (Yeung et al. 2001; Goldstein and Butler 2009; Linnane et al. 2010; Singh et al. 2018).
An increase in the dispersal properties of the LC during hatching, alongside its increased strength in summer, while the CC decreased, is likely to have given rise to alternative pathways for puerulus reaching the nearshore after the low settlement period. It is important to note that with only minimal field observations, little is definitively known about the movements of P. cygnus during its larval duration, in particular, cues for metamorphosis into puerulus and likely swimming speeds (Jeffs et al. 2005; Phillips and McWilliam 2009). This is a limitation of our study, alongside the assumptions of larval dispersal modelling. In particular, biological processes such as variation in productivity and predation were not factors in this analysis aside from a timing constraint (40 days) to leave the shelf in Run H. However, we aimed to highlight the possible influence that oceanographic conditions could have on planktonic and weakly swimming particles before, during and after this period of unusually low settlement.
Our results have added to the growing body of knowledge surrounding the influence of the LC and coastal currents on P. cygnus larvae and settlement patterns represented by PI (Caputi and Brown 1993; Pearce and Phillips 1988; Kolbusz et al. 2021b). Our research indicated that localised influences of the LC interacting with CC can positively and negatively affect the successful settlement of P. cygnus puerulus. However, this depends on where and when it has reached the shelf and whether it is in the early or late portion of the season. This research, therefore, calls for greater monitoring of the CC strength when considering larval stages of species over the summer months along the western coast of Australia, for example, characterisation of the spatial scale of the LC and CC from satellite images in austral winter and summer. Importantly, larval dispersal modelling should consider it necessary for the CC to be resolved.
Conclusions
Our findings have significant implications on fishery management and understanding of the transport mechanism of the early life of the western rock lobster. The southward-flowing Leeuwin Current (LC) and the northward-flowing Capes Current (CC) on the shelf interact over summer to generate different pathways for puerulus to reach nearshore areas. Before the unusually low puerulus settlement period, the LC was weaker, therefore not impeding the CC northward shelf flow. This allowed pueruli to be transported more directly to suitable habitats. During and after the period of low settlement, puerulus reached nearshore areas more successfully from further north, particularly in the southern latitudes of the fishery. The same occurred afterwards, although to a lesser extent. These differences were highlighted statistically and through trajectory mapping. We have shown that it is not so much the strength of the LC but somewhat nuanced pathways influencing the CC in summer, that characterises puerulus settlement into the fishery.
Supplementary material
Supplementary material is available online.
Data availability
Monthly, daily and hourly mean 3-D currents, temperature and salinity together with mean water levels available for the ozROMS simulations described in this paper are available through contacting the first author.
Conflicts of interest
The authors declare that they have no conflicts of interest.
Declaration of funding
J. Kolbusz was supported by the UWA University Postgraduate Award and Australia Research Training Program scholarships.
Acknowledgements
We acknowledge the developers of Ichthyop and their feedback. Access to the supercomputing facilities of the Pawsey Centre was enabled through the partner allocation scheme. The authors thank Kye Adams for his assistance with statistics and Danielle Su for her guidance in the initial stages of learning Ichthyop.
References
Anderson, MJ (2001). A new method for non-parametric multivariate analysis of variance. Austral Ecology 26, 32–46.Anderson MJ, Gorley RN, Clarke KR (2008) ‘PERMANOVA for PRIMER: Guide to Software and Statistical Methods.’ (PRIMER-E: Plymouth, UK)
Andrews, JC (1977). Eddy structure and the West Australian current. Deep-Sea Research 24, 1133–1148.
| Eddy structure and the West Australian current.Crossref | GoogleScholarGoogle Scholar |
Benthuysen, J, Feng, M, and Zhong, L (2014). Spatial patterns of warming off Western Australia during the 2011 Ningaloo Niño: quantifying impacts of remote and local forcing. Continental Shelf Research 91, 232–246.
| Spatial patterns of warming off Western Australia during the 2011 Ningaloo Niño: quantifying impacts of remote and local forcing.Crossref | GoogleScholarGoogle Scholar |
Boening, C, Willis, JK, Landerer, FW, Nerem, RS, and Fasullo, J (2012). The 2011 La Niña: so strong, the oceans fell. Geophysical Research Letters 39, 1–5.
| The 2011 La Niña: so strong, the oceans fell.Crossref | GoogleScholarGoogle Scholar |
Booth, J, and Phillips, B (1994). Early life history of spiny lobster. Crustaceana 66, 295–317.
| Early life history of spiny lobster.Crossref | GoogleScholarGoogle Scholar |
Braine SJ, Rimmer DW, Phillips BF (1979) An illustrated key and notes on the phyllosoma stages of the western rock lobster Panulirus cygnus George. Report 102. (CSIRO Division of Fisheries and Oceanography: Sydney, NSW, Australia) Available at http://www.cmar.csiro.au/e-print/open/CMReport_102.pdf
Briones-Fourzán, P, Candela, J, and Lozano-Álvarez, E (2008). 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 53, 970–985.
| Postlarval settlement of the spiny lobster Panulirus argus along the Caribbean coast of Mexico: patterns, influence of physical factors, and possible sources of origin.Crossref | GoogleScholarGoogle Scholar |
Brown R (2009) ‘Western Rock Lobster Low Puerulus Settlement Risk Assessment Workshop’, 1–2 April 2009, Perth, WA, Australia. Fisheries Occasional Publication number 71. (Department of Fisheries: Perth, WA, Australia)
Butler, M, Paris, C, Goldstein, J, Matsuda, H, and Cowen, R (2011). Behavior constrains the dispersal of long-lived spiny lobster larvae. Marine Ecology Progress Series 422, 223–237.
| Behavior constrains the dispersal of long-lived spiny lobster larvae.Crossref | GoogleScholarGoogle Scholar |
Caputi, N (2008). 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 17, 147–152.
| 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.Crossref | GoogleScholarGoogle Scholar |
Caputi, N, and Brown, R (1993). The effect of environment on puerulus settlement of the western rock lobster (Panulirus cygnus) in Western Australia. Fisheries Oceanography 2, 1–10.
| The effect of environment on puerulus settlement of the western rock lobster (Panulirus cygnus) in Western Australia.Crossref | GoogleScholarGoogle Scholar |
Caputi, N, Chubb, C, and Pearce, A (2001). Environmental effects on recruitment of the western rock lobster, Panulirus cygnus. Marine and Freshwater Research 52, 1167–1174.
| Environmental effects on recruitment of the western rock lobster, Panulirus cygnus.Crossref | GoogleScholarGoogle Scholar |
Caputi, N, Feng, M, Denham, A, de Lestang, S, Penn, J, Slawinski, D, Pearce, A, and How, J (2018). Optimising an oceanographic-larval model for assessment of the puerulus settlement of the western rock lobster, Panulirus cygnus, in Western Australia. Bulletin of Marine Science 94, 775–800.
| Optimising an oceanographic-larval model for assessment of the puerulus settlement of the western rock lobster, Panulirus cygnus, in Western Australia.Crossref | GoogleScholarGoogle Scholar |
Cetina-Heredia, P, Roughan, M, van Sebille, E, Feng, M, and Coleman, MA (2015). Strengthened currents override the effect of warming on lobster larval dispersal and survival. Global Change Biology 21, 4377–4386.
| Strengthened currents override the effect of warming on lobster larval dispersal and survival.Crossref | GoogleScholarGoogle Scholar |
Cetina-Heredia, P, Roughan, M, van Sebille, E, Keating, S, and Brassington, GB (2019). Retention and leakage of water by mesoscale eddies in the East Australian Current system. Journal of Geophysical Research: Oceans 124, 2485–2500.
| Retention and leakage of water by mesoscale eddies in the East Australian Current system.Crossref | GoogleScholarGoogle Scholar |
Chiswell, S, and Booth, J (2008). Sources and sinks of larval settlement in Jasus edwardsii around New Zealand: where do larvae come from and where do they go? Marine Ecology Progress Series 354, 201–217.
| Sources and sinks of larval settlement in Jasus edwardsii around New Zealand: where do larvae come from and where do they go?Crossref | GoogleScholarGoogle Scholar |
Chiswell, SM, Wilkin, J, Booth, JD, and Stanton, B (2003). Trans-Tasman Sea larval transport: Is Australia a source for New Zealand rock lobsters? Marine Ecology Progress Series 247, 173–182.
| Trans-Tasman Sea larval transport: Is Australia a source for New Zealand rock lobsters?Crossref | GoogleScholarGoogle Scholar |
Chittleborough, RG, and Thomas, LR (1969). Larval ecology of the Western Australian Marine Crayfish, with notes upon other panulirid larvae from the Eastern Indian Ocean. Marine and Freshwater Research 20, 199–223.
| Larval ecology of the Western Australian Marine Crayfish, with notes upon other panulirid larvae from the Eastern Indian Ocean.Crossref | GoogleScholarGoogle Scholar |
Clavel-Henry, M, Solé, J, Kristiansen, T, Bahamon, N, Rotllant, G, and Company, JB (2020). Modeled buoyancy of eggs and larvae of the deep-sea shrimp Aristeus antennatus (Crustacea: Decapoda) in the northwestern Mediterranean Sea. PLoS ONE 15, e0223396.
| Modeled buoyancy of eggs and larvae of the deep-sea shrimp Aristeus antennatus (Crustacea: Decapoda) in the northwestern Mediterranean Sea.Crossref | GoogleScholarGoogle Scholar |
Cosoli, S, Pattiaratchi, C, and Hetzel, Y (2020). High-frequency radar observations of surface circulation features along the south-western australian coast. Journal of Marine Science and Engineering 8, 97.
| High-frequency radar observations of surface circulation features along the south-western australian coast.Crossref | GoogleScholarGoogle Scholar |
Cowen, RK, and Sponaugle, S (2009). Larval dispersal and marine population connectivity. Annual Review of Marine Science 1, 443–466.
| Larval dispersal and marine population connectivity.Crossref | GoogleScholarGoogle Scholar |
Cresswell, GR, Boland, FM, Peterson, JL, and Wells, GS (1989). Continental shelf currents near the Abrolhos Islands, Western Australia. Marine and Freshwater Research 40, 113–128.
| Continental shelf currents near the Abrolhos Islands, Western Australia.Crossref | GoogleScholarGoogle Scholar |
de Lestang S, Caputi N, How J, Melville-Smith R, Thomson A, Stephenson P (2012) Stock assessment for the west coast rock lobster fishery. Fisheries Research Report 217, Department of Fisheries, Western Australia.
de Lestang, S, Caputi, N, Feng, M, Denham, A, Penn, J, Slawinski, D, Pearce, A, and How, J (2015). What caused seven consecutive years of low puerulus settlement in the western rock lobster fishery of Western Australia? ICES Journal of Marine Science 72, 49–58.
| What caused seven consecutive years of low puerulus settlement in the western rock lobster fishery of Western Australia?Crossref | GoogleScholarGoogle Scholar |
Espinosa-Magaña, AF, Briones-Fourzán, P, Jeffs, A, and Lozano-Álvarez, E (2018). Energy cost of the onshore transport of postlarvae of the Caribbean spiny lobster, Panulirus argus. Bulletin of Marine Science 94, 801–819.
| Energy cost of the onshore transport of postlarvae of the Caribbean spiny lobster, Panulirus argus.Crossref | GoogleScholarGoogle Scholar |
Everett, JD, van Sebille, E, Taylor, MD, Suthers, IM, Setio, C, Cetina-Heredia, P, and Smith, JA (2017). Dispersal of Eastern King Prawn larvae in a western boundary current: new insights from particle tracking. Fisheries Oceanography 26, 513–525.
| Dispersal of Eastern King Prawn larvae in a western boundary current: new insights from particle tracking.Crossref | GoogleScholarGoogle Scholar |
Fang, F, and Morrow, R (2003). Evolution, movement and decay of warm-core Leeuwin Current eddies. Deep-Sea Research – II. Topical Studies in Oceanography 50, 2245–2261.
| Evolution, movement and decay of warm-core Leeuwin Current eddies.Crossref | GoogleScholarGoogle Scholar |
Feng, M, Meyers, G, Pearce, A, and Wijffels, S (2003). Annual and interannual variations of the Leeuwin Current at 32°S. Journal of Geophysical Research 108, 19–39.
| Annual and interannual variations of the Leeuwin Current at 32°S.Crossref | GoogleScholarGoogle Scholar |
Feng, M, Wijffels, S, Godfrey, S, and Meyers, G (2005). Do eddies play a role in the momentum balance of the Leeuwin Current? Journal of Physical Oceanography 35, 964–975.
| Do eddies play a role in the momentum balance of the Leeuwin Current?Crossref | GoogleScholarGoogle Scholar |
Feng, M, Slawinski, D, Beckley, LE, and Keesing, JK (2010). Retention and dispersal of shelf waters influenced by interactions of ocean boundary current and coastal geography. Marine and Freshwater Research 61, 1259–1267.
| Retention and dispersal of shelf waters influenced by interactions of ocean boundary current and coastal geography.Crossref | GoogleScholarGoogle Scholar |
Feng, M, Caputi, N, Penn, J, Slawinski, D, de Lestang, S, Weller, E, and Pearce, A (2011). Ocean circulation, stokes drift, and connectivity of western rock lobster (Panulirus cygnus) population. Canadian Journal of Fisheries and Aquatic Sciences 68, 1182–1196.
| Ocean circulation, stokes drift, and connectivity of western rock lobster (Panulirus cygnus) population.Crossref | GoogleScholarGoogle Scholar |
Fitzgibbon, QP, Jeffs, AG, and Battaglene, SC (2014). The Achilles heel for spiny lobsters: the energetics of the non-feeding post-larval stage. Fish and Fisheries 15, 312–326.
| The Achilles heel for spiny lobsters: the energetics of the non-feeding post-larval stage.Crossref | GoogleScholarGoogle Scholar |
Gersbach, GH, Pattiaratchi, CB, Ivey, GN, and Cresswell, GR (1999). Upwelling on the south-west coast of Australia: source of the Capes Current. Continental Shelf Research 19, 363–400.
| Upwelling on the south-west coast of Australia: source of the Capes Current.Crossref | GoogleScholarGoogle Scholar |
Goldstein, J, and Butler, MJ (2009). Behavioral enhancement of onshore transport by postlarval Caribbean spiny lobster (Panulirus argus). Limnology and Oceanography 54, 1669–1678.
| Behavioral enhancement of onshore transport by postlarval Caribbean spiny lobster (Panulirus argus).Crossref | GoogleScholarGoogle Scholar |
Greenwood, JE, Feng, M, and Waite, AM (2007). A one-dimensional simulation of biological production in two contrasting mesoscale eddies in the south eastern Indian Ocean. Deep-Sea Research – II. Topical Studies in Oceanography 54, 1029–1044.
| A one-dimensional simulation of biological production in two contrasting mesoscale eddies in the south eastern Indian Ocean.Crossref | GoogleScholarGoogle Scholar |
Griffin, DA, Wilkin, JL, Chubb, CF, Pearce, AF, and Caputi, N (2001). Ocean currents and the larval phase of Australian western rock lobster, Panulirus cygnus. Marine and Freshwater Research 52, 1187–1199.
| Ocean currents and the larval phase of Australian western rock lobster, Panulirus cygnus.Crossref | GoogleScholarGoogle Scholar |
Hamon, BV (1965). Geostrophic currents in the southeastern Indian Ocean. Australian Journal of Marine and Freshwater Research 16, 255–271.
| Geostrophic currents in the southeastern Indian Ocean.Crossref | GoogleScholarGoogle Scholar |
Hawkins, S, Allcock, A, Bates, A, Firth, L, Smith, I, Swearer, S, and Todd, P (2019). Oceanography and Marine Biology: Preface. Oceanography and Marine Biology: An Annual Review 57, 325–356.
Hill, AE, Brown, J, and Fernand, L (1996). The western Irish Sea gyre: a retention system for Norway lobster (Nephrops norvegicus)? Oceanologica Acta 19, 357–368.
Hinojosa, IA, Green, BS, Gardner, C, Hesse, J, Stanley, JA, and Jeffs, AG (2016). Reef sound as an orientation cue for shoreward migration by pueruli of the rock lobster, Jasus edwardsii. PLoS ONE 11, e0157862.
| Reef sound as an orientation cue for shoreward migration by pueruli of the rock lobster, Jasus edwardsii.Crossref | GoogleScholarGoogle Scholar |
Hinojosa, IA, Gardner, C, Green, BS, Jeffs, A, Leon, R, and Linnane, A (2017). Differing environmental drivers of settlement across the range of southern rock lobster (Jasus edwardsii) suggest resilience of the fishery to climate change. Fisheries Oceanography 26, 49–64.
| Differing environmental drivers of settlement across the range of southern rock lobster (Jasus edwardsii) suggest resilience of the fishery to climate change.Crossref | GoogleScholarGoogle Scholar |
Huang, Z, and Feng, M (2015). Remotely sensed spatial and temporal variability of the Leeuwin Current using MODIS data. Remote Sensing of Environment 166, 214–232.
| Remotely sensed spatial and temporal variability of the Leeuwin Current using MODIS data.Crossref | GoogleScholarGoogle Scholar |
Jeffs, AG, Montgomery, JC, and Tindle, CT (2005). How do spiny lobster post-larvae find the coast? New Zealand Journal of Marine and Freshwater Research 39, 605–617.
| How do spiny lobster post-larvae find the coast?Crossref | GoogleScholarGoogle Scholar |
Jones, S, Inall, M, Porter, M, Graham, JA, and Cottier, F (2020). Storm-driven across-shelf oceanic flows into coastal waters. Ocean Science 16, 389–403.
| Storm-driven across-shelf oceanic flows into coastal waters.Crossref | GoogleScholarGoogle Scholar |
Kolbusz, J, de Lestang, S, Langlois, T, and Pattiaratchi, C (2021a). Changes in Panulirus cygnus settlement along Western Australia using a long time series. Frontiers in Marine Science 8, 628912.
| Changes in Panulirus cygnus settlement along Western Australia using a long time series.Crossref | GoogleScholarGoogle Scholar |
Kolbusz, J, Langlois, T, Pattiaratchi, C, and de Lestang, S (2021b). Using an oceanographic model to investigate the mystery of the missing puerulus. Biogeosciences 19, 517–539.
| Using an oceanographic model to investigate the mystery of the missing puerulus.Crossref | GoogleScholarGoogle Scholar |
Kvile, KØ, Romagnoni, G, Dagestad, K-F, Langangen, Ø, and Kristiansen, T (2018). Sensitivity of modelled North Sea cod larvae transport to vertical behaviour, ocean model resolution and interannual variation in ocean dynamics. ICES Journal of Marine Science 75, 2413–2424.
| Sensitivity of modelled North Sea cod larvae transport to vertical behaviour, ocean model resolution and interannual variation in ocean dynamics.Crossref | GoogleScholarGoogle Scholar |
Lemmens, JWTJ (1994). Biochemical evidence for absence of feeding in puerulus larvae of the Western rock lobster Panulirus cygnus (Decapoda: Palinuridae). Marine Biology 118, 383–391.
| Biochemical evidence for absence of feeding in puerulus larvae of the Western rock lobster Panulirus cygnus (Decapoda: Palinuridae).Crossref | GoogleScholarGoogle Scholar |
Lett, C, Verley, P, Mullon, C, Parada, C, Brochier, T, Penven, P, and Blanke, B (2008). A Lagrangian tool for modelling ichthyoplankton dynamics. Environmental Modelling & Software 23, 1210–1214.
| A Lagrangian tool for modelling ichthyoplankton dynamics.Crossref | GoogleScholarGoogle Scholar |
Liddy, GC, Phillips, BF, and Maguire, GB (2004). Effects of temperature and food density on the survival and growth of early stage phyllosoma of the western rock lobster, Panulirus cygnus. Aquaculture 242, 207–215.
| Effects of temperature and food density on the survival and growth of early stage phyllosoma of the western rock lobster, Panulirus cygnus.Crossref | GoogleScholarGoogle Scholar |
Linnane, A, Gardner, C, Hobday, D, Punt, A, McGarvey, R, Feenstra, J, Matthews, J, and Green, B (2010). Evidence of large-scale spatial declines in recruitment patterns of southern rock lobster Jasus edwardsii, across south-eastern Australia. Fisheries Research 105, 163–171.
| Evidence of large-scale spatial declines in recruitment patterns of southern rock lobster Jasus edwardsii, across south-eastern Australia.Crossref | GoogleScholarGoogle Scholar |
Melville-Smith, R, De Lestang, S, and Thomson, AW (2009). Spatial and temporal changes in egg production in the western rock lobster (Panulirus cygnus) fishery. New Zealand Journal of Marine and Freshwater Research 43, 151–161.
| Spatial and temporal changes in egg production in the western rock lobster (Panulirus cygnus) fishery.Crossref | GoogleScholarGoogle Scholar |
Meuleners, MJ, Ivey, GN, and Pattiaratchi, CB (2008). A numerical study of the eddying characteristics of the Leeuwin Current System. Deep-Sea Research – I. Oceanographic Research Papers 55, 261–276.
| A numerical study of the eddying characteristics of the Leeuwin Current System.Crossref | GoogleScholarGoogle Scholar |
Moore, TS, Matear, RJ, Marra, J, and Clementson, L (2007). Phytoplankton variability off the Western Australian Coast: mesoscale eddies and their role in cross-shelf exchange. Deep-Sea Research – II. Topical Studies in Oceanography 54, 943–960.
| Phytoplankton variability off the Western Australian Coast: mesoscale eddies and their role in cross-shelf exchange.Crossref | GoogleScholarGoogle Scholar |
Mullaney, TJ, and Suthers, IM (2013). Entrainment and retention of the coastal larval fish assemblage by a short-lived, submesoscale, frontal eddy of the East Australian Current. Limnology and Oceanography 58, 1546–1556.
| Entrainment and retention of the coastal larval fish assemblage by a short-lived, submesoscale, frontal eddy of the East Australian Current.Crossref | GoogleScholarGoogle Scholar |
Nolasco, R, Gomes, I, Peteiro, L, Albuquerque, R, Luna, T, Dubert, J, Swearer, SE, and Queiroga, H (2018). Independent estimates of marine population connectivity are more concordant when accounting for uncertainties in larval origins. Scientific Reports 8, 2641.
| Independent estimates of marine population connectivity are more concordant when accounting for uncertainties in larval origins.Crossref | GoogleScholarGoogle Scholar |
Paris C, Vaz A, Helgers J, Wood S, Ross R (2017) Connectivity Modeling System User’s Guide. Available at https://github.com/beatrixparis/connectivity-modeling-system/blob/master/User-Guide-v2.pdf
Pattiaratchi, C, and Buchan, S (1991). Implications of long-term climate change for the Leeuwin Current. Journal of the Royal Society of Western Australia 74, 133–140.
Pattiaratchi C, Hetzel Y (2020) Sea surface temperature variability: state and trends of Australia’s oceans report. Integrated Marine Observing System, Hobart, Tas., Australia.
Pattiaratchi C, Siji P (2020) Variability in ocean currents around Australia: State and trends of Australia’s oceans report. Integrated Marine Observing System.
Pattiaratchi, C, and Woo, M (2009). The mean state of the Leeuwin Current system between North West Cape and Cape Leeuwin. Journal of the Royal Society of Western Australia 92, 221–241.
Pattiaratchi, C, Hegge, B, Gould, J, and Eliot, I (1997). Impact of sea-breeze activity on nearshore and foreshore processes in southwestern Australia. Continental Shelf Research 17, 1539–1560.
| Impact of sea-breeze activity on nearshore and foreshore processes in southwestern Australia.Crossref | GoogleScholarGoogle Scholar |
Pearce, AF (1991). Eastern boundary currents of the southern hemisphere. Journal of the Royal Society of Western Australia 74, 35–45.
Pearce, AF, and Griffiths, RW (1991). The mesoscale structure of the Leeuwin Current: a comparison of laboratory models and satellite imagery. Journal of Geophysical Research 96, 16739–16757.
| The mesoscale structure of the Leeuwin Current: a comparison of laboratory models and satellite imagery.Crossref | GoogleScholarGoogle Scholar |
Pearce, A, and Pattiaratchi, C (1999). The Capes Current: a summer countercurrent flowing past Cape Leeuwin and Cape Naturaliste, Western Australia. Continental Shelf Research 19, 401–420.
| The Capes Current: a summer countercurrent flowing past Cape Leeuwin and Cape Naturaliste, Western Australia.Crossref | GoogleScholarGoogle Scholar |
Pearce, AF, and Phillips, BF (1988). Enso events, the leeuwin current, and larval recruitment of the western rock lobster. ICES Journal of Marine Science 45, 13–21.
| Enso events, the leeuwin current, and larval recruitment of the western rock lobster.Crossref | GoogleScholarGoogle Scholar |
Peliz, A, Marchesiello, P, Dubert, J, Marta-Almeida, M, Roy, C, and Queiroga, H (2007). A study of crab larvae dispersal on the Western Iberian Shelf: physical processes. Journal of Marine Systems 68, 215–236.
| A study of crab larvae dispersal on the Western Iberian Shelf: physical processes.Crossref | GoogleScholarGoogle Scholar |
Phillips, BF, and McWilliam, PS (2009). Spiny lobster development: where does successful metamorphosis to the puerulus occur? A review. Reviews in Fish Biology and Fisheries 19, 193–215.
| Spiny lobster development: where does successful metamorphosis to the puerulus occur? A review.Crossref | GoogleScholarGoogle Scholar |
Phillips, BF, and Olsen, L (1975). Swimming behaviour of the puerulus larvae of the western rock lobster. Marine and Freshwater Research 26, 415–417.
| Swimming behaviour of the puerulus larvae of the western rock lobster.Crossref | GoogleScholarGoogle Scholar |
Phillips, BF, Rimmer, DW, and Reid, DD (1978). Ecological investigations of the late-stage phyllosoma and puerulus larvae of the western rock lobster Panulirus longipes cygnus. Marine Biology 45, 347–357.
| Ecological investigations of the late-stage phyllosoma and puerulus larvae of the western rock lobster Panulirus longipes cygnus.Crossref | GoogleScholarGoogle Scholar |
Phillips, BF, Brown, PA, Rimmer, DW, and Braine, SJ (1981). Description, distribution and abundance of late larval stages of the Scyllaridae (slipper lobsters) in the south-eastern Indian Ocean. Australian Journal of Marine and Freshwater Research 32, 417–437.
| Description, distribution and abundance of late larval stages of the Scyllaridae (slipper lobsters) in the south-eastern Indian Ocean.Crossref | GoogleScholarGoogle Scholar |
Pires, RFT, Peliz, Á, Pan, M, and dos Santos, A (2020). ‘There and back again’ − how decapod megalopae find the way home: a modelling exercise for Pachygrapsus marmoratus. Progress in Oceanography 184, 102331.
| ‘There and back again’ − how decapod megalopae find the way home: a modelling exercise for Pachygrapsus marmoratus.Crossref | GoogleScholarGoogle Scholar |
Rafiq, S, Pattiaratchi, C, and Janeković, I (2020). Dynamics of the land–sea breeze system and the surface current response in South-West Australia. Journal of Marine Science and Engineering 8, 931.
| Dynamics of the land–sea breeze system and the surface current response in South-West Australia.Crossref | GoogleScholarGoogle Scholar |
Rennie, S, Pattiaratchi, C, and Mccauley, R (2007). Eddy formation through the interaction between the Leeuwin Current, Leeuwin Undercurrent and topography. Deep-Sea Research – I. Oceanographic Research Papers 54, 818–836.
| Eddy formation through the interaction between the Leeuwin Current, Leeuwin Undercurrent and topography.Crossref | GoogleScholarGoogle Scholar |
Rimmer, DW, and Phillips, BF (1979). Diurnal migration and vertical distribution of phyllosoma larvae of the western rock lobster Panulirus cygnus. Marine Biology 54, 109–124.
| Diurnal migration and vertical distribution of phyllosoma larvae of the western rock lobster Panulirus cygnus.Crossref | GoogleScholarGoogle Scholar |
Robson, NA, Hetzel, Y, Whiting, S, Wijeratne, S, Pattiaratchi, CB, Withers, P, and Thums, M (2017). Use of particle tracking to determine optimal release dates and locations for rehabilitated neonate sea turtles. Frontiers in Marine Science 4, 173.
| Use of particle tracking to determine optimal release dates and locations for rehabilitated neonate sea turtles.Crossref | GoogleScholarGoogle Scholar |
Roughan, M, Macdonald, HS, Baird, ME, and Glasby, TM (2011). Modelling coastal connectivity in a Western Boundary Current: seasonal and inter-annual variability. Deep-Sea Research – II. Topical Studies in Oceanography 58, 628–644.
| Modelling coastal connectivity in a Western Boundary Current: seasonal and inter-annual variability.Crossref | GoogleScholarGoogle Scholar |
Saunders, MI, Thompson, PA, Jeffs, AG, Säwström, C, Sachlikidis, N, Beckley, LE, and Waite, AM (2012). Fussy feeders: phyllosoma larvae of the Western Rocklobster (Panulirus cygnus) demonstrate prey preference. PLoS ONE 7, e36580.
| Fussy feeders: phyllosoma larvae of the Western Rocklobster (Panulirus cygnus) demonstrate prey preference.Crossref | GoogleScholarGoogle Scholar |
Säwström, C, Beckley, LE, Saunders, MI, Thompson, PA, and Waite, AM (2014). The zooplankton prey field for rock lobster phyllosoma larvae in relation to oceanographic features of the south-eastern Indian Ocean. Journal of Plankton Research 36, 1003–1016.
| The zooplankton prey field for rock lobster phyllosoma larvae in relation to oceanographic features of the south-eastern Indian Ocean.Crossref | GoogleScholarGoogle Scholar |
Schiller, A, Oke, PR, Brassington, G, Entel, M, Fiedler, R, Griffin, DA, and Mansbridge, JV (2008). Eddy-resolving ocean circulation in the Asian–Australian region inferred from an ocean reanalysis effort. Progress in Oceanography 76, 334–365.
| Eddy-resolving ocean circulation in the Asian–Australian region inferred from an ocean reanalysis effort.Crossref | GoogleScholarGoogle Scholar |
Schilling, HT, Everett, JD, Smith, JA, Stewart, J, Hughes, JM, Roughan, M, Kerry, C, and Suthers, IM (2020). Multiple spawning events promote increased larval dispersal of a predatory fish in a western boundary current. Fisheries Oceanography 29, 309–323.
| Multiple spawning events promote increased larval dispersal of a predatory fish in a western boundary current.Crossref | GoogleScholarGoogle Scholar |
Singh, SP, Groeneveld, JC, Hart-Davis, MG, Backeberg, BC, and Willows-Munro, S (2018). Seascape genetics of the spiny lobster Panulirus homarus in the Western Indian Ocean: understanding how oceanographic features shape the genetic structure of species with high larval dispersal potential. Ecology and Evolution 8, 12221–12237.
| Seascape genetics of the spiny lobster Panulirus homarus in the Western Indian Ocean: understanding how oceanographic features shape the genetic structure of species with high larval dispersal potential.Crossref | GoogleScholarGoogle Scholar |
Smith, RL, Huyer, A, Godfrey, S, and Church, JA (1991). The Leeuwin Current off Western Australia, 1986–1987. Journal of Physical Oceanography 21, 323–345.
| The Leeuwin Current off Western Australia, 1986–1987.Crossref | GoogleScholarGoogle Scholar |
Stanley, JA, Hesse, J, Hinojosa, IA, and Jeffs, AG (2015). Inducers of settlement and moulting in post-larval spiny lobster. Oecologia 178, 685–697.
| Inducers of settlement and moulting in post-larval spiny lobster.Crossref | GoogleScholarGoogle Scholar |
Steele, RW, and Neuheimer, AB (2022). Assessing the ability of the growing degree-day metric to explain variation in size-at-age and duration-to-moult of lobsters and crabs. Canadian Journal of Fisheries and Aquatic Sciences 79, 850–860.
| Assessing the ability of the growing degree-day metric to explain variation in size-at-age and duration-to-moult of lobsters and crabs.Crossref | GoogleScholarGoogle Scholar |
Telcik, N, and Pattiaratchi, C (2014). Influence of Northwest Cloudbands on Southwest Australian Rainfall. Journal of Climatology 2014, 671394.
| Influence of Northwest Cloudbands on Southwest Australian Rainfall.Crossref | GoogleScholarGoogle Scholar |
van Sebille, E, Griffies, S, Abernathey, R, Adams, T, Berloff, P, Biastoch, A, Blanke, B, Chassignet, E, Cheng, Y, Cotter, C, Deleersnijder, E, Döös, K, Drake, H, Drijfhout, S, Gary, S, Heemink, A, Kjellsson, J, Koszalka, I, Lange, M, et al. (2018). Lagrangian ocean analysis: fundamentals and practices. Ocean Modelling 121, 49–75.
| Lagrangian ocean analysis: fundamentals and practices.Crossref | GoogleScholarGoogle Scholar |
Waite, AM, Raes, E, Beckley, LE, Thompson, PA, Griffin, D, Saunders, M, Säwström, C, O’Rorke, R, Wang, M, Landrum, J, and Jeffs, A (2019). Production and ecosystem structure in cold-core vs. warm-core eddies: implications for the zooplankton isoscape and rock lobster larvae. Limnology and Oceanography 64, 2405–2423.
| Production and ecosystem structure in cold-core vs. warm-core eddies: implications for the zooplankton isoscape and rock lobster larvae.Crossref | GoogleScholarGoogle Scholar |
Wang, M, O’Rorke, R, Waite, AM, Beckley, LE, Thompson, P, and Jeffs, AG (2015). Condition of larvae of western rock lobster (Panulirus cygnus) in cyclonic and anticyclonic eddies of the Leeuwin Current off Western Australia. Marine and Freshwater Research 66, 1158–1167.
| Condition of larvae of western rock lobster (Panulirus cygnus) in cyclonic and anticyclonic eddies of the Leeuwin Current off Western Australia.Crossref | GoogleScholarGoogle Scholar |
Wijeratne, S, Pattiaratchi, C, and Proctor, R (2018). Estimates of surface and subsurface boundary current transport around Australia. Journal of Geophysical Research: Oceans 123, 3444–3466.
| Estimates of surface and subsurface boundary current transport around Australia.Crossref | GoogleScholarGoogle Scholar |
Wilkin, JL, and Jeffs, AG (2011). 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 1, 163–175.
| 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?Crossref | GoogleScholarGoogle Scholar |
Woo, M, and Pattiaratchi, C (2008). Hydrography and water masses off the western Australian coast. Deep-Sea Research – I. Oceanographic Research Papers 55, 1090–1104.
| Hydrography and water masses off the western Australian coast.Crossref | GoogleScholarGoogle Scholar |
Woodings, LN, Murphy, NP, Jeffs, A, Suthers, IM, Liggins, G, and Strugnell, JM (2019). Distribution of Palinuridae and Scyllaridae phyllosoma larvae within the East Australian current: a climate change hot spot. Marine and Freshwater Research 70, 1020–1033.
| Distribution of Palinuridae and Scyllaridae phyllosoma larvae within the East Australian current: a climate change hot spot.Crossref | GoogleScholarGoogle Scholar |
Yeung, C, Jones, DL, Criales, MM, Jackson, TL, and Richards, WJ (2001). Influence of coastal eddies and counter-currents on the influx of spiny lobster, Panulirus argus, postlarvae into Florida Bay. Marine And Freshwater Research 52, 1217–1232.
| Influence of coastal eddies and counter-currents on the influx of spiny lobster, Panulirus argus, postlarvae into Florida Bay.Crossref | GoogleScholarGoogle Scholar |