Moving toward stock units identification based on spatial population structure of marine species in the Mediterranean Sea and adjacent waters through multidisciplinary and holistic approaches for the sustainability of fisheries resources

Investigating stock identity of marine species in a multidisciplinary holistic approach can reveal patterns of complex spatial population structure and signatures of potential local adaptation. The population structure of common sole (Solea solea) in the Mediterranean Sea was delineated using genomic and otolith data, including single nucleotide polymorphisms (SNPs) markers and otolith data. SNPs were correlated with environmental and spatial variables to evaluate the impact of these features on the actual genetic population structure. Integrated holistic approach was applied to combine the tracers with different spatio-temporal scales. SNPs data was also used to illustrate the population structure of European hake (Merluccius merluccius) within the Alboran Sea, extending into the neighboring Mediterranean Sea and Atlantic Ocean. The aim was to identify patterns of neutral and potential adaptive genetic variation by applying seascape genomic framework. Results from both genetic and otolith data suggested significant divergence among putative populations of common sole, confirming a clear separation between Western, Adriatic Sea and Eastern Mediterranean Sea. Evidence of fine-scale population structure in the Western Mediterranean Sea was observed at outlier loci level and in the Adriatic. Our study not only indicates that separation among Mediterranean sole population is led primarily by neutral processes, but it also suggests the presence of local adaptation influenced by environmental and spatial factors. The holistic approach by considering the spatio-temporal scales of variation confirmed that the same pattern of separation between these geographical sites is currently occurring and has occurred for many generations. Results showed the occurrence of population structure in Merluccius merluccius by detecting westward–eastward differentiation among populations and distinct subgroups at a fine geographical scale using outlier SNPs. These results enhance the knowledge of the population structure of commercially relevant species to support the application of spatial stock assessment models, including a redefinition of fishery management units.

Spatial and temporal variability and ecological processes in the epibenthic assemblages of the northern Adriatic Sea

Several coralligenous reefs occur in the soft bottoms of the northern Adriatic continental shelf. Mediterranean coralligenous habitats are characterised by high species diversity and are intrinsically valuable for their biological diversity and for the ecological processes they support. The conservation and management of these habitats require quantifying spatial and temporal variability of their benthic assemblages. This PhD thesis aims to give a relevant contribution to the knowledge of the structure and dynamics of the epibenthic assemblages on the coralligenous subtidal reefs occurring in the northern Adriatic Sea. The epibenthic assemblages showed a spatial variation larger compared to temporal changes, with a temporal persistence of reef-forming organisms. Assemblages spatial heterogeneity has been related to morphological features and geographical location of the reefs, together with variation in the hydrological conditions. Manipulative experiments help to understand the ecological processes structuring the benthic assemblages and maintaining their diversity. In this regards a short and long term experiment on colonization patterns of artificial substrata over a 3-year period has been performed in three reefs, corresponding to the three main types of assemblages detected in the previous study. The first colonisers, largely depending by the different larval supply, played a key role in determining the heterogeneity of the assemblages in the early stage of colonisation. Lateral invasion, from the surrounding assemblages, was the driver in structuring the mature assemblages. These complex colonisation dynamics explained the high heterogeneity of the assemblages dwelling on the northern Adriatic biogenic reefs. The buildup of these coralligenous reefs mainly depends by the bioconstruction-erosion processes that has been analysed through a field experiment. Bioconstruction, largely due to serpulid polychaetes, prevailed on erosion processes and occurred at similar rates in all sites. Similarly, the total energy contents in the benthic communities do not differ among sites, despite being provided by different species. Therefore, we can hypothesise that both bioconstruction processes and energetic storage may be limited by the availability of resources. Finally the major contribution of the zoobenthos compared to the phytobenthos to the total energetic content of assemblages suggests that the energy flow in these benthic habitats is primarily supported by planktonic food web trough the filter feeding invertebrates.

Diet of bluefin tuna (Thunnus thynnus) and swordfish (Xhipias gladius) in the Adriatic Sea

The present issue analyses bluefin tuna (Thynnus thunnus) and swordfish (Xiphias glaudis) diet, caught by professional long-line fishing in the middle Adriatic Sea (Pomo pit). These species represent apex predators in pelagic environment that may play key roles in determining food web structure and ecosystem dynamics. The studies about their feedings habits, based upon stomach contents analysis, are important for the comprehension of biological and ecological interaction. Over the years, many studies have been performed on the diet of tuna and swordfish in the Mediterranean Sea. This research is based on a fairly wide number of analyzed stomach contents, in comparison with the previous ones. In this work, the analysis of 340 stomach contents of bluefin tuna caught by long-line in the central Adriatic sea confirms in general the opportunistic behaviour of this species. Finding support the hypothesis that Adriatic tuna and swordfish chase their food over a wide bathymetric zone and probably near the surface at night. No indication of food preference respect to size of predator or sample season are found. It seems that the two species are able to cohabit because their trophic niche are not overlapped, changing during the time and the vertical and horizontal space.

Demersal communities in the Mediterranean Sea: a case study of Triglidae (Osteichthyes, Scorpaeniformes) on the conservation and sustainable use of marine resources

An appropriate management of fisheries resources can only be achieved with the continuous supply of information on the structure and biology of populations, in order to predict the temporal fluctuations. This study supports the importance of investigating the bio-ecology of increasingly exploited and poorly known species, such as gurnards (Osteichthyes, Triglidae) from Adriatic Sea (Mediterranean), to quantify their ecological role into marine community. It also focuses on investigate inter and intra-specific structuring factor of Adriatic population. These objectives were achieved by: 1) investigating aspects of the population dynamics; 2) studying the feeding biology through the examination of stomach contents; 3) using sagittal otoliths as potential marker of species life cycle; 4) getting preliminary data on mDNA phylogeny. Gurnards showed a specie-specific “critical size” coinciding with the start of sexual maturity, the tendency to migrate to greater depths, a change of diet from crustaceans to fish and an increase of variety of food items eaten. Distribution of prey items, predator size range and depth distribution were the main dimensions that influence the breadth of trophic niche and the relative difference amongst Adriatic gurnards. Several feeding preferences were individuated and a possible impact among bigger-size gurnards and other commercial fishes (anchovy, gadoids) and Crustacea (such as mantis prawn and shrimps) were to be necessary considered. Otolith studies showed that gurnard species have a very fast growth despite other results in other areas; intra-specific differences and the increase in the variability of otolith shape, sulcus acusticus shape, S:O ratios, sulcus acusticus external crystals arrangement were shown between juveniles and adults and were linked to growth (individual genetic factors) and to environmental conditions (e.g. depth and trophic niche distribution). In order to facilitate correct biological interpretation of data, molecular data were obtained for comparing morphological distance to genetic ones.

Toward an understanding of the atmospheric forcing uncertainties in the Mediterranean Sea

The thesis has extensively investigated for the first time the statistical distributions of atmospheric surface variables and heat fluxes for the Mediterranean Sea. After retrieving a 30-year atmospheric analysis dataset, we have captured the spatial patterns of the probability distribution of the relevant atmospheric variables for ocean atmospheric forcing: wind components (U,V), wind amplitude, air temperature (T2M), dewpoint temperature (D2M) and mean sea-level pressure (MSL-P). The study reveals that a two-parameter PDF is not a good fit for T2M, D2M, MSL-P and wind components (U,V) and a three parameter skew-normal PDF is better suited. Such distribution captures properly the data asymmetric tails (skewness). After removing the large seasonal cycle, we show the quality of the fit and the geographic structure of the PDF parameters. It is found that the PDF parameters vary between different regions, in particular the shape (connected to the asymmetric tails) and the scale (connected to the spread of the distribution) parameters cluster around two or more values, probably connected to the different dynamics that produces the surface atmospheric fields in the Mediterranean basin. Moreover, using the atmospheric variables, we have computed the air-sea heat fluxes for a 20-years period and estimated the net heat budget over the Mediterranean Sea. Interestingly, the higher resolution analysis dataset provides a negative heat budget of –3 W/m2 which is within the acceptable range for the Mediterranean Sea heat budget closure. The lower resolution atmospheric reanalysis dataset(ERA5) does not satisfy the heat budget closure problem pointing out that a minimal resolution of the atmospheric forcing is crucial for the Mediterranean Sea dynamics. The PDF framework developed in this thesis will be the basis for a future ensemble forecasting system that will use the statistical distributions to create perturbations of the atmospheric ocean forcing.