Habitat Connectivity of Fish in Temperate Shallow-Water Seascapes

Movements of organisms comprise a fundamental aspect of coastal habitat connectivity. Determining the distribution and co-existence of habitat specialists and generalists in shallow-water seascapes ...


INTRODUCTION
Shallow-water habitats are vital for healthy coastal areas and contribute globally to fisheries productivity and maintenance of biodiversity (Stål et al., 2008;Ramos et al., 2015;Nordlund et al., 2017). Fundamental to ecology is an understanding of assemblage composition and interactions. Recent efforts have shifted toward examining connectivity between habitats (and habitat patches) within the seascape mosaic, in an effort to broaden applicability to marine spatial planning and comprehensive management, (e.g., Gillanders et al., 2003;Unsworth et al., 2009;Yeager et al., 2011Yeager et al., , 2016Caldwell and Gergel, 2013).
Fish are capable of utilizing multiple habitats, and are therefore influenced by the complexity and structure of the seascape as a whole (Pittman et al., 2007;Grober-Dunsmore et al., 2009;Staveley et al., 2017). A tremendous effort has been made toward describing and comparing distribution patterns of fish in different shallow-water habitats (e.g., Dean et al., 2000;Mathieson et al., 2000;Able, 2005;Dorenbosch et al., 2007;Reiss et al., 2010), such as the significant role of seagrass meadows as nursery grounds (reviewed by Heck et al., 2003). However, in temperate waters the strength of habitat linkages and the relative importance of multiple habitats from a seascape perspective (considering ecological functions in relation to spatial patterns) has, as yet, been given sparse attention (see Perry et al., 2017;Staveley et al., 2017). Fish connect habitats by the exchange of biomass and energy via ontogenetic migrations, larval dispersal and daily movement. Many coastal marine fish species use different habitats during different life stages (Gillanders et al., 2003). Indeed, many vulnerable species, such as Gadus morhua, use key habitats (e.g., seagrass meadows) during at least part of their life cycle (Musick et al., 2000;Orth et al., 2006). Studies have shown that fish species richness and abundance are higher for juveniles and subadults in seagrass meadows compared to larger adults (Gullström et al., 2008;Bertelli and Unsworth, 2013). Additionally, research from tropical regions indicates that seagrass meadows close to other habitats, such as mangroves and coral reefs, have a positive influence on nursery species (i.e., fish species that use specific habitats exclusively during juvenile life stages), and are essential in shaping fish assemblage structure (Gullström et al., 2008;Berkström et al., 2013b). Given that mobile species use multiple coastal habitats, it is important to understand patterns and processes at a seascape level.
As human activity fragments, and even destroys, important coastal habitats (Baden et al., 2003(Baden et al., , 2012Pihl et al., 2006;Crook et al., 2015;Nagelkerken et al., 2015), defining critical areas and nursery grounds for fish is imperative. This is of importance given the fact that species possess varying degrees of mobility. Therefore, it can be expected that stationary species suffer more negative consequences than migrating species as a result of habitat fragmentation and destruction (Caldwell and Gergel, 2013). Maintaining the integrity, diversity, and connectivity of shallow-water habitats is important for healthy fish communities. Due to their dependence on these shallowwater biomes during various life-stages, commercially and economically profitable fisheries/fish assemblages will certainly benefit from the preservation of these habitats (Pihl et al., 2006;Stål et al., 2008;Baden et al., 2012;Bertelli and Unsworth, 2013;Ramos et al., 2015).
Research has found that nursery species prefer the clearer water in seagrass meadows compared to the turbid water found in unvegetated sandy areas . The authors argue that this pattern is likely a result of the protection provided by the structural complexity of seagrass meadows rather than the reduced visibility due to turbidity. In addition, the food supply in vegetated habitats is important for food webs within the shallow-water environment (Valentine and Heck, 1999;Baden et al., 2010). Heck et al. (2003) posit that the existence of structure (regardless of type) is critical for shallowwater habitats to function as nursery grounds. Therefore, habitat structure and linkages between areas of complex construction, rather than specific habitats themselves, may drive improved species abundances seen in vegetated areas (vs. unvegetated) (Orth et al., 1984;Hemminga and Duarte, 2000;Wennhage and Pihl, 2002;Pittman et al., 2004;Gullström et al., 2008Gullström et al., , 2011. Nagelkerken et al. (2013) proposed the "seascape nursery" as a conceptual model that defines a mosaic of coastal habitats as an interlinked entity. While some research in tropical regions has evaluated the seascape nursery concept and the strength of connectivity between shallow-water areas (e.g., Nagelkerken et al., 2000;Dorenbosch et al., 2007;Gullström et al., 2012;Berkström et al., 2013a,b), this has seldom been investigated in temperate regions. In temperate shallow coastal waters, fishhabitat associations have been studied in various biomes such as seagrass meadows and rocky bottoms (e.g., Wennhage and Pihl, 2002;Jackson et al., 2006;Pihl et al., 2007;Hutchinson et al., 2013) though the focus on habitat connectivity has often been overlooked.
The primary focus of this study was to examine the nature and strength of habitat usage and connectivity in shallow-water coastal seascapes by comparing fish assemblage composition in three different (but adjacent) key habitats: seagrass meadows, rocky bottoms covered by macroalgae, and unvegetated areas. This was done using field observations at multiple sites on the temperate Swedish Skagerrak coast. Using habitat preferences and the life stage of individuals, we hypothesized high levels of similarity in terms of species composition and distribution patterns of fish between adjoining habitats in shallow-water environments, thereby demonstrating connectivity between neighboring habitats in shallow coastal seascapes.

Study Location
The study took place on the Swedish Skagerrak coast (58 • 00 ′ N−58 • 5 ′ 6 ′ N, 11 • 00 ′ E−11 • 67 ′ E) in June 2014 (Figure 1). This area is characterized by low tidal fluctuations with an average of ∼30 cm daily, though it can oscillate as much as 2 m depending on the strength of the winds (Johannesson, 1989). The archipelago is a fjord-like system consisting of many rocky islands and islets. The area is a productive transitional zone connecting the oceanic waters of the North Sea with the low saline waters of the Kattegat region in the south, an area influenced by the brackish waters of the Baltic Sea. Within the study area, the surface water salinity varies from 15 to 25 (Baden et al., 2012), though it has been reported to be as high as 33 (Björk and Nordberg, 2003).
We examined six seascapes each with a diameter of 600 m, all of which had been previously mapped with the use of a drop video system (SeaDrop camera, SeaViewer Cameras, Florida, USA) . We use the term "seascape" following the definition by Pittman et al. (2011) stating that seascapes are "wholly or partially submerged marine landscapes." Each seascape contained the three distinct focal coastal habitats necessary for the current study: (i) seagrass (Zostera marina) meadows, (ii) rocky bottoms covered by macroalgae (rock-macroalgae), and (iii) unvegetated soft bottom areas (unvegetated). The minimum area considered a distinct habitat used for placing the cameras was approximately 850 m 2 . For each seascape, an unbaited stereo-video system was deployed in each focal habitat.

Camera Surveys and Assessment of Fish Assemblages
Using Remote Underwater Video (RUV) for marine research is a relatively new approach that has grown in popularity because (1) it is a non-extractive method, as opposed to more destructive sample methods employing fishing gear, and (2) it can mitigate problems faced by diver surveys (Underwater Visual Census-UVC), such as disturbance and expense (Langlois et al., 2010). RUV has great potential for studying movements of individuals and connectivity through the use of stereo video cameras. With stereo video, two cameras are mounted on a calibrated frame and then synchronized to record the same object simultaneously (Mallet and Pelletier, 2014). These recordings give the observer highly accurate three-dimensional images, while collecting data on both flora and fauna unobtrusively (Harvey et al., 2003). This ability to observe without disturbance makes this method a valuable tool for studying mobile species such as fish. While the remote underwater camera is a suitable method for studying fish communities, and is well designed for habitats such as rocky areas that are difficult to study with traditional net methods (Cappo et al., 2004;McIlwain et al., 2011), it also has its limitations.
Though RUV methods are able to capture fish behavior undisturbed, they are limited by water visibility; there must be suitable light for recording, and water clarity must be adequate for data collection. Furthermore, it is very difficult to determine if the video is recording a new fish entering the field of view or if it is the same fish previously recorded. For this reason, double counting is an inherent limitation of the method. In order to standardize the methodology, all fish that were out of the field of view for at least 3 s were considered a new fish and therefore counted. Though this is a limitation of the stereo-video method, it would limit all videos equally in the current study and therefore not influence the overall results.
Through calibration and the use of the EventMeasure (www. seagis.com.au) analysis software, data was collected on species identification and abundance information, as well as the length of individuals (Harvey et al., 2001(Harvey et al., , 2002(Harvey et al., , 2003. All fish were identified to the lowest taxonomic level possible, counted, and measured. Subsequently, abundances were calculated (m −2 ) based on fish activity. This was defined as the number of instances fish entered the field of view divided by the recording area, which therefore corrects for the variation in visibility among sites (Hammar et al., 2015;Perry et al., 2017). Field of view is defined by the maximum distance at which a recorded fish is identified for each specific film. To obtain length information fish must be measurable in both camera fields of view with their image overlapping in the two (stereo) images simultaneously. Hence, length measurements were difficult to obtain in highly vegetated habitats where certain parts of the fish could be obstructed. Additionally, cryptic species with little to no movement can be difficult to measure; their ability to camouflage themselves with their surroundings may obscure the edge of the fish in relation to the habitat. However, when fish length data were obtainable, adult and juvenile life stages were determined using their length at maturity (Froese and Pauly, 2015). When specific maturity data were unavailable, a method commonly used to determine life stage was employed, by which individuals 1/3 of the maximum recorded length (according to FishBase) are considered to be juveniles (Dorenbosch et al., 2004(Dorenbosch et al., , 2007Unsworth et al., 2009).
The camera system was deployed from a boat and placed in the center of a habitat aided by a snorkeler, angled to record each habitat properly (upward toward the rock-macroalgae area, down toward the unvegetated soft bottom and straight forward/slightly up in the seagrass meadow), ensuring that no vegetation directly obstructed the field of view. Two GoPro R HERO2 cameras mounted on a calibrated frame with 60 cm base separation were used for data collection. Specific calibration (Hammar et al., 2013) and additional methodological details  can be found in previous investigations. Once positioned, the cameras were left to record for 50 min of continuous filming, with an additional 5 min at the start of all videos to exclude disturbance effects. In total, 900 min of video were analyzed by the same observer. The camera was positioned at depths ranging from approximately 1.5 m to 3 m and the field of view ranged between approximately 0.5 and 3.5 m 2 (see Appendix A in Supplementary Material).
All videos were analyzed for fish community information within each of the three habitats. Additionally, species were grouped into guilds based on habitat preference in order to elucidate similarities and differences in fish assemblage structure between various coastal shallow-water habitats, to understand species connectivity between habitats. Guilds were selected based on information from FishBase (Froese and Pauly, 2015) and other sources (Elliott and DeWailly, 1995;Pihl and Wennhage, 2002;Pihl et al., 2006;Perry et al., 2017;Staveley et al., 2017). The guild groupings were as follows: stationary species (SS), shallowwater generalists (SWG), occasional shallow-water visitors (OSV) and juvenile migrants (JM). (1) Stationary species do not actively leave a specific habitat, though they may be found in various habitats and transported via water movement to different areas at early life-stage phases.
(2) Shallow-water generalists regularly move between coastal shallow-water habitats but typically not outside the coastal areas. (3) OSV utilize certain shallow-water habitats occasionally as well as other habitats, such as deepwater areas or fresh water rivers. (4) Juvenile migrants use shallow-water habitats in their early life stages before migrating to nearby habitats as adults (or subadults) but may also return occasionally as adults for feeding. Fish within the Pleuronectidae family were not included in the habitat preference guild grouping for this study. We were unable to identify the individuals to species level and different species have different habitat preferences.

Data Analysis
Differences in fish abundance by habitat, number of species and diversity (Shannon-Wiener diversity index) were analyzed by randomized block design analysis of variance (ANOVA) using R (v. 3.2.0). Prior to analyses, data were log 10 (x+1)-transformed in order to validate statistical assumptions. Pairwise post-hoc Tukey's tests were conducted to determine significant differences between habitats. Additionally, multivariate analyses were conducted on fish assemblage data using the Primer software (v. 6, Plymouth Marine Laboratory). Given the extreme variation in species abundances (because of some highly abundant species), presence/absence data were used for the Sorensen similarity measures. One-way analysis of similarity (ANOSIM) was carried out to test for differences in fish assemblage structure between habitats, and patterns were visualized using nonmetric multidimensional scaling (nMDS) ordinations. Moreover, similarity of percentages (SIMPER) analyses were conducted to determine which species were driving the differences in fish assemblage structure among habitats.
It should be noted that Atlantic herring, Clupea harengus, was removed from the analysis when abundance data were utilized, given that this species was found in extremely high numbers compared to all other species identified (or shown with/without C.harengus for comparative purposes). When only presence/absence data were considered C. harengus was included. The reason for removal was to avoid the extremely high abundance muting all other fish abundance results, however it was included in presence/absence calculations because we felt it important to show that it is occasionally a part of the shallowwater habitats in the area. Atlantic herring are a coastal pelagic fish exhibiting schooling behavior and are therefore expected in all types of coastal habitats, as well as in the open ocean, and when seen are found in very high densities (Reid et al., 1999).

Fish Assemblage Compositions
A total of 11,744 fish were recorded, including 15 taxa (13 to species level) within 9 families. Clupeidae was the most abundant family followed by Gobiidae and Labridae, while Gobiidae exhibited the highest species richness (Table 1). At species level, C. harengus was by far the most abundant species throughout all habitats, followed by Gobiusculus flavescens (Table 1). Total fish abundance was significantly higher in seagrass meadows and rock-macroalgae habitats, compared to unvegetated soft bottom areas (S-U p < 0.001; R-U p = 0.002), while seagrass meadows and rock-macroalgae habitat abundances did not differ from each other [(S-R p = 0.217) Block ANOVA; Habitat: F (2, 10) = 24.482, p < 0.001; Site: F (5, 10) = 2.644, p = 0.090] (Figure 2). Interestingly, this is the case despite the fact that the total area of the unvegetated habitat, on average over the six seascapes, was the largest, followed by the area of seagrass, while the smallest coverage area was that of the rock-macroalgae habitat (Figure 2). Although the abundance differed significantly between some of the habitats, neither the Shannon-Wiener diversity index (excluding the abundance of C. harengus) [Block ANOVA; Habitat: F (2, 10) = 0.602, p = 0.567; Site: F (5, 10) = 0.453, TABLE 1 | List of fish taxa identified from camera surveys in each of the three studied habitats (seagrass meadows, rocky bottoms covered by macroalgae and unvegetated soft bottom areas), their habitat preference guilds and the mean activity per m 2 (proxy for density) per 50 min sampling period for each fish taxon as well as the percent of each species compared to the total abundance of fish within each habitat.
The assemblage structure of fish (Figure 3) varied significantly among the three habitats (ANOSIM; Global R = 0.435, p < 0.001), with seagrass-and rock-macroalgae-habitats differing FIGURE 3 | Non-metric multidimensional scaling (nMDS) based on Bray-Curtis similarity coordinates for presence/absence of fish species separated into the three studied shallow-water habitats. Site numbers are labeled (compare with Figure 1); note that two unvegetated habitat sites overlap.
from the unvegetated areas (Seagrass vs. Unvegetated: R = 0.594, p < 0.01; Rock-macroalgae vs. Unvegetated: R = 0.635, p < 0.01), while the two vegetated habitats did not differ from each other (Seagrass vs. Rock-macroalgae: R = 0.026, p = 0.37). When evaluating the assemblage structure based on presence/absence data in order to control for the disproportionately large abundance of C. harengus, the SIMPER analysis showed that the taxon contributing most to the dissimilarities between the rock-macroalgae-and unvegetated habitats was Pomatoschistus spp., which was only present in unvegetated areas, whereas G. flavescens was found more often in the rock-macroalgae habitat ( Table 2). Gobius niger, on the other hand, was only found in unvegetated soft bottom habitats (Table 2). Similarly, these were also the same species driving statistical differences between seagrass meadows and unvegetated areas ( Table 3). Species such as C. harengus, A. minuta, G. flavescens, and C. rupestris were, however, observed in all three habitats.

HABITAT PREFERENCE GUILDS
When reviewing all camera data the habitat preference guild with the highest abundance was the OSV with approximately 90-98% of the total abundance in all three habitats when C. harengus was included. Excluding C. harengus, the OSV yielded an approximate abundance less than half a percent in both the unvegetated-and rock-macroalgae habitats. When omitting C. harengus, the seagrass-and rock-macroalgae habitats were dominated by the SS guild, while the SWG was only about 20% of the total abundance. However, the opposite was found in the unvegetated habitat as the SWG was most abundant followed by the SS. Taxa from the JM guild were only seen within the seagrass habitat (Figure 4).

Length and Life Stage of Fish
The average length of fish did not differ significantly among habitats [ANOVA; F (2, 5) = 2.065, p = 0.178], though the fish within the unvegetated areas were, on average, slightly larger than those in the seagrass, followed by the rock-macroalgae habitats (5.7 cm ± 0.9 SE, 5.3 cm ± 0.3 SE, 4.9 cm ± 0.2 SE, respectively). It must be noted that the number of measurable fish was not equal to the number of identified fish as not all fish are seen within both cameras' field of view (refer to Methods section for the discussion of methodological limitations). For the unvegetated soft bottom habitat, 2% (n = 108/6008) of the observed fish were measurable, 4% (n = 78/2048) within the seagrass meadow and 16% (n = 582/3688) from the rock-macroalgae habitat. While length, and therefore life stage information, was estimated on an individual basis, the limited number of measurable fish made analysis of averages necessary. From these measurable fish an evaluation of the proportion of juveniles was performed using ANOVA and it was shown that there was no significant difference between habitats [F (2, 5) = 0.634, p = 0.551], with approximately 75% juveniles comprising all three habitats (Figure 5).

DISCUSSION
The concept that the complexity and interconnectedness of habitats (rather than evaluating each specific habitat individually), is of importance for those organisms dependent thereon, has been gaining momentum (Heck et al., 2003). In this study, we tried to understand in more detail how fish assemblages within three adjacent shallow-water habitats overlap. We did this by assessing variation between habitat types, as well as explicitly focusing on habitat variation within different coastal seascapes, to see whether differences exist between habitats and/or seascapes. In support of our hypothesis, we found that there was high similarity in species composition between adjoining habitats. Specifically, we show that neither the number of species nor species diversity differed among the three shallowwater habitats, with many specimens of the same species found in all three habitat types implying a degree of connectivity between habitats. This indicates that the fish community is similar in adjoining habitats of temperate shallow-water seascapes, though structurally complex habitats (i.e., seagrass meadows and macroalgae beds) had significantly higher abundances of fish compared to unvegetated areas. Also, interestingly, the proportion of juveniles found within each habitat did not differ between the three habitats.  The mean fish abundance was highest within the seagrass habitat, a result initially appearing to contrast with those found previously on the Swedish west coast by Stål et al. (2007), where rocky areas contained higher abundances than soft bottom areas. However, the soft bottom areas examined in the study by Stål et al. (2007) included depths (6-10 m) not evaluated in the current study and where seagrass does not grow, and thus is the most probable cause for the contradictory findings. While seagrass beds had the highest fish abundance in our study, they did not differ significantly from the rock-macroalgae habitat. These findings emphasize the need to consider the shallow-water seascape as an integrally-linked coastal area important for species assemblages, rather than an area of isolated homogenous habitats important for particular species. This is consistent with recent results from the tropics showing that both seagrass beds and macroalgae habitats are important for juvenile fish (Tano et al., 2017).
Comparable fish communities within shallow-water habitats on the Swedish west coast have been demonstrated by Pihl and Wennhage (2002). However, the perspective of a continuous and interconnected shallow-water seascape has not previously been highlighted as in the current research. Here, we evaluated all habitats within the same seascape together (600 m diameter), whereas previous studies utilized random habitat sampling methods. In our design and analyses all three habitats were adjacent and therefore part of the same contiguous shallow-water seascape. Many of the species from the current study were found in at least two of the three shallow-water habitats examined, with a few species found in all three habitats. Clearly, the species within this study utilize multiple habitats within the shallowwater seascape, indicating that species may be moving between the adjacent habitats. However, some species, such as those from the stationary species guild, may be habitat-specific at the individual level as opposed to the species level. Similar findings of a substantial overlap of species within multiple coastal habitats have been found in different areas of the world Franco et al., 2006;Unsworth et al., 2009;Berkström et al., 2013b;Lilley and Unsworth, 2014), including temperate Swedish waters Wennhage and Pihl, 2002).

Habitat Connectivity
Our results strengthen the postulation that habitat connectivity is strong in temperate shallow-water seascapes. Different species, such as C. harengus, G. flavescens and C. rupestris, were observed within all three focal habitats, which made up three of the four habitat preference guilds assessed (OSV, SS, and SWG, respectively). These data, when combined with life-stage history information from the literature, show highly connected seascapes where all shallow-water habitats evaluated are important for these species as well as the juvenile migrant (JM) G. morhua (Figure 6). Clupea harengus is a species known for its schooling behavior utilizing coastal pelagic waters as juveniles as well as the open ocean during all life stages (Reid et al., 1999). As adults, specimens of C. harengus are found to make vertical migrations between deep and shallow waters The stationary species G. flavescens is ecologically associated with both macroalgae-and seagrass habitats as adults, while different parts of the water column are used at different life stages in response to predator avoidance (Folkestad, 2005;Froese and Pauly, 2015). Juveniles have been found in seagrass meadows  as the eggs are demersal; after hatching the larvae are typically seen in the sublittoral zone (Folkestad, 2005). While G. morhua was only seen in the seagrass habitat in the current study, previous studies have shown that juveniles may be associated with rockyand unvegetated areas as well . It is possible that our results are due to a limitation of the camera method, which could potentially be limited by species behavior; G. morhua can respond strongly to changes in light, and therefore are more active at dawn and dusk (Staveley unpublished data). Adults are known to perform vertical migrations from deep-to shallow-water areas and the eggs and larval stage have been seen in the pelagic zone (Petitgas et al., 2013). For C. rupestris, while seen in all habitats within the shallow-water seascape, literature suggests that juveniles and adults are typically associated with seagrass, and with rocky areas Staveley et al., 2017). Clearly, the species in Figure 6, among others in the current study, are associated with many habitats, both coastal and offshore, showing strong connectivity patterns between various habitats.
The species observed in multiple habitats represented the habitat preference guilds occasional shallow-water visitors (OSV), stationary species (SS), and shallow-water generalists (SWG), while juvenile migrants (JM) were identified only in seagrass meadows. These findings are similar to a study in tropical waters by Kendall et al. (2011), who found that when evaluating the seascape influence on fish mobility guilds, all guilds (residents, mobile fish and transient fish) were influenced by the surrounding seascape, including seagrass meadows, sand/soft bottom and hard bottom areas. In this study, the vegetated habitats consisted of more than 75% of specimens belonging to the SS guild, while half of the specimens in unvegetated areas represented the SWG guild. This may be attributable to the fact that stationary species are reliant on specific habitat to provide them with food and protection, and therefore structurally complex areas offer more opportunity to hide (e.g., while foraging) (Heck et al., 2003;Nagelkerken and van der Velde, 2004). In contrast, unvegetated areas may largely be used for temporary passage between nearby habitats or by cryptic, sand-colored species such as Pomatoschistus spp. and some bottom-dwelling flatfish (Froese and Pauly, 2015). It is also possible that the stationary species are found in the unvegetated area as a result of passive transport via water current movement.

Fish Assemblage
The SIMPER analysis indicated that similar species contributed to the variation in seagrass meadows and rock-macroalgae habitats compared to unvegetated areas. The differences in fish assemblage structure between the structurally complex habitats and the unvegetated sites were mainly driven by the shallowwater generalists Pomatoschistus spp. and Gobius niger, and the stationary species G. flavescens, with G. flavescens found in the more structurally complex habitats while Pomatoschistus spp. and G. niger were only found in the unvegetated areas. Interestingly, G. niger has been reported to have very strong habitat selection shifts based on the presence of predators, avoiding seagrass meadows completely when certain predator species are present and only inhabiting unvegetated soft bottom areas (Kruschel and Schultz, 2011), illustrating the complexity of trophic interactions. FIGURE 6 | Theoretical diagrams illustrating the habitat connectivity based on ontogenetic migrations of four species found within our study representing the four habitat preference guilds, based on literature as described in the text. (A) Clupea harengus, an occasional shallow-water visitor (OSV) habitat usage. Eggs found in deep water with larvae located in the pelagic zone (Reid et al., 1999;Petitgas et al., 2013). Juveniles and adults in large schools utilizing coastal waters and the open ocean (Reid et al., 1999). (B) Stationary species (SS) Gobiusculus flavescens; adults linked to seagrass and rock-macroalgae areas (Folkestad, 2005;Froese and Pauly, 2015), juveniles in seagrass habitats , eggs demersal and larvae in the sublittoral zone (Folkestad, 2005). (C) Juvenile migrant (JM) Gadus morhua; adults perform vertical migrations, juveniles utilize multiple shallow-water habitats (Staveley et al., unpublished data;Wennhage and Pihl, 2002), and larvae and eggs are in the pelagic zone (Petitgas et al., 2013). (D) Shallow-water generalist (SWG) Ctenolabrus rupestris. Juveniles and adults are associated with seagrass and rock-macroalgae Staveley et al., 2017), eggs are laid in the rock-macroalgae habitats (Darwall et al., 1992) and larvae in the pelagic zone (Kullander et al., 2012). Fish illustrations of G. flavescens and C. rupestris by Maj Persson. Other symbols courtesy of the Integration and Application Network (IAN).
While the abundance of different species differed among certain shallow-water habitats, the average size of fish within each habitat did not differ significantly. All habitats consisted of approximately 75% juveniles, indicating an important nursery function. These results are very strong evidence in support of the holistic seascape nursery concept initially posited by Nagelkerken et al. (2013). Although RUV is limited with regard to the number of fish that can be measured, its results are consistent across all seascapes and habitats. The number of measurable fish obtained is likely sufficient to draw relevant conclusions from, such as similarities in the life stage proportions among habitats.

CONCLUSIONS
Shallow-water habitats are integrally linked through nutrient flow and the movements of mobile organisms such as fish, and should therefore be studied and analyzed with this holistic perspective in mind (Ng et al., 2013). Here, we found support for the idea that coastal habitats, while differing in the abundance of fish contained within each habitat, all support juvenile fish assemblages. Our results, together with life history information of specific species, show that habitats within the shallow-water seascape are connected via species movements and ontogenetic shifts. Type of habitat was important for the abundance of fish but not for the diversity or number of species. Therefore, we suggest that the different habitats of shallow-water seascapes are of similar importance for fish assemblages and as a seascape nursery. Thus, management should mimic the ecological system and focus on a holistic-, heterogeneous-and interconnectedseascape approach, rather than take a single-habitat or singlespecies perspective.
AUTHOR CONTRIBUTIONS DP, TS, and MG all contributed to the design of the study, analysis and interpretation of the data. DP and TS were responsible for data acquisition. DP drafted the work and DP, TS, and MG participated in all manuscript revisions. DP, TS, and MG all gave final approval of the submitted manuscript and all agree to be accountable for every aspect of the work.

ACKNOWLEDGMENTS
We wish to thank Said Mgeleka and Mathew Silas for their help in assisting with all the field work as well as Linus Hammar for assistance during a portion of the field work. We are also grateful to Maria Asplund for her generosity in loaning equipment as well as the Sven Lovén Centre for Marine Sciences (Kristineberg) facility. Thank you Earl and Alex Perry for reviewing and editing the language of the manuscript. We gratefully acknowledge the support from the Swedish Research Council Formas (grant number 2011-1640). Additional support was graciously provided by a scholarship from the Royal Swedish Academy of Sciences.