By Emmanis Dorval Jenny McDaniel Southwest Fisheries Science Center Dianna Porzio

1. Age and Growth of Pacific Sardine in California During a Period of Stock Recovery and Geographical Expansion By Emmanis Dorval Jenny McDaniel Southwest Fisheries Science Center Dianna Porzio California Department of Fish and Game

2. Pacific SardineSardinop sagax

3. Background Patterns of individual fish growth rate • After population decline: • Following low population level due to exploitation (or to other events), individual fish growth rate is expected to increase (Le Cren et al. 1972, Botsford 1981). • Following low population level, individual fish growth rate may be several times higher than normal growth rate (i.e., pre-decline or pre-exploitation level ), Le Cren et al. (1972).

4. Background Patterns of individual fish growth rate • After population rebuilding: • Individual fish growth rate is generally expected to decline toward normal growth rate (i.e., pre-decline or pre-exploitation level ). • However, the extent of this decline may depend on the strength of regulating factors that are in play; and/or on how much the environment/habitat of this species been modified compared to pre-decline or pre-exploitation level.

5. Background Patterns of individual sardine growth rate • In the historical fishery (1937-1957): • Length-at-first annulus formation: 101 – 131 mm (Marr 1960) • During the recovery of the stock (1990s – 2000s) • Butler et al. (1996) • Hill et al. (2007, 2009, 2011)

6. Background Patterns of individual sardine growth rate P. sardine growth curve Butler et (1996) Fish samples Daily Egg Production Method (DEPM): April - May 1994 Ensenada fishery landings Von Bertallanfy growth model L∞= 205.4 ± 1.6 k= 1.19 ± 0.04 t0= 0 (fixed) Birthdate: January 1

7. Background Patterns of individual sardine growth rate P. Sardine growth curve Hill et al. (2011) Fish samples Daily Egg Production Method (DEPM): April – May: 2004-2010 Ensenada/California fishery landings: 1993 - 2010 Pacific Northwest fishery landings: 1999-2009 Birthdate: July 1 • L0.5yr= 114 mm • L∞= 242 mm • k= 0.364

8. Background • Sardine growth rate: current issues • Density -dependent changes in growth rates or reproduction have not been identified nor evaluated (Hill et al. 2011). • Combining fishery-dependent and fishery-independent data may not be the most reliable way to detect density-dependent. • A single growth model is derived in stock assessment model and applied to all fisheries (i.e., Ensenada and/or California, Pacific Northwest). • There are no known aged-fish: it is difficult to estimate age-reading bias and imprecision from multiple readings of otoliths.

9. Research Objectives • Develop methods to estimate and compare growth rates in the California spawning stock: • During the recovery of the stock in California (1980s-1990s). • During the recovery and expansion of the stock from California into the Pacific Northwest and British Columbia (1990s – present). • Determine how fish migration and fishing development may affect sardine growth rate in the California spawning stock.

10. Methods • Sample collection • DEPM April survey: 2004-2010 • California: San Diego – San Francisco • Age estimation • Annual increment from whole otolith, i.e. unpolished otolith, Yaremko (1996) • Birthdate assumption: July 1 • Age-readers: 3

11. Methods

12. Methods • Growth model • Von Bertalanffy growth model (Beverton and Holt, 1954) • Maximum likelihood estimation of parameters • Mixed-effect model • Lt~ Age | Reader • L∞ (random factor ) • k(fixed factor) • t0(random factor)

13. Methods Age-1+ sardine biomass Growth comparison between: Year-classes: 1994-2003 Year-classes: 2004-2008

14. Preliminary results Butler et al. 1996 Least-squares based model (Reader best ages) L∞= 205.4 mm k= 1.19 t0= 0 (fixed)

15. Preliminary results Least-squares based model (Reader best ages) Mixed -effect model Length ~Age | Reader L∞= 328 mm k= 0.14 t0= -3.13 L∞=283 mm ± 9 k=0.22 ± 0.01 t0= -3.54 ± 0.5

16. Preliminary results Mixed effect model Length ~Age | Reader Mixed -effect model Length ~ Age | Reader L∞=294 mm ± 10 k=0.19 ± 0.01 t0= -4.66 ± 1.31 L∞=263 mm ± 10 k=0.26 ± 0.01 t0= -3.11± 0.68

17. Summary • Compared to previous study, i.e. Butler et al. (1996), sardine growth rate in the spawning stock appears to have decreased. However, early age and growth data may need to be re-analyzed so that comparison can be based on most recent age-reading and statistical methods. • Growth rate of sardine in the 1994-2003 and in the 2004-2008 year-classes are significantly different. However, comparison of these year-classes may need to be based on data that include similar age-classes.

18. Summary • We will continue to age our historical sardine otolith samples to address these concerns. • We will include age-reading errors in future models, so that all model assumptions can be met.

19. Acknowledgement This study is an ongoing research collaboration project between the California Department of Fish and Game (CDFG) and the Southwest Fisheries Science Center (SWFSC). We thank all CDFG staff that participated in the age-reading of sardine otoliths, specially Valerie Taylor. We tank all SWFSC staff that participated in the collection and management of sardine otolith samples, specially Beverly Macewicz.