- Open Access
Identifying the key biophysical drivers, connectivity outcomes, and metapopulation consequences of larval dispersal in the sea
Movement Ecologyvolume 3, Article number: 17 (2015)
Population connectivity, which is essential for the persistence of benthic marine metapopulations, depends on how life history traits and the environment interact to influence larval production, dispersal and survival. Although we have made significant advances in our understanding of the spatial and temporal dynamics of these individual processes, developing an approach that integrates the entire population connectivity process from reproduction, through dispersal, and to the recruitment of individuals has been difficult.
We present a population connectivity modelling framework and diagnostic approach for quantifying the impact of i) life histories, ii) demographics, iii) larval dispersal, and iv) the physical seascape, on the structure of connectivity and metapopulation dynamics. We illustrate this approach using the subtidal rocky reef ecosystem of Port Phillip Bay, were we provide a broadly-applicable framework of population connectivity and quantitative methodology for evaluating the relative importance of individual factors in determining local and system outcomes.
The spatial characteristics of marine population connectivity are primarily influenced by larval mortality, the duration of the pelagic larval stage, and the settlement competency characteristics, with significant variability imposed by the geographic setting and the timing of larval release. The relative influence and the direction and strength of the main effects were strongly consistent among 10 connectivity-based metrics.
These important intrinsic factors (mortality, length of the pelagic larval stage, and the extent of the precompetency window) and the spatial and temporal variability represent key research priorities for advancing our understanding of the connectivity process and metapopulation outcomes.
Understanding the causes and consequences of dispersal is a foundational goal in population and community ecology, and evolution. Dispersal, or the exchange of individuals among natal and non-natal sites, is the primary process that ‘connects’ populations, with important impacts on local demography, landscape-wide population dynamics and gene flow. Although there are a diversity of evolutionary drivers of dispersal, such as reducing kin competition and inbreeding , ‘bet hedging’ offspring success through spatial-temporal variation in reproduction , and escaping unfavourable or ephemeral local conditions , it is the more proximate drivers of dispersal at ecological time scales which influence the selective pressures on dispersal traits . At ecological time scales, population connectivity is critical for population growth [5, 6] and system persistence , aids in the local-scale recovery or rescue from severe disturbances , and plays a key role in driving metapopulation dynamics . Connectivity is also believed to be important in determining how species will cope in a changing climate by allowing species’ ranges to expand or constrict in response to environmental shifts .
The importance of population connectivity has also fuelled efforts to integrate this process into conservation and management, particularly in relation to population persistence and viability , reserve design and spatial management strategies , and the alignment between management efforts and ecological processes . Unfortunately, identifying quantitative conservation objectives with respect to connectivity remains difficult, particularly in regards to the specific process, landscape feature, or population outcome being targeted for conservation . For example, targeting strong source habitat patches (self-sustaining subpopulations that are net exporters of individuals) would lead to different conservation outcomes than if one were to target locally persistent patches, critical stepping-stones, or sites receiving the highest diversity of settlers . Clearly, a more holistic understanding of the drivers of population connectivity and the population and system-wide implications is needed.
Recent advances in movement ecology, population biology, and land/seascape ecology have improved our understanding of this biological-physical process. The study of the mechanisms of dispersal movements per se has elevated the importance of both intrinsic (biological) and extrinsic (e.g., environmental) drivers [16, 17] and the role of individual-based decisions in determining dispersal outcomes [18, 19]. This work appears to be coalescing into the study of dispersal syndromes [20, 21], or the patterns of covariance between dispersal potential and suites of life history traits based on shared evolutionary or environmental histories . Concurrently, the exploration of the patterns and implications of dispersal has often taken a spatially-explicit or landscape ecological approach, which has increased our ability to quantify the impact of matrix structure and habitat topology on connectivity and how it interacts with life history traits in determining broad-scale emergent patterns of metapopulation dynamics . The mechanistic approaches to studying population connectivity have not been well integrated with landscape ecological approaches, largely due to the differences in spatiotemporal scales. Unifying these often disparate approaches, however, would provide a more holistic quantitative framework for investigating the factors that influence dispersal-driven connectivity and their population level consequences.
Studying dispersal-driven population connectivity in benthic marine species with complex life cycles is particularly challenging due to the strong influence of currents, the age, size, and behavioural complexities of the dispersing individuals, and the spatiotemporal scales (and variability) of the process . Despite these difficulties, technological advances in larval tagging, computation power and model sophistication have enabled significant progress in estimating the scales of connectivity and identifying several key drivers for a number of taxa . At local scales, the proportion of total larvae released that ultimately recruit back to the natal population, termed local larval retention [26, 27], is essential for determining demographically meaningful estimates of population replenishment  and quantifying a population’s dependence on subsidies from non-local sources for population persistence . Local retention is driven by local-scale hydrodynamics or ‘sticky water’ , early-stage larval behaviour , and aspects of the local habitat structure , with estimates as high as 20 to 30 % in some systems [33, 34]. Together, local retention and the amount of immigrating larvae arriving from upstream sources determine the relative dependencies on natal and non-natal larvae to population growth and persistence. Although measuring this mixture between natal and non-natal recruits is becoming easier with genetic  and otolith-based  approaches, estimates must be viewed in the context of the local population size, reproductive output, and local demographic rates to determine relevant recruitment rates .
Although individual biophysical parameters are important in determining connectivity outcomes, such as larval mortality , larval behaviour  and sensing , reproductive output , duration of the pelagic stage , and local-scale ocean physics [39–41], we have very little understanding of the relative importance and interactions of these parameters in any given system or for any particular species (but see [40, 32]). To move beyond one-at-a-time empirical evaluations of parameter importance, we are largely dependent on models to develop a more comprehensive understanding of this complex system across scales .
Here we have taken a system-level perspective and define metapopulation connectivity as the aggregate process integrating natal dispersal, post-settlement survival and reproduction (i.e., recruitment), in both natal and non-natal sites. This process-based conceptualisation includes four-stages of population connectivity (Fig. 1) and provides a clear framework for investigating the primary intrinsic and extrinsic drivers of marine larval dispersal  and the local to population-wide consequences. This connectivity definition is consistent with the recent marine literature  and genetic descriptions , and incorporates the three phases of dispersal common in the movement ecology field [45, 20, 46, 18]. Guided by this framework, we have used a well-validated, high-resolution, and three-dimensional biophysical model of marine larval dispersal to gain a better understanding of the drivers of marine population connectivity and its metapopulation implications. Specifically, we have three primary aims: i) to present a process-based conceptual framework of marine population connectivity and their intrinsic and extrinsic drivers, ii) to quantify the relative importance of individual parameters in determining population connectivity outcomes, and iii) to identify important knowledge gaps and prioritise research questions to improve our understanding of population connectivity.
Population connectivity framework
Due to the nearly passive dispersal qualities of the early developmental stages, the strong influence and dynamics of the physical environment (i.e., currents), and the potential mobility and sensing capacity of late-stage dispersers, the existing mechanistic models of dispersal are inadequate. As a result, we have developed a framework describing the four stages of population connectivity, incorporating the three stages of dispersal [47, 20] appropriate for both vector-mediated passive dispersal  and active dispersers, and includes the final stage of recruitment (post-settlement survival to reproduction), which is critical in determining ‘realised’ connectivity  and important in evolutionary models of dispersal . This framework illustrates the key intrinsic and extrinsic drivers acting upon each stage of connectivity and captures the complex bio-physical and contextual interdependencies of this process characteristic of marine and aquatic environments (Fig. 1).
The first stage in the connectivity process is the initiation of emigration (‘departure’ in , ‘initiation’ in ) in which some quantity of gametes, spores, or larvae are released from the parent. Once released, this stage is followed by the transport and movement stage [‘transience’ in 20, ‘transport’ in 46] where the disperser’s trajectory is determined both by the potential advection and turbulence of currents, and the motility and behaviour of individuals, often extending from days to weeks, with wide variability among taxa . Settlement marks the end of the dispersal period [‘termination’ in 46, ‘settlement’ in 20] in which dispersers actively settle to some suitable habitat patch, either within the natal source site or in a non-natal location. Individuals successfully settling into viable habitat enter the final stage of recruitment, in which some may survive and mature to reproduce, thereby contributing to subpopulation demographics and gene flow. Together, these four stages of population connectivity represent the unique biophysical processes determining the connectivity of subpopulations.
The parameters important in emigration are related to reproductive output [or vector seed load in 46], and include the fecundity, abundance, and fertilisation success of the parents . Reproductive output is strongly context-dependent and influenced by extrinsic drivers such as the quality, quantity, and spatial structure of the natal habitat, and their effect on individual parents [phenotype-dependence, 20]. The transport and movement stage depends on mortality , larval development , individual sensing and motility , and the extrinsic role of currents  and spatial habitat structure . Settlement can be as biophysically complex and governed by intrinsic (settlement competency window, behaviour, motility, and sensing) and extrinsic factors (habitat quality and structure and currents) and is likely context- (phenotype-environment mismatch, ) and condition-dependent [54, 55]. Similarly, once settled, recruitment into the adult stage is determined by individual growth and survival to maturation, which are influenced by habitat quality, competition, and individual condition. This conceptual model of population connectivity effectively partitions the primary intrinsic and extrinsic parameters among the key life stages enabling the relative importance of each parameter and stage to be quantified.
Test case: bay-wide marine population connectivity
Port Phillip Bay (PPB) is a large (~2,000 km2) semi-enclosed temperate marine system in Victoria, Australia, and is ideally suited for exploring the biophysics of marine population connectivity. One of the most prominent and economically important habitat features in the bay are the sub-tidal rocky reefs (Fig. 2). These reefs harbour the vast majority of the Bay’s biodiversity  and are important to commercial and recreational fisheries . The rocky reef habitat is restricted to the shallow periphery of the Bay with each discrete reef isolated by a matrix of unconsolidated mud and sand. For the purpose of this study, we used a subset of eight individual reefs (Fig. 2) to explore the relative importance of individual biological and physical drivers to population connectivity.
The hydrodynamics of PPB have been well studied and accurately modelled since the mid 1980′s  accompanied by a long history of in situ measurements (salinity, temperature, etc.) and ecosystem monitoring , resulting in an accurate and well validated representation of the Bay’s dynamics. For this study, we used a PPB-wide 3-dimensional hydrodynamic model (400 m horizontal resolution, 8 vertical layers, and hourly time-steps, ) incorporating wind, sea level, temperature, air pressure, tides, and solar radiation forcing . The hydrodynamics, together with high-resolution rocky reef habitat data , represent PPB’s physical domain, within which the simulated individual dispersers interact.
An individual-based model of dispersal (modified from ) was used to simulate the movement and settlement of dispersers throughout PPB. Although this numerical approach has been used successfully in quantifying passive fish dispersal with simplified biology (e.g., [64, 65]), we have made several novel advancements to more accurately represent the reality of key biological traits typical of the larval life histories of marine invertebrates and fishes. First, we have implemented a larval growth scheme in which all biological parameters related to development are capable of changing through time, allowing individuals to change behaviour and sensing capacity, as well as become less susceptible to physiological stress (resulting in decreased mortality over time), for example. Second, we have added a suite of behavioural schemes to match the known and hypothesised movement patterns (e.g., diel vertical migration, homing behaviour) of marine taxa . All key biological parameters (Table 1) important to population connectivity and their function in the growth and behavioural schemes are described below.
At the initiation of emigration, the reproductive output (RO) is determined by the fecundity and abundance of the modelled taxon, and the habitat characteristics of the natal site. While individuals move through the environment (Transport and movement stage, Fig. 1), their trajectories and success are governed by mortality, growth and development, and individual-based behavioural parameters. Mortality was modelled as a Weibull function with a rate (A) and shape (B) parameter, capturing a full range of empirically-based mortality functions (e.g., ) due to processes such as natural mortality and starvation. This function has the capacity to represent the exponential decay function common to most dispersal models (if B = 1, then the Weibull reduces to the exponential function), and the flexibility to represent a ‘fat tail’ dispersal kernel, which is perhaps more appropriate for many taxa . Individual growth and development is controlled through four parameters: 1) precompetency period (PreP) in which individuals are not physiologically capable of settlement; 2) competency rate (ComR) describing the developmental transition to competency; 3) development period (DevP) as a proportion of the precompetency period within which the individuals are passive dispersers before the onset of active behaviour; and 4) the initial fall velocity (iFV) describing the buoyancy of the individuals during the early development period. The flexibility in this growth scheme has the capacity to represent a broad range of taxa, including drifting seaweeds, slowly developing fish larvae and quickly developing invertebrate larvae (e.g., species with non-feeding larvae). Following the developmental period, individuals have the capacity to move, implemented through three potential behavioural strategies: 1) passive strategy where the individuals are transported by currents only; 2) benthic seeking strategy where individuals actively and constantly swim (with speed, Sp*0.05) to a specific target depth (TD); and 3) diel vertical migration where individuals actively swim (at speed, Sp*0.05) to a target depth (TD) only during daylight hours. Concurrently with these vertical swimming strategies, all non-passive individuals have the capacity to sense and swim (at speed, Sp) to nearby habitat patches at a given detection limit or homing distance (HmD) from individual reefs, to simulate the typical distances over which larvae have the potential to detect and orientate to benthic settlement habitat using auditory and/or olfactory cues . Finally, following successful transport and movement through the environment, competent individuals may settle to suitable reef habitat patches and a proportion of these survive to transition to the final stage of recruitment. Post-settlement mortality rates in benthic marine organisms are highly variable, but can be as high as 96 % within the first 24 h . As such, post-settlement survival was determined by a simple survival probability (Sr), allocated randomly across patches to represent (unknown) recruitment costs .
Simulations were completed for each of eight selected reef habitat patches to explore the impact of geographic location on connectivity outcomes. In addition, simulations were initiated on three separate dates (rls, Table 1) to capture temporal variability across three seasons and representing potential spawning dates for different taxa, while maintaining computational feasibility. In this way, we were able to quantify the relative impact of both intrinsic (e.g., growth, behaviour) and extrinsic (currents, habitat quality) factors, and their geographic and temporal signatures, on 10 different ecologically-relevant metrics of population connectivity that span three scales and a range of questions. At the local scale we calculated the proportion of individuals released at initiation that recruit back to the natal habitat patch, termed local retention (LR), as well as the proportion of successfully recruiting individuals that originated from the focal patch, or self-recruitment (SR). In addition, we calculated the diversity of successful settlers (H’) to each habitat patch, which may influence population persistence and long-term resilience . At the scale of downstream connectivity, we quantified several distance-based measures of population connectivity, including the median geographic distance displaced by individuals (mdG) and the maximum distance (mxG), as well as the total proportion of successful settlers to all downstream patches (S). The downstream contributions from each source (i.e., source strength) was quantified by counting the number of downstream linkages (dC) and calculating the weighted degree centrality (Cw), a network-based measure characterising the strength and evenness of downstream linkages . Finally, at the eight-patch metapopulation scale, we calculated the metapopulation growth rate [λM, 6] and the metapopulation capacity [λmax, 5] (see Table 2 for additional details on the ecological significance of each metric).
To quantify the main effects and all interactions, a global sensitivity analysis (GSA) was used [69, 70]. Due to the computational requirements and complexity of the PPB connectivity model, we performed a non-parametric regression tree GSA on a meta-model [69, 71] derived from a full suite of input parameters to quantify parameter importance (R package CompModSA with ‘sensitivity’ function). To aid in the interpretation of these importance values, and to effectively quantify the main effect and direction of influence of parameters on model output, we paired this with a generalized linear regression (GLM) analysis . For the GLM, we calculated the main effects on the standardized data using the identity link function. Sensitivities where visualized by plotting the effect of one standard deviation change in each model parameter on the response . For computational feasibility, we ran the connectivity model for a suite of parameter values generated using an optimum Latin Hypercube Sample (LHS) scheme , assuming each parameter has a uniform distribution. The LHS scheme ensures that the entire multivariate range is sampled and the full behaviour of the model explored. The PPB connectivity model was used to simulate 21,600 scenarios, representing 900 unique parameter combinations (each containing values for all 13 parameters) for each of eight habitat patches across three release dates (seasons). The complete ensemble of simulations consisted of 2,583 million dispersers being tracked (~130,000 individuals per simulation) for the full sensitivity analysis. The parameter importance was based on the total sensitivities on the model output variables (Table 2) and a bootstrap technique was applied to characterise uncertainty.
The mean relative influence and 95 % confidence limits were plotted for parameters using the LHS (each with 900 parameter sets) for all eight reefs and three release dates (Figs. 3, 4, 5 and 6). The mean variance explained (R2) in the recursive partitioning regression analysis for each response variable and the global mean relative influence of each parameter is presented in Table 3, and plotted as vertical grey bars in Figs. 3, 4, 5 and 6. Of the 10 connectivity response variables, four are presented: local retention (Fig. 3), median geographic distance (Fig. 4), weighted downstream degree centrality (Fig. 5), and metapopulation growth rate (Fig. 6), with the remaining six in the Additional file 1. The importance analysis was completed across all simulations for each vertical behaviour strategy (Behav, Table 1) independently to quantify any behavioural-specific responses. The relative parameter importance across all behaviour strategies were strongly consistent and therefore all results presented do not separate benthic-seeking from diel vertical migration. The relative importance of larval behaviour is then quantified in the swimming speed (Sp), target depth (TD), and homing distance (HmD) parameters.
Across all response variables, high variability resulted from both the geographic context and timing of initiation. Variation driven by the geographic location of the eight reefs is evident as the horizontal displacement when comparing individual bars (reefs) across the vertical extent of each parameter’s row. Similarly, the variation resulting from different release times is seen as the total spread across colours (release dates) along each horizontal plane (unique reefs). Importance values for all simulations are shown, and those that are not significant are plotted as points along the y-axis and some points and lines may fall on top of each other. Generally, across the majority of response variables, larval mortality rate (A) and the length of the pelagic larval stage (PLD) were the most influential parameters (Figs. 3, 4, 5 and 6, Table 3). For the metapopulation growth rate (Fig. 6), reproductive output (RO) outweighed PLD in importance. At the next level of importance, the extent of the precompetency period (PreP) was consistently identified as being an important parameter across all response variables. A suite of parameters displayed highly variable levels of importance across most response variables, including the initial fall velocity (iFV), and swimming speed (Sp). Three parameters were consistently of low importance (HmD, DevP, and TD), with the remaining parameters having minimal influence on population connectivity (ComR, K, B).
Focusing on the local patch scale, retention of individuals within the natal site (Fig. 3, Table 3) increased primarily with a decrease in PLD, a decrease in the precompetency period (PreP), and a decrease in larval mortality (A), with levels below 10 % across all reefs (Additional file 2: Table S1). At local scales, the relative influences of these parameters were somewhat more idiosyncratic when settler diversity (H) and levels of self-recruitment (SR) were considered. Distance-based measures of downstream connectivity consistently ranked PLD as the most influential parameter (distance increased with increasing PLD), followed by the initial fall velocity (iFV), swim speed (Sp), and precompetency period (PreP). The distribution of downstream connectivity (number of connections, dC, and the degree centrality, CW) was determined primarily by larval mortality (A), followed by the precompetency period (PreP), swim speed (Sp), PLD, and to a lesser extent, reproductive output (RO). The proportion of successful settlers was driven by mortality rate (A), followed by the precompetency period (PreP), with swimming speed (Sp), PLD, and homing distance (HmD) all with modest influences. At the metapopulation scale, larval mortality rate (A), the precompetency period (PreP), and the reproductive output (RO) were the most influential drivers, followed by the PLD and swim speed (Sp).
Across the majority of population connectivity response variables the following six parameters consistently had low relative influences: larval mortality shape parameter (B), rate of transition to competency (ComR), initial development time (DevP), bio-physical cohesion (K), the target depth (TD), and homing distance (HmD). The strong influence of geographic location or seascape context was clear when comparing values in the response variables across all reefs (Additional file 2: Table S1).
The unique characteristics of connectivity in marine and aquatic environments, such as the early developmental stages of dispersers and the strongly advective environment, require a more biologically-physically balanced and comprehensive framework to study the importance of individual parameters to this complex process. Here, we have suggested an alternative conceptualisation of connectivity which integrates recent advances across systems and taxa (Fig. 1) and provided insights from applying this framework to the analysis of potential connectivity in Port Phillip Bay, Australia. Although the outcomes presented are from a large semi-enclosed bay, the relative importance of bio-physical parameters would be expected to scale up to broader scales and be relevant for other marine systems. Indeed, the importance of PLD, mortality and behaviour are consistent with many other study systems [40, 38, 74, 36, 75, 19, 76], and the importance of geographic and temporal context is only beginning to be recognised [77, 78].
In particular, the strong temporal and geographic variability evident across all parameters and connectivity metrics (Figs. 3, 4, 5 and 6) may be both a unique consequence of dispersal in the ocean, and relatively unrecognized in the study of population connectivity . The geographic location of individual reefs (Fig. 2) and the surrounding hydrodynamic environment had a large impact on connectivity outcomes (Additional file 2: Table S1), yet strong coherence remained in the relative importance of bio-physical parameters. For example, reef #16 is most isolated (low SR, Additional file 2: Table S1), whereas the nearest neighbours, reefs #18 and #12, are strong sources of dispersing individuals to many patches throughout the Bay (high S, mdG, and mxG). Despite the semi-enclosed nature of the bay, the system is generally well connected for many modelled parameter combinations (mean mxG across all reefs is 0.81, 81 % of total possible downstream connections are made on average). Similarly, the timing of spawning (i.e., initiation) had considerable impact on connectivity outcomes, evident in the wide horizontal spread in importance values across many parameters and connectivity outcomes (Figs. 3, 4 and 5). This strong spatial and temporal variability resulting from the geographic and temporal context suggests predicting demographic consequences or local-scale patterns in marine environments from generalised ‘dispersal syndromes’ [20, 21] may be challenging for benthic species with pelagic larvae where the environment can exert considerable influence on dispersal patterns . Our results suggest that even when species exhibit stereotypical behaviours during early development, local oceanographic conditions can interact with such behaviours to result in fundamentally different dispersal outcomes among locations and release times.
Despite this spatiotemporal variability, broad patterns in the dominant drivers of connectivity outcomes were apparent and strongly consistent. Across the majority of response variables (Table 2), three parameters were consistently the most influential on population connectivity across scales: 1) the mortality rate during the dispersal phase, 2) the maximum duration of the pelagic larval stage (PLD), and 3) the relative duration of the pre-competency and competency windows of dispersers. The importance and main effect of mortality rate (Figs. 3, 4, 5 and 6) on connectivity is intuitive and matches empirical data (e.g., ), theoretical studies (e.g., ), and expectations based on dispersal syndromes (e.g., ). Although the importance of PLD in marine systems has been well recognized [50, 38], the main effect in several connectivity outcomes is less so. Here, an increase in PLD increased the connectivity distance (mdG, mxG) and number of connections (dC), yet decreased total settlement (S), local retention (LR) and self-recruitment (SR), and resulted in negative impacts on metapopulaton growth rate and capacity (λM, λmax, respectively). Increasing the time spent dispersing effectively increases the likelihood of making some long-distance connections, while at the same time, increases the likelihood of being lost at sea resulting in larval wastage. Quantifying these mixed effects of PLD on connectivity outcomes could help disentangle the equally mixed results present in the empirically-based, albeit indirect, correlations between PLD and genetic distance or gene flow .
The relative length of the precompetency period (PreP), a trait unique to marine and aquatic taxa, was the third most important parameter. Increasing the precompetency period decreases the potential settlement window, effectively increasing the geographic distance (mdG and mxG) individuals travel while decreasing the proportion that settled locally (LR, SR). These parameters have been highlighted in the past as being relatively important [32, 40], yet the scale-dependence and effect highlighted here are new. Delaying settlement competency until later in the transport stage causes more individuals to be swept away from their natal habitat patch (lower LR), many not finding suitable habitat (decreases dC, CW) and effectively decreasing metapopulation growth (λM, λmax). Although relatively little is known of the development transition phase of many marine species (but for corals see [52, 81]), it is growth dependent and sensitive to temperature , therefore having obvious implications under future warming scenarios [83, 10].
A suite of parameters were identified as having an intermediate level of influence, often with inconsistencies in the strength and direction of effect on connectivity (right panels in Figs. 3, 4, 5 and 6). These parameters include the initial fall velocity (iFV), swimming capacity (Sp), distance at which individuals can sense suitable reef habitat (HmD), and the reproductive output (RO). The remaining parameters explored had consistently low influence on the connectivity outcomes (B, ComR, K, and TD). Across these low and intermittently important parameters, there was a high level of variability in the direction of impact, making it difficult to identify direct causal relationship.
Considering the proposed four stages of marine population connectivity (Fig. 1), our sensitivity analysis suggests that the settlement stage (PLD, competency, and habitat structure) and the transport and movement stage (mortality and currents) are the most critical drivers of connectivity outcomes in Port Phillip Bay. The initiation of the emigration stage is important in metapopulation-wide measures, which are largely driven by reproductive output (fecundity and abundance of habitat patches). Increased reproductive output also leads to a greater magnitude in connection strengths throughout, although this is not reflected in the relative measures presented here. Recruitment, the final stage of population connectivity, had relatively little influence on connectivity outcomes, compared to the other stages. This was somewhat surprising as post-settlement mortality was allowed to vary between 0 and 100 % on a patch-by-patch basis per simulation. Implicit in our approach to post-settlement mortality is the assumption that the specific cause of increased mortality is acting at, or below, the patch-scale (e.g., poor habitat quality, increased predators) and is not spatially autocorrelated, such as some disturbance events (e.g., storms and urchin outbreaks). Although these added complexities were beyond the scope of the present study, further research is needed to elucidate how environmental heterogeneity and anthropogenic disturbances may influence recruitment and connectivity outcomes.
Implications for future research and management
Marine larval dispersal is biophysically complex and our understanding of this process and its consequences to local and regional population dynamics and management is still limited, despite considerable research. Our findings identify several areas where future research should be targeted.
Better estimates of the key intrinsic drivers of dispersal
Larval mortality was the most important intrinsic parameter influencing dispersal outcomes at all scales. Unfortunately, empirical measures of larval mortality under natural conditions are scarce and most marine dispersal models use guestimates based on the few published studies available  and often assume these rates are invariant . Ecologists and larval biologists have made great strides in recent decades on amassing knowledge of larval durations and the timing of competency to settle for a diversity of marine taxa (e.g., ). Efforts need to be redirected to increase our understanding of mortality, particularly in how rates change as a function of age, size, or condition.
Ground-truthing biophysical models of marine dispersal and population connectivity
Consistent with other research, dispersal outcome are highly variable across space and time. Such spatio-temporal variability makes it difficult to gain a full understanding of the dynamics of marine populations from empirical studies, which are often logistically and financially constrained to snapshots of the dispersal process in space and time. This, in combination with greater computational capability and finer-resolved hydrodynamic models, has led to an increasing focus and reliance on predictions from biophysical models. Although validation of hydrodynamic models is virtually a requirement, corroborating dispersal predictions from such models is considerably more challenging. To our knowledge, no larval dispersal model has been ground-truthed with empirical estimates of dispersal (although this has been achieved in comparisons of generalized ocean circulation with estimates of gene flow among populations (e.g., )). This greatly limits the confidence in modelled estimates of dispersal and connectivity. Our modelling framework provides a mechanism for model parameter tuning and validation in comparison to empirical estimates. Once ground-truthed, simulations can be run across the full spectrum of spatio-temporal variability to generate more realistic estimates of connectivity outcomes and their impacts on metapopulation dynamics.
Applying connectivity models for marine management
Our framework has direct applicability to many marine and aquatic systems, and may assist in gaining a more holistic and integrated view of dispersal-based connectivity to aid in management. Connectivity is a multifaceted process and using a holistic framework to assess the primary drivers of dispersal and population connectivity will lead to greater insight into where best to target management efforts. We suggest that the metapopulation-based metrics, such as λM and λmax, should be used to help identify locations that are important contributors to overall growth rates or capacity to recover from disturbance. Identifying such keystone populations, which are likely to be well-connected populations that are important sources, has been a successful approach for informing conservation efforts in terrestrial ecosystems . We suggest that, with realistic marine population connectivity data, identifying keystone populations will be equally informative in managing natural resources in marine ecosystems.
Here, we have presented clear evidence on the relative importance of the transport and settlement stages to marine population connectivity, the key intrinsic drivers of larval mortality, the length of the pelagic larval phase, and the settlement competency characteristics, and the influence of the extrinsic factors of habitat geography and currents. Gaining a better understanding of these drivers and how they vary across species will greatly enhance our ability to predict contemporary connectivity patterns, will aid in our study of the evolution of larval dispersal strategies (e.g., ), and may guide a proactive approach to understanding species’ potential for adaptation to habitat and climate change to better inform marine environmental management.
Gandon S. Kin competition, the cost of inbreeding and the evolution of dispersal. J Theor Biol. 1999;200(4):345–64. doi:10.1006/jtbi.1999.0994.
Levin SA, Cohen D, Hastings A. Dispersal strategies in patchy environments. Theor Popul Biol. 1984;26(2):165–91.
Travis JMJ, Dytham C. Habitat persistence, habitat availability and the evolution of dispersal. Proc R Soc B-Biol Sci. 1999;266(1420):723–8.
Starrfelt J, Kokko H. The theory of dispersal under multiple influences. In: Clobert J, Baguette M, Benton TG, Bullock JM, editors. Dispersal ecology and evolution. Oxford: Oxford University Press; 2012. p. 19–28.
Hanski I, Ovaskainen O. The metapopulation capacity of a fragmented landscape. Nature. 2000;404(6779):755–8.
Figueira WF, Crowder LB. Defining patch contribution in source-sink metapopulations: the importance of including dispersal and its relevance to marine systems. Popul Ecol. 2006;48(3):215–24. doi:10.1007/s10144-006-0265-0.
Hastings A, Botsford LW. A simple persistence condition for structured populations. Ecol Lett. 2006;9(7):846–52.
Thrush SF, Hewitt JE, Lohrer AM, Chiaroni LD. When small changes matter: the role of cross-scale interactions between habitat and ecological connectivity in recovery. Ecol Appl. 2013;23(1):226–38.
Hanski I. Single-species metapopulation dynamics: concepts, models and observations. Biol J Linnean Soc. 1991;42(1–2):17–38. doi:10.1111/j.1095-8312.1991.tb00549.x.
Munday PL, Jones GP, Pratchett MS, Williams AJ. Climate change and the future for coral reef fishes. Fish Fish. 2008;9(3):261–85. doi:10.1111/j.1467-2979.2008.00281.x.
Moilanen A, Franco AMA, Eary RI, Fox R, Wintle B, Thomas CD. Prioritizing multiple-use landscapes for conservation: methods for large multi-species planning problems. Proc R Soc B-Biol Sci. 2005;272(1575):1885–91.
Gaines SD, Gaylord B, Largier JL. Avoiding current oversights in marine reserve design. Ecol Appl. 2003;13(1):S32–46.
Treml EA, Halpin PN. Marine population connectivity identifies ecological neighbors for conservation planning in the Coral Triangle. Conserv Lett. 2012;5(6):441–9. doi:10.1111/j.1755-263X.2012.00260.x.
Magris RA, Pressey RL, Weeks R, Ban NC. Integrating connectivity and climate change into marine conservation planning. Biol Conserv. 2014;170:207–21. doi:10.1016/j.biocon.2013.12.032.
Beger M, Selkoe KA, Treml E, Barber PH, von der Heyden S, Crandall ED, et al. Evolving coral reef conservation with genetic information. Bull Mar Sci. 2014;90(1):159–85. doi:10.5343/bms.2012.1106.
Bowler DE, Benton TG. Causes and consequences of animal dispersal strategies: relating individual behaviour to spatial dynamics. Biol Rev. 2005;80(2):205–25. doi:10.1017/s1464793104006645.
Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, Saltz D, et al. A movement ecology paradigm for unifying organismal movement research. Proc Natl Acad Sci U S A. 2008;105(49):19052–9. doi:10.1073/pnas.0800375105.
Bonte D, Van Dyck H, Bullock JM, Coulon A, Delgado M, Gibbs M, et al. Costs of dispersal. Biol Rev. 2012;87(2):290–312. doi:10.1111/j.1469-185X.2011.00201.x.
Fossette S, Gleiss AC, Chalumeau J, Bastian T, Armstrong CD, Vandenabeele S, et al. Current-oriented swimming by Jellyfish and its role in bloom maintenance. Curr Biol. 2015;25(3):342–7. doi:10.1016/j.cub.2014.11.050.
Clobert J, Le Galliard JF, Cote J, Meylan S, Massot M. Informed dispersal, heterogeneity in animal dispersal syndromes and the dynamics of spatially structured populations. Ecol Lett. 2009;12(3):197–209. doi:10.1111/j.1461-0248.2008.01267.x.
Ronce O, Clobert J. Dispersal syndromes. In: Clobert J, Baguette M, Benton TG, Bullock JM, editors. Dispersal ecology and evolution. Oxford: Oxford University Press; 2012. p. 119–38.
Stevens VM, Whitmee S, Le Galliard JF, Clobert J, Bohning-Gaese K, Bonte D, et al. A comparative analysis of dispersal syndromes in terrestrial and semi-terrestrial animals. Ecol Lett. 2014;17(8):1039–52. doi:10.1111/ele.12303.
Ricketts TH. The matrix matters: effective isolation in fragmented landscapes. Am Nat. 2001;158(1):87–99. doi:10.1086/320863.
Jones GP, Almany GR, Russ GR, Sale PF, Steneck RS, van Oppen MJH, et al. Larval retention and connectivity among populations of corals and reef fishes: history, advances and challenges. Workshop on Connectivity and Resilience Sustaining Coral Reefs during the coming Century. Townsville, AUSTRALIA: Springer; 2009.
Kool JT, Moilanen A, Treml EA. Population connectivity: recent advances and new perspectives. Landsc Ecol. 2013;28(2):165–85. doi:10.1007/s10980-012-9819-z.
Cowen RK, Paris CB, Olson D, Fortuna JL. The role of long distance dispersal versus local retention in replenishing marine populations. Gulf and Caribbean Research. 2003;14(2):129–37.
Swearer SE, Caselle JE, Lea DW, Warner RR. Larval retention and recruitment in an island population of a coral-reef fish. Nature. 1999;402(6763):799–802. doi:10.1038/45533.
Cowen RK, Paris CB, Srinivasan A. Scaling of connectivity in marine populations. Science. 2006;311(5760):522–7.
Hastings A, Botsford LW. Persistence of spatial populations depends on returning home. Proc Natl Acad Sci U S A. 2006;103(15):6067–72.
Andutta FP, Kingsford MJ, Wolanski E. ‘Sticky water’ enables the retention of larvae in a reef mosaic. Estuar Coast Shelf Sci. 2012;101:54–63. doi:10.1016/j.ecss.2012.02.013.
Paris CB, Cherubin LM, Cowen RK. Surfing, spinning, or diving from reef to reef: effects on population connectivity. Mar Ecol Prog Ser. 2007;347:285–300. doi:10.3354/meps06985.
Treml EA, Roberts JJ, Chao Y, Halpin PN, Possingham HP, Riginos C. Reproductive output and duration of the pelagic larval stage determine seascape-wide connectivity of marine populations. Integr Comp Biol. 2012;52(4):525–37. doi:10.1093/icb/ics101.
Jones GP, Planes S, Thorrold SR. Coral reef fish larvae settle close to home. Curr Biol. 2005;15(14):1314–8. doi:10.1016/j.cub.2005.06.061.
Saenz-Agudelo P, Jones GP, Thorrold SR, Planes S. Connectivity dominates larval replenishment in a coastal reef fish metapopulation. Proc R Soc B-Biol Sci. 2011;278(1720):2954–61. doi:10.1098/rspb.2010.2780.
Saenz-Agudelo S, Jones GP, Thorrold SR, Planes S. Estimating connectivity in marine populations: an empirical evaluation of assignment tests and parentage analysis under different gene flow scenarios. Mol Ecol. 2009;18(8):1765–76.
Cowen RK, Lwiza KMM, Sponaugle S, Paris CB, Olson DB. Connectivity of marine populations: open or closed? Science. 2000;287(5454):857–9.
Gerlach G, Atema J, Kingsford MJ, Black KP, Miller-Sims V. Smelling home can prevent dispersal of reef fish larvae. Proc Natl Acad Sci U S A. 2007;104(3):858–63.
Shanks AL. Pelagic larval duration and dispersal distance revisited. Biol Bull. 2009;216(3):373–85.
Largier JL. Considerations in estimating larval dispersal distances from oceanographic data. Ecol Appl. 2003;13(1):S71–89.
Wolanski E, Kingsford MJ. Oceanographic and behavioural assumptions in models of the fate of coral and coral reef fish larvae. J R Soc Interface. 2014;11(98). doi:10.1098/rsif.2014.0209.
Putman NF, He RY. Tracking the long-distance dispersal of marine organisms: sensitivity to ocean model resolution. J R Soc Interface. 2013;10(81). doi:10.1098/rsif.2012.0979.
Levin LA. Recent progress in understanding larval dispersal: new directions and digressions. Integr Comp Biol. 2006;46(3):282–97. doi:10.1093/icb/icj024.
Cowen RK, Sponaugle S. Larval dispersal and marine population connectivity. Ann Rev Mar Sci. 2009;1(1):443–66.
Lowe WH, Allendorf FW. What can genetics tell us about population connectivity? Mol Ecol. 2010;19(15):3038–51. doi:10.1111/j.1365-294X.2010.04688.x.
Buoro M, Carlson SM. Life-history syndromes: integrating dispersal through space and time. Ecol Lett. 2014;17(6):756–67. doi:10.1111/ele.12275.
Nathan R, Schurr FM, Spiegel O, Steinitz O, Trakhtenbrot A, Tsoar A. Mechanisms of long-distance seed dispersal. Trends Ecol Evol. 2008;23(11):638–47. doi:10.1016/j.tree.2008.08.003.
Ronce O. How does it feel to be like a rolling stone? Ten questions about dispersal evolution. Annual Review of Ecology Evolution and Systematics. Annual Review of Ecology Evolution and Systematics. Palo Alto: Annual Reviews; 2007. p. 231–53.
Shima JS, Swearer SE. The legacy of dispersal: larval experience shapes persistence later in the life of a reef fish. J Anim Ecol. 2010;79(6):1308–14. doi:10.1111/j.1365-2656.2010.01733.x.
Travis JMJ, Mustin K, Barton KA, Benton TG, Clobert J, Delgado MM, et al. Modelling dispersal: an eco-evolutionary framework incorporating emigration, movement, settlement behaviour and the multiple costs involved. Methods Ecol Evol. 2012;3(4):628–41. doi:10.1111/j.2041-210X.2012.00193.x.
Kinlan BP, Gaines SD. Propagule dispersal in marine and terrestrial environments: a community perspective. Ecology. 2003;84(8):2007–20.
Sponaugle S, Cowen RK, Shanks A, Morgan SG, Leis JM, Pineda JS, et al. Predicting self-recruitment in marine populations: biophysical correlates and mechanisms. Bull Mar Sci. 2002;70(1):341–75.
Connolly SR, Baird AH. Estimating dispersal potential for marine larvae: dynamic models applied to scleractinian corals. Ecology. 2010;91(12):3572–83. doi:10.1890/10-0143.1.
Marshall DJ, Monro K, Bode M, Keough MJ, Swearer S. Phenotype-environment mismatches reduce connectivity in the sea. Ecol Lett. 2010;13(1):128–40.
Burgess SC, Treml EA, Marshall DJ. How do dispersal costs and habitat selection influence realized population connectivity? Ecology. 2012;93(6):1378–87.
Pechenik JA, Wendt DE, Jarrett JN. Metamorphosis is not a new beginning. Bioscience. 1998;48(11):901–10. doi:10.2307/1313294.
Edgar GJ, Stuart-Smith RD. Ecological effects of marine protected areas on rocky reef communities-a continental-scale analysis. Mar Ecol Prog Ser. 2009;388:51–62. doi:10.3354/meps08149.
Jenkins GP, Wheatley MJ. The influence of habitat structure on nearshore fish assemblages in a southern Australian embayment: Comparison of shallow seagrass, reef-algal and unvegetated sand habitats, with emphasis on their importance to recruitment. J Exp Mar Biol Ecol. 1998;221(2):147–72.
Black K, Hatton D, Rosenberg M. Locally and externally-driven dynamics of a large semi-enclosed bay in Southern Australia. J Coast Res. 1993;9(2):509–38. doi:10.2307/4298106.
Lee RS, Black KP, Bosserel C, Greer D. Present and future prolonged drought impacts on a large temperate embayment: Port Phillip Bay, Australia. Ocean Dyn. 2012;62(6):907–22. doi:10.1007/s10236-012-0538-4.
Black K. The hydrodynamic model 3DD and support software. New Zealand: Department of Earth Sciences, University of Waikato. Contract No.: 19. 1995.
Black KP, Moran PJ, Hammond LS. Numerical-models show coral reefs can be self-seeding. Mar Ecol Prog Ser. 1991;74(1):1–11.
Blake S, Ball D. Victorian marine habitat database: seagrass mapping of Port Phillip Bay. In: The State of Victoria DoNRaE, editor. Marine and freshwater resources institute. Queenscliff: Geospatial Systems Section, Marine and Freshwater Resources Institute; 2001.
Black K. Lagrangian dispersal and sediment transport model POL3DD. New Zealand: Department of Earth Sciences, University of Waikato. Contract No.: 21. 1996.
Jenkins GP, Black KP, Hamer PA. Determination of spawning areas and larval advection pathways for King George whiting in southeastern Australia using otolith microstructure and hydrodynamic modelling. I. Victoria. Mar Ecol Prog Ser. 2000;199:231–42. doi:10.3354/meps199231.
Jenkins GP, Black KP, Wheatley MJ, Hatton DN. Temporal and spatial variability in recruitment of a temperate, seagrass-associated fish is largely determined by physical processes in the pre- and post-settlement phases. Mar Ecol Prog Ser. 1997;148(1–3):23–35. doi:10.3354/meps148023.
Leis JM. Are larvae of demersal fishes plankton or nekton? In: Southward AJSDW, editor. Advances in Marine Biology, Vol 51. Advances in Marine Biology, 2006. p. 57–141. http://www.elsevier.com/books/advances-in-marine-biology/sims/978-0-12-026152-9
Doherty PJ, Dufour V, Galzin R, Hixon MA, Meekan MG, Planes S. High mortality during settlement is a population bottleneck for a tropical surgeonfish. Ecology. 2004;85(9):2422–8. doi:10.1890/04-0366.
Opsahl T, Agneessens F, Skvoretz J. Node centrality in weighted networks: generalizing degree and shortest paths. Soc Networks. 2010;32(3):245–51.
Makler-Pick V, Gal G, Gorfine M, Hipsey MR, Carmel Y. Sensitivity analysis for complex ecological models – a new approach. Environ Model Software. 2011;26(2):124–34.
Saltelli A, Scott EM, Chan K. Sensitivity analysis. In: Andrea S, Karen C, Marian Scott E, editors. Wiley series in probability and statistics. Chichester: Wiley, c2000; 2000.
De’ath G. Boosted trees for ecological modeling and prediction. Ecology. 2007;88(1):243–51. doi:10.1890/0012-9658(2007)88[243:BTFEMA]2.0.CO;2.
Coutts S, Yokomizo H. Meta-models as a straightforward approach to the sensitivity analysis of complex models. Popul Ecol. 2014;56(1):7–19. doi:10.1007/s10144-013-0422-1.
Stocki R. A method to improve design reliability using optimal Latin hypercube sampling. Comput Assist Mech Eng Sci. 2005;12(4):393–411.
Cowen RK, Gawarkiewicz G, Pineda J, Thorrold SR, Werner FE. Population connectivity in marine systems: an overview oceanography. Special Issue on Marine Population Connectivity. 2007;20(3):14–21.
Putman NF, Scott R, Verley P, Marsh R, Hays GC. Natal site and offshore swimming influence fitness and long-distance ocean transport in young sea turtles. Mar Biol. 2012;159(10):2117–26. doi:10.1007/s00227-012-1995-5.
Scott R, Marsh R, Hays GC. A little movement orientated to the geomagnetic field makes a big difference in strong flows. Mar Biol. 2012;159(3):481–8. doi:10.1007/s00227-011-1825-1.
Carson HS, López-Duarte PC, Rasmussen L, Wang D, Levin LA. Reproductive timing alters population connectivity in marine metapopulations. Current biology: CB. 2010;20(21):1926–31.
Treml EA, Roberts J, Halpin PN, Possingham HP, Riginos C. The emergent geography of biophysical dispersal barriers across the Indo-West Pacific. Divers Distrib. 2015;21(4):465–76. doi:10.1111/ddi.12307.
Siegel DA, Mitarai S, Costello CJ, Gaines SD, Kendall BE, Warner RR et al. The stochastic nature of larval connectivity among nearshore marine populations. Proc Natl Acad Sci. 2008:0802544105. doi:10.1073/pnas.0802544105.
Weersing K, Toonen RJ. Population genetics, larval dispersal, and connectivity in marine systems. Mar Ecol Prog Ser. 2009;393:1–12. doi:10.3354/meps08287.
Figueiredo J, Baird AH, Connolly SR. Synthesizing larval competence dynamics and reef-scale retention reveals a high potential for self-recruitment in corals. Ecology. 2012. doi:10.1890/12-0767.1.
O’Connor MI, Bruno JF, Gaines SD, Halpern BS, Lester SE, Kinlan BP, et al. Temperature control of larval dispersal and the implications for marine ecology, evolution, and conservation 10.1073/pnas.0603422104. Proc Natl Acad Sci U S A. 2007;104(4):1266–71.
Gerber LR, Mancha-Cisneros MDM, O’Connor MI, Selig ER. Climate change impacts on connectivity in the ocean: Implications for conservation. Ecosphere. 2014;5(3):art33. doi:10.1890/ES13-00336.1.
Houde ED. Comparative growth, mortality, and energetics of marine fish larvae: temperature and implied latitudinal effects. Fish Bull. 1989;87(3):471–95.
Crandall ED, Treml EA, Barber PH. Coalescent and biophysical models of stepping-stone gene flow in neritid snails. Mol Ecol. 2012;21(22):5579–98. doi:10.1111/mec.12031.
Pulliam HR, Danielson BJ. Sources, sinks, and habitat selection - a landscape perspective on population-dynamics. Am Nat. 1991;137:S50–66.
We thank Randall Lee and Dougal Greer for hydrodynamic modelling and data assistance, David Ball for PPB reef habitat data, and the Swearer Lab, Graeme Hays and two anonymous reviewers for constructive comments on earlier drafts. Funding for this work was provided by the Victoria State Government, Department of Environment and Primary Industries.
The authors declare that they have no competing interests.
ET, JF and SS designed the research and approach and KB wrote the dispersal model. ET developed the spatial model and completed the analysis. ET, JF, and SS interpreted the data, ET and SS wrote the manuscript. All authors read and approved the final manuscript.
Sensitivity results for select output parameters are shown as two-panel plots of the influence of model parameters (y-axis) on remaining model output.
Median values for all reef-based response variables are reported per reef patch and summarized across all reefs in system.