Projecting marine species range shifts from only temperature can mask climate vulnerability

Climate change is causing range shifts in many marine species, with implications for biodiversity and fisheries. Previous research has mainly focused on how species' ranges will respond to changing ocean temperatures, without accounting for other environmental covariates that could affect future distribution patterns. Here, we integrate habitat suitability modeling approaches, a high‐resolution global climate model projection, and detailed fishery‐independent and ‐dependent faunal datasets from one of the most extensively monitored marine ecosystems—the U.S. Northeast Shelf. We project the responses of 125 species in this region to climate‐driven changes in multiple oceanographic factors (e.g., ocean temperature, salinity, sea surface height) and seabed characteristics (i.e., rugosity and depth). Comparing model outputs based on ocean temperature and seabed characteristics to those that also incorporated salinity and sea surface height (proxies for primary productivity and ocean circulation features), we explored how an emphasis on ocean temperature in projecting species' range shifts can impact assessments of species' climate vulnerability. We found that multifactor habitat suitability models performed better in explaining and predicting species historical distribution patterns than temperature‐based models. We also found that multifactor models provided more concerning assessments of species' future distribution patterns than temperature‐based models, projecting that species' ranges will largely shift northward and become more contracted and fragmented over time. Our results suggest that using ocean temperature as a primary determinant of range shifts can significantly alter projections, masking species' climate vulnerability, and potentially forestalling proactive management.     

Offshore Habitat Preference of Overwintering Juvenile and Adult Black Sea Bass,Centropristis striata,and the Relationship to Year-Class Success

Black sea bass (Centropristis striata) migrations are believed to play a role in overwinter survival and connectivity between juvenile and adult populations. This study investigated oceanographic drivers of winter habitat choice and regional differences between populations of juvenile and adult black sea bass. Trends in cohort strength, as a result of juvenile survival, were also identified. Oceanographic and fisheries survey data were analyzed using generalized additive models. Among the oceanographic variables investigated, salinity was the main driver in habitat selection with an optimal range of 33–35 practical salinity units (PSU) for both juveniles and adults. Preferred temperature ranges varied between juveniles and adults, but held a similar minimum preference of >8°C. Salinity and temperature ranges also differed by regions north and south of Hudson Canyon. Shelf water volume had less of an effect than temperature or salinity, but showed an overall negative relationship with survey catch. The effect of winter conditions on juvenile abundance was also observed across state and federal survey index trends. A lack of correlation observed among surveys in the fall paired with a strong correlation in the spring identifies the winter period as a factor determining year-class strength of new recruits to the population. A rank order analysis of spring indices identified three of the largest year classes occurring during years with reduced shelf water volumes, warmer winter shelf waters, and a 34 PSU isohaline aligned farther inshore. While greater catches of black sea bass in the northwest Atlantic Ocean remain south of Hudson Canyon, the species’ range has expanded north in recent years.     

Climate and habitat configuration limit range expansion and patterns of dispersal in a non‐native lizard

Invasive species are one of the main causes of biodiversity loss worldwide. As introduced, populations increase in abundance and geographical range, so does the potential for negative impacts on native communities. As such, there is a need to better understand the processes driving range expansion as species become established in recipient landscapes. Through an investigation into capacity for population growth and range expansion of introduced populations of a non‐native lizard (Podarcis muralis), we aimed to demonstrate how multi‐scale factors influence spatial spread, population growth, and invasion potential in introduced species. We collated location records of P. muralis presence in England, UK through data collected from field surveys and a citizen science campaign. We used these data as input for presence‐background models to predict areas of climate suitability at a national‐scale (5 km resolution), and fine‐scale habitat suitability at the local scale (2 m resolution). We then integrated local models into an individual‐based modeling platform to simulate population dynamics and forecast range expansion for 10 populations in heterogeneous landscapes. National‐scale models indicated climate suitability has restricted the species to the southern parts of the UK, primarily by a latitudinal cline in overwintering conditions. Patterns of population growth and range expansion were related to differences in local landscape configuration and heterogeneity. Growth curves suggest populations could be in the early stages of exponential growth. However, annual rates of range expansion are predicted to be low (5–16 m). We conclude that extensive nationwide range expansion through secondary introduction is likely to be restricted by currently unsuitable climate beyond southern regions of the UK. However, exponential growth of local populations in habitats providing transport pathways is likely to increase opportunities for regional expansion. The broad habitat niche of P. muralis, coupled with configuration of habitat patches in the landscape, allows populations to increase locally with minimal dispersal.     

Fine‐scale foraging habitat selection by two diving central place foragers in the Northeast Atlantic

Habitat selection and spatial usage are important components of animal behavior influencing fitness and population dynamic. Understanding the animal–habitat relationship is crucial in ecology, particularly in developing strategies for wildlife management and conservation. As this relationship is governed by environmental features and intra‐ and interspecific interactions, habitat selection of a population may vary locally between its core and edges. This is particularly true for central place foragers such as gray and harbor seals, where, in the Northeast Atlantic, the availability of habitat and prey around colonies vary at local scale. Here, we study how foraging habitat selection may vary locally under the influence of physical habitat features. Using GPS/GSM tags deployed at different gray and harbor seals’ colonies, we investigated spatial patterns and foraging habitat selection by comparing trip characteristics and home‐range similarities and fitting GAMMs to seal foraging locations and environmental data. To highlight the importance of modeling habitat selection at local scale, we fitted individual models to colonies as well as a global model. The global model suffered from issues of homogenization, while colony models showed that foraging habitat selection differed markedly between regions for both species. Despite being capable of undertaking far‐ranging trips, both gray and harbor seals selected their foraging habitat depending on local availability, mainly based on distance from the last haul‐out and bathymetry. Distance from shore and tidal current also influenced habitat preferences. Results suggest that local conditions have a strong influence on population spatial ecology, highlighting the relevance of processes occurring at fine geographical scale consistent with management within regional units.     

Thermal ecology and baseline energetic requirements of a large‐bodied ectotherm suggest resilience to climate change

1. Most studies on how rising temperatures will impact terrestrial ectotherms have focused on single populations or multiple sympatric species. Addressing the thermal and energetic implications of climatic variation on multiple allopatric populations of a species will help us better understand how a species may be impacted by altered climates.
2. We used eight years of thermal and behavioral data collected from four populations of Pacific rattlesnakes (Crotalus oreganus) living in climatically distinct habitat types (inland and coastal) to determine the field‐active and laboratory‐preferred body temperatures, thermoregulatory metrics, and maintenance energetic requirements of snakes from each population.
3. Physical models showed that thermal quality was best at coastal sites, but inland snakes thermoregulated more accurately despite being in more thermally constrained environments. Projected increases of 1 and 2°C in ambient temperature result in an increase in overall thermal quality at both coastal and inland sites.
4. Population differences in modeled standard metabolic rate estimates were driven by body size and not field‐active body temperature, with inland snakes requiring 1.6× more food annually than coastal snakes.
5. All snakes thermoregulated with high accuracy, suggesting that small increases in ambient temperature are unlikely to impact the maintenance energetic requirements of individual snakes and that some species of large‐bodied reptiles may be robust to modest thermal perturbations under conservative climate change predictions.

​     

An open spatial capture–recapture model for estimating density,movement, and population dynamics from line‐transect surveys

The purpose of many wildlife population studies is to estimate density, movement, or demographic parameters. Linking these parameters to covariates, such as habitat features, provides additional ecological insight and can be used to make predictions for management purposes. Line‐transect surveys, combined with distance sampling methods, are often used to estimate density at discrete points in time, whereas capture–recapture methods are used to estimate movement and other demographic parameters. Recently, open population spatial capture–recapture models have been developed, which simultaneously estimate density and demographic parameters, but have been made available only for data collected from a fixed array of detectors and have not incorporated the effects of habitat covariates. We developed a spatial capture–recapture model that can be applied to line‐transect survey data by modeling detection probability in a manner analogous to distance sampling. We extend this model to a) estimate demographic parameters using an open population framework and b) model variation in density and space use as a function of habitat covariates. The model is illustrated using simulated data and aerial line‐transect survey data for North Atlantic right whales in the southeastern United States, which also demonstrates the ability to integrate data from multiple survey platforms and accommodate differences between strata or demographic groups. When individuals detected from line‐transect surveys can be uniquely identified, our model can be used to simultaneously make inference on factors that influence spatial and temporal variation in density, movement, and population dynamics.     

Fish Species and Community Distributions as Proxies for Seafloor Habitat Distributions: The Stellwagen Bank National Marine Sanctuary Example (Northwest Atlantic, Gulf Of Maine)

Defining the habitats of fishes and associated fauna on outer continental shelves is problematic given the paucity of data on the actual types and distributions of seafloor habitats. However many regions have good data on the distributions of fishes from resource surveys or catch statistics because of the economic importance of the fisheries. Fish distribution data (species or communities) have been used as a proxy for the distribution of habitats to develop precautionary conservation strategies for habitat protection (e.g., marine protected areas, fishing gear restrictions). In this study we assessed the relationships between the distributions of fish communities and species derived from trawl survey data with the spatial distribution of sediment types determined by sampling and acoustic reflectance derived from multibeam sonar surveys in Stellwagen Bank National Marine Sanctuary. Fish communities were correlated with reflectance values but all communities did not occur in unique sediment types. This suggests that use of community distributions as proxies for habitats should include the caveat that a greater number of communities within an area could indicate a greater range of habitat types. Single species distributions showed relationships between abundance and reflectance values. Trawl catches with low abundances had wide variations in reflectance values while those with high abundances had narrower ranges indicating habitat affinities. Significant non-random frequency-dependent relationships were observed for 17 of 20 species although only 12 of 20 species had significant relationships based on rank correlation. These results suggest that species distributions based on trawl survey data can be used as proxies for the distribution of seafloor habitats. Species with known habitat associations can be used to infer habitat requirements of co-occurring species and can be used to identify a range of habitat types.     

Plastic adjustments of biparental care behavior across embryonic development under elevated temperature in a marine ectotherm

Phenotypic plasticity in parental care investment allows organisms to promptly respond to rapid environmental changes by potentially benefiting offspring survival and thus parental fitness. To date, a knowledge gap exists on whether plasticity in parental care behaviors can mediate responses to climate change in marine ectotherms. Here, we assessed the plasticity of parental care investment under elevated temperatures in a gonochoric marine annelid with biparental care, Ophryotrocha labronica, and investigated its role in maintaining the reproductive success of this species in a warming ocean. We measured the time individuals spent carrying out parental care activities across three phases of embryonic development, as well as the hatching success of the offspring as a proxy for reproductive success, at control (24℃) and elevated (27℃) temperature conditions. Under elevated temperature, we observed: (a) a significant decrease in total parental care activity, underpinned by a decreased in male and simultaneous parental care activity, in the late stage of embryonic development; and (b) a reduction in hatching success that was however not significantly related to changes in parental care activity levels. These findings, along with the observed unaltered somatic growth of parents and decreased brood size, suggest that potential cost‐benefit trade‐offs between offspring survival (i.e., immediate fitness) and parents'' somatic condition (i.e., longer‐term fitness potential) may occur under ongoing ocean warming. Finally, our results suggest that plasticity in parental care behavior is a mechanism able to partially mitigate the negative effects of temperature‐dependent impacts.     

Evaluating new technology for biodiversity monitoring: Are drone surveys biased?

Drones and machine learning‐based automated detection methods are being used by ecologists to conduct wildlife surveys with increasing frequency. When traditional survey methods have been evaluated, a range of factors have been found to influence detection probabilities, including individual differences among conspecific animals, which can thus introduce biases into survey counts. There has been no such evaluation of drone‐based surveys using automated detection in a natural setting. This is important to establish since any biases in counts made using these methods will need to be accounted for, to provide accurate data and improve decision‐making for threatened species. In this study, a rare opportunity to survey a ground‐truthed, individually marked population of 48 koalas in their natural habitat allowed for direct comparison of the factors impacting detection probability in both ground observation and drone surveys with manual and automated detection. We found that sex and host tree preferences impacted detection in ground surveys and in manual analysis of drone imagery with female koalas likely to be under‐represented, and koalas higher in taller trees detected less frequently when present. Tree species composition of a forest stand also impacted on detections. In contrast, none of these factors impacted on automated detection. This suggests that the combination of drone‐captured imagery and machine learning does not suffer from the same biases that affect conventional ground surveys. This provides further evidence that drones and machine learning are promising tools for gathering reliable detection data to better inform the management of threatened populations.     

Considering sampling bias in close‐kin mark–recapture abundance estimates of Atlantic salmon

Genetic methods for the estimation of population size can be powerful alternatives to conventional methods. Close‐kin mark–recapture (CKMR) is based on the principles of conventional mark–recapture, but instead of being physically marked, individuals are marked through their close kin. The aim of this study was to evaluate the potential of CKMR for the estimation of spawner abundance in Atlantic salmon and how age, sex, spatial, and temporal sampling bias may affect CKMR estimates. Spawner abundance in a wild population was estimated from genetic samples of adults returning in 2018 and of their potential offspring collected in 2019. Adult samples were obtained in two ways. First, adults were sampled and released alive in the breeding habitat during spawning surveys. Second, genetic samples were collected from out‐migrating smolts PIT‐tagged in 2017 and registered when returning as adults in 2018. CKMR estimates based on adult samples collected during spawning surveys were somewhat higher than conventional counts. Uncertainty was small (CV < 0.15), due to the detection of a high number of parent–offspring pairs. Sampling of adults was age‐ and size‐biased and correction for those biases resulted in moderate changes in the CKMR estimate. Juvenile dispersal was limited, but spatially balanced sampling of adults rendered CKMR estimates robust to spatially biased sampling of juveniles. CKMR estimates based on returning PIT‐tagged adults were approximately twice as high as estimates based on samples collected during spawning surveys. We suggest that estimates based on PIT‐tagged fish reflect the total abundance of adults entering the river, while estimates based on samples collected during spawning surveys reflect the abundance of adults present in the breeding habitat at the time of spawning. Our study showed that CKMR can be used to estimate spawner abundance in Atlantic salmon, with a moderate sampling effort, but a carefully designed sampling regime is required.     

Age‐specific habitat preference,carrying capacity,and landscape structure determine the response of population spatial variability to fishing‐driven age truncation

1. Understanding the mechanisms underlying spatial variability of exploited fish is critical for the sustainable management of fish stocks. Empirical studies suggest that size‐selective fishing can elevate fish population spatial variability (i.e., more heterogeneous distribution) through age truncation, making the population less resilient to changing environment. However, species differ in how their spatial variability responds to age truncation and the underlying mechanisms remain unclear.
2. We hypothesize that age‐specific habitat preference, together with environmental carrying capacity and landscape structure, determines the response of population spatial variability to fishing‐induced age truncation. To test these hypotheses, we design an individual‐based model of an age‐structured fish population on a two‐dimensional landscape under size‐selective fishing. Individual fish reproduces and survives, and moves between habitats according to age‐specific habitat preference and density‐dependent habitat selection.
3. Population spatial variability elevates with increasing age truncation, and the response is stronger for populations with stronger age‐specific habitat preference. On a gradient landscape, reducing carrying capacity elevates the relative importance of density dependence in habitat selection, which weakens the response of spatial variability to age truncation for populations with strong age‐specific habitat preference. On a fragmented landscape, both populations with strong and weak age‐specific habitat preferences are restricted at local optimal habitats, and reducing carrying capacity weakens the responses of spatial variability to age truncation for both populations.
4. Synthesis and applications. We demonstrate that to track and predict the changes in population spatial variability under exploitation, it is essential to consider the interactive effects of age‐specific habitat preference, carrying capacity, and landscape structure. To improve spatial management in fisheries, it is crucial to enhance empirical and theoretical developments in the methodology to quantify age‐specific habitat preference of marine fish, and to understand how climatic change influences carrying capacity and landscape continuity.

     

Autumn larval cold tolerance does not predict the northern range limit of a widespread butterfly species

Climate change is driving range shifts, and a lack of cold tolerance is hypothesized to constrain insect range expansion at poleward latitudes. However, few, if any, studies have tested this hypothesis during autumn when organisms are subjected to sporadic low‐temperature exposure but may not have become cold‐tolerant yet. In this study, we integrated organismal thermal tolerance measures into species distribution models for larvae of the Giant Swallowtail butterfly, Papilio cresphontes (Lepidoptera: Papilionidae), living at the northern edge of its actively expanding range. Cold hardiness of field‐collected larvae was determined using three common metrics of cold‐induced physiological thresholds: the supercooling point, critical thermal minimum, and survival following cold exposure. P. cresphontes larvae were determined to be tolerant of chilling but generally die at temperatures below their SCP, suggesting they are chill‐tolerant or modestly freeze‐avoidant. Using this information, we examined the importance of low temperatures at a broad scale, by comparing species distribution models of P. cresphontes based only on environmental data derived from other sources to models that also included the cold tolerance parameters generated experimentally. Our modeling revealed that growing degree‐days and precipitation best predicted the distribution of P. cresphontes, while the cold tolerance variables did not explain much variation in habitat suitability. As such, the modeling results were consistent with our experimental results: Low temperatures in autumn are unlikely to limit the distribution of P. cresphontes. Understanding the factors that limit species distributions is key to predicting how climate change will drive species range shifts.     

Spatial,Temporal, and Habitat-Related Variation in Abundance of Pelagic Fishes in the Gulf of Mexico: Potential Implications of the Deepwater Horizon Oil Spill

Time-series data collected over a four-year period were used to characterize patterns of abundance for pelagic fishes in the northern Gulf of Mexico (GoM) before (2007–2009) and after (2010) the Deepwater Horizon oil spill. Four numerically dominant pelagic species (blackfin tuna, blue marlin, dolphinfish, and sailfish) were included in our assessment, and larval density of each species was lower in 2010 than any of the three years prior to the oil spill, although larval abundance in 2010 was often statistically similar to other years surveyed. To assess potential overlap between suitable habitat of pelagic fish larvae and surface oil, generalized additive models (GAMs) were developed to evaluate the influence of ocean conditions on the abundance of larvae from 2007–2009. Explanatory variables from GAMs were then linked to environmental data from 2010 to predict the probability of occurrence for each species. The spatial extent of surface oil overlapped with early life habitat of each species, possibly indicating that the availability of high quality habitat was affected by the DH oil spill. Shifts in the distribution of spawning adults is another factor known to influence the abundance of larvae, and the spatial occurrence of a model pelagic predator (blue marlin) was characterized over the same four-year period using electronic tags. The spatial extent of oil coincided with areas used by adult blue marlin from 2007–2009, and the occurrence of blue marlin in areas impacted by the DH oil spill was lower in 2010 relative to pre-spill years.     

Spatial structure of reproductive success infers mechanisms of ungulate invasion in Nearctic boreal landscapes

1. Landscape change is a key driver of biodiversity declines due to habitat loss and fragmentation, but spatially shifting resources can also facilitate range expansion and invasion. Invasive populations are reproductively successful, and landscape change may buoy this success.
2. We show how modeling the spatial structure of reproductive success can elucidate the mechanisms of range shifts and sustained invasions for mammalian species with attendant young. We use an example of white‐tailed deer (deer; Odocoileus virginianus) expansion in the Nearctic boreal forest, a North American phenomenon implicated in severe declines of threatened woodland caribou (Rangifer tarandus).
3. We hypothesized that deer reproductive success is linked to forage subsidies provided by extensive landscape change via resource extraction. We measured deer occurrence using data from 62 camera traps in northern Alberta, Canada, over three years. We weighed support for multiple competing hypotheses about deer reproductive success using multistate occupancy models and generalized linear models in an AIC‐based model selection framework.
4. Spatial patterns of reproductive success were best explained by features associated with petroleum exploration and extraction, which offer early‐seral vegetation resource subsidies. Effect sizes of anthropogenic features eclipsed natural heterogeneity by two orders of magnitude. We conclude that anthropogenic early‐seral forage subsidies support high springtime reproductive success, mitigating or exceeding winter losses, maintaining populations.
5. Synthesis and Applications. Modeling spatial structuring in reproductive success can become a key goal of remote camera‐based global networks, yielding ecological insights into mechanisms of invasion and range shifts to inform effective decision‐making for global biodiversity conservation.

     

Modeling the Distribution of Geodia Sponges and Sponge Grounds in the Northwest Atlantic

Deep-sea sponge grounds provide structurally complex habitat for fish and invertebrates and enhance local biodiversity. They are also vulnerable to bottom-contact fisheries and prime candidates for Vulnerable Marine Ecosystem designation and related conservation action. This study uses species distribution modeling, based on presence and absence observations of Geodia spp. and sponge grounds derived from research trawl catches, as well as spatially continuous data on the physical and biological ocean environment derived from satellite data and oceanographic models, to model the distribution of Geodia sponges and sponge grounds in the Northwest Atlantic. Most models produce excellent fits with validation data although fits are reduced when models are extrapolated to new areas, especially when oceanographic regimes differ between areas. Depth and minimum bottom salinity were important predictors in most models, and a Geodia spp. minimum bottom salinity tolerance threshold in the 34.3-34.8 psu range was hypothesized on the basis of model structure. The models indicated two currently unsampled regions within the study area, the deeper parts of Baffin Bay and the Newfoundland and Labrador slopes, where future sponge grounds are most likely to be found.     

Straight Line Foraging in Yellow-Eyed Penguins: New Insights into Cascading Fisheries Effects and Orientation Capabilities of Marine Predators

Free-ranging marine predators rarely search for prey along straight lines because dynamic ocean processes usually require complex search strategies. If linear movement patterns occur they are usually associated with travelling events or migratory behaviour. However, recent fine scale tracking of flying seabirds has revealed straight-line movements while birds followed fishing vessels. Unlike flying seabirds, penguins are not known to target and follow fishing vessels. Yet yellow-eyed penguins from New Zealand often exhibit directed movement patterns while searching for prey at the seafloor, a behaviour that seems to contradict common movement ecology theories. While deploying GPS dive loggers on yellow-eyed penguins from the Otago Peninsula we found that the birds frequently followed straight lines for several kilometres with little horizontal deviation. In several cases individuals swam up and down the same line, while some of the lines were followed by more than one individual. Using a remote operated vehicle (ROV) we found a highly visible furrow on the seafloor most likely caused by an otter board of a demersal fish trawl, which ran in a straight line exactly matching the trajectory of a recent line identified from penguin tracks. We noted high abundances of benthic scavengers associated with fisheries-related bottom disturbance. While our data demonstrate the acute way-finding capabilities of benthic foraging yellow-eyed penguins, they also highlight how hidden cascading effects of coastal fisheries may alter behaviour and potentially even population dynamics of marine predators, an often overlooked fact in the examination of fisheries’ impacts.     

Incorporating capture heterogeneity in the estimation of autoregressive coefficients of animal population dynamics using capture–recapture data

Population dynamic models combine density dependence and environmental effects. Ignoring sampling uncertainty might lead to biased estimation of the strength of density dependence. This is typically addressed using state‐space model approaches, which integrate sampling error and population process estimates. Such models seldom include an explicit link between the sampling procedures and the true abundance, which is common in capture–recapture settings. However, many of the models proposed to estimate abundance in the presence of capture heterogeneity lead to incomplete likelihood functions and cannot be straightforwardly included in state‐space models. We assessed the importance of estimating sampling error explicitly by taking an intermediate approach between ignoring uncertainty in abundance estimates and fully specified state‐space models for density‐dependence estimation based on autoregressive processes. First, we estimated individual capture probabilities based on a heterogeneity model for a closed population, using a conditional multinomial likelihood, followed by a Horvitz–Thompson estimate for abundance. Second, we estimated coefficients of autoregressive models for the log abundance. Inference was performed using the methodology of integrated nested Laplace approximation (INLA). We performed an extensive simulation study to compare our approach with estimates disregarding capture history information, and using R‐package VGAM, for different parameter specifications. The methods were then applied to a real data set of gray‐sided voles Myodes rufocanus from Northern Norway. We found that density‐dependence estimation was improved when explicitly modeling sampling error in scenarios with low process variances, in which differences in coverage reached up to 8% in estimating the coefficients of the autoregressive processes. In this case, the bias also increased assuming a Poisson distribution in the observational model. For high process variances, the differences between methods were small and it appeared less important to model heterogeneity.     

Mutualism promotes site selection in a large marine planktivore

1. Mutualism is a form of symbiosis whereby both parties benefit from the relationship. An example is cleaning symbiosis, which has been observed in terrestrial and marine environments. The most recognized form of marine cleaning symbiosis is that of cleaner fishes and their clients.
2. Cleaner species set up cleaning stations on the reef, and other species seek out their services. However, it is not well understood how the presence of cleaning stations influence movements of large highly mobile species. We examined the role of cleaning stations as a driver of movement and habitat use in a mobile client species.
3. Here, we used a combination of passive acoustic telemetry and in‐water surveys to investigate cleaning station attendance by the reef manta ray Mobula alfredi. We employed a novel approach in the form of a fine‐scale acoustic receiver array set up around a known cleaning area and tagged 42 rays. Within the array, we mapped structural features, surveyed the distribution of cleaner wrasse, and observed the habitat use of the rays.
4. We found manta ray space use was significantly associated with blue‐streak cleaner wrasse Labroides dimidiatus distribution and hard coral substrate. Cleaning interactions dominated their habitat use at this site, taking precedence over other life history traits such as feeding and courtship.
5. This study has demonstrated that cleaning symbiosis is a driver for highly mobile, and otherwise pelagic, species to visit inshore reef environments. We suggest that targeted and long‐term use of specific cleaning stations reflects manta rays having a long‐term memory and cognitive map of some shallow reef environments where quality cleaning is provided. We hypothesize that animals prefer cleaning sites in proximity to productive foraging regions.

     

Climatic change drives dynamic source–sink relationships in marine species with high dispersal potential

While there is now strong evidence that many factors can shape dispersal, the mechanisms influencing connectivity patterns are species‐specific and remain largely unknown for many species with a high dispersal potential. The rock lobsters Jasus tristani and Jasus paulensis have a long pelagic larval duration (up to 20 months) and inhabit seamounts and islands in the southern Atlantic and Indian Oceans, respectively. We used a multidisciplinary approach to assess the genetic relationships between J. tristani and J. paulensis, investigate historic and contemporary gene flow, and inform fisheries management. Using 17,256 neutral single nucleotide polymorphisms we found low but significant genetic differentiation. We show that patterns of connectivity changed over time in accordance with climatic fluctuations. Historic migration estimates showed stronger connectivity from the Indian to the Atlantic Ocean (influenced by the Agulhas Leakage). In contrast, the individual‐based model coupled with contemporary migration estimates inferred from genetic data showed stronger inter‐ocean connectivity in the opposite direction from the Atlantic to the Indian Ocean driven by the Subtropical Front. We suggest that the J. tristani and J. paulensis historical distribution might have extended further north (when water temperatures were lower) resulting in larval dispersal between the ocean basis being more influenced by the Agulhas Leakage than the Subtropical Front. As water temperatures in the region increase in accordance with anthropogenic climate change, a southern shift in the distribution range of J. tristani and J. paulensis could further reduce larval transport from the Indian to the Atlantic Ocean, adding complexity to fisheries management.