Drivers of Megabenthic Community Structure in One of the World’s Deepest Silled-Fjords, Sognefjord (Western Norway)

The Sognefjord is the longest (205 km) and deepest (1308 m) fjord in Norway, and the second-longest in the world. Coast-fjord exchange in Sognefjord is limited by a seaward sill at 170 m water depth, which causes a clear stratification between water masses as the dense oxygen-poor basin water mixes slowly with the well-oxygenated water directly above from the coastal ocean. Due to the homogeneity and limited variability in the deep-water, the deep slopes of Sognefjord represent the ideal setting to study how abiotic factors influence the deep-water benthic community structure. During the summer of 2017, two remotely operated vehicle (ROV) video transects were performed to compare the megabenthic community behind the sill (water depth: 1230 to 55 m; transect length: 1.39 km; distance from sill: ∼17 km) and within the central fjord (water depth: 1155–85 m; transect length: 2.43 km; distance from sill: ∼79 km). Accompanying conductivity–temperature–depth (CTD) deployments were made to measure the in situ abiotic factors and nutrient concentrations at each transect location, while the substrate characteristics (percent cover of soft and hard exposed substrate) were documented from the video footage. Here, Sognefjord’s megabenthic community composition, distribution, and species richness were analyzed in relation to abiotic factors (e.g., depth, salinity, dissolved oxygen, chlorophyll a concentration, and percent cover of hard and soft substrata) within the fjord. Basin communities were homogeneous and characterized by sponges, echinoderms, and crustaceans, whereas the shallower regions were dominated by mobile scavengers. Contrary to other fjord-based studies, species richness and diversity were stable in the fjord basin and decreased with proximity to the sill, decreasing water depth, and at the boundary between intermediate and basin water. The findings demonstrate that highly stratified fjords support stable communities in their basins; however, further research is needed to investigate the influence water mass dynamics have on silled-fjord megafauna communities.


INTRODUCTION
Deep fjords are valuable study areas because they allow easy access to habitats that share similarities with continental shelf or deep-sea communities found in the open ocean (Bernd, 1993;Sweetman and Witte, 2008;Storesund et al., 2017). Their accessibility allows researchers to study the influence of abiotic factors on the shelf or deep-sea community ecology whilst reducing the limitations of cost and transportation that is often problematic for deep-sea research.
It is well documented that fjord communities are influenced by both the coast-fjord and/or depth gradients (Buhl-Mortensen and Høisaeter, 1993;Holte et al., 2004;Włodarska-Kowalczuk and Pearson, 2004;Storesund et al., 2017;Molina et al., 2019). The interaction between the ocean water and the fjord system helps carry seawater into the fjord, which aids in the distribution of fauna, organic nutrients, and inorganic material (Buhl-Mortensen and Høisaeter, 1993;Holte et al., 2004). However, within silled-fjords, the exchange of seawater between the coastal and fjord systems is reduced. The sill height, fjord topography, and freshwater input influences the transport of nutrients, pelagic larvae, organic matter, and dissolved oxygen, where benthic communities toward the inner fjord regions are negatively impacted due to the decreased access to resources (e.g., dissolved oxygen, organic carbon, nutrients) (Buhl-Mortensen and Høisaeter, 1993;Blanchard et al., 2010;Storesund et al., 2017). Density stratification between upper and bottom water masses is often observed and the inner basin(s) below the sill depth becomes isolated from the adjacent coastal system. This stratification leads to relatively stable temperatures and salinities within the basins (Renaud et al., 2007;Drewnik et al., 2016;Molina et al., 2019). With limited mixing within the water column, dissolved oxygen levels tend to decrease with depth, distance from the sill, and over time until a renewal event occurs, whereby more oxygenated ocean water mixes with fjord water. In general, faunal diversity and species richness is seen to decrease with increasing distance from the coastal regions and increasing water depth (Buhl-Mortensen and Høisaeter, 1993;Holte et al., 2004;Molina et al., 2019).
Deep-water stagnation is thought have a major influence on fjord basin communities, where lower oxygen concentrations can negatively impact the community structure and species composition (Blanchard et al., 2010;Molina et al., 2019), resulting in lower species richness with lower oxygen levels (Buhl-Mortensen and Høisaeter, 1993). In periods of hypoxic and anoxic conditions, where the oxygen concentration is <2.1 mL L −1 , defaunation of macrofauna and changes in faunal assemblages have been observed within the basin (Holte et al., 2005;Molina et al., 2019). In extreme cases of deoxygenation and strong stratification, a rise of acidification within fjord basins can occur (Jantzen et al., 2013), particularly if water exchange with adjacent coastal systems is insufficient.
Sognefjord is Norway's longest and deepest fjord. It is host to numerous towns and villages and has become a sought-out destination for many cruise ships, where thousands of tourists visit the fjord each year. There is concern on how these anthropogenic influences impact the fjord's marine habitat (Manzetti and Stenersen, 2010). Numerous studies in recent years have been conducted primarily on the microbial community (Poremba and Jeskulke, 1995;Storesund et al., 2017) or the influence of phytodetrital pulses on the macrofaunal community (Witte et al., 2003). Despite the accessibility of Sognefjord, the epibenthic megafauna community has been poorly studied, especially in recent years (Bernd, 1993). Bernd (1993) found that the central fjord and shallower adjacent side-fjord, Høyangsfjord, were primarily dominated by softbottom dwelling burrowing decapods (Munida sarsi and Munida tenuimana), holothurians (Bathyplotes natans and Parastichopus tremulus), and anthozoans. Høyangsfjord was found to have a lower megafaunal density compared to the main fjord. However, besides depth and general observations of substrate type, no thorough analysis was conducted on the influence of abiotic variables on the megafaunal community, especially over a horizontal or vertical gradient. In recent years several expeditions by the University of Bergen and the Institute for Marine Research in Norway, have taken place in the main-as well as the sidefjords, and knowledge about the fauna and benthic communities in Sognefjord is expected to increase rapidly when new data are made available over the coming years (Buhl-Mortensen et al., 2017;H. Glenner, personal communication).
In the present study, we used visual data collected with the remotely operated vehicle (ROV) AEgir 6000 to investigate the benthic megafaunal community diversity and distribution in Sognefjord and the influence of abiotic variables on the community structure at two locations in the fjord. It is expected that the communities would be less diverse and dense with increasing distance from the sill due to the extremely stable conditions identified by Storesund et al. (2017) (see the section "Study Area" for site description). For this study, we had three main objectives: (1) characterize the benthic community structure based on their depth and distance from the sill, (2) detect any response in community structure and diversity to changes in water mass characteristics above and below the sill, and (3) examine the relationships between environmental conditions and community composition and diversity (using generalized linear models, or GLMs).

Study Area
Sognefjord is located on the western Norwegian coast, extends to about 205 km and has a maximum water depth of 1308 m (Poremba and Jeskulke, 1995;Storesund et al., 2017). The fjord has a 3-layer water column structure (Svendsen, 2006;Storesund et al., 2017): the top layer is brackish water (salinity ≤ 33 psu), formed by the mixture of freshwater runoff and seawater, and moves out of the fjord; the intermediate layer (salinity between 33 and 35 psu) is well-oxygenated, owing to exchanges with the Norwegian Coastal Current above the sill depth, and hosts compensatory flows that may be in-fjord or out-fjord; and the bottom, or basin, water below the sill depth (salinity > 35 psu after deep-water renewal) originates from Atlantic water and becomes gradually less dense between renewal events because of diffusion and mixing with the layer above, driven by tidal forcing. Whilst the shallow sill (170 m) at the mouth of the fjord allows for some exchange of water between the fjord system and the adjacent coastal water, the mixing is fairly reduced and strong stratification occurs because of the influx of terrestrial runoff (Storesund et al., 2017). Renewal of the basin water occurs approximately every 8 years (Buhl-Mortensen et al., 2020).
Conditions within the fjord basin below sill depth are fairly stable and homogeneous. Water temperature is consistently around 7.4 • C and the salinity is generally greater than 35.0 psu (Witte et al., 2003;Storesund et al., 2017). Oxygen concentrations have been found to decrease with increasing distance from the sill and with depth, where Storesund et al. (2017) found the inner fjord and lower basin water to have oxygen concentrations below 4.5 mL L −1 in May and November 2012.

Data Collection
Two video transects were performed with the work-class ROV AEgir 6000 in Sognefjord in July 2017 during a SponGES cruise with the R.V. G.O. Sars (Figure 1). Dive 1 (D1) was conducted within the central fjord (1155 -85 m; 61 • 6 N, 6 • 39 E) and Dive 2 (D2) was located near the sill (1230 -55 m; 61 • 3 N, 5 • 22 E), respectively. D1 and D2 had approximate transect lengths of 2.43 km and 1.39 km, and were located approximately 79 km and 17 km away from the sill, respectively. The approximate fjord side slopes for D1 and D2 are 43 • and 44 • inclined to the horizontal plane. Physical samples were collected by the ROV during the transects to help confirm identifications of fauna.
One ship-based conductivity-temperature-depth (CTD) sensor package cast was conducted for each transect to profile temperature ( • C), salinity (psu), dissolved oxygen (mL L −1 ) and chlorophyll a concentration (µg L −1 ) throughout the water column. The bottom depths of the CTD casts corresponding to D1 and D2 were 1017 m and 1223 m, respectively. In addition, water samples were collected (using a rosette water sampler, on which the CTD package was mounted) from the different water masses for nutrient analysis. From the video footage, the percent cover of exposed hard substrate and soft sediment was estimated for each image analyzed.

Nutrient Content Analysis
Seawater samples for the analysis of inorganic dissolved nutrients (Si,PO 4 ,NH 4 , NO 3 , and NO 2 ) were collected with the CTDrosette from selected depths. Sample depths were selected based on the profile of the CTD downcast, whereby samples were collected from five different depths. Subsamples were collected in 50 mL syringes, which were rinsed three times with water from the niskin bottles of the CTD rosette before being filled. After sampling on deck, samples were filtered through 0.2 µm filters and instantly sub-sampled into two vials, one of which was used for samples of ortho-phosphate (PO 4 ), ammonium (NH 4 ), nitrate (NO 3 ), and nitrite (NO 2 ), stored at −20 • C and the other for silicate (Si), stored at 4 • C. Nutrients were analyzed at NIOZ with a QuAAtro Gas Segmented Continuous Flow Analyzer. Measurements were made simultaneously on four channels for PO 4 (Murphy and Riley, 1962), NH 4 (Helder and De Vries, 1979), NO 2 , and NO 3 (Grasshoff et al., 1983). Si was measured during a separate run (Strickland and Parsons, 1968). All measurements were calibrated with standards diluted in low nutrient seawater. For a detailed description of the sampling procedure we refer to Roberts et al. (2018).

Video Annotation
Still images were extracted from the videos approximately every 30 s to ensure there was no spatial overlap between images during analysis. Due to fluctuations in ROV altitude and changes in turbidity or topography, some areas of the transects were not suitable for analysis. Images were excluded if they contained any of the following characteristics: (1) image area obscured by part of the ROV or suspended material, (2) ROV was collecting physical specimens, (3) ROV was too far from the substrate (>10 m), (4) camera angle was not seabed-facing, (5) image contained poor light visibility, (6) image was blurred, and (7) overlapping image area. Videos were rescanned during image annotation to help identify individuals that were difficult to decipher in the stills alone, and in some cases, new stills were extracted if a more suitable image was present ±5 s from the original still. This was mostly the case for images that were blurred or contained poor light visibility and more suitable images of the same area were available within 5 s of the original image extracted. Parallel lasers of known separation (which project spots onto the seabed in order to determine image scale) were only active for the first hour of D1. It was therefore not possible to determine area (m 2 ) and density (individuals m −2 ) for the survey and megafauna were enumerated based on individuals per image.
ImageJ (version 1.52) was used to annotate the extracted imagery. A virtual grid with 496 cells was overlaid on each image in ImageJ. The grid size was selected because the cells completely overlapped the images. Percent cover of substrate type was estimated by counting the number of grid cells that contained that particular substrate type, then calculating the percentage of the total number of grid cells represented by that value. The grid resolution was selected such that the cells were small enough to minimize the number containing multiple substrate types as a proportion of total number of cells (and thereby increase the precision in the percent cover estimates) and large enough that the gridlines were not so dense as to obscure the fauna present in the image. In cases where cells contained multiple substrate types, the cell was characterized as the substrate that covered more than half of the cell. All epibenthic megafauna individuals that were easily visible within the imagery were counted and identified to the lowest taxonomic level. Due to fluctuations in altitude during video transects, some fauna had to be identified based on gross morphology (e.g., white sponge 1, yellow sponge 2, etc.). Taxa were classified as rare if represented by three or less individuals.

Data Preparation
Based on the description of the Sognefjord water mass structure by Storesund et al. (2017), the biotic and abiotic data were assigned to three depth zones to identify any depth-related changes in the benthic community structure: "Above Sill" (≤170 m), "Intermediate" (>170 to ≤300 m), and "Basin" (>300 m). No data from the surface brackish top-layer (1-10 m) was included. Due to differences in sampling frequency between the CTD profiling and image annotation, the abiotic variables were interpolated at 10 m depth intervals in Rstudio (version 1.2.5; RStudio Team, 2019). Nutrient composition was sampled with low frequency and was thereby excluded from statistical analysis. The megafauna abundances were summed into 10 m depth intervals. To account for missing data in the biotic dataset, certain depth intervals were removed from the abiotic dataset prior to statistical analysis. All multivariate statistical analysis was conducted in Primer (version 7) unless otherwise specified.

Environmental Variables
The environmental variables of D1 (central) and D2 (nearsill) were plotted against each other (e.g., temperatures from D1 vs. temperatures, at corresponding depths, from D2) to identify any notable differences in abiotic conditions between dives (Supplementary Figure S1). Furthermore, a correlation matrix was generated in Rstudio with package "corrplot" (version 0.84; Wei and Simko, 2017) to identify which abiotic variables were correlated (Supplementary Figure S2). Temperature had a strong negative correlation with salinity and positive correlation with dissolved oxygen (rho > 0.9), and was thereby dropped from further analysis. While depth was significantly correlated with the majority of the abiotic variables (p-value < 0.05) and had a moderately strong positive correlation with salinity and negative correlation with dissolved oxygen (rho > 0.7), it was selected to remain since it often acts as a proxy for other abiotic variables that were not measured in the present study. The remaining variables were normalized in Primer prior to multivariate analysis due to the different units used for each variable. A principal component analysis (PCA) of the selected abiotic variables was used to examine the environmental conditions within each depth zone per dive.

Community Composition and Diversity
Due to the lack of lasers throughout most of the dives, the abundances were converted to presence-absence data. A Sørensen similarity matrix between the two dives was calculated on the presence-absence dataset. Non-metric multidimensional scaling (nMDS) plots were constructed for each dive to identify differences between the community structure within each depth zone. An analysis of similarities (ANOSIM) was generated to identify significant differences between the depth zones and dives. SIMPER was used to determine which taxa were considered indicator organisms for each depth zone.
Diversity indices such as species richness and Shannon-Wiener diversity were calculated from the megafauna presenceabsence data for the depth zones in both dives to compare the changes in species richness and diversity over the vertical and horizontal gradients. In SPSS (version 25), a Levene's test of homogeneity of variances was used to determine if the diversity indices were homogeneous prior to running a one-way analysis of variance (ANOVA) to determine if there were significant differences for each index. A Tukey honestly significant difference (HSD) post hoc test on the diversity indices was used to identify which depth zones were significantly different for each dive.
To examine which environmental variables best explained the variance in species richness and Shannon-Wiener diversity indices, GLMs were generated in Rstudio. The residual deviance was larger than the residual degrees of freedom, therefore a quasi-Poisson error was fitted to the GLMs to account for overdispersion (Zuur et al., 2009). Depth, salinity, dissolved oxygen, chlorophyll a, percent cover of exposed hard substrate and soft sediment were included in the GLMs.

Environmental Conditions
In general, the water was slightly warmer, more saline, and more oxygenated in the water column below the sill at D2 (near-sill), relative to the same depths in D1 (central). The top layer of water was made up of warm brackish water (Figure 2), and a sub-surface chlorophyll maximum occurred at the halocline between the brackish surface water and intermediate water (approximately at 20 to 30 m). At approximately 80 to 100 m, dissolved oxygen decreased. Dissolved oxygen levels recovered in the intermediate water just below the sill depth. There was a gradual transition between the intermediate water mass and basin water until about 300 m, where a drop in temperature and dissolved oxygen levels occurred. The basin water was fairly homogeneous and there was not much difference in the water properties between the dives, where temperature was around 7.5 • C, salinity at 35.06 psu, and dissolved oxygen around 4.2 mL L −1 . D1 was frequently covered in soft sediment with exposed hard substrate patches throughout the transect, whereas D2 had distinct regions of solely soft sediment (e.g., bottom of fjord basin) and exposed hard substrate (e.g., along slopes and cliffs). Within the Intermediate Zone, small boulders and rocks became more frequent, and the sediment type became more coarse.

Dive
Depth waters at both sites. Nutrient concentrations did not show large differences between the two sites. In D1, from the suspended particulates observed in the video footage, predominately vertical settling appeared to occur. In these regions, sessile fauna and vertical rock walls were covered with a fine layer of particulate matter. However, D2 had higher observed turbidity throughout the dive compared to D1, with the settling of particulate matter apparent and some evidence of horizontal flow based on the position of the feeding apparatus of filter feeders and particulate direction changes.
In the PCA ordination, the first two principal components (PC) represented approximately 76% of the environmental variability within the two dives combined (47.1% and 28.6% for PC1 and PC2, respectively) (Figure 3). The PCA ordination showed that the images within the Basin Zones of both D1 and D2 were clearly distinguishable from the Intermediate and Above Sill Zones along the axes of PC1 and PC2. Dissolved oxygen (Eigenvector = 0.503), salinity (Eigenvector = −0.498), and depth (Eigenvector = −0.490) had the strongest influence on PC1. Percent cover of soft sediment (Eigenvector = −0.623) and exposed hard substrate (Eigenvector = 0.613) most strongly influenced PC2.

Sognefjord Megafauna Composition
In total, D1 had 79 taxa with a total of 11557 individuals and D2 had 89 taxa with a total of 10615 individuals ( Table 2). After rare taxa were excluded from the dataset, there was a total of 22105 individuals and 72 taxa identified within 511 images, where D1 had 11528 individuals and 57 taxa and D2 had 10577 individuals and 63 taxa. Porifera made up the majority of the taxa (24), followed by Echinodermata (16), and Cnidaria (10)

Community Trends and Observations
Within the Basin Zone, exposed hard substrate was often covered with polychaete tubes, the bivalve A. excavata, cnidarians, and encrusting sponges (Figure 4). Acesta excavata was present in high densities when observed on vertical rocky walls, often with juveniles and dense accumulations of encrusting polychaetes nearby. The octocoral Anthomastus grandiflorus (Verril, 1878) (Figure 4b) and large glass sponges Asconema aff. foliatum (Fristedt, 1887) (Figure 4e) occurred only on vertical rock walls. In D1, A. aff. foliatum was covered by a fine layer of suspended particulate matter (Figure 4e).
In regions with mixed substrate types (e.g., exposed hard substrate and soft sediment), Psolus squamatus was often observed concentrated at breaks in the slope or protruding surfaces. Phakellia spp., Phakellia ventilabrum (Linnaeus, 1767), and Axinella infundibuliformis (Linnaeus, 1759) were commonly positioned along slopes, aligned with the direction of observed horizontal particle flow (Figure 4m). Large anemones like Bolocera tuediae (Johnston, 1832) were observed residing on exposed hard substrate walls with soft sediment surrounding the substrate.
In the Intermediate and Above Sill Zone of D2, it should also be noted that there was a sudden occurrence of Echinoidea 1 in large quantities (max. 859 individuals per image) from 240 m to 60 m water depth (Figure 4u). For D2, the Above Sill Zone had a higher proportion of echinoderms present (Figures 4v,y). Numerous fish species, including Chimaera monstrosa (Linnaeus, 1758), Sebastes viviparus (Krøyer, 1845) (Figure 4v), and Coryphaenoides rupestris (Gunnerus, 1765), were observed within the upper regions of D1 and D2. Anthropogenic waste was found in the D1 Above Sill Zone (Figure 4w).

Community Structure
The community composition of the non-rare taxa showed significant differences between the depth zones and dives (ANOSIM Global R: 0.261, p = 0.001) (Figure 5). For D1, the nMDS plots and pairwise ANOSIM test indicated that the Basin and Intermediate Zones shared a similar community composition (p > 0.05). All other zones showed significant differences in community composition.
The SIMPER analysis revealed that for both dives, all zones except for D2-Above Sill had P. squamatus, Hymedesmia sp., and White Encrusting Sponge 1 within the top five contributing taxa ( Table 3). Taxa from the Echinodermata were more common in D2 than D1. For both dives, the Basin and Above Sill Total percentage (%) is the percentage of total abundance for each respective taxon. Asterisks denotes morphotaxa that were only observed within one dive, where * = Dive 1 and ** = Dive 2. Highlighted taxa are rare (<3 individuals) and dropped from statistical analyses. Figure 4 plate reference indicates which Figure 4 image plate that particular taxa can be found within. Bolded taxa are rare (<3 individuals) and dropped from statistical analyses.
Frontiers in Marine Science | www.frontiersin.org FIGURE 4 | Example of taxa and megafauna communities observed in Sognefjord, Norway. Panels (a-g) represents Basin Zone taxa on exposed hard substrate. Panels (h-m) displays Basin taxa in environments with both hard and soft substrate present. Panels (n-r) shows taxa common in soft sediment habitats. Panels (s-u) displays taxa within the Intermediate Zone. Panels (v-y) shows the Above Sill Zone. Refer to Table 2 for the taxa identifications.
FIGURE 5 | Non-metric MDS ordination plots visualizing the Sørensen resemblances of the megafauna presence-absence data at different depth zones between Dive 1 (stress = 0.23; left) and Dive 2 (stress = 0.19; right). The distance between the points relates to the similarity of community composition at each depth, whereby the closer the points the more similar the community.

Diversity of Sognefjord's Megafauna Community
The diversity indices for the D2 (near-sill) zones were consistently lower than those of the corresponding zones in D1 (central) (Figure 6). For both dives, the diversity indices for the Above Sill Zones were lower than those of the deeper zones. The diversity indices all passed the Levene's test of homogeneity (p > 0.05), and the one-way ANOVA indicated that there were significant differences (p < 0.001) in the indices for D1 and D2 and their respective depth zones. The species richness and diversity were statistically significantly different between the D1 and D2 respective Basin Zones (Tukey HSD: p < 0.01). There was a significant difference in species richness between the Basin and Above Sill Zones for D1 (Tukey HSD: p = 0.03) and trend toward significant difference between the Basin and Above Sill Zones for D2 (Tukey HSD: p = 0.069). For the Shannon-Wiener diversity index, there was a significant difference between the Basin and Above Sill Zones for D2 (Tukey HSD: p = 0.005), and a significant difference between the Intermediate and Above Sill Zones (Tukey HSD: p = 0.047).

Environmental Influence on the Megafauna Community
Salinity and dissolved oxygen were the most influential variables on the diversity indices when considered separately, as revealed by GLMs (Table 4 and Supplementary Figure S3), and depth had little influence alone. For species richness, the combination of depth and dissolved oxygen explained 10.85% of the deviance within the dataset. For diversity, the combination of depth, salinity, and dissolved oxygen explained 19.01% of the deviance in the dataset.

DISCUSSION
The present study provides a more recent overview of the Sognefjord megabenthic community composition than Bernd (1993) and focuses on the abiotic conditions more than was recently published by Buhl-Mortensen et al. (2017, 2020. The observations were similar to the findings from Bernd, where the deeper regions were characterized by sponges, holothurians and Munida tenuimana; however, in addition to sponges and holothurians, the shallower regions, particularly above the sill, had a higher abundance of urchins and anemones present. Sognefjord shares some of the same fauna elements found in Hardangerfjord (Buhl-Mortensen and Buhl-Mortensen, 2014;Buhl-Mortensen et al., 2020). Both fjords are dominated by Munida sp., Parastichopus tremulus, Psolus squamatus, and Phakellia species. However, as stated by Buhl-Mortensen et al. (2020), and as was observed in this study, many of the taxa observed in Sognefjord are not present in Hardangerfjord.
As is common for fjord systems, many of the taxa observed are continental slope or deep-water species. For example, Anthomastus grandiflorus (Verril, 1878), which is considered a deep-sea species with a distribution of 457-1760 m (Molodtsova et al., 2008), was observed in small clusters on vertical rock walls at depths below 540 m, and Kophobelemnon stelliferum and Virgularia mirabilis (Müller, 1776) were both observed in low quantities in soft bottom regions at depths below 630 and 500 m, respectively. Coryphaenoides rupestris was observed in the Basin Zone of D2, a deep-water fish that has been found to have isolated subpopulations within Sognefjord (Delaval et al., 2018). A peculiar finding was the presence of very large specimens (up to 140 cm long) of the hexactinellid (glass) sponge Asconema aff. foliatum on vertical cliffs below 800 m depth at both dive sites, representing a group of species normally confined to deep and cold waters on the outer shelf off Northern Norway, or along the Reykjanes or Arctic Mid-Ocean Ridges (e.g., Tabachnick and Menshenina, 2007;Maldonado et al., 2016;Roberts et al., 2018). Bolded taxa have the highest similarity (SIM)/standard deviation values (< 1.5). Asterisks (*) indicate taxa that were present in the respective depth zone for both dives.
Now we address each of this study's objectives in turn.

Community Patterns With Depth and Distance From the Sill
In general, much of the same taxa composition was observed in both dives for depth zones below the sill depth. The largest difference in megabenthic community composition was found between the deepest and shallowest zones for both dives, and similar trends have been observed in other surveys (Starmans et al., 1999;Sswat et al., 2015;Molina et al., 2019). In the present study, the Basin Zone was characterized more by sessile fauna (e.g., P. squamatus, Acesta excavata, Hymenodiscus coronata, and sponges) and M. tenuimana, whereas the Above Sill Zone was more dominated by echinoderms and anemones. Contrary to numerous fjord and shelf-based studies (Buhl-Mortensen and Høisaeter, 1993;Holte et al., 2004;Webb et al., 2009;Włodarska-Kowalczuk et al., 2012;Sswat et al., 2015;Gasbarro et al., 2018;Molina et al., 2019), we find that communities at the shallowest depths (Above Sill Zone) and in closer proximity to the sill (D2) have the lowest number of species and diversity. Buhl-Mortensen et al. (2017, 2020 observed a similar trend in relation to the proximity to sill, where the outer region in Sognefjord had lower species richness compared to the middle region (Buhl-Mortensen et al., 2017, 2020. However, Buhl-Mortensen et al. (2020) observed a trend of decreasing species richness with increasing depth, which was not observed in the present study. The trends observed in the present study are more consistent with shelf community patterns observed by Starmans et al. (1999), where shallower regions contained a lower number of highly abundant taxa than the deeper stations and diversity increased with increasing water depth. This reduction in species richness and diversity in D2 (near-sill) and the areas above the sill could be driven by changes in water mass characteristics or increased physical stress on the Percentage (%) explained is the percentage of null deviance in the data explained by the model.

Influence of Water Mass Properties and Sill Depth
The basin waters of both dives were fairly homogeneous (Storesund et al., 2017) and likely contributed to the homogeneity observed in the species composition in the deeper regions. Water mass properties (e.g., temperature, salinity, dissolved oxygen) play a significant role in megabenthic community composition (Buhl-Mortensen and Høisaeter, 1993;Williams et al., 2010: Meyer et al., 2015, which appears to be the case for Sognefjord as well. The changes in species composition appear to gradually occur around the transition between the basin and intermediate water masses, which is at approximately 300 m (Storesund et al., 2017), and more clearly near the sill depth. Studies have shown that sills affect water mass dynamics in ways that are critical to the structuring of benthic communities (Strømgren, 1970;Rüggeberg et al., 2011). As Buhl-Mortensen and Høisaeter (1993) stated, the environment in fjord basins is influenced by the sill depth. With shallow sills, organic matter becomes trapped within the inner fjord below the sill depth and is not flushed out readily by the adjacent coastal water (Klitgaard-Kristensen and Buhl-Mortensen, 1999). As such, it is possible that organic input from renewal events and terrestrial sources (e.g., rivers, runoff, snowmelt, etc.) accumulates and has longer residence times in fjord basins (relative to shallower waters), providing food and nutrients to the benthic fauna below the sill depth. In a recent study of the Sognefjord by Buhl-Mortensen et al. (2020), the authors observed continuous detritus cover on sloping bedrocks at depths greater than 400 m and terrestrial organic material mixed in with the basin's soft sediment. The observed higher species richness and presence of deposit-feeding holothurians (Bathyplotes natans and Mesothuria intestinalis) and suspensionfeeding Hymenodiscus coronata in the Basin Zone's soft bottom regions indicate availability of organic matter to the basin floor (Roberts and Moore, 1997;Flach et al., 1998;Amaro et al., 2015).
The vertical-falling particulate matter along the rocky walls observed in D1 is likely an important food source for many of the filter-and suspension-feeders (e.g., encrusting sponges, Asconema aff. foliatum, encrusting polychaetes, and Acesta exacta) residing on the vertical rock walls or under overhangs. Areas of flow acceleration owing to irregular topography (e.g., ridges, peaks, and other elevated substrate) experience increased particle fluxes and are also likely important for suspension feeders (e.g., Psolus squamatus, Phakellia spp., Phakellia ventilabrum, and Axinella infundibuliformis) (Flach et al., 1998;Buhl-Mortensen et al., 2020). As is common in fjord environments, it is likely that the quality of food is lower in the basin and inner fjord compared to regions nearer to the sill (Klitgaard-Kristensen and Buhl-Mortensen, 1999). However, the higher species richness and presence of suspension-and deposit-feeders within the basin suggests the fauna may be adapted to the low quality of food, or that this is compensated by the stability of the basin environment. It is clear that a more rigorous study should be conducted to quantify and assess the quality of the organic matter supplied to the basin communities.

Environmental Dynamics
Depth, salinity, and dissolved oxygen were highlighted as important variables for the diversity indices. Depth acts as a proxy for other factors and it is likely that parameters which were not accounted for in the present study (e.g., food availability, particulate organic matter, localized hydrodynamics, pollution) are also influencing the patterns observed (Jones et al., 2007;Webb et al., 2009;Williams et al., 2010). Dissolved oxygen and percentage cover of substrate type varied most between the two dives, both of which are known to be critical for many benthic habitats (Holte et al., 2005;Williams et al., 2010;Sswat et al., 2015). The availability of hard substrate is important for sessile invertebrates (Williams et al., 2010;Buhl-Mortensen et al., 2012), and in this study, regions with exposed hard substrate were often covered with sponges, serpulid worm tubes, bivalves, and holothurians, similar to observations made by Gasbarro et al. (2018). However, for the Sognefjord megafauna community, there was not much difference in diversity and species richness between percent cover of hard substrate, soft bottom or mixed substrates. Therefore, it is possible that other factors like environmental dynamics or food availability is driving the patterns observed.
The megafauna communities at the mouth of the fjord (D2) showed lower diversity and species richness compared to central fjord (D1) communities. The central fjord environment is more stable than that of the fjord mouth, which is subjected to greater temporal variability (Storesund et al., 2017) due to exchange with the coastal ocean. Differences between the two dives are likely to be largely a result of horizontal environmental gradients along the fjord set up by biogeochemical and physical processes. For example, dissolved oxygen concentrations at the interface between the basin and intermediate water and at the sill depth differed slightly between dives (Supplementary Figure S1), the water being more oxygenated toward the sill (D2) because of the influence of coastal water. Dissolved oxygen concentrations at these depths are likely diminished with distance up-fjord by diffusion to and entrainment of vertically adjacent, less oxygenated waters and by the cumulative effects of (bacterial) respiration with distance from the sill (Storesund et al., 2017).
Areas with high environmental disturbance are characterized by an increase in mobile species, decrease in sessile fauna, and overall lower diversity (Włodarska-Kowalczuk et al., 2005;Jones et al., 2006;Webb et al., 2009;Hughes et al., 2010;Włodarska-Kowalczuk et al., 2012), as was observed in D2 and regions above the sill depth. The higher turbidity observed in D2 may have impacted the fjord benthic communities. Sessile suspensionfeeding invertebrates are at a risk of smothering in areas with high quantities of suspended material in the water column (Jones et al., 2006;Kutti et al., 2015;Meyer et al., 2015). Fauna that are not limited by such conditions can persist (Rygg, 1985;Włodarska-Kowalczuk et al., 2005Gasbarro et al., 2018), sometimes in high abundances, which could contribute to the increased abundance of echinoderms and reduced occurrences of sponges and sessile holothurians in the shallower regions of the fjord. Buhl-Mortensen et al. (2020) also noted that the shallower and silled regions of the fjord have relatively strong currents, whereas, the bottom currents in the basin were weak. This supports the general picture of gradients in environmental variability and stress within the fjord.

Future Implications
The environmental conditions in Sognefjord are affected by anthropogenic influences, such as cruise ships, fish farms, hydroelectric stations, and pollution (Manzetti and Stenersen, 2010). There is limited information concerning how such influences impact the Sognefjord community, though fish stocks have seen a considerable reduction (see Manzetti and Stenersen, 2010) and the shellfish community showed increased diarrhetic shellfish poisoning toxins with increased distance into the fjord (Ramstad et al., 2001).
A study by Rygg (1985) found that fjords with high pollutant concentrations were characterized by opportunistic species. Additional organic input from anthropogenic sources like fish farming or nutrient runoff may lead to hypoxic conditions in the fjord basin (Levin et al., 2009;Johansen et al., 2018). Johansen et al. (2018) predicted that increased organic matter within Norwegian fjord basins will lead to a dominance of deposit feeders, while the presence of suspension feeders will decline. Similar to the findings of Rygg (1985), Johansen et al. (2018) also found a shift in the community structure toward opportunistic species as a result of oxygen depletion and increased temperatures. Coastal water temperatures have been rising (Aure, 2016), which has led to increased temperatures within fjord basins (Johansen et al., 2018), resulting in changes in stratification and reduced oxygen supply. Long-time series temperature and organic input data are not readily available for Sognefjord, but if its environmental conditions follow the trajectories of other Norwegian fjords it is possible to predict a similar shift toward more opportunistic species. The present study does not include any temporal replication and the impact and future implications of anthropogenic-derived environmental change on the system is largely unknown and requires more research.

CONCLUSION
This study provides a recent overview of Sognefjord's megabenthic community near the sill of the fjord and the central fjord. Megafauna community composition was homogeneous within the fjord basin; however, species richness and diversity declined with proximity to the sill and with decreasing water depth, particularly at the boundary between basin and intermediate water and at the sill depth. The fjord basin was characterized by Psolus squamatus, Munida tenuimana, Phakellia sp., Acesta excavata, and encrusting sponges. At shallower depths, the fjord was dominated by echinoderms, particularly in the dive closest to the sill. It is clear that more research is needed to understand the influence of shallow sills and water mass structure on fjord communities, as this study shows these features are important to Sognefjord's megabenthic communities. The clear stratification occurring between the basin water and intermediate water within Sognefjord would make it well suited for future surveys designed to monitor a wider range of environmental conditions or to understand future scenarios with stratification changes or deoxygenation.

DATA AVAILABILITY STATEMENT
The datasets generated for this study can be found at PANGAEA: https://doi.pangaea.de/10.1594/PANGAEA.914801.

AUTHOR CONTRIBUTIONS
HR, ER, and FM collected the video footage, CTD casts, and nutrient data for the study. HM annotated the video footage. HR provided fauna identifications. ER interpreted the in situ abiotic conditions and calculated the slope and distance from the sill for the transects. HM performed the statistical analysis. FM provided the sampling design, analysis and description for the nutrient contents. HM wrote the manuscript. All the authors read, edited, and approved the final manuscript.

FUNDING
The work leading to this publication has received funding from the European Union's Horizon 2020 Research and Innovation Programme through the SponGES project (Grant Agreement No. 679849). This document reflects only the authors' view and the Executive Agency for Small and Medium-sized Enterprises (EASME) is not responsible for any use that may be made of the information it contains. FM is supported by the Innovational Research Incentives Scheme of the Netherlands Organisation for Scientific Research (NWOVIDI Grant No. 016.161.360).

ACKNOWLEDGMENTS
The video footage and CTD casts were collected in 2017 on the RV G.O. Sars during a SponGES cruise, therefore, the crew of the RV G.O. Sars and the ROV AEgir 6000 as well as participating SponGES team are thanked for their contribution to this project. EMODnet Bathymetry Consortium (2018) is acknowledged for the use of high-resolution bathymetry map for Sognefjord. The work presented here is dedicated to the memory of our friend and mentor HR, who spent his life improving global understanding of deep-sea sponges.