Validating a biophysical dispersal model with the early life-history traits of common sole (Solea solea L.)

Larval dispersal and juvenile survival are crucial in determining variation in recruitment, stock size and adult distribution of commercially important fish. This study investigates the dispersal of early-life stages of common sole (Solea solea L.) in the southern North Sea, both empirically and through modeling. Age at different life-history events of juvenile flatfish sampled along the coasts of Belgium, the Netherlands and the United Kingdom in 2013, 2014 and 2016, was determined through the counting of daily growth rings in the otoliths. Juveniles captured between August and October were estimated to be on average 140 days old with an average pelagic larval duration of 34 days. The hatching period was estimated between early April and mid-May followed by arrival and settlement in the nurseries between May and mid-June. Growth rates were higher off the Belgian coast than in the other nursery areas, especially in 2013, possibly due to a post-settlement differentiation. Empirical pelagic larval duration and settlement distributions were compared with the Larvae&Co larval dispersal model, which combines local hydrodynamics in the North Sea with sole larval behavior. Yearly predicted and observed settlement matched partially, but the model estimated a longer pelagic phase. The observations fitted even better with the modelled average (1995–2015) distribution curves. Aberrant results for the small juvenile sole sampled along the UK coast in March 2016, led to the hypothesis of a winter disruption in the deposition of daily growth rings, potentially related to starvation and lower food availability. The similarities between measured and modelled distribution curves cross-validated both types of estimations and accredited daily ageing of juveniles as a useful method to calibrate biophysical models and to understand early-life history of fish, both important tools in support of efficient fisheries management strategies.

Amphidromous but endemic: Population connectivity of Rhinogobius gigas (Teleostei: Gobioidei)

Rhinogobius gigas is an amphidromous fish endemic to eastern Taiwan. Fishes with the diadromous behavior are expected to have a broader distribution range and higher genetic homogeneity despite that some amphidromous fishes with limited distribution are observed and R. gigas is an additional exception with a limited distribution range. Rhinogobius gigas has been documented to be retained inshore near the river plume with a short pelagic larval duration of 30–40 days, which may account for the endemism of this species. The short marine larval stage of R. gigas may imply a population genetic structure and the aim of the present study is to test whether the population genetic structure is present in R. gigas. To test the population genetic structure, fragments of mitochondrial displacement loop and cytochrome c oxidase subunit I were sequenced to provide molecular inference for genetic structure among populations. Sixty-nine haplotypes were identified among 191 R. gigas from 10 populations of eastern Taiwan and the mean haplotype and nucleotide diversities for all samples were 0.956 and 0.0024, respectively, implying a bottleneck followed by a recent population expansion further supported by Fu’s Fs (-26.6; p < 0.001) and Tajima’s D (-1.5; p = 0.037) values. The phylogenetic analysis revealed lack of genetic structure and the bush-like median joining network without commonly shared haplotypes supports the same scenario. The genetic homogeneity is probably due to the amphidromous life history providing the opportunity for passive larval transportation among the rivers through coastal currents in eastern Taiwan. The endemism to eastern Taiwan may be a consequence of complicated interactions among short pelagic larval duration, interspecific competition and coastal currents.

Networks of marine protected areas under fluctuating connectivity: The importance of within-MPA dynamics

The design of marine protected areas (MPAs) has been optimized under assumptions of spatially and temporally homogeneous larval dispersal, despite complex spatiotemporal patterns displayed by ocean currents. Here we studied the effect of dispersal variability on the effectiveness of MPA networks across scales. We adopted a nested approach integrating the dynamics of both within and among MPA connectivity into a stochastic metapopulation model and first derived metapopulation persistence (required reproductive effort) and stability over MPA networks by partitioning within and among MPA contributions in relation to the spatial resolution of within-MPA connectivity. We applied this framework over a range of dispersal traits (spawning time and pelagic larval duration) and MPA network configurations, based on simulated biophysical connectivity along the northeast Pacific coast. Our results show how within-MPA dynamics affect predictions based on parameters of MPA networks such as MPA size, spacing, and pelagic larval duration. Increasing within-MPA spatial resolution predicted increasing population persistence and stability independently of other network properties. High-resolution within-MPA dynamics also predicted a negative relationship between species persistence and MPA spacing while that relationship was non-monotonic under low-resolution within-MPA dynamics. Our analysis also resolved the role of pelagic larval duration for scaling up within-MPA dynamics to MPA networks: species with short larval duration led to increasing network stability with MPA spacing while the opposite was observed for species with long larval duration. Our study stresses the importance of integrating fluctuating larval connectivity, both within and among MPAs, and more specifically suggest the benefit of small and nearby MPAs under increasing ocean variability.

Facultative amphidromy and pelagic larval duration plasticity of Rhinogobius formosanus (Teleostei, Gobioidei)

Rhinogobius formosanus Oshima, 1919 has long been considered an amphidromous goby. However, a landlocked population recently found in the Jingualiao Creek upstream of the Feitsui Reservoir in Taipei suggests that R. formosanus may complete its life in the river. This study aims to verify the habitat use of the landlocked population of R. formosanus collected from the Feitsui Reservoir and an amphidromous population collected in Malian Creek using otolith Sr:Ca ratio analysis. The hypothesis that early life history varies between the landlocked and migratory gobies was also tested. Genetic analyses show that the Feitsui Reservoir and Malian Creek populations are not genetically different. Rhinogobius formosanus from Malian Creek showed high-to-low otolith Sr:Ca ratios suggesting that these specimens spent a planktonic larval stage in the sea followed by a freshwater life at later stages. In contrast, R. formosanus from the Feitsui Reservoir showed constant lower otolith Sr:Ca ratios, implying a landlocked life history of fish in the creek upstream of the reservoir. In addition, the analysis of growth increments showed a longer pelagic larval duration for the fish in the Malian Creek (58.8 days) than those in the Feitsui Reservoir (38.8). Variation of pelagic larval duration in two genetically homogenous populations implies acclimatization to the reservoir by the landlocked gobies. This study shows that R. formosanus, like some other congeners, is capable of adapting to a freshwater landlocked environment in its early developmental stage and supports the hypothesis that landlocked populations may have a shorter pelagic larval duration.

Pelagic larval duration, growth rate, and population genetic structure of the tidepool snake morayUropterygius micropterusaround the southern Ryukyu Islands, Taiwan, and the central Philippines

The relationships between pelagic larval duration (PLD) and geographic distribution patterns or population genetic structures of fishes remain obscure and highly variable among species. To further understand the early life history of the tidepool snake morayUropterygius micropterusand the potential relationship between PLD and population genetic structure of this species, otolith microstructure and population genetics based on concatenated mtDNA sequence (cytochromeband cytochrome oxidase subunit I, 1,336 bp) were analyzed for 195 specimens collected from eight locations around the southern Ryukyu Islands, Taiwan, and the central Philippines. Eels with longer PLD and lower otolith growth rates were observed at relatively higher latitudes with lower water temperatures (54.6 ± 7.7 days and 1.28 ± 0.16 µm day−1on Ishigaki Island, Japan, vs. 43.9 ± 4.9 days and 1.60 ± 0.19 µm day−1on Badian, the Philippines), suggesting that leptocephali grew faster and had shortened pelagic periods in warmer waters. Meanwhile, the eels along the southwest coast of Taiwan had relatively longer PLD (57.9 ± 10.5 days), which might be associated with the more complex ocean current systems compared to their counterparts collected along the east coast of Taiwan (52.6 ± 8.0 days). However, the southwestern and eastern Taiwan groups had similar otolith growth rates (1.33 ± 0.19 µm day−1vs. 1.36 ± 0.16 µm day−1). Despite the intergroup variation in PLD, genetic analysis revealed fluent gene flow among the tidepool snake morays in the study regions, implying that intraspecies PLD variation had a weak effect on genetic structure. The leptocephalus stage might have ensured the widespread gene flow among the study areas and leptocephalus growth was likely influenced by regional water temperature.

Sensitivity and robustness of larval connectivity diagnostics obtained from Lagrangian Flow Networks

Lagrangian Flow Network (LFN) is a modelling framework in which ocean sub-areas are represented as nodes in a network interconnected by links representing transport of propagules (eggs and larvae) by currents. We asses the sensitivity and robustness of four LFN-derived connectivity metrics measuring retention and exchange. The most relevant parameters are tested over large ranges and a wide region with contrasting hydrodynamics: density of released particles, node size (spatial scale of discretization), Pelagic Larval Duration (PLD) and spawning modality. We find a minimum density of released particles that guarantees reliable values for most of the metrics examined. We also find that node size has a nontrivial influence on them. Connectivity estimates for long PLDs are more robust against biological uncertainties (PLD and spawning date) than for short PLDs. For mass-spawners releasing propagules over short periods (≈ 2-10 days), daily release must be simulated to properly consider connectivity fluctuations due to variable currents. In contrast, average connectivity estimates for species that spawn repeatedly over longer durations (few weeks to few months) remain robust even using longer periodicity (5-10 days). Our results have implications to design connectivity experiments with particle-tracking models and to evaluate the reliability of their results.