Résumé: Understanding larval connectivity patterns is critical for marine spatial planning, particularly for designing marine protected areas and managing fisheries. Patterns of larval dispersal and connectivity can be inferred from numerical transport models at large spatial and temporal scales. We assess model-based connectivity patterns between seamounts of the Southwestern Indian Ocean (SWIO) and the coastal ecosystems of Mauritius, La Reunion, Madagascar, Mozambique and South Africa, with emphasis on three shallow seamounts (La Pemuse [LP], MAD-Ridge [MR] and Walters Shoal [WS]). Using drifter trajectory and a Lagrangian model of ichthyo-plankton dispersal, we show that larvae can undertake very long dispersion, with larval distances increasing with pelagic larval duration (PLD). There are three groups of greater connectivity: the region between the eastern coast of Madagascar, Mauritius and La Reunion islands; the seamounts of the South West Indian Ridge; and the pair formed by WS and a nearby un-named seamount. Connectivity between these three groups is evident only for the longest PLD examined (360 d). Connectivity from seamounts to coastal ecosystems is weak, with a maximum of 2% of larvae originating from seamounts reaching coastal ecosystems. Local retention at the three focal seamounts (LP, MR and WS) peaks at about 11% for the shortest PLD considered (15 d) at the most retentive seamount (WS) and decreases sharply with increasing PLD. Information on PLD and age of larvae collected at MR and LP are used to assess their putative origin. These larvae are likely self-recruits but it is also plausible that they immigrate from nearby coastal sites, i.e. the southern coast of Madagascar for MR and the islands of La Reunion and Mauritius for LP.

Résumé: Understanding how ocean currents impact the distribution and connectivity of marine species, provides vital information for the effective conservation management of migratory marine animals. Here, we used a combination of molecular genetics and ocean drift simulations to investigate the spatial ecology of juvenile green turtle (Chelonia mydas) developmental habitats, and assess the role of ocean currents in driving the dispersal of green turtle hatchlings. We analyzed mitochondrial (mt)DNA sequenced from 358 juvenile green turtles, and from eight developmental areas located throughout the Southwest Indian Ocean (SWIO). A mixed stock analysis (MSA) was applied to estimate the level of connectivity between developmental sites and published genetic data from 38 known genetic stocks. The MSA showed that the juvenile turtles at all sites originated almost exclusively from the three known SWIO stocks, with a clear shift in stock contributions between sites in the South and Central Areas. The results from the genetic analysis could largely be explained by regional current patterns, as shown by the results of passive numerical drift simulations linking breeding sites to developmental areas utilized by juvenile green turtles. Integrating genetic and oceanographic data helps researchers to better understand how marine species interact with ocean currents at different stages of their lifecycle, and provides the scientific basis for effective conservation management.

Résumé: This work is part of the MADRidge Project special issue which aims to describe pelagic ecosystems in the vicinity of three prominent shallow seamounts in the South West Indian Ocean: one here named MAD-Ridge (240 m below the surface) plus Walters Shoal (18 m) on the Madagascar Ridge, and La Perouse (60 m) on the abyssal plain east of Madagascar. The three span latitudes 20 degrees S and 33 degrees S, some 1500 km. The study provides the background oceanography for the once-off, multidisciplinary snapshot cruise studies around the seamounts. As life on seamounts is determined by factors such as summit depth, proximity to the light layers of the ocean, and the ambient circulation, a first description of regional spatial-field climatologies (16-22 years) and monthly along-ridge gradients of surface wind (driving force), water column properties of sea surface temperature, mixed layer depth, chlorophyll-a and eddy kinetic energy, plus ocean currents is provided. Being relevant to many applications in the study domain, these properties in particular reveal contrasting environments along the Madagascar Ridge and between the three seamounts that should drive biological differences. Relative to the other two seamounts, MAD-Ridge is in the more extreme situation, being at the end of the East Madagascar Current, where it experiences sturdy, albeit variable, currents and the frequent passing of mesoscale eddies.

Résumé: The blue shark (Prionace glauca) is the most frequently captured shark in pelagic oceanic fisheries, especially pelagic longlines targeting swordfish and/or tunas. As part of cooperative scientific efforts for fisheries and biological data collection, information from fishery observers, scientific projects and surveys, and from recreational fisheries from several nations in the Atlantic and Indian Oceans was compiled. Data sets included information on location, size and sex, in a total of 478,220 blue shark records collected between 1966 and 2014. Sizes ranged from 36 to 394cm fork length. Considerable variability was observed in the size distribution by region and season in both oceans. Larger blue sharks tend to occur in equatorial and tropical regions, and smaller specimens in higher latitudes in temperate waters. Differences in sex ratios were also detected spatially and seasonally. Nursery areas in the Atlantic seem to occur in the temperate south-east off South Africa and Namibia, in the south-west off southern Brazil and Uruguay, and in the north-east off the Iberian Peninsula and the Azores. Parturition may occur in the tropical north-east off West Africa. In the Indian Ocean, nursery areas also seem to occur in temperate waters, especially in the south-west Indian Ocean off South Africa, and in the south-east off south-western Australia. The distributional patterns presented in this study provide a better understanding of how blue sharks segregate by size and sex, spatially and temporally, and improve the scientific advice to help adopt more informed and efficient management and conservation measures for this cosmopolitan species.

Résumé: We looked at how three ecological tracers may influence the characterization and interpretation of trophic structures in a tropical marine system, with a view to informing tracer(s) selection in future trophic ecology studies. We compared the trophic structures described by stable isotope compositions (carbon and nitrogen), the total mercury concentration (THg) and levels of essential fatty acids (EFA) at both the individual and species level. Analyses were undertaken on muscle tissue samples from fish and crustacean species caught in the waters surrounding the Seychelles. The carbon isotope composition (delta C-13) correlated to the proportion of arachidonic acid (ARA), whereas the nitrogen isotope composition (delta N-15) correlated to the proportion of docosahexaenoic acid (DHA) and THg. At the individual level, trophic position obtained with these three last tracers are similar. In ' contrast, the eicosapentaenoic acid (EPA) was not clearly correlated to any of the tracers. At the species level, the use of EFA (ARA and DHA), as compared to stable isotopes, resulted in slight structural modifications, mainly in the middle trophic levels. For example, the EFA overestimated the trophic positions of Thunnus alalunga and Etelis coruscans but underestimated the trophic positions of other snappers and groupers. While ARA mainly originates from coastal/benthic areas, DHA is conserved throughout the food web and may be used as a proxy indicator of trophic position. However, metabolic disparities can affect ecological tracers and in turn, distort the trophic structures derived from their results. This is especially true for species with close trophic ecologies. Despite these caveats, we think that analysing at the individual level the wealth of ARA, DHA and THg data that has already been obtained through earlier nutrition or food security studies would enhance our understanding of trophic structures.

Résumé: The reproductive biology of albacore tuna, Thunnus alalunga, in the western Indian Ocean was examined through analysis of the sex ratio, spawning season, length-at-maturity (L50), spawning frequency and fecundity. From 2013 to 2015, a total of 923 female and 867 male albacore were sampled. A bias in sex ratio was found in favor of females with fork length (LF) < 100 cm. Using histological analyses and gonadosomatic index, spawning was found to occur between 10°S and 30°S, mainly to the east of Madagascar from October to January. Large females contributed more to reproduction through their longer spawning period compared to small individuals. The L50 (mean ± standard error) of female albacore was estimated at 85.3 ± 0.7 cm LF. Albacore spawn on average every 2.2 days within the spawning region and spawning months, from November to January. Batch fecundity ranged between 0.26 and 2.09 million oocytes and the relative batch fecundity (mean ± standard deviation) was estimated at 53.4 ± 23.2 oocytes g-1 of somatic-gutted weight. The study provides new information on the reproductive development and classification of albacore in the western Indian Ocean. The reproductive parameters will reduce uncertainty in current stock assessment models which will eventually assist the fishery to be sustainable for future generations.

Résumé: The lipid composition of somatic and reproductive tissues was determined for female skipjack tuna Katsuwonus pelamis caught in the western Indian Ocean between latitude 10 degrees N and 20 degrees S and longitude 40 degrees and 70 degrees E. The highest total lipid (TL) contents were in the liver and gonads, with white muscle levels approximately three-fold lower. Three lipid classes dominated: triacylglycerols (TAG), sterol esters and wax esters (SE-WE) and phospholipids (PL). Collectively, these accounted for between 70 and 80% of TLs. Changes in lipid concentrations were evaluated over the maturation cycle. Immature fish had the lowest gonad and liver TL levels; concentrations of TL, TAG, SE-WE and PL accumulated from immature to mature (spawning-capable) phase, reflecting sustained vitellogenic activity of the liver and a transfer of lipids to developing oocytes from the onset of vitellogenesis. Gonado-somatic and hepato-somatic indices were positively correlated with each other and positively related to TL in the gonads and liver. Fulton's condition index and lipid concentrations in muscle did not vary significantly over the maturation cycle; fat content in the main storage tissues was undepleted as the ovary developed. Hence, K. pelamis apparently supports reproduction directly from food intake over the breeding season. In the gonads, reserve lipids (SE-WE and TAG) and sterols were related to batch fecundity but this was not the case for somatic and hepatic tissues. These results suggest that K. pelamis utilizes an income breeding strategy.

Résumé: The Indian Ocean Tuna Tagging Program provided a unique opportunity to collect demographic data on the key commercially targeted tropical tuna species in the Indian Ocean. In this paper, we focused on estimating growth rates for one of these species, yellowfin (Thunnus albacares). Whilst most growth studies only draw on one data source, in this study we use a range of data sources: individual growth rates derived from yellowfin that were tagged and recaptured, direct age estimates obtained through otolith readings, and length-frequency data collected from the purse seine fishery between 2000 and 2010. To combine these data sources, we used an integrated Bayesian model that allowed us to account for the process and measurement errors associated with each data set. Our results indicate that the gradual addition of each data type improved the model's parameter estimations. The Bayesian framework was useful, as it allowed us to account for uncertainties associated with age estimates and to provide additional information on some parameters (e.g., asymptotic length). Our results support the existence of a complex growth pattern for Indian Ocean yellowfin, with two distinct growth phases between the immature and mature life stages. Such complex growth patterns, however, require additional information on absolute age of fish and transition rates between growth stanzas. This type of information is not available from the data. We suggest that bioenergetic models may address this current data gap. This modeling approach explicitly considers the allocation of metabolic energy in tuna and may offer a way to understand the underlying mechanisms that drive the observed growth patterns.

Résumé: Diversity in the by-catch communities from the pelagic ecosystem in the tropical tuna purse seine fishery has been poorly studied. This study uses different biodiversity measures to compare drifting fish aggregating devices (FADs) and Free School sets (sets made on schools of tuna) of the Western Indian Ocean. Data was collected from observer programs carried out by the European Union between 2003 and 2010 on board Spanish and French fleets. Alpha (species diversity of a particular area) and Beta diversity (difference in species composition between different areas) was analyzed to assess differences in the number of species, abundances and the species composition between areas and fishing modes. Generalized additive models were undertaken to explore which geographical/environmental variables explain the distribution of species richness index and Shannon diversity index in both fishing modes. Results showed that by-catch species in FAD communities may be used as observatories of surface pelagic biodiversity in combination with Free School communities. FAD communities were more diverse with higher number of species (74 species) and evenly distributed than Free School communities (56 species). However, environmental variables played a more important role in Free School communities. Somalia area and Mozambique Channel were the areas with highest biodiversity rates in both fishing modes. This work contributed for the future implementation of the EAFM to manage the pelagic ecosystem in a holistic and more integrated way.

Résumé: In order to optimize frontal detection in sea surface temperature fields at 4 km resolution, a combined statistical and expert-based approach is applied to test different spatial smoothing of the data prior to the detection process. Fronts are usually detected at 1 km resolution using the histogram-based, single image edge detection (SIED) algorithm developed by Cayula and Cornillon in 1992, with a standard preliminary smoothing using a median filter and a 3 × 3 pixel kernel. Here, detections are performed in three study regions (off Morocco, the Mozambique Channel, and north-western Australia) and across the Indian Ocean basin using the combination of multiple windows (CMW) method developed by Nieto, Demarcq and McClatchie in 2012 which improves on the original Cayula and Cornillon algorithm. Detections at 4 km and 1 km of resolution are compared. Fronts are divided in two intensity classes (“weak” and “strong”) according to their thermal gradient. A preliminary smoothing is applied prior to the detection using different convolutions: three type of filters (median, average and Gaussian) combined with four kernel sizes (3 × 3, 5 × 5, 7 × 7, and 9 × 9 pixels) and three detection window sizes (16 × 16, 24 × 24 and 32 × 32 pixels) to test the effect of these smoothing combinations on reducing the background noise of the data and therefore on improving the frontal detection. The performance of the combinations on 4 km data are evaluated using two criteria: detection efficiency and front length. We find that the optimal combination of preliminary smoothing parameters in enhancing detection efficiency and preserving front length includes a median filter, a 16 × 16 pixel window size, and a 5 × 5 pixel kernel for strong fronts and a 7 × 7 pixel kernel for weak fronts. Results show an improvement in detection performance (from largest to smallest window size) of 71% for strong fronts and 120% for weak fronts. Despite the small window used (16 × 16 pixels), the length of the fronts has been preserved relative to that found with 1 km data. This optimal preliminary smoothing and the CMW detection algorithm on 4 km sea surface temperature data are then used to describe the spatial distribution of the monthly frequencies of occurrence for both strong and weak fronts across the Indian Ocean basin. In general strong fronts are observed in coastal areas whereas weak fronts, with some seasonal exceptions, are mainly located in the open ocean. This study shows that adequate noise reduction done by a preliminary smoothing of the data considerably improves the frontal detection efficiency as well as the global quality of the results. Consequently, the use of 4 km data enables frontal detections similar to 1 km data (using a standard median 3 × 3 convolution) in terms of detectability, length and location. This method, using 4 km data is easily applicable to large regions or at the global scale with far less constraints of data manipulation and processing time relative to 1 km data.

Résumé: Climate-induced changes in the physical, chemical, and biological environment are expected to increasingly stress marine ecosystems, with important consequences for fisheries exploitation. Here, we use the APECOSM-E numerical model (Apex Predator ECOSystem Model – Estimation) to evaluate the future impacts of climate change on the physiology, spatial distribution, and abundance of skipjack tuna, the worldwide most fished species of tropical tuna. The main novelties of our approach lie in the mechanistic link between environmental factors, metabolic rates, and behavioral responses and in the fully three dimensional representation of habitat and population abundance. Physical and biogeochemical fields used to force the model are provided by the last generation of the IPSL-CM5 Earth System Model run from 1990 to 2100 under a &8216;business-as-usual&8217; scenario (RCP8.5). Our simulations show significant changes in the spatial distribution of skipjack tuna suitable habitat, as well as in their population abundance. The model projects deterioration of skipjack habitat in most tropical waters and an improvement of habitat at higher latitudes. The primary driver of habitat changes is ocean warming, followed by food density changes. Our projections show an increase of global skipjack biomass between 2010 and 2050 followed by a marked decrease between 2050 and 2095. Spawning rates are consistent with population trends, showing that spawning depends primarily on the adult biomass. On the other hand, growth rates display very smooth temporal changes, suggesting that the ability of skipjack to keep high metabolic rates in the changing environment is generally effective. Uncertainties related to our model spatial resolution, to the lack or simplification of key processes and to the climate forcings are discussed.

Résumé: Effective use of spatial management in the pelagic realm presents special challenges due to high fish and fisher mobility, limited knowledge and significant governance challenges. The tropical Indian Ocean provides an ideal case study for testing our ability to apply existing data sources to assessing impacts of spatial management on tuna fisheries because of several recent controversial spatial closures. We review the scientific underpinnings of pelagic MPA effects, spatio-temporal patterns of Indian Ocean tuna catch, by catch and fish movements, and the consequences of these for the efficacy of spatial management for Indian Ocean tropical tuna fisheries. The tropical Indian Ocean is characterized by strong environmental fluctuations, regular seasonal variability in catch, large observed tuna displacement distances, relatively uniform catch-per-unit-effort and bycatch rates over space, and high fisher mobility, all of which suggest significant variability and movement in tropical tuna fisheries that are simply not well adapted to static spatial closures. One possible exception to this overall conclusion would be a large time/area closure east of Somalia. If closed for a significant fraction of the year it could reduce purse-seine bycatch and juvenile tuna catch. Dynamic closures following fish migratory patterns are possible, but more focused information on fish movements will be needed for effective implementation. Fortunately, several recent improvements in conventional fishery management and reporting will likely enhance our ability to evaluate spatial and non-spatial management options in the near future, particularly as pertaining to bycatch species.

Résumé: We analyzed stomach contents and muscle isotopic composition (delta C-13, delta N-15) of yellowstripe goatfish Mulloidichthys flavolineatus from a coral reef to investigate size-related dietary changes and resource or habitat partitioning. Juveniles (< 12 cm total length [TL]), young adults (12 <= TL < 17 cm) and adults (>= 17 cm TL) showed a high diet overlap, especially between juveniles and young adults. According to stomach contents analysis, M. flavolineatus widens its prey spectrum with increasing size from a common prey pool that includes polychaetes, tanaids and harpacticoid copepods. We observed a significant increase in delta C-13 values (from -17.1 +/- 0.5 parts per thousand for juveniles to -10.7 +/- 0.8 parts per thousand for adults), which were correlated to fish size. Adults (delta N-15 mean of 11.1 +/- 1.8 parts per thousand) were one trophic level above juveniles and young adults (7.7 +/- 0.5 and 7.4 +/- 0.5 parts per thousand, respectively). These patterns of isotopic changes confirmed ontogenetic dietary shifts. However, trophodynamics can be influenced by physiological factors such as growth and sexual maturity. M flavolineatus shift from a pelagic to a macrobenthic diet, which is equilibrated at the adult stage. Results from combined stomach contents (prey volume) and stable isotope analyses suggested a dominant contribution of polychaetes (macrofauna), possibly through selective feeding. Conversely, on the basis of prey volume and stable isotope data, meiofauna did not feature significantly in the diet despite their high abundance in stomachs. Two adult groups were distinguished based on their delta N-15 values (11.9 +/- 0.8 and 7.8 +/- 0.6 parts per thousand), indicating possible stage-specific partitioning in habitat use inside the reef.