Bay - Delta Science Conference November 15 - 17, 2016 Flora Cordoleani, Noble Hendrix, Eric Danner and Steve Lindley
Bay-Delta Science ConferenceNovember 15-17, 2016
Flora Cordoleani, Noble Hendrix, Eric Danner and Steve Lindley
Historic vs current distribution of spring-run Chinook
Only 3 out of 19 historic independent populations of CV spring-run Chinook salmon are extant: Mill, Deer, and Butte creeks
Represent only the Northern Sierra Nevada diversity group
Listed as threatened under the federal Endangered Species Act (ESA) since 1999.
“The status of the CV spring-run Chinook salmon ESU has probably improved on balance sincethe 2010 status review, through 2014 […].”
“The recent declines of many of the dependent populations, high pre-spawn and egg mortality, and uncertain juvenile survival during the 2012 to 2015 drought, ocean conditions, as well as the level of straying of FRFH spring-run Chinook salmon to other CV spring-run Chinook salmon populations are all causes for concern for the long-term viability of the CV spring-run Chinook salmon ESU.” [Johnson and Lindley, SR viability report (2016) and NOAA-NMFS 5 year status review report (2016)]
Central Valley spring-run Chinook viability status
Objectives
1. Understand clearly the dynamics of Central Valley spring-run Chinook salmon in the freshwater and the ocean
2. Identify the environmental factors influencing changes in abundance of spring-run Chinook salmon populations
3. Predicting possible impacts of future water management and climate change scenarios on their dynamics
CV spring-run Chinook salmon life cycle
. ·-_ .... _ --=---- -.e; ...... --
' I <j eu'TM ~vrMa ~ ' •111\J" ~""'' t IV .... ull
• VOI.O I V ... Il
·~--
CV spring-run life cycle model
rua r
)-" {:~. 1£ hn-..
, .. 01 =
Js.c.IMr •1\0chll
22
All\ 1 1
76 parameters Period of model simulation: 1985 – 2008Temporal Resolution
• Annual for ocean stages• Monthly for freshwater stages
Spatial Resolution• Regional depiction of rearing habitat types into Tributaries, Sutter Bypass,
Sac. River, Yolo Bypass, Delta, and BayModel validation by fitting simulated adult abundance to historical adult
escapement abundance (Grand Tab)
Model Structure
MarAprMayJunJulAugSepOctNovDecJanFebMarAprMayJunJulAugSepOctNovDecJanFeb
MarAprMayJunJulAugSepOctNovDecJanFebMarAprMayJunJulAugSepOctNovDecJanFeb
MarAprMayJunJulAugSepOctNov
Adult enter tributaries
Spawn
Age 3
Juvenile emergence
YoY out migration
Yearling out migration
Age 2 exit ocean
Age 2
Spawning(Age 2)
Age 4
Age 3 exit ocean
Spawning(Age 3)
Age 4 exit ocean
Spawning(Age 4)
Age 1
Spring-run aging convention
Juvenile salmon
survival during rearing and
outmigration to the Ocean
Young of the Year vs Yearling
Different juvenile rearing/migration strategy for spring-run Chinook
1. Young of the Year that rear for several months and migrate in the spring
2. Yearling that stays an entire year in the natal reaches before migrating to the Ocean
Tidal Fry disperse instantaneously after emergence
TidalFry = PTF * Fry
PTF = Proportion of Tidal fry
Density independent migration of fry
Butte,Mill/Deer Cr. Fry
Tidal Fry Migrant Fry
Creek Fry
Sac River Fry
Floodplain Fry
Bay Fry
Delta Fry
Sutter Fry
Ni,t+1 = Si 1−m Ni,t1+ ⁄Si 1−m Ni,t Ki,t
and Mi,t = Si * Ni,t – Ni,t+1
Ni,t+1 = resident fry abundanceMi,t = migrant fry abundanceSi = fry survivalm = migration rate without density dependenceKi,t = rearing capacity of habitat i
Density dependent migration of fry
Butte,Mill/Deer Cr. Fry
Tidal Fry Migrant Fry
Habitat type Variable Habitat quality Variable rangeMainstem Velocity High <= 0.15 m/s
Low > 0.15 m/sDepth High > 0.2 m, <= 1 m
Low <= 0.2 m, > 1 mDelta Channel type High Blind channels
Low Mainstem, distributaries, open water
Depth High > 0.2 m, <= 1.5 mLow <= 0.2 m, > 1.5 m
Cover High VegetatedLow Not vegetated
Bay Shoreline type High Beaches, marshes, vegetated banks, tidal flats
Low Riprap, structures, rocky shores, exposed habitats
Depth High > 0.2 m, <= 1.5 mLow <= 0.2 m, > 1.5 m
Salinity High <= 10 pptLow > 10 ppt
Rearing Capacity estimate (C. Greene, NWFSC)
Survival of rearing fry in the Delta
Use Newman (2003) survival rate relationship:
logit(SDelta,t)= Brearing * XRearing,i,t
Xrearing = Flow, Temperature, Exports, DCC
Survival of sm0lt migrating to the Ocean
Survival rate in the Sutter Bypass based on acoustic tagging study:
logit(St)= B0 + B1 * Flow
Survival rate through the Delta from:1. ePTM simulations [Sridharan, V., and Byrne, B.]2. Newman equations
200 400 600 800 1000
0.2
0.4
0.6
0.8
Flow experienced
Surv
ival
Early ocean survival
Early ocean survival of smolts depends on ocean conditions in the Gulf of Farallones and the fish rearing origin
logit(Si) = (B0,i + B0-add,i )+ (B1,I + B1-add,i) * OPI
OPI = Ocean productivity indexB0-add,I = poor habitat interceptB1-add,I = poor habitat slope
Survival of adult during holding period
Significant adult pre-spawning mortality events in 2002 and 2003 have been reported for Butte Creek population
Summer pre-spawing survival expressed as a function of water temperature [Thompson et al. (2012)] :
Sps,t = 11+𝑒𝑒−𝑏𝑏1−𝑏𝑏2𝑇𝑇
T = Temperature in holding habitat
Mill/Deer Cr. Model sensitivity analysis
ePTMNewman
0.0
0.2
0.4
0.6
0.8
1.0
Sens
itiv_
avg.
MD
C
logi
t.s.e
ggs.
MD
Clo
git.p
.tf.M
DC
logi
t.s.tf
.MD
CS
aclo
git.p
.tf.M
DC
Sac
logi
t.s.tf
.MD
CS
utte
rlo
git.s
.tf.M
DC
fplo
git.s
.tf.M
DC
delta
logi
t.s.re
ar.M
DC
trib
logi
t.s.re
ar.M
DC
sutte
rlo
git.s
.rear
.MD
Criv
erlo
git.s
.rear
.MD
Cyo
lolo
git.s
.rear
.MD
Cba
ylo
git.m
logi
t.s.m
fblo
git.p
.yea
rl.M
DC
logi
t.s.y
earl.
MD
Clo
git.p
.yea
rlmig
r.MD
Clo
git.s
.15.
MD
Clo
git.s
.16a
.MD
Clo
git.s
.sm
olt.b
aylo
git.s
.17a
logi
t.s.1
8alo
git.s
.19a
logi
t.s.2
1b0
.22_
27b1
.22_
27b0
.22_
27.a
ddb1
.22_
27.a
ddlo
git.s
.28
logi
t.s.2
9lo
git.m
.2lo
git.s
.spa
wng
logi
t.s.s
paw
ng.y
earl
logi
t.s.3
1lo
git.s
.32
logi
t.m.3
logi
t.s.3
4lo
git.s
.35
Age
2.fe
mra
teA
ge3.
fem
rate
Age
4.fe
mra
te
SI = dist(N_adult_+/−10%,N_adult_fixed)/N_adult_fixeddist(Par_+/−10%, Par_fixed) /Par_fixed
01
23
45
Sens
itiv_
avg.
MD
C
logi
t.s.e
ggs.
MD
Clo
git.p
.tf.M
DC
logi
t.s.tf
.MD
CS
aclo
git.p
.tf.M
DC
Sac
logi
t.s.tf
.MD
CS
utte
rlo
git.s
.tf.M
DC
fplo
git.s
.tf.M
DC
delta
logi
t.s.re
ar.M
DC
trib
logi
t.s.re
ar.M
DC
sutte
rlo
git.s
.rear
.MD
Criv
erlo
git.s
.rear
.MD
Cyo
lolo
git.s
.rear
.MD
Cba
ylo
git.m
logi
t.s.m
fblo
git.p
.yea
rl.M
DC
logi
t.s.y
earl.
MD
Clo
git.p
.yea
rlmig
r.MD
Clo
git.s
.15.
MD
Clo
git.s
.16a
.MD
Clo
git.s
.sm
olt.b
aylo
git.s
.17a
logi
t.s.1
8alo
git.s
.19a
logi
t.s.2
1b0
.22_
27b1
.22_
27b0
.22_
27.a
ddb1
.22_
27.a
ddlo
git.s
.28
logi
t.s.2
9lo
git.m
.2lo
git.s
.spa
wng
logi
t.s.s
paw
ng.y
earl
logi
t.s.3
1lo
git.s
.32
logi
t.m.3
logi
t.s.3
4lo
git.s
.35
Age
2.fe
mra
teA
ge3.
fem
rate
Age
4.fe
mra
te
Mill/Deer Cr. model sensitivity analysis
d I
'""'
- -::-----::::---,1 10 .' ~,_ ,._ I 9 l, 4 ~L fiM!or SHy . ,._, ~ ~
' J SUIIIIf 11 Sutltr Ovp•u GH) J-'1 ovpus RFry Smolt I Sutter Bypa s
In Oceah
® r--- --t
ePTMNewman
0.0
0.5
1.0
1.5
2.0
2.5
3.0
Sens
itiv_
avg.
BC
logi
t.s.e
ggs.
BC
logi
t.p.tf
.BC
logi
t.s.tf
.BC
Sut
ter
logi
t.p.tf
.BC
Sut
ter
logi
t.s.tf
.BC
fplo
git.s
.tf.B
Cde
ltalo
git.s
.rear
.BC
trib
logi
t.s.re
ar.B
Csu
tter
logi
t.s.re
ar.B
Cyo
lolo
git.s
.rear
.BC
bay
logi
t.mlo
git.s
.mfb
logi
t.p.y
earl.
BC
logi
t.s.y
earl.
BC
logi
t.p.y
earlm
igr.B
Clo
git.s
.15.
BC
logi
t.s.1
6a.B
Clo
git.s
.sm
olt.b
aylo
git.s
.18a
logi
t.s.1
9alo
git.s
.21
b0.2
2_27
b1.2
2_27
b0.2
2_27
.add
b1.2
2_27
.add
logi
t.s.2
8lo
git.s
.29
logi
t.m.2
logi
t.s.s
paw
nglo
git.s
.spa
wng
.yea
rllo
git.s
.31
logi
t.s.3
2lo
git.m
.3lo
git.s
.34
logi
t.s.3
5A
ge2.
fem
rate
Age
3.fe
mra
teA
ge4.
fem
rate
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
Sens
itiv_
avg.
BC
logi
t.s.e
ggs.
BC
logi
t.p.tf
.BC
logi
t.s.tf
.BC
Sut
ter
logi
t.p.tf
.BC
Sut
ter
logi
t.s.tf
.BC
fplo
git.s
.tf.B
Cde
ltalo
git.s
.rear
.BC
trib
logi
t.s.re
ar.B
Csu
tter
logi
t.s.re
ar.B
Cyo
lolo
git.s
.rear
.BC
bay
logi
t.mlo
git.s
.mfb
logi
t.p.y
earl.
BC
logi
t.s.y
earl.
BC
logi
t.p.y
earlm
igr.B
Clo
git.s
.15.
BC
logi
t.s.1
6a.B
Clo
git.s
.sm
olt.b
aylo
git.s
.18a
logi
t.s.1
9alo
git.s
.21
b0.2
2_27
b1.2
2_27
b0.2
2_27
.add
b1.2
2_27
.add
logi
t.s.2
8lo
git.s
.29
logi
t.m.2
logi
t.s.s
paw
nglo
git.s
.spa
wng
.yea
rllo
git.s
.31
logi
t.s.3
2lo
git.m
.3lo
git.s
.34
logi
t.s.3
5A
ge2.
fem
rate
Age
3.fe
mra
teA
ge4.
fem
rate
Butte Cr. Model sensitivity analysis
Butte cr. model sensitivity analysis
I I r~~~ :
~·,··-"-- --~--{AfeJS,:..•vn '.. ... -
---... 39
-------- AJ~A,pnwM" ~
Butte Cr
Time
Estim
ated
# o
f Spa
wne
rs
1985 1990 1995 2000 2005
050
0015
000
2500
0 Mill/Deer Cr.
TimeEs
timat
ed #
of S
paw
ners
1985 1990 1995 2000 2005
010
0020
0030
0040
0050
00
Model simulations
Next Steps
Refine parameter values and finish model calibration- Rearing capacity in Tributaries and Sutter Bypass- Proportion of tidal fry vs rearing fry vs yearling - Evaluate relationship between egg survival and temperature in spawning
habitat
Use model for inference in evaluating water management and climate change scenarios
- Effect of increased temperature in spawning habitat?- Sutter Bypass flooding scenarios - Delta water management scenarios
0.00
0.10
0.20
0.30
obsv
spre
d
logi
t.s.e
ggs.
MD
Clo
git.p
.tf.M
DC
logi
t.s.tf
.MD
CSa
clo
git.p
.tf.M
DC
Sac
logi
t.s.tf
.MD
CSu
tter
logi
t.s.tf
.MD
Cfp
logi
t.s.tf
.MD
Cde
ltalo
git.s
.rear
.MD
Ctri
blo
git.s
.rear
.MD
Csu
tter
logi
t.s.re
ar.M
DC
river
logi
t.s.re
ar.M
DC
yolo
logi
t.s.re
ar.M
DC
bay
logi
t.mlo
git.s
.mfb
logi
t.p.y
earl.
MD
Clo
git.s
.yea
rl.M
DC
logi
t.p.y
earlm
igr.M
DC
logi
t.s.1
5.M
DC
logi
t.s.s
mol
t.bay
logi
t.s.1
9alo
git.s
.21
b0.2
2_27
b1.2
2_27
b0.2
2_27
.add
b1.2
2_27
.add
logi
t.s.2
8lo
git.s
.29
logi
t.m.2
logi
t.s.s
paw
nglo
git.s
.spa
wng
.yea
rllo
git.s
.31
logi
t.s.3
2lo
git.m
.3lo
git.s
.34
logi
t.s.3
5Ag
e2.fe
mra
teAg
e3.fe
mra
teAg
e4.fe
mra
te
0.00
0.10
0.20
0.30
obsv
spre
d
logi
t.s.e
ggs.
BClo
git.p
.tf.B
Clo
git.s
.tf.B
CSu
tter
logi
t.p.tf
.BC
Sutte
rlo
git.s
.tf.B
Cfp
logi
t.s.tf
.BC
delta
logi
t.s.re
ar.B
Ctri
blo
git.s
.rear
.BC
sutte
rlo
git.s
.rear
.BC
yolo
logi
t.s.re
ar.B
Cba
ylo
git.m
logi
t.s.m
fblo
git.p
.yea
rl.BC
logi
t.s.y
earl.
BClo
git.p
.yea
rlmig
r.BC
logi
t.s.1
5.BC
logi
t.s.s
mol
t.bay
logi
t.s.1
9alo
git.s
.21
b0.2
2_27
b1.2
2_27
b0.2
2_27
.add
b1.2
2_27
.add
logi
t.s.2
8lo
git.s
.29
logi
t.m.2
logi
t.s.s
paw
nglo
git.s
.spa
wng
.yea
rllo
git.s
.31
logi
t.s.3
2lo
git.m
.3lo
git.s
.34
logi
t.s.3
5Ag
e2.fe
mra
teAg
e3.fe
mra
teAg
e4.fe
mra
te
0.0
0.2
0.4
0.6
0.8
obsv
spre
d
logi
t.s.e
ggs.
BClo
git.p
.tf.B
Clo
git.s
.tf.B
CSu
tter
logi
t.p.tf
.BC
Sutte
rlo
git.s
.tf.B
Cfp
logi
t.s.tf
.BC
delta
logi
t.s.re
ar.B
Ctri
blo
git.s
.rear
.BC
sutte
rlo
git.s
.rear
.BC
yolo
logi
t.s.re
ar.B
Cba
ylo
git.m
logi
t.s.m
fblo
git.p
.yea
rl.BC
logi
t.s.y
earl.
BClo
git.p
.yea
rlmig
r.BC
logi
t.s.1
5.BC
logi
t.s.s
mol
t.bay
logi
t.s.1
9alo
git.s
.21
b0.2
2_27
b1.2
2_27
b0.2
2_27
.add
b1.2
2_27
.add
logi
t.s.2
8lo
git.s
.29
logi
t.m.2
logi
t.s.s
paw
nglo
git.s
.spa
wng
.yea
rllo
git.s
.31
logi
t.s.3
2lo
git.m
.3lo
git.s
.34
logi
t.s.3
5Ag
e2.fe
mra
teAg
e3.fe
mra
teAg
e4.fe
mra
te
0.0
0.2
0.4
0.6
obsv
spre
d
logi
t.s.e
ggs.
MD
Clo
git.p
.tf.M
DC
logi
t.s.tf
.MD
CSa
clo
git.p
.tf.M
DC
Sac
logi
t.s.tf
.MD
CSu
tter
logi
t.s.tf
.MD
Cfp
logi
t.s.tf
.MD
Cde
ltalo
git.s
.rear
.MD
Ctri
blo
git.s
.rear
.MD
Csu
tter
logi
t.s.re
ar.M
DC
river
logi
t.s.re
ar.M
DC
yolo
logi
t.s.re
ar.M
DC
bay
logi
t.mlo
git.s
.mfb
logi
t.p.y
earl.
MD
Clo
git.s
.yea
rl.M
DC
logi
t.p.y
earlm
igr.M
DC
logi
t.s.1
5.M
DC
logi
t.s.s
mol
t.bay
logi
t.s.1
9alo
git.s
.21
b0.2
2_27
b1.2
2_27
b0.2
2_27
.add
b1.2
2_27
.add
logi
t.s.2
8lo
git.s
.29
logi
t.m.2
logi
t.s.s
paw
nglo
git.s
.spa
wng
.yea
rllo
git.s
.31
logi
t.s.3
2lo
git.m
.3lo
git.s
.34
logi
t.s.3
5Ag
e2.fe
mra
teAg
e3.fe
mra
teAg
e4.fe
mra
te