Are there distribution patterns and population structure differences among demersal fish species in relation to Antarctic benthic communities? A case study in the Weddell Sea

Despite the general belief that the Southern Ocean harbors low fish biodiversity, the Weddell Sea hosts one of the richest fish communities in the region. Parallelly, the Weddell Sea is also known for the presence of dense and diverse macrobenthos. Most macrobenthic invertebrates, such as gorgonians, sponges and bryozoans, are considered ecosystem engineers as they generate a three-dimensional structure that increases habitat heterogeneity. This structural complexity serves as a refuge against predators as well as a nursery ground for many organisms, including fish species. By analyzing video transects recorded by a Remotely Operated Vehicle, we investigated density, spatial distribution and size-frequency of populations of the demersal fish species inhabiting macrobenthic communities in the southernmost part of the Weddell Sea. We also attempted to unveil whether there is any relationship between benthic and fish communities and substrate, as well as some fish behavioral patterns. The dominance of juveniles in the surveyed fish assemblages provides evidence that, at this life stage, some fish species appear to be positively associated with complex benthic communities conformed by bryozoans, sponges and gorgonians which are more common in sand matrix with sparse rocks substrates. Moreover, about 37% of all specimens recorded were resting on benthic invertebrates or were using them to hide, implying that Antarctic benthic communities might offer suitable habitat. As such, it can be concluded that there was an apparent relationship between certain species of fish and the different benthic communities, yet the exact triggers and/or factors behind such an association remain partially elusive.


Introduction
The Southern Ocean is the water body occupying a total area of 34.8 million km 2 comprised between the coast of the Antarctic continent and the Antarctic Polar Front (also known as Antarctic Convergence, located at about 50°S in the Atlantic and Indian sectors and 60°S in the Pacific sector) (Hofmann 1985). It is characterized by low water temperature and seasonal primary production due to seasonal ice cover, which has remained practically unaltered over the last 20 million years (Clarke and Crame 1989). Yet, despite its large dimension and long evolutionary history, fish biodiversity is less diverse than expected (Eastman 2005).
Fish species with benthic habits are the most diverse and abundant fish group in Antarctic waters, representing more than 60% of all Antarctic fish species, with the Notothenioidei suborder being the major component in terms of biomass and number of demersal species (DeWitt 1971;Eastman and Hubold 1999;Eastman et al. 2013). In fact, most of the notothenioid species are endemic to the Southern Ocean (Andriashev 1965;Kock 1992). This high endemism has occurred due to the Antarctic Polar Front, which acts as a zoogeographic boundary, as well as the constant and severe conditions of the Southern Ocean waters (Eastman and Clarke 1998;Peck 2018). The Notothenioidei currently includes 9 families and 140 species, from which 5 families and 110 species are exclusively present in Antarctic waters, while only 4 families and 30 species are non-Antarctic (Eastman and Eakin 2021; Daane and Detrich 2022).
Despite the aforementioned generalization that the Southern Ocean presents lower fish diversity than expected, the Weddell Sea stands out as a hotspot of fish diversity within the area, harboring the highest number of fish species within the Southern Ocean (Hubold 1991). Within the High-Antarctic Zone (sensu Kock 1992), the Weddell Sea is characterized by a wide shelf as deep as more than 1000 m (Ekau 1990;Anderson 1999). It also shows the most durable seaice coverage, which plays a fundamental role in physical and biological processes, such as the formation of the coldest water masses and productivity (Clarke and Ackley 1984;Peck 2018). Furthermore, the bottom of the High-Antarctic Zone stands out due to the presence of diverse and dense microbenthic communities, known as the Antarctic Marine Animal Forests (AMAF) (Arntz et al. 1994;Gutt et al. 2004;Gili et al. 2006). Most of the macrobenthic sessile organisms in these communities, such as gorgonians, sponges and bryozoans, are considered ecosystem engineers since they generate a three-dimensional structure that increases habitat heterogeneity (Jones et al. 1994;Gili and Coma 1998;Rossi et al. 2017). This increase in structural complexity provides a wide array of niches in the areas those species occur, which in turn might serve as a refuge against predators or as nursery grounds for many other species, including fish species (Gutt et al. 1994;Gratwicke and Speight 2005;Buhl-Mortensen et al. 2010;La Mesa et al. 2019). In this regard, previous studies have noted that some fish species settle on elevated benthic structures such as sponges, a behavior considered characteristic for certain species (Ekau and Gutt 1991;Gutt and Ekau 1996). Yet, while macrobenthic sessile organisms might increase the diversity and density of other taxa, their presence does not always lead to an increase in diversity, as the degree of habitat modification depends on the density and size of the engineer species (Jones et al. 1994;Cerrano et al. 2010).
In the past, it was only possible to study marine ecosystems through destructive methods but, thanks to technological advances, that has changed. Technological advances and the development of remotely operated vehicles (ROVs), manned submersibles, autonomous underwater vehicles (AUVs) and ocean floor observation systems (OFOS) have made it possible to go further in understanding the marine environment, enabling an in situ overview of the ecosystem (Gutt and Ekau 1996;Santín et al. 2018;La Mesa et al. 2022). While catch sampling is a destructive technique that normally results in a disordered sampling of animals, underwater photography and video are conservative methods that significantly minimize impacts on sessile organisms. Image methods allow for direct observation, as well as a much greater study area in terms of extension and depth, compared to using destructive methods (Ninio et al. 2003;Matarrese et al. 2004;Rossi et al. 2008). As so, underwater video and photo technology has been used for studies on biomass (Mühlenhardt-Siegel et al. 1988;Pabis et al. 2011), animal behavior (Ekau and Gutt 1991;La Mesa et al. 2019) demography and spatial distribution of benthic (Segelken-Voigt et al. 2016;Ambroso et al. 2017) and motile species (Ekau and Gutt 1991;Lorance and Trenkel 2006;Amsler et al. 2015;La Mesa et al. 2022). Moreover, some studies have provided the same results on the distribution and abundance of fish species by comparing catch sampling and image techniques, validating the usage of the later for fish community studies (Ekau and Gutt 1991;Gutt et al. 1994;La Mesa et al. 2019). Yet, despite the possibilities offered by imaging techniques, few studies have focused on the macrobenthic species that accompany fish assemblages in different study areas (Gutt et al. 1994;Gutt and Ekau 1996;La Mesa et al. 2019).
Hence, the present study aims to fill the said knowledge gap and respond to the hypothesis of a possible relationship between complex benthic communities and greater biodiversity of Antarctic fishes. In order to achieve this, three main objectives have been pursued through video analyses: (1) identification of demersal fish species, their density, spatial distribution and size-frequency distribution populations; (2) description of fish assemblages; and (3) a study of the relationship between benthic communities and the different fish species, including any behavioral patterns observed.

Study area
The study area was located in front of the Filchner-Ronne Ice Shelf (FRIS), in the southernmost part of the Weddell Sea ( Fig. 1), which has, until now, been poorly studied due to heavy sea-ice conditions and consequently difficult access. Therefore, any new information about this region is of special relevance. This area is characterized by cold and dense water masses with a cyclonic clockwise flow (Deacon 1979;Daae et al. 2020). Near the bottom, the Weddell Sea Bottom Water (WSBW) can be found, which is slightly warmer than the Weddell Sea Deep Water (WSDW). Above this, there is the Warm Deep Water (WDW), which occupies the layer of 200-1200 m depth and its temperature is between − 1 and -2 °C (CTD data). In addition to this, the continental shelf includes two water masses that make up the Dense Shelf Water (DSW): The Ice Shelf Water (ISW) and the High-Salinity Shelf Water (HSSW). They originate from the seaice production (Haid and Timmermann 2013) and, in turn, contribute to the formation of the Weddell Sea Deep and Bottom Waters (Gordon et al. 1993;Orsi et al. 1999). All of the conditions of the shelf water masses have remained practically stable since the 1980s (Janout et al. 2021). The entire area is mainly characterized by soft bottoms with different proportions of sand, gravel and boulders (Anderson et al. 1980;Haase 1986).

Video recording
A total of six video transects were analyzed, covering a total area of 1,845.9 m 2 between 251 and 361 m depth, as part of the multidisciplinary scientific cruise PS82 (ANT XXIX/9) expedition on board the RV Polarstern, conducted from December 2013 to March 2014.
To study the composition and distribution of Antarctic fishes and their relationships with invertebrate benthic macrofauna, an inspection-class ROV (Ocean Modules V8Sii) was deployed at six stations in the area of the Filchner Trough (Fig. 1). The videos were acquired with a high definition (HD) video camera (Kongsberg oe14-502) recording with a viewing angle of 29° × 45° (vertical × horizontal) and equipped with two lasers with a 4 cm separation that serves as a reference scale, used for the spatial density and fish size analyses. Furthermore, the ROV was equipped with ultra-short baseline (USBL) positioning data. All of the ROV video material pertaining to this study is available from the data publisher PANGAEA® at www. panga ea. de (Owsianowski et al. 2017).

Video analysis
Video transects were analyzed using Premiere Pro CC version 12.1.2 (Adobe Inc.) software and following the methodology described in Gori et al. (2011). Each track was edited before the analysis in order to remove pauses, loops or video sequences where the image quality was not suitable for the study, for example, when there was sediment resuspension. After editing, every benthic organism and fish observed along the transect was annotated, within a section of 0.3 m width around the central position of the laser beams. A reference time was assigned to each benthic organism and fish was when it was observed and their position along the transect was recorded with geographical coordinates. In addition, following the same analysis procedure, substrate type and depth was annotated along each transect. Substrate types were classified into three categories, namely: sand matrix with rocks, gravel matrix, and rock matrix. Seabed slope was also analyzed, and it was consistently found to be 0.

Data analysis
Once fish species and benthic organisms had been identified and georeferenced, it was decided to take an area of 50 m 2 as a sampling unit (0.3 m width and 166.5 m length) to analyze the data statically. Since the target fish of the study are demersal, but the benthic communities are mainly composed of sessile animals and could present a highly patchy distribution, the sampling unit chosen to study both of them is 50 m 2 . Each sampling unit was characterized by a type of macrobenthic community and substrate. The slope of all the transects was zero, and the sea bottom was horizontal. This factor was therefore not considered to characterize the sampling units, and not used for statistical analysis either. Depth was also ruled out as an influencing factor in this study because an ADONIS test was carried out using the vegan package (Oksanen et al. 2013) from RStudio (RStudio Team 2019), and no influence of depth was found on the studied fish assemblages. The study area corresponds to a unique water mass (WDW), and that could be the reason Fig. 1 Study area. Location of the video transects in the study area in front of the Filchner Ronne Ice Shelf in the Weddell Sea why depth did not appear to be an influent factor in this case. Regarding benthic classification, four benthic communities were defined by the most abundant engineer species: bryozoan community (B), when more than 85% of the macrobenthic organisms were bryozoans and no other macrobenthic phylum reached 10% coverage by itself; bryozoan and sponge community (BS), when bryozoan represented more than 60% and sponges more than 15% of the macrobenthos; sponge and bryozoan community (SB), which was dominated by sponges, with any inverse distribution (≥ 60% and ≥ 15%, respectively); and lastly, bryozoan, gorgonian and sponge community (BGS), where bryozoans represented more than 60%, gorgonians more than 20%, and sponges more than 10% (Fig. 2). Besides the types of benthic communities, two kinds of bottom substrate were defined: sand matrix with sparse rocks and gravel matrix. Each sampling unit usually presented a combination of both, but was defined by the most abundant (with a percentage of 80% or more).

Density and spatial distribution and size-frequency of fish populations
For statistical analysis, all of the transects were divided into sampling units of 50 m 2 (0.3 m width and 166.5 m length). A total of 37 sampling units were obtained from six transects, composed of four sampling units of transect 49, four more of transect 81, eight of transect 86, seven of transect 128, seven sampling units of transect 136, and seven more of transect 170, respectively. In order to study the distribution for the most representative species, a geographically-referenced map was generated using QGIS version 3.16 (QGIS. org 2022), in which both fish density and the corresponding type of benthic community of the sampling unit were represented. In addition, the species composition was examined and their density (ind. 100 m −2 ) were quantified.
The size-frequency of fish populations was computed only for those species with more than 13 individuals documented. The ROV lasers were set at a known distance from each other (30 cm) and provided a scale for the measurement of fish individuals. When a fish appeared in the video transect, a photo was taken when the lasers were parallel to the individual, which was then analyzed for measurements by means of the software Macnification version 2.0.5 (Schols and Lorson 2008) In order to compare species with different size population, size-frequency histograms were performed for each species by means of the ggplot2 package (Wickham 2016) from the R software platform (RStudio Team 2019). In the High Antarctic, more than 90% of notothenioids do not reach 45 cm total length, and other fishes like artedidraconids are typically less than 15 cm. The variability of maximum fish sizes between observed fish species is the reason why the sizes of fishes are classified in small ranges of 3 cm.

Fish assemblage composition
Aiming to find some assemblage present in the study area, abundances of fish species were represented by means of a dendrogram and non-metric multi-dimensional scaling (nMDS). Before the analysis, the abundance data were transformed using the square-root method to normalize them, and then distances between pairs of sampling units were calculated using the Bray-Curtis dissimilarity index. To obtain the dendrogram, the hclust function of vegan package (Oksanen et al. 2013) was used from the RStudio platform (RStudio Team 2019). In order to elaborate the nMDS, the metaMDS function was applied, again available via the vegan package (Oksanen et al. 2013). Three sampling units were discarded from the statistical analysis, having no fish because of the analysis requirements. Dendrograms and nMDS are only representations of how sampling units are grouped, as an exploratory method to graphically visualize if there might be some kind of relationship among the data. The ADONIS test was applied using the vegan package (Oksanen et al. 2013) to determine any statistical differences among fish assemblages with RStudio software (RStudio Team 2019). To carry out these analyses, individuals classified at the genus or family level have not been considered, because they could be different species and, therefore, could have different substrate and benthic community preferences, which could introduce errors in the interpretation of the results.

Relationship with environmental features
The distribution of fish species, like all other animals, is determined by environmental features. In this sense, in the present study, the environmental factors determining fish distribution were explored by canonical correspondence analysis (CCA), which is a multivariate analysis. CCA allows the identification of probable relationships between fish species and environmental factors (Greenacre 2013). The environmental features considered were depth, substrate, and benthic community. Sampling units without fish were excluded because of the statistical test requirements. For this reason, only a total of 34 sampling units were analyzed. Additionally, the individuals not classified at species level were deleted for this analysis, for the same reason as in the dendrogram analysis. Another statistical test to find out how fish species are related to environmental factors was Indicator Value Analysis (IndVal), which was described by Dufrêne and Legendre (1997). The analysis presents a table of indicator species that are quite abundant and significantly characteristic in a type of habitat and with a p-value of less than 0.05. In this study, this was tested with the different types of benthic communities and substrates previously defined, using the indval function of the labdsv package (Roberts and Roberts 2016) of RStudio (RStudio Team 2019). In addition, the richness and diversity of fish species were quantified by means of the Shannon Diversity Index, which was calculated for the different types of benthic communities and substrates. This is represented in two box plots carried out via the ggplot package (Wickham 2016) of RStudio (RStudio Team 2019).

Density, spatial distribution and size-frequency of fish populations
A total of six stations of video transects of the seabed were conducted in the Weddell Sea from 251 to 361 m water depth, covering an area of 1.845,9 m 2 . Thanks to the high resolution of the camera, the majority of the fishes observed were identified at species level. The rest of them were identified within a family or genus category. In these cases, the high swimming speed or the incomplete vision of the body prevent identification with the highest level. Overall, 12 fish genera and 14 fish species were identified. A total of 414 specimens were recorded that could be identified to species level (87.2%), and the rest of the individuals have been identified to genus (11.3%) and family level (1.4%). Icefish of the genus Chionodraco were not identified to the species level, as we were not able to clearly distinguish between them. The major diversity and density of fish were found at stations 128, 136 and 170 (Tables 1 and 2). Fishes recorded belonged to four families, such as Artedidraconidae, Bathydraconidae, Channichthyidae and Nototheniidae. The family Nototheniidae was the most abundant, including in decreasing order of abundance Trematomus scotti (Boulenger 1907) (30.7%), T. lepidorhinus (Pappenheim 1911) (13.8%), and T. eulepidotus (Regan 1914) (10.4%). The Channichtyidae was the second most abundant and speciose family, encompassing at least six different species, followed by the Bathydraconidae with four species, whereas the Artedidraconidae family was represented by a single specimen of Pogonophryne scotti (Regan 1914). In terms of spatial distribution, only the six most abundant species were represented in maps (Fig. 60). Chionodraco spp. and Chaenodraco wilsoni (Regan 1914) were observed only in northern stations, while T.scotti, T. eulepidotus, T. lepidorhinus and Pagetopsis maculatus (Barsukov and Permitin 1958) were distributed both in the northern and southern stations. The size-frequency distributions are summarized in Fig. 3. The size population of the fish species was described only when the number of observations of each species amounted to more than 13 individuals (Fig. 4). All nototheniids belonged to Trematomus genus and reached a maximum total length (TL) of 21 cm, with a modal size of 9 cm. This pattern of size distribution can also be observed with the specimens of Pagetopsis, whereas Cryodraco antarcticus (Dollo 1900) and Chionodraco spp. showed the largest sizes reaching more than 30 cm TL. In Chionodraco spp. it should be highlighted the absence of individuals smaller than 18 cm TL.

Fish assemblage composition
To find out the fish assemblage composition, a dendrogram was performed representing the sample clustering based on a Bray-Curtis dissimilarity matrix (Fig. 5). Cluster analysis separated five assemblages (at 60% of dissimilarity), with the 8, 9 and 13 sampling units being considered as outliers.
Most of the sample units classified as BGS benthic community (red color in the cluster, Fig. 5) are grouped in the representation of the cluster, and all of them are characterized by sand as the most abundant substrate. The second group contains all the SB sampling units (grey color, Fig. 5), and the majority of them correspond to the matrix of gravel substrate. The last group shows the majority of the BS sampling units (green color, Fig. 5) and this grouping can be divided into two depending on the substrate. As a result, fish distribution appears to be related benthic communities and substrate. The non-metric multi-dimensional scaling (nMDS) was another analysis to pursue some relation between benthic communities, substrate and fish distribution (Fig. 6). The ordination of the sampling units that correspond to the SB Fig. 3 Geographical distribution in the study area of the most common fish species represented based on density and benthic community. Bryozoan community (B), bryozoan and sponge community (BS), sponge and bryozoan community (SB) and bryozoan, gorgonian and sponge community (BGS) community are clearly grouped; the same occurs with the BGS community, and some grouping formed by the sampling units of the BS community could also be seen. Consequently, there could be a correlation between fish biodiversity and density and these types of benthic communities. On the other hand, despite some sampling units of the B community also being grouped, others are not. Regarding the substrate, there is a spatial ordination of the different types of substrate that proves that there also seems to be a relationship between this environmental factor and the distribution and abundance of fish species. All this coincides with the result obtained from the dendrogram. In addition, the ADONIS test confirmed that benthic communities are different (p-value < 0.01) ( Table 3). Comparing benthic communities in pairs, results show that all of them were different from each other, except for the B and BS communities (p-value = 0.238) ( Table 3). The two categories of substrate, sand matrix with some rocks and gravel matrix, also appear to be different from each other (p-value < 0.001).

Relationship with environmental features
Regarding geographical distribution and benthic communities' preferences of the studied fish species (Fig. 3; Table 4), P. scotti and Gymnodraco acuticeps (Boulenger 1902) were observed exclusively restricted to the BS community, whereas Pagetopsis spp. and Trematomus loennbergii (Regan 1913) were recorded within the BGS community. Gerlachea australis (Dollo 1900), Racovitzia glacialis (Dollo 1900), Chionodraco spp., C. antarcticus, Prionodraco evansii (Regan, 1914) and some non-identified individuals belonging to the Channichthyidae family were associated to the BS and BGS communities. Pagetopsis macropterus (Boulenger 1907) and C. wilsoni occurred within three benthic communities B, BS and BGS, whereas the remaining species were observed everywhere. The highest densities were found within the BS community, where T. scotti density was 8 ind. m −2 , and also in the BGS community, where the maximum density of C. wilsoni was 6.75 ind. m −2 , but the latter was the benthic community with the highest total density of fish species. All species were observed on sampling units with a matrix of sand with rocks, yet some species were also found on sampling units defined as having gravel bottoms, but with low densities, such as P. maculatus and P. macropterus. Conversely, T. lepidorhinus is the only species that showed higher density on gravel bottoms, yet it occurred on both types of substrate.
As observed in the results obtained with the dendrogram and nMDS, it seems that fish assemblage has a relation with benthic communities and substrate. To delve further into this, a CCA was performed that shows the relations between fish species with the benthic community and substrate (Fig. 7). Sand substrate with rocks seems to have a strong relationship with the vast majority of identified species, except for T. lepidorhinus. Regarding the benthic communities, most fish species are positively related to BGS communities, such as G. australis, T. loennbergii, C. wilsoni and C. antarcticus. Differently, G. acuticeps and P. scotti prefer the BS community, and Cygnograco mawsoni (Waite 1916) seems to be more associated with B community. The CCA analysis is complemented by a statistic analysis, the Indicator Value Analysis (IndVal), which shows, with a p-value of less than 0.05, that T. lepidorrhinus is characteristic of the SB community, Chionodraco spp. are common in the BS community, and the C. wilsoni, C. antarcticus, T. scotti, G. australis and T. loennbergii species are characteristic in the BGS community. Moreover, this same test shows, with a p-value of < 0.05, that T. lepidorrhinus is typical on gravel bottoms, while T. scotti, T. eulepidotus, P. maculatus and C. wilsoni are more characteristic on sandy ground.
In agreement with our previous results, Shannon's index was highest in BGS community and lowest in SB community (Fig. 8). Regarding substrate, sand matrix with sparse rocks shows the highest Shannon values, so this type of substrate appears to favor fish biodiversity and abundance. In addition, it should be said that the most complex benthic communities studied (BGS) were observed only in sampling units with sand matrix with rocks, so a relationship might also be deduced. Conversely, in sampling units with a gravel matrix, simpler benthic communities have been observed. In fact, most of them were B communities. Indeed, the only three fishless sampling units that were excluded from the dendrogram and CCA analyses are predominantly characterized by a gravel substrate and B community. Therefore, the results seem to be consistent. The biodiversity indexes were calculated considering all the fishes identified at species and genus level. Following this, one genus could be represented by more than one species, so the biodiversity index calculated could be the same or less than real biodiversity index.

Discussion
The analysis of the images has allowed the identification of a total of four fish families: Artedidraconidae, Bathydraconidae, Channichthyidae and Nototheniidae of the perciform suborder Notothenioidei. These families are the most abundant in the High-Antarctic Zone, especially on the continental shelf of the Weddell Sea (Ekau 1990). Nevertheless, other families occur in the High-Antarctic Zone, like Liparididae and Zoarcidae, which are less abundant but equally diverse as those mentioned, were not observed in this survey (Anderson 1990). In this sense, some fishes may be so elusive to the presence of submersible devices, to the point of becoming undetectable (Katsanevakis et al. 2012). As a result of this circumstance, the number of observations in this study could be underestimated. Throughout the study area, the fish assemblage was mainly characterized by representatives of the family Nototheniidae, being Trematomus the most abundant genus, followed by the Channichthyidae family, whose most representative genera were Pagetopsis and Chionodraco (Table 2), in accordance with data from previous studies in the southeastern Weddell Sea and other regions of the High-Antarctic Zone (DeWitt 1971;Ekau 1990;La Mesa et al. 2019). This contrasts with the Antarctic Peninsula that belongs to Low-Antarctic Zone also called Seasonal Pack-Ice Zone, where the Nototheniops spp. and Notothenia spp. were predominant (Targett 1981;Kock and Stransky 2000).
Matching the size-frequency distributions of fishes with size ( Fig. 4) at sexual maturity reported elsewhere, it was possible to infer the population structure of the species  Duhamel et al. 1993;La Mesa et al. 2008), so the population surveyed consisted exclusively of juveniles and subadults (Fig. 4). Conversely, the population of T. scotti was composed of both juveniles and adults, as they reach sexual maturity at 12-13 cm TL (Duhamel et al. 1993). In this sense, although belonging to the same genus, different species might present contrasting population demographics, highlighting the importance of accurate identification of the specimens. Among channichthyids, P. macropterus and P. maculatus attain sexual maturity at 18-19 cm TL, whereas for C. wilsoni it is at 23-25 cm TL (Duhamel et al. 1993;Kock et al. 2008), so their populations consisted solely of juveniles and, in the case of C. wilsoni, with the additional presence of just a few adults (Fig. 4). Similarly, all the specimens of C. antarcticus were juveniles, as they reach sexual maturity at a size larger than 35 cm (Kock and Jones 2002), which were not observed in the present study (Fig. 4). These results are consistent with previous studies (La Mesa et al. 2020), confirming that early juveniles show a pelagic distribution and feed on krill, mysids and fish larvae, while adults are demersal and more sedentary, and hunt prey following the "sit and wait" feeding strategy (Ekau and Gutt 1991;Kock et al. 2013). Conversely, the populations of Chionodraco spp. included both juveniles and adults as C. myersi (DeWitt and Tyler 1960) and C. hamatus (Lönnberg, 1905) attain sexual maturity at 25-30 cm and 35 cm TL, respectively) (Duhamel et al. 1993;La Mesa et al. 2003). In summary, T. eulepidotus, T. lepidorhinus, T. scotti, C. wilsoni, P. maculatus and P. macropterus, six of eight species whose size population has been analyzed, exhibited a population dominated by juvenile stages. These species have in common that they reflect some dependence on the water column, since much of their diet is made up of krill or other euphausiids (Targett 1981;La Mesa et al. 2004). These results lead to hypothesize that adult individuals might show more active behavior for prey in the water column, while juveniles belonging to these species might prefer the most complex benthic communities formed by sponges, bryozoans and gorgonians, which provide a suitable place to hide from predators until they reach a larger size.
In this sense, observation in situ by means of video transects also allowed us to assess the relationships between fish and benthic communities, as well as to evaluate their habitat preferences and specific behaviors. As has been mentioned previously, benthic organisms such as sponges, gorgonians and bryozoans offer protection, food and a place for breeding and nursery for different notothenioid fishes (Ekau 1990;Ekau and Gutt 1991). In agreement with previous studies (Ekau and Gutt 1991;Gutt and Ekau 1996;La Mesa et al. 2019) T. lepidorhinus prefers benthic communities characterized by populations of sponges and bryozoans. In fact, during the course of this study T. lepidorhinus had often been observed resting on or hiding inside volcano sponges (more than 60% of observations for this species). Similarly, several species (T. scotti, T. loennbergii, C. antarcticus, C. wilsoni and G. australis) were common where bryozoans, gorgonians and sponges are abundant (Table 4; Fig. 7). Some authors explain this preference by the three-dimensional structure that these engineer animals offer for protection against predators, being an advantage for the "sit-observe and hide" strategy of some fish species during their juvenile stage (Gutt and Ekau 1996;La Mesa et al. 2019). Further, these kinds of complex benthic communities offer suitable habitats for sedentary species that use the same strategy but for hunting. They lie in wait hidden and hunt on the prowl (Ekau and Gutt 1991). Conversely, was C. mawsoni the only fish species which showed a preferential relationship for bryozoan communities and more open areas without highgrowing sponges (Table 4; Fig. 7) (Ekau and Gutt 1991). Regarding substrate preference observed, fish diversity rose in sandy areas with rocks where the most complex benthic communities were common and consequently, the variety of ecological niches were greater (Fig. 8). Conversely, fish biodiversity decreased in gravel areas hosting simpler benthic communities (Fig. 8). In this sense, the different kinds of substrates and benthic communities studied make that fish assemblage vary depending on them. Consistent with their benthic and sedentary behaviors, about 37% of all specimens recorded were resting on benthic invertebrates such as sponges and bryozoans or were using the tridimensional structures created by the macrobenthos to hide. The rest of the observed fishes were found swimming a few centimeters above the substrate. It should be also mentioned that C. wilsoni was observed seven times out of 31 sightings guarding eggs on nets composed of flat drop stones, with a behavior previously observed in the west Antarctic Peninsula (Ziegler et al. 2017) and the Weddell Sea (La Mesa et al. 2019;Kock et al. 2008). For this reason, in order to ward off potential predators such as starfish and other fish species, nesting and parental care are common in demersal species like C. wilsoni (Kock et al. 2006).
To conclude, the distribution and abundance of fish species depend on a large number of factors not considered in this research, but which are still important, such as temperature and food availability (Ekau and Gutt 1991). Further, ice-scouring also represents a determining factor in the distribution of species at a small scale in High-Antarctic waters (Brenner et al. 2001). This study proves that seems to be a close relationship between the different species of fish observed and the different benthic communities defined by the predominant engineer organisms, and consequently, by their structural complexity. Fish species take advantage of the structures of benthic organisms to hide from predators, like T. lepidorhinus, or to hunt on the prowl as C. antarcticus does. Despite there not being a clear correlation with all of the fish species, biodiversity is indeed greater in sandy areas with rocks where complex benthic communities are more common.