<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Archiving and Interchange DTD v2.3 20070202//EN" "archivearticle.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="methods-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mar. Sci.</journal-id>
<journal-title>Frontiers in Marine Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mar. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-7745</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmars.2023.1219283</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Marine Science</subject>
<subj-group>
<subject>Methods</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Evaluating the impacts of reduced sampling density in a systematic fisheries-independent survey design</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>DeFilippo</surname>
<given-names>Lukas</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2140233"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Kotwicki</surname>
<given-names>Stan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2333054"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Barnett</surname>
<given-names>Lewis</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1231170"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Richar</surname>
<given-names>Jon</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Litzow</surname>
<given-names>Michael A.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Stockhausen</surname>
<given-names>William T.</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2390993"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Palof</surname>
<given-names>Katie</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2385031"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Resource Assessment and Conservation Engineering Division, Alaska Fisheries Science Center, National Oceanic and Atmospheric Administration (NOAA)</institution>, <addr-line>Seattle, WA</addr-line>, <country>United States</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Resource Assessment and Conservation Engineering Division, Alaska Fisheries Science Center, National Oceanic and Atmospheric Administration (NOAA)</institution>, <addr-line>Kodiak, AK</addr-line>, <country>United States</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Resource Ecology and Fisheries Management Division, Alaska Fisheries Science Center, National Oceanic and Atmospheric Administration (NOAA)</institution>, <addr-line>Seattle, WA</addr-line>, <country>United States</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Alaska Department of Fish and Game, Division of Commercial Fisheries</institution>, <addr-line>Juneau, AK</addr-line>, <country>United States</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Nathan Bacheler, Beaufort Laboratory, Southeast Fisheries Science Center (NOAA), United States</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Sarah Roberts, University of North Carolina at Chapel Hill, United States; Matthew Vincent, Southeast Fisheries Science Center (NOAA), United States</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Lukas DeFilippo, <email xlink:href="mailto:lukas.defilippo@noaa.gov">lukas.defilippo@noaa.gov</email>
</p>
</fn>
<fn fn-type="present-address" id="fn003">
<p>&#x2020;Present address: Lukas DeFilippo, Ecosystem Monitoring and Assessment Program, Alaska Fisheries Science Center, NOAA, Seattle, WA, United States</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>25</day>
<month>08</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>10</volume>
<elocation-id>1219283</elocation-id>
<history>
<date date-type="received">
<day>08</day>
<month>05</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>07</day>
<month>08</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 DeFilippo, Kotwicki, Barnett, Richar, Litzow, Stockhausen and Palof</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>DeFilippo, Kotwicki, Barnett, Richar, Litzow, Stockhausen and Palof</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Fisheries-independent surveys provide critical data products used to estimate stock status and inform management decisions. While it can be possible to redistribute sampling effort to improve survey efficiency and address changing monitoring needs in the face of unforeseen challenges, it is important to assess the consequences of such changes. Here, we present an approach that relies on existing survey data and simulations to evaluate the impacts of strategic reductions in survey sampling effort. We apply this approach to assess the potential effects of reducing high density sampling near St. Matthew Island and the Pribilof Islands in the NOAA eastern Bering Sea (EBS) bottom trawl survey. These areas contain high density &#x201c;corner stations&#x201d; that were implemented for finer-scale monitoring of associated blue king crab stocks (<italic>Paralithodes platypus</italic>) which historically supported commercial fisheries but have since declined and are seldom eligible for harvest. We investigate the effects of removing these corner stations on survey data quality for focal <italic>P. platypus</italic> stocks and other crab and groundfish species monitored by the EBS survey. We find that removing the St. Matthew and Pribilof Islands corner stations has negligible effects on data quality for most stocks, except for those whose distributions are concentrated in these areas. However, the data quality for such stocks was relatively low even with higher density sampling, and corner station removal had only minor effects on stock assessment outcomes. The analysis we present here provides a generic approach for evaluating strategic reductions in sampling effort for systematic survey designs and can be applied by scientists and managers facing similar decisions elsewhere.</p>
</abstract>
<kwd-group>
<kwd>survey design</kwd>
<kwd>stock assessment</kwd>
<kwd>groundfish</kwd>
<kwd>crab</kwd>
<kwd>spatiotemporal model</kwd>
</kwd-group>
<counts>
<fig-count count="7"/>
<table-count count="4"/>
<equation-count count="9"/>
<ref-count count="70"/>
<page-count count="16"/>
<word-count count="8820"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Marine Fisheries, Aquaculture and Living Resources</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Fisheries-independent surveys provide some of the most important data used to inform stock assessments and management decisions (<xref ref-type="bibr" rid="B21">Gunderson, 1993</xref>; <xref ref-type="bibr" rid="B5">Chen et&#xa0;al., 2003</xref>). Indices of abundance generated from survey data are typically assumed to be proportional to the true abundance of target populations, a relationship which is expressed via the catchability parameter (<italic>q</italic>) in assessment models (<xref ref-type="bibr" rid="B48">Quinn and Deriso, 1999</xref>). The relationship between the survey index and population size is assumed to be constant in most cases, violation of which can lead to spurious estimates of stock dynamics (<xref ref-type="bibr" rid="B47">Pope and Shepherd, 1985</xref>; <xref ref-type="bibr" rid="B67">Wilberg et&#xa0;al., 2009</xref>; <xref ref-type="bibr" rid="B64">Thorson et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B30">Kotwicki et&#xa0;al., 2014</xref>). Consequently, changes to survey designs and operations that affect catchability may disrupt the stationarity of this relationship and confound true population trends if not accounted for (<xref ref-type="bibr" rid="B20">Godo, 1994</xref>; <xref ref-type="bibr" rid="B28">Kimura and Somerton, 2006</xref>; <xref ref-type="bibr" rid="B2">Cadigan et&#xa0;al., 2022</xref>). Changes in gear selectivity and effort allocation across space may also shift the composition of ages and sizes sampled, potentially affecting estimates of population age and size structure (<xref ref-type="bibr" rid="B40">Ono et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B31">Kotwicki et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B39">O&#x2019;Leary et&#xa0;al., 2020</xref>). As such, to the extent possible, fisheries monitoring agencies typically avoid altering survey designs and sampling procedures in order to minimize the impacts such changes could have on data products and population estimates.</p>
<p>While long-term time series of fisheries-independent survey data with minimal interruptions or modifications are desirable (<xref ref-type="bibr" rid="B20">Godo, 1994</xref>; <xref ref-type="bibr" rid="B55">Stompe et&#xa0;al., 2020</xref>), circumstances may arise where changes are unavoidable or beneficial. For instance, sections of historically surveyed areas may become inaccessible due to marine protected area designations or wind energy development (<xref ref-type="bibr" rid="B15">Field et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B35">Methratta et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B22">Hare et&#xa0;al., 2022</xref>). Similarly, hazardous weather, vessel breakdowns, staffing and budget shortages, as well as pandemics, international conflicts, and economic turmoil can disrupt monitoring agencies&#x2019; ability to complete surveys (<xref ref-type="bibr" rid="B24">ICES, 2020</xref>; <xref ref-type="bibr" rid="B50">Santora et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B25">ICES, 2023</xref>). Changes may also offer advantages, such as modernizing fishing gear or improving the efficiency of sampling designs (<xref ref-type="bibr" rid="B1">Brown et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B41">Oyafuso et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B42">Oyafuso et&#xa0;al., 2022</xref>). Moreover, marine ecosystems are dynamic, and updating survey designs may be necessary to adapt to shifting species distributions, abundance trends, and management priorities (<xref ref-type="bibr" rid="B13">Dulvy et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B46">Pinsky and Mantua, 2014</xref>; <xref ref-type="bibr" rid="B34">Maureaud et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B11">DeFilippo et&#xa0;al., 2023</xref>). While some fisheries-independent survey designs can readily accommodate changes in sampling density (<italic>e.g.</italic>, stratified random), others are less flexible (<italic>e.g.</italic>, systematic) (<xref ref-type="bibr" rid="B7">Cochran, 1977</xref>) and may require more detailed evaluation of the consequences of such changes.</p>
<p>The U.S. National Oceanic and Atmospheric Administration (NOAA) eastern Bering Sea (EBS) bottom trawl survey provides critical data in support of stock assessment and management for some of the world&#x2019;s most commercially valuable fish and crab stocks. The EBS survey follows a stratified systematic design composed of 350 rectangular grid cells which are each sampled as part of annual surveys (<xref ref-type="bibr" rid="B32">Lauth et&#xa0;al., 2019</xref>). However, near St. Matthew Island and the Pribilof Islands the corners of the grid cells are sampled in addition to the centers. The rationale for finer-scale sampling near St. Matthew and the Pribilof Islands has been to improve monitoring and data products for blue king crab (<italic>Paralithodes platypus</italic>) stocks in these areas, which historically supported valuable commercial fisheries. Blue king crab exhibit a sparse and patchy distribution, resulting in highly uncertain abundance estimates. Sampling these corner stations was instituted to increase the probability of encountering concentrated patches of blue king crab and reduce uncertainty in survey data products. However, both the Pribilof and St. Matthew stocks of blue king crab have declined substantially and are now closed to fishing with little sign of rebuilding to harvestable levels in the near future (<xref ref-type="bibr" rid="B43">Palof et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B53">Stockhausen, 2021a</xref>). Consequently, it is unclear if the effort and funds required to continue sampling the corner stations is justifiable, or if resources might be better allocated to other priorities, such as improving data quality for other species.</p>
<p>Here, we present an approach for evaluating the consequences of reduced sampling effort in systematic survey designs and apply it to estimate the effects of removing the St. Matthew and Pribilof Island corner stations from the EBS survey. While changes to survey fishing gear and sampling protocols can be resolved by paired fishing and intercalibration studies (<italic>e.g.</italic>, <xref ref-type="bibr" rid="B3">Cadigan and Dowden, 2010</xref>; <xref ref-type="bibr" rid="B36">Miller, 2013</xref>; <xref ref-type="bibr" rid="B31">Kotwicki et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B2">Cadigan et&#xa0;al., 2022</xref>), there is less guidance for evaluating changes in sampling effort (but see <xref ref-type="bibr" rid="B70">Zimmermann and Enberg, 2017</xref>; <xref ref-type="bibr" rid="B24">ICES, 2020</xref>; <xref ref-type="bibr" rid="B25">ICES, 2023</xref>), particularly for systematic survey designs. Using existing survey time-series and simulation analyses, our approach estimates the impacts of effort reduction on the precision and accuracy of survey data products, as well as stock assessment output and biological and management reference points. As fisheries management agencies are tasked with monitoring changing marine ecosystems under static or declining budgets, tools for addressing effort reduction and reallocation decisions will be essential for optimizing survey efficiency and ensuring reliable data products (<xref ref-type="bibr" rid="B24">ICES, 2020</xref>; <xref ref-type="bibr" rid="B25">ICES, 2023</xref>). The approach we present here is generic and flexible and can be widely applied to other species and regions by scientists and managers facing survey effort reduction decisions.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Case study background</title>
<p>The NOAA EBS bottom trawl survey occurs southeast of the U.S. &#x2013; Russian international maritime boundary from Bristol Bay and the Alaska Peninsula to the south, to north of Nunivak and St. Matthew Island (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>). Occurring annually from May to early August, the EBS survey samples a fixed set of 350 rectangular grid cells as part of a systematic design with a minimum grid resolution of 37.04 km<sup>2</sup> and a maximum depth of 200m (<xref ref-type="bibr" rid="B32">Lauth et&#xa0;al., 2019</xref>). The EBS survey begins in Bristol Bay and proceeds west using two chartered commercial vessels. Each vessel tows a standard 83&#x2019;-112&#x2019; eastern otter trawl with 10cm mesh for a duration of 30 minutes at a target speed of 3 knots (<xref ref-type="bibr" rid="B32">Lauth et&#xa0;al., 2019</xref>). Survey effort is measured as the area swept by the trawl gear, which is calculated as the product of the distance fished (measured with a GPS and a bottom contact sensor) and net width (measured by an acoustic sensor). Catches of commercially important fishes are identified to species and weighed, with a subset selected for length measurements and otolith extraction following protocols outlined in <xref ref-type="bibr" rid="B32">Lauth et&#xa0;al. (2019)</xref> and <xref ref-type="bibr" rid="B51">Stauffer (2004)</xref>. Catches of commercially important crab species are sexed, assigned a shell condition and maturity status as per classifications in <xref ref-type="bibr" rid="B27">Jadamec et&#xa0;al. (1999)</xref>, and carapace width (<italic>Chionoecetes</italic> spp.) or carapace length (<italic>Lithodes</italic>/<italic>Paralithodes</italic> spp.) is measured to the nearest 0.1&#xa0;mm using either Vernier (prior to 2016) or digital (2016 and later) calipers.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Bathymetric map of the southeastern Bering Sea shelf and centroids of the NOAA EBS bottom trawl survey grid. The St. Matthew&#x2019;s and Pribilof Islands corner stations are shown in purple, while the standard grid center stations are shown in blue.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1219283-g001.tif"/>
</fig>
<p>Surrounding St. Matthew and the Pribilof Islands, the corners of the EBS survey grid cells are sampled in addition to the centers, leading to higher density sampling in these areas (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>). There are 26 of these &#x201c;corner stations&#x201d; in total, which, at an average rate of 4 stations/day, require roughly one week of vessel time to sample. This additional sampling effort was initiated in 1981 for the Pribilof Islands, and 1983 for St. Matthew Island in order to improve the quality of data products for blue king crab stocks in these areas. Blue king crab exhibit a &#x201c;contagious&#x201d; distribution, in which individuals tend to be encountered in sparse, concentrated patches. Moreover, blue king crab (particularly females) are often found on rocky substrate and thus are difficult to sample by the EBS bottom trawl survey which primarily occurs in soft bottom habitat (<xref ref-type="bibr" rid="B69">Zacher et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B65">V&#xf8;lstad, 2012</xref>). As a result, bottom trawl survey data products for blue king crab are highly uncertain (<xref ref-type="bibr" rid="B53">Stockhausen, 2021a</xref>; <xref ref-type="bibr" rid="B44">Palof, 2022</xref>). Sampling at a finer spatial resolution via the addition of the corner stations was intended to reduce uncertainty by increasing both the probability of encountering patches of blue king crab, and overall sample sizes. However, following pronounced declines in abundance, the St. Matthew blue king crab stock was declared overfished and closed to harvest in 1999. After a 10-year rebuilding plan, fishing resumed in 2009/2010 but after three years of modest harvest, the fishery was closed again in 2013/2014 due to declining survey abundance estimates and concerns about the productivity of the stock. Fishing resumed in 2014/2015 but fishery catches were poor and harvest has remained closed since 2016; the stock was once again declared overfished in 2018 when it entered into the current rebuilding plan (<xref ref-type="bibr" rid="B44">Palof, 2022</xref>). The Pribilof blue king crab stock has been closed to fishing since 1999, was declared overfished in 2002 and subsequently put on a rebuilding plan (<xref ref-type="bibr" rid="B53">Stockhausen, 2021a</xref>).</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Analysis overview</title>
<p>Our approach to evaluating the effects of survey effort reduction proceeds in three general stages. In the first stage (Empirical Analysis, section 2.2.1), we quantify the retrospective effects of station removal on core survey data products (<italic>i.e.</italic> biomass, compositional estimates) that are used as inputs in stock assessment models. This stage of the analysis is focused on (1) identifying qualitative differences in survey data products as a result of station removal (<italic>i.e.</italic> changes in estimated stock biomass trajectories), and (2) evaluating effects on the uncertainty/precision of the survey data products (<italic>i.e.</italic> changes in the estimated variance of survey biomass indices). This first stage of the analysis uses historical survey data, generating survey data products without any corner station data and comparing them to estimates produced with the data from these stations. The second stage of our analysis focuses on assessing the effects of station removal on the accuracy of the survey data products using simulations (section 2.2.2. Simulation Analysis). Specifically, new survey datasets are simulated from spatiotemporal operating models and we examine how removing the corner stations from these simulated datasets affects the accuracy of estimated biomass indices relative to the true values specified in the operating models. In the third stage of our approach, we propagate the effects of station removal into the stock assessment models (section 2.2.3. Stock Assessment Analysis). This is accomplished in our case study by fitting stock assessment models to the biomass and compositional survey data products generated without corner station data and comparing the output to that from model fits that included the corner station data. Because integrated stock assessment models are complex and time consuming to fit, we do not perform this exercise for the full suite of stocks under consideration in our analysis, but select a representative subset to evaluate. Given that the corner stations were specifically intended to improve data quality for crab, our case study focuses more heavily on crab stocks, with a similar but limited set of analyses performed for groundfish.</p>
<sec id="s2_2_1">
<label>2.2.1</label>
<title>Empirical analysis</title>
<p>For the first stage of our analyses, we retrospectively evaluated the effects of removing corner station data from the existing survey time-series on derived data products. To do so, we withheld all data collected from the 26 corner stations from 1982-2019 and produced design-based (<italic>i.e.</italic> area-swept catch-per-unit effort expansions) and model-based (<italic>i.e.</italic> standardized indices produced via spatiotemporal models) estimates using only data from the remaining stations. While the corner stations were implemented to improve data quality specifically for blue king crab, it is important to understand the effects of their removal on estimates for other species that occur in these areas as well. As such, in addition to both the St. Matthew and Pribilof Islands blue king crab stocks, we evaluated the effects of corner station removal on Pribilof Islands red king crab (<italic>P. camtschaticus</italic>), snow crab (<italic>Chionoecetes opilio</italic>) and Tanner crab (<italic>C. bairdi</italic>) biomass indices. For all crab biomass estimates, we focused on key size classes for each stock that are particularly important from a management standpoint. Specifically, we examined biomass estimates for the GE103 (GE=greater than or equal to) size class (103 &#x2013; 189&#xa0;mm) of male Tanner crab, the GE95 size class (95mm &#x2013; 178&#xa0;mm) of male snow crab, the GE120 size class of male Pribilof blue and red king crab (120 &#x2013; 173&#xa0;mm and 120 &#x2013; 209&#xa0;mm respectively), and both the GE90 (90 &#x2013; 173&#xa0;mm) and GE105 (105&#x2013;173 mm) size classes of St. Matthew blue king crab. Additionally, we evaluated the effects of corner station removal on data products for ten groundfish species that inhabit these areas: Pacific cod (<italic>Gadus macrocephalus</italic>), walleye pollock (<italic>G. chalcogrammus</italic>), yellow Irish lord (<italic>Hemilepidotus jordani</italic>), Alaska plaice (<italic>Pleuronectes quadrituberculatus</italic>), northern rock sole (<italic>Lepidopsetta polyxystra</italic>), yellowfin sole (<italic>Limanda aspera</italic>), Bering flounder (<italic>Hippoglossoides robustus</italic>), flathead sole (<italic>H. elassodon</italic>), Alaska skate (<italic>Bathyraja parmifera</italic>), and Bering skate (<italic>B. interrupta</italic>). While improving data quality for groundfish stocks was not a consideration in implementing the corner stations as part of the EBS survey, groundfish are nonetheless sampled at these locations and it is useful to consider impacts on these stocks&#x2019; data as well.</p>
<p>Design-based biomass indices for all species in our study were produced following the methods outlined in <xref ref-type="bibr" rid="B66">Wakabayashi et&#xa0;al. (1985)</xref> and using standard protocols for EBS groundfish and crab index production (<xref ref-type="bibr" rid="B32">Lauth et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B69">Zacher et&#xa0;al., 2020</xref>). For the crab stocks in our analysis, we also produced design-based size compositional estimates with and without the corner station data; these results are presented in Appendix A.</p>
<p>Model-based estimators are increasingly being used for index standardization, as they can often improve the precision of indices from fisheries-independent survey data (<xref ref-type="bibr" rid="B63">Thorson et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B4">Cao et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B62">Thorson and Haltuch, 2019</xref>). Model-based estimates are currently included in stock assessments for some of the groundfish species considered in this study, but none of the crab stocks. However, by leveraging spatial autocorrelation to extrapolate to unsampled areas, some model-based estimators may be more robust to reductions in survey effort. As such, we compared model-based biomass indices produced with and without the corner station data to explore how these estimates responded to station removal. Model-based estimates were generated using spatiotemporal models built via the Vector Autoregressive Spatiotemporal (VAST) package (release number 3.9.0) (<xref ref-type="bibr" rid="B60">Thorson and Barnett, 2017</xref>) in R-4.1.2 (<xref ref-type="bibr" rid="B49">R Core Team, 2021</xref>). The specification of the spatiotemporal models used here is standard for abundance index production, and further details can be found in <xref ref-type="bibr" rid="B61">Thorson (2019)</xref>. Briefly, VAST models are an extension of a generalized linear mixed model (GLMM) that estimates dependent variable(s) via two linear predictors and a link function. Variation in the response variable(s) over space and time is partitioned into three components: (1) temporal variation (<inline-formula>
<mml:math display="inline" id="im1">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>, which represents changes from year-to-year that are equal across all locations, (2) spatial patterns (<inline-formula>
<mml:math display="inline" id="im2">
<mml:mi>&#x3c9;</mml:mi>
</mml:math>
</inline-formula>), which correspond to variation over space that is constant over time (<italic>i.e</italic>. long-term habitat associations), and (3) spatiotemporal variation (<inline-formula>
<mml:math display="inline" id="im3">
<mml:mi>&#x3f5;</mml:mi>
</mml:math>
</inline-formula>), which represents changes from year-to-year that are expressed differently across locations. To account for zero-inflated and skewed distributions, a Poisson-link delta modeling approach is used with two estimated linear predictors, <inline-formula>
<mml:math display="inline" id="im4">
<mml:mi>n</mml:mi>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im5">
<mml:mi>w</mml:mi>
</mml:math>
</inline-formula>, which represent expected numerical density and biomass-per-individual, respectively, such that <inline-formula>
<mml:math display="inline" id="im6">
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> gives the expected biomass density (<inline-formula>
<mml:math display="inline" id="im7">
<mml:mrow>
<mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) of sample (survey haul) <inline-formula>
<mml:math display="inline" id="im8">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> (<xref ref-type="bibr" rid="B61">Thorson, 2019</xref>):</p>
<disp-formula>
<label>(1)</label>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:mi>log</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>

<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mi>log</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>

<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>

</mml:mrow>
</mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im9">
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im10">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the location and year associated with sample <inline-formula>
<mml:math display="inline" id="im11">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula>. The annual intercepts (<inline-formula>
<mml:math display="inline" id="im12">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) were specified as fixed effects independent among years, and the spatial variation terms (<inline-formula>
<mml:math display="inline" id="im13">
<mml:mi>&#x3c9;</mml:mi>
</mml:math>
</inline-formula>) were estimated as random effects following a multivariate normal distribution:</p>
<disp-formula>
<label>(2)</label>
<mml:math display="block" id="M3">
<mml:mrow>
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x223c;</mml:mo>
<mml:mtext>MVN</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mstyle mathvariant="bold" mathsize="normal">
<mml:mi>R</mml:mi>
</mml:mstyle>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3b7;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im14">
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the marginal spatial variance and <inline-formula>
<mml:math display="inline" id="im15">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3b7;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the correlation matrix among locations (<inline-formula>
<mml:math display="inline" id="im16">
<mml:mi>s</mml:mi>
</mml:math>
</inline-formula>) which is modeled as a Mat&#xe9;rn function with decorrelation distance of <inline-formula>
<mml:math display="inline" id="im17">
<mml:mi>&#x3b7;</mml:mi>
</mml:math>
</inline-formula> and a transformation matrix that allows for geometric anisotropy such that decorrelation distance varies with cardinal direction (<xref ref-type="bibr" rid="B63">Thorson et&#xa0;al., 2015</xref>). Spatiotemporal effects were specified similarly:</p>
<disp-formula>
<label>(3)</label>
<mml:math display="block" id="M4">
<mml:mrow>
<mml:mstyle mathvariant="bold" mathsize="normal">
<mml:mi>&#x3b5;</mml:mi>
</mml:mstyle>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x223c;</mml:mo>
<mml:mtext>MVN</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c3;</mml:mi>
<mml:mtext>&#x3b5;</mml:mtext>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mstyle mathvariant="bold" mathsize="normal">
<mml:mi>R</mml:mi>
</mml:mstyle>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>&#x3b7;</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im18">
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the marginal spatiotemporal variance. The predicted density of individuals for each sample <inline-formula>
<mml:math display="inline" id="im19">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> follows a Poisson process with expectation <inline-formula>
<mml:math display="inline" id="im20">
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> such that the encounter probability (<inline-formula>
<mml:math display="inline" id="im21">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) is defined as:</p>
<disp-formula>
<label>(4)</label>
<mml:math display="block" id="M5">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>exp</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mi>n</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im22">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> is the area swept for bottom trawl sample <inline-formula>
<mml:math display="inline" id="im23">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula>. For any years in which a given species exhibited a 100% encounter rate across the EBS survey (<italic>i.e.</italic> at least one individual encountered at every station) the encounter probability (<inline-formula>
<mml:math display="inline" id="im24">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) was fixed at one. The positive catch rate <inline-formula>
<mml:math display="inline" id="im25">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> for sample <inline-formula>
<mml:math display="inline" id="im26">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> is obtained from the numerical density (<inline-formula>
<mml:math display="inline" id="im27">
<mml:mi>n</mml:mi>
</mml:math>
</inline-formula>) and average biomass per individual (<inline-formula>
<mml:math display="inline" id="im28">
<mml:mi>w</mml:mi>
</mml:math>
</inline-formula>) as:</p>
<disp-formula>
<label>(5)</label>
<mml:math display="block" id="M6">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Given the predicted encounter probability <inline-formula>
<mml:math display="inline" id="im29">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and positive catch rate <inline-formula>
<mml:math display="inline" id="im30">
<mml:mrow>
<mml:mtext>&#xa0;</mml:mtext>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, the probability distribution of the biomass <inline-formula>
<mml:math display="inline" id="im31">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> for sample <inline-formula>
<mml:math display="inline" id="im32">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> was specified as:</p>
<disp-formula>
<label>(6)</label>
<mml:math display="block" id="M7">
<mml:mrow>
<mml:mtext>Pr</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi>B</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>B</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2022;</mml:mo>
<mml:mi>G</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mo>|</mml:mo>
<mml:msubsup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>B</mml:mi>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im33">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the probability associated with a biomass of zero, and <inline-formula>
<mml:math display="inline" id="im34">
<mml:mrow>
<mml:mtext>Gamma</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mo>|</mml:mo>
<mml:msubsup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the probability of biomass <inline-formula>
<mml:math display="inline" id="im35">
<mml:mi>B</mml:mi>
</mml:math>
</inline-formula> given the expected encounter probability <inline-formula>
<mml:math display="inline" id="im36">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and positive catch rate <inline-formula>
<mml:math display="inline" id="im37">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> with estimated Gamma shape and scale parameters <inline-formula>
<mml:math display="inline" id="im38">
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im39">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>The model form described here was fitted individually to each of the ten species of groundfish considered and each of the six crab stocks/size classes, with the exceptions of Pribilof blue and red king crab. The full model configuration described above failed to converge for these latter two stocks, likely due to a limited spatial distribution and low encounter rates. As such, we removed the second spatiotemporal term <inline-formula>
<mml:math display="inline" id="im40">
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> from the Pribilof blue king crab model and the second temporal (<inline-formula>
<mml:math display="inline" id="im41">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and spatial (<inline-formula>
<mml:math display="inline" id="im42">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>), and both spatiotemporal terms (<inline-formula>
<mml:math display="inline" id="im43">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) from the Pribilof red king crab model. Additionally, to account for years in which no Pribilof red king crab were observed, the first temporal intercept (<inline-formula>
<mml:math display="inline" id="im44">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) was specified as an independent annual random effect rather than as a fixed effect.</p>
</sec>
<sec id="s2_2_2">
<label>2.2.2</label>
<title>Simulation analysis</title>
<p>The empirical analysis described above was focused on assessing the (1) qualitative impacts (<italic>i.e.</italic> changes in biomass trends) and (2) effects on precision (<italic>i.e.</italic> coefficients of variation (CVs) of biomass estimates) of corner station removal on existing survey data products. However, it is not possible from these analyses to evaluate the effects of corner station removal on the accuracy of survey data products. To do so, we conducted a simulation analysis using the VAST model fits for each species as the basis for spatiotemporal operating models (OMs). New fixed and random effects were simulated from the joint precision matrix of the spatiotemporal model fits (<italic>i.e.</italic> conditional on the original data), and new data were then simulated conditional upon these new fixed and random effects via parametric bootstrapping (<italic>e.g.</italic>, <xref ref-type="bibr" rid="B59">Thorson et&#xa0;al., 2021</xref>). New observations were simulated at each location and year in which sampling occurred in the original survey data set (including both positive observations, and observations of zero biomass). Data were simulated from model fits that included the corner station data and thus included simulated data points at the corner station locations. These simulated corner station data were then either withheld or retained as design and model-based biomass estimates were obtained using the simulated data. The accuracy of these estimates was measured using the log accuracy ratio (LAR):</p>
<disp-formula>
<label>(7)</label>
<mml:math display="block" id="M8">
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>A</mml:mi>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mtext>log</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>B</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im45">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>B</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the estimated (design or model-based) biomass index from year <inline-formula>
<mml:math display="inline" id="im46">
<mml:mi>t</mml:mi>
</mml:math>
</inline-formula> and simulation replicate <inline-formula>
<mml:math display="inline" id="im47">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula>, and <inline-formula>
<mml:math display="inline" id="im48">
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the true biomass value specified in the OM. LAR was summarized across years and simulation replicates via the median symmetric accuracy (MSA):</p>
<disp-formula>
<label>(8)</label>
<mml:math display="block" id="M9">
<mml:mrow>
<mml:mtext>MSA</mml:mtext>
<mml:mo>=</mml:mo>
<mml:mtext>exp</mml:mtext>
<mml:mo stretchy="false">(</mml:mo>
<mml:mtext>Median</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>|</mml:mo>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>A</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mo>|</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</disp-formula>
<p>LAR and MSA offer a number of advantages over other accuracy metrics (<italic>e.g.</italic>, mean absolute percentage error (MAPE), root mean square error (RMSE)), including scale independence, robustness to outliers, symmetry, and interpretability (<xref ref-type="bibr" rid="B37">Morley et&#xa0;al., 2018</xref>). LAR and MSA are interpretable such that values of zero represent perfect accuracy, and larger values indicate progressively worse accuracy.</p>
<p>The VAST operating models for Pribilof blue and red king crab and St. Matthew blue king crab were poorly conditioned, likely due to the low abundance levels, restricted ranges, and infrequent encounters for these stocks. Consequently, simulated datasets for Pribilof and St. Matthew king crab generated with the same procedures used for other stocks exhibited unrealistic distributions, with frequent positive encounters outside the Pribilof or St. Matthew management areas. As such, a modified simulation procedure was used for these stocks in which new data were simulated conditional upon the estimated fixed and random effects rather than simulating new fixed and random effects. Additionally, simulated datasets for St. Matthew and Pribilof king crab stocks were generated from operating models with spatial extents that were limited to the St. Matthew and Pribilof management areas respectively, thereby constraining the simulated observations within these boundaries. Given the difficulty of conditioning the operating models, we recommend that the results of the simulation analyses for St. Matthew&#x2019;s blue and Pribilof blue and red king crab be interpreted cautiously.</p>
</sec>
<sec id="s2_2_3">
<label>2.2.3</label>
<title>Stock assessment analysis</title>
<p>In addition to examining the effects of corner station removal on the precision and accuracy of survey data products themselves, we also investigated how these effects propagated into stock assessment output. For this phase of the analysis, we focus solely on crab stocks as these were the primary consideration for the implementation of the corner stations. We fitted stock assessment models using survey estimates of biomass (Methods: 2.2.2. Empirical analysis) and size composition (Appendix A) produced without the corner stations. We compared the resulting assessment model predictions and estimated biological and management reference points to those from model fits to data products that included corner station data. For this analysis we focused on Tanner crab and St. Matthew blue king crab. We selected these two stocks to investigate the effects of corner station removal on assessment output for stocks for which the corner stations are (St. Matthew blue king crab) and are not (Tanner crab) focal areas of the stocks&#x2019; distributions.</p>
<p>The Tanner crab assessment model is a stage/size-based population dynamics model that incorporates sex (male, female), shell condition (new shell, old shell), and maturity (immature, mature) as different categories into which the overall stock is divided on a size-specific basis (<xref ref-type="bibr" rid="B54">Stockhausen, 2021b</xref>). The model is fit using a penalized maximum likelihood approach to the design-based survey biomass and size composition time-series, molt increment data, retained catch biomass and size composition time series from the directed fishery, bycatch data (biomass and size composition time series) from the directed fishery (sub-legal males, all females), and bycatch data (biomass and size composition time series) from several other crab and groundfish fisheries. Management quantities are subsequently derived from the maximum likelihood solution using spawner-per-recruit proxies for F<sub>MSY</sub> and B<sub>MSY</sub>. For this analysis, the 2021 assessment model was fitted using the design-based survey biomass and size composition time series estimated without the corner stations and compared to existing assessment output produced with the corner stations. Model results for the estimated sex/maturity-specific population biomass time series, recruitment time series, management quantities, and rates of natural mortality from the two scenarios were compared.</p>
<p>The St. Matthew&#x2019;s blue king crab stock assessment is a simpler form of the Tanner crab model, using only mature male crab in the size/stage structured model (<xref ref-type="bibr" rid="B44">Palof, 2022</xref>). The stock is modeled via the Generalized Modeling for Crustacean Stocks (GMACS) framework in which it is fit using a penalized maximum likelihood approach to design-based survey biomass and size composition time series, retained catch biomass and size composition time series from the directed fishery, bycatch data (biomass and size composition time series) from the directed fishery, and bycatch data (biomass) from several other crab and groundfish fisheries. Estimates of life history parameters, such as natural mortality and growth for this stock are borrowed from other well studied king crab stocks, such as Bristol Bay red king crab. Management quantities are approximated from the long term average of mature male biomass, as directed by the Bering Sea and Aleutian Island (BSAI) crab fishery management plan. For this analysis, the 2021 model was fit using the design-based EBS bottom trawl survey biomass and size composition time series with and without corner stations. Resulting time series of mature male biomass and recruitment, as well as management quantities were compared between model fits.</p>
</sec>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Empirical analyses</title>
<p>For the six stocks/size classes of crab that we examined, corner station removal had little qualitative impact on either design or model-based biomass estimates. Biomass trends were generally coherent with one another regardless of corner station inclusion, with some notable transient discrepancies for St. Matthew blue king crab and Pribilof blue and red king crab (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). Corner station removal had little effect on the precision of design and model-based estimates for both snow crab and Tanner crab (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>; <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>), with increases in average CV (&#x394; CV) of&lt;0.02 for these species. More substantial declines in precision were observed for the other crab stocks, particularly for design-based biomass estimates (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>; <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). For the two size classes of St. Matthew blue king crab, average design-based CVs estimates increased by ~ 0.1, while model-based CVs were much less affected by removing the corner stations (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). The average design-based CV increased by 0.185 without the corner station data for Pribilof blue king crab, and 0.071 for Pribilof red king crab, while average the model-based CV for these stocks increased by 0.143 and 0.077 (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). While the precision of model-based estimates was generally more robust to corner station removal compared to that of design-based estimates, we note that this was not the case for Pribilof red king crab; the average model-based CV for this stock increased by more than its design-based counterpart as a result of corner station removal, although the model-based CVs themselves were smaller in all scenarios (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Design and model-based survey biomass estimates (left) and associated CVs (right) produced with and without corner station data for crab stocks/size classes.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1219283-g002.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Mean coefficients of variation (CV) of design and model-based biomass estimates for crab species produced with versus without the corner stations.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">&#xa0;<break/>Species/Size Class</th>
<th valign="top" align="center">Estimator</th>
<th valign="top" align="center">Mean CV (corners)</th>
<th valign="top" align="center">Mean CV&#xa0;(no corners)</th>
<th valign="top" align="center">&#x394; CV</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Tanner GE103</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.174</td>
<td valign="top" align="center">0.190</td>
<td valign="top" align="center">0.016</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.065</td>
<td valign="top" align="center">0.069</td>
<td valign="top" align="center">0.004</td>
</tr>
<tr>
<td valign="top" align="left">Snow GE95</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.147</td>
<td valign="top" align="center">0.158</td>
<td valign="top" align="center">0.011</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.064</td>
<td valign="top" align="center">0.065</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Pribilof blue king GE120</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.503</td>
<td valign="top" align="center">0.688</td>
<td valign="top" align="center">0.185</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.352</td>
<td valign="top" align="center">0.495</td>
<td valign="top" align="center">0.143</td>
</tr>
<tr>
<td valign="top" align="left">Pribilof red king GE120</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.648</td>
<td valign="top" align="center">0.719</td>
<td valign="top" align="center">0.071</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.318</td>
<td valign="top" align="center">0.395</td>
<td valign="top" align="center">0.077</td>
</tr>
<tr>
<td valign="top" align="left">St. Matthew blue king GE90</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.320</td>
<td valign="top" align="center">0.415</td>
<td valign="top" align="center">0.095</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.189</td>
<td valign="top" align="center">0.184</td>
<td valign="top" align="center">-0.005</td>
</tr>
<tr>
<td valign="top" align="left">St. Matthew blue King GE105</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.312</td>
<td valign="top" align="center">0.412</td>
<td valign="top" align="center">0.1</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.190</td>
<td valign="top" align="center">0.212</td>
<td valign="top" align="center">0.022</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>&#x394; CV represents the increase in mean CV between estimates produced without the corner stations relative to estimates produced with the corner stations.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>For nine of the ten groundfish species that we examined, there was no qualitative effect of corner station removal on the trend or scale of either design or model-based biomass estimates (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). Similarly, for most of the groundfish species we considered, differences in the precision of biomass estimates produced with and without corner station data were negligible (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>; <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). The precision of model-based groundfish biomass estimates was generally more robust to corner station removal than their design-based counterparts (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). The only groundfish species for which corner station removal caused a substantial increase in biomass CVs was yellow Irish lord (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>; <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). Average design-based CV for this species increased by 0.096, while model-based CVs were more robust and increased by 0.05 without the corner station data. However, it is worth noting that even with the corner station data included, the biomass CVs for yellow Irish lord were large in many years (<italic>e.g.</italic>, 0.3-0.6), and were the highest of any of the groundfish species we examined (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Design and model-based survey biomass estimates (left) and associated CVs (right) produced with and without corner station data for groundfish species. Note that Bering skate and Alaska skate were not distinguished on the NOAA EBS bottom trawl survey until 1996, hence why the time-series for those species do not begin until 1996.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1219283-g003.tif"/>
</fig>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Mean coefficients of variation (CV) of design and model-based biomass estimates for groundfish species produced with versus without the corner stations.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">&#xa0;Species</th>
<th valign="top" align="center">Estimator</th>
<th valign="top" align="center">Mean CV (corners)</th>
<th valign="top" align="center">Mean CV (no corners)</th>
<th valign="top" align="center">&#x394; CV</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Bering skate</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.165</td>
<td valign="top" align="center">0.167</td>
<td valign="top" align="center">0.002</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.118</td>
<td valign="top" align="center">0.119</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Alaska skate</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.061</td>
<td valign="top" align="center">0.062</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.042</td>
<td valign="top" align="center">0.043</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Flathead sole</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.110</td>
<td valign="top" align="center">0.118</td>
<td valign="top" align="center">0.008</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.061</td>
<td valign="top" align="center">0.060</td>
<td valign="top" align="center">-0.001</td>
</tr>
<tr>
<td valign="top" align="left">Bering flounder</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.169</td>
<td valign="top" align="center">0.173</td>
<td valign="top" align="center">0.004</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.123</td>
<td valign="top" align="center">0.125</td>
<td valign="top" align="center">0.002</td>
</tr>
<tr>
<td valign="top" align="left">Yellowfin sole</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.080</td>
<td valign="top" align="center">0.080</td>
<td valign="top" align="center">&lt;0.001</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.048</td>
<td valign="top" align="center">0.046</td>
<td valign="top" align="center">-0.002</td>
</tr>
<tr>
<td valign="top" align="left">Northern rock sole</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.088</td>
<td valign="top" align="center">0.094</td>
<td valign="top" align="center">0.006</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.060</td>
<td valign="top" align="center">0.059</td>
<td valign="top" align="center">-0.001</td>
</tr>
<tr>
<td valign="top" align="left">Alaska plaice</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.115</td>
<td valign="top" align="center">0.116</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.057</td>
<td valign="top" align="center">0.057</td>
<td valign="top" align="center">&lt;0.001</td>
</tr>
<tr>
<td valign="top" align="left">Yellow Irish lord</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.350</td>
<td valign="top" align="center">0.446</td>
<td valign="top" align="center">0.096</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.151</td>
<td valign="top" align="center">0.201</td>
<td valign="top" align="center">0.05</td>
</tr>
<tr>
<td valign="top" align="left">Pacific cod</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.085</td>
<td valign="top" align="center">0.089</td>
<td valign="top" align="center">0.004</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.047</td>
<td valign="top" align="center">0.048</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Walleye pollock</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.106</td>
<td valign="top" align="center">0.109</td>
<td valign="top" align="center">0.003</td>
</tr>
<tr>
<td valign="top" align="center">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.066</td>
<td valign="top" align="center">0.068</td>
<td valign="top" align="center">0.002</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>&#x394; CV represents the increase in mean CV between estimates produced without the corner stations relative to estimates produced with the corner stations.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Simulation analyses</title>
<p>For snow and Tanner crab, the accuracy of design and model-based estimates was relatively robust to removal of the corner station data (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>; <xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>). Conversely, the accuracy of biomass estimates for Pribilof blue and red king crab and St. Matthew blue king crab showed more substantial declines resulting from corner station removal (<xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>). The accuracy of model-based biomass estimates for the Pribilof and St. Matthew&#x2019;s king crab stocks was considerably more robust to corner station removal than were their design-based counterparts (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>; <xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>). For instance, removing simulated corner station data for Pribilof red king crab resulted in declines in accuracy that were an order of magnitude lower for the model-based (&#x394; MSA=0.170) compared to design-based (&#x394; MSA=2.106) estimates (<xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Accuracy of design and model-based survey biomass estimates with and without corner station data for crab stocks/size classes. Time-series on the left show the log accuracy ratio (LAR) over time for design and model-based biomass estimates with and without corner station data. The median LAR across simulation replicates over time is shown as solid/dashed lines, while the 95% distribution across simulation replicates is shown as shaded boundaries. Histograms on the right show the median symmetric accuracy (MSA) across years and simulation replicates. A LAR/MSA value of 0 indicates perfect accuracy while larger values indicate progressively worse accuracy.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1219283-g004.tif"/>
</fig>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Median symmetric accuracy (MSA) of design and model-based biomass estimates for crab species/size classes produced with versus without the corner stations.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">&#xa0;Species/Size class</th>
<th valign="top" align="center">Estimator</th>
<th valign="top" align="center">MSA (corners)</th>
<th valign="top" align="center">MSA (no corners)</th>
<th valign="top" align="center">&#x394; MSA</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Tanner GE103</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.080</td>
<td valign="top" align="center">0.088</td>
<td valign="top" align="center">0.008</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.065</td>
<td valign="top" align="center">0.069</td>
<td valign="top" align="center">0.004</td>
</tr>
<tr>
<td valign="top" align="left">Snow GE95</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.085</td>
<td valign="top" align="center">0.086</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.064</td>
<td valign="top" align="center">0.067</td>
<td valign="top" align="center">0.003</td>
</tr>
<tr>
<td valign="top" align="left">Pribilof blue king GE120</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.314</td>
<td valign="top" align="center">1.075</td>
<td valign="top" align="center">0.761</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.232</td>
<td valign="top" align="center">0.350</td>
<td valign="top" align="center">0.118</td>
</tr>
<tr>
<td valign="top" align="left">Pribilof red king GE120</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.680</td>
<td valign="top" align="center">2.786</td>
<td valign="top" align="center">2.106</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.469</td>
<td valign="top" align="center">0.639</td>
<td valign="top" align="center">0.170</td>
</tr>
<tr>
<td valign="top" align="left">St. Matthew blue king GE90</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.183</td>
<td valign="top" align="center">0.499</td>
<td valign="top" align="center">0.316</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.163</td>
<td valign="top" align="center">0.197</td>
<td valign="top" align="center">0.034</td>
</tr>
<tr>
<td valign="top" align="left">St. Matthew blue King GE105</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.176</td>
<td valign="top" align="center">0.286</td>
<td valign="top" align="center">0.110</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.155</td>
<td valign="top" align="center">0.192</td>
<td valign="top" align="center">0.037</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>&#x394; MSA refers to the difference in MSA between estimates produced with versus without the corner stations.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>For nine of the ten groundfish species we examined, the accuracy of both design and model-based biomass estimates was robust to corner station removal (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>; <xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>). There were modest declines (~13%) in the accuracy of design-based estimates for northern rock sole arising from corner station removal (<xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>). The only species for which substantial declines (~20-30%) in accuracy occurred as a result of removing the corner station data was yellow Irish lord (<xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>). However, the accuracy of both design and model-based biomass estimates for yellow Irish lord was limited even with the corner stations included, and the lowest of any groundfish species that we considered (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>; <xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Accuracy of design and model-based survey biomass estimates with and without corner station data for groundfishes. Time-series on the left show the log accuracy ratio (LAR) over time for design and model-based biomass estimates with and without corner station data. The median LAR across simulation replicates over time is shown as solid/dashed lines, while the 95% distribution across simulation replicates is shown as shaded boundaries. Histograms on the right show the median symmetric accuracy (MSA) across years and simulation replicates. A LAR/MSA value of 0 indicates perfect accuracy while larger values indicate progressively worse accuracy.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1219283-g005.tif"/>
</fig>
<table-wrap id="T4" position="float">
<label>Table&#xa0;4</label>
<caption>
<p>Median symmetric accuracy (MSA) of design and model-based biomass estimates for groundfish species produced with versus without the corner stations.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Species</th>
<th valign="top" align="center">Estimator</th>
<th valign="top" align="center">MSA (corners)</th>
<th valign="top" align="center">MSA (no corners)</th>
<th valign="top" align="center">&#x394; MSA</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Bering skate</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.101</td>
<td valign="top" align="center">0.100</td>
<td valign="top" align="center">-0.001</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.091</td>
<td valign="top" align="center">0.092</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Alaska skate</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.036</td>
<td valign="top" align="center">0.035</td>
<td valign="top" align="center">-0.001</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.033</td>
<td valign="top" align="center">0.033</td>
<td valign="top" align="center">0</td>
</tr>
<tr>
<td valign="top" align="left">Flathead sole</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.071</td>
<td valign="top" align="center">0.073</td>
<td valign="top" align="center">0.002</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.061</td>
<td valign="top" align="center">0.062</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Bering flounder</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.126</td>
<td valign="top" align="center">0.128</td>
<td valign="top" align="center">0.002</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.106</td>
<td valign="top" align="center">0.107</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">Yellowfin sole</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.061</td>
<td valign="top" align="center">0.061</td>
<td valign="top" align="center">0</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.042</td>
<td valign="top" align="center">0.041</td>
<td valign="top" align="center">-0.001</td>
</tr>
<tr>
<td valign="top" align="left">Northern rock sole</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.084</td>
<td valign="top" align="center">0.095</td>
<td valign="top" align="center">.011</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.053</td>
<td valign="top" align="center">0.055</td>
<td valign="top" align="center">0.002</td>
</tr>
<tr>
<td valign="top" align="left">Alaska plaice</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.071</td>
<td valign="top" align="center">0.071</td>
<td valign="top" align="center">0</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.056</td>
<td valign="top" align="center">0.056</td>
<td valign="top" align="center">0</td>
</tr>
<tr>
<td valign="top" align="left">Yellow Irish lord</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.182</td>
<td valign="top" align="center">0.238</td>
<td valign="top" align="center">0.056</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.138</td>
<td valign="top" align="center">0.170</td>
<td valign="top" align="center">0.032</td>
</tr>
<tr>
<td valign="top" align="left">Pacific cod</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.047</td>
<td valign="top" align="center">0.049</td>
<td valign="top" align="center">0.002</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.040</td>
<td valign="top" align="center">0.042</td>
<td valign="top" align="center">0.002</td>
</tr>
<tr>
<td valign="top" align="left">Walleye pollock</td>
<td valign="top" align="center">Design-based</td>
<td valign="top" align="center">0.077</td>
<td valign="top" align="center">0.078</td>
<td valign="top" align="center">0.001</td>
</tr>
<tr>
<td valign="top" align="left">&#xa0;</td>
<td valign="top" align="center">Model-based</td>
<td valign="top" align="center">0.065</td>
<td valign="top" align="center">0.065</td>
<td valign="top" align="center">0</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>&#x394; MSA refers to the difference in MSA between estimates produced with versus without the corner stations.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Stock assessment analyses</title>
<p>Predictions of Tanner crab biomass over time were nearly identical with and without the corner station data for immature and mature males and females (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>). Similarly, estimates of Tanner crab recruitment over time were virtually unchanged (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>), and average recruitment was only marginally greater without the corner station data (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>). Biomass-related reference points were also slightly greater in model runs without the corner station data (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>). The biomass at which maximum spawning potential occurs (B<sub>100</sub>) was estimated to be ~3.5% greater without the corner station data, while estimates of current and projected biomass were ~5% and ~6% greater respectively without the corner station data. Estimates of fishing mortality rates associated with maximum sustained yield (F<sub>MSY</sub>) and the overfishing limit (F<sub>OFL</sub>) were nearly identical regardless of whether the corner station data were included (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>), while estimates of MSY and the overfishing limit (OFL) themselves were ~1.5% and ~4.4% greater without the corner station data. Estimated natural mortality rates for both mature males and females were slightly lower in model runs without the corner station data (<xref ref-type="supplementary-material" rid="SM1">
<bold>Figure S1</bold>
</xref>), likely contributing to the slightly higher estimates of biomass-related reference points.</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Tanner crab stock assessment model output. Left-hand panels show model predictions with (yellow) and without (teal) corner stations for recruits (top row) immature females (second row), mature females (third row), immature males (fourth row), and mature males (fifth row). Barplots on the right show estimated tanner crab reference points produced with (yellow) and without (teal) corner stations. Average recruitment is shown in the top left barplot. Biomass-related reference points are shown in the top-right barplots, including B<sub>100</sub>, current biomass (B<sub>cur</sub>), and projected biomass (B<sub>prj</sub>). Estimates of fishing mortality rates associated with MSY (F<sub>MSY</sub>) and the OFL (F<sub>OFL</sub>) are shown in the bottom barplots. Estimates of MSY and OFL themselves are shown in the bottom right barplot.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1219283-g006.tif"/>
</fig>
<p>Stock assessment model runs for St. Matthew blue king crab completed without the corner station data resulted in biomass estimates that were slightly lower than those produced with the corner station data (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>). However, the biomass time-series of these two scenarios were highly correlated with one another (lag-zero cross correlation coefficient = 0.998), exhibiting nearly identical trends (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>). Estimates of B<sub>MSY</sub> and mature male biomass (MMB) produced without the corner stations were 324 tons and 166 tons lower respectively compared to model fits that included corner station data (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>). However, the estimated ratio of current biomass relative to B<sub>MSY</sub> (B/B<sub>MSY</sub>) was identical regardless of corner station inclusion (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>). Model predictions of recruitment over time were similar regardless of corner station inclusion, though average recruitment was slightly lower in the absence of the corner station data (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>). The estimated overfishing limit from model runs without the corner station data was 0.01 tons lower compared to model runs that included corner station data (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>).</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>St. Matthew blue king crab stock assessment model output. Left-hand panels show model predictions with (yellow) and without (teal) corner stations for recruits (top row) mature male biomass (second row). Barplots on the right show estimated reference points produced with (yellow) and without (teal) corner stations. The biomass at which maximum sustained yield (B<sub>MSY</sub>) is shown in the top left barplot, and the ratio of current biomass to B<sub>MSY</sub> (B/B<sub>MSY</sub>) is shown in the top right barplot. The current level of mature male biomass is shown in the bottom left barplot, and the estimated overfishing limit (OFL) is shown in the bottom right barplot.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-10-1219283-g007.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>The crab stocks for which corner station removal had the greatest impacts on data quality were those whose distributions are concentrated around St. Matthew and the Pribilof Islands. However, it is important to note that the uncertainty in these stocks&#x2019; data products was often substantial even with the corner station data included. The difficulty of producing high quality biomass estimates for the St. Matthew and Pribilof Island king crab stocks may be partially due to these species&#x2019; associations with hard-bottom habitats that the EBS bottom trawl survey cannot sample effectively (<italic>e.g.</italic>, female blue king crab; <xref ref-type="bibr" rid="B65">V&#xf8;lstad, 2012</xref>; <xref ref-type="bibr" rid="B69">Zacher et&#xa0;al., 2020</xref>), and/or limited ranges and low abundance levels (<italic>e.g.</italic>, Pribilof red king crab; <xref ref-type="bibr" rid="B69">Zacher et&#xa0;al., 2020</xref>). At such levels of uncertainty, the contributions these data make to stock assessment models are limited and further declines in precision from corner station removal may not represent a meaningful erosion of their information content. For instance, the annual survey biomass estimates for Pribilof blue and red king crab are considered too variable and uncertain for use in stock assessment models even with the corner station data, such that these assessments rely on <italic>post hoc</italic> smoothing of the survey biomass time series (<xref ref-type="bibr" rid="B53">Stockhausen, 2021a</xref>; <xref ref-type="bibr" rid="B56">Szuwalski, 2022</xref>). In some cases the declines in data precision and accuracy arising from corner station removal were mitigated to some extent by using a model-based versus design-based estimator.</p>
<p>The impacts of corner station removal on both the precision and accuracy of biomass estimates were minor for all groundfish species that we investigated except for yellow Irish lord, which exhibits high concentrations around the Pribilof Islands (<xref ref-type="supplementary-material" rid="SM1">
<bold>Figure S2</bold>
</xref>). The other groundfish species we examined exhibit broader distributions, either due to lack of association with the St. Matthew or Pribilof Island habitat areas or because their distributions are more strongly driven by dynamic environmental conditions than by static habitat characteristics. For instance, Alaska skate are widely distributed throughout the outer EBS shelf area (<xref ref-type="supplementary-material" rid="SM1">
<bold>Figure S3</bold>
</xref>), and yellowfin sole are broadly associated with shallow, inshore areas of the shelf (<xref ref-type="supplementary-material" rid="SM1">
<bold>Figure S4</bold>
</xref>). Moreover, many groundfishes are mobile and alter their distributions in response to environmental conditions. For example, walleye pollock and Pacific cod actively avoid a mass of &lt;2&#xb0; C subsurface water (the cold pool) that occurs on the EBS shelf (<xref ref-type="bibr" rid="B68">Wyllie-Echeverria and Wooster, 1998</xref>; <xref ref-type="bibr" rid="B29">Kotwicki and Lauth, 2013</xref>). As such, static habitat associations are likely a less important component of these species&#x2019; distributions. Importantly, with the exception of yellow Irish lord (which does not have a stock assessment and is of no commercial importance) it is unlikely that the removal of corner stations from the EBS survey would meaningfully impact estimates of stock status or management recommendations for any of the groundfish species we examined.</p>
<p>Our analysis of the effects of corner station removal on stock assessment outcomes showed negligible effects for Tanner crab. The Tanner crab assessment model predictions were nearly identical regardless of whether or not the corner station data were included. The only appreciable changes in model output arising from corner station removal were marginally greater estimates of biomass-related reference points and the OFL, likely attributable to the slightly lower estimates of natural mortality produced without the corner station data. Such minor changes in model predictions and reference points seem unlikely to affect decision-making for this stock, indicating that assessment and management of Tanner crab would be robust to removal of the corner stations.</p>
<p>Not surprisingly, we found somewhat greater effects of corner station removal on stock assessment output for St. Matthew blue king crab. Assessment model predictions made without the corner station for this stock indicated slightly lower mature male biomass and recruitment over time compared to baseline model output that included the corner stations. However, the trends in biomass and recruitment were nearly identical regardless of corner station inclusion. Despite the minor differences between time-series produced with versus without the corner stations, the coherence of trends across data scenarios suggests that ability to detect stock recovery would not be impaired without the corner station data. Similarly, the estimated ratio of current biomass relative to B<sub>MSY</sub> (B/B<sub>MSY</sub>), the reference point used for tracking stock rebuilding, was not affected by removal of the corner station data. As such, despite the substantial effects of corner station removal on survey data products for St. Matthew blue king crab, it does not appear that there is a correspondingly appreciable effect on assessment outcomes for this stock, at least with respect to monitoring stock rebuilding. This may be partially due to the fact that the EBS bottom trawl survey data is not the only source of fisheries-independent data that are used in the St. Matthew blue king crab stock assessment, the other being the Alaska Department of Fish and Game (ADFG) pot survey (<xref ref-type="bibr" rid="B44">Palof, 2022</xref>).</p>
<p>While our analysis focuses on the effects of survey effort reduction on stock assessment inputs (survey data products) and outputs (biological/management reference points), there are other objectives that can be important to fisheries and ecosystem management. For instance, there may be benefits to finer-scale sampling of ecologically important areas &#x2013; beyond those to stock assessment outcomes &#x2013; in improving understanding of crucial habitats. Indeed, the Pribilof Islands are important habitat for many fish, invertebrate, seabird and marine mammal species (<xref ref-type="bibr" rid="B10">Craighead and Oppenheim, 1985</xref>; <xref ref-type="bibr" rid="B18">Gentry, 1998</xref>; <xref ref-type="bibr" rid="B14">Ferrero et&#xa0;al., 2000</xref>) and one of the most productive regions in the Bering Sea (<xref ref-type="bibr" rid="B8">Cooney and Coyle, 1982</xref>; <xref ref-type="bibr" rid="B9">Coyle and Cooney, 1993</xref>), with oceanographic and ecological dynamics that are distinct from the rest of the EBS shelf (<xref ref-type="bibr" rid="B23">Hunt et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B6">Ciannelli et&#xa0;al., 2004</xref>). However, it is difficult to quantitatively assess the value of finer-scale survey observations for understanding ecological and oceanographic processes. Nonetheless, objectives beyond utility to stock assessment such as ecosystem considerations may be important to managers, and it is important to note that such concerns cannot be assessed using the approach we present here.</p>
<p>We found that the stocks for which corner station removal caused substantial declines in data quality were also those with the lowest quality data to begin with. This result suggests that while sampling the corner stations may improve data quality for some stocks, the extent of improvement may not be sufficient to produce satisfactorily informative data products. The limited data quality for Pribilof and St. Matthew king crab and yellow Irish lord may simply reflect the fact that the EBS bottom trawl survey is not optimized for sampling certain species, a challenge that naturally arises in fisheries-independent survey designs with multispecies objectives (<xref ref-type="bibr" rid="B7">Cochran, 1977</xref>; <xref ref-type="bibr" rid="B20">Godo, 1994</xref>; <xref ref-type="bibr" rid="B41">Oyafuso et&#xa0;al., 2021</xref>). Data quality for these stocks may be more effectively improved by alternative approaches beyond increased bottom trawl sampling density, such as expanded collection of other forms of fisheries-independent data (<italic>e.g.</italic>, pot surveys for king crab; <xref ref-type="bibr" rid="B19">Gish and Vanek, 2010</xref>; <xref ref-type="bibr" rid="B44">Palof, 2022</xref>), or building capacity for generating and using model-based indices (<xref ref-type="bibr" rid="B61">Thorson, 2019</xref>). While we found that model-based biomass estimates were more robust to corner station removal than their design-based counterparts, adequate model performance was difficult to achieve in some cases (<italic>e.g.</italic>, St. Matthew blue king crab, Pribilof red and blue king crab) and stock assessments for Bering Sea crab still rely entirely on design-based abundance and compositional estimates (<xref ref-type="bibr" rid="B52">Stockhausen, 2014</xref>; <xref ref-type="bibr" rid="B57">Szuwalski et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B58">Szuwalski and Turnock, 2016</xref>; <xref ref-type="bibr" rid="B43">Palof et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B53">Stockhausen, 2021a</xref>; <xref ref-type="bibr" rid="B54">Stockhausen, 2021b</xref>; <xref ref-type="bibr" rid="B56">Szuwalski, 2022</xref>). Continued innovation in model-based index production methods for crab, as well as expanded capacity to use model-based indices in crab assessment models may lead to improved inference on stock status and mitigate adverse consequences from potential survey effort reduction. Model-based indices can also incorporate environmental covariates to facilitate extrapolation to unsampled areas (<xref ref-type="bibr" rid="B61">Thorson, 2019</xref>; <xref ref-type="bibr" rid="B38">O'Leary et&#xa0;al., 2022</xref>) which may also help offset the impacts of reductions in sampling effort.</p>
<p>While stock assessment models are sensitive to the precision of input data, reductions in data quality beyond a certain point may not exert a corresponding effect on their information content. This is due to the integrated design of many contemporary stock assessment models (including those for Tanner and St. Matthew blue king crab), in which multiple data sources contribute to parameter estimates via a joint likelihood (<xref ref-type="bibr" rid="B16">Fournier and Archibald, 1982</xref>; <xref ref-type="bibr" rid="B33">Maunder and Punt, 2013</xref>). In such integrated models, the greater the uncertainty in one data source, the less information it will contribute to the model, whether by explicit data weighting procedures (<italic>e.g.</italic>, <xref ref-type="bibr" rid="B17">Francis, 2011</xref>) or via the joint likelihood itself (<italic>e.g.</italic>, <xref ref-type="bibr" rid="B12">DeFilippo et&#xa0;al., 2021</xref>). In our analysis, this effect is best demonstrated by St. Matthew blue king crab, for which corner station removal had a minor effect on assessment model outputs despite exerting a large influence on the input survey data products. This result emphasizes the importance of propagating the impacts of survey effort reduction through stock assessment models to understand the impacts of proposed survey design changes. Indeed, contemporary stock assessment models are complex, with multiple contributing (and often conflicting) sources of information that can interact to affect parameter estimates in unpredictable ways (<xref ref-type="bibr" rid="B33">Maunder and Punt, 2013</xref>; <xref ref-type="bibr" rid="B26">Ichinokawa et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B45">Peterson et&#xa0;al., 2021</xref>).</p>
<p>The approach we demonstrate here is flexible and can be used to evaluate the potential effects of removing any set of existing stations from a systematic fisheries-independent survey. The corner stations we considered for removal in our analysis were selected because they represent an area of higher sampling density than the rest of the survey grid and because the stocks that motivated their implementation have remained below harvestable levels for a prolonged period of time. However, there are other situations in which a set of stations may be identified <italic>a priori</italic> as candidates for removal. For instance, a particular area within a survey&#x2019;s boundaries may be under consideration for wind energy development or MPA designation, or become inaccessible due to political (<italic>e.g.</italic>, international strife) or logistical (<italic>e.g.</italic>, prohibitive fuel costs for reaching remote areas) reasons (<xref ref-type="bibr" rid="B15">Field et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B24">ICES, 2020</xref>; <xref ref-type="bibr" rid="B35">Methratta et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B22">Hare et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B25">ICES, 2023</xref>). In such scenarios, scientists could follow the general approach outlined here to quantify the expected impacts of removing the stations in question on data quality and stock assessment outcomes. These results could be used to understand and predict the impacts of the effort reduction on scientific inference, and/or be used to develop calibration factors to correct for these effects in the survey data time series. Conversely, there may be a need to reduce sampling effort more generally (<italic>e.g.</italic>, due to budgetary and staffing shortages), without a specific subset of stations indicated as leading options for removal (<xref ref-type="bibr" rid="B24">ICES, 2020</xref>; <xref ref-type="bibr" rid="B25">ICES, 2023</xref>). In response to such a need to reduce sampling effort without obvious candidates for station removal, scientists may use our approach to evaluate a range of alternative scenarios (<italic>e.g.</italic>, thinning the sampling density of the entire survey grid by a specified fraction, removing stations from various areas based on expert opinion), and compare the impacts of each effort reduction scenario on data quality and stock assessment output to inform decision-making.</p>
<p>While evaluating the effects of station removal on all survey data products and stock assessment outcomes is optimal, scientists may prioritize certain stocks to focus effort on (<italic>e.g.</italic>, crab in our case study) and conduct more limited investigations on others for which expected impacts are lower (<italic>e.g.</italic>, groundfish in our case study). As shifting ecosystems and species distributions place growing demands on fisheries-independent monitoring amid rising costs and limited institutional budgets, unavoidable reductions in survey effort are likely to become increasingly common (<xref ref-type="bibr" rid="B24">ICES, 2020</xref>; <xref ref-type="bibr" rid="B25">ICES, 2023</xref>). Consequently, it is imperative that scientists can evaluate the effects of survey effort reductions to inform decision-making on such actions and understand their consequences. The analysis we present here demonstrates a generic set of steps for evaluating reductions in sampling density for systematic survey designs that can be applied in similar situations elsewhere to provide a quantitative basis for decision-making.</p>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>Publicly available datasets were analyzed in this study. This data can be found here: <ext-link ext-link-type="uri" xlink:href="https://www.fisheries.noaa.gov/foss/f?p=215%3A28">https://www.fisheries.noaa.gov/foss/f?p=215%3A28</ext-link>.</p>
</sec>
<sec id="s6" sec-type="ethics-statement">
<title>Ethics statement</title>
<p>The manuscript presents research on animals that do not require ethical approval for their study.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>SK and LD conceived the study. LD designed and executed the analyses and wrote the paper. WS and KP conducted and helped interpret the stock assessment model runs. JR and ML contributed crab data and input on their use and interpretation. LB, LD, and SK contributed input on interpreting study results. All authors contributed to manuscript writing. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s8" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s9" sec-type="disclaimer">
<title>Publisher's note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s10" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmars.2023.1219283/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmars.2023.1219283/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet_1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brown</surname> <given-names>R. W.</given-names>
</name>
<name>
<surname>Fogarty</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Legault</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Miller</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Nordahl</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Politis</surname> <given-names>P.</given-names>
</name>
<etal/>
</person-group>. (<year>2007</year>). <article-title>Survey transition and calibration of bottom trawl surveys along the northeastern continental shelf of the United States</article-title>. <source>Int. Council Explor. Sea CM.</source>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cadigan</surname> <given-names>N. G.</given-names>
</name>
<name>
<surname>Dowden</surname> <given-names>J. J.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Statistical inference about the relative efficiency of a new survey protocol, based on paired-tow survey calibration data</article-title>. <source>Fishery Bull.</source> <volume>108</volume> (<issue>1</issue>).</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cadigan</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Yin</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Beno&#xee;t</surname> <given-names>H. P.</given-names>
</name>
<name>
<surname>Walsh</surname> <given-names>S. J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A nonparametric-monotone regression model and robust estimation for paired-tow bottom-trawl survey comparative fishing data</article-title>. <source>Fisheries Res.</source> <volume>254</volume>, <fpage>106422</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fishres.2022.106422</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Richards</surname> <given-names>R. A.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Spatiotemporal index standardization improves the stock assessment of northern shrimp in the Gulf of Maine</article-title>. <source>Can. J. Fisheries Aquat. Sci.</source> <volume>74</volume> (<issue>11</issue>), <fpage>1781</fpage>&#x2013;<lpage>1793</lpage>. doi: <pub-id pub-id-type="doi">10.1139/cjfas-2016-0137</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Stergiou</surname> <given-names>K. I.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Impacts of data quantity on fisheries stock assessment</article-title>. <source>Aquat. Sci.</source> <volume>65</volume> (<issue>1</issue>), <fpage>92</fpage>&#x2013;<lpage>98</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s000270300008</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ciannelli</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Robson</surname> <given-names>B. W.</given-names>
</name>
<name>
<surname>Francis</surname> <given-names>R. C.</given-names>
</name>
<name>
<surname>Aydin</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Brodeur</surname> <given-names>R. D.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Boundaries of open marine ecosystems: an application to the Pribilof Archipelago, southeast Bering Sea</article-title>. <source>Ecol. Appl.</source> <volume>14</volume> (<issue>3</issue>), <fpage>942</fpage>&#x2013;<lpage>953</lpage>. doi: <pub-id pub-id-type="doi">10.1890/03-5016</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Cochran</surname> <given-names>W. G.</given-names>
</name>
</person-group> (<year>1977</year>). <source>Sampling techniques</source> (Hoboken, New Jersey: <publisher-name>John Wiley &amp; Sons</publisher-name>).</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cooney</surname> <given-names>R. T.</given-names>
</name>
<name>
<surname>Coyle</surname> <given-names>K. O.</given-names>
</name>
</person-group> (<year>1982</year>). <article-title>Trophic implications of cross-shelf copepod distributions in the southeastern Bering Sea</article-title>. <source>Mar. Biol.</source> <volume>70</volume> (<issue>2</issue>), <fpage>187</fpage>&#x2013;<lpage>196</lpage>. doi: <pub-id pub-id-type="doi">10.1007/BF00397684</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Coyle</surname> <given-names>K. O.</given-names>
</name>
<name>
<surname>Cooney</surname> <given-names>R. T.</given-names>
</name>
</person-group> (<year>1993</year>). <article-title>Water column sound scattering and hydrography around the Pribilof Islands, Bering Sea</article-title>. <source>Continental Shelf Res.</source> <volume>13</volume> (<issue>7</issue>), <fpage>803</fpage>&#x2013;<lpage>827</lpage>. doi: <pub-id pub-id-type="doi">10.1016/0278-4343(93)90028-V</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Craighead</surname> <given-names>F. L.</given-names>
</name>
<name>
<surname>Oppenheim</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>1985</year>). &#x201c;<article-title>Population estimates and temporal trends of Pribilof Island seabirds</article-title>,&#x201d; in <source>Outer Continental Shelf Environmental Assessment Program (OCSEAP)  (Anchorage Alaska:U.S. Department of Commerce) Final Report</source>, vol. <volume>30</volume>. , <fpage>307</fpage>&#x2013;<lpage>356</lpage>.</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>DeFilippo</surname> <given-names>L. B.</given-names>
</name>
<name>
<surname>Buehrens</surname> <given-names>T. W.</given-names>
</name>
<name>
<surname>Scheuerell</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Kendall</surname> <given-names>N. W.</given-names>
</name>
<name>
<surname>Schindler</surname> <given-names>D. E.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Improving short-term recruitment forecasts for coho salmon using a spatiotemporal integrated population model</article-title>. <source>Fisheries Res.</source> <volume>242</volume>, <fpage>106014</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fishres.2021.106014</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>DeFilippo</surname> <given-names>L. B.</given-names>
</name>
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>O'Leary</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Kotwicki</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Hoff</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Ianelli</surname> <given-names>J. N.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Characterizing dominant patterns of spatiotemporal variation for a transboundary groundfish assemblage</article-title>. <source>Fisheries Oceanography.</source> doi: <pub-id pub-id-type="doi">10.1111/fog.12651</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dulvy</surname> <given-names>N. K.</given-names>
</name>
<name>
<surname>Rogers</surname> <given-names>S. I.</given-names>
</name>
<name>
<surname>Jennings</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Stelzenm&#xfc;ller</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Dye</surname> <given-names>S. R.</given-names>
</name>
<name>
<surname>Skjoldal</surname> <given-names>H. R.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Climate change and deepening of the North Sea fish assemblage: a biotic indicator of warming seas</article-title>. <source>J. Appl. Ecol.</source> <volume>45</volume> (<issue>4</issue>), <fpage>1029</fpage>&#x2013;<lpage>1039</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1365-2664.2008.01488.x</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ferrero</surname> <given-names>R. C.</given-names>
</name>
<name>
<surname>DeMaster</surname> <given-names>D. P.</given-names>
</name>
<name>
<surname>Hill</surname> <given-names>P. S.</given-names>
</name>
<name>
<surname>Muto</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Lopez</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Brix</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Alaska marine mammal stock assessment</article-title>. <source>NOAA technical memorandum NMFS-AFSC</source>. <fpage>p 119</fpage>.</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Field</surname> <given-names>J. C.</given-names>
</name>
<name>
<surname>Punt</surname> <given-names>A. E.</given-names>
</name>
<name>
<surname>Methot</surname> <given-names>R. D.</given-names>
</name>
<name>
<surname>Thomson</surname> <given-names>C. J.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Does MPA mean &#x2018;Major Problem for Assessments'? Considering the consequences of place-based management systems</article-title>. <source>Fish Fisheries</source> <volume>7</volume> (<issue>4</issue>), <fpage>284</fpage>&#x2013;<lpage>302</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1467-2979.2006.00226.x</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fournier</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Archibald</surname> <given-names>C. P.</given-names>
</name>
</person-group> (<year>1982</year>). <article-title>A general theory for analyzing catch at age data</article-title>. <source>Can. J. Fisheries Aquat. Sci.</source> <volume>39</volume> (<issue>8</issue>), <fpage>1195</fpage>&#x2013;<lpage>1207</lpage>. doi: <pub-id pub-id-type="doi">10.1139/f82-157</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Francis</surname> <given-names>R. C.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Data weighting in statistical fisheries stock assessment models</article-title>. <source>Can. J. Fisheries Aquat. Sci.</source> <volume>68</volume> (<issue>6</issue>), <fpage>1124</fpage>&#x2013;<lpage>1138</lpage>. doi: <pub-id pub-id-type="doi">10.1139/f2011-025</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Gentry</surname> <given-names>R. L.</given-names>
</name>
</person-group> (<year>1998</year>). <source>Behavior and ecology of the northern fur seal</source> Vol. <volume>382</volume> (Princeton, New Jersey: <publisher-name>Princeton University Press</publisher-name>).</citation>
</ref>
<ref id="B19">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Gish</surname> <given-names>R. K.</given-names>
</name>
<name>
<surname>Vanek</surname> <given-names>V. A.</given-names>
</name>
</person-group> (<year>2010</year>). <source>Project operational plan for the 2010 St. Matthew Island blue king crab survey</source> (<publisher-loc>Kodiak</publisher-loc>: <publisher-name>Alaska Department of Fish and Game, Regional Information Report 4K10-12</publisher-name>).</citation>
</ref>
<ref id="B20">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Godo</surname> <given-names>O. R.</given-names>
</name>
</person-group> (<year>1994</year>). &#x201c;<article-title>Factors affecting the reliability of groundfish abundance estimates from bottom trawl surveys</article-title>,&#x201d; in <source>Marine fish behavior in capture and abundance estimation</source>. (Hoboken New Jersey:Wiley and Sons)</citation>
</ref>
<ref id="B21">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Gunderson</surname> <given-names>D. R.</given-names>
</name>
</person-group> (<year>1993</year>). <source>Surveys of fisheries resources</source> (Hoboken, New Jersey: <publisher-name>John Wiley &amp; Sons</publisher-name>).</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hare</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Blythe</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Ford</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Hooker</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Jensen</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Lipsky</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>NOAA fisheries and BOEM federal survey mitigation implementation strategy - northeast U.S. Region</article-title>. <source>NOAA Technical Memorandum</source> <volume>292</volume>, <fpage>33</fpage>. doi: <pub-id pub-id-type="doi">10.25923/jqse-x746</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hunt</surname> <given-names>G. L.</given-names>
<suffix>Jr.</suffix>
</name>
<name>
<surname>Stabeno</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Strom</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Napp</surname> <given-names>J. M.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Patterns of spatial and temporal variation in the marine ecosystem of the southeastern Bering Sea, with special reference to the Pribilof Domain</article-title>. <source>Deep Sea Res. Part II: Topical Stud. Oceanography</source> <volume>55</volume> (<issue>16-17</issue>), <fpage>1919</fpage>&#x2013;<lpage>1944</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.dsr2.2008.04.032</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<collab>ICES</collab>
</person-group> (<year>2020</year>). <article-title>Workshop on unavoidable survey effort reduction (WKUSER)</article-title>. <source>ICES Sci. Rep.</source> <volume>2</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.17895/ices.pub.7453</pub-id>. 72.92pp(ICES, 2020).</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<collab>ICES</collab>
</person-group> (<year>2023</year>). <article-title>Workshop on unavoidable survey effort reduction 2 (WKUSER2)</article-title>. <source>ICES Sci. Rep</source>. doi:&#xa0;<pub-id pub-id-type="doi">10.17895/ices.pub.22086845.v1</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ichinokawa</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Okamura</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Takeuchi</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Data conflict caused by model mis-specification of selectivity in an integrated stock assessment model and its potential effects on stock status estimation</article-title>. <source>Fisheries Res.</source> <volume>158</volume>, <fpage>147</fpage>&#x2013;<lpage>157</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fishres.2014.02.003</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Jadamec</surname> <given-names>L. S.</given-names>
</name>
<name>
<surname>Donaldson</surname> <given-names>W. E.</given-names>
</name>
<name>
<surname>Cullenberg</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>1999</year>). <source>Biological field techniques for Chionoecetes crabs</source> (Fairbanks, AK: <publisher-name>University of Alaska Sea Grant Program Report 99-02</publisher-name>). 80 p.</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kimura</surname> <given-names>D. K.</given-names>
</name>
<name>
<surname>Somerton</surname> <given-names>D. A.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Review of statistical aspects of survey sampling for marine fisheries</article-title>. <source>Rev. Fisheries Sci.</source> <volume>14</volume> (<issue>3</issue>), <fpage>245</fpage>&#x2013;<lpage>283</lpage>. doi: <pub-id pub-id-type="doi">10.1080/10641260600621761</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kotwicki</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Ianelli</surname> <given-names>J. N.</given-names>
</name>
<name>
<surname>Punt</surname> <given-names>A. E.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Correcting density-dependent effects in abundance estimates from bottom-trawl surveys</article-title>. <source>ICES J. Mar. Sci.</source> <volume>71</volume> (<issue>5</issue>), <fpage>1107</fpage>&#x2013;<lpage>1116</lpage>. doi: <pub-id pub-id-type="doi">10.1093/icesjms/fst208</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kotwicki</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Lauth</surname> <given-names>R. R.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Detecting temporal trends and environmentally-driven changes in the spatial distribution of bottom fishes and crabs on the eastern Bering Sea shelf</article-title>. <source>Deep Sea Res. Part II: Topical Stud. Oceanography</source> <volume>94</volume>, <fpage>231</fpage>&#x2013;<lpage>243</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.dsr2.2013.03.017</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kotwicki</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Lauth</surname> <given-names>R. R.</given-names>
</name>
<name>
<surname>Williams</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Goodman</surname> <given-names>S. E.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Selectivity ratio: a useful tool for comparing size selectivity of multiple survey gears</article-title>. <source>Fisheries Res.</source> <volume>191</volume>, <fpage>76</fpage>&#x2013;<lpage>86</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fishres.2017.02.012</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lauth</surname> <given-names>R. R.</given-names>
</name>
<name>
<surname>Dawson</surname> <given-names>E. J.</given-names>
</name>
<name>
<surname>Conner</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Results of the 2017 eastern and northern bering sea continental shelf bottom trawl survey of groundfish and invertebrate fauna</article-title>. <source>Natl. Mar. Fish. Serv.</source>, <volume>396</volume>, <fpage>1</fpage>&#x2013;<lpage>265</lpage>. doi: <pub-id pub-id-type="doi">10.25923/h118-nw41</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maunder</surname> <given-names>M. N.</given-names>
</name>
<name>
<surname>Punt</surname> <given-names>A. E.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>A review of integrated analysis in fisheries stock assessment</article-title>. <source>Fisheries Res.</source> <volume>142</volume>, <fpage>61</fpage>&#x2013;<lpage>74</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fishres.2012.07.025</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maureaud</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Frelat</surname> <given-names>R.</given-names>
</name>
<name>
<surname>P&#xe9;cuchet</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Shackell</surname> <given-names>N.</given-names>
</name>
<name>
<surname>M&#xe9;rigot</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Pinsky</surname> <given-names>M. L.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Are we ready to track climate-driven shifts in marine species across international boundaries?-A global survey of scientific bottom trawl data</article-title>. <source>Global Change Biol.</source> <volume>27</volume> (<issue>2</issue>), <fpage>220</fpage>&#x2013;<lpage>236</lpage>. doi: <pub-id pub-id-type="doi">10.1111/gcb.15404</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Methratta</surname> <given-names>E. T.</given-names>
</name>
<name>
<surname>Hawkins</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Hooker</surname> <given-names>B. R.</given-names>
</name>
<name>
<surname>Lipsky</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Hare</surname> <given-names>J. A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Offshore wind development in the Northeast US Shelf large marine ecosystem</article-title>. <source>Oceanography</source> <volume>33</volume> (<issue>4</issue>), <fpage>16</fpage>&#x2013;<lpage>27</lpage>. doi: <pub-id pub-id-type="doi">10.5670/oceanog.2020.402</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Miller</surname> <given-names>T. J.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>A comparison of hierarchical models for relative catch efficiency based on paired-gear data for US Northwest Atlantic fish stocks</article-title>. <source>Can. J. Fisheries Aquat. Sci.</source> <volume>70</volume> (<issue>9</issue>), <fpage>1306</fpage>&#x2013;<lpage>1316</lpage>. doi: <pub-id pub-id-type="doi">10.1139/cjfas-2013-0136</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Morley</surname> <given-names>S. K.</given-names>
</name>
<name>
<surname>Brito</surname> <given-names>T. V.</given-names>
</name>
<name>
<surname>Welling</surname> <given-names>D. T.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Measures of model performance based on the log accuracy ratio</article-title>. <source>Space Weather</source> <volume>16</volume> (<issue>1</issue>), <fpage>69</fpage>&#x2013;<lpage>88</lpage>. doi: <pub-id pub-id-type="doi">10.1002/2017SW001669</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>O'Leary</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>DeFilippo</surname> <given-names>L. B.</given-names>
</name>
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Kotwicki</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Hoff</surname> <given-names>G. R.</given-names>
</name>
<name>
<surname>Kulik</surname> <given-names>V. V.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Understanding transboundary stocks&#x2019; availability by combining multiple fisheries-independent surveys and oceanographic conditions in spatiotemporal models</article-title>. <source>ICES J. Mar. Sci.</source> <volume>79</volume> (<issue>4</issue>), <fpage>1063</fpage>&#x2013;<lpage>1074</lpage>. doi: <pub-id pub-id-type="doi">10.1093/icesjms/fsac046</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>O&#x2019;Leary</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Ianelli</surname> <given-names>J. N.</given-names>
</name>
<name>
<surname>Kotwicki</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Adapting to climate-driven distribution shifts using model-based indices and age composition from multiple surveys in the walleye pollock (Gadus chalcogrammus) stock assessment</article-title>. <source>Fisheries Oceanography</source> <volume>29</volume> (<issue>6</issue>), <fpage>541</fpage>&#x2013;<lpage>557</lpage>. doi: <pub-id pub-id-type="doi">10.1111/fog.12494</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ono</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Licandeo</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Muradian</surname> <given-names>M. L.</given-names>
</name>
<name>
<surname>Cunningham</surname> <given-names>C. J.</given-names>
</name>
<name>
<surname>Anderson</surname> <given-names>S. C.</given-names>
</name>
<name>
<surname>Hurtado-Ferro</surname> <given-names>F.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>The importance of length and age composition data in statistical age-structured models for marine species</article-title>. <source>ICES J. Mar. Sci.</source> <volume>72</volume> (<issue>1</issue>), <fpage>31</fpage>&#x2013;<lpage>43</lpage>. doi: <pub-id pub-id-type="doi">10.1093/icesjms/fsu007</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oyafuso</surname> <given-names>Z. S.</given-names>
</name>
<name>
<surname>Barnett</surname> <given-names>L. A.</given-names>
</name>
<name>
<surname>Kotwicki</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Incorporating spatiotemporal variability in multispecies survey design optimization addresses trade-offs in uncertainty</article-title>. <source>ICES J. Mar. Sci.</source> <volume>78</volume> (<issue>4</issue>), <fpage>1288</fpage>&#x2013;<lpage>1300</lpage>. doi: <pub-id pub-id-type="doi">10.1093/icesjms/fsab038</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Oyafuso</surname> <given-names>Z. S.</given-names>
</name>
<name>
<surname>Barnett</surname> <given-names>L. A. K.</given-names>
</name>
<name>
<surname>Siple</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Kotwicki</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2022</year>). <source>A flexible approach to optimizing the Gulf of Alaska groundfish bottom trawl survey design for abundance estimation</source> (Seattle, WA: <publisher-name>U.S. Dep. Commer., NOAA Tech. Memo. NMFS-AFSC-434</publisher-name>). 142 p.</citation>
</ref>
<ref id="B44">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Palof</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2022</year>). <source>Saint matthew island blue king crab stock assesssment</source> (<publisher-name>Anchorage, AK: Alaska Department of Fish and Game</publisher-name>).</citation>
</ref>
<ref id="B43">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Palof</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Zheng</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Ianelli</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2020</year>). &#x201c;<article-title>Saint matthew island blue king crab stock assessment, (2020)</article-title>,&#x201d; in <source>Stock assessment and fishery evaluation report for the king and Tanner crab fisheries of the Bering Sea and Aleutian Islands regions</source> (<publisher-loc>Alaska</publisher-loc>: <publisher-name>North Pacific Fishery Management Council, Anchorage</publisher-name>). 2020 Final Crab SAFE.</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peterson</surname> <given-names>C. D.</given-names>
</name>
<name>
<surname>Courtney</surname> <given-names>D. L.</given-names>
</name>
<name>
<surname>Cort&#xe9;s</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Latour</surname> <given-names>R. J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Reconciling conflicting survey indices of abundance prior to stock assessment</article-title>. <source>ICES J. Mar. Sci.</source> <volume>78</volume> (<issue>9</issue>), <fpage>3101</fpage>&#x2013;<lpage>3120</lpage>. doi: <pub-id pub-id-type="doi">10.1093/icesjms/fsab179</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pinsky</surname> <given-names>M. L.</given-names>
</name>
<name>
<surname>Mantua</surname> <given-names>N. J.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Emerging adaptation approaches for climate-ready fisheries management</article-title>. <source>Oceanography</source> <volume>27</volume> (<issue>4</issue>), <fpage>146</fpage>&#x2013;<lpage>159</lpage>. doi: <pub-id pub-id-type="doi">10.5670/oceanog.2014.93</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pope</surname> <given-names>J. G.</given-names>
</name>
<name>
<surname>Shepherd</surname> <given-names>J. G.</given-names>
</name>
</person-group> (<year>1985</year>). <article-title>A comparison of the performance of various methods for tuning VPAs using effort data</article-title>. <source>Ices J. Mar. Sci.</source> <volume>42</volume> (<issue>2</issue>), <fpage>129</fpage>&#x2013;<lpage>151</lpage>. doi: <pub-id pub-id-type="doi">10.1093/icesjms/42.2.129</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Quinn</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Deriso</surname> <given-names>R. B.</given-names>
</name>
</person-group> (<year>1999</year>). <source>Quantitative fish dynamics</source> (<publisher-name>New York, New York: Oxford university Press</publisher-name>).</citation>
</ref>
<ref id="B49">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>R Core Team</collab>
</person-group> (<year>2021</year>). <source>R: A language and environment for statistical computing</source> (<publisher-loc>Vienna, Austria</publisher-loc>: <publisher-name>R foundation for Statistical Computing</publisher-name>). Available at: <uri xlink:href="https://www.R-project.org/">https://www.R-project.org/</uri>.</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Santora</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Rogers</surname> <given-names>T. L.</given-names>
</name>
<name>
<surname>Cimino</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Sakuma</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Hanson</surname> <given-names>K. D.</given-names>
</name>
<name>
<surname>Dick</surname> <given-names>E. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Diverse integrated ecosystem approach overcomes pandemic-related fisheries monitoring challenges</article-title>. <source>Nat. Commun.</source> <volume>12</volume> (<issue>1</issue>), <fpage>1</fpage>&#x2013;<lpage>10</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41467-021-26484-5</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Stauffer</surname> <given-names>G. D.</given-names>
</name>
</person-group> (<year>2004</year>). <source>NOAA protocols for groundfish bottom trawl surveys of the nation's fishery resources, March 16, 2003</source>. (Seattle WA:U.S. Department of Commerce)</citation>
</ref>
<ref id="B52">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Stockhausen</surname> <given-names>W. T.</given-names>
</name>
</person-group> (<year>2014</year>). <source>2014 stock assessment and fishery evaluation report for the tanner crab fisheries of the bering sea and aleutian islands regions</source>. (Anchorage AK:North Pacific Fishery Management Council)</citation>
</ref>
<ref id="B53">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Stockhausen</surname> <given-names>W. T.</given-names>
</name>
</person-group> (<year>2021</year>a). <source>2021 stock assessment and fishery evaluation report for the pribilof islands blue king crab fisheries of the bering sea and aleutian islands regions</source>. (Anchorage AK:North Pacific Fishery Management Council)</citation>
</ref>
<ref id="B54">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Stockhausen</surname> <given-names>W. T.</given-names>
</name>
</person-group> (<year>2021</year>b). <source>2021 stock assessment and fishery evaluation report for the tanner crab fisheries of the bering sea and aleutian islands regions</source> (Anchorage, AK: <publisher-name>North Pacific Fishery Management Council</publisher-name>). Available at: <uri xlink:href="https://meetings.npfmc.org/CommentReview/DownloadFile?p=55d56b05-4f74-4c5d-9372-28ced4a2bfc2.pdf&amp;fileName=C1%20Tanner%20Crab%20SAFE.pdf">https://meetings.npfmc.org/CommentReview/DownloadFile?p=55d56b05-4f74-4c5d-9372-28ced4a2bfc2.pdf&amp;fileName=C1%20Tanner%20Crab%20SAFE.pdf</uri>.</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stompe</surname> <given-names>D. K.</given-names>
</name>
<name>
<surname>Moyle</surname> <given-names>P. B.</given-names>
</name>
<name>
<surname>Kruger</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Durand</surname> <given-names>J. R.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Comparing and integrating fish surveys in the San Francisco Estuary: why diverse long-term monitoring programs are important</article-title>. <source>San Francisco Estuary Watershed Sci.</source> <volume>18</volume> (<issue>2</issue>). doi: <pub-id pub-id-type="doi">10.15447/sfews.2020v18iss2art4</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Szuwalski</surname> <given-names>C. S.</given-names>
</name>
</person-group> (<year>2022</year>) <source>2022 assessment for Pribilof Islands red king crab</source>. Available at: <uri xlink:href="https://meetings.npfmc.org/CommentReview/DownloadFile?p=4783b2a2-e8a2-4efe-bdfb-37f89216c1b3.pdf&amp;fileName=PIRKC%20SAFE%202022%20Final.pdf">https://meetings.npfmc.org/CommentReview/DownloadFile?p=4783b2a2-e8a2-4efe-bdfb-37f89216c1b3.pdf&amp;fileName=PIRKC%20SAFE%202022%20Final.pdf</uri>.</citation>
</ref>
<ref id="B57">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Szuwalski</surname> <given-names>C. S.</given-names>
</name>
<name>
<surname>Foy</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Turnock</surname> <given-names>B. J.</given-names>
</name>
</person-group> (<year>2014</year>). &#x201c;<article-title>2014 Stock assessment and fishery evaluation report for the Pribilof Islands red king crab fishery of the Bering Sea and Aleutian Islands Regions</article-title>,&#x201d; in <source>Stock assessment and fishery evaluation report for the king and tanner crab fisheries of the bering sea and aleutian islands regions</source> (<publisher-loc>605 W. 4th Avenue, no. 306, Anchorage, AK 99501</publisher-loc>: <publisher-name>North Pacific Fishery Management Council</publisher-name>), <fpage>546</fpage>&#x2013;<lpage>605</lpage>.</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Szuwalski</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Turnock</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>A stock assessment for eastern Bering Sea snow crab</article-title> <fpage>167</fpage>&#x2013;<lpage>251</lpage>. Stock Assessment and Fishery Evaluation Report for the King and Tanner Crab Fisheries of the Bering Sea and Aleutian Islands Regions. <source>North Pacific Fishery Management Council, Stock Assessment and Fishery Evaluation Report</source></citation>
</ref>
<ref id="B61">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Guidance for decisions using the Vector Autoregressive Spatio-Temporal (VAST) package in stock, ecosystem, habitat and climate assessments</article-title>. <source>Fisheries Res.</source> <volume>210</volume>, <fpage>143</fpage>&#x2013;<lpage>161</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fishres.2018.10.013</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Barbeaux</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Goethel</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Kearney</surname> <given-names>K. A.</given-names>
</name>
<name>
<surname>Laman</surname> <given-names>E. A.</given-names>
</name>
<name>
<surname>Nielsen</surname> <given-names>J. K.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Estimating fine-scale movement rates and habitat preferences using multiple data sources</article-title>. <source>Fish Fisheries</source> <volume>22</volume> (<issue>6</issue>), <fpage>1359</fpage>&#x2013;<lpage>1376</lpage>. doi: <pub-id pub-id-type="doi">10.1111/faf.12592</pub-id>
</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Barnett</surname> <given-names>L. A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Comparing estimates of abundance trends and distribution shifts using single-and multispecies models of fishes and biogenic habitat</article-title>. <source>ICES J. Mar. Sci.</source> <volume>74</volume> (<issue>5</issue>), <fpage>1311</fpage>&#x2013;<lpage>1321</lpage>. doi: <pub-id pub-id-type="doi">10.1093/icesjms/fsw193</pub-id>
</citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Clarke</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Stewart</surname> <given-names>I. J.</given-names>
</name>
<name>
<surname>Punt</surname> <given-names>A. E.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>The implications of spatially varying catchability on bottom trawl surveys of fish abundance: a proposed solution involving underwater vehicles</article-title>. <source>Can. J. Fisheries Aquat. Sci.</source> <volume>70</volume> (<issue>2</issue>), <fpage>294</fpage>&#x2013;<lpage>306</lpage>. doi: <pub-id pub-id-type="doi">10.1139/cjfas-2012-0330</pub-id>
</citation>
</ref>
<ref id="B62">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Haltuch</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Spatiotemporal analysis of compositional data: increased precision and improved workflow using model-based inputs to stock assessment</article-title>. <source>Can. J. Fisheries Aquat. Sci.</source> <volume>76</volume> (<issue>3</issue>), <fpage>401</fpage>&#x2013;<lpage>414</lpage>. doi: <pub-id pub-id-type="doi">10.1139/cjfas-2018-0015</pub-id>
</citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Shelton</surname> <given-names>A. O.</given-names>
</name>
<name>
<surname>Ward</surname> <given-names>E. J.</given-names>
</name>
<name>
<surname>Skaug</surname> <given-names>H. J.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Geostatistical delta-generalized linear mixed models improve precision for estimated abundance indices for West Coast groundfishes</article-title>. <source>ICES J. Mar. Sci.</source> <volume>72</volume> (<issue>5</issue>), <fpage>1297</fpage>&#x2013;<lpage>1310</lpage>. doi: <pub-id pub-id-type="doi">10.1093/icesjms/fsu243</pub-id>
</citation>
</ref>
<ref id="B65">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>V&#xf8;lstad</surname> <given-names>J. H.</given-names>
</name>
</person-group> (<year>2012</year>). <source>CIE independent peer review report</source>. (Bergen Norway:Center for Independent Experts) <uri xlink:href="https://www.st.nmfs.noaa.gov/science-quality-assurance/cie-peer-reviews/cie-review-2012">https://www.st.nmfs.noaa.gov/science-quality-assurance/cie-peer-reviews/cie-review-2012</uri>.</citation>
</ref>
<ref id="B66">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Wakabayashi</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Bakkala</surname> <given-names>R. G.</given-names>
</name>
<name>
<surname>Alton.</surname> <given-names>M. S.</given-names>
</name>
</person-group> (<year>1985</year>). &#x201c;<article-title>Methods of the U.S.-Japan demersal trawl surveys</article-title>,&#x201d; in <source>Results of cooperative U.S.-Japan groundfish investigations in the Bering Sea during May-August 1979</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Bakkala</surname> <given-names>R. G.</given-names>
</name>
<name>
<surname>Wakabayashi</surname> <given-names>K.</given-names>
</name>
</person-group> (<publisher-name>Vancouver, CA: Int. North Pac. Fish. Comm. Bull. 44</publisher-name>), <fpage>7</fpage>&#x2013;<lpage>29</lpage>.</citation>
</ref>
<ref id="B67">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wilberg</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Thorson</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Linton</surname> <given-names>B. C.</given-names>
</name>
<name>
<surname>Berkson</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Incorporating time-varying catchability into population dynamic stock assessment models</article-title>. <source>Rev. Fisheries Sci.</source> <volume>18</volume> (<issue>1</issue>), <fpage>7</fpage>&#x2013;<lpage>24</lpage>. doi: <pub-id pub-id-type="doi">10.1080/10641260903294647</pub-id>
</citation>
</ref>
<ref id="B68">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wyllie-Echeverria</surname> <given-names>T. I. N. A.</given-names>
</name>
<name>
<surname>Wooster</surname> <given-names>W. S.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>Year-to-year variations in Bering Sea ice cover and some consequences for fish distributions</article-title>. <source>Fisheries Oceanography</source> <volume>7</volume> (<issue>2</issue>), <fpage>159</fpage>&#x2013;<lpage>170</lpage>. doi: <pub-id pub-id-type="doi">10.1046/j.1365-2419.1998.00058.x</pub-id>
</citation>
</ref>
<ref id="B69">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Zacher</surname> <given-names>L. S.</given-names>
</name>
<name>
<surname>Richar</surname> <given-names>J. I.</given-names>
</name>
<name>
<surname>Foy</surname> <given-names>R. J.</given-names>
</name>
</person-group> (<year>2020</year>). <source>The 2019 eastern and northern Bering Sea continental shelf trawl surveys: Results for commercial crab species</source> (<publisher-name>Anchorage, AK: U.S. Department of Commerce NOAA Technical Memorandum NMFS-AFSC-400</publisher-name>). 234 p.</citation>
</ref>
<ref id="B70">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zimmermann</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Enberg</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Can less be more? Effects of reduced frequency of surveys and stock assessments</article-title>. <source>ICES J. Mar. Sci.</source> <volume>74</volume> (<issue>1</issue>), <fpage>56</fpage>&#x2013;<lpage>68</lpage>. doi: <pub-id pub-id-type="doi">10.1093/icesjms/fsw134</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>