Estimating dispersal distance in the deep sea: challenges and applications to marine reserves

Hilário, Ana; Metaxas, Anna; Gaudron, Sylvie M.; Howell, Kerry L.; Mercier, Annie; Mestre, Nélia C.; Ross, Rebecca E.; Thurnherr, Andreas M.; Young, Craig

doi:10.3389/fmars.2015.00006

REVIEW article

Front. Mar. Sci., 13 February 2015

Sec. Deep-Sea Environments and Ecology

Volume 2 - 2015 | https://doi.org/10.3389/fmars.2015.00006

Estimating dispersal distance in the deep sea: challenges and applications to marine reserves

$\r\nAna Hilrio$ Ana Hilário¹

Anna Metaxas²^*

Andreas M. Thurnherr⁷

Craig Young⁸

¹Departamento de Biologia & Centre for Environmental and Marine Studies (CESAM), Universidade de Aveiro, Aveiro, Portugal
²Department of Oceanography, Dalhousie University, Halifax, NS, Canada
³UMR8187 Laboratoire d'Océanologie et de Géosciences, Station Marine de Wimereux, Wimereux, France
⁴Marine Biology and Ecology Research Centre, Marine Institute, Plymouth University, Plymouth, UK
⁵Department of Ocean Sciences, Memorial University, St. John's, NL, Canada
⁶Centro de Investigação Marinha e Ambiental, Universidade do Algarve, Faro, Portugal
⁷Lamont-Doherty Earth Observatory, The Earth Institute of Columbia University, Palisades, NY, USA
⁸Oregon Institute of Marine Biology, University of Oregon, Charleston, OR, USA

Population connectivity refers to the exchange of individuals among populations: it affects gene flow, regulates population size and function, and mitigates recovery from natural or anthropogenic disturbances. Many populations in the deep sea are spatially fragmented, and will become more so with increasing resource exploitation. Understanding population connectivity is critical for spatial management. For most benthic species, connectivity is achieved by the planktonic larval stage, and larval dispersal is, in turn, regulated by complex interactions between biological and oceanographic processes. Coupled biophysical models, incorporating ocean circulation and biological traits, such as planktonic larval duration (PLD), have been used to estimate population connectivity and generate spatial management plans in coastal and shallow waters. In the deep sea, knowledge gaps in both the physical and biological components are delaying the effective use of this approach. Here, we review the current efforts in conservation in the deep sea and evaluate (1) the relevance of using larval dispersal in the design of marine protected areas and (2) the application of biophysical models in the study of population connectivity. Within biophysical models, PLD can be used to estimate dispersal distance. We propose that a PLD that guarantees a minimum dispersal distance for a wide range of species should be used in the planning of marine protected areas in the deep sea. Based on a review of data on species found at depths >200 m, a PLD of 35 and 69 days ensures a minimum distance for 50 and 75%, respectively, of eurybathic and deep-sea species. We note that more data are required to enhance accuracy and address the high variability in PLD between and within taxonomic groups, limiting generalizations that are often appealing to decision-makers. Given the imminent expansion of resource exploitation in the deep sea, data relevant to spatial management are needed urgently.

Introduction

The deep sea, although the largest biome on the planet, remained unexplored until the late 19th century (Tyler, 2003) when the cosmopolitan occurrence of deep-sea fauna was established for the first time. Another 100 years of research and technological developments were necessary before the habitat heterogeneity (Ramirez-Llodra et al., 2010), high biodiversity (e.g., Grassle and Maciolek, 1992), and contribution to global ocean processes (Jahnke, 1996) of the deep sea challenged prevailing views and raised new ecological questions (Danovaro et al., 2014). In the last few decades, demand for deep-sea products, such as those from fishing, hydrocarbon extraction, and mining, has been rapidly expanding (Ramirez-Llodra et al., 2011; Thurber et al., 2014), and deep-sea ecologists are asked to provide solutions for the mitigation of exploitation impacts.

Although the deep seafloor includes some of the largest contiguous features of the planet, such as the abyssal plains and the sedimentary slopes of the continental margins (Ramirez-Llodra et al., 2010), many deep-sea populations are spatially fragmented, and may become more so as a consequence of human disturbance during resource exploitation and extraction. Concurrently, one of the main challenges of deep-sea ecology is the elucidation of the processes that lead to connectivity among spatially isolated populations, which would ultimately regulate their persistence and recovery after disturbance (reviewed in Cowen et al., 2007). Differences in population connectivity contribute greatly to the spatiotemporal patterns in the distribution of organisms and must be considered when defining spatial management strategies, including in the design of marine protected areas (MPAs) (Gaines et al., 2003).

Many marine benthic species have complex life cycles that include a pelagic larval stage and sessile/sedentary adults (Thorson, 1964). For these species, the main process that connects populations is larval transport; consequently, the factors that regulate larval dispersal and population connectivity have been receiving increased attention. Although Thorson (1950) predicted that deep-sea animals should brood their young or have direct development, recent analyses show that pelagic lecithotrophy, rather than brooding, is the main reproductive mode in the deep sea (reviewed in Young, 2003). Different approaches have been used to evaluate population connectivity by larval dispersal including: (i) measuring the hydrodynamic and biological processes involved in larval transport (e.g., Jackson et al., 2010; Thurnherr et al., 2011; Mullineaux et al., 2013); and (ii) deriving larval origins and dispersal pathways using genetic or geochemical markers (Levin, 2006; Cowen and Sponaugle, 2009), or a combination of the two (Foster et al., 2012). To achieve a mechanistic understanding of larval transport, the interaction of biological and physical processes on different spatial and temporal scales needs to be elucidated (Cowen et al., 2007). Because sampling over all relevant scales is presently not possible, numerical models that incorporate both physical dynamics and biological traits are increasingly being used to quantify larval transport and assess its role in regulating population connectivity (e.g., Cowen et al., 2006; Siegel et al., 2008; Nolasco et al., 2013). In coastal and shallow areas, such coupled biophysical models have provided information of relevance to decision-makers in determining the spatial arrangement of marine reserves (e.g., Guizien et al., 2012; Treml et al., 2012). However, in the deep sea, this field of research is still in its infancy and fewer than a handful of studies have explored this approach (Lavelle et al., 2010; McGillicuddy et al., 2010; Yearsley and Sigwart, 2011; Young et al., 2012).

Biological parameterization of the biophysical models presents a challenge even in well studied shallow-water systems (Metaxas and Saunders, 2009). Components, such as diel and ontogenetic vertical migration, buoyancy of embryos, mortality, food availability, developmental rate, and physiological tolerances, can play an important role in dispersal patterns, and should be incorporated in biophysical models. However, most of these factors have not been studied for many species, particularly in the deep sea. In contrast, the pelagic or planktonic larval duration (PLD), commonly defined as the developmental period of a species in the water column, has been estimated for a relatively large number of marine fishes and invertebrates (Shanks et al., 2003; Shanks, 2009) and is perhaps the most often cited biological variable potentially affecting population connectivity (Sponaugle et al., 2002). Although the validity of the long-standing hypothesis that species with long larval duration also have greater dispersal potential remains equivocal (Eckert, 2003; Siegel et al., 2003; Weersing and Toonen, 2009), PLD can be used to set an upper bound on dispersal distance (Selkoe and Toonen, 2011).

In this study, we firstly provide an overview of the current efforts in conservation associated with resource extraction in deep-sea seafloor habitats and evaluate the relevance of larval dispersal in the design of marine reserves. Secondly, we review the application of biophysical models to the study of population connectivity in the deep sea, and provide an evaluation of their performance. Lastly, we assess the extent of the current knowledge on PLD for deep-sea species as one of the main biological components consistently included in biophysical models. We then compare the available estimates from the deep sea with those from the closest taxonomic relatives that live in shallow systems to determine whether PLD is taxonomically conserved. Genetic tools have also been used in many studies to estimate dispersal distances and genetic connectivity, which we do not review here. A separate synthesis is underway, focused specifically on genetic-based estimates of deep-sea dispersal distances, including analyses of how different life-history factors may affect these estimates and that compare these estimates with those in shallow water (A. Baco et al in prep). These two syntheses are parallel products of the Population Connectivity working group of INDEEP (International Network of Scientific Investigations of Deep-Sea Ecosystems; www.indeep-project.org/). The ultimate goal is to provide recommendations for obtaining accurate estimates of larval dispersal and population connectivity that can be used on the spatial management of different deep-sea habitats.

Using Larval Dispersal in the Design of Marine Reserves in the Deep Sea

With the depletion of mineral and biological resources on land and in coastal waters, resource extraction has been extending into the water column and the seafloor of the deep sea. Oil and gas have long been extracted offshore, in waters >200 m in depth. The Deepwater Horizon blowout in the Gulf of Mexico in 2010 was the largest oil spill in history (in terms of amount of oil spilled) and one that occurred directly on the seafloor of the deep sea. It impacted deep-water coral communities as far as 22 km away from the accident site at depths of 1850–1950 m (Fisher et al., 2014). The presence of some of these communities was unknown until surveys were conducted after the oil spill. The continuous expansion of oil and gas exploration onto the continental margins all around the globe is greatly enhancing the threat for similar accidents and impacts. Deep-water fishing has been occurring since the late 1950s but developed into a commercial industry in the last 40 years. The impact of bottom trawling on both deep-sea fish and benthic communities has been highlighted by a number of studies (Koslow et al., 2000; Bailey et al., 2009), and growing concern has resulted in recent proposals for a ban on deep-sea bottom trawling in European waters. An emerging potential pressure on the seafloor is through the development of deep-ocean industrial mining, which is rapidly gaining momentum. Deep-sea mining will potentially encompass polymetallic nodules in the abyssal plains, deposits of seafloor massive polymetallic sulfides (SMS) from hydrothermal vents, cobalt crusts from seamounts, among others. As of mid-2014, the International Seabed Authority (ISA) had granted 16 exploration contracts in the Atlantic, Pacific and Indian Oceans, covering all three types of resources (ISBA, 2014). However, deep-sea mining will occur both in national and international jurisdictions, and the laws and regulations that will apply to the industry are currently still under development.

Although the extraction methods can potentially be highly destructive, the spatial and temporal scales of their impact are not known. Additionally, both the biological communities and the drivers that regulate these communities are mostly unknown, particularly in the abyssal plains, but even at some of the hydrothermal vents and seamounts currently being targeted for exploration. A high probability of endemicity at some of these locations further underscores the potential impact of unfettered anthropogenic activities on these largely undescribed ecosystems. For these reasons, a recent call for the precautionary approach in the management of human activities in the deep sea includes plans for protection of the ecosystems, research to increase knowledge, and governance collaboration across sectors (Mengerink et al., 2014).

The strategies being considered for spatial management and protection in the deep sea are based on our existing practices from shallow waters: MPAs or other marine reserves partially or fully restricting certain human activities (Halpern, 2003). The spatial arrangements of marine reserves can vary from a single reserve to a network of many reserves within a habitat or region, but the target is usually a subset of habitats (or species) in a region (Hastings and Botsford, 2003; Sale et al., 2005). The selection of the subset of habitats to be protected is based on how representative they are within the region of interest, their uniqueness or rarity, their vulnerability to potential threats or presumed slow recovery from disturbance, whether they support high biodiversity or high productivity, or whether they are a key habitat for a particular stage in the life-history of species, particularly if the latter is threatened or endangered (CBD, 2007). Many of these selection criteria apply to habitats in the deep sea, particularly their vulnerability and presumed slow recovery from perturbation and high biodiversity (e.g., the abyssal plains), uniqueness (e.g., hydrothermal vents) and importance to certain life-history stages (e.g., seamounts). For a network of marine reserves or MPAs, additional criteria apply, such as maximizing connectivity between individual MPAs and maintaining viable populations across the network. These last two criteria are closely linked, particularly in spatially fragmented populations where persistence of a population will depend on either sufficiently large local replenishment in a single patch or, in its absence, sufficiently strong connectivity among patches in a network (Burgess et al., 2014). Since larval dispersal influences population connectivity (Cowen and Sponaugle, 2009), knowledge of the magnitude and pathways of dispersal can be critical elements in the design of effective marine reserves.

The concept of marine reserves in the deep sea in the face of potentially heightened exploitation is increasingly gaining support. An example from areas under international jurisdiction include an environmental management plan for the Clarion-Clipperton Zone (CCZ), generated by the International Seabed Authority (ISA) in which the Authority acknowledges its responsibility to afford effective protection of the environment from harmful effects of prospecting, exploitation and exploration activities (ISBA, 2011). The conservation objectives were to maintain regional biodiversity, ecosystem structure and ecosystem function across the CCZ, manage the CCZ consistently with the principles of integrated ecosystem-based management, and enable the preservation of representative and unique marine ecosystems. Based on environmental and ecological data, which included (presumed) faunal dispersal capabilities and distances, the ISA recommended the allocation of 9 areas of environmental interest, within each of 9 biogeographic subregions, each 400 × 400 km (a 200 × 200 km core area surrounded by a 100-km buffer zone) (ISBA, 2011; Wedding et al., 2013). The size of each area was presumed to be sufficient for each to maintain viable populations of species potentially restricted to a sub-region. Assessing the viability of a population requires the combined estimates of larval retention, reproductive output of the population, and population size at minimum (Burgess et al., 2014); this information largely did not exist when ISA made its recommendations. Further examples of marine reserves in the deep sea include the OSPAR network of MPAs in the North East Atlantic (OSPAR, 2003) as well as bottom trawl closures on the Mid Atlantic Ridge and various seamounts for the protection of vulnerable marine ecosystems (NEAFC, 2011). Within national jurisdiction, a few MPAs (or a national monument in the case of the US) have been established for the protection, at least in part, of deep-sea hydrothermal vents and they include the Endeavor Segment (Canada), the Marianas Trench region (USA), the mid-Atlantic Ridge off the Azores (Portugal), and the Guaymas Basin and Eastern Pacific Rise (Mexico). In addition, particularly within Europe and the USA, a number of MPAs have been established for the protection of vulnerable deep-sea habitats, principally cold water corals and deep-sea sponges. Generally, the conservation objectives and management plans of these align with those agreed upon in the Convention of Biological Diversity. As in the example with the CCZ, for all of the national MPAs, the information required to assess population connectivity and viability was weak to non-existent.

Recommendations on the spatial management through marine reserves have also been made by the scientific community directly. For example, the Dinard workshop, attended by stakeholders from various sectors and 14 countries provided a clear set of guidelines for setting up reserves in chemosynthetic environments (Ardron et al., 2011; Van Dover et al., 2012). The proposed design followed the same criteria as recommended for shallow water MPAs, including “ensuring connectivity” among reserves (Van Dover et al., 2012). Clark et al. (2014) assessed how each of the criteria for the selection of Ecologically and Biologically Significant Areas (EBSAs), described above in the context of MPAs, can be applied to deep-sea ecosystems, and provided a test case on seamounts in the Pacific Ocean. Other studies have addressed the management requirements and have made recommendations for conservation specifically with respect to mining SMS deposits, cobalt-rich crust regions on seamounts and manganese nodules in abyssal plains (Boschen et al., 2013; Wedding et al., 2013; Schlacher et al., 2014). Boschen et al. (2013) recommended the establishment of “preservation reference zones” during SMS mining, including upstream set-asides that can supply colonizing larvae, in addition to preserving an unimpacted section of the population. Wedding et al. (2013) provided a systematic framework for conservation in abyssal plains, including the incorporation of design principles utilized in shallow water, such as ecosystem-based management and networks of MPAs. The concepts of realized dispersal distances and the size of each MPA in the network, particularly as they relate to population viability, were addressed in the framework (Wedding et al., 2013).

While spatial planning is gaining attention in the context of deep-sea resource extraction, the data to support decisions are scarce. Classification systems of deep-sea habitat are being developed on which the criteria of representation of habitats within MPA networks may then be assessed (Howell, 2010; Howell et al., 2010), and species-area relationships have been used to inform baseline conservation targets for the deep North East Atlantic (Foster et al., 2013). Although recent science-based studies have started to address their relevance to the conservation of potentially vulnerable ecosystems in the deep sea (e.g., Rengstorf et al., 2013; Ross and Howell, 2013; Jackson et al., 2014; Nakajima et al., 2014), studies designed specifically to collect relevant data are still lacking. In addition to the attention being recent, the logistical constraints in collecting data from a remote environment, such as the deep sea, are great. Data on larval dispersal and population connectivity that are purported to be relevant in the design of marine reserves are particularly difficult to obtain, even in the more accessible shallow-water ecosystems (Burgess et al., 2014). Numerical models are one promising approach allowing the calculation of dispersal matrices under different scenarios, and their performance can be progressively improved with gaining biological and physical information.

Applications of Biophysical Models to Larval Dispersal

Lagrangian particle tracking methods, traditionally employed by atmospheric scientists and oceanographers, can be used to simulate the release of passive particles to track the fate of the advected drifters in the ocean. These passive drifters can be used to represent theoretical larvae in order to identify likely dispersal pathways, highlighting the oceanographic mechanisms and barriers to dispersal (Werner et al., 2007). Particles can also be given “behavior,” simulating swimming abilities such as diel vertical migration or ontogenetic buoyancy properties, to adjust predictions where dispersal is likely not passive (Levin, 2006; Werner et al., 2007).

Many of the biological parameters (e.g., planktonic larval duration, larval buoyancy, mortality over time, vertical migration, settlement probability, settlement behavior) included in the biophysical models cannot be estimated at this time for marine benthic species (in shallow or deep waters) (Metaxas and Saunders, 2009). Where these data are available, they may alter estimates of connectivity (Cowen et al., 2000). Many tracer parameters require data derived from biological traits, such as number of particles released (fecundity), particle release depth (spawning location), and particle tracking time (planktonic larval duration). The accuracy in these tracer characteristics along with rate of particle loss from the system (mortality) and “behavior” is the challenge that biologists face when they attempt to predict and validate dispersal pathways (Metaxas and Saunders, 2009).

Biophysical models can be applied over a variety of time scales (Levin, 2006) and studies based on contemporary time scales can provide insight into current metapopulation management and demographic dynamics (Cowen and Sponaugle, 2009). Modeling can be used to define the average route and distance of dispersal paths from a release site (Cowen et al., 2007; Cowen and Sponaugle, 2009; Kool et al., 2013) and provide retention estimates and, thus, inform MPA networks (Paris and Cowen, 2004; Treml and Halpin, 2012). Most modeling studies are used to explore the effect of physical drivers on dispersal (Martins et al., 2010; Blanke et al., 2012; Soria et al., 2012; Young et al., 2012) but, in a system where knowledge of life-history traits is lacking, null models can test the effects of behavior on dispersal (Paris et al., 2007, 2009; Carr et al., 2008; Sundelöf and Jonsson, 2012), the relationship between PLD and dispersal ability (Siegel et al., 2003; Young et al., 2012), and the role of other abiotic/biotic factors in dispersal (Ayata et al., 2010; Martins et al., 2010; McGillicuddy et al., 2010; Treml and Halpin, 2012). The ability to run a model both forwards and backwards in time also allows for the prediction of both sources and sinks of propagules (Brickman et al., 2009). Results from these modeling efforts can then be compared to empirical data for biological validation (e.g., Foster et al., 2012), while also being used to constrain dispersal estimates within an ecological timeframe.

Because of the paucity of data on life histories of deep-sea fauna as well as deep-sea circulation, accurate modeling and precise validation are not feasible at this point. Basic models, which exclude life-history data can provide estimates of the bounds of dispersal ranges for future validation and hypothesis generation. Studies on population connectivity in the deep sea may benefit from biophysical modeling more than in shallow environments, given the inherent barriers of accessibility, scale and expense associated with the collection of samples. However, sampling is still required for validation of models and predictions.

Choosing Physical Models for Dispersal Studies in the Deep Sea

Horizontal dispersal of planktonic propagules, such as larvae, is primarily passive, i.e., the greatest component of displacement is through advection by the oceanic velocity field. If the velocities are known across all scales of interest, dispersal reduces to a problem of advection. If, as is usually the case, knowledge of the velocity field is incomplete, the effects of the unknown velocities must be parameterized somehow. Often, the unresolved velocities are modeled as random-walk processes (e.g., Berg, 1993), which cause down-gradient fluxes proportional to the spatial gradients (Fickian diffusion), leading to advection-diffusion models.

There is a hierarchy of techniques that has been used to infer bounds on planktonic dispersal in the deep sea using advection-diffusion models, with the simplest ones being either purely advective or diffusive. The simplest advective model uses a representative “mean” velocity together with PLDs to estimate a dispersal distance (e.g., McClain and Hardy, 2010). As this model ignores both spatial and temporal variability in the velocity field, the relevance of the resulting estimates is restricted to temporal and spatial scales over which the circulation can be considered steady and homogeneous. On time scales that typically range from weeks to months and even years, horizontal dispersal across most of the deep ocean is either dominated or strongly affected by eddy diffusion (Speer et al., 2003), implying that advection-by-mean-flow models are not appropriate. Another simple technique that has sometimes been used to estimate dispersal distances is based on progressive vector diagram (PVDs) derived from Eulerian measurements (e.g., Marsh et al., 2001). While temporal variability of the velocity field is included in this method, spatial variability is not, restricting the relevance of the resulting estimates to the scales of the processes that dominate the velocities, such as eddies, Rossby waves, equatorial jets, and boundary currents. In particular near sloping topography (continental slopes, seamounts, mid-ocean ridges, etc.), the spatial scales of subinertial oceanic flows are often on the order of kilometers (e.g., Brink, 1995; Cannon and Pashinski, 1997; Stahr and Sanford, 1999; Thurnherr and Richards, 2001; McGillicuddy et al., 2010; Thurnherr et al., 2011), limiting the use of PVDs to temporal and spatial scales of days and 10 s of kilometers at most.

Given the difficulty of observing the velocity field in the deep ocean across a wide range of spatial and temporal scales, numerical ocean general circulation models (OGCMs) are often the only viable option for obtaining the velocities that are required to study larval dispersal. However, it is unlikely that all physical processes will be captured in a single OGCM: their greatest inadequacies are in relation to the processes occurring at the fine spatial and temporal scales which are of greatest relevance to a larva (Metaxas and Saunders, 2009). For example, Lacroix et al. (2009) suggested that a 3-km grid is required to resolve a 20-km eddy indicating the intensity of data coverage required to resolve sub-mesoscale and small scale processes. Because of the small scales involved, modeling the circulations near topography often requires dedicated regional models with high spatial and temporal resolution (e.g., Proehl et al., 2005; Mitarai et al., 2009; Lavelle et al., 2010; McGillicuddy et al., 2010), which can be hard to source. As even the highest-resolution regional models cannot resolve the small scales associated with mechanical turbulence, the effects of sub-gridscale processes on dispersal should be parameterized, in particular when vertical dispersal is to be investigated (e.g., Proehl et al., 2005).

Modeling of vertical dispersal requires sufficient vertical resolution in the velocity fields. The vertical resolution of most “general-purpose” models is typically quite coarse below the thermocline, where the vertical gradients of temperature and salinity, and the corresponding diffusive fluxes, tend to be small. Such models are less suitable for simulating vertical dispersal of tracers and propagules that are associated with strong vertical gradients, such as larvae released at the seabed, which set the diffusive vertical fluxes between adjacent grid cells or isopycnal layers. Within a grid cell or layer, diffusive vertical dispersal in a numerical model is instantaneous due to the standard assumption that any particles contained within a cell are distributed uniformly across its volume. In reality, it may take ~a year for a tracer sheet to diffuse across a vertical distance of 100 m in the deep ocean, away from the immediate vicinity of topography (Ledwell and Watson, 1991; Ledwell et al., 1993, 2011). As a result, simulated vertical dispersal of propagules can be much more rapid than in the real ocean, even in models with accurate diffusive fluxes of heat, salt, oxygen, nutrients, etc. It is noted that fine vertical resolution is typically also required to simulate the small scales associated with the topographic flows discussed above. Such inaccuracies in this advective/diffusive vertical dispersion parameter will be exacerbated if any biological parameters are also used to modify vertical position within a model (e.g., buoyancy or vertical swimming behavior), so care must be used to minimize or acknowledge the error here.

Another important consequence of unresolved sub-gridscale processes in numerical models is that validation with velocity measurements is difficult, as the model velocities represent spatial averages over grid cells. As a result, direct comparisons between observed and modeled velocity time series often show sizable differences, even for high-resolution regional circulation models that have considerable skill in predicting dispersal as validated, for example, with tracer-release experiments (e.g., Proehl et al., 2005; Lavelle et al., 2010; McGillicuddy et al., 2010). In general, Lagrangian data from tracer and dye release experiments, and from float and drifter trajectories are more suitable for validating the dispersal characteristics of a model than Eulerian velocity measurements, because they integrate the effects of all processes affecting dispersal, regardless of their scales. While Lagrangian experiments are expensive and difficult to carry out, especially in the deep ocean, there are data available from previous and on-going deep tracer-release experiments (e.g., Ledwell and Watson, 1991; Ledwell et al., 1993, 2000, 2011; Jackson et al., 2010) and float studies (e.g., Hautala and Riser, 1993; Hogg and Owens, 1999; Argo data http://www.argo.ucsd.edu/) that can be used to validate the dispersal characteristics of large-scale circulation models, at least in some regions of the deep ocean.

The choice of an appropriate OGCM, including subgridscale parameterizations, is key in the process of model parameterization. It should be recognized that OGCMs are generally not designed explicitly to estimate larval dispersal. Consequently, a model designed (and validated) to represent global thermohaline circulation may perform less well within particular regions (Fossette et al., 2012). Further complications arise in the choice of particle tracer. Online particle tracers run natively within the OGCM utilizing the full resolution model output to infer advection and diffusion, but access and computational restrictions can become prohibitive for repeated runs (North et al., 2009; Fossette et al., 2012). Offline particle tracer models (Supplementary Table 1) use outputs from OGCMs, with both the offline models and OGCM outputs being more accessible. However, outputs from temporally and spatially highly resolved OGCMs are often averaged to lower resolutions to reduce the required storage capacity. This averaging can in turn reduce the resolution of captured hydrographic phenomena, e.g., de-trending tides and smoothing eddies, potentially resulting in erroneous trajectory predictions (Putman and He, 2013). Sensitivity analyses can be very informative in terms of the limitations and predictive ability of an OGCM/particle tracer coupling, and can assist in model choices and discourage “black box” model usage (Simons et al., 2013). Coupled with study specific validation, sensitivity analysis is an advisable step prior to settling upon model choice and asking questions of dispersal (North et al., 2009).

Parameterizing the Biology for Biophysical Models in the Deep Sea: The Role of Planktonic Larval Duration

Most marine benthic species exhibit a biphasic life cycle, which includes a pelagic larva, but there are exceptions such as pericarid crustaceans and nematodes that have direct development. Hence, models of the distribution of benthic organisms typically incorporate species-specific biological parameters that account for this potentially dispersive larval phase. The most frequently utilized is the PLD during which larvae are susceptible to physical mixing and advection (Sponaugle et al., 2002; Treml et al., 2008). In spite of its reference to larval development, values of PLD provided in the literature often encompass the entire development between egg release and settlement, although such data would be best described as planktonic propagule duration (PPD; embryonic + larval phase). Because embryos (to late gastrula) may differ from larvae in physical properties (e.g., shape, buoyancy) and swimming capacity, a true PLD (restricted to the larval phase) may be distinguished from PPD, where such data are available (e.g., Brooke and Young, 2009; Selkoe and Toonen, 2011; Mercier et al., 2013). A clearer distinction enables the inclusion of passive vs. active dispersive phases (e.g., egg/embryos vs. larvae), as well as transient planktonic phases in species that undergo parental care or demersal development for a portion of the embryonic or larval phases (e.g., certain gastropods, polychaetes and anthozoans).

Taken as the length of the planktonic phase, PLD has long been a central variable of biophysical models (Sponaugle et al., 2002; Lett et al., 2010; Liggins et al., 2013). The simplest models assume that pelagic propagules are passive and that population connectivity is therefore inversely related to PLD (e.g., in reef fishes; Roberts, 1997). However, the strength of the relationship between PLD and dispersal is being debated (Paulay and Meyer, 2006; Strathmann, 2007; Shanks, 2009; Weersing and Toonen, 2009; Selkoe and Toonen, 2011; Mercier et al., 2013). The emerging view is that dispersal is not only determined by the length of the planktonic phase, but also by circulation processes (e.g., Watson et al., 2010) and larval behavior (e.g., Metaxas and Saunders, 2009; Shanks, 2009; Butler et al., 2011). Hence, the use of PLD in combination with other biological variables is now gaining favor in designing biophysical models of species recruitment and population connectivity (Sponaugle et al., 2002; Levin, 2006; Fiksen et al., 2007; Gilbert et al., 2010; Domingues et al., 2012; Treml et al., 2012; Kough et al., 2013; Nicolle et al., 2013; Nolasco et al., 2013). Fine predictions tend to include many biotic variables. For instance, the LARVAHS model proved to be effective at estimating recruitment success in clams, emphasizing the role of PLD, as well as habitat suitability, larval swimming behavior, wind patterns (seasons), spawning ground location and tidal phase at spawning (Bidegain et al., 2013). Nevertheless, simpler models can be relatively robust. A recent study of invertebrate and fish larvae showed that PLD and depth distribution explained 80% of total variation in dispersal distance, whereas spawning season, and geographic and annual variations in circulation had only marginal effects (Corell et al., 2012). Conversely, differences in reproductive seasons were determined to drive opposite source-sink dynamics in two congeneric mussel species (Carson et al., 2010). Testing various idealized larval behaviors also supported the role of vertical swimming/migration during planktonic development as a key determinant of nearshore settlement site (Drake et al., 2013).

While it may be desirable, the inclusion of several biotic variables, particularly behavioral traits, is generally more difficult for deep-sea than for shallow-water species. Even obtaining reliable estimates of PLDs can present a challenge; however, the coupling of PLD with oceanographic data can provide estimates of the upper bounds of dispersal distances (Young et al., 2012). To date, PLDs have been estimated for deep-sea species using four different methods: (1) larval culture in the laboratory, which presents several challenges in terms of maintaining appropriate rearing conditions, or in the field which may not allow for the completion of the life cycle; (2) computation of PLD from metabolic rates and available energy stores, which is valid only for lecithotrophic larvae; this approach also requires knowledge of the relationship between temperature and metabolic efficiency that for most deep-sea species can only be assumed; (3) tracking of larval cohorts in the plankton, an approach only possible for species with discrete spawning periods; and (4) calculation based on the timings of settlement and spawning times, which also may require back-calculation of settlement time using juvenile growth rates, seldom known for deep-sea species.

Analysis of Existing Data on Planktonic Larval durations in the Deep Sea

We assessed the current knowledge on planktonic larval durations (PLD) for deep-sea species and compared the available estimates with those from shallow-water species, using published observational and experimental PLD values for 305 species belonging to seven marine benthic phyla (Table 1, Supplementary Data): Cnidaria (12), Annelida (25), Sipuncula (1), Mollusca (31), Arthropoda (68), Echinodermata (167), Chordata (1). When multiple PLD values were available in the literature for the same species, only the minimum and maximum values were kept for analyses and used to calculate the median PLD; when only one value was available in the literature that value was used as minimum, maximum and median in all analyses.

TABLE 1

Table 1. Number of species for which PLD has been estimated.

Species were classified according to their bathymetric distribution into shallow (0–200 m), eurybathic (0–>200 m) and deep (>200 m). In total, PLD estimates were available for 212 shallow, 72 eurybathic and 21 deep-sea species (Table 1). Eurybathic—living below 200 m and deep-sea species were further categorized according to their habitat (Figure 1). The relatively high number of species from the sedimentary slope (excluding other specific sub-habitats found on slopes such as cold-water corals, cold seeps; N = 75), for which the PLD has been investigated is the result of numerous reproductive studies on echinoderms with a eurybathic distribution (N = 52); echinoderms are the best-studied group of deep-sea animals in terms of reproduction (Young, 2003) and are the taxonomic group for which many PLD values are available. Apart from these, PLD from species living below 200 m depth has been investigated mostly for polychaetes, molluscs, and crustaceans from chemosynthesis-based habitats (hydrothermal vents and cold seeps, N = 12) since questions related to how these insular and ephemeral habitats are maintained and new sites colonized by larvae have been of much interest in recent decades (Tyler and Young, 1999; Metaxas and Kelly, 2010). We examined differences in PLD among shallow, eurybathic and deep species (available values pooled within each of these three categories) with a Kruskal-Wallis test, followed by Dunn's multiple comparisons using the statistical software GraphPad Prism (version 6.0).

FIGURE 1

Figure 1. Number of deep-sea and eurybathic species for which PLD has been estimated in different deep-sea habitats (SS, sedimentary slope; Sm, seamounts; CWC, cold water coral; CS, cold seeps; FF, food falls; HV, hydrothermal vents).

The ranges in minimum, maximum and median PLD were quite wide for all three bathymetric distributions (Table 2). Statistically significant differences were found between shallow and both eurybathic and deep species; PLD values of eurybathic and deep-sea species were not significantly different from one another (Table 3). Overall, shallow-water species present shorter PLD than eurybathic and deep species (Figure 2, Table 2).

TABLE 2

Table 2. Descriptive statistics of the available minimum, maximum and median PLD values (days) for shallow (N = 212); eurybathic (N = 72) and deep species (N = 20; the value of 660 days determined for Sclerasterias tanneri was not included).

TABLE 3

Table 3. Results of the Kruskal-Wallis and Dunn's multiple comparison tests used to investigate differences of minimum, maximum and median values between shallow, eurybathic and deep species.

FIGURE 2

Figure 2. Mean and standard deviation of minimum (left), maximum (middle) and median (right) PLD values for shallow, eurybathic and deep-sea species. The value of 660 days determined for Sclerasterias tanneri was not included.

PLD estimates have been included in biophysical models to set boundaries on dispersal distances of individual species. However, in the context of spatial planning, connectivity is often not considered for single species, but rather between different areas. Models can estimate the durations that would be required to connect particular metapopulations; however, in most cases there is insufficient data on the distribution of deep-sea species to use this methodology. An alternative approach would be to incorporate in the models a PLD value that guarantees a minimum dispersal distance for a wide range of species. Based on the existing data, we propose that a PLD of 35 and 69 days ensures a minimum distance for 50 and 75%, respectively, of the eurybathic and deep-sea species in our study (Figure 3). Despite its potential utility, the paucity of data points and the high variability in PLD between and within taxonomic groups underscore the limitations of our proposal.

FIGURE 3

Figure 3. Histogram of frequency of minimum PLD estimates available for shallow, eurybathic and deep-sea species. Median, percentile 25 and 75 are shown for PLD values of eurybathic and deep-sea species pooled together. The value of 660 days determined for Sclerasterias tanneri was not included.

The current knowledge on the PLD of deep-sea species is scarce and unevenly distributed between habitats and taxonomic groups (Figure 4), limiting generalizations that are often appealing to decision-makers. More data on the larval ecology is undoubtedly necessary to develop effective conservation strategies for deep-sea habitats and species. Nevertheless, the few biophysical models integrating estimates of larval duration have already generated important information to understand dispersal and connectivity in the deep sea (e.g., Marsh et al., 2001; Yearsley and Sigwart, 2011; Young et al., 2012).

FIGURE 4

Figure 4. Minimum (left) and maximum (right) PLD values available for shallow, eurybathic and deep-sea species of the different invertebrate phyla.

Using Biophysical Models to Predict Dispersal in the Deep Sea

Only a handful of studies have attempted to model larval dispersal in the deep sea, and most of these have been limited by the availability of reliable estimates of biological parameters, including planktonic larval duration, spawning season, and dispersal depth. Consequently, most modeling exercises to date have been experiments that attempt to identify the biological values required to produce a given distributional or genetic pattern. In one of the earliest such studies, Chevaldonné et al. (1997) used a particle flux model to estimate the magnitude of dispersal required by larvae of the vent worm Alvinella pompejana to reach known vent sites on the East Pacific Rise (EPR). Short PLDs and lecithotrophic development were assumed based on circumstantial evidence, and the resulting simulations did not agree with an observed absence of genetic structure over large spatial scales. To explain the absence of genetic structure, the authors invoked spatially variable and frequent geological events, rather than the possibility that dispersal times were longer than assumed. Subsequent embryological studies with the same species (Pradillon et al., 2001, 2005) suggested a mechanism of developmental arrest at cold temperatures that could easily reconcile the genetic data with model predictions and confirming that lecithotrophic development does not necessarily result in short PLD (e.g., Shilling and Manahan, 1994).

Using progressive vector models based on Eulerian current measurements during a single year, Marsh et al. (2001) predicted the dispersal potential of the vent tubeworm Riftia pachyptila on the EPR. Reliable estimates of PLD were based on laboratory and in situ larval rearing, as well as metabolic measurements predicting the depletion rate of internal lipid stores. The models predicted peak dispersal distances at ~103 km, with high cumulative mortality by the terminal (presumably competent) larval stages. In a subsequent study, Brooke and Young (2009) showed that R. pachyptila disperse as unciliated embryos for the first 3 weeks of their 5-week development. If this observation is superimposed on the mortality schedule estimated by the simulations of Marsh et al. (2001), one must conclude that more than half of the individuals are lost from the vent system before they even become larvae. Brooke and Young (2009) also showed that ontogenetic migration by this species is limited by larval tolerance to the warm temperatures and low pressures that prevail in the shallower depths of the water column. The progressive vector approach of Marsh et al. (2001) was extended from 9°N to 13°N on the EPR by Mullineaux et al. (2002) who estimated dispersal distances for theoretical larvae with both shorter and longer PLDs than those determined empirically for Riftia pachyptila.

The dispersal of simulated larvae from hydrothermal vents as either passive particles or “balloonists” that migrate ontogenetically to and from the upper water column was modeled by McGillicuddy et al. (2010). The model was both driven and validated with data from moored current meters and particle trajectories based on larval releases every 12 min in all seasons of the year. The model assumed a 30-day dispersal period as predicted by Marsh et al. (2001) for Riftia pachyptila. This modeling exercise yielded several important conclusions, including the observation that larvae released at the crest of the EPR dispersed greater distances near the sea floor than higher in the water column. This result suggests that ontogenetic migration is not a viable strategy for increasing dispersal in this system and is in agreement with the physiological tolerances reported by Brooke and Young (2009). McGillicuddy et al. (2010) explained this pattern with reference to observed vigorous along-flank flows and high flow variability near the ridge axis, but were not able to test these hypotheses empirically with their moored current meter data. Another important conclusion was that the likelihood of successful settlement varied significantly with time of year. Larvae “released” in early summer were much less likely to find themselves near a suitable settlement site than ones released in winter. Lastly, the model showed that the period of precompetency may play a highly significant role in successful dispersal. As predicted by Jackson and Strathmann (1981), larvae that spend a higher percentage of their time being able to undergo metamorphosis will be more likely to encounter suitable habitat during that period. In the case of R. pachyptila, however, the developmental observations (Marsh et al., 2001; Brooke and Young, 2009) indicate that metamorphic competency is not attained until the very end of larval life.

In non-chemosynthetic habitats of the deep sea, Yearsley and Sigwart (2011) used a biophysical model to estimate connectivity among known populations of deep-sea chitons in archipelagos of the Southwest Pacific. Because there were no biological parameters available for any of the species examined, the models were based on several assumptions inferred from shallow-water species with known larval development. Because all known chitons from shallow water have lecithotrophic development, it was reasonably assumed that deep-water species would have the same developmental mode. Planktonic larval duration is unknown for any deep-sea chiton, and this important parameter was extrapolated from known shallow-water PLDs at relatively high temperatures. The assumption that the temperature/PLD curve is the same for shallow- and deep-water species is untested for this group and, based on observations in echinoderms and other phyla, may be unwarranted (Young, unpublished data). It was also assumed that dispersal occurred at approximately the same depth as where the adults are found. The observed dispersal kernels were likely too small to maintain connectivity among known populations within a single generation, suggesting the existence of additional undiscovered metapopulations that serve as stepping stones. The accuracy of the biological parameters used is questionable in the absence of empirical data, but because the authors intentionally overestimated dispersal by adding in a pre-larval dispersal phase of 50 days, the major conclusion may be robust.

A recent study used Lagrangian (LTRANS) modeling with regional ocean models system (ROMS) to predict dispersal trajectories of larvae with known PLD, either obtained from larval rearing studies or estimated indirectly from other types of biological data (Young et al., 2012). Dispersal trajectories were estimated for seven cold seep and non-seep species at each of two depths, originating from known adult locations off Barbados, in the Tongue of the Ocean Bahamas, and at seeps in the Gulf of Mexico. Planktonic larval durations ranged from 21 days in cold-seep tubeworms (Lamellibrachia luymesi) to nearly 2 years in bipinnaria larvae of the asteroid Sclerasterias tanneri. For those species with long PLDs, the models showed significantly greater dispersal at 100 m depth than at 500 m depth, although the actual ontogenetic movements are unknown for all species. For some species, simulations were run more than once a year, revealing highly significant temporal effects in dispersal distance and underscoring the importance of spawning time. Data on PLD had been measured for two Bahamian echinoids, Stylocidaris lineata and Cidaris blakei (Bennett et al., 2012), and an extended larval life had been documented for the sea star, which provided a conservative estimate of dispersal time. PLD also had been estimated for the tubeworm (Young et al., 1996) on the basis of larval rearing. For the other species, PLD had to be estimated from indirect methods. Based on data on the timing of spawning and seasonal recruitment, and on in situ growth rates of juveniles of Bathymodiolus childressi (Arellano and Young, 2009), settlement times of individuals were back-calculated. The number of days since the previous spawning peak to settlement was then determined as the PLD. A similar method was used for the sipunculan Phascolosoma turnerae, known to have a strongly seasonal spawning peak and for which a growth curve had been generated for laboratory-reared juveniles (Rice et al., 2012). Two species of seep molluscs, Bathymodiolus childressi and Bathynerita naticoidea, were collected in the water column during some seasons of the year. Because gametogenic cycles had already been documented for both species, it was possible to infer the approximate ages of the plankton-collected larvae to assist with the estimation of PLD.

The study of Young et al. (2012) demonstrates the value of incorporating various types of biological measurements into the estimation of dispersal distances. A major shortfall of this and all other deep-sea dispersal models is the absence of information on dispersal depth as determined by actual vertical distributions in the plankton (Arellano et al., 2014), or estimated by swimming speed, egg buoyancy, direction of swimming, and physiological tolerances.

Main Data Gaps for Estimating Connectivity in the Deep Sea

Estimating connectivity in marine ecosystems requires the understanding of the biological and physical processes regulating larval dispersal, settlement and recruitment, which is hindered by the small size of larvae coupled with the vast and complex fluid environment. In the deep sea, the problem is magnified because of inherent barriers of accessibility and sampling constraints. Biophysical modeling approaches, an established technique in marine connectivity research, which incorporate key physical dynamics and biological traits, provide a partial solution to this problem. In the deep sea, however, knowledge gaps in both the physical and biological components are delaying the effective use of finely tuned biophysical models.

The physical components of biophysical models are limited by scale and computational capacity. Large scale processes are well understood and parameterized within model equations, but sub-mesoscale ocean physics is an ongoing area of research. In any case, the phenomena at the small scales relevant to a larva are prohibitively expensive to parameterize within a large spatial scale model. For this reason, a sub-gridscale parameterisation is usually considered an adequate enough approximation for most purposes although this is difficult to estimate in itself. There are additional conflicts of resolution, such as poorly represented topography, and consequential flow modifications, that result from low resolution in bathymetry data. The best high-resolution models rely upon high quality bathymetry which is costly to acquire over large spatial scales, particularly at great depth and distance from shore. Accurate model performance requires the assimilation of and validation by high resolution observational data over large temporal and spatial scales.

Vertical velocity remains the most elusive component of both observed and modeled velocity fields, often reduced to secondary calculation in line with the conservation of energy. Without improved understanding of vertical velocities, the potential for passive vertical migration of larvae also remains elusive. In the end, numerical models should be taken for what they are—simplified approximations of reality offered as a best guess of average hydrographic conditions in the area concerned.

Biological processes that control larval dispersal include the reproductive effort of adults, which determines the timing and number of larvae in the water column, and larval development and behavior. Both these components determine how larvae interact with the oceanic circulation and influence the timing, distance and trajectory of larvae among habitats. Reproductive biology has only been studied for a small fraction of deep-sea species and most of the utilized knowledge on parental investment, egg size and fecundity has been extrapolated from a limited number of samples. Because time-series analyses are rare, measurements on reproductive synchrony and periodicity remains elusive for most deep-sea species. As a result, information on the initial factors controlling larval dispersal, “how many” and “when” larvae enter the water column is largely unavailable.

“How long” larvae spend in the water column is defined by the planktonic larval duration (PLD), a fundamental parameter in dispersal models that is unknown for most deep-sea species and highly variable among those for which it has been estimated. One of the biggest obstacles in acquiring accurate PLD estimates is the difficulty in culturing embryos and larvae in the laboratory. In addition, it remains difficult to assess whether and how pressure and temperature conditions may affect PLDs determined in culture. Studies on early life-history stages are scarce and hampered by sampling difficulties and the general lack of facilities available for long-term maintenance under appropriate conditions.

Lastly, “where” larvae are positioned in the water column is largely determined by their swimming behavior. Most larvae are poor horizontal swimmers, but they can alter their vertical positions actively through vertical swimming and/or passively through differential buoyancy; for deep-sea larvae, changes in vertical position may result in major changes in temperature and pressure, in turn resulting in major consequences in terms metabolic and feeding rates, and other vital processes. Further, because of the vertical structure in current velocity different dispersal pathways will also be affected. Ontogenetic changes in anatomical features can indirectly affect larval behavior, motility (speed) within the water column and nutritional reserves (or feeding regime), in turn influencing vertical positioning. However, detailed larval development characterization has only been achieved for few deep-sea species. Larval behavior and development studies are not only constrained by the difficulty in rearing deep-sea larvae in the laboratory, but also by the difficulty in identifying field collected larvae.

Recent technological advances can help expand our knowledge of reproductive and larval biology for deep-sea species. For example, increasing availability and development of new technologies, including pressure and temperature-controlled sampling vessels and holding facilities (e.g., Pradillon et al., 2004; Mestre et al., 2009, 2013; Ravaux et al., 2013) can facilitate the quantification of PLD and swimming behavior. Moreover, the development of sampling systems associated with larval identification by high-throughput molecular techniques allows direct observation of distribution of deep-sea larvae in the water column.

While many gaps exist in our ability to collect empirical data and use predictive models on the factors that regulate dispersal in the deep sea, anthropogenic pressure on this habitat for resource extraction is rapidly increasing. For most of the deep-sea regions currently under threat of major disturbances associated with fishing and mining, our limited understanding of the resident species and communities and the mechanism that regulate them can compromise our ability to manage them sustainably. Despite the cost in filling our knowledge gaps, there is an urgent need to do so.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fmars.2015.00006/abstract

References

Ardron, J., Arnaud-Haond, S., Beaudoin, Y., Bezaury, J., Billet, D., Boland, G., et al. (2011). Environmental Management of Deep-Sea Chemosynthetic Ecosystems: Justification of and Considerations for a Spatially Based Approach. Technical Report of the International Seabed Authority. Avalible online at: https://hal.archives-ouvertes.fr/hal-00667838

Arellano, S. M., Van Gaest, A. L., Johnson, S. B., Vrijenhoek, R. C., and Young, C. M. (2014). Larvae from deep-sea methane seeps disperse in surface waters. Proc. R. Soc. B 281:20133276. doi: 10.1098/rspb.2013.3276