Relative Inﬂuence of Environmental Factors on Biodiversity and Behavioural Traits of a Rare Mesopelagic Fish, Trachipterus trachypterus (Gmelin, 1789), in a Continental Shelf Front of the Mediterranean Sea

: Coastal environments can be inﬂuenced by water body masses with particular physical, chemical, and biological properties that create favourable conditions for the development of unique planktonic communities. In this study, we investigated a continental shelf front at Ponza Island (Tyrrhenian Sea) and discussed its diversity and complexity in relation to major environmental parameters. Moon phase and current direction were found to play a signiﬁcant role in shaping species abundance and behaviour. During in situ observations, we also provided the ﬁrst data on the behaviour of juveniles of a rare mesopelagic species, Trachipterus trachypterus , suggesting the occurrence of Batesian mimicry.


Introduction
Open water environments are highly dynamic ecosystems [1], with ecological processes spreading over a wide range of spatial and temporal scales [2]. At a meso-and sub-mesoscale (~1-100 km), these processes emerge as dominant structuring regimes for populations of marine organisms in the short time span (days to months) [3,4], encompassing key biological dynamics, such as bloom lifetime or behavioural switches of marine predators [5]. These ecological processes are particularly relevant for drifting organisms (plankton) that are advected and dispersed by the ocean currents [6]. Within this complex scenario, sampling probabilities for uncommon and little-known species can be higher than those reported in all other marine environments, in which these species are typically rare and occur at low abundances. This may result in uncertainty about the real species abundances and their ecological traits. As a result, modern biodiversity surveillance programs often investigate the occurrence of species mostly present at moderate-to-high abundance level [7,8].

Study Area and Sampling Design
The study was performed in the framework of the AQUATILIS EXPEDITION at Ponza Island (Tyrrhenian Sea, Italy) from the 8th to the 28th of April 2018. Due to its bathymetric profile, the area is exposed to prominent surface and upwelling filaments, resulting in a complex hydrodynamism that favours the accumulation of planktonic species in patched areas, according to the prevalent weather forcing [29]. Diving sites were close to "Secca delle Formiche" (40.884947 • N, 12.978043 • E), a group of emerged rocks, located 0.8 nautical miles (nm) SE from the port of Ponza. The bank rapidly sinks from a depth of 5 m to over 500 m (Figure 1), representing a suitable place to study the upwelling of deep mesopelagic species.
Previous surveys conducted in April 2016 and 2017 highlighted a remarkable difference in plankton abundance and biodiversity in the study area according to time; in particular, observation during night dives showed a more complex plankton biodiversity, if compared with day dives (authors' personal note). As a consequence, in the present study, scuba night dives have been planned to start one hour after sunset. Two different groups of three divers each equipped with Nikon D850 with 50 mm f2 Zeiss Milvus macro lens for photography, Panasonic Lumix GH5 with Panasonic Lumix G 8 mm f3.5 fish eye and Panasonic Lumix 15 mm f1.7 lenses for video and tools (fine mesh hand nets, 250 mL plastic cans) for sampling collection, performed alternate dives. In order to increase the image quality, all the cameras were equipped with a buoyancy stabilization system. Sampling depth varied between 0.5 and 10 m, freely drifting within the surface current. To ensure taxon identity, 4 specimens of T. trachypterus were collected and analyzed for morphological traits. Behaviour was recorded by filming all the specimens detected. Observations followed the same schedule, starting one hour after sunset, freely drifting in the first 10 m of the water column. Surveys were conducted daily between 21:00 and 24:00 h. Previous surveys conducted in April 2016 and 2017 highlighted a remarkable difference in plankton abundance and biodiversity in the study area according to time; in particular, observation during night dives showed a more complex plankton biodiversity, if compared with day dives (authors' personal note). As a consequence, in the present study, scuba night dives have been planned to start one hour after sunset. Two different groups of three divers each equipped with Nikon D850 with 50mm f2 Zeiss Milvus macro lens for photography, Panasonic Lumix GH5 with Panasonic Lumix G 8 mm f3.5 fish eye and Panasonic Lumix 15 mm f1.7 lenses for video and tools (fine mesh hand nets, 250 mL plastic cans) for sampling collection, performed alternate dives. In order to increase the image quality, all the cameras were equipped with a buoyancy stabilization system. Sampling depth varied between 0.5 and 10 m, freely drifting within the surface current. To ensure taxon identity, 4 specimens of T. trachypterus were collected and analyzed for morphological traits. Behaviour was recorded by filming all the specimens detected. Observations followed the same schedule, starting one hour after sunset, freely drifting in the first 10 m of the water column. Surveys were conducted daily between 21:00 and 24:00 h.

Environmental Data Collection
Main abiotic environmental conditions were inferred from numerical models and in situ measurements (current direction). To this aim, wind records from available meteorological station and numerical results of Surface Sea Temperature (SST) and Salinity (SS) were collected considering both surveying period and long-time series. The wind records were downloaded from the Gaeta Port Authority and the Mareographic network websites (www.mareografico.it; www.portidiroma.it). Ponza's records extend only for a 4-year period (2011-2015), while more than 10-year records are available from Gaeta (2007-2018) (Supplementary Materials, Figure S1). Despite the distance between stations, wind roses can be considered comparable [29,30]. During surveys, we used Gaeta Meteo Station period to analyse wind conditions. The SST and SS datasets were downloaded from the ''Copernicus'' data portal (http://marine.copernicus.eu, accessed on the 21 st of November 2018). SST and SS from 1987 to present were analysed, producing long-term temporal mean distribution of variables during April months and mean conditions during the surveyed period.

Environmental Data Collection
Main abiotic environmental conditions were inferred from numerical models and in situ measurements (current direction). To this aim, wind records from available meteorological station and numerical results of Surface Sea Temperature (SST) and Salinity (SS) were collected considering both surveying period and long-time series. The wind records were downloaded from the Gaeta Port Authority and the Mareographic network websites (www.mareografico.it; www.portidiroma.it). Ponza's records extend only for a 4-year period (2011-2015), while more than 10-year records are available from Gaeta (2007-2018) (Supplementary Materials, Figure S1). Despite the distance between stations, wind roses can be considered comparable [29,30]. During surveys, we used Gaeta Meteo Station period to analyse wind conditions. The SST and SS datasets were downloaded from the "Copernicus" data portal (http://marine.copernicus.eu, accessed on the 21st of November 2018). SST and SS from 1987 to present were analysed, producing long-term temporal mean distribution of variables during April months and mean conditions during the surveyed period.
Moonlight is known to affect zooplankton and micro-nekton dynamics, as well as their predators in tropical, subtropical, and Arctic waters [31][32][33][34][35][36]. Moonlight also drives vertical migrations to bathypelagic depths (~1000 m) involving a cascade response, with deeper-living organisms responding to movements of shallower-living, surface-synchronized populations [37]. According to this evidence, the main macro-zooplanktonic species composition was recorded and correlated to moon phases and moonlight intensity.

Data Analysis
Non-parametric multivariate analyses were performed to assess the relation between the species observed and the environmental parameters over the survey time period. To analyze the environmental conditions (SST, moon phase, wind direction, and current direction), the parameters were normalized and standardized to carry out a Principal Component Analysis (PCA) with Statistical Package for Social Science (SPSS, Version 26.0, IBM Corp© 2019). Once the significant environmental patterns were established, standardized data of the observed species were used to investigate the variability of these patterns. To assess variations in the species diversity, data of the thirty-five taxa were analyzed using distance-based permutational multivariate analysis of variance (PERMANOVA) [38,39] with the method of permutation of the residuals under a reduced model, according to a 2-factor nested design: moon phase and current direction. The homogeneity of multivariate dispersions was tested with the PERMDISP routine, mainly to estimate dispersion at the moon phase level. Non-parametric multidimensional scaling ordination (nMDS) was used to examine the behaviour of the species grouped according to the current direction pattern. The Similarity Percentage analysis procedure (SIMPER) was used to identify the main species that characterized each current direction pattern. These analyses were performed on Primer v.7 with PERMANOVA+ software [39,40]. A PCA with supplementary variables was performed with CANOCO 5 [41] in order to investigate the species response to each variable. Finally, the hierarchical cluster analysis was applied to estimate the number of clusters using Primer v.7 [40].

Oceanographic Landscape
Available time-lapse records of Ponza Island and Gaeta showed three main directions (North-East, South-East and West-Northwest) of wind. Between the 12th and the 18th of April 2018, a moderate north-west wind dominated, while in the other days, weak northern flows (8th-11th of April) and eastern flows (19th-28th of April) were recorded (Supplementary Materials, Figure S2B). The SST mean distribution of April showed a quite uniform temperature at a regional scale, with variations of less than 1 • C (16.8-17.6 • C). Between the 9th and the 18th of April 2018, the SST was generally cooler than the mean long-term map, with a weak NW-SE increasing gradient. In correspondence of the study area, the SST reached 15.3-15.5 • C. After the 18th of April, a sensitive increase (more than 2 • C) in SST occurred at a regional scale (Supplementary Materials, Figure S2A). SS mean distribution map (Supplementary Materials, Figure S2C) shows the effect of continental waters (Volturno River) which reduced SS in the nearshore. In the study site, the SS reached 38 psu on the 9th-11th of April 2018, slowly decreasing until the end of the survey (Supplementary Materials, Figure S2C).

Assemblages
A total of 35 families of meso-and macroplankton taxa has been recorded over the whole studied period, together with the associated environmental factors (Supplementary Materials, Table S1). The species list included pelagic predators (Forskaliidae, Ommastrephidae, Pelagidae, and Phronimidae) as well as planktivorous species (Alciopydae, Fritillariidae, Leucotheidae, Mysidae, and Salpidae) representing exploiting food characteristics of the current fronts. Details of the trophic and reproductive behaviour of these species are reported in the Supplementary Materials (Plate I-III).

Role of Environmental Factors
The PCA model computed with SPSS explained more than 90% of the variability for the considered environmental variables ( Figure 2). Based on these results and the PERMANOVA, two environmental patterns were retained as main factors: the moon phase (p-value = 0.0011) and the current direction (p-value = 0.0015). The SST did not play a significant role in species composition and diversity, while showing high variability (range 14.5-19 • C) in the surveyed period.
Moonlight seems to rule the occurrence and abundance of plankton in the surface layer of the water column, with a clear increase in the occurrence of 12 taxa, consequently to the extremes of the moon phase (see Supplementary Materials, Figure S3). With the exception of Cestidae, all the other taxa included in this group showed a negative correlation to moonlight, as for Nereidae and Sepiolidae, which was also correlated to the record of mating behaviour (see Supplementary Materials, Plate I-II). Juvenile and larval stages of benthic and benthopelagic species, as Pleuronectiformes, Phycidae, Anguilliformes, and Scyllaridae, also negatively responded to moonlight. environmental patterns were retained as main factors: the moon phase (p-value = 0.0011) and the current direction (p-value = 0.0015). The SST did not play a significant role in species composition and diversity, while showing high variability (range 14.5-19 °C) in the surveyed period. Moonlight seems to rule the occurrence and abundance of plankton in the surface layer of the water column, with a clear increase in the occurrence of 12 taxa, consequently to the extremes of the moon phase (see Supplementary Materials, Figure S3). With the exception of Cestidae, all the other taxa included in this group showed a negative correlation to moonlight, as for Nereidae and Sepiolidae, which was also correlated to the record of mating behaviour (see Supplementary Materials, Plate I-II). Juvenile and larval stages of benthic and benthopelagic species, as Pleuronectiformes, Phycidae, Anguilliformes, and Scyllaridae, also negatively responded to moonlight.
According to the observed pattern, the most abundant taxon during the main Eastern flow direction was represented by Forskaliidae, while in the least North-Western flow the most abundant taxon was represented by Alciopidae.
The PCA model performed with CANOCO 5 shows correlations among moon phase, current directions, and taxa. It is possible to observe how most taxa took the opposite direction as the moon phase progressed and how the main current blows eastward (Figure 3a). Arrows indicate the direction of the taxa at increasing values of the environmental variables, with longer vectors representing a greater range of variation in the observed values ( Figure 3b). A stronger correlation is associated to more acute angles between vectors and the axes. The x axis explains 57.6% of the variation among parameters at each taxon, whereas the y axis explains 47.1% (Figure 3a, b). According to the observed pattern, the most abundant taxon during the main Eastern flow direction was represented by Forskaliidae, while in the least North-Western flow the most abundant taxon was represented by Alciopidae.
The PCA model performed with CANOCO 5 shows correlations among moon phase, current directions, and taxa. It is possible to observe how most taxa took the opposite direction as the moon phase progressed and how the main current blows eastward (Figure 3a). Arrows indicate the direction of the taxa at increasing values of the environmental variables, with longer vectors representing a greater range of variation in the observed values (Figure 3b). A stronger correlation is associated to more acute angles between vectors and the axes. The x axis explains 57.6% of the variation among parameters at each taxon, whereas the y axis explains 47.1% (Figure 3a,b).  The hierarchical cluster analysis performed with Primer v7 shows 5 main clusters based on the taxa presence and the current direction as environmental factor (Figure 4).  The hierarchical cluster analysis performed with Primer v7 shows 5 main clusters based on the taxa presence and the current direction as environmental factor (Figure 4). According to the statistical analyses on the relation between environmental constrains and species diversity, the contribution of the eastern currents in shaping plankton communities is well supported (see Figure 4). Lowest level of plankton diversity is related to north western stream, which can be considered local, as the diving site results covered by the island from this direction (0.8 nm off the bank). Species diversity mainly reflects typical pelagic assemblages, whereas the occurrence of emergent zooplankton (Amphipoda, Nereidae) and more benthic species (Octopodidae) can be considered mainly influenced by the contribution of the nearby shallow water banks ("Secca delle Formiche") and affected by moon phase rather than current direction. The occurrence of other larval and juvenile stages, as for Anguilliformes, Phycidae, Pleuronectiformes, and Scyllaridae (Figure 4), showed a clear convergence in the contribution of current direction ( Figure 4) and moon phase (Supplementary Materials, Figure S3), which may be related to a shared behavioural form of detection avoidance. According to the statistical analyses on the relation between environmental constrains and species diversity, the contribution of the eastern currents in shaping plankton communities is well supported (see Figure 4). Lowest level of plankton diversity is related to north western stream, which can be considered local, as the diving site results covered by the island from this direction (0.8 nm off the bank). Species diversity mainly reflects typical pelagic assemblages, whereas the occurrence of emergent zooplankton (Amphipoda, Nereidae) and more benthic species (Octopodidae) can be considered mainly influenced by the contribution of the nearby shallow water banks ("Secca delle Formiche") and affected by moon phase rather than current direction. The occurrence of other larval and juvenile stages, as for Anguilliformes, Phycidae, Pleuronectiformes, and Scyllaridae (Figure 4), showed a clear convergence in the contribution of current direction ( Figure 4) and moon phase (Supplementary Materials, Figure S3), which may be related to a shared behavioural form of detection avoidance.

Behavioral Traits of T. trachypterus
Over the whole study period, a total of 18 juveniles of T. trachypterus were observed, photographed and filmed during SCUBA dives at night between the 16th and the 28th of April 2018. Prolonged observations (>2 min) of the specimens were possible only for 10 fish, given the high avoidance behaviour of the species. Juveniles were observed swimming between 0.5 and 7 m below the surface, with the body obliquely or almost vertically pointed with the head towards the surface. Locomotion was performed through undulating movements of the long dorsal fin only, except for the anterior prolonged dorsal fin rays that showed limited oscillating movements (Supplementary Materials, Video 1-10). A faster escaping movement was displayed using the whole body, as showed in Supplementary Materials, Video 8. Pelvics and caudal fin remained spread, but substantially motionless or with small and occasional "snap" movements. While dorsal and pectoral fin membranes were entirely translucent, those of the wide pelvic and caudal fins showed orange blotches, with rays protruding far over the fin membrane. In some cases, juveniles oriented the dorsal or abdominal surface towards the observer, greatly limiting the area exposed to the view of the observer itself ( Figure 5C,D; Supplementary Materials, Video 5-7, 9-10). At the same time, caudal and pelvic fins were maintained orthogonal to the observer, with a significant torsion of the caudal fin up to 90 • angle ( Figure 5C). This latter behaviour was displayed for a prolonged time.
pectoral fin membranes were entirely translucent, those of the wide pelvic and caudal fins showed orange blotches, with rays protruding far over the fin membrane. In some cases, juveniles oriented the dorsal or abdominal surface towards the observer, greatly limiting the area exposed to the view of the observer itself ( Figure 5C,D; Supplementary Materials, [9][10]. At the same time, caudal and pelvic fins were maintained orthogonal to the observer, with a significant torsion of the caudal fin up to 90° angle ( Figure 5C). This latter behaviour was displayed for a prolonged time. Fish were also observed rapidly protruding their jaws, probably preying on small planktonic organisms (Supplementary Materials, Video 4), that were very abundant in the surface water column in which specimens were observed.
Four specimens were caught by hand net and preserved in alcohol for laboratory analysis. All the main morphometric parameters are reported in Supplementary Materials, Table S2. Fish were also observed rapidly protruding their jaws, probably preying on small planktonic organisms (Supplementary Materials, Video 4), that were very abundant in the surface water column in which specimens were observed.

Discussion
Four specimens were caught by hand net and preserved in alcohol for laboratory analysis. All the main morphometric parameters are reported in Supplementary Materials, Table S2.

Ponza Island Coastal front and Pelagic Biodiversity
Biological-physical interactions structure the variability of the marine environment at a wide range of spatial and temporal scales, affecting population dynamics and trophic interactions [42,43]. Such example is a front, a physical interface with gradients of water properties that include temperature, salinity, and turbidity [44,45]. Fronts occur across the world oceans, ranging from basin-scale features to small river plumes, and can be persistent or ephemeral [44,46]. Frontal zones are often associated with enhanced biomass, and may serve as an important foraging grounds by aggregating species from multiple trophic levels [47][48][49], being key habitats for successful energy transfer through food webs. In situ observations and modelling studies suggest that mesoscale and sub-mesoscale processes may affect biodiversity in the Mediterranean Sea, especially where coastal morphology and intense wind stress lead to the formation of energetic filaments [50,51], contributing to the dispersal of coastal inputs towards the open sea, along with plankton.
Ponza Island is located in the western sector of the Pontine volcanic archipelago and originated during the Plio-Pleistocene on the outer margin of the Latium continental shelf. The island slope rapidly sinks to deep waters bottoms, especially on the east and south-west side (Figure 1). Winds are known to support the organic enrichment process of offshore waters, bringing the contribution of coastal freshwater input. Available time lapse records of Ponza Island and Gaeta showed three prevalent wind directions, the eastern one of which can be considered the most relevant, both in terms of intensity and frequency. This is clearly shown by the SS map presented in Supplementary Materials, Figures S1C and S2. The overlap of deep water pools, eutrophic waters, and wind driven currents may lead to an organic enrichment in the study site, due to both the upwelling and the contribution of coastal freshwater waters, supporting an increase in plankton biodiversity. Vertical transport has a well known effect in increasing the biodiversity of the photic layers only in restricted areas, as, for example, where convection is sufficiently deep, in a small number of frontal regions and in the few upwelling sites. Coastal inputs have thus an important role in sustaining food webs in the whole Mediterranean Sea [52].
Despite the data acquired in this expedition they cannot be considered exhaustive to completely depict these processes (both in terms of temporal and sampling extension), they show intriguing insights into the relation between local hydrodynamism and plankton diversity. This latter conclusion is enhanced by the record of a number of rare species, some of which appear to be far less rare than thought. An example is the case of the larval stages of the tripod fish Bathypterois sp., for which, to the best of our knowledge, we provide the first available photo in environment; see Supplementary Materials, Plate IV). Despite the short period of the investigation conducted in this work, over the annual cycle, mesozooplankton abundance in offshore waters of the Tyrrhenian Sea oscillated within a narrow range and revealed lower seasonal variability [53,54], with peaks occurring in April.
Moonlight has been documented to affect the distribution of plankton biodiversity, with cyclic patterns of vertical movement synchronised with variations of irradiance in all aquatic habitats [55]. Species abundances are higher during new moon and lower during full moon, suggesting that predation pressure from planktivorous species may affect plankton abundance [55]. Moon phase is also used by species to synchronize reproduction [55].
Occurrence of species during the studied period showed two different scenarios according to the moon phase. Whereas the abundance of most taxa [23] seemed not to be directly correlated only with moon cycle, for 12 taxa, there was a clear correlation between their occurrence and moon phases (Supplementary Materials, Figure S3

Ecological Insight on the Uncommon Mesopelagic Fish Trachipterus trachypterus
Trachipterus trachypterus is considered rare and with low population densities; however, the mesopelagic habits of the species could lead to an underestimation of its abundance [28]. This fish is only occasionally caught with professional fishing gear, and usually as one or two specimens. Most of the specimens are caught by longline or trawl at depth of 370-700 m [21,24,28,56]; in other cases, they were caught by hand, or observed near the surface or stranded [20,28]. In the present study, the species was observed and caught in exceptionally large numbers. The presence of several juveniles of T. trachypterus in the study area could be related to the favorable hydrological and ecological conditions typical of banks and nearby areas. Banks are important areas for fish biodiversity and often host poorly known and rare fish species [57]. Uncommon species, such as T. trachypterus, could be good indicators of environmental changes; however, data about the ecology and biology of this species remain scarce.
In ecological communities, it is generally believed that the principal survival option available to marine fish larvae is a fast growth and a quick enhancement of swimming performances [58][59][60]. The occurrence of T. trachypterus displayed a circa-lunar cycle, in contrast with what observed with juveniles of other fish species. A more in depth analysis showed that specimens collected referred to two different size ranges, with the bigger ones (229-157 mm) sampled during almost full moon nights (Supplementary Materials, Table S2, specimens 1 and 2) and the smallerones (91-58 mm) during darkest nights (see also Supplementary Materials, Table S1). Despite the fact that this correlation must be supported by further observations, it is reasonable that the species behaviour changes according to the size, and that it occurs in the surface layer of the pelagic domain as a consequence of hiding or planktotrophic behaviour.
The larval fish phenotype also may have a strong impact on predation through dissembling or active defenses against predators. Diverse and complex morphologies can be found within the same family [61,62], suggesting the occurrence of different predation strategies, even if visual predation by larger fishes is likely the dominant selective force.
Batesian mimicry is a more widespread survival strategy than previously thought, especially within marine communities [63], as demonstrated by the convergence of numbers of distantly related fish taxa on the production of complex features (especially in larval stages and juveniles), like extremely long fin rays (Gadiformes, Lampridae, Lophiidae, Pleuronectiformes, and Serranidae), external guts (Gadiformes, Myctophidae), and even false eyes on stalks (Myctophidae, Stomiidae) [62,63]. Within marine pelagic environments, the evolution of visual Batesian mimicry in larval fishes must meet three fundamental conditions: (1) visual predation must be a strong source of mortality; (2) there must be a relatively common and functional model that the mimic resembles through morphology or behavior, and (3) mimics with only a slight resemblance to the model receives some degree of protection [64,65].
Within this expedition, we had the chance to observe, for the first time, some peculiar behavioural traits of juveniles of T. trachypterus which may be linked to Batesian mimicry.
The almost vertical position taken by juveniles of T. trachypterus during swimming and the dorsal-fin based locomotion recorded during our study is in agreement with the swimming behaviour previously reported for juveniles of the same species [66]. A similar swimming behaviour was recently reported also for juveniles of T. arcticus [67]. On the other hand, new data emerge from our study: 1-the "dorsal/ventral surface exhibition", in which the juvenile fish oriented its dorsal/ventral surface towards the observer; 2-the potential Batesian mimicry. We hypothesise that these two behaviours are related. Dorsal-orienting swimming behaviour may represent a strategy to hide the general shape of the fish and, at the same time, the extension and orientation of caudal and pelvic fins to the potential predator may enhance its appearance as a noxious (or not palatable) organism. The small occasional snap movements of these fins further recall a jellyfish species as a model for the Batesian mimicry (e.g., Pelagia noctiluca). Noteworthy, the elongated rays of the anterior part of the dorsal fin and the small orange expansions on these rays are similar in appearance to the cormidia of siphonophores. These peculiar morphological traits seem not to be conserved in the adult and therefore can be considered a functional adaptation to reduce predation risks during the juvenile epipelagic stage of the species. Indeed, this has already been suggested for juveniles and larvae of other mesopelagic fish species, in which particularly long and delicate fins (often in association with specific behaviours and colors) play an important role in reducing mortality due to predation [63,[68][69][70].

Conclusions
Despite the numerous investigations of the last decades, the emerging picture of plankton dynamics in the Mediterranean Sea is still largely unexplored. Many factors, such as the structure of the seafloor and continental water inflows, may represent core drivers in the distribution of plankton diversity, structuring complex and dynamic trophic interactions between coastal and pelagic environments. In this research, we show how the vertical transport of coastal nutrients close to the continental shelf boundary sustains a high planktonic diversity and the occurrence of rare mesopelagic fish species. Furthermore, for a total of 12 taxa, a clear correlation between their presence and moon phases is reported.
The concept behind AQUATILIS field-based approach, along with the use of modern filming equipment, allowed us to describe, for the first time, the behaviour of the uncommon mesopelagic fish T. trachypterus, depicting the presence of Batesian mimicry in juveniles, and to observe and photograph rare and little known species, such as the larva of Bathypterois sp. Our results support the role of field studies as a fundamental instrument improving our knowledge on complex ecosystems, as those inhabited by the mesopelagic communities.
Supplementary Materials: The following are available online at http://www.mdpi.com/2077-1312/8/8/581/s1: Table S1-Planktonic taxa surveyed at the diving site; Figure S1-(a) and (b) show respectively the wind rose of Ponza (2011)(2012)(2013)(2014)(2015) and Gaeta (2007-present). Map (c) and (d) show respectively the mean distribution of SST and Salinity from 1987 to present; Figure S2-In the panel (a) is presented in descending order the distribution of SST during 9th-11th of April, during 12th-18th of April and during 19th-28th of April. In the panel (b) the polar scatter of wind indicates the direction (polar orientation) and intensity of wind (x-axis). In the panel (c) is presented in descending order the distribution of Salinity during 9th-11th of April, during 12th-18th of April and during 19th-28th of April; Figure S3-Taxa occurrence according to moon phase; Figure S4. The four sampled juvenile specimens of Trachipterus trachypterus from the central Tyrrhenian Sea (A-C), measurements are reported in Table S2; Table S2. Morphometric measurements of the four juvenile specimens of T. trachypterus reported in Figure S4