Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
ICES Journal of Marine Science (2011), 68(6), 1131–1137. doi:10.1093/icesjms/fsr033 Regional differences in climate factors controlling chum and pink salmon abundance Masa-aki Fukuwaka 1*, Toshiki Kaga 2‡, and Tomonori Azumaya1 1 Hokkaido National Fisheries Research Institute (HNFRI), Fisheries Research Agency (FRA), 116 Katsurakoi, Kushiro 085-0802, Japan National Salmon Resources Center, FRA, 2-2 Nakanoshima, Toyohira-ku, Sapporo 062-0922, Japan 2 *Corresponding Author: tel: +81 154 92 1715; fax: +81 154 91 9355; e-mail: [email protected]. Present address: HNFRI, FRA, 116 Katsurakoi, Kushiro 085-0802, Japan. ‡ Fukuwaka, M., Kaga, T., and Azumaya, T. 2011. Regional differences in climate factors controlling chum and pink salmon abundance. – ICES Journal of Marine Science, 68: 1131 – 1137. Received 17 June 2010; accepted 20 February 2011; advance access publication 11 April 2011. Chum and pink salmon abundances vary on a decadal time-scale. We examined the relationship between large-scale climate indices (CIs), regional climate factors (RFs), and rates of change in regional catches (RCs) of chum and pink salmon in five regions of the North Pacific. Correlation coefficients of RCs with RFs were larger than those of RCs with CIs, although the correlation coefficient of particular variables varied among regions. Climate affected salmon stocks as indicated by significant relationships with various terrestrial and ocean climate factors on a regional scale. These results suggest that no single CI or RF controls salmon abundance in all regions; however, global climate changes could affect regional climate directly and regional salmon abundance indirectly. A warming trend in the North Pacific might affect the long-term change in salmon abundance. The mechanisms controlling regional salmon abundance must be understood better to forecast successfully future conditions for Pacific salmon stocks, because the response of salmon stocks to global climate change varies among regions. Keywords: chum salmon, climate change, large-scale climate indices, pink salmon, regional process, regional stock dynamics. Introduction Pacific salmon (Oncorhynchus spp.) abundance, as indexed by commercial catch, has changed on a decadal time-scale with peaks in the 1930–1940s and 1990–2000s (Irvine and Fukuwaka, 2011). Decadal changes in salmon catch followed large-scale climate changes that were indicated by ocean climate indices (CIs; Beamish and Bouillon, 1993; Mantua et al., 1997). However, salmon catch can vary among regions. For example, along the west coast of North America, salmon catches in southern regions varied inversely with catches in northern regions (Hare et al., 1999). Catch trends for Asian pink salmon (O. gorbuscha) have been the opposite of those for Alaskan stocks (Nagasawa, 2000). Salmon abundance appears to be affected by large-scale climate change, although regional variation exists in the ways that salmon abundances respond. Spatial patterns of recruitment covariation for North American pink, chum (O. keta), and sockeye salmon (O. nerka), indicated that regional processes during freshwater and early marine stages might be more important for controlling salmon abundance than ocean basin-scale processes (Peterman et al., 1998; Pyper et al., 2001, 2002). Salmon survival and/or production have been correlated with winter air temperature (AT), discharge and flooding, sea surface temperature (SST) and salinity, and upwelling in coastal environments (Wickett, 1958; Levanidov, 1964; Blackbourn, 1990; Holtby et al., 1990; Mueter et al., 2002a; Morita et al., 2006). This is consistent with the concept that early mortality and growth are important in determining the # 2011 abundance of individual year classes (Parker, 1962, 1968; Bradford, 1995; Beamish and Mahnken, 2001). Chum and pink salmon are the two most abundant Pacific salmon species, whose catches account for 67– 91% of the total Pacific salmon catch (Irvine and Fukuwaka, 2011). These species have a short freshwater life, migrate to the sea soon after their emergence from gravel, and have similar early ocean lives in estuaries and subsequent ocean migration (Quinn, 2005). Therefore, their population dynamics might be affected similarly by climate. The objectives of this study were to clarify how climate change affects chum and pink salmon differently in various regions, and to determine whether recent climate change favours chum and pink salmon. We examined the relationship between rates of generational change at regional scales of chum and pink salmon abundances, large-scale CIs, and regional climate factors (RFs) around the North Pacific. Material and methods We analysed the relationship between the rates of change in regional catches (RCs) of chum and pink salmon, large-scale CIs, and RFs. Because salmon catch data were available by nation or region, we set five regions around the North Pacific: Japan, Russia, the US Northwest and British Columbia (NWBC), southeastern Alaska (SEAK), and western and central Alaska (WCAK). The study period included the 1925–2001 brood years of chum and pink salmon for which complete sets of catch and climate data were available. International Council for the Exploration of the Sea. Published by Oxford Journals. All rights reserved. For Permissions, please email: [email protected] 1132 M. Fukuwaka et al. Salmon catch and climate data A time-series of Pacific salmon catches since 1925 was compiled by the Working Group on Stock Assessment, North Pacific Anadromous Fish Commission (Eggers et al., 2003; Irvine et al., 2009). We used RCs for chum and pink salmon, except for high seas catches, because those catches included fish from mixed stocks. Catch is not the same as abundance, but catch trends can be used as indices of abundance if changes in fishing effort are small (Morita et al., 2006; Irvine and Fukuwaka, 2011). For RCs, we used the rate of change in catch per salmon generation, which should index population growth rates, assuming that fishing effort did not change between generations: RC = ln Ct+T , Ct where Ct is the salmon catch in numbers at year t, and T is a generation time, i.e. four for chum and two for pink salmon. For CIs, we used the Aleutian Low Pressure (ALPI), Pacific Decadal Oscillation (PDO) and Victoria (VI) indices, which were closely related to the climate and ecosystems of the North Pacific Ocean (Beamish et al., 1997; Mantua et al., 1997; Bond et al., 2003). We obtained the annual ALPI from the Pacific Biological Station, Fisheries and Oceans Canada (http://www. pac.dfo-mpo.gc.ca/science/species-especes/climatology-ie/ cori-irco/alpi/index-eng.html). For monthly PDO and VI, we used the first and second principal components of SST in the North Pacific recalculated by S. McKinnell (pers. comm.). Recalculated PDO correlated highly with the original PDO (r ¼ 0.811, p , 0.001; http://jisao. washington.edu/pdo/). The PDO and VI indices were averaged over the four seasons, i.e. winter, January –March; spring, April –June; summer, July –September; and autumn, October– December. For RFs, we used precipitation, AT on the land surface and SST. Monthly precipitation and AT data on a 0.58 latitude by 0.58 longitude grid until 2002 were available from the Climatic Research Unit, University of East Anglia (Mitchell and Jones, 2005; CRU TS 2.1 available at http://www.cru.uea.ac.uk/). Monthly SST data on a 18 latitude by 18 longitude grid were obtained from the Kobe Collection, Japanese Meteorological Agency. These data were averaged annually for the four seasons (i.e. winter, January –March; spring, April –June; summer, July –September; and autumn, October–December) and over the same regions for RCs, which were in the ranges 38 –468N 138–1468E for Japan; 46 –588N 135–1658E and 58 –638N 140–1808E for Russia; 45 – 558N 120– 1338W for NWBC; 55– 618N 130–1448W for SEAK; and 55 –628N 144–1658W for WCAK (Figure 1). Statistical analyses To examine relationships among RCs, CIs, and RFs, we used Pearson’s correlation coefficient with t-tests after adjusting the degree of freedom as appropriate, considering temporal autocorrelation between variables (Pyper and Peterman, 1998). We analysed the relationships among RCs for brood year t, CIs and RFs in seasons during freshwater and early ocean life, i.e. autumn at year t and winter, spring, and summer at year t + 1. To interpret relationships among variables, we used structural equation modelling (SEM) with unobserved latent variables (recently reviewed by Grace et al., 2010). SEM is the process of developing and evaluating structural equation models that can test theoretical models consisting of multiple cause-and-effect relationships with an observed dataset and evaluate causal effects through the study of direct and indirect path relations. SEM includes hypothesized linear relationships based on a theoretical understanding of the system of interest. Different model configurations reflecting alternative hypotheses are fitted to the data and evaluated. Path coefficients of the model demonstrate the strength of the estimated effect of a single variable on the response. Overall, model fit can be tested using the x 2-test. Computations were performed using the sem package for the statistical package R (version 2.9.2; Fox et al., 2009; R Development Core Team, 2009). To test whether large-scale ocean climate affects RCs in broad areas of the North Pacific, we analysed the relationships between CIs and RCs in all five regions using SEM. In this analysis, we included latent variables for decadal climate changes and for broad-area RCs, because of the correlation between CIs and RCs. The initial SEM model included all paths from latent decadal climate changes to observed CIs, latent decadal climate changes to latent broad-area RCs, observed CIs to latent broadarea RCs, and latent decadal climate changes to observed RCs. To select the best model, we used the backward elimination of these paths based on the Bayesian information criterion (BIC). By region and by species, to evaluate the indirect and direct effects of large-scale climate changes, we conducted an SEM Figure 1. Map displaying the study regions around the North Pacific. Rectangles indicate the range of a region for collecting data on climate factors. 1133 Climate controls of pink and chum salmon with variables for each RC, CIs, latent decadal climate changes, and RFs. RFs used in this analysis were chosen based on their significance in correlation analysis with the targeted RC. The starting model of the SEM process included all paths from CIs, latent decadal climate changes to RFs, and an RC. To select the best model, we again used the backward elimination of paths based on BIC. To evaluate the indirect effect of large-scale climate through regional climate to regional salmon dynamics, we used a compound path coefficient, which is a product of path coefficients along an indirect linkage from an independent variable to a dependent variable. A compound path coefficient indicates the strength of effect of an independent variable on a dependent variable through an indirect compound path, and the sign of the coefficient indicates the direction of the effect (i.e. positive or negative). In the analyses, because the correlation analysis included a total of 21 tests between an RC and separated independent variables (five variables by four seasons and ALPI), this multiple correlation test increased the possibility of type I error, where at least one independent variable would be incorrectly detected as being correlated significantly. An SEM for an RC could include all independent variables in a test simultaneously and avoid such a multiple-test problem. However, CIs and RFs might correlate with each other and potentially be confounding. For help in interpreting the results of SEMs considering multicollinearity, correlation matrices of CIs and RFs for SEMs are presented in the Supplementary material. Results Commercial catches of chum and pink salmon in the North Pacific peaked in the 1930–1940s and in the 1990–2000s (Figure 2). Recent catches have been at historically high levels. Pink salmon catches varied from 200 million to ,100 million fish, fluctuating annually with a 2-year cycle. Because the pink salmon has a 2-year generation time, a dominant brood line tends to occur in either odd- or even-numbered years. Chum salmon catches were high in the 1930–1940s and 1990–2000s at 50 and 100 million fish, respectively, and low in the 1960–1970s at 20 million fish. No apparent cycles exist in the total catch of chum salmon. Significant correlations were detected between RCs and CIs/ RFs, except for chum salmon stocks off Japan (Table 1). The strongest correlation was observed between RCs and RFs, but season and RFs varied among regions. In six of nine salmon stocks, the strongest correlation was observed between RC and seasonal precipitation, but the sign of the correlation was inconsistent. In the other three stocks, the strongest correlation was observed between RC and seasonal thermal environment (i.e. AT or SST). This indicates that regional salmon catches were not driven by a common RF or CI, but were affected by different RFs in different regions. However, RCs of both species correlated with a common RF in NWBC (i.e. spring precipitation) and in SEAK (i.e. summer precipitation), which suggests that similar mechanisms might drive both chum and pink salmon catches on a regional scale. The estimated effects of large-scale CIs on a given RC were weak (Figure 3). RCs were combined into two groups, i.e. southern stocks, including both species in Japan, NWBC, and SEAK, and northern stocks, including both species in WCAK and chum salmon in Russia. Because CIs were correlated with each other, common temporal patterns were extracted from CIs using two latent variables, i.e. a combined PDO index and a combined VI. The ALPI was correlated with both PDO and VI indices. A path from a combined PDO index to southern stocks remained in the final model, but was non-significant. No paths from any CI to northern stocks remained in the final model. Climate effects were inconsistent among stocks, even in southern or northern areas of the North Pacific (Table 2). Some Table 1. Significant Pearson’s correlation coefficient (r) between RC of chum and pink salmon and large-scale CIs or RFs in five regions around the North Pacific: Japan, Russia, NWBC, SEAK, and WCAK during the 1925– 2001 brood years (i.e. 77 years, the maximum d.f. is 75). Region Japan Russia Species Chum Pink Chum Pink NWBC Chum Pink SEAK Chum Pink WCAK Chum Pink Figure 2. Commercial catch of chum and pink salmon in number of fish around the North Pacific. CI/RF – Winter AT Winter precipitation Summer VI Summer AT Summer SST Spring precipitation Spring precipitation Summer PDO Spring PDO Summer SST Summer precipitation Spring AT Spring PDO Summer precipitation Summer precipitation Summer SST Spring SST Spring AT Adjusted d.f. p-value 0.256 20.359 75 47 ,0.05 ,0.05 0.339 20.290 20.274 20.262 49 73 65 75 ,0.05 ,0.05 ,0.05 ,0.05 20.339 75 ,0.01 20.303 20.277 20.258 0.318 59 75 66 64 ,0.05 ,0.05 ,0.05 ,0.01 0.315 0.290 0.332 47 47 68 ,0.05 ,0.05 ,0.01 20.231 75 ,0.05 0.367 0.286 0.266 75 75 71 ,0.01 ,0.05 ,0.05 r Statistical significance of the correlation was tested by Pyper and Peterman’s (1998) t-test adjusting d.f. –, no CIs or RFs with p , 0.05. 1134 M. Fukuwaka et al. The signs of indirect and direct path coefficients of CIs to RCs were also inconsistent: that of ALPI was positive in two stocks, that of PDO was negative in two stocks, but positive in one stock, and that of VI was negative in four stocks and positive in one stock. The inconsistent effects of large-scale climate change indicated that no single CI or RF controls salmon abundance in all regions of the North Pacific. This demonstrates the importance of using RFs when predicting the effects of climate change on regional patterns of salmon stock survival and productivity. Discussion How do climate change effects on chum and pink salmon vary among regions? Effects of large-scale climate change Figure 3. Path diagram indicating relationships among large-scale CIs and RCs of chum and pink salmon around the North Pacific: Japan, Russia, NWBC, SEAK, WCAK. A rectangle indicates an observed variable and an ellipse indicates a latent variable. An arrow indicates a significant path with a solid line or an insignificant path with a broken line. Numbers beside arrows illustrate path coefficients: ***p ≤ 0.001, **p ≤ 0.01, *p ≤ 0.05, n.s.p . 0.05. Errors for variables are omitted in this figure. Table 2. Compound path coefficients of indirect and direct effects from large-scale CIs to rates of change in RC of chum and pink salmon through RFs around the North Pacific. Region Russia Species Chum WCAK Pink Chum Pink RF Winter precipitation 2 Summer AT 2 Summer SST Summer SST Japan Pink NWBC Chum Pink SEAK Chum Pink Winter AT 2 Spring precipitation Spring precipitation 2 Summer precipitation Spring AT Spring AT Summer precipitation CI ALPI Compound path 0.025 Summer VI – – Summer PDO Summer VI 0.126 – – 0.140 20.045 Winter VI ALPI Spring PDO 20.145 0.053 20.040 Spring PDO 20.044 Summer PDO Spring VI 20.174 20.045 Spring PDO Autumn VI Spring VI 0.078 20.042 20.063 Dash in the RF column indicates a direct effect from a CI. Northern stocks (see the Results section) are displayed above the horizontal line and southern stocks are displayed below the line. CIs and RFs that correlated with RCs in Table 1 were deleted in the SEM processes, because of correlations with other independent variables (Table 2; see also Supplementary material). In seven of the nine salmon stocks, large-scale climate change indexed by CIs indirectly affected RCs through RFs. Direct effects of large-scale climate change were also detected in three regional salmon stocks. Environmental changes in early life could affect the population dynamics of chum and pink salmon stocks strongly. Chum and pink salmon have anadromous life histories, where fish live in a freshwater habitat during the embryonic, larval, and juvenile stages and in a marine habitat during juvenile to immature stages (Quinn, 2005). Mortality in the freshwater life stage is higher than during the marine life stage for Pacific salmon (Bradford, 1995). In the marine life stage, juvenile mortality in the early estuarine and coastal periods is much higher than in the later oceanic period for chum and pink salmon (Parker, 1968; Fukuwaka et al., 2010). The range in variation of freshwater and coastal environments could differ greatly among regions of the North Pacific. Therefore, relationships between the population dynamics of salmon and freshwater and coastal environments could also differ among regions. Responses of such regional environmental factors and salmon population dynamics to global-scale climate change are expected to differ among regions. The results of this study, however, indicate that CIs are related directly and indirectly to regional chum and pink salmon stocks. In terrestrial systems, CIs often predict population processes better than local weather patterns. Most measures of local weather used by researchers fail to capture complex associations between bad weather and population processes, such as high rainfall, strong wind, or low temperature over a short period, and pulsed and lagged mortality of Soay sheep (Ovis aries; Hallett et al., 2004). The frequency of extreme weather events could be related to CIs. In marine systems, because local environmental events might continue for longer durations and affect broader areas, the effects of local environmental change on a fish population might be detectable. The recruitment of two commercial pelagic fish, the Japanese sardine (Sardinops melanostictus) and the chub mackerel (Scomber japonicus), was correlated with SST at spawning grounds in the Kuroshio Current region and also with CIs (Yatsu et al., 2005). CIs could affect terrestrial and ocean RFs, which could control the population processes of chum and pink salmon, because of their anadromous nature. Effects of regional-scale climate The differences in catch fluctuations between southern and northern stocks could suggest the importance of the thermal environment to salmon population dynamics. In the Gulf of Alaska, temporal changes in salmon catch in southern areas were opposite to those in northern areas (Hare et al., 1999). This might be explained by opposite responses of salmon recruitment to coastal SST (Mueter et al., 2002b). On a smaller geographical scale along the Japanese coast, the return rate of hatchery-released 1135 Climate controls of pink and chum salmon chum salmon was correlated negatively with SST along the Japan Sea coast of Honshu near the southern limit of chum salmon distribution, but correlated positively with SST along the Okhotsk Sea coast and Nemuro Strait off Hokkaido in northern Japan (Fukuwaka and Suzuki, 2000; Saito and Nagasawa, 2009). The thermal environment might determine local environmental factors controlling regional salmon stocks, such as SST or AT in northern stocks, whereas precipitation and other factors might control southern stocks. Changes in freshwater, estuarine, or coastal environments caused by precipitation were other important factors affecting chum and pink salmon populations. However, the results of this study demonstrated that the effect of precipitation was inconsistent and could act either positively or negatively on chum and pink salmon stocks. During their freshwater life, river discharge within a normal range correlated positively, but extremely high discharge or flooding correlated negatively with Pacific salmon production (Vernon, 1958; Wickett, 1958; Thorne and Ames, 1987; Beamish et al., 1994). Precipitation also affects coastal marine environments, because river discharge supplies freshwater and nutrients from terrestrial systems to the surface waters of coastal seas. During their early marine life, chum salmon juveniles stay within and may utilize riverine plumes as nursery areas (Fukuwaka and Suzuki, 1998). However, the effects of surface sea salinity have been found to be inconsistent, i.e. negative for the total catch of Fraser River pink salmon, but positive for the survival of several chum and pink salmon in British Columbia and Washington State (Vernon, 1958; Blackbourn, 1990). The mechanisms with which precipitation affects salmon populations might differ among stocks and further investigation is needed to clarify the effects of regional aquatic environments on salmon stocks. The relationship between environmental factors and salmon population dynamics could also differ among regions. CIs or RFs controlling chum and pink salmon catches were different among regions. Some RFs correlated with RCs were consistent with climate factors reported in previous studies on regional salmon abundances, i.e. winter AT for Japanese pink, spring to summer precipitation for NWBC chum and pink, and summer coastal SST for WCAK chum and pink (Blackbourn, 1990; Mueter et al., 2002a; Morita et al., 2006). Although the environmental conditions indispensable for salmon populations might be similar among regional stocks, limiting environmental factors probably differ among regional stocks, because the ranges of variable environmental factors differ among regions. Local adaptations could also cause the differences in environmental conditions required by salmon stocks. Such local adaptations in salmonids can develop rapidly, even on a contemporary time-scale (reviewed by Kinnison and Hendry, 2004). Factors masking the relationship between climate changes and salmon dynamics Human activities, such as artificial propagation, urbanization, or water utilization, might affect or mask the relationship between environmental factors and salmon population dynamics. Almost all Japanese chum salmon stocks have been maintained by hatchery release (Hiroi, 1998). Here, we used the rate of change in catch as an index of salmon population growth. Because the numbers of fish released from hatcheries are controlled artificially, the adult – offspring relationship could be violated even if the survival rate after release is affected by climate. Hatchery fish are incubated and reared artificially; they then migrate to the sea soon after their release from hatcheries (Mayama et al., 1982). Hatchery fish might be affected less by freshwater environments than wild fish and improvements in hatchery techniques might not be related to environmental fluctuation. Therefore, we might have failed to detect effects of climate factors on Japanese chum salmon because of improvements in hatchery techniques (Mayama, 1985). Another masking factor could be an artefact of this analysis. We assumed a fixed generation time for chum salmon (i.e. 4 years). However, the age at maturity for chum salmon ranges from 3 to 6 years, with 4-year-old fish predominant (60 –90%; Salo, 1991). Where the proportion of 4-year-old fish was lower, the relationship between CIs or RFs and RCs might be less detectable. Environmental factors controlling salmon population dynamics can differ among regional populations, because of differences in the ranges of environmental fluctuations, environmental requirements by salmon populations, and human influences on salmon populations and/or their environment. However, in this study, we might have failed to detect a relationship between environmental factors and salmon population dynamics for regional populations where hatchery activity is extensive or the age at maturity differs greatly from our assumption. Does recent climate change favour chum and pink salmon? Interdecadal climate changes as quantified by the CIs examined in this study might not be related directly to the recent large catches of chum and pink salmon. Previous reports have pointed out that ALPI and PDO correlated with local zooplankton production, coastal and terrestrial climate, and also with salmon abundance (Beamish and Bouillon, 1993; Mantua et al., 1997). Since the publication of these reports, although chum and pink salmon catches have remained at high levels (Figure 2), ALPI changed greatly in the early 1990s and PDO changed greatly in the early 2000s. In this study, no single CI was found to affect changes in salmon abundance in a similar way across regions. Although we examined a limited number of CIs, we failed to find any interdecadal climate changes favouring chum and pink salmon population dynamics in the North Pacific. Alternatively, warming trends might have contributed to the recent increases in salmon catch in some regions. The major common trend in the Pacific climate has been a persistent SST warming tendency since 1900 (Schwing et al., 2010). In intracentennial trends, SST increased rapidly in the late 1980s to 1990s in the Oyashio and Kuroshio current systems in the western North Pacific and in the 1970s in the Gulf of Alaska and California Current system in the eastern North Pacific. Radchenko et al. (2007) demonstrated that the increase in Russian pink salmon catch was coincident with an increase in the ocean heat content for the 0 –700-m layer. The current results demonstrated that temperature and SST correlated positively with RCs in some regions of the North Pacific. However, in future, rapid warming might affect salmon catches negatively. In southern and warm areas of salmon distribution, chum and pink salmon survival or productivity correlated negatively with SST (Fukuwaka and Suzuki, 2000; Mueter et al., 2002b). The relationships between salmon growth and distribution and water temperature have been shown to be bell-shaped (Morita et al., 2010a, b). The somatic growth and number of chum and pink salmon decreased at temperatures higher than the optimum. In this study, we 1136 analysed linear relationships between RCs and environmental factors using correlation analyses and SEM. If the true response function of an RC to AT or SST is bell-shaped, the positive effects of warming in the past might become negative effects in future. Conclusions Responses of chum and pink salmon stocks to global climate change vary among regions. Interdecadal climate changes might not be related to recent large increases in catches of chum and pink salmon. Because the thermal environment plays an important role in salmon population dynamics, a warming trend in the North Pacific might affect long-term changes in salmon abundance. In addition, global warming could affect not only the thermal environment, but also other freshwater and coastal environmental conditions for Pacific salmon (Bryant, 2009). The mechanisms controlling regional salmon abundance should be better understood to forecast successfully future conditions for Pacific salmon stocks, because responses of salmon stocks to global climate change vary among regions. Supplementary material The following supplementary material is available at the ICES JMS online version of the paper. Table S1 is a correlation matrix for large-scale climate indices used in the SEM for Figure 3; Table S2 is the correlation coefficient between large-scale climate indices (CIs) and regional climate factors (RFs; above the horizontal line) and correlation matrix for RFs (below the horizontal line) used in the SEM for Russian chum and pink salmon in Table 2; Table S3 is the correlation coefficient between large-scale climate indices (CIs) and regional climate factors (RFs; above the horizontal line) and correlation matrix for RFs (below the horizontal line) used in the SEM for chum and pink salmon in the western and central Alaska in Table 2; Table S4 is the correlation coefficient between large-scale climate indices (CIs) and winter air temperature (winter AT) used in the SEM for Japanese pink salmon in Table 2; Table S5 is the correlation coefficient between large-scale climate indices (CIs) and regional climate factors (RFs; above the horizontal line) and correlation matrix for RFs (below the horizontal line) used in the SEM for chum and pink salmon in the US Northwest and British Columbia in Table 2; Table S6 is the correlation coefficient between large-scale climate indices (CIs) and regional climate factors (RFs; above the horizontal line) and correlation matrix for RFs (below the horizontal line) used in the SEM for chum and pink salmon in the Southeast Alaska in Table 2. Acknowledgements We thank members of the NPAFC Working Group on Stock Assessment for compiling salmon catch data. We also thank Dr Skip McKinnell for providing the PDO and VI indices, and Drs McKinnell and James Irvine for their valuable comments on an earlier version of the manuscript. This study was supported financially by the Promotion Programme of International Fisheries Resources Survey from the Fisheries Agency of Japan. References Beamish, R. J., and Bouillon, D. R. 1993. Pacific salmon production trends in relation to climate. Canadian Journal of Fisheries and Aquatic Sciences, 50: 1002– 1016. M. Fukuwaka et al. Beamish, R. J., and Mahnken, C. 2001. A critical size and period hypothesis to explain natural regulation of salmon abundance and the linkage to climate and climate change. Progress in Oceanography, 49: 423– 437. Beamish, R. J., Neville, C. E., and Cass, A. J. 1997. Production of Fraser River sockeye salmon (Oncorhynchus nerka) in relation to decadal-scale changes in the climate and the ocean. Canadian Journal of Fisheries and Aquatic Sciences, 54: 543– 554. Beamish, R. J., Neville, C-E. M., Thomson, B. L., Harrison, P. J., and St John, M. 1994. A relationship between Fraser River discharge and interannual production of Pacific salmon (Oncorhynchus spp.) and Pacific herring (Clupea pallasi) in the Strait of Georgia. Canadian Journal of Fisheries and Aquatic Sciences, 51: 2843– 2855. Blackbourn, D. J. 1990. Comparison of release size and environmental data with the marine survival rates of some wild and enhanced stocks of pink and chum salmon in British Columbia and Washington State. In 14th Northeast Pacific Pink and Chum Salmon Workshop, pp. 82– 87. Ed. by P. A. Knudsen. Washington State Department of Fisheries, Olympia, WA. Bond, N. A., Overland, J. E., Spillane, M., and Stabeno, P. 2003. Recent shifts in the state of the North Pacific. Geophysical Research Letters, 30: 2183 – 2186. Bradford, M. J. 1995. Comparative review of Pacific salmon survival rates. Canadian Journal of Fisheries and Aquatic Sciences, 52: 1327– 1338. Bryant, M. D. 2009. Global climate change and potential effects on Pacific salmonids in freshwater ecosystems of southeast Alaska. Climate Change, 95: 169 – 193. Eggers, D. M., Irvine, J. R., Fukuwaka, M., and Karpenko, V. I. 2003. Catch trends and status of North Pacific Salmon. NPAFC Document No. 723, Rev. 3. 35 pp. http://www.npafc.org. Fox, J., Kramer, A., and Friendly, M. 2009. Sem: Structural Equation Models. R Package version 0.9-19. http://CRAN.R-project.org/ package=sem. Fukuwaka, M., Sato, S., Yamamura, O., Sakai, O., Nagasawa, T., Nishimura, A., and Azumaya, T. 2010. Biomass and mortality of chum salmon in the pelagic Bering Sea. Marine Ecology Progress Series, 403: 219 – 230. Fukuwaka, M., and Suzuki, T. 1998. Role of a riverine plume as a nursery area for chum salmon Oncorhynchus keta. Marine Ecology Progress Series, 173: 289– 297. Fukuwaka, M., and Suzuki, T. 2000. Density-dependence of chum salmon in coastal waters of the Japan Sea. North Pacific Anadromous Fish Commission Bulletin, 2: 75 – 81. Grace, J. B., Anderson, T. M., Olfe, H., and Scheiner, S. M. 2010. On the specification of structural equation models for ecological systems. Ecological Monographs, 80: 67 –87. Hallet, T. B., Coulson, T., Pilkington, J. G., Clutton-Brock, T. H., Pemberton, J. M., and Grenfell, B. T. 2004. Why large-scale climate indices seem to predict ecological processes better than local weather. Nature, 430: 71– 75. Hare, S. R., Mantua, N. J., and Francis, R. C. 1999. Inverse production regimes: Alaska and West Coast Pacific salmon. Fisheries, 24: 6– 14. Hiroi, O. 1998. Historical trends of salmon fisheries and stock conditions in Japan. North Pacific Anadromous Fish Commission Bulletin, 1: 23 – 27. Holtby, L. B., Andersen, B. C., and Kadowaki, R. K. 1990. Importance of smolt size and early ocean growth to interannual variability in marine survival of coho salmon (Oncorhynchus kisutch). Canadian Journal of Fisheries and Aquatic Sciences, 47: 2181– 2194. Irvine, J. R., and Fukuwaka, M. 2011. Pacific salmon abundance trends and climate change. ICES Journal of Marine Science, 68: 1122–1131. Irvine, J. R., Fukuwaka, M., Kaga, T., Park, J. H., Seong, K. B., Kang, S., Karpenko, V., et al. 2009. Pacific salmon status and abundance trends. NPAFC Document 1199, Rev. 1. 153 pp. http://www. npafc.org. Climate controls of pink and chum salmon Kinnison, M. T., and Hendry, A. P. 2004. From macro- to microevolution: tempo and mode in salmonid evolution. In Evolution Illuminated: Salmon and Their Relatives, pp. 208– 231. Ed. by A. P. Hendry, and S. C. Stearns. Oxford University Press, New York. Levanidov, V. Ya. 1964. Salmon population trends in the Amur basin and means of maintaining the stocks. In Lososevoe Khozyaistvo Dal’nego Vostoka. Ed. by E. N. Pavlovskii. Izdatel’stvo Nauka, Moscow, USSR. 381 pp. (Translation from Russian to Japanese; Soren Hokuyo Gyogyo Kankei Bunken-syu, 66: 49 –80.). Mantua, N. J., Hare, S. R., Zhang, Y., Wallace, J. M., and Francis, R. C. 1997. A Pacific-interdecadal climate oscillation with impacts on salmon production. Bulletin of the American Meteorological Society, 78: 1069– 1079. Mayama, H. 1985. Technical innovations in chum salmon enhancement with special reference to fry condition and timing of release. In Proceedings of the Eleventh U.S. – Japan Meeting on Aquaculture, Salmon Enhancement, Tokyo, Japan, 19– 20 October 1982, pp. 83– 86. Ed. by C. J. Sindermann. U.S. Department of Commerce, NOAA Technical Report NMFS 27. Mayama, H., Kato, H., Seki, J., and Shimizu, I. 1982. Studies on the chum salmon released in the Ishikari River System—I. On the seaward migration and inshore distributions of liberated fry in 1979. Scientific Reports of the Hokkaido Salmon Hatchery, 36: 1 – 7. Mitchell, T. D., and Jones, P. D. 2005. An improved method of constructing a database of monthly climate observations and associated high-resolution grids. International Journal of Climatology, 25: 693– 712. Morita, K., Fukuwaka, M., and Tanimata, N. 2010b. Age-related thermal habitat use by Pacific salmon Oncorhynchus spp. Journal of Fish Biology, 77: 1024– 1029. Morita, K., Fukuwaka, M., Tanimata, N., and Yamamura, O. 2010a. Size-dependent thermal preferences in a pelagic fish. Oikos, 119: 1265– 1272. Morita, K., Morita, S. H., and Fukuwaka, M. 2006. Population dynamics of Japanese pink salmon (Oncorhynchus gorbuscha): are recent increases explained by hatchery programs or climatic variations? Canadian Journal of Fisheries and Aquatic Sciences, 63: 55 – 62. Mueter, F. J., Peterman, R. M., and Pyper, B. J. 2002b. Opposite effects of ocean temperature on survival rates of 120 stocks of Pacific salmon (Oncorhynchus spp.) in northern and southern areas. Canadian Journal of Fisheries and Aquatic Sciences, 59: 456 – 463. Mueter, F. J., Ware, D. M., and Peterman, R. M. 2002a. Spatial correlation patterns in coastal environmental variables and survival rates of salmon in the north-east Pacific Ocean. Fisheries Oceanography, 11: 205– 218. Nagasawa, K. 2000. Long-term changes in the climate and ocean environment in the Okhotsk Sea and western North Pacific and abundance and body weight of East Sakhalin pink salmon (Oncorhynchus gorbuscha). North Pacific Anadromous Fish Commission Bulletin, 2: 203 – 211. Parker, R. R. 1962. Estimations of ocean mortality rates for Pacific salmon (Oncorhynchus). Journal of the Fisheries Research Board of Canada, 19: 561 –589. 1137 Parker, R. R. 1968. Marine mortality schedules of pink salmon of the Bella Coola River, central British Columbia. Journal of the Fisheries Research Board of Canada, 25: 757– 794. Peterman, R. M., Pyper, B. J., Lapointe, M. F., Adkison, M. D., and Walters, C. J. 1998. Patterns of covariation in survival rates of British Columbia and Alaskan sockeye salmon (Oncorhynchus nerka) stocks. Canadian Journal of Fisheries and Aquatic Sciences, 55: 2503– 2517. Pyper, B. J., Mueter, F. J., and Peterman, R. M. 2002. Spatial covariation in survival rates of northeast Pacific chum salmon. Transactions of the American Fisheries Society, 131: 343 – 363. Pyper, B. J., Mueter, F. J., Peterman, R. M., Blackbourn, D. J., and Wood, C. C. 2001. Spatial covariation on survival rates of northeast Pacific pink salmon (Oncorhynchus gorbuscha). Canadian Journal of Fisheries and Aquatic Sciences, 58: 1501– 1515. Pyper, B. J., and Peterman, R. M. 1998. Comparison of methods to account for autocorrelation in correlation analyses of fish data. Canadian Journal of Fisheries and Aquatic Sciences, 55: 2127– 2140. Quinn, T. P. 2005. The Behavior and Ecology of Pacific Salmon and Trout. University of Washington Press, Seattle, WA. Radchenko, V. I., Temnykh, O. S., and Lapko, V. V. 2007. Trends in abundance and biological characteristics in pink salmon (Oncorhynchus gorbuscha) in the North Pacific Ocean. North Pacific Anadromous Fish Commission Bulletin, 4: 7– 21. R Development Core Team. 2009. R: a Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org. Saito, T., and Nagasawa, K. 2009. Regional synchrony in return rates of chum salmon (Oncorhynchus keta) in Japan in relation to coastal temperature and size at release. Fisheries Research, 95: 14– 27. Salo, E. O. 1991. Life history of chum salmon (Oncorhynchus keta). In Pacific Salmon Life Histories, pp. 231 – 309. Ed. by C. Groot, and L. Margolis. UBC Press, Vancouver, BC, Canada. Schwing, F. B., Mendelssohn, R., Bograd, S. J., Overland, J. E., Wang, M., and Ito, S. 2010. Climate change, teleconnection patterns, and regional processes forcing marine populations in the Pacific. Journal of Marine Systems, 79: 245– 257. Thorne, R. E., and Ames, J. J. 1987. A note on variability of marine survival of sockeye salmon (Oncorhynchus nerka) and effects of flooding on spawning success. Canadian Journal of Fisheries and Aquatic Sciences, 44: 1791– 1795. Vernon, E. H. 1958. An examination of factors affecting the abundance of pink salmon in the Fraser River. International Pacific Salmon Fisheries Commission Progress Report. International Pacific Salmon Fisheries Commission, New Westminster, BC, Canada. Wickett, W. P. 1958. Review of certain environmental factors affecting the production of pink and chum salmon. Journal of the Fisheries Research Board of Canada, 15: 1103– 1126. Yatsu, A., Watanabe, T., Ishida, M., Sugisaki, H., and Jacobson, L. D. 2005. Environmental effects on recruitment and productivity of Japanese sardine Sardinops melanostictus and chub mackerel Scomber japonicus with recommendations for management. Fisheries Oceanography, 14: 263– 278.