How to select networks of marine protected areas for multiple species with different dispersal strategies

Indicators to assess temporal variability in marine connectivity processes: A semi-theoretical approach

Oceanographic connectivity in an effective network of protected areas is crucial for restoring and stabilising marine populations. However, temporal variability in connectivity is rarely considered as a criterion in designing and evaluating marine conservation planning. In this study, indicators were defined to characterise the temporal variability in occurrence, flux, and frequency of connectivity in a northwestern Mediterranean Sea area. Indicators were tested on semi-theoretically-estimated connections provided by the runs of a passive particle transport model in a climatological year and in three years between 2006-2020, showing large deviation from the climatological year. The indicators allowed comparing the temporal variability in connectivity of four zones, highlighted differences in connectivity due to their locations and the mesoscale hydrodynamics, and identified areas that require further investigation. The three indicators also showed that the temporal variability in connectivity was influenced by the duration and depth of particle transport, although no consistent pattern was observed in the indicator variations of the compared zones. Provided that specific objectives will be given when parameterising transport models (i.e., selection of focus species and time period), indicators of temporal variability in connectivity have potential to support spatial conservation planning, prioritise the protection of marine resources, and measure the effectiveness of Marine Protected Areas, in line with a long-term vision of ocean management.

Unlocking the secret life of blue mussels: Exploring connectivity in the Skagerrak through biophysical modeling and population genomics

Knowledge of functional dispersal barriers in the marine environment can be used to inform a wide variety of management actions, such as marine spatial planning, restoration efforts, fisheries regulations, and invasive species management. Locations and causes of dispersal barriers can be studied through various methods, including movement tracking, biophysical modeling, demographic models, and genetics. Combining methods illustrating potential dispersal, such as biophysical modeling, with realized dispersal through, e.g., genetic connectivity estimates, provides particularly useful information for teasing apart potential causes of observed barriers. In this study, we focus on blue mussels (Mytilus edulis) in the Skagerrak-a marginal sea connected to the North Sea in Northern Europe-and combine biophysical models of larval dispersal with genomic data to infer locations and causes of dispersal barriers in the area. Results from both methods agree; patterns of ocean currents are a major structuring factor in the area. We find a complex pattern of source-sink dynamics with several dispersal barriers and show that some areas can be isolated despite an overall high dispersal capability. Finally, we translate our finding into management advice that can be used to sustainably manage this ecologically and economically important species in the future.

Ecological connectivity of the marine protected area network in the Baltic Sea, Kattegat and Skagerrak: Current knowledge and management needs

Marine protected areas (MPAs) have become a key component of conservation and fisheries management to alleviate anthropogenic pressures. For MPA networks to efficiently promote persistence and recovery of populations, ecological connectivity, i.e. dispersal and movement of organisms and material across ecosystems, needs to be taken into account. To improve the ecological coherence of MPA networks, there is hence a need to evaluate the connectivity of species spreading through active migration and passive dispersal. We reviewed knowledge on ecological connectivity in the Baltic Sea, Kattegat and Skagerrak in the northeast Atlantic and present available information on species-specific dispersal and migration distances. Studies on genetic connectivity are summarised and discussed in relation to dispersal-based analyses. Threats to ecological connectivity, limiting dispersal of populations and lowering the resilience to environmental change, were examined. Additionally, a review of studies evaluating the ecological coherence of MPA networks in the Baltic Sea, Kattegat and Skagerrak was performed, and suggestions for future evaluations to meet management needs are presented.

Introduction to the theme issue 'Species' ranges in the face of changing environments'

Understanding where, when and how species' ranges will be modified is both a fundamental problem and essential to predicting how spatio-temporal environmental changes in abiotic and biotic factors impact biodiversity. Notably, different species may respond disparately to similar environmental changes: some species may overcome an environmental change only with difficulty or not at all, while other species may readily overcome the same change. Ranges may contract, expand or move. The drivers and consequences of this variability in species' responses remain puzzling. Importantly, changes in a species' range creates feedbacks to the environmental conditions, populations and communities in its previous and current range, rendering population genetic, population dynamic and community processes inextricably linked. Understanding these links is critical in guiding biodiversity management and conservation efforts. This theme issue presents current thinking about the factors and mechanisms that limit and/or modify species' ranges. It also outlines different approaches to detect changes in species' distributions, and illustrates cases of range modifications in several taxa. Overall, this theme issue highlights the urgency of understanding species' ranges but shows that we are only just beginning to disentangle the processes involved. One way forward is to unite ecology with evolutionary biology and empirical with modelling approaches. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Reconciling seascape genetics and fisheries science in three codistributed flatfishes

Uncertainty hampers innovative mixed-fisheries management by the scales at which connectivity dynamics are relevant to management objectives. The spatial scale of sustainable stock management is species-specific and depends on ecology, life history and population connectivity. One valuable approach to understand these spatial scales is to determine to what extent population genetic structure correlates with the oceanographic environment. Here, we compare the level of genetic connectivity in three codistributed and commercially exploited demersal flatfish species living in the North East Atlantic Ocean. Population genetic structure was analysed based on 14, 14 and 10 neutral DNA microsatellite markers for turbot, brill and sole, respectively. We then used redundancy analysis (RDA) to attribute the genetic variation to spatial (geographical location), temporal (sampling year) and oceanographic (water column characteristics) components. The genetic structure of turbot was composed of three clusters and correlated with variation in the depth of the pycnocline, in addition to spatial factors. The genetic structure of brill was homogenous, but correlated with average annual stratification and spatial factors. In sole, the genetic structure was composed of three clusters, but was only linked to a temporal factor. We explored whether the management of data poor commercial fisheries, such as in brill and turbot, might benefit from population-specific information. We conclude that the management of fish stocks has to consider species-specific genetic structures and may benefit from the documentation of the genetic seascape and life-history traits.

Integrating genetics, biophysical, and demographic insights identifies critical sites for seagrass conservation

The eelgrass Zostera marina is an important foundation species of coastal areas in the Northern Hemisphere, but is continuing to decline, despite management actions. The development of new management tools is therefore urgent in order to prioritize limited resources for protecting meadows most vulnerable to local extinctions and identifying most valuable present and historic meadows to protect and restore, respectively. We assessed 377 eelgrass meadows along the complex coastlines of two fjord regions on the Swedish west coast-one is currently healthy and the other is substantially degraded. Shoot dispersal for all meadows was assessed with Lagrangian biophysical modeling (scale: 100-1,000 m) and used for barrier analysis and clustering; a subset (n = 22) was also assessed with population genetic methods (20 microsatellites) including diversity, structure, and network connectivity. Both approaches were in very good agreement, resulting in seven subpopulation groupings or management units (MUs). The MUs correspond to a spatial scale appropriate for coastal management of "waterbodies" used in the European Water Framework Directive. Adding demographic modeling based on the genetic and biophysical data as a third approach, we are able to assess past, present, and future metapopulation dynamics to identify especially vulnerable and valuable meadows. In a further application, we show how the biophysical approach, using eigenvalue perturbation theory (EPT) and distribution records from the 1980s, can be used to identify lost meadows where restoration would best benefit the present metapopulation. The combination of methods, presented here as a toolbox, allows the assessment of different temporal and spatial scales at the same time, as well as ranking of specific meadows according to key genetic, demographic and ecological metrics. It could be applied to any species or region, and we exemplify its versatility as a management guide for eelgrass along the Swedish west coast.

Incorporating the Connectivity Timescale in Metapopulation Partitioning

The often complex spatial patterns of propagule dispersal across a metapopulation present a challenge for species management, motivating efforts to represent the connectivity in simpler but meaningful ways. The reduction of complexity may be achieved by partitioning the metapopulation into groups of highly connected patches called "subpopulations." To have relevance for management, these subunits must be defined from ecological or evolutionary principles. The probabilities of dispersal-mediated propagule interchange between sites, commonly organized into a connectivity matrix, entail a timescale that is usually ignored in subpopulation analyses, limiting their utility and possibly leading to misinterpretation and wrong management decisions. Recognition of the essentially dynamical role played by metapopulation connectivity naturally leads to the incorporation of the generational timescale into the partitioning analysis. An algorithm is proposed to determine the subpopulations-both their cardinality and their composition-as a function of the generational timescale and of a limiting probability of connection, illustrated with a novel empirical estimate of mesopelagic connectivity. The proposed framework allows the unambiguous determination of the timescales corresponding to dispersal barriers and the identification of effective ecological units across the spectrum of management-relevant time horizons.

Spatial genetic structure in a crustacean herbivore highlights the need for local considerations in Baltic Sea biodiversity management

Incorporating species' eco-evolutionary responses to human-caused disturbances remains a challenge in marine management efforts. A prerequisite is knowledge of geographic structure and scale of genetic diversity and connectivity-the so-called seascape genetic patterns. The Baltic Sea is an excellent model system for studies linking seascape genetics with effects of anthropogenic stress. However, seascape genetic patterns in this area are only described for a few species and are completely unknown for invertebrate herbivores, which constitute a critical part of the ecosystem. This information is crucial for sustainable management, particularly under future scenarios of rapid environmental change. Here, we investigate the population genetic structure among 31 locations throughout the Baltic Sea, of which 45% were located in marine protected areas, in one of the most important herbivores of this region, the isopod crustacean Idotea balthica, using an array of 33,774 genome-wide SNP markers derived from 2b-RAD sequencing. In addition, we generate a biophysical connectivity matrix for I. balthica from a combination of oceanographic current models and estimated life history traits. We find population structure on scales of hundreds of kilometers across the Baltic Sea, where genomic patterns in most cases closely match biophysical connectivity, indicating passive transport with oceanographic currents as an important mean of dispersal in this species. We also find a reduced genetic diversity in terms of heterozygosity along the main salinity gradient of the Baltic Sea, suggesting periods of low population size. Our results provide crucial information for the management of a key ecosystem species under expected changes in temperature and salinity following global climate change in a marine coastal area.

Seascape genetics and biophysical connectivity modelling support conservation of the seagrass Zostera marina in the Skagerrak-Kattegat region of the eastern North Sea

Maintaining and enabling evolutionary processes within meta-populations are critical to resistance, resilience and adaptive potential. Knowledge about which populations act as sources or sinks, and the direction of gene flow, can help to focus conservation efforts more effectively and forecast how populations might respond to future anthropogenic and environmental pressures. As a foundation species and habitat provider, Zostera marina (eelgrass) is of critical importance to ecosystem functions including fisheries. Here, we estimate connectivity of Z. marina in the Skagerrak-Kattegat region of the North Sea based on genetic and biophysical modelling. Genetic diversity, population structure and migration were analysed at 23 locations using 20 microsatellite loci and a suite of analytical approaches. Oceanographic connectivity was analysed using Lagrangian dispersal simulations based on contemporary and historical distribution data dating back to the late 19th century. Population clusters, barriers and networks of connectivity were found to be very similar based on either genetic or oceanographic analyses. A single-generation model of dispersal was not realistic, whereas multigeneration models that integrate stepping-stone dispersal and extant and historic distribution data were able to capture and model genetic connectivity patterns well. Passive rafting of flowering shoots along oceanographic currents is the main driver of gene flow at this spatial-temporal scale, and extant genetic connectivity strongly reflects the "ghost of dispersal past" sensu Benzie, 1999. The identification of distinct clusters, connectivity hotspots and areas where connectivity has become limited over the last century is critical information for spatial management, conservation and restoration of eelgrass.

The Baltic Sea as a time machine for the future coastal ocean

Coastal global oceans are expected to undergo drastic changes driven by climate change and increasing anthropogenic pressures in coming decades. Predicting specific future conditions and assessing the best management strategies to maintain ecosystem integrity and sustainable resource use are difficult, because of multiple interacting pressures, uncertain projections, and a lack of test cases for management. We argue that the Baltic Sea can serve as a time machine to study consequences and mitigation of future coastal perturbations, due to its unique combination of an early history of multistressor disturbance and ecosystem deterioration and early implementation of cross-border environmental management to address these problems. The Baltic Sea also stands out in providing a strong scientific foundation and accessibility to long-term data series that provide a unique opportunity to assess the efficacy of management actions to address the breakdown of ecosystem functions. Trend reversals such as the return of top predators, recovering fish stocks, and reduced input of nutrient and harmful substances could be achieved only by implementing an international, cooperative governance structure transcending its complex multistate policy setting, with integrated management of watershed and sea. The Baltic Sea also demonstrates how rapidly progressing global pressures, particularly warming of Baltic waters and the surrounding catchment area, can offset the efficacy of current management approaches. This situation calls for management that is (i) conservative to provide a buffer against regionally unmanageable global perturbations, (ii) adaptive to react to new management challenges, and, ultimately, (iii) multisectorial and integrative to address conflicts associated with economic trade-offs.

Designing connected marine reserves in the face of global warming

Marine reserves are widely used to protect species important for conservation and fisheries and to help maintain ecological processes that sustain their populations, including recruitment and dispersal. Achieving these goals requires well-connected networks of marine reserves that maximize larval connectivity, thus allowing exchanges between populations and recolonization after local disturbances. However, global warming can disrupt connectivity by shortening potential dispersal pathways through changes in larval physiology. These changes can compromise the performance of marine reserve networks, thus requiring adjusting their design to account for ocean warming. To date, empirical approaches to marine prioritization have not considered larval connectivity as affected by global warming. Here, we develop a framework for designing marine reserve networks that integrates graph theory and changes in larval connectivity due to potential reductions in planktonic larval duration (PLD) associated with ocean warming, given current socioeconomic constraints. Using the Gulf of California as case study, we assess the benefits and costs of adjusting networks to account for connectivity, with and without ocean warming. We compare reserve networks designed to achieve representation of species and ecosystems with networks designed to also maximize connectivity under current and future ocean-warming scenarios. Our results indicate that current larval connectivity could be reduced significantly under ocean warming because of shortened PLDs. Given the potential changes in connectivity, we show that our graph-theoretical approach based on centrality (eigenvector and distance-weighted fragmentation) of habitat patches can help design better-connected marine reserve networks for the future with equivalent costs. We found that maintaining dispersal connectivity incidentally through representation-only reserve design is unlikely, particularly in regions with strong asymmetric patterns of dispersal connectivity. Our results support previous studies suggesting that, given potential reductions in PLD due to ocean warming, future marine reserve networks would require more and/or larger reserves in closer proximity to maintain larval connectivity.

Sensitivity of marine protected area network connectivity to atmospheric variability

International efforts are underway to establish well-connected systems of marine protected areas (MPAs) covering at least 10% of the ocean by 2020. But the nature and dynamics of ocean ecosystem connectivity are poorly understood, with unresolved effects of climate variability. We used 40-year runs of a particle tracking model to examine the sensitivity of an MPA network for habitat-forming cold-water corals in the northeast Atlantic to changes in larval dispersal driven by atmospheric cycles and larval behaviour. Trajectories of Lophelia pertusa larvae were strongly correlated to the North Atlantic Oscillation (NAO), the dominant pattern of interannual atmospheric circulation variability over the northeast Atlantic. Variability in trajectories significantly altered network connectivity and source-sink dynamics, with positive phase NAO conditions producing a well-connected but asymmetrical network connected from west to east. Negative phase NAO produced reduced connectivity, but notably some larvae tracked westward-flowing currents towards coral populations on the mid-Atlantic ridge. Graph theoretical metrics demonstrate critical roles played by seamounts and offshore banks in larval supply and maintaining connectivity across the network. Larval longevity and behaviour mediated dispersal and connectivity, with shorter lived and passive larvae associated with reduced connectivity. We conclude that the existing MPA network is vulnerable to atmospheric-driven changes in ocean circulation.

Sensitivity of marine protected area network connectivity to atmospheric variability

International efforts are underway to establish well-connected systems of marine protected areas (MPAs) covering at least 10% of the ocean by 2020. But the nature and dynamics of ocean ecosystem connectivity are poorly understood, with unresolved effects of climate variability. We used 40-year runs of a particle tracking model to examine the sensitivity of an MPA network for habitat-forming cold-water corals in the northeast Atlantic to changes in larval dispersal driven by atmospheric cycles and larval behaviour. Trajectories of Lophelia pertusa larvae were strongly correlated to the North Atlantic Oscillation (NAO), the dominant pattern of interannual atmospheric circulation variability over the northeast Atlantic. Variability in trajectories significantly altered network connectivity and source-sink dynamics, with positive phase NAO conditions producing a well-connected but asymmetrical network connected from west to east. Negative phase NAO produced reduced connectivity, but notably some larvae tracked westward-flowing currents towards coral populations on the mid-Atlantic ridge. Graph theoretical metrics demonstrate critical roles played by seamounts and offshore banks in larval supply and maintaining connectivity across the network. Larval longevity and behaviour mediated dispersal and connectivity, with shorter lived and passive larvae associated with reduced connectivity. We conclude that the existing MPA network is vulnerable to atmospheric-driven changes in ocean circulation.