Open Access
How to translate text using browser tools
1 January 2015 Quantification of Habitat and Community Relationships among Nearshore Temperate Fishes Through Analysis of Drop Camera Video
Ryan R. Easton
Author Affiliations +
Abstract

Temperate nearshore reefs along the Pacific coast of North America are highly valuable to commercial and recreational fisheries yet comprise a small fraction of the seabed. Monitoring fisheries resources in this region is difficult; high-relief structural complexity and adverse sea conditions have led to a paucity of information on temperate reef species assemblage patterns. Reliable, inexpensive tools and methods for monitoring are needed, as many traditional tools are both logistically complicated and expensive, limiting the frequency of their implementation over a large scale. Video drop cameras of varying designs have previously been employed to estimate fish abundance and distribution. We surveyed a nearshore rocky reef off the northern Oregon coast with a video lander (a video camera mounted on a landing platform so it can be dropped to the seafloor) over the spring and winter of 2011. We designed a 272-point systematic grid to document the species assemblage and the distribution and habitat associations of the reef species, including two overfished rockfishes: Canary Rockfish Sebastes pinniger and Yelloweye Rockfish Sebastes ruberrimus. Species assemblages differed significantly across the reef by depth and by season for the outer part of the reef. Well-defined habitat associations existed for many species; Canary Rockfish were associated with complex moderate-relief habitat types such as large boulders and small boulders, while Yelloweye Rockfish were associated with high-relief habitats like vertical walls. Species associations were evaluated pairwise to identify nearshore complexes. We compared our site with five exploratory reef sites off the central Oregon coast and found that nearshore reefs differed from our site, while offshore reefs were more similar. Video landers provide a solution to the need for increased sampling of temperate reef systems that are subject to difficult conditions and can contribute to habitat mapping, fish abundance indices, and fish assemblage information for monitoring and management of fisheries resources.

Temperate nearshore reefs along the Pacific coast of North America represent areas of valuable economic resources for commercial and recreational fisherman (Williams and Ralston 2002; Fox et al. 2004; Gunderson et al. 2008). In Oregon's territorial sea (out 5.6 km), nearshore rocky reefs make up only a small fraction of the total area (∼7%), with the remaining region comprised predominately of sand and unconsolidated sediments (D. Fox, Oregon Department of Fish and Wildlife, personal communication). However, these reefs constitute much of the Essential Fish Habitat designated for many pelagic and demersal fishes, which currently inhabit the nearshore region (PFMC 2005; PFMC 2011).

Monitoring fisheries resources and habitat in temperate reefs is difficult due to their high-relief structural complexity, adverse sea conditions that are common in temperate regions, and depths that often exceed those safe for visual surveys by scuba (Adams et al. 1995; Williams et al. 2010). Therefore, there is a paucity of information on species assemblage patterns on temperate reefs at different times of the year and along depth gradients. More information is needed to determine how the reefs function as critical habitat and to refine the fine-scale habitat associations of the fish that utilize those habitats (Gunderson et al. 2008). Annual trawl surveys conducted by the National Marine Fisheries Service currently cover the continental shelf from Cape Flattery, Washington, to the USA—Mexico border, but this survey does not come close to shore and trawls are unable to adequately catch fish that typically reside among rock escarpments and boulders (Zimmermann 2003; Cordue 2007). There have been repeated calls for more comprehensive sampling, particularly for reef-associated species that are considered to be below or near overfishing thresholds (Yoklavich et al. 2007; Williams et al. 2010).

Visual survey tools, such as remotely operated vehicles (ROVs) and human-occupied vehicles, are useful to survey untrawlable rocky habitats. These methods collect valuable information regarding the distribution, relative abundance, and species—habitat associations of various fish species, further aiding Essential Fish Habitat designation (Stein et al. 1992; Krieger 1993; Adams et al. 1995; Johnson et al. 2003; Yoklavich et al. 2007). However, the expense and expertise required for these methods, as well as confounding depth and sea conditions, can make them prohibitive to employ over a large scale within shallow (<70 m), highly productive nearshore waters. Recently, in the nearshore waters of the U.S. West Coast, there have been efforts to establish Marine Protected Areas and Marine Reserves, as well as to expand comprehensive multibeam seafloor mapping. Reliable, inexpensive tools and methods to monitor the effects of protected areas on local fish stocks are needed, as traditional tools (ROVs, autonomous underwater vehicles, and human-occupied vehicles) are both logistically complicated and expensive, generally limiting the frequency of their implementation. Video lander systems (a video camera mounted on a landing platform so it can be dropped to the seafloor) may be the solution to the increased sampling needed to ground-truth habitat maps and determine the relative abundance and distribution of nearshore fisheries resources over broad areas and in winter months.

In northeastern Pacific Ocean waters, highly variable and changing conditions represent obstacles to those researchers who are looking to obtain nonextractive visual data on nearshore Pacific rocky-reef fishes. Video drop cameras of varying designs have previously been employed as noninvasive tools for estimating fish abundance and distribution. Baited underwater video stations and baited remote underwater video stations have been shown to be effective in estimating the relative abundance of many fish that are mobile or solitary and that have low population sizes or avoid other visual survey methods (Ellis and DeMartini 1995; Priede and Merrett 1996; Willis and Babcock 2000; Cappo et al. 2004; Harvey et al. 2007; Stobart et al. 2007; Hannah and Blume 2012; Wakefield et al. 2013). Nonbaited underwater photo and video lander platforms have been effective in capturing accurate and repeatable fish and habitat data throughout a wide range of depths and habitat types without artificially attracting fish with bait, which would bias potential species-habitat associations (Gledhill et al. 1996; Roberts et al. 2005; Hannah and Blume 2012).

The objectives of this research were multifaceted, with a primary aim to determine the ability of a low-cost drop camera system to comprehensively survey a temperate nearshore rocky reef. We documented the species assemblage and the distribution and habitat associations of nearshore Pacific rocky-reef fishes in spring and winter to describe the distribution of key fished species, including two overfished rockfishes that are under intensive “stock rebuilding plans” by the Pacific Marine Fisheries Council. The video lander was evaluated as a survey tool for monitoring protected areas in nearshore temperate reef complexes, while concurrently assessing the ability to be used as a comprehensive ground-truthing tool to identify habitat types that are currently used by multibeam sonar surveys. Finally, we analyzed fish community composition with a recently developed quantitative method that measures pairwise species co-occurrence (Stone and Roberts 1990; Ulrich and Gotelli 2010; Groundfish Management Team 2013).

METHODS

The video lander we used is an autonomous underwater video system designed and built for use in high-relief rocky habitat by the Oregon Department of Fish and Wildlife's Marine Resources Program (Hannah and Blume 2012). The video lander is composed of an aluminum tube frame, a Deep Sea Power and Light (DSPL) Multi-SeaCam 2060, dual DSPL LED Ritelites (850 lm, 3,000 K), and an aluminum pressure housing containing two 13.2V rechargeable NiMH battery packs, a controller board, a Sony TRV-11 digital camcorder (recording video received from the DSPL Multi-SeaCam 2060) recording onto 60-min Mini DVC cassette tapes, and either a depth activated pressure switch used in waters deeper than 18 m or a push activated switch for shallower depths (Figure 1). The sacrificial base is designed so that, if stuck in rocky habitat, the lander can release from the base and rotate around multiple attachment points to maximize the retrieval probability in high-relief rocky habitat. The digital video footage from each lander drop was transferred from the original Sony DVC 60-min cassettes into digital format on a personal desktop computer with Adobe Premiere Pro through a firewire-connected Sony GV-HD700 portable video recorder deck. The video lander was deployed unbaited to avoid drawing in fish from other habitat types near the sampling point.

We selected the Three Arch Rocks rocky-reef complex located off Oceanside, Oregon, approximately 11 km south of the entrance to Tillamook Bay, for its broad depth range and known species diversity (Figure 2). The structure of the reef is a horseshoe pattern, running approximately 5 km east—west and 2 km north—south (Figure 3). Three Arch Rocks reef has a broad depth range, from surface to approximately 75 m as it runs east to west, and is known to support a high diversity of marine species (E. Schindler, Oregon Department of Fish and Wildlife, personal communication).

The video lander was initially deployed on the Three Arch Rocks reef on 12 separate days between April 17 and June 28, 2011. A systematic grid consisting of 272 individual drop points, spaced 175 m apart, was designed to maximize coverage of the reef structure while capturing all possible habitat types throughout the reef's entire depth range (Figure 3). The grid spacing was designed to maintain independence, while reducing the possibility for double counting individual fish during any given sampling day (Matthews 1990a, 1990b, 1992; Pacunski and Palsson 2002). The entire grid was blocked into seven regions prior to the survey, with the goal of completing at least one blocked section each day. The video lander platform used for the spring survey was again employed for winter sampling in December 2011, with the addition of a set of 10-cm paired scaling lasers to better quantify substrate grain size. Two separate attempts were made to complete the grid between December 1 and December 9, 2011; however, due to poor underwater visibility conditions we were only able to successfully survey a contiguous block of 70 drops comprising the outer quarter of the grid.

FIGURE 1.

The video lander platform utilized at Three Arch Rocks reef. Displayed in the photo are the Deep Sea Power and Light (DSPL) Multi-Sea-Cam 2060 (1), dual DSPL LED Mini-Sealites (2), pressure tube containing batteries and Sony TRV-11 digital camcorder (3), sacrificial (breakaway) base (4), and steel-rod weight bar (5), with arrows showing the break-away connection points.

f01_87.jpg

FIGURE 2.

Three Arch Rocks rocky-reef study area off of Oceanside, Oregon, located approximately 11 km south of the entrance to Tillamook Bay. The study area covers approximately 15 km2.

f02_87.jpg

The video lander was deployed for a fixed duration during daylight hours at each sampling location, following the protocols in Hannah and Blume (2012). Each video sample consisted of 5 min of recorded bottom time, beginning at the estimated time of the lander reaching the seafloor and ending when retrieval began. Five minutes of recorded bottom time allowed for sufficient sediment settling, as well as capturing a maximum count of species present. Completed video tapes were reviewed aboard the vessel to determine underwater visibility and if any drops needed to be repeated based on low water clarity, camera orientation, or visual obstruction.

FIGURE 3.

The video lander systematic survey grid of the Three Arch Rocks reef completed in the spring (April—June) of 2011. Each point represents an individual drop site within the grid by calendar day (175-m spacing).

f03_87.jpg

During the spring survey, 415 individual drops were completed over 12 boat-days, providing over 48.5 h of video footage. Of the 12 total sampling days, 7 d showed good bottom visibility, 1 d showed moderate bottom visibility, and 4 d showed little to no underwater visibility on the bottom, where turbidity obscured the view to the point that neither habitat nor fish were discernable, thus requiring resampling of these sites. A total of 143 drops had to be repeated due to poor underwater visibility or an obstructed view. The final sample size for analysis was 272 drops, a single drop for each of the 272 sampling locations. We used the drop with the highest score for visibility and view in cases when an individual site required resampling because habitat type or species were indiscernible due to high turbidity or marine snow. Acceptable weather for the winter survey occurred between December 6 and 9, 2011, and yielded a total of 108 usable drops across the Three Arch Rocks reef grid. Of these 108 usable drops, a continuous block of 70 drops comprising the outer reef section were used for species composition comparison with the spring survey results.

Following the field deployments, videos were reviewed in the laboratory by the primary author. Video review consisted of two separate components. The initial review was used to describe camera visibility and view, as well as topographic relief (Table 1). Topographic relief was defined as flat, low, or high depending on the observed habitat type at each drop location. We did not have a way to accurately measure distance sampled by the camera but roughly classified visibility conditions on a scale from 0 to 2 (poor, medium, good) following the criteria in Table 1 (Hannah and Blume 2012). Primary habitat (dominant habitat type in the camera's view) and secondary habitat (second most abundant habitat feature in view) were classified for each drop based on the habitat criteria shown in Table 2 (Hannah and Blume 2012). Drops that had view and visibility scores of 0 were excluded from analysis. The second review was used to identify fish observed on the video, which were identified to the lowest taxonomic level possible, usually to species. The maximum count of individuals from each species in any single frame (MaxN) was recorded for each drop to eliminate the potential for double counting of individuals (Harvey et al. 2007). Fish observed as the video lander was being retrieved were not included in the maximum count. Habitat observations made from the review of video lander footage were compared with a habitat classification map of the Three Arch Rocks region developed and provided by the Active Tectonics and Seafloor Mapping Lab (ATSML) at Oregon State University. The ATSML habitat maps are developed from multibeam sonar scans that collect data on a 4-m × 4-m (16 m2) grid pixel size. To develop a habitat map, these grid boxes are smoothed into 10-m × 10-m (100 m2) mapping units for which the dominant habitat type is displayed in a habitat box.

TABLE 1.

Criteria used to classify relief, underwater visibility, and view when reviewing video lander footage from Three Arch Rocks reef. Video footage had to receive at least a 1 in the Visibility or View categories to be used in further analysis.

t01_87.gif

We investigated the species—habitat associations of the 9 most abundant fish species observed on the Three Arch Rocks reef during the spring survey and the 10 most abundant fish species observed during the winter survey. These included the following: Black Rockfish Sebastes melanops, Blue Rockfish Sebastes mystinus, Canary Rockfish Sebastes pinniger, Quillback Rockfish Sebastes maliger, Yelloweye Rockfish Sebastes ruberrimus, Yellowtail Rockfish Sebastes flavidus, Lingcod Ophiodon elongatus, Kelp Greenling Hexagrammos decagrammus, Spotted Ratfish Hydrolagus colliei (winter only), and Pile Perch Damalichthys vacca. The number of individual fish species was totaled for each drop location to determine species richness. A Fisher's exact test was used to compare presence-absence data for each of the most abundant fish species on the reef to both of the primary and secondary habitat types identified. Unidentified adult fish were excluded from all analyses; however, unidentified juvenile rockfish Sebastes spp. were included in the analyses and treated as their own category.

TABLE 2.

Habitat criteria used to classify the primary and secondary habitat types observed at the Three Arch Rocks reef complex from video lander survey footage (Hannah and Blume 2012).

t02_87.gif

Differences in species assemblages and relative abundance (maximum count per drop in a single frame; MaxN) between sections of the reef were analyzed with a pairwise one-way analysis of similarity (ANOSIM), adjusted for multiple comparisons with a step-down sequential Bonferroni correction, using the software package Paleontological Statistics 2.15. The degree of difference in the species composition was then measured using the Bray—Curtis dissimilarity index (BCDI) in the SIMPER routine (Bray and Curtis 1957; Hammer et al. 2001; Hannah and Blume 2012). Every drop, including those with no observations, was included in both analyses. The relative abundance of each species observed among sections of the study site was compared using a nonparametric Wilcoxon test, as has been previously used with video lander data (Hannah and Blume 2012). Additionally, the species assemblage at Three Arch Rocks reef was compared with the species assemblages of five exploratory sites surveyed by a video lander and presented in Hannah and Blume (2012). This assemblage comparison was again performed using pairwise one-way ANOSIMs, adjusted for multiple comparisons using a step-down sequential Bonferroni correction in Paleontological Statistics 2.15.

FIGURE 4.

Total distribution of habitat types (see Table 2) across the Three Arch Rocks reef as observed by the video lander (April through June 2011) overlaid on a habitat classification map (≥ 100-m2 mapping unit patch size) developed by Oregon State University's Active Tectonics and Seafloor Mapping Lab (ATSML) for the state waters mapping project. Video-lander-observed habitat classifications are shown as a divided circle with primary habitat type on the left and secondary habitat type on the right; if only one habitat type was observed the circle is shown as a contiguous color. (For simplification, the video lander habitat classification legend shows only the reference color of the primary habitat type observed).

f04_87.jpg

Pairwise species co-occurrence was evaluated using the checkerboard score (C-score) metric to provide a single score of co-occurrence of a pair of species using presence-absence data (Stone and Roberts 1990; Ulrich and Gotelli 2010; Groundfish Management Team 2013):

e01_87.gif

where Ki = the number of occurrences of species i, Kj = the number of occurrences of species j, and Sij = the number of co-occurrences of species i and j. The C-score analysis provides a normalized value of 0–1, where 1 indicates perfect segregation between the two species and 0 indicates complete overlap. The Groundfish Management Team considers C-scores above 0.70 as a strong indication that the two species are segregated and scores of 0.30 and below as the two species exhibiting a high degree of overlap.

RESULTS

Habitat characterization over the entire survey grid revealed that sand was the most abundant habitat type for both primary (41.2%) and secondary (35.3%) habitat. Bedrock outcrop was the second most frequent primary habitat type (21.7%), with the high-relief habitat types (large boulder [4.0%], vertical wall [5.9%], and crevice [3.7%]) and moderate-to-low-relief habitat types (cobble [9.2%], gravel-pebble [5.9%], and small boulder [8.4%]) registering lower in overall frequency (Figure 4). We compared the habitat classifications characterized from the review of video lander footage with habitat maps developed by Oregon State University's ATSML for the state waters mapping project. We found an 80.1% agreement of primary habitat type between our habitat classifications and the ATSML classifications at the lowest resolution of the two classifications. Differences in agreement between the two methodologies is primarily due to scale; the video lander provides a view over a relatively small area, while habitat maps generated from multibeam sonar and backscatter data generally blend multiple similar habitat types into more uniform classifications. This habitat smoothing inherently leads to a loss of fine-scale habitat resolution in multibeam sonar habitat mapping, which is information the video lander is able to provide. Discrepancies along the edges between habitat types may indicate important areas for resampling.

TABLE 3.

Habitat associations (based on P-values obtained from Fisher's exact test) of the 10 most abundant species observed by the video lander during both the spring (April—June 2011) and winter (December 2011) surveys of the Three Arch Rocks reef complex. Kelp Greenling is the only fish with clear sexual dimorphism that could be observed in the videos. Significant P-values are color-coded by survey period: yellow = spring significant positive association, blue = winter significant positive association, green = spring and winter significant positive association, and red = significant negative association. If a significant P-value was identified during both survey periods, the “less significant” of the two seasons is presented. The abbreviations for habitat type are as follows: bedrock outcrop (BR), large boulder (LB), small boulder (SB), crevice (CR), vertical wall (VW), cobble (CO), gravel—pebble (GP), and sand (SA).

t03a_87.gif

Continued.

t03b_87.gif

Species Abundance and Distribution

While many (46% in spring, 36% in winter) of the drops performed over the course of the survey did not have any fish observed, the majority of the sites surveyed by the video lander had one or more fish species present. Over the course of the spring survey, the nine most abundant species observed by the video lander showed distinct habitat and depth associations, while three other identified fishes (Copper Rockfish Sebastes caurinus, Cabezon Scorpaenichthys marmoratus, and unidentified juvenile rockfishes) were not observed with enough regularity to identify any significant associations (Table 3).

In total, 745 individual rockfish were observed, among 939 total individual fish of all species, as well as 17 Dungeness crab Metacarcinus magister (Table 4). Over the course of the spring survey, Black Rockfish was the most abundant species overall (34 stations, 277 individuals), with Canary Rockfish being the second most abundant (41 stations, 225 individuals) (Table 4). The most abundant and frequently observed demersal rockfish species was Yelloweye Rockfish (22 stations, 27 individuals) (Table 4). Kelp Greenlings were the most frequently observed fish other than rockfish on the reef (67 stations, 90 individuals), while Lingcod were second (48 stations, 61 individuals) (Table 4). These two hexagrammid fishes were also the species with the broadest habitat associations and depth distributions.

During the winter survey, pelagic schooling rockfish (Black Rockfish, Blue Rockfish, and Yellowtail Rockfish) were the most abundant group of fishes observed on the reef, with Yellowtail Rockfish (14 stations, 92 individuals) being the single most abundant and frequently observed of the three (Table 4). Canary Rockfish (19 stations, 40 individuals) was the second most abundant species overall, excluding Northern Anchovy Engraulis mordax, and were observed in the greatest frequency (Table 4). The most abundant and frequently observed demersal rockfish was Yelloweye Rockfish (10 stations, 15 individuals), followed by Quillback Rockfish (4 stations, 5 individuals) (Table 4).

In winter, Kelp Greenlings (16 stations, 17 individuals) exhibited a broad distribution across depth and habitat, while Lingcod appeared noticeably absent (5 stations, 6 individuals) (Table 4). This absence of Lingcod may be due to spawning migration, as Lingcod are known to nest in shallow-water habitats in winter months (Matthews 1992; O'Connell 1993; Martell et al. 2000). Spotted Ratfish, which were not observed during the spring survey, were observed at 11 stations during December over multiple habitat types. Additionally, schools of varying size of Northern Anchovy, also not observed during the spring survey, were observed at 16 drop stations across the reef in December (Table 4). In total, 213 identified rockfish were observed, among 273 total fish (excluding Northern Anchovy) and one Dungeness crab (Table 4). Some of the species showed consistency in sighting locations between spring and winter surveys. For instance, of the 10 sites where Yelloweye Rockfish were observed in December, 5 were the same as the spring survey, while Kelp Greenlings were observed again at 8 of 16 total sites. It is impossible to know if these were the same individual fish, but consistency does provide strong evidence for habitat association.

TABLE 4.

Maximum count (n), number of drop sites each species was observed at (Drops), relative abundance (%; n/total number of fish observed), and rank abundance of observed fish taxa on the Three Arch Rocks reef over both the spring (April—June 2011) and winter (December 2011) surveys; NA denotes lack of observation.

t04_87.gif

FIGURE 5.

Distribution of (A) Canary Rockfish and (B) Yelloweye Rockfish across the Three Arch Rocks reef as observed by the video lander over both the spring (April—June 2011) and winter (December 2011) surveys combined. Slight differences in position location between spring and winter are due to the current shifting the final recorded drop location between seasons.

f05_87.jpg

FIGURE 6.

Spring (April—June 2011) species composition of Three Arch Rocks reef divided into three depth categories; inner (89 drops), middle (95 drops), and outside (88 drops) across the reef based on video lander video analysis (arrows represent the reef sections compared for species composition and abundance with corresponding Bray—Curtis dissimilarity index [BCDI] values and P-values from pairwise one-way ANOSIMs); RF = rockfish.

f06_87.jpg

TABLE 5.

Normalized C-scores (pairwise associations) of the 10 most abundant species observed by the video lander during the spring survey. A C-score of 0 indicates the total overlap of two species, while a C-score of 1 would indicate that the two species are never found together. In general, C-scores less than 0.3 (dark gray) are thought to signify species pairs that commonly associate, scores between 0.3 and 0.7 (light gray) indicate moderate association, while scores over 0.7 indicate very low association (unshaded).

t05_87.gif

Habitat associations were often significant but varied for some species between winter and spring sampling (Table 3). Nearly all of the species we observed showed negative correlations with sand habitat. Schooling pelagic rockfish (Black Rockfish, Blue Rockfish, and Yellowtail Rockfish) all showed a qualitative relationship with the reef structure, exhibiting a distribution pattern that mirrored the shape of the reef. Yelloweye Rockfish and Canary Rockfish, both currently managed as overfished stocks on the West Coast, exhibited significant associations with a variety of habitat types across the Three Arch Rocks reef (Table 3). During each survey period, Canary Rockfish had the broadest distribution across depths and habitat types of any of the observed rockfish species (Figure 5A). Yelloweye Rockfish on the other hand, while exhibiting significant relationships with high-vertical-relief habitat types (Table 4), showed a more restricted distribution than Canary Rockfish, with the vast majority of observations occurring on the outer third portion of the reef (Figure 5B).

Reef Fish Community Composition

We conducted an analysis of the species composition across the reef to investigate how species composition changed with depth and location. The survey grid was divided into three segments moving east to west across the reef: inside (89 drops), middle (95 drops), and outside (88 drops) (Figure 6). Comparison of the habitat composition across the three reef sections (inside, middle, outside) showed that only the inner and outer sections differed significantly (BCDI = 75.88; ANOSIM: P < 0.0427). This result was primarily driven by differences in sand, bedrock, and small boulder (secondary habitat) between these two regions of the reef. The species assemblages of the three sections all differed significantly in composition (Figure 6). The inside section, dominated primarily by Black Rockfish and Kelp Greenling, had the lowest overall species abundance and richness. This section also differed significantly in species composition from the middle (BCDI = 68.82; ANOSIM: P < 0.0029) and the outside (BCDI = 71.38; ANOSIM: P < 0.0001) sections of the reef. The middle section was dominated by Black Rockfish, followed by Blue Rockfish and Canary Rockfish. The middle section also differed significantly from the outside section (BCDI = 79.11; ANOSIM: P < 0.0334), which was dominated by Canary Rockfish, but also had the highest abundance of Yellowtail Rockfish, Yelloweye Rockfish, and Quillback Rockfish (Figure 6).

A comparison of the species composition between the spring and winter surveys, restricted to the contiguous outer-most 70 stations of the Three Arch Rocks reef survey that were sampled in both seasons, revealed significant differences in species composition (BCDI = 76.41; ANOSIM: P < 0.0155). This difference was driven by the winter presence of Spotted Ratfish at the reef and the overall lower winter abundance of Canary Rockfish (Wilcoxon test: P < 0.0168), Kelp Greenling (P < 0.0241), and Lingcod (P < 0.0456).

We investigated species correlations using the technique utilized by the Northwest Fisheries Science Center to identify species that commonly associate with others (C-score; Groundfish Management Team 2013) (Table 5). Black Rockfish and Blue Rockfish were commonly associated, with C-scores of 0.15 and 0.07 in the spring and winter, respectively. During the winter survey, strong pairwise associations were observed for Yelloweye Rockfish and Quillback Rockfish, Blue Rockfish and Yellowtail Rockfish, and Yelloweye Rockfish and Canary Rockfish. Other species pairs exhibited intermediate scores, likely due to low sample sizes (Table 5).

Nearshore Species Assemblage Variation on the Oregon Coast

The presence or absence and relative abundance of nearshore temperate fishes is likely to be quite variable from reef to reef due to differences in environmental conditions and fishing pressure, as well as year-to-year or seasonal movements and recruitment events (Gunderson et al. 2008). Three Arch Rocks reef is unique in its structure in that it runs east to west, spanning a wide range of depths, giving it a very diverse assemblage of temperate Pacific reef fishes. We compared our results from the spring survey with those from five exploratory sites, four nearshore (Cape Perpetua, East Siletz, West Siletz, and Seal Rocks) and one offshore (Stonewall Bank), off the central Oregon coast (presented in Hannah and Blume 2012) and found some interesting differences (Table A.1 in the appendix). The species assemblage at the Three Arch Rocks reef was significantly different from that of Cape Perpetua, approximately 150 km south (BCDI = 86.02; ANOSIM: P < 0.0011). The Cape Perpetua reef, which runs north to south, showed greater relative abundances of Canary Rockfish and juvenile rockfish but lower relative abundances of Kelp Greenling and Black Rockfish during the July survey (Table A.1 in the appendix). The differences were less conclusive for the February survey at Cape Perpetua (BCDI = 81.07; ANOSIM: P < 0.064), but this was likely due to the smaller numbers of fish observed in winter months. The Eastern Siletz reef showed the greatest disparity in species assemblage in comparison with the Three Arch Rocks reef (BCDI = 88.88; ANOSIM: P < 0.0001) due to the high relative abundance of Blue Rockfish and Kelp Greenling at the Eastern Siletz site (Table A.1 in the appendix). Similar results were observed when comparing with the Seal Rocks site (BCDI = 85.75; ANOSIM: P < 0.0002) due to the greater relative abundance of Black Rockfish, Canary Rockfish, Kelp Greenling, and juvenile rockfish observed at Seal Rocks (Table A.1 in the appendix). In contrast, the Western Siletz reef site was more similar to Three Arch Rocks (BCDI = 78.62; ANOSIM: P < 0.0721) (Table A.1 in the appendix). Similar results were observed when comparing the species assemblage to Stonewall Bank (BCDI = 70.63; ANOSIM: P < 0.0574), due primarily to the absence of Black Rockfish at Stonewall Bank and the relative greater abundance of Canary Rockfish and juvenile rockfish there but also the greater relative abundance of Lingcod at Three Arch Rocks (Table A.1 in the appendix).

DISCUSSION

Assessing and monitoring the distribution and abundance of fishes inhabiting rocky-reef structures along the West Coast of North America is a significant challenge. Typically, scientists must take a dynamic approach, utilizing a myriad of sources of data and a variety of data collection techniques to identify the relative abundance of species and their associated habitat types across geographic regions (Francis 1986; Parker et al. 2000). In order for resource surveys to be effective, a sizeable amount of planning, funding, and personnel from multiple agencies is often needed. However, even with substantial time and effort, the results of these surveys can still be highly variable, with poor representation of some habitat types in areas that are difficult to sample (Krieger 1993; Jagielo et al. 2003). Our study provides a comprehensive look at how video drop camera survey data can be analyzed to contribute to our understanding of nearshore fishes and their habitat associations.

The limitations of the most commonly used survey gear types (i.e., bottom trawl) are well documented and understood both by those who employ the gear and those who utilize the data in stock assessments (Adams et al. 1995; Williams and Ralston 2002). The question, however, is how to integrate different sources of data and promote the use of novel sampling methods that can overcome the shortcomings of currently used sampling methods for nearshore rocky-reef species. Previous work has demonstrated that surveys utilizing various direct video observation techniques, such as video landers, baited underwater video stations, and scuba, are effective at surveying nearshore, shallow-water fish populations (Watson et al. 2005; Langlois et al. 2010; Hannah and Blume 2012). Additionally, video data analyzed from direct observations has been shown to provide relative abundance estimates of a variety of fish species through nonextractive, fishery-independent means (Gledhill et al. 2006; Yoklavich et al. 2007; Coleman et al. 2011; Merritt et al. 2011). Much of this work, however, has been performed in areas such as Hawaii, southern California, and the Gulf of Mexico— regions which are not inhibited to the same extent by the survey challenges present in the highly productive and often turbulent waters of the northeastern Pacific Ocean. Our research demonstrates that a video lander platform can comprehensively survey a temperate nearshore reef under difficult conditions, collecting broad-scale fish assemblage, relative abundance, and habitat data. This illustrates the opportunity for video landers to be developed into a more widely utilized, cost-effective survey tool to monitor nearshore rocky reefs along the Pacific coast of North America and other poorly surveyed systems.

Visual surveys, however, are not without their own limitations. Visibility is the single most important factor when it comes to a successful video survey, be it with the video lander, scuba divers, an ROV, or a submersible. The video lander has an advantage over other visual methods in its simplicity and ease of use and deployment, allowing minimal time and effort for data collection relative to conventional ROVs or submersibles (Stein et al. 1992; Krieger and Ito 1999; Johnson et al. 2003; Yoklavich et al. 2007; Pacunski et al. 2008). The relative simplicity, small size, and ease of use allows for rapid deployment of the video lander when weather conditions become favorable, as well as minimal cost to abort a survey when conditions prove unsuitable for sampling.

The data from this study provide information on habitat associations for overfished and rarely surveyed species in both spring and winter seasons that can further contribute to stock assessment and spatial management. The species—habitat associations of many demersal species determined from the spring survey were reinforced by the winter survey (e.g., Yelloweye Rockfish), while others were expanded upon (e.g., Quillback Rockfish) (Table 3). However, due to the low sample size of these species, more drops will be required to refine and substantiate these results. Additionally, the ability of the video lander to accurately identify species-habitat associations for pelagic schooling rockfish (Black Rockfish, Blue Rockfish, and Yellowtail Rockfish) may be limited and requires further investigation. Furthermore, these species have reduced interactions with the benthic substrate compared with demersal species (i.e., Yelloweye Rockfish, Quillback Rockfish), with much of their time spent schooling in the water column, potentially decreasing the capacity of the video lander to survey these species.

As with previous studies, the video lander exhibited limitations in its ability to identify and count flatfishes and very cryptic species like Cabezon due to the oblique angle of the camera and the standard definition resolution of the camera system (Hannah and Blume 2012). Additionally, we were unable to identify juvenile rockfish and other small fish to species due to the low resolution of standard-definition video, as well as some fish being too distant from the camera, an issue that will be improved with advances in high-definition and stereo video camera systems (Hannah and Blume 2014). The extent to which the video lander attracts or repels fish, as well as the limited and highly variable size of the area viewed, represent potential sampling biases which are currently not quantified. Despite this, extensive research has shown maximum counts from video surveys to be an accurate, although conservative, index of relative abundance (Willis and Babcock 2000; Stoner et al. 2008; Merritt et al. 2011). Other research has suggested using a mean count of fish observed in single snapshots over the entire course of the video in place of a single maximum count to improve the estimate of true abundance (Schobernd et al. 2014). Continued analysis of these methods of estimating relative abundance will be needed to further evaluate the video lander as an effective survey tool.

The C-score analysis of the video lander datasets offers another dimension with which to examine the distribution of nearshore species across the reef. Previous C-score analysis focused on deepwater slope rockfish and “other” roundfish found in trawl samples enumerated by the West Coast Groundfish Observer Program and the Alaska Fisheries Science Center (Groundfish Management Team 2013). These are robust datasets, both spatially and temporally, but are lacking in nearshore species like those observed by the video lander at Three Arch Rocks. Because C-scores are calculated from presence— absence data, the video lander provides ideal high-resolution data amenable to C-score analysis, enabling it to contribute to an additional assessment need. The video lander provides data on both overlap and segregation of species pairs, allowing for potential targeted management actions at the species, or species pair, level. Further surveys are required to expand on the results as many of the observed species had low total counts; however, preliminary analysis shows the versatility and robustness of video lander data for this type of analysis. Our analysis of C-scores for key groundfish species compliments the analyses performed by the Groundfish Management Team by providing results for nearshore, high-relief habitats that are not sampled in existing surveys.

Our comparison of observations at Three Arch Rocks reef to those of other experimental survey sites off the central Oregon coast displays the uniqueness and complexity of the species assemblage at different sites (Table A.1 in the appendix). While these locations are within 150 km of Three Arch Rocks reef, the east-to-west distribution of Three Arch Rocks reef creates a larger depth range and distance from shore, with a reef structure that is home to both shallow (Blue Rockfish and Black Rockfish) and deeper-water (Yelloweye Rockfish and Canary Rockfish) species. This combination of depths and species diversity at Three Arch Rocks reef highlights the ecological diversity of this area and the need for additional surveys to characterize nearshore species assemblages and habitats in Oregon.

The ability to survey nearshore reefs off the Pacific coast is sporadic, and the weather windows are generally short. The video lander therefore is an ideal tool to use in all seasons because of its rapid, intensive survey capability combined with a short preparation and implementation schedule. The winter survey yielded interesting results, including the first video observations of Spotted Ratfish at Three Arch Rocks reef. While visibility issues plagued the winter survey, it provided new insight into this nearshore rocky-reef environment in winter. Given acceptable ocean conditions, the video lander has the capability to perform at the same level in the winter as it did in the spring.

Finally, the video lander may serve as a critical tool for evaluations of Marine Protected Areas and Marine Reserves, which have recently been employed as a conservation tool in Oregon and elsewhere in the United States. Nonextractive monitoring methods are needed, as many of these areas are closed to fishing, even for scientific study. The high spatial coverage, low operational and logistical cost, and nonextractive nature of video lander surveys make them an ideal tool for Marine Protected Area or Marine Reserve assessment, monitoring, and habitat ground-truthing. Further refinement of video quality, processing, and distance-area sampling will make video landers essential components of nearshore reef monitoring.

ACKNOWLEDGMENTS

Steve and Ray Dana provided the commercial passenger fishing vessel Blue Water Too as a sampling platform for the spring survey. Al Pazar provided the research vessel Pacific Surveyor as a sampling platform for the winter survey. Chris Goldfinger of the ATSML at Oregon State University provided funding for the winter survey, while Chris Romsos of ATSML provided habitat maps and technical support. Matthew T.O. Blume from the Marine Resources Program in the Oregon Department of Fish and Wildlife assisted with field sampling during the spring survey and provided technical support throughout. This manuscript was improved through the input of W. Waldo Wakefield at the National Oceanic and Atmospheric Administration's Northwest Fisheries Science Center in Newport, Oregon, Kevin Thompson and Scarlett Arbuckle at Oregon State University, and two anonymous reviewers.

REFERENCES

1.

P. B. Adams , J. L. Butler , C. H. Baxter , T. E. Laidig , K. A. Dahlin , and W. Wakefield . 1995. Population estimates of Pacific coast groundfishes from video transects and swept-area trawls. U.S. National Marine Fisheries Service Fishery Bulletin 93:446–455. Google Scholar

2.

J. R. Bray , and J. T. Curtis . 1957. An ordination of the upland forest communities of southern Wisconsin. Ecological Monographs 27:325–349. Google Scholar

3.

M. Cappo , P. Speare , and G. De'ath . 2004. Comparison of baited remote underwater video stations (BRUVS) and prawn (shrimp) trawls for assessments of fish biodiversity in inter-reefal areas of the Great Barrier Reef Marine Park. Journal of Experimental Marine Biology and Ecology 302:123–152. Google Scholar

4.

F. C. Coleman , K. M. Scanlon , and C. C. Koenig . 2011. Groupers on the edge: shelf edge spawning habitat in and around marine reserves of the northeastern Gulf of Mexico. Professional Geographer 63:456–474. Google Scholar

5.

P. L. Cordue 2007. A note on non-random error structure in trawl survey abundance indices. ICES Journal of Marine Science 64:1333–1337. Google Scholar

6.

D. M. Ellis , and E. E. DeMartini . 1995. Evaluation of a video camera technique for indexing abundances of juvenile Pink Snapper Pristipomoides filamentosus, and other Hawaiian insular shelf fishes. U.S. National Marine Fisheries Service Fishery Bulletin 93:67–77. Google Scholar

7.

D. Fox , A. Merems , M. Amend , H. Weeks , C. Romsos , and M. Appy . 2004. Comparative characterization of two nearshore rocky reef areas: a high-use recreational fishing reef vs. an unfished reef. Oregon Department of Fish and Wildlife, Newport. Google Scholar

8.

R. C. Francis 1986. Two fisheries biology problems in West Coast groundfish management. North American Journal of Fisheries Management 6:453–462. Google Scholar

9.

C. T. Gledhill , G. W. Ingram Jr , K. R. Rademacher , P. Felts , B. Trigg , and M. S. Pascagoula . 2006. NOAA Fisheries reef fish video surveys: yearly indices of abundance for Gag (Mycteroperca microlepis), SEDAR10-DW-12. Southeast Data Assessment and Review, North Charleston, South Carolina. Google Scholar

10.

C. T. Gledhill , J. Lyczkowski-Shultz , K. Rademacher , E. Kargard , G. Crist , and M. A. Grace . 1996. Evaluation of video and acoustic index methods for assessing reef-fish populations. ICES Journal of Marine Science 53:483–485. Google Scholar

11.

Groundfish Management Team. 2013. Groundfish management team report on methods and results that may be used to evaluate alternatives for stock complex reorganization. Pacific Fishery Management Council, June 2013 Briefing Book, Agenda Item F.8.b, Portland. Google Scholar

12.

D. R. Gunderson , A. M. Parma , R. Hilborn , J. M. Cope , D. L. Fluharty , M. L. Miller , R. D. Vetter , S. S. Heppell , and H. G. Greene . 2008. The challenge of managing nearshore rocky reef resources. Fisheries 33:172–179. Google Scholar

13.

O. Hammer , D. A. T. Harper , and P. D. Ryan . 2001. PAST: Paleontological statistics package for education and data analysis. Paleontologia Electrobica [online serial] 4:article 4. Google Scholar

14.

R. W. Hannah , and M. T. O. Blume . 2012. Tests of an experimental unbaited video lander as a marine fish survey tool for high-relief deepwater rocky reefs. Journal of Experimental Marine Biology and Ecology 430–431:1–9. Google Scholar

15.

R. W. Hannah , and M. T. O. Blume . 2014. The influence of bait and stereo video on the performance of a video lander as a survey tool for marine demersal reef fishes in Oregon waters. Marine and Coastal Fisheries: Dynamics, Management, and Ecosystem Science [online serial] 6:181–189. Google Scholar

16.

E. S. Harvey , M. Cappo , J. J. Butler , N. Hall , and G. A. Kendrick . 2007. Bait attraction affects the performance of remote underwater video stations in assessment of demersal fish community structure. Marine Ecology Progress Series 350:245–254. Google Scholar

17.

T. Jagielo , A. Hoffmann , J. Tagart , and M. Zimmermann . 2003. Demersal groundfish densities in trawlable and untrawlable habitats off Washington: implications for the estimation of habitat bias in trawl surveys. U.S. National Marine Fisheries Service Fishery Bulletin 101:545–565. Google Scholar

18.

S. W. Johnson , M. L. Murphy , and D. J. Csepp . 2003. Distribution, habitat, and behavior of rockfishes, Sebastes spp., in nearshore waters of southeastern Alaska: observations from a remotely operated vehicle. Environmental Biology of Fishes 66:259–270. Google Scholar

19.

K. J. Krieger 1993. Distribution and abundance of rockfish determined from a submersible and by bottom trawling. U.S. National Marine Fisheries Service Fishery Bulletin 91:87–96. Google Scholar

20.

K. J. Krieger , and D. H. Ito . 1999. Distribution and abundance of Shortraker Rockfish, Sebastes borealis, and Rougheye Rockfish, S. aleutianus, determined from a manned submersible. U.S. National Marine Fisheries Service Fishery Bulletin 97:264–272. Google Scholar

21.

T. J. Langlois , E. S. Harvey , B. Fitzpatrick , J. J. Meeuwig , G. Shedrawi , and D. L. Watson . 2010. Cost-efficient sampling of fish assemblages: comparison of baited video stations and diver video transects. Aquatic Biology 9:155–168. Google Scholar

22.

S. J. D. Martell , C. J. Walters , and S. S. Wallace . 2000. The use of marine protected areas for conservation of Lingcod (Ophiodon elongatus). Bulletin of Marine Science 66:729–743. Google Scholar

23.

K. R. Matthews 1990a. An experimental study of the habitat preferences and movement patterns of Copper, Quillback, and Brown rockfishes (Sebastes spp.). Environmental Biology of Fishes 29:161–178. Google Scholar

24.

K. R. Matthews 1990b. A telemetric study of the home ranges and homing routes of Copper and Quillback rockfishes on shallow rocky reefs. Canadian Journal of Zoology 68:2243–2250. Google Scholar

25.

K. R. Matthews 1992. A telemetric study of the home ranges and homing routes of Lingcod Ophiodon elongatus on shallow rocky reefs off Vancouver Island, British Columbia. U.S. National Marine Fisheries Service Fishery Bulletin 90:784–790. Google Scholar

26.

D. Merritt , M. K. Donovan , C. Kelley , L. Waterhouse , M. Parke , K. Wong , and J. C. Drazen . 2011. BotCam: a baited camera system for nonextractive monitoring of bottomfish species. U.S. National Marine Fisheries Service Fishery Bulletin 109:56–67. Google Scholar

27.

V. M. O'Connell 1993. Submersible observations of Lingcod, Ophiodon elongatus, nesting below 30 m off Sitka, Alaska. Marine Fisheries Review 55:19–24. Google Scholar

28.

R. E. Pacunski , and W. A. Palsson . 2002. Macro- and micro-habitat relationships of adult and sub-adult rockfish, Lingcod, and Kelp Greenling in Puget Sound. Washington Department of Fish and Wildlife, Olympia. Google Scholar

29.

R. E. Pacunski , W. A. Palsson , H. G. Greene , and D. Gunderson . 2008. Conducting visual surveys with a small ROV in shallow water. Pages 109–128 in J. R. Reynolds and H. G. Greene , editors. Marine habitat mapping technology for Alaska. Alaska Sea Grant for North Pacific Research Board, Fairbanks. Google Scholar

30.

S. J. Parker , S. A. Berkeley , J. T. Golden , D. R. Gunderson , J. Heifetz , M. A. Hixon , R. Larson , B. M. Leaman , M. S. Love , J. A. Musick , V. M. O'Connell , S. Ralston , H. J. Weeks , and M. M. Yoklavich . 2000. Management of Pacific rockfish. Fisheries 25(3):22–30. Google Scholar

31.

PFMC (Pacific Fishery Management Council). 2005. Pacific coast groundfish management plan for the California, Oregon, and Washington groundfish fishery. Appendix B, part 1. Assessment methodology for groundfish essential fish habitat. PFMC, Portland, Oregon. Google Scholar

32.

PFMC (Pacific Fishery Management Council). 2011. Pacific coast groundfish fishery management plan for the California, Oregon, and Washington groundfish fishery. PFMC, Portland, Oregon. Google Scholar

33.

I. G. Priede and N. R. Merrett . 1996. Estimation of abundance of abyssal demersal fishes; a comparison of data from trawls and baited cameras. Journal of Fish Biology 49:207–216. Google Scholar

34.

J. M. Roberts , O. C. Peppe , L. A. Dodds , D. J. Mercer , W. T. Thomson , J. D. Gage , D. T. Meldrum , A. Freiwald , and J. M. Roberts . 2005. Monitoring environmental variability around cold-water coral reefs: the use of a benthic photolander and the potential of seafloor observatories. Pages 483–502 in A. Freiwald , editor. Cold-water corals and ecosystems. Springer Berlin-Heidelberg, Berlin. Google Scholar

35.

Z. H. Schobernd , N. M. Bacheier , and P. B. Conn . 2014. Examining the utility of alternative video monitoring metrics for indexing reef fish abundance. Canadian Journal of Fisheries and Aquatic Sciences 71:464–471. Google Scholar

36.

D. L. Stein , B. N. Tissot , M. A. Hixon , and W. H. Barss . 1992. Fish-habitat associations on a deep reef at the edge of the Oregon continental shelf. U.S. National Marine Fisheries Service Fishery Bulletin 90:540–551. Google Scholar

37.

B. Stobart , J. A. García-Charton , C. Espejo , E. Rochel , R. Goñi , O. Reñones , A. Herrero , R. Crec'hriou , S. Polti , C. Marcos , S. Planes , and A. Pérez-Ruzafa . 2007. A baited underwater video technique to assess shallow-water Mediterranean fish assemblages: methodological evaluation. Journal of Experimental Marine Biology and Ecology 345:158–174. Google Scholar

38.

L. Stone , and A. Roberts . 1990. The checkerboard score and species distributions. Oecologia 85:74–79. Google Scholar

39.

A. W. Stoner , B. J. Laurel , and T. P. Hurst . 2008. Using a baited camera to assess relative abundance of juvenile Pacific Cod: field and laboratory trials. Journal of Experimental Marine Biology and Ecology 354:202–211. Google Scholar

40.

W. Ulrich , and N. J. Gotelli . 2010. Null model analysis of species associations using abundance data. Ecology 91:3384–3397. Google Scholar

41.

C. B. Wakefield , P. D. Lewis , T. B. Coutts , D. V. Fairclough , and T. J. Langlois . 2013. Fish assemblages associated with natural and anthropogenically-modified habitats in a marine embayment: comparison of baited videos and opera-house traps. PLoS (Public Library of Science) ONE [online serial] 8:e59959. Google Scholar

42.

D. L. Watson , E. S. Harvey , M. J. Anderson , and G. A. Kendrick . 2005. A comparison of temperate reef fish assemblages recorded by three underwater stereo-video techniques. Marine Biology 148:415–425. Google Scholar

43.

E. H. Williams , and S. Ralston . 2002. Distribution and co-occurrence of rockfishes (family: Sebastidae) over trawlable shelf and slope habitats of California and southern Oregon. U.S. National Marine Fisheries Service Fishery Bulletin 100:836–855. Google Scholar

44.

K. Williams , C. N. Rooper , and R. Towler . 2010. Use of stereo camera systems for assessment of rockfish abundance in untrawlable areas and for recording Pollock behavior during midwater trawls. U.S. National Marine Fisheries Service Fishery Bulletin 108:352–362. Google Scholar

45.

T. J. Willis , and R. C. Babcock . 2000. A baited underwater video system for the determination of relative density of carnivorous reef fish. Marine and Freshwater Research 51:755–763. Google Scholar

46.

M. M. Yoklavich , M. S. Love , and K. A. Forney . 2007. A fishery-independent assessment of an overfished rockfish stock, Cowcod (Sebastes levis), using direct observations from an occupied submersible. Canadian Journal of Fisheries and Aquatic Sciences 64:1795–1804. Google Scholar

47.

M. Zimmermann 2003. Calculation of untrawlable areas within the boundaries of a bottom trawl survey. Canadian Journal of Fisheries and Aquatic Sciences 60:657–669. Google Scholar

Appendices

APPENDIX: Nearshore Species Assemblage Variation on the Oregon Coast

Table A.1.

Comparison of video lander survey data, including the total numbers of fish observed (summed maximum counts across stations of the 25 most abundant species between the five study locations), by species or group and survey area (n denotes the number of stations sampled). All data not from the Three Arch Rocks reef was obtained from Hannah and Blume 2012.

tA01_87.gif

Notes

[1] This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The moral rights of the named author(s) have been asserted.

© Ryan R. Easton, Selina S. Heppell, Robert W. Hannah
Ryan R. Easton "Quantification of Habitat and Community Relationships among Nearshore Temperate Fishes Through Analysis of Drop Camera Video," Marine and Coastal Fisheries: Dynamics, Management, and Ecosystem Science 7(7), 87-102, (1 January 2015). https://doi.org/10.1080/19425120.2015.1007184
Received: 18 October 2014; Accepted: 3 January 2015; Published: 1 January 2015
Back to Top