Model-derived dispersal pathways from multiple source populations explain variability of invertebrate larval supply

Connectivity modelling identifies sources and sinks of coral recruitment within reef clusters

Connectivity aids the recovery of populations following disturbances, such as coral bleaching and tropical cyclones. Coral larval connectivity is a function of physical connectivity and larval behaviour. In this study, we used OceanParcels, a particle tracking simulator, with 2D and 3D velocity outputs from a high resolution hydrodynamic-biogeochemical marine model (RECOM) to simulate the dispersal and settlement of larvae from broadcast spawning Acropora corals in the Moore Reef cluster, northern Great Barrier Reef, following the annual spawning events in 2015, 2016 and 2017. 3D velocity simulations showed 19.40-68.80% more links and sinks than those of 2D simulations. Although the patterns of connectivity among sites vary over days and years, coral larvae consistently dispersed from east to west in the cluster domain, with some sites consistently acting as sources or sinks for local larval recruitment. Results can inform coral reef intervention plans for climate change, such as the design of marine protected areas and the deployment of proposed interventions within reef clusters. For example, the wider benefits of interventions (e.g., deployment of heat adapted corals) may be optimised when deployed at locations that are a source of larvae to others within comparable habitats across the reef cluster.

Climatic resilience: Marine heatwaves do not influence the variations of green crab (Carcinus maenas) megalopae supply patterns to a Western Iberian estuary

Extreme climatic events like marine heatwaves (MHWs) are becoming more frequent, intense, and longer lasting all around the world. The consequences of these anomalously warm periods are devastating for marine ecosystems. Still, little is known about these extreme events off the western Iberia coast. Here we analyzed MHW events occurring from 1982 to 2020 on the Aveiro coast, western Iberia coast of Portugal. A total of 79 events were detected for the region, with an average duration of 15.8 days, and a mean intensity of 1.9 °C ± 0.4 °C above the 90th percentile of sea surface temperatures (SST) for the region. The maximum intensity of the events has increased by 0.5 °C over the last decade. The relation between SST, and therefore, MHW events, the North Atlantic Oscillation index (NAO), and the regional Iberian Upwelling Index (UI) was identified. The intense upwelling of the region seems to mitigate the duration of warming conditions, resulting in shorter MHW events. Furthermore, the impacts of SST and MHW events on the supply patterns of Carcinus maenas megalopae were examined, utilizing daily data from 2002, 2006-2009, 2012, and 2013, collected at the entrance of Ria de Aveiro. Cross-correlations were employed to assess the effect of SST on megalopae supply, while ordinary least square cumulative sums were used to identify variations over time. The influence of SST on supply was noticed with a 5-to-11-day lag, but this relation changed over the years. Contrary to our hypothesis, we found no evidence supporting a diminishment in megalopae supply due to MHW events. These elusive findings, coupled with the apparent lack of influence of these extreme events, highlight the relatively weak intensity and brief duration of the MHW events in the region, coupled with the high thermal tolerance of these species.

Reef larval recruitment in response to seascape dynamics in the SW Atlantic

Advances in satellite observation have improved our capacity to track changes in the ocean with numerous ecological and conservation applications, which are yet under-explored for coastal ecology. In this study, we assessed the spatio-temporal dynamics in invertebrate larval recruitment and the Seascape Pelagic Habitat Classification, a satellite remote-sensing product developed by the Marine Biodiversity Observation Network (MBON) and delivered by the US National Oceanic and Atmospheric Administration to monitor biodiversity globally. Our ultimate goal was to identify and predict changes in coastal benthic assemblages at tropical reefs in the SW Atlantic based on integrated pelagic conditions, testing the use of MBON Seascape categorization. Our results revealed that the pelagic Seascapes correlated with monthly and seasonal variations in recruitment rates and assemblage composition. Recruitment was strongly influenced by subtropical Seascapes and was reduced by the presence of warm waters with high-nutrient contents and phytoplankton blooms, which are likely to affect reef communities in the long term. Recruitment modeling indicates that Seascapes may be more efficient than sea surface temperature in predicting benthic larval dynamics. Based on historical Seascape patterns, we identified seven events that may have impacted benthic recruitment in this region during the last decades. These findings provide new insights into the application of novel satellite remote-sensing Seascape categorizations in benthic ecology and evidence how reef larval supply in the SW Atlantic could be impacted by recent and future ocean changes.

Fatty acid ratio analysis identifies changes in competent meroplanktonic larvae sampled over different supply events

Planktonic communities are a cornerstone of ocean food webs. Early benthic performance of meroplanktonic organisms is shaped by their life stages in planktonic communities. Fatty acid profiles of marine invertebrates are a good indicator of their nutritional state and allow inferring how dietary regimes experienced during larval pelagic life may drive their pre- and post-metamorphosis performance. Fatty acid profiles of Carcinus maenas megalopae were analysed during four larval supply events in two consecutive years to better understand the variability in their nutritional state at settlement. The logratio analysis of fatty acids showed differences between the four larval supply events, with five ratios explaining 83.1% of the variance. The ratios that contributed to separate larval supply events presented a combination of essential, de novo synthetized and diet origin fatty acids (e.g., phytanate/20:4 n-6, 16:0/18:2 n-4). The high fatty acid signature dispersion found within the same supply event suggests that larvae settling at Ria de Aveiro (Portugal) developed through different planktonic feeding zones and experienced contrasting feeding regimes. The fatty acid profile of megalopae demonstrated a high contribution of diatoms, flagellates and bacteria in the larval diet of C. maenas. The present study demonstrated differences between supply events, although a high variability of larval phenotypes was recorded within the same supply event.

Independent estimates of marine population connectivity are more concordant when accounting for uncertainties in larval origins

Marine larval dispersal is a complex biophysical process that depends on the effects of species biology and oceanography, leading to logistical difficulties in estimating connectivity among populations of marine animals with biphasic life cycles. To address this challenge, the application of multiple methodological approaches has been advocated, in order to increase confidence in estimates of population connectivity. However, studies seldom account for sources of uncertainty associated with each method, which undermines a direct comparative approach. In the present study we explicitly account for the statistical uncertainty in observed connectivity matrices derived from elemental chemistry of larval mussel shells, and compare these to predictions from a biophysical model of dispersal. To do this we manipulate the observed connectivity matrix by applying different confidence levels to the assignment of recruits to source populations, while concurrently modelling the intrinsic misclassification rate of larvae to known sources. We demonstrate that the correlation between the observed and modelled matrices increases as the number of observed recruits classified as unknowns approximates the observed larval misclassification rate. Using this approach, we show that unprecedented levels of concordance in connectivity estimates (r = 0.96) can be achieved, and at spatial scales (20-40 km) that are ecologically relevant.

Integrated Monitoring of Mola mola Behaviour in Space and Time

Over the last decade, ocean sunfish movements have been monitored worldwide using various satellite tracking methods. This study reports the near-real time monitoring of fine-scale (< 10 m) behaviour of sunfish. The study was conducted in southern Portugal in May 2014 and involved satellite tags and underwater and surface robotic vehicles to measure both the movements and the contextual environment of the fish. A total of four individuals were tracked using custom-made GPS satellite tags providing geolocation estimates of fine-scale resolution. These accurate positions further informed sunfish areas of restricted search (ARS), which were directly correlated to steep thermal frontal zones. Simultaneously, and for two different occasions, an Autonomous Underwater Vehicle (AUV) video-recorded the path of the tracked fish and detected buoyant particles in the water column. Importantly, the densities of these particles were also directly correlated to steep thermal gradients. Thus, both sunfish foraging behaviour (ARS) and possibly prey densities, were found to be influenced by analogous environmental conditions. In addition, the dynamic structure of the water transited by the tracked individuals was described by a Lagrangian modelling approach. The model informed the distribution of zooplankton in the region, both horizontally and in the water column, and the resultant simulated densities positively correlated with sunfish ARS behaviour estimator (r_s = 0.184, p<0.001). The model also revealed that tracked fish opportunistically displace with respect to subsurface current flow. Thus, we show how physical forcing and current structure provide a rationale for a predator’s fine-scale behaviour observed over a two weeks in May 2014.

Influence of Biological Factors on Connectivity Patterns for Concholepas concholepas (loco) in Chile

In marine benthic ecosystems, larval connectivity is a major process influencing the maintenance and distribution of invertebrate populations. Larval connectivity is a complex process to study as it is determined by several interacting factors. Here we use an individual-based, biophysical model, to disentangle the effects of such factors, namely larval vertical migration, larval growth, larval mortality, adults fecundity, and habitat availability, for the marine gastropod Concholepas concholepas (loco) in Chile. Lower transport success and higher dispersal distances are observed including larval vertical migration in the model. We find an overall decrease in larval transport success to settlement areas from northern to southern Chile. This spatial gradient results from the combination of current direction and intensity, seawater temperature, and available habitat. From our simulated connectivity patterns we then identify subpopulations of loco along the Chilean coast, which could serve as a basis for spatial management of this resource in the future.

Behaviorally Mediated Larval Transport in Upwelling Systems

Highly advective upwelling systems along the western margins of continents are widely believed to transport larvae far offshore in surface currents resulting in larval wastage, limited recruitment, and increased population connectivity. However, suites of larval behaviors effectively mediate interspecific differences in the extent of cross-shelf migrations between nearshore adult habitats and offshore larval habitats. Interspecific differences in behavior determining whether larvae complete development in estuaries or migrate to the continental shelf are evident in large estuaries, but they sometimes may be disrupted by turbulent tidal flow or the absence of a low-salinity cue in shallow, low-flow estuaries, which are widespread in upwelling systems. Larvae of most species on the continental shelf complete development in the coastal boundary layer of reduced flow, whereas other species migrate to the mid- or outer shelf depending on how much time is spent in surface currents. These migrations are maintained across latitudinal differences in the strength and persistence of upwelling, in upwelling jets at headlands, over upwelling-relaxation cycles, and among years of varying upwelling intensity. Incorporating larval behaviors into numerical models demonstrates that larvae recruit closer to home and in higher numbers than when larvae disperse passively or remain in surface currents.

Overview on the European green crab Carcinus spp. (Portunidae, Decapoda), one of the most famous marine invaders and ecotoxicological models

Green crabs (Carcinus, Portunidae) include two species native to Europe--Carcinus aestuarii (Mediterranean species) and Carcinus maenas (Atlantic species). These small shore crabs (maximal length carapace, approximately 10 cm) show rapid growth, high fecundity, and long planktonic larval stages that facilitate broad dispersion. Carcinus spp. have a high tolerance to fluctuations of environmental factors including oxygen, salinity, temperature, xenobiotic compounds, and others. Shipping of Carcinus spp. over the past centuries has resulted in its invasions of America, Asia, and Australia. Classified as one of the world's 100 worst invaders by the International Union for Conservation of Nature, Carcinus spp. are the most widely distributed intertidal crabs in the world. Their voracious predatory activity makes them strong interactors in local communities, and they are recognized as a model for invasiveness in marine systems as well as a sentinel species in ecotoxicology. This review shows an exhaustive analysis of the literature on the life cycle, diversity, physiological tolerance, genomic investigations, ecotoxicological use, historical invasion, control programs, and putative economical valorization of shore crabs.

Role of upwelling on larval dispersal and productivity of gooseneck barnacle populations in the Cantabrian Sea: management implications

The effect of coastal upwelling on the recruitment and connectivity of coastal marine populations has rarely been characterized to a level of detail to be included into sound fishery management strategies. The gooseneck barnacle (Pollicipes pollicipes) fishery at the Cantabrian Coast (Northern Spain) is located at the fringes of the NW Spanish Upwelling system. This fishery is being co-managed through a fine-scale, interspersed set of protected rocks where each rock receives a distinct level of protection. Such interspersion is potentially beneficial, but the extent to which such spacing is consistent with mean larval dispersal distances is as yet unknown. We have simulated the spread of gooseneck barnacle larvae in the Central Cantabrian Coast using a high-resolution time-series of current profiles measured at a nearshore location. During a year of high upwelling activity (2009), theoretical recruitment success was 94% with peak recruitment predicted 56 km west of the emission point. However, for a year of low upwelling activity (2011) theoretical recruitment success dropped to 15.4% and peak recruitment was expected 13 km east of the emission point. This is consistent with a positive correlation between catch rates and the Integrated Upwelling Index, using a 4-year lag to allow recruits to reach commercial size. Furthermore, a net long-term westward larval transport was estimated by means of mitochondrial cytochrome c oxidase subunit I (COI) sequences for five populations in the Cantabrian Sea. Our results call into question the role of long distance dispersal, driven by the mesoscale processes in the area, in gooseneck barnacle populations and point to the prevalent role of small-scale, asymmetric connectivity more consistent with the typical scale of the co-management process in this fishery.