INTRODUCTIONTop

Environmental changes in the northern Gulf of California (NGC) are thought to negatively impact the primary productivity of the region, with adverse consequences on the development of eggs and larvae of numerous marine species (invertebrates and fishes, Álvarez-Borrego 1983Álvarez-Borrego S. 1983. The Gulf of California. In: Ketchum B.H. (ed.), Estuaries and Enclosed Seas, Elsevier, New York, pp. 427-449., Silber 1990Silber G. 1990. Occurrence and Distribution of the Vaquita Phocoena sinus in the Northern Gulf of California. Fish. Bull. 88: 339-346.). Another factor considered to negatively affect marine populations is fishing activity by industrial shrimp and shark fleets, through retention of juveniles and pre-adults. This activity also impacts species of particular conservation interest in the region, such as totoaba (Totoaba macdonaldi) and marine vaquita (Phocoena sinus), due to their status as endangered species (Diario Oficial de la Federación de México 2010Diario Oficial de la Federación de México. 2010. NOM-059-SEMARNAT-2010. Protección ambiental-Especies nativas de México de flora y fauna silvestres-Categorías de riesgo. Lista de especies en riesgo. SEMARNAT, México. 78 pp.). The NGC and the Colorado River Delta currently constitute a natural reserve protected by the Mexican Federal Government and included since 1995 in the National Programme of Natural Protected Areas (Fig. 1).

The relationship between oceanic patterns and variability in fisheries has been documented previously (Hare and Mantua 2002Hare S., Mantua N. 2002. An historical narrative on the Pacific Decadal Oscillation, interdecadal climate variability and ecosystem impacts. Rep. 20th NE Pacific Pink and Chum Workshop, Seattle, WA, March 22, 2001.). In the Gulf of California, fluctuations in sardine fishery are related to intradecadal events such as the El Niño Southern Oscillation (ENSO) (Lluch-Belda et al. 1986Lluch-Belda D., Magallon F., Schwartzlose R. 1986. Large fluctuations in the sardine fishery in the Gulf of California: possible causes. CalCOFI Rep. 27: 136-140.). Lluch-Belda et al. (1992)Lluch-Belda D., Schwartzlose R., Serra R., et al. 1992. Sardine and anchovy regime fluctuations of abundance in four regions of the world oceans: a workshop report. Fish. Oceanogr. 1: 339-347. also observed that the Pacific Decadal Oscillation (PDO) index, has a significant effect on the sardine populations (Arreguín-Sánchez and Martínez-Aguilar 2004Arreguín-Sánchez F., Martínez-Aguilar S. 2004. Manejo adaptativo de la pesquería de sardina, Sardinops caeruleus, del Golfo de California. In: Quiñonez Velázquez C., Elorduy Garay J.F. (eds), Ambiente y pesquería de pelágicos menores en el noroeste de México, IPN, La Paz BCS, México, pp. 1-16. ). Similar relationships occur for other species in the Gulf of California, including fish and shrimp (Lluch-Cota et al. 2007Lluch-Cota S.E., Aragón-Noriega E.A., Arreguín-Sánchez F., et al. 2007. The Gulf of California: review of ecosystem status and sustainability challenges. Prog. Oceanogr. 73: 1-26., Arreguín-Sánchez et al. 2017Arreguín Sánchez F., Del Monte L.P., Zetina-Rejón M., et al. 2017. The Gulf of California Large Marine Ecosystem: Fisheries and other natural resources. Env. Dev. 22: 71-77.), primary producers (Espinosa-Carreón et al. 2004Espinosa-Carreón T., Strub P., Beier E., et al. 2004. Seasonal and interannual variability of satellite-derived chlorophyll pigment, surface height, and temperature off Baja California. J. Geo. Res. 109: 1-20.), totoaba (Lercari and Chávez 2007Lercari D., Chávez E. 2007. Possible causes related to historic stock depletion of the totoaba, Totoaba macdonaldi (Perciformes: Sciaenidae), endemic to the Gulf of California. Fish. Res. 86: 136-142.), giant squid (Nevárez Martínez et al. 2002Nevárez Martínez M., Rivera Parra G., Morales Bojórquez E., et al. 2002. The jumbo squid (Dosidicus gigas) fishery of the Gulf of California and its relation to environmental variability. Inv. Mar. 30: 193-194.) and pelicans (Sydeman et al. 2001Sydeman W., Hester M., Thayer J., et al. 2001. Climate change, reproductive performance and diet composition of marine birds in the southern California Current system, 1969-1997. Prog. Oceanogr. 49: 309-329.). These reports suggest interactions between the physical environment and the ecosystem that affect species throughout their ontogeny (Lluch-Cota et al. 2007Lluch-Cota S.E., Aragón-Noriega E.A., Arreguín-Sánchez F., et al. 2007. The Gulf of California: review of ecosystem status and sustainability challenges. Prog. Oceanogr. 73: 1-26.).

It also has been documented that large-scale oceanic climate indices such as the PDO index and the Multivariate ENSO Index (MEI) (Wolter 2018Wolter K. 2018. MEI Index (Database) from National Centers for Environmental Information. National Oceanic and Atmospheric Administration.) have relevant effects on fisheries in the Gulf of California (Castro-Ortiz and Lluch-Belda 2008Castro-Ortiz J., Lluch-Belda D. 2008. Impacts of interanual environmental variation on the shrimp fishery off the Gulf of California. CalCOFI Rep. 49: 183-190.). The Pacific Ocean has experienced an interannual-to-interdecadal recurring pattern of ocean-atmosphere climate variability centred over the mid-latitude Pacific (Zhang et al. 1997Zhang Y., Wallace J.M., Battisti D.S. 1997. ENSO-like interdecadal variability: 1900-93. J. Clim. 10: 1004-1020., Newman et al. 2003Newman M., Compo G.P., Alexander M.A. 2003. ENSO-forced variability of the Pacific decadal oscillation. J. Clim. 16: 3853-3857., Yu et al. 2017Yu L., Zhong S., Winkler J.A., et al. 2017. Possible connections of the opposite trends in Arctic and Antarctic sea-ice cover. Sci. Rep. 7: 45804.), which has significant influences on the fishery (Litz et al. 2011Litz M.N., Phillips A.J., Brodeur R.D., et al. 2011. Seasonal occurrences of Humboldt squid (Dosidicus gigas) in the northern California Current System. CalCOFI Rep. 52: 97-108., Newman et al. 2016Newman M., Alexander M.A., Ault T.R., et al. 2016. The Pacific decadal oscillation, revisited. J. Clim. 29: 4399-4427.). The PDO is characterized by the two phases of the sea surface temperature (SST) anomaly: positive (warm) and negative (cool) phases. During the positive phase, the water becomes cooler in the northwestern and central Pacific Ocean and warmer in the eastern Pacific Ocean. The SST anomaly pattern is reversed during the negative phase: warmer in the northwestern Pacific Ocean and cooler in the eastern Pacific Ocean (Mantua and Hare 2002Mantua N., Hare S. 2002. The Pacific decadal oscillation. J. Oceanogr. 58: 35-44.). On other hand, the Pacific–North American (PNA) pattern affects coastal sea and continental surface air temperatures, as well as streamflow in major west coast river systems and terrestrial precipitation variability in North and Central America (Yu and Lin 2019Yu B., Lin H. 2019. Modification of the wintertime Pacific-North American pattern related North American climate anomalies by the Asian-Bering-North American teleconnection. Clim. Dyn. 55: 313-328.).

The Gulf of California is a marine system with a complex environment. One of the greatest challenges for resource managers is to understand how much of the observed change in a population’s abundance can be explained by the environment. Understanding these relationships will enhance the probability of success in managing populations for both fishing and conservation. This is especially true if fishing and conservation efforts can be adapted based on observed variability and if species show specific responses to environmental factors in relation to life histories. Therefore, the focus of the present study is to estimate the role of environmental variables in explaining observed variability of the relative abundance of several commercial species in the NGC.

METHODSTop

Catch per unit effort (CPUE) (kilogram/boat, k/b) was used as an indicator for the relative abundance of commercial species (Spare and Venema 1998Spare P. Venema S. 1998. Introduction to Tropical Fish Stock Assessment. FAO Fish. Tech. Pap. 306. FAO, Rome, Italy.). The time series of daily catch and fishing effort for commercial species were obtained from logbooks of the Regional Fisheries Offices (1995 to 2007). Species catches were listed by common commercial names which may include one or several species, such as chano (Chanos chanos, Micropogonias megalops), corvina (Cynoscion othonopterus, C. parvipinnis, C. xanthulus, C. reticulatus, C. nobilis), blue shrimp (Litopenaeus stylirostris), Pacific sierra (Scomberomorus sierra, S. concolor), ray (Rhinobatos productus, Dasyatis brevis, Gymnura marmorata, Urolophus maculatus, U. halleri, Myliobatis californica, Rinopthera steindachneri, Narcine entemedor), shark (Heterodontus francisci, Heterodontus mexicanus, Rhizoprinodon longurio, Mustelus spp., M. lunulatus, Sphyearna lewini, Squatina californica) and totoaba (Totoaba macdonaldi). For the totoaba, the catch time series included 1929 to 1971, when the fishery was active (Arvizu and Chávez 1972Arvizu J., Chávez H. 1972. Sinopsis sobre la biología de la totoaba Cynoscion macdonaldi Gilbert, 1890. FAO Fish. Synop. 108: 1-26. BCS, México.) and the fishing effort corresponded to a reconstructed time series (Lercari and Chávez 2007Lercari D., Chávez E. 2007. Possible causes related to historic stock depletion of the totoaba, Totoaba macdonaldi (Perciformes: Sciaenidae), endemic to the Gulf of California. Fish. Res. 86: 136-142.).

The selection of environmental variables affecting species considered in this study was based on information from scientific publications: time series of chlorophyll and primary productivity (Espinosa-Carreón et al. 2004Espinosa-Carreón T., Strub P., Beier E., et al. 2004. Seasonal and interannual variability of satellite-derived chlorophyll pigment, surface height, and temperature off Baja California. J. Geo. Res. 109: 1-20.), SST (Aragón-Noriega 2007Aragón-Noriega E. 2007. Coupling the reproductive period of blue shrimp Litopenaeus stylirostris Stimpson, 1874 (Decapoda: Penaeidae) and sea surface temperature in the Gulf of California. Rev. Biol. Mar. Oceanogr. 42: 167-175.), ENSO (Lluch-Cota et al. 2007Lluch-Cota S.E., Aragón-Noriega E.A., Arreguín-Sánchez F., et al. 2007. The Gulf of California: review of ecosystem status and sustainability challenges. Prog. Oceanogr. 73: 1-26., 2010Lluch-Cota S., Pares-Sierra A., Magaña-Rueda V., et al. 2010. Changing climate in the Gulf of California. Prog. Oceanogr. 87: 114-126.), the Colorado River discharge (Rowell et al. 2005Rowell K. Flessa K., Dettman D., et al. 2005. The importance of Colorado River flow to nursery habitats of the Gulf corvina (Cynoscion othonopterus). Can. J. Fish. Aqua. Sci. 62: 2874-2885., Pérez-Arvizu et al. 2009Pérez-Arvizu E., Aragón-Noriega E., Carreon L. 2009. Response of the shrimp population in the Upper Gulf of California to fluctuations in discharges of the Colorado River. Crustaceana 82: 615-625., Aragón-Noriega y Calderón-Aguilera 2000Aragón-Noriega E., Calderón-Aguilera L.E. 2000. Does damming of the Colorado River affect the nursery area of blue shrimp Litopenaeus stylirostris (Decapoda: Penaeidae) in the Upper Gulf California. Rev. Biol. Trop. 48: 867-871.), precipitation (Valdez-Muñoz et al. 2010Valdez-Muñoz C., Aragón-Noriega E.A., Ortega-Rubio A., et al. 2010. Distribución y abundancia de juveniles de totoaba, Totoaba macdonaldi y la salinidad del hábitat de crianza. Interciencia 35: 136-139.), the PDO index and the MEI (Lluch-Cota et al. 2007Lluch-Cota S.E., Aragón-Noriega E.A., Arreguín-Sánchez F., et al. 2007. The Gulf of California: review of ecosystem status and sustainability challenges. Prog. Oceanogr. 73: 1-26., 2010Lluch-Cota S., Pares-Sierra A., Magaña-Rueda V., et al. 2010. Changing climate in the Gulf of California. Prog. Oceanogr. 87: 114-126.) and the PNA pattern, (Cabello-Pasini et al. 2003Cabello-Pasini A., Muñiz-Salazar R., Ward D.H. 2003. Annual variations of biomass and photosynthesis in Zostera marina at its southern end of distribution in the North Pacific. Aq. Bot. 76: 31-47.), and the sources are shown in Appendix 1. Generalized additive models (GAMs) were used to identify variables that explain changes in CPUE because they allowed several environmental variables acting on stock abundances (CPUE) to be incorporated during a single time. The resulting GAMs for each group were chosen considering the lower AIC value (Akaike 1973Akaike H. 1973. Maximum likelihood identification of Gaussian autoregressive moving average models. Biometrika 60: 255-265.), and were also validated by the generalized cross validation estimator (GCV), used as a measure of goodness of fit period (Wood and Augustin 2002Wood S., Augustin N. 2002. GAMs with integrated model selection using penalized regression spline and applications to environmental modelling. Ecol. Model. 157: 157-177., Venables and Dichmont 2004Venables W. Dichmont C. 2004. GLMs, GAMs and GLMMs: an overview of theory for applications in fisheries research. Fish. Res. 70: 319-337.). The GAM model was based on R script (R Core Team 2018R Core Team. 2018. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.) version 3.5.2, using the mgcv package version 1.7-5.

Broad-scale climate indices not only affect weather and climate variability worldwide but have also been demonstrated to considerably influence ecological processes (Stenseth et al. 2003Stenseth N.C., Ottersen G., Hurrell J.W., et al. 2003. Studying climate effects on ecology through the use of climate indices: the North Atlantic Oscillation, El Niño Southern Oscillation and beyond. Proc. R. Soc. London B 270: 2087-2096.). To identify possible biological effects related more strongly to long-scale climate patterns than seasonal factors, and to avoid masking biotic response by seasonal effects, CPUE time series were transformed using natural logarithm (x+1), and were also standardized by z-transformation to facilitate visualization (Blanchard et al. 2010Blanchard J.L., Coll M., Trenkel V.M., et al. 2010. Trend analysis of indicators: a comparison of recent changes in the status of marine ecosystems around the world. ICES J. Mar. Sci. 67: 732-744.). Then we conducted a seasonal and trend decomposition with repeated locally weighted scatterplot smoothing (LOESS), using the stl function in the R package (Zuur et al. 2007Zuur A.F., Ieno E.N. Smith G.M. 2007. Analyzing Ecological Data. In: Gail M., Krickeberg K., Samet A., et al., Statistics for Biology and Health. Springer, New York. 672 pp.). We used a large span width of periodic signal (blue shrimp, chano and sierra, 6 months; sharks and corvina, 5 months; rays, 8 months; totoaba, 2 years) to provide an overview of long-term trends. We also used the trend component and verified for uncorrelated noise processes.

De-seasonalized significant environmental and climate time series derived from GAM models were used to test cross-correlation analysis between them and CPUE of all commercial groups, using the ccf function in R (Walton et al. 2017Walton S., Nunn A., Probst W., et al. 2017. Do fish go with the flow? The effects of periodic and episodic flow pulses on 0+ fish biomass in a constrained lowland river. Ecohydrology 10: e1777.). This function is an estimate of the association between two variables (CPUE and the physical environmental variables) at the significant temporal leads and lags (George et al. 2005George E.P.B., Gwilym M.J., Gregory C.R. 2005. Time Series Analysis: Forecasting and Control. Prentice-Hall, New Jersey.). Temporal lags in the relationship between predictors and the response variable were accounted for, because it has been reported that climate indices influence recruitment success of fish species several months and even years later (Kashkooli et al. 2017Kashkooli O.B., Gröger J., Núñez-Riboni I. 2017. Qualitative assessment of climate-driven ecological shifts in the Caspian Sea. PLoS ONE 12: e0176892.). Positive time lags represent a positive correlation between driver and response and could indicate causality. Negative correlation would be represented by time lags opposed in the driver and response variables (Walton et al. 2017Walton S., Nunn A., Probst W., et al. 2017. Do fish go with the flow? The effects of periodic and episodic flow pulses on 0+ fish biomass in a constrained lowland river. Ecohydrology 10: e1777.).

RESULTSTop

For the predictive model, seven groups of exploited species were considered, including nine fish species (including the totoaba), sixteen elasmobranchs and one crustacean species. Table 1 shows model parameters and the significant environmental variables that corresponded to deviances in abundance for different species. Figure 2 shows the observed CPUE values of the species analysed, as well as the resulting fittings through the GAM for each population. The GAM coefficients and the respective standard errors for all populations analysed in this study are presented in Appendix 2. According to the results obtained, SST was statistically significant with respect to the variability in CPUE, except for the totoaba, for which the Colorado River discharge (CRD) explained the greatest variability in abundance. Other environmental variables that showed significant effects on several species were the PDO and the CRD.

For both sharks and rays the statistically significant environmental variable that explained most of the deviance in CPUE was SST, with 23.69% explained for sharks and 52.10% for rays (Table 1). Based on these results, it was further observed that the highest values of relative abundance (based on CPUE) occurred within an SST range of 22°C to 29.5°C for sharks and 16°C to 20°C for rays (Fig. 3A and B, respectively). The PDO index explained 18.36% of the deviance for sharks and 5.10% for rays. For both, greater abundance was observed during the cold phases of the PDO (Fig. 3C and D, respectively). The higher abundances for both sharks and rays occurred between the cold to neutral phases (Fig. 4A and B, respectively). The total deviance explained using the GAM model is given in Table 1, while Figure 2A and B show observed vs estimated CPUE values.

For blue shrimp, the environmental variable with the most significant impact was SST, which contributed 32.7% of the explained deviance, followed by the MEI, which contributed 12.4%. For SST, the greatest abundance of shrimp was observed within the range of 24°C and 30°C; for the MEI, higher CPUE values were observed when the anomaly remains within a range of ±1°C (Fig. 3E and F, respectively). During El Niño event periods, low CPUE values were recorded. The total deviance explained by the GAM model was 45.10% (Table 1). Figure 2C shows the observed CPUE values vs. the CPUE values estimated by the GAM method, and Figure 4C shows the standardized values of the CPUE and the significant environmental variables related to the blue shrimp.

The most significant environmental variable affecting the Pacific sierra was SST, contributing 67.10% of the explained deviance for variability in CPUE, followed by the PNA index, contributing 7.7% (Table 1). For SST, the highest abundance values were observed within the range of 22°C to 27°C (Fig. 3G), coinciding with the negative phases of PNA pattern (Fig. 3H). The GAM model explained a total deviance of 74.80% of the variability in serranid abundance, including both SST and the PNA pattern (Fig. 2D). Figure 4D shows the standardized values of CPUE and the significant environmental variables related to the Pacific sierra.

For sciaenids, including both chanos and corvina, SST explained most of the variability in CPUE. SST explained 63.10% of the deviance in CPUE for chanos and 67.30% for corvina. The results further indicate a positive correlation between SST (up to 25°C) and changes in the abundance of chanos and a negative correlation for SST (values >21°C). Maximum abundance occurred in the range of SST from 19°C to 24°C (Fig. 3I). For the corvina, a negative correlation between abundance and SST was identified for some ranges. CPUE of the corvina showed no changes at temperatures below 18.5°C. When SST was greater than 19°C, the relationship between this and the CPUE of the corvina showed an inverse relationship (Fig. 3J). The CRD was another environmental variable that significantly affected the chanos and contributed 13.2% to the explained deviance of CPUE variability (Fig. 3K); for the corvina, the CRD contributed 13.9%. The CRD showed a positive correlation for the chano, but for the corvina, for discharges lower than 9 m³ s^–1 of water of Colorado River, CPUE of the corvina showed an inverse relationship with the CRD; for discharges higher than 9 m³ s^–1, the abundance of the corvina showed a direct relationship with the CRD (Fig. 3L). The PDO index was also statistically significant for chanos, with 7.7% of the explained deviance (Fig. 3M), while PDO contributed 5.9% of the explained deviance for corvina (Fig. 3N). In both cases, the negative phase of the PDO resulted in favourable conditions for the abundance of both chanos and corvinas. Figure 4E and F illustrate the standardized values of CPUE and the environmental variables for the chanos and corvina, respectively. Finally, the GAMs resulted in a total explained deviance of 84% for chanos and 87.10% for corvina (Table 1, Fig. 2E and F).

For totoaba, the GAM suggests that the CRD explained most of the variability in CPUE, with 42.0% of the deviance, followed by the PDO with 15.2%. The total explained variation in CPUE was 57.20% (Table 1, Fig. 2G). The CRD exhibited a significant and positive relationship with abundance of totoaba at flows lower than 400 m³ s^–1. The highest abundances were observed when the CRD flows ranged between 200 and 500 m³ s^–1 (Fig. 3O). The PDO showed a trend that was similar to that of CPUE (Fig. 3P). These results suggest that the totoaba stock responds positively to the warm phases of the PDO. Figure 4G shows the standardized observed values of CPUE and the significant environmental variables related to GAM for totoaba.

We cross-correlated the CPUE of elasmobranchs, shrimp, sierra and sciaenids from the NGC with SST, the CRD and climate indices. Cross-correlations (Table 2) indicated that elasmobranch CPUE significantly decreased with SST with no time lag. Shark CPUE was positively correlated with SST anomalies three years earlier (lag of –3 year). Shrimp CPUE was also negatively correlated with SST anomalies and with no time lag, while the Pacific sierra was positively correlated with a time lag of –2 years. (Table 2). Abundance of chano and corvina were negatively correlated with SST, with a three- and one-year time lag, respectively (Table 2).

The long-scale climate PDO index was negatively correlated with elasmobranch CPUE, but with no time lag for sharks and a one-year time lag for rays (Table 2). The PDO was significant for sciaenids: negatively for chano and corvina, and positively for totoaba. Cross-correlation shows a two-year time lag for corvina and no time lag for chano. For totoaba the positive correlation also indicates a two-year time lag with the PDO (Table 2). Blue shrimp was the only commercial species that showed a significant correlation with the MEI. It was positive but showed a five-year time lag, while for the PNA pattern the Pacific sierra was the only one that significantly correlated, being positive and with two-year time lag (Table 2). Finally, the Colorado River discharge was positively correlated with all sciaenids, but with different time lags: 1 year for chano, –1 year for corvina and no time lag for totoaba. (Table 2).

DISCUSSIONTop

According to our results, six of the seven groups under study demonstrated a high correlation between abundance and SST, as shown by several authors for sharks in the Mexican Pacific (Soriano Velázquez et al. 2006Soriano-Velásquez S., Acal Sánchez D., Castillo-Géniz J., et al. 2006. Tiburón del Golfo de Tehuantepec. In: Arreguín Sánchez F., Beléndez-Moreno L.F.J., et al. (eds). Sustentabilidad y Pesca Responsable en México. INA-Pesca, SAGARPA, México, pp. 323-360.). This correlation support reports indicating that the SST plays an important role in seasonal migration for several elasmobranch species in coastal waters (Talent 1985Talent L.G. 1985. The occurrence, seasonal distribution, and repro-ductive condition of elasmobranch fishes in Elkhorn Slough, California. Calif. Fish Game 71: 210-219., Wallman and Bennett 2006Wallman H.L., Bennett W.A. 2006. Effects of parturition and feeding on thermal preference of Atlantic stingray, Dasyatis sabina (Lesueur). Environ. Biol. Fish. 75: 259-267.).

For rays, fisheries in the NGC have increased in recent years, but little research has been published regarding seasonal changes in abundance related to temperature (Cudney-Bueno and Turk-Boyer 1998Cudney-Bueno R., Turk-Boyer R. 1998. Pescando entre mareas del Alto Golfo de California: una guía sobre la pesca artesanal, su gente, y sus propuestas de manejo. CEDO, Pto. Peñasco, Sonora, México, 166 pp.). In both cases, considering the species mobility capacities, the literature suggests that the relationships with SST are related to habitat preferences.

For the blue shrimp, SST is the most significant environmental variable affecting stock abundance in the NGC, coinciding with fishing areas and greater abundances, and particularly with the reproductive process (Aragón-Noriega 2007Aragón-Noriega E. 2007. Coupling the reproductive period of blue shrimp Litopenaeus stylirostris Stimpson, 1874 (Decapoda: Penaeidae) and sea surface temperature in the Gulf of California. Rev. Biol. Mar. Oceanogr. 42: 167-175.).

For the above three study cases, cross-correlations suggest that temperature effects occurred at seasonal level (higher cross-correlations at time t=0), probably associated with ENSO effects and with mesoscale variability (McClatchie 2014McClatchie S. 2014. Scales of variability relevant to fisheries in the southern California Current System. In: Regional fisheries oceanography of the California Current System, Springer. pp. 77-138.).

With respect to the Pacific sierra (Scomberomorus sierra), Medina-Gómez (2006)Medina-Gómez S. 2006. Edad y crecimiento de la sierra del Pacífico Scomberomorus sierra (Jordan y Starks 1895) en el Golfo de California, México. PhD thesis, Instituto Politécnico Nacional, CICIMAR, La Paz, Baja California Sur, México. and Valdovinos-Jacobo (2006)Valdovinos-Jacobo L. 2006. Edad, crecimiento y mortalidad de la sierra del golfo Scomberomorus concolor (Lockington, 1879) en el Golfo de California. M.Sc. thesis, CICIMAR-IPN, México, 60 pp. suggested that migratory behaviour related to feeding and reproduction is governed by changes in temperature, which favoured fishing activities (from October to May). Higher Pacific sierra CPUE coincided with the minimum average values of SST (colder conditions), and the cross-correlation of –2 years with SST is probably related to recruitment success, as shown in the fishery. In contrast, the behaviour of S. concolor, which is an endemic species in the NGC, is governed by the dynamics of the Colorado River Delta, particularly in relation to the reproductive process (Valdovinos-Jacobo 2006Valdovinos-Jacobo L. 2006. Edad, crecimiento y mortalidad de la sierra del golfo Scomberomorus concolor (Lockington, 1879) en el Golfo de California. M.Sc. thesis, CICIMAR-IPN, México, 60 pp.), suggesting a synchrony with warm waters favouring larval development (Moser et al. 1974Moser H.G., Ahlstrom E.H., Kramer D., et al. 1974. Distribution and abundance of fish eggs and larvae in the Gulf of California. CalCOFI Rep. 17: 112-128.).

It is well known that sciaenid fish such as chano (milkfish) and corvina (croaker) inhabit coastal waters, estuaries, rivers and deltas mainly during the reproductive season (Chao 1995Chao L. 1995. Sciaenidae. Corvinas, barbiches, bombaches, corvinatas, corvinetas, corvinillas, lambes, pescadillas, roncachos, verrugatos. In: Fischer W., Krupp F., et al. (eds), Guía FAO para la identificación de especies para los fines de la pesca. Pacifico centro-oriental. FAO, Roma, pp. 1427-1518.). In the NGC, the milkfish shows spawning reproductive aggregations coinciding with warmer waters (CONANP 2007Comisión Nacional de Áreas Naturales Protegidas (CONANP). 2007. Programa de Conservación y Manejo Reserva de la Biosfera Alto Golfo de California y Delta del Río Colorado. CONANP, México DF, 323 pp.), as shown in Figure 3I, while the lag observed by cross-correlation suggests that SST is related to adults during the reproductive process (Table 2). In contrast, corvina (Cynoscion othonopterus), which is widely distributed in the Gulf of California, shows an inverse relationship with SST (Fig. 3J), coinciding with that reported by Chao (1995)Chao L. 1995. Sciaenidae. Corvinas, barbiches, bombaches, corvinatas, corvinetas, corvinillas, lambes, pescadillas, roncachos, verrugatos. In: Fischer W., Krupp F., et al. (eds), Guía FAO para la identificación de especies para los fines de la pesca. Pacifico centro-oriental. FAO, Roma, pp. 1427-1518., while cross-correlation indicates higher value with only a one-year time lag (Table 2), probably reflecting the wide distributions of the species and habitat preference during the reproductive migration period in spring (Pérez-Valencia et al. 2012Pérez-Valencia S., Gorostieta-Monjaraz M., Castañeda-Fernández de Lara V., et al. 2012. Manifestación de Impacto Ambiental para la Pesca Ribereña Responsable en la Reserva de la Biosfera Alto Golfo de California y Delta del Río Colorado: Costa Oeste. CEDO, A.C., Puerto Peñasco, Sonora, México, 221 pp.).

The GAM model indicated that shark and ray catches in the NGC significantly decreased in the same year as a high PDO (and increased after a low PDO anomaly) (Fig. 3C and D respectively). The influence of the PDO index during the shark fishery season could be explained by the fact that SST decreases because of the influence of cold weather and strong winter winds in the northern Pacific Ocean (Castro-Ortiz and Lluch-Belda 2008Castro-Ortiz J., Lluch-Belda D. 2008. Impacts of interanual environmental variation on the shrimp fishery off the Gulf of California. CalCOFI Rep. 49: 183-190.). Differences between the high CPUE for ray in spring and summer found by Ramirez-Amaro et al. (2013)Ramirez-Amaro S.R., Cartamil D., Galvan-Magaña F., et al. 2013. The artisanal elasmobranch fishery of the Pacific coast of Baja California Sur, Mexico, management implications. Sci. Mar. 77: 473-487. and our results suggest that considerable fishing effort may be opportunistically directed.

The influence of the PDO in the dynamics of the NGC has been described by several authors (e.g. Hare and Mantua 2002Hare S., Mantua N. 2002. An historical narrative on the Pacific Decadal Oscillation, interdecadal climate variability and ecosystem impacts. Rep. 20th NE Pacific Pink and Chum Workshop, Seattle, WA, March 22, 2001.), who pointed out that the PDO’s positive phase is characterized by a higher SST in the Alaskan current in the summer months. These events increased the CRD (through the drainage of dams) during the periods 1979-1981 and 1983-1987 (Lavín and Sánchez 1999Lavín M., Sánchez S. 1999. On how the Colorado River affected the hydrography of the upper Gulf of California. Cont. Shelf Res. 19: 1545-1560., Rowell et al. 2005Rowell K. Flessa K., Dettman D., et al. 2005. The importance of Colorado River flow to nursery habitats of the Gulf corvina (Cynoscion othonopterus). Can. J. Fish. Aqua. Sci. 62: 2874-2885.) (Fig. 5), and were associated with the increase in abundance of marine species in the NGC, including the corvina, whose abundance has been increasing since 1992 (Rowell et al. 2005Rowell K. Flessa K., Dettman D., et al. 2005. The importance of Colorado River flow to nursery habitats of the Gulf corvina (Cynoscion othonopterus). Can. J. Fish. Aqua. Sci. 62: 2874-2885.). GAM models for chano and corvina show a negative correlation between CPUE and PDO anomalies (Fig. 3M and N).

Totoaba stock responds directly to the warm and cold phases of the PDO, showing a significant delay with PDO anomalies of –2 years. (Table 2). In this regards, Fischer et al. (1995)Fischer W., Krupp F., Schneider C., et al. 1995. Guía FAO para la identificación de especies para los fines de la pesca. P.C.O., FAO, Roma I: 1-646. indicate that wind patterns and upwellings affect reproduction (Fig. 3P). According to Cury and Roy (1989)Cury P., Roy C. 1989. Optimal environmental window and pelagic fish recruitment success in upwelling areas. Can. J. Fish. Aquat. Sci. 46: 670-680., the synchrony between reproductive success and upwelling intensity can promote the success or failure of fish larva survival. In general, totoaba CPUE coincided with large-scale annual PDO and CRD trends (Fig. 4G).

Castro-Ortiz and Lluch-Belda (2008)Castro-Ortiz J., Lluch-Belda D. 2008. Impacts of interanual environmental variation on the shrimp fishery off the Gulf of California. CalCOFI Rep. 49: 183-190. reported the influence of the PDO on the declining of blue shrimp catch between 1988 and 1991, related to the cold period (Lavín and Sánchez 1999Lavín M., Sánchez S. 1999. On how the Colorado River affected the hydrography of the upper Gulf of California. Cont. Shelf Res. 19: 1545-1560.), and a steadily rising trend that peaked in 1996-1997. Our results suggest that longer blue shrimp reproductive periods occur when MEI anomalies are within the range of ±1 (called the neutral phase, Figs 3F and 4C), decreasing with MEI anomalies greater than ±1. Though blue shrimp CPUE is positively correlated with the MEI with a time lag of –5 year (Table 2), there is no evidence of the process behind this correlation.

High catches of sierra coincide with the negative anomaly of the PNA pattern, associated with colder conditions. It has been suggested that the effect of the PNA pattern occurs through changes in the displacement of air masses. and therefore changes in temperature and precipitation (Wallace and Gutzler 1981Wallace J., Gutzler D. 1981. Teleconnections in the geopotential height field during the Northern Hemisphere winter. Mon. Weather Rev. 109: 784-812.), and is frequently associated with El Niño-La Niña events (Wise 2012Wise E.K. 2012. Hydroclimatology of the US Intermountain West. Prog. Phys. Geogr. 36: 458-479.). The GAM model suggests that higher stock abundance values occurred during the negative and neutral phases of the PNA pattern, as evidenced by adult fishes having a greater presence in cold temperatures (Fig. 3H).

Cudney-Bueno and Turk-Boyer (1998)Cudney-Bueno R., Turk-Boyer R. 1998. Pescando entre mareas del Alto Golfo de California: una guía sobre la pesca artesanal, su gente, y sus propuestas de manejo. CEDO, Pto. Peñasco, Sonora, México, 166 pp. reported that corvina abundance seemed to increase with CRD increments during the period 1979-1988. This increase also corresponded to a decreasing abundance of totoaba, suggesting a species replacement process because of their similar feeding and reproductive niches. Our results in GAM models prove a positive correlation between sciaenid CPUE and CRD (Fig. 3K, L and O respectively), but with different time lags between species (Table 2). Correlations found between CRD and CPUE in sciaenids in the NGC are explained by the fact that chano, corvina and totoaba are endemic species that use the delta region for reproductive and breeding aggregations because they require estuarine conditions. The importance of CRD as an essentially environmental factor have been well documented (Gillanders and Kingsford 2002Gillanders B., Kingsford M. 2002. Impact of changes in flow of freshwater on estuarine and open coastal habitats and the associated organisms. Oceanogr. Mar. Biol. Annu. Rev. 40: 233-309.), while Rowell et al. (2005Rowell K. Flessa K., Dettman D., et al. 2005. The importance of Colorado River flow to nursery habitats of the Gulf corvina (Cynoscion othonopterus). Can. J. Fish. Aqua. Sci. 62: 2874-2885., 2008)Rowell K., True C., Flessa K., et al. 2008. Fish without water: validation and application of the δ18O in Totoaba macdonaldi otoliths. Cien. Mar. 34: 55-68. suggested that higher CRD would increase spawning and nursery habitats, thus benefitting recruitment. Differences in years’ time lags could be explained as particular adaptive processes according to life histories.

Our analysis of totoaba only cover the period in which totoaba fishing was permitted (1925 to 1975). During this period, the CRD represented an environmental variable with a strong and significant effect on stock abundance. At present, the totoaba is considered a critically endangered species by the IUCN red list (Findley 2010Findley L. 2010. Totoaba macdonaldi. The IUCN Red List of Threatened Species 2010: e.T22003A9346099.), and the probable causes mentioned are changes in the CRD, spawning and nursery habitat losses, commercial overfishing and illegal captures, all of them interfering with the success of recovery management measures (Cisneros-Mata et al. 1995Cisneros-Mata M., Montemayor-López G., Román-Rodríguez M. 1995. Life history and conservation of Totoaba macdonaldi. Conserv. Biol. 9: 806-814., Rowell et al. 2008Rowell K., True C., Flessa K., et al. 2008. Fish without water: validation and application of the δ¹⁸O in Totoaba macdonaldi otoliths. Cien. Mar. 34: 55-68.). Since 1943, the CRD has been considerably reduced, coinciding with a very marked decrease in the abundance of totoaba. Unfortunately, present stock regulations (permanent closure) do not permit further investigation into the current effects of climate variables (due to the absence of abundance data), since it is expected that the habitat has changed over the last 30 to 40 years.

CONCLUSIONSTop

The NGC is a highly dynamic ecosystem influenced by several environmental variables that affect species differently in accordance with their life history. However, with the exception of the totoaba, temperature appears as a key environmental variable affecting stock abundance. PDO influences elasmobranches and sciaenid fishes, suggesting effects on species with high mobility, while the CRD affects sciaenid fishes, suggesting, in accordance with the literature, an impact on the reproductive habitat and/or behaviour. Because of its status as a critically endangered species, results are especially relevant for the totoaba. In general terms, the time lags observed between environmental variables and stock abundance seem to be explained through the traits or processes of the life histories of the species and their adaptation to habitat conditions. Our results demonstrate that environmental variables can explain some changes in stock abundance and must be taken into account in the expectations of the management policies of exploited stocks, particularly those aimed at stock recovery of endangered species.

ACKNOWLEDGEMENTSTop

The authors would like to acknowledge the support received through the projects CONACyT 221705 and SIP-IPN 20180929. FAS also thanks the PIA (18714), EDI and COFAA programmes, TMRB thankfully acknowledges fellowships from CONACyT 353507/250618, the PhD programme of the Instituto Politécnico Nacional (CICIMAR) and the BEIFI programme. Dr Juan Antonio de Anda Montañez is thanked for his valuable advice on the GAM statistical method and the most suitable guidelines to follow. The authors also appreciate the constructive comments provided by anonymous referees and by Dr Vinod Sasidharan, which improved the original manuscript.

REFERENCESTop

Akaike H. 1973. Maximum likelihood identification of Gaussian autoregressive moving average models. Biometrika 60: 255-265.
https://doi.org/10.1093/biomet/60.2.255

Álvarez-Borrego S. 1983. The Gulf of California. In: Ketchum B.H. (ed.), Estuaries and Enclosed Seas, Elsevier, New York, pp. 427-449.

Aragón-Noriega E. 2007. Coupling the reproductive period of blue shrimp Litopenaeus stylirostris Stimpson, 1874 (Decapoda: Penaeidae) and sea surface temperature in the Gulf of California. Rev. Biol. Mar. Oceanogr. 42: 167-175.
https://doi.org/10.4067/S0718-19572007000200004

Aragón-Noriega E., Calderón-Aguilera L.E. 2000. Does damming of the Colorado River affect the nursery area of blue shrimp Litopenaeus stylirostris (Decapoda: Penaeidae) in the Upper Gulf California. Rev. Biol. Trop. 48: 867-871.

Arreguín Sánchez F., Del Monte L.P., Zetina-Rejón M., et al. 2017. The Gulf of California Large Marine Ecosystem: Fisheries and other natural resources. Env. Dev. 22: 71-77.
https://doi.org/10.1016/j.envdev.2017.03.002

Arreguín-Sánchez F., Martínez-Aguilar S. 2004. Manejo adaptativo de la pesquería de sardina, Sardinops caeruleus, del Golfo de California. In: Quiñonez Velázquez C., Elorduy Garay J.F. (eds), Ambiente y pesquería de pelágicos menores en el noroeste de México, IPN, La Paz BCS, México, pp. 1-16.

Arvizu J., Chávez H. 1972. Sinopsis sobre la biología de la totoaba Cynoscion macdonaldi Gilbert, 1890. FAO Fish. Synop. 108: 1-26. BCS, México.

Blanchard J.L., Coll M., Trenkel V.M., et al. 2010. Trend analysis of indicators: a comparison of recent changes in the status of marine ecosystems around the world. ICES J. Mar. Sci. 67: 732-744.
https://doi.org/10.1093/icesjms/fsp282

Cabello-Pasini A., Muñiz-Salazar R., Ward D.H. 2003. Annual variations of biomass and photosynthesis in Zostera marina at its southern end of distribution in the North Pacific. Aq. Bot. 76: 31-47.
https://doi.org/10.1016/S0304-3770(03)00012-3

Castro-Ortiz J., Lluch-Belda D. 2008. Impacts of interanual environmental variation on the shrimp fishery off the Gulf of California. CalCOFI Rep. 49: 183-190.

Chao L. 1995. Sciaenidae. Corvinas, barbiches, bombaches, corvinatas, corvinetas, corvinillas, lambes, pescadillas, roncachos, verrugatos. In: Fischer W., Krupp F., et al. (eds), Guía FAO para la identificación de especies para los fines de la pesca. Pacifico centro-oriental. FAO, Roma, pp. 1427-1518.

Cisneros-Mata M., Montemayor-López G., Román-Rodríguez M. 1995. Life history and conservation of Totoaba macdonaldi. Conserv. Biol. 9: 806-814.
https://doi.org/10.1046/j.1523-1739.1995.09040806.x

Comisión Nacional de Áreas Naturales Protegidas (CONANP). 2007. Programa de Conservación y Manejo Reserva de la Biosfera Alto Golfo de California y Delta del Río Colorado. CONANP, México DF, 323 pp.

Cudney-Bueno R., Turk-Boyer R. 1998. Pescando entre mareas del Alto Golfo de California: una guía sobre la pesca artesanal, su gente, y sus propuestas de manejo. CEDO, Pto. Peñasco, Sonora, México, 166 pp.

Cury P., Roy C. 1989. Optimal environmental window and pelagic fish recruitment success in upwelling areas. Can. J. Fish. Aquat. Sci. 46: 670-680.
https://doi.org/10.1139/f89-086

Diario Oficial de la Federación de México. 2010. NOM-059-SEMARNAT-2010. Protección ambiental-Especies nativas de México de flora y fauna silvestres-Categorías de riesgo. Lista de especies en riesgo. SEMARNAT, México. 78 pp.

Espinosa-Carreón T., Strub P., Beier E., et al. 2004. Seasonal and interannual variability of satellite-derived chlorophyll pigment, surface height, and temperature off Baja California. J. Geo. Res. 109: 1-20.
https://doi.org/10.1029/2003JC002105

Findley L. 2010. Totoaba macdonaldi. The IUCN Red List of Threatened Species 2010: e.T22003A9346099.
http://doi.org/10.2305/IUCN.UK.2010-3.RLTS.T22003A9346099.en

Fischer W., Krupp F., Schneider C., et al. 1995. Guía FAO para la identificación de especies para los fines de la pesca. P.C.O., FAO, Roma I: 1-646.

George E.P.B., Gwilym M.J., Gregory C.R. 2005. Time Series Analysis: Forecasting and Control. Prentice-Hall, New Jersey.

Gillanders B., Kingsford M. 2002. Impact of changes in flow of freshwater on estuarine and open coastal habitats and the associated organisms. Oceanogr. Mar. Biol. Annu. Rev. 40: 233-309.

Hare S., Mantua N. 2002. An historical narrative on the Pacific Decadal Oscillation, interdecadal climate variability and ecosystem impacts. Rep. 20th NE Pacific Pink and Chum Workshop, Seattle, WA, March 22, 2001.

Kashkooli O.B., Gröger J., Núñez-Riboni I. 2017. Qualitative assessment of climate-driven ecological shifts in the Caspian Sea. PLoS ONE 12: e0176892.
https://doi.org/10.1371/journal.pone.0176892

Lavín M., Sánchez S. 1999. On how the Colorado River affected the hydrography of the upper Gulf of California. Cont. Shelf Res. 19: 1545-1560.
https://doi.org/10.1016/S0278-4343(99)00030-8

Lercari D., Chávez E. 2007. Possible causes related to historic stock depletion of the totoaba, Totoaba macdonaldi (Perciformes: Sciaenidae), endemic to the Gulf of California. Fish. Res. 86: 136-142.
https://doi.org/10.1016/j.fishres.2007.05.010

Litz M.N., Phillips A.J., Brodeur R.D., et al. 2011. Seasonal occurrences of Humboldt squid (Dosidicus gigas) in the northern California Current System. CalCOFI Rep. 52: 97-108.

Lluch-Belda D., Magallon F., Schwartzlose R. 1986. Large fluctuations in the sardine fishery in the Gulf of California: possible causes. CalCOFI Rep. 27: 136-140.

Lluch-Belda D., Schwartzlose R., Serra R., et al. 1992. Sardine and anchovy regime fluctuations of abundance in four regions of the world oceans: a workshop report. Fish. Oceanogr. 1: 339-347.
https://doi.org/10.1111/j.1365-2419.1992.tb00006.x

Lluch-Cota S., Pares-Sierra A., Magaña-Rueda V., et al. 2010. Changing climate in the Gulf of California. Prog. Oceanogr. 87: 114-126.
https://doi.org/10.1016/j.pocean.2010.09.007

Lluch-Cota S.E., Aragón-Noriega E.A., Arreguín-Sánchez F., et al. 2007. The Gulf of California: review of ecosystem status and sustainability challenges. Prog. Oceanogr. 73: 1-26.
https://doi.org/10.1016/j.pocean.2007.01.013

Mantua N., Hare S. 2002. The Pacific decadal oscillation. J. Oceanogr. 58: 35-44.
https://doi.org/10.1023/A:1015820616384

McClatchie S. 2014. Scales of variability relevant to fisheries in the southern California Current System. In: Regional fisheries oceanography of the California Current System, Springer. pp. 77-138.
https://doi.org/10.1007/978-94-007-7223-6_4

Medina-Gómez S. 2006. Edad y crecimiento de la sierra del Pacífico Scomberomorus sierra (Jordan y Starks 1895) en el Golfo de California, México. PhD thesis, Instituto Politécnico Nacional, CICIMAR, La Paz, Baja California Sur, México.

Moser H.G., Ahlstrom E.H., Kramer D., et al. 1974. Distribution and abundance of fish eggs and larvae in the Gulf of California. CalCOFI Rep. 17: 112-128.

Nevárez Martínez M., Rivera Parra G., Morales Bojórquez E., et al. 2002. The jumbo squid (Dosidicus gigas) fishery of the Gulf of California and its relation to environmental variability. Inv. Mar. 30: 193-194.
https://doi.org/10.4067/S0717-71782002030100082

Newman M., Compo G.P., Alexander M.A. 2003. ENSO-forced variability of the Pacific decadal oscillation. J. Clim. 16: 3853-3857.
https://doi.org/10.1175/1520-0442(2003)016<3853:EVOTPD>2.0.CO;2

Newman M., Alexander M.A., Ault T.R., et al. 2016. The Pacific decadal oscillation, revisited. J. Clim. 29: 4399-4427.
https://doi.org/10.1175/JCLI-D-15-0508.1

Pérez-Arvizu E., Aragón-Noriega E., Carreon L. 2009. Response of the shrimp population in the Upper Gulf of California to fluctuations in discharges of the Colorado River. Crustaceana 82: 615-625.
https://doi.org/10.1163/156854009X404789

Pérez-Valencia S., Gorostieta-Monjaraz M., Castañeda-Fernández de Lara V., et al. 2012. Manifestación de Impacto Ambiental para la Pesca Ribereña Responsable en la Reserva de la Biosfera Alto Golfo de California y Delta del Río Colorado: Costa Oeste. CEDO, A.C., Puerto Peñasco, Sonora, México, 221 pp.

R Core Team. 2018. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
http://www.R-project.org/

Ramirez-Amaro S.R., Cartamil D., Galvan-Magaña F., et al. 2013. The artisanal elasmobranch fishery of the Pacific coast of Baja California Sur, Mexico, management implications. Sci. Mar. 77: 473-487.
https://doi.org/10.3989/scimar.03817.05A

Rowell K. Flessa K., Dettman D., et al. 2005. The importance of Colorado River flow to nursery habitats of the Gulf corvina (Cynoscion othonopterus). Can. J. Fish. Aqua. Sci. 62: 2874-2885.
https://doi.org/10.1139/f05-193

Rowell K., True C., Flessa K., et al. 2008. Fish without water: validation and application of the δ¹⁸O in Totoaba macdonaldi otoliths. Cien. Mar. 34: 55-68.
https://doi.org/10.7773/cm.v34i1.1256

Silber G. 1990. Occurrence and Distribution of the Vaquita Phocoena sinus in the Northern Gulf of California. Fish. Bull. 88: 339-346.

Smith T.M., Reynolds R.W., Peterson T.C., et al. 2008. Improvements to NOAA’s historical merged land-ocean surface temperature analysis (1880-2006). J. Clim. 21: 2283-2296.
https://doi.org/10.1175/2007JCLI2100.1

Soriano-Velásquez S., Acal Sánchez D., Castillo-Géniz J., et al. 2006. Tiburón del Golfo de Tehuantepec. In: Arreguín Sánchez F., Beléndez-Moreno L.F.J., et al. (eds). Sustentabilidad y Pesca Responsable en México. INA-Pesca, SAGARPA, México, pp. 323-360.

Spare P. Venema S. 1998. Introduction to Tropical Fish Stock Assessment. FAO Fish. Tech. Pap. 306. FAO, Rome, Italy.

Stenseth N.C., Ottersen G., Hurrell J.W., et al. 2003. Studying climate effects on ecology through the use of climate indices: the North Atlantic Oscillation, El Niño Southern Oscillation and beyond. Proc. R. Soc. London B 270: 2087-2096.
https://doi.org/10.1098/rspb.2003.2415

Sydeman W., Hester M., Thayer J., et al. 2001. Climate change, reproductive performance and diet composition of marine birds in the southern California Current system, 1969-1997. Prog. Oceanogr. 49: 309-329.
https://doi.org/10.1016/S0079-6611(01)00028-3

Talent L.G. 1985. The occurrence, seasonal distribution, and repro-ductive condition of elasmobranch fishes in Elkhorn Slough, California. Calif. Fish Game 71: 210-219.

Valdez-Muñoz C., Aragón-Noriega E.A., Ortega-Rubio A., et al. 2010. Distribución y abundancia de juveniles de totoaba, Totoaba macdonaldi y la salinidad del hábitat de crianza. Interciencia 35: 136-139.

Valdovinos-Jacobo L. 2006. Edad, crecimiento y mortalidad de la sierra del golfo Scomberomorus concolor (Lockington, 1879) en el Golfo de California. M.Sc. thesis, CICIMAR-IPN, México, 60 pp.

Venables W. Dichmont C. 2004. GLMs, GAMs and GLMMs: an overview of theory for applications in fisheries research. Fish. Res. 70: 319-337.
https://doi.org/10.1016/j.fishres.2004.08.011

Wallace J., Gutzler D. 1981. Teleconnections in the geopotential height field during the Northern Hemisphere winter. Mon. Weather Rev. 109: 784-812.
https://doi.org/10.1175/1520-0493(1981)109<0784:TITGHF>2.0.CO;2

Wallman H.L., Bennett W.A. 2006. Effects of parturition and feeding on thermal preference of Atlantic stingray, Dasyatis sabina (Lesueur). Environ. Biol. Fish. 75: 259-267.
https://doi.org/10.1007/s10641-006-0025-1

Walton S., Nunn A., Probst W., et al. 2017. Do fish go with the flow? The effects of periodic and episodic flow pulses on 0+ fish biomass in a constrained lowland river. Ecohydrology 10: e1777.
https://doi.org/10.1002/eco.1777

Wise E.K. 2012. Hydroclimatology of the US Intermountain West. Prog. Phys. Geogr. 36: 458-479.
https://doi.org/10.1177/0309133312446538

Wolter K. 2018. MEI Index (Database) from National Centers for Environmental Information. National Oceanic and Atmospheric Administration.
https://www.esrl.noaa.gov/psd/enso/mei/

Wood S., Augustin N. 2002. GAMs with integrated model selection using penalized regression spline and applications to environmental modelling. Ecol. Model. 157: 157-177.
https://doi.org/10.1016/S0304-3800(0 2)00193-X

Yu B., Lin H. 2019. Modification of the wintertime Pacific-North American pattern related North American climate anomalies by the Asian-Bering-North American teleconnection. Clim. Dyn. 55: 313-328.
https://doi.org/10.1007/s00382-018-4586-4

Yu L., Zhong S., Winkler J.A., et al. 2017. Possible connections of the opposite trends in Arctic and Antarctic sea-ice cover. Sci. Rep. 7: 45804.
https://doi.org/10.1038/srep45804

Zhang Y., Wallace J.M., Battisti D.S. 1997. ENSO-like interdecadal variability: 1900-93. J. Clim. 10: 1004-1020.
https://doi.org/10.1175/1520-0442(1997)010<1004:ELIV>2.0.CO;2

Zuur A.F., Ieno E.N. Smith G.M. 2007. Analyzing Ecological Data. In: Gail M., Krickeberg K., Samet A., et al., Statistics for Biology and Health. Springer, New York. 672 pp.
https://doi.org/10.1007/978-0-387-45972-1