-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathESP_Pacific_cod.Rmd
127 lines (106 loc) · 10.5 KB
/
ESP_Pacific_cod.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
---
title: ESP Pacific cod morphometric condition
author:
- affiliation: RACE
description: Research Fisheries Biologist
email: Sean.Rohan@NOAA.gov
name: Sean Rohan
output: word_document
fontsize: 12pt
bibliography: EBS_references.bib
csl: fish-and-fisheries.csl
addr:
l1: 7600 Sand Point Way NE
l2: NMFS RACE Division, Groundfish Assessment Program
l3: Seattle, WA 98115
---
# Eastern Bering Sea Pacific cod
## Adults (≥46 cm FL)
*Status and Trends:* In 2022, the condition of adult Pacific cod in the EBS was neutral (within one standard deviation of the time series mean), which continues the trend of neutral morphometric condition since 2018. The neutral condition in recent years (2018–2022) represents an increase from the three prior years with negative to neutral condition from 2015–2017. This year, morphometric condition indicator was estimated using VAST [@Thorson2019; @Gruss2020a] instead of the stratum biomass weighted length-weight residual method from previous years. Trends in the historical and new indicator are similar based on the strong correlation between the historical and new indicator (r = 0.87), although there are notable differences for some years. The year with the lowest condition for the old indicator was 1999 (a ‘cold’ year with an early survey start), with an anomaly greater than two standard deviations from the mean. Based on the new VAST relative condition indicator, 1999 was a neutral year and the year with the lowest condition was 2012. Despite these differences, new indicator trends generally match the trend from the old indicator with condition increasing from 1999 to 2003, decreasing from 2003 to 2006, then fluctuating around neutral from 2007 to 2022, aside from negative (>1 standard deviation below the mean) years in 2009, 2012, 2015, and 2017.
## Juveniles (<46 cm FL)
*Status and Trends:* In 2022, the condition of juvenile Pacific cod in the EBS was neutral (within one standard deviation of the time series mean), which continues the trend of neutral morphometric condition since 2017. This year, the morphometric condition indicator was estimated using VAST [@Thorson2019; @Gruss2020a] instead of the stratum biomass weighted length-weight residual method from previous years. Trends in the historical and new indicator are similar based on the strong correlation between the historical and new indicator (r = 0.91), although there are noteable differences in some years. Specifically, 2017 was a negative year using the old indicator (~1 standard deviation below the mean) but a neutral year with the new indicator, negative years in 2009 and 2012 are still negative but the anomaly is larger, and the anomaly in 1999 decreased from 3.2 standard deviations below them mean to 1.8 standard deviations below the mean (a ‘cold’ year with an early survey start). Despite these differences, new indicator trends generally match the trend from the old indicator with condition increasing from 1999 to 2004, decreasing from 2005 to 2010, then fluctuating around neutral from 2011 to 2022, aside from a negative year in 2012 and positive year in 2016.
## Influential Factors
*Influential Factors:* Many factors contribute to variation in morphometric condition so it is unclear which specific factors contributed to neutral condition of adult Pacific cod in the EBS in 2022. Factors that may contribute to variation in morphometric condition include environmental conditions that affect prey quality and temperature-dependent metabolic rates, survey timing, stomach fullness of individual fish, fish migration patterns, and the distribution of samples within survey strata. Temperature is an important factor that can influence the morphometric condition of Pacific cod by influencing metabolic rates, prey availability, and prey quality. Historically in the eastern Bering Sea (EBS), ‘cold’ years (with a small cold pool) were associated with negative morphometric condition (e.g., 1999, 2012) and warm years (e.g., 2002-2005) were associated with positive morphometric condition. However, during recent (2018–2021) exceptionally warm recent years, the morphometric condition of Pacific cod was neutral for adult and juvenile Pacific cod and this trend continued into the average temperature year in 2022. Temperature can negatively affect growth rates if prey resources are insufficient to make up for increased metabolic demand. In GOA, elevated temperatures during the 2014–2016 marine heatwave in the Gulf of Alaska were associated with lower growth rates of Pacific cod and lower morphometric condition in 2015 (adults and juveniles combined), likely because diminished prey resources during the heatwave were insufficient to make up for increased metabolic demand [@Barbeaux2020]. Additional information about the groundfish morphometric condition indicator and factors that can influence estimates of morphometric condition are described in the EBS Groundfish Morphometric Condition contribution in the 2022 Eastern Bering Sea Ecosystem Status Report (Rohan et al., _In prep_).
```{r steup_data, include=FALSE}
library(akfishcondition)
akfishcondition:::PCOD_ESP$FULL_REGION_EBS %>%
dplyr::filter(common_name == "Pacific cod (juvenile)") %>%
dplyr::mutate(indicator_name = "Summer_Pacific_Cod_Condition_Juvenile_EBS_Survey",
data_value = vast_relative_condition) %>%
dplyr::select(year, indicator_name, data_value) %>%
write.csv(file = here::here("output", "ESP_Summer_Pacific_Cod_Condition_Juvenile_EBS_Survey.csv"),
row.names = FALSE)
akfishcondition:::PCOD_ESP$FULL_REGION_EBS %>%
dplyr::filter(common_name == "Pacific cod (adult)") %>%
dplyr::mutate(indicator_name = "Summer_Pacific_Cod_Condition_Adult_EBS_Survey",
data_value = vast_relative_condition) %>%
dplyr::select(year, indicator_name, data_value) %>%
write.csv(file = here::here("output", "ESP_Summer_Pacific_Cod_Condition_Adult_EBS_Survey.csv"),
row.names = FALSE)
print(paste0("Last Update: ", akfishcondition:::EBS_INDICATOR$LAST_UPDATE))
```
```{r ebs_anomaly, include=TRUE, echo=FALSE, fig.height=6,fig.width=12, fig.cap="\\label{fig:figs}Figure 1. VAST relative condition for adult and juvenile Pacific cod collected during AFSC/RACE GAP standard summer bottom trawl surveys of the eastern Bering Sea shelf, 1999 to 2022. The dash in the blue boxes denote the mean for that year, the blue box denotes one standard error, and the lines on the boxes denote two standard error. Horizontal lines on each plot represent the historical mean, dashed lines denote one standard deviation, and dotted lines denote two standard deviations.", message=FALSE, warning=FALSE}
fig1_rmd <- plot_anomaly_timeseries(x = akfishcondition::PCOD_ESP$FULL_REGION_EBS,
region = "BS",
fill_color = "#0085CA",
var_y_name = "vast_relative_condition",
var_y_se_name = "vast_relative_condition_se",
var_x_name = "year",
y_title = "VAST relative condition",
plot_type = "box",
set_intercept = 1,
format_for = "rmd")
fig1_png <- plot_anomaly_timeseries(x = akfishcondition::PCOD_ESP$FULL_REGION_EBS,
region = "BS",
fill_color = "#0085CA",
var_y_name = "vast_relative_condition",
var_y_se_name = "vast_relative_condition_se",
var_x_name = "year",
y_title = "VAST relative condition",
plot_type = "box",
set_intercept = 1,
format_for = "png")
print(fig1_rmd + theme_condition_index())
```
```{r ebs_anomaly_png, include=FALSE, echo=FALSE, message=FALSE, warning=FALSE}
png(here::here("plots", "BS", "ESP_EBS_PCOD_condition.png"),width=6,height=3,units="in",res=600)
print(fig1_png + theme_blue_strip())
dev.off()
```
```{r ebs_vast_sbw, include=TRUE, echo=FALSE, fig.height=6,fig.width=12,fig.cap="\\label{fig:figs}Figure 3. Time series of VAST relative condition and length-weight residual condition anomalies for adult and juvenile Pacific cod in the eastern Bering Sea. Triangles denote the length-weight residual, while points denote the VAST relative condition. Lines represent the historical mean, dashed lines denote one standard deviation, and dotted lines denote two standard deviations. The Pearson correlation coefficient (_r_) is shown at the bottom right of each panel.", message=FALSE, warning=FALSE}
fig2_rmd <- plot_two_timeseries(x_1 = akfishcondition::PCOD_ESP$FULL_REGION_EBS,
x_2 = akfishcondition::PCOD_ESP$FULL_REGION_EBS,
region = "BS",
series_name_1 = "Length-weight residual (Historical)",
series_name_2 = "VAST relative condition (New)",
var_y_name_1 = "mean_wt_resid",
var_x_name_1 = "year",
var_y_name_2 = "vast_relative_condition",
var_x_name_2 = "year",
y_title = "Condition anomaly (Z-score)",
scale_y = TRUE,
year_break = 2020,
fill_colors = c("#001743", "#0085CA"),
format_for = "rmd")
fig2_png <- plot_two_timeseries(x_1 = akfishcondition::PCOD_ESP$FULL_REGION_EBS,
x_2 = akfishcondition::PCOD_ESP$FULL_REGION_EBS,
region = "BS",
series_name_1 = "Length-weight residual (Historical)",
series_name_2 = "VAST relative condition (New)",
var_y_name_1 = "mean_wt_resid",
var_x_name_1 = "year",
var_y_name_2 = "vast_relative_condition",
var_x_name_2 = "year",
y_title = "Condition anomaly (Z-score)",
scale_y = TRUE,
year_break = 2020,
fill_colors = c("#001743", "#0085CA"),
format_for = "png")
print(fig2_rmd + theme_condition_index() + theme(legend.position = "bottom"))
```
```{r ebs_vast_sbw_png, include=FALSE, echo=FALSE, message=FALSE, warning=FALSE}
png(file = here::here("plots", "BS", "ESP_PCOD_VAST_vs_SBW_time_series.png"), width = 6, height = 3, units = "in", res = 600)
print(fig2_png + theme_blue_strip() + theme(legend.position = "bottom"))
dev.off()
```
# References