SciELO - Scientific Electronic Library Online

 
vol.49Rasgos de la vida temprana de la dorada común de dos bandas, Diplodus vulgaris (Perciformes: Sparidae), que habita en las aguas poco profundas de la isla de Gökçeada, TurquíaComposición y biomasa fitoplanctónica bajo condiciones oligotróficas en la cuenca de Guaymas (golfo de California) índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Ciencias marinas

versión impresa ISSN 0185-3880

Cienc. mar vol.49  Ensenada ene./dic. 2023  Epub 08-Dic-2023

https://doi.org/10.7773/cm.y2023.3260 

Articles

Temporal oceanographic variation using satellite imagery data in the central Mexican Pacific convergence zone

Carlos Vladimir Pérez-de Silva1 
http://orcid.org/0000-0002-8200-7411

Amílcar Leví Cupul-Magaña1  * 
http://orcid.org/0000-0002-6455-1253

Fabián Alejandro Rodríguez-Zaragoza2 
http://orcid.org/0000-0002-0066-4275

Alma Paola Rodríguez-Troncoso1 
http://orcid.org/0000-0001-6243-7679

1Laboratorio de Ecología Marina, Departamento de Ciencias Biológicas, Centro Universitario de la Costa, Universidad de Guadalajara, 48280 Puerto Vallarta, Jalisco, Mexico.

2Laboratorio de Ecología Molecular, Microbiología y Taxonomía, Departamento de Ecología, Centro Universitario de Ciencias Biológicas y Agropecuarias, Universidad de Guadalajara, 45200 Zapopan, Jalisco, Mexico.


Abstract

The oceanographic variation of the central Mexican Pacific (CMP) is mainly driven by the California Current and the Mexican Coastal Current, which generate a convergence zone. Little information is available regarding changes in oceanographic variables at the mesoscale level in the CMP. Therefore, this study characterized 6 oceanographic variables (sea surface temperature, chlorophyll concentration [Chl], primary productivity [PP], diffuse attenuation coefficient [K490], and particulate inorganic/organic carbon concentration [PIC, POC]) in the CMP from 2010 to 2017 and their relationships with El Niño/Southern Oscillation (ENSO). The variables were standardized to monthly pixel values of 0.08 latitude degrees, and the study area covered 48,846.48 km2. Friedman tests were used to compare the temporal variation in the variables, while Spearman correlations were used to evaluate the relationship between each variable and the Multivariate ENSO Index (MEI). A cross-correlation analysis was performed to determine the temporal lag between the oceanographic variables and the MEI. The cyclicity of the variation in the CMP was determined by spectral analysis. All variables showed significant differences between months and years. Two seasons defined by temperature were also detected: a cold season (December-June), in which high values of these variables were observed, and a warm season (July-November), in which low values of these variables were observed. No variables were correlated with the MEI; however, a 4-month time lag was identified between the variables and the MEI. The cyclicity of the variables corresponded to the cold and warm seasons. The cold phase of ENSO increased the values of PP, Chl, and K490 up to 4 times compared to those of other years. Taken together, the observed variation in oceanographic conditions makes the CMP one of the most dynamic coastal regions of the Mexican Pacific.

Key words: Mexican coast; oceanographic variability; ENSO

Resumen

La variación oceanográfica del Pacífico central mexicano (PCM) está influenciada por la corriente de California y la corriente costera Mexicana, que generan una zona de convergencia. Existe poca información acerca de los cambios en las variables oceanográficas a nivel de mesoescala en el PCM. Por lo tanto, este estudio caracterizó 6 variables oceanográficas (temperatura superficial del mar, concentración de clorofila [Cl], productividad primaria [PP], coeficiente de atenuación difusa [K490] y concentración de carbono inorgánico/orgánico particulado [CIP, COP]) en el PCM de 2010 a 2017 y su relación con El Niño/Oscilación del Sur (ENOS). Las variables se estandarizaron a valores mensuales en píxeles de 0.08 grados latitudinales, y el área de estudio cubrió 48,846.48 km2. Para comparar la variación temporal en las variables, se utilizó la prueba de Friedman, y la relación entre cada variable y el Índice Multivariado de ENOS (IME) se obtuvo con una correlación de Spearman. Para establecer el desfase temporal entre las variables y el IME, se realizó una correlación cruzada. La ciclicidad de la variación en el PCM se determinó mediante un análisis espectral. Todas las variables mostraron diferencias significativas entre meses y años. Se detectaron 2 temporadas definidas por la temperatura, una temporada fría (diciembre-junio), con valores de las variables altos, y una temporada cálida (julio-noviembre), con valores de las variables bajos. Ninguna variable mostró correlación con el IME; sin embargo, se determinó un lapso de respuesta de 4 meses entre estas y el IME. La ciclicidad de las variables corresponde a la estación fría y la estación cálida. La fase fría del ENOS incrementó los valores de PP, Cl y K490 hasta 4 veces más en comparación con aquellos observados en otros años. La variación observada en las condiciones oceanográficas convierte al PCM en una de las regiones costeras más dinámicas del Pacífico mexicano.

Palabras clave: costa mexicana; variabilidad oceanográfica; ENOS

Introduction

The central Mexican Pacific (CMP) constitutes a marine convergence zone that is mainly influenced by the water masses of the California Current (CC) and Mexican Coastal Current (MCC) (Lavín and Marinone 2003, López-Sandoval et al. 2009b). Changes in temperature, primary productivity, and light attenuation can influence the biological attributes of marine regions (Espinosa-Carreón and Valdez-Holguín 2007, Ribas-Ribas et al. 2011, Speers et al. 2016). Particulate organic and inorganic carbon (POC and PIC, respectively) are both indicators of CO2production, and variations in their content can alter the air-sea flux of CO2at the surface and carbon cycle in the water column (Gafar et al. 2019, Hopkins et al. 2019). Variations in energy flow, temperature, water quality, and carbon availability can affect biodiversity, coastal ecosystems, and ecosystem services (IPCC 2002). Moreover, the convergence of 2 biogeographical regions (Spalding et al. 2007) in the CMP highlights spatiotemporal transitions in the area from both oceanographic and biological perspectives.

Seasonal, annual, and interannual environmental conditions are influenced by mesoscale processes, such as El Niño/Southern Oscillation (ENSO) events and oceanographic patterns, especially over transition areas like the CMP (Maul et al. 1984, Brodeur et al. 2003, Ballance et al. 2006, Palacios et al. 2006, Robertson and Cramer 2009). This affects species richness and community composition as well as oceanographic processes like upwelling (Harley et al. 2006, Carpenter et al. 2008).

ENSO events are the primary source of interannual variation in the CMP (Lavín and Marinone 2003), and the organisms present within this transition zone have adapted to annual cyclical changes in their environment (Green and Fisher 2004). However, anomalous changes in environmental conditions can negatively affect their growth (e.g., reduced carbonate assimilation by calcifying organisms like corals) and reproduction (e.g., reduced fecundity in various fish species; Brander 2010, Wild et al. 2011, Rodríguez-Troncoso et al. 2016). Therefore, the recent changes in the intensity and periodicity of ENSO events may threaten the organisms of the CMP (Palacios et al. 2006). These unexpected changes often occur over short time periods (i.e., a few months), which are insufficient for organisms to adequately respond. In particular, important life cycle stages (e.g., breeding, feeding, and growth) can be severely affected by sudden environmental change (Donelson et al. 2010).

In the CMP, the warm phase of ENSO (El Niño) is characterized by the intrusion of warm southern waters from the central Pacific that increase sea surface temperature and deepen the thermocline and nutricline, followed by a reduction in primary productivity in shallow, coastal waters (McPhaden and Zhang 2002, Páez-Osuna et al. 2016). Regional variation coupled with the recently enhanced ENSO effects have generated a zone of oceanographic instability that affects resident species (Magaña-Rueda et al. 1997, Reyes-Coca and Troncoso-Gaytán 2001, Lee and McPhaden 2010). Within the Gulf of California (GC), oceanographic variability and ENSO effects have been well documented (Alvarez-Borrego et al. 1979, Álvarez-Borrego and Lara-Lara 1991, Lavín and Marinone 2003, Espinosa-Carreón and Valdez-Holguín 2007, Ledesma-Vázquez et al. 2009, Mercado-Santana et al. 2017); however, these phenomena remain understudied in areas outside the GC, such as the CMP (Filonov et al. 2000, López-Sandoval et al. 2009b), even when these areas harbor species of both commercial and scientific interest (Moncayo-Estrada et al. 2006).

Currently, no recent studies have described changes in the oceanographic conditions of the CMP, which are increasing in both intensity and frequency due to climate change (IPCC 2002, di Lorenzo et al. 2010, Lee and McPhaden 2010). Understanding how these changes occur within the CMP is important to evaluate their effects on resident species. Therefore, the aim of this study was to describe the variation in 6 environmental variables within the CMP over a 7-year period, which included 2 intense ENSO events.

Materials and methods

Study area

The CMP comprises the coasts of Michoacán, Colima, Jalisco, and Nayarit. The study area (17-23° N, 108-102° W) covered 48,846.48 km2of oceanic and coastal zones. The CMP is characterized by a narrow continental shelf (~100 km) that further narrows to the south and by depths >2,000 m throughout the region, with the deepest areas reaching 5,000 m. In the northern portion of the study area off Nayarit and Jalisco, the shelf is shallow (<200 m), whereas in the southern area off Colima and Michoacán, the shelf is narrow (<10 km) with deep waters near the coast. Finally, this area includes 3 insular systems, Marias Islands, Isabel Island, and Marietas Islands, which are located on the continental platform of Nayarit in the northern CMP (Reyes-Bonilla 2003).

The CMP is influenced nearly year-round by northwesterly winds (Fiedler and Talley 2006, Pantoja et al. 2012) that increase the number of eddies and meanders, which in turn promote an increase in productivity near the coast (Lavín et al. 2006). In the summer, the CMP is influenced by the poleward MCC, which is characterized by high temperature and nutrient-poor water (Badan-Dangon 1998, Fiedler and Talley 2006). During winter, the CMP is influenced by the CC, which is characterized by cold water with high nutrient and oxygen concentrations (Cepeda-Morales et al. 2009). In addition, the cyclonic season in the CMP generates rainfall from June to October (García-Oliva et al. 1991, 2002).

Data collection

To characterize the temporal oceanographic variation in the CMP, variables that have a major influence on the distribution of resident species were evaluated. These variables were sea surface temperature (SST), primary productivity (PP), chlorophyll concentration (Chl), and water clarity based on the diffuse attenuation coefficient (K490). In addition, the availability of both particulate inorganic carbon (PIC) and particulate organic carbon (POC) were included, as these determine the distribution of zones with high CO2concentrations, which in turn affect the zooplankton distributions throughout the water column (Gafar et al. 2019, Hopkins et al. 2019).

The numerical values of the oceanographic variables for the period 2010 to 2017 were obtained through the ERDDAP platform of the National Oceanic and Atmospheric Administration (NOAA; Simons and Mendelssohn 2012, https://coastwatch.pfeg.noaa.gov/erddap/index.html). Data corresponding to the aforementioned variables collected by the Aqua/MODIS sensor with monthly periodicity and a minimum pixel value of 0.025° latitude were used in the analyses. An atmospheric correction was applied to the data to obtain a measurement of the outflow radiance of the water. These radiance values were processed for each of the environmental variables using the OC3M algorithm developed by the National Aeronautics and Space Administration (NASA). Prior to the analyses, all oceanographic variables were standardized to a pixel resolution of 0.08 × 0.08 latitude degrees (8.45 km2, 5,625.00 pixels per image) and to mean monthly values (total data sets = 94). Finally, the monthly Multivariate ENSO Index (MEI) was obtained from the NOAA physical sciences laboratory (Wolter and Timlin 2011) and used to perform the statistical procedures described below.

Data analyses

As the data did not comply with the assumptions of parametric statistics, nonparametric tests were used. To compare temporal variation among the variables, 2 Friedman tests (nonparametric repeated-measures analysis of variance) were conducted for each variable using 2 sources of primary variation: month-to-month for intrannual variation and year-to-year for interannual variation. To determine the differences between months and years, ana posteriorinonparametric Tukey test was performed. A comparison of anomalies standardized to the mean was conducted based on the monthly average value of each variable and the corresponding yearly average to identify differences with respect to the general trend of each variable.

The variables were standardized to zero mean and were used to establish the relationships between the variables based on nonparametric Spearman correlations, the monthly values of each oceanographic measurement, and the monthly MEI values. In addition, a cross-correlation test was performed to identify the lag between the highest value of each oceanographic variable and the MEI. Finally, a spectral analysis was performed to determine the frequency of variation of the time series corresponding to the 6 variables. All data analyses and graphics were performed using the packages ‘stats’ (statistical analysis) and ‘ggplot2’ and ‘gplots’ (graphing) (R Core Team 2021).

Results

Temporal characterization and anomaly analysis

Sea surface temperature

The mean SST was 28 °C (SD = 2.08, Fig. 1a), with maximum and minimum monthly values of 31.19 °C in July 2014 and 23.34 °C in February 2011, respectively (Fig. 1a). Differences among monthly SST were observed in all months (Table 1) except in February and March, which showed the lowest values (Fig. 2a). The results showed that the region is characterized by a cold season from December to May and a warm season from July to November.

Figure 1.  Mean monthly values of 6 oceanographic variables from 2010 to 2017 in the central Mexican Pacific. Blue shade represents standard deviation values. (a) Sea surface temperature (SST, °C), (b) primary productivity (PP, mg C·m−2·d−1), (c) diffuse attenuation coefficient (K490, m-1), (d) chlorophyll concentration (Chl, mg·m−3), (e) particulate inorganic carbon (PIC, mol·m−3), and (f) particulate organic carbon (POC, mg·m−3). 

Table 1.  Friedman test results showing interannual and monthly differences SST: sea surface temperature; PP: primary productivity; K490: diffuse attenuation coefficient; Chl: chlorophyll concentration; PIC: particulate inorganic carbon concentration; and POC: particulate organic carbon concentration. 

SST PP K490
Factor Friedman χ² Pvalue Friedman χ² Pvalue Friedman χ² Pvalue
Months 38.361 <0.001 44.720 <0.001 46.889 <0.001
Year 79.840 <0.002 60.231 <0.001 59.212 <0.001
Chl PIC POC
Factor Friedman χ² Pvalue Friedman χ² Pvalue Friedman χ² Pvalue
Months 23.750 <0.001 18.303 0.010 43.667 <0.001
Year 36.667 <0.001 23.374 0.010 58.000 <0.001

Figure 2.  Monthly mean values of 6 oceanographic variables from 2010 to 2017 in the central Mexican Pacific. Blue shade represents standard error values. (a) Sea surface temperature (SST, °C), (b) primary productivity (PP, mg C·m−2·d−1), (c) diffuse attenuation coefficient (K490, m-1), (d) chlorophyll concentration (Chl, mg·m−3), (e) particulate inorganic carbon (PIC, mol·m−3), (f) Particulate organic carbon (POC, mg·m−3). Ave: monthly average value of the variable constructed from the years of observations. Names and months with the same letter are similar according to Tukey a posteriori tests. 

Differences were also observed among years (Table 1), with 2011 and 2015 being the coldest and warmest years, respectively (Fig. 3a). Cold years (defined as having an average value lower than the average value of all years) corresponded to 2010-2013, whereas warm years (defined as having an average value higher than the average value of all years) were 2014 and 2015 (Fig. 3a). The analysis of the annual averages with respect to the average SST value for each year indicated that 2010-2012 exhibited the highest proportion of negative variation based on the average monthly values. In contrast, 2014-2016 exhibited the highest proportion of positive variation in SST based on the average monthly values (Fig. 2a).

Figure 3.  Mean annual values of 6 oceanographic variables from 2010 to 2017 in the central Mexican Pacific. Blue shade represents standard error values. (a) Sea surface temperature (SST, °C), (b) primary productivity (PP, mg C·m−2·d−1), (c) diffuse attenuation coefficient (K490, m-1), (d) chlorophyll concentration (Chl, mg·m−3), (e) particulate inorganic carbon (PIC, mol·m−3), (f) particulate organic carbon (POC, mg·m−3). Names and years with the same letter are similar according to Tukey a posteriori tests. 

Primary productivity

The average PP value in the study area was 818.41 mg C·m−2·d−1(SD = 2,511.62), with the highest value observed in February 2011 (4,375.07 mg C·m−2·d−1) and the lowest value observed in November 2016 (319.28 mg C·m−2·d−1) (Fig. 1b). The monthly pattern highlighted an increase in PP during February 2011 that was 5 times greater than the overall average (Fig. 1b). PP showed differences among months, with the highest and lowest values corresponding to January-May and June-December, respectively (Fig. 2b).

When comparing differences among years, 2011 showed the highest PP concentration (1,630.00 mg C·m−2·d−1), which was twice the average of all years. The a posteriori tests revealed significant differences among all years except for 2013, 2014, and 2017 (Fig. 3b). For this variable, the years with the highest proportion of positive monthly anomalies with respect to the mean for all years corresponded to the period of 2010-2013 (Fig. 2b), whereas 2016 predominantly exhibited negative anomalies.

Diffuse attenuation coefficient

Overall, the mean K490 value in the study area was 0.060 m−1(SD = 0.120). The highest and lowest K490 values were observed in April 2011 (0.232 m−1) and July 2012 (0.035 m−1) (Fig. 1c). Monthly comparisons resulted in significant differences (Table 1) and revealed 2 periods: one from January to May with values above the monthly average for all years and the other from June to December with values below the monthly average for all years (Table 1, Fig. 2c).

Differences among years were also observed, with 2010-2012 exhibiting values higher than the total average (Fig. 3c). The largest positive anomalies in K490 were observed from 2010 to 2013 (Fig. 2c). In addition, the negative anomalies in K490 were found closer to the total average than were the positive K490 anomalies, with 2016 exhibiting the greatest negative anomaly value.

Chlorophyll concentration

The mean Chl concentration in the study area was 0.793 mg Chl·m−3(SD = 4.300), with the highest and lowest values observed in May 2012 (2.507 mg Chl·m−3) and June 2013 (0.124 mg Chl·m−3), respectively (Fig. 1d). Monthly comparisons revealed 2 periods: one from February to May characterized by high Chl concentrations and the other from June to December characterized by a decrease in the Chl concentration (Fig. 2d, Table 1). Differences in the Chl concentration among years were less marked than those observed among months, and differences across years were mainly found in 2010, 2011, 2012, and 2014, all of which exhibited the highest mean Chl concentrations (Fig. 3d). The positive Chl anomalies were observed during 2010-2013, whereas the greatest negative Chl anomalies were observed in 2016 (Fig. 2d).

Concentration of particulate inorganic carbon

The mean PIC concentration in the study area was 0.077 mol·m−3(SD = 0.260), with the highest and lowest concentrations observed in April 2011 (0.035 mol·m−3) and December 2012 (0.017 mol·m−3), respectively (Fig. 1e). Pairwise comparisons resulted in significant differences among monthly PIC concentrations. Two peaks were observed each year: one from August to October and the other from April to May (Fig. 2e).

Annual differences were mainly explained by 2011 and 2012, which showed the highest PIC concentrations, and by 2017, which showed the lowest PIC concentrations (Fig. 3e). As with the other variables, the greatest positive anomalies were observed from 2010 to 2013, although no particular pattern was observed with negative monthly anomalies over the study period (Fig. 2e).

Concentration of particulate organic carbon

The mean POC concentration in the study area was 90.52 mg·m−3(SD = 94.38), with February 2011 exhibiting the highest value (252.92 mg·m−3) and July 2016 exhibiting the lowest value (50.78 mg·m−3). Among all years, 2011 and 2012 showed the widest ranges between maximum and minimum monthly values (Fig. 1f). In addition, significant differences were present among months, with higher and lower POC concentrations observed from January to May and June to December, respectively (Fig. 2f). However, no significant differences were observed among the months with either the highest or lowest POC concentrations.

In contrast, significant differences were also observed among years, with 2011 and 2016 showing the highest and lowest POC concentrations, respectively. Overall, the highest POC concentrations were observed from 2010 to 2012, whereas the lowest POC concentrations were observed from 2013 to 2017 (Fig. 3f). Positive POC anomalies were present in 2011 and 2012, whereas 2016 was characterized by negative anomalies (Fig. 2f).

Multivariate El Niño/Southern Oscillation Index

Our analyses indicated that the MEI varied among and within years. From August 2010 to March 2011, the MEI values were lower than 1.5 SD. Thus, this period can be considered a cold phase of an ENSO event (La Niña). In contrast, from May 2015 to May 2016, the MEI values were higher than 1.5 SD, which indicate that this period of time constituted a warm phase of an ENSO event (El Niño). When comparing both the cold and warm phases, the El Niño years exhibited higher variations in both time and MEI values than those of the La Niña years (Fig. 4).

Figure 4.  Multivariate El Niño/Southern Oscillation (ENSO) Index (MEI) values and their intensity during 2010-2017 (data from Wolter and Timlin 2011). 

Relationships among variables

Among all variables, the highest correlations were observed between SST and PP, followed by the correlations between K490 and Chl and K490 and POC (Fig. 5). Both the MEI and PIC showed the lowest correlations with the other variables (Fig. 5). The cross-correlation analysis of SST with MEI revealed a time lag of 3 months between the highest MEI values and the highest SST values (r= 0.041, Table 2). Moreover, the highest MEI values were associated with time lags greater than 3 months in all comparisons (Table 2).

Figure 5.  Correlogram of the environmental variables and Multivariate El Niño/Southern Oscillation Index (MEI). The top right shows the results and significance of the Spearman correlations. The lower left shows the correlation diagrams between the variables. SST: sea surface temperature; PP: primary productivity; K490: diffuse attenuation coefficient; Chl: chlorophyll concentration; PIC: particulate inorganic carbon concentration; and POC: particulate organic carbon concentration. 

Table 2.  Temporal cross-correlation analysis. The first value represents the time lag (months), and the second value represents the correlation value between the higher average values of 6 oceanographic variables. SST: sea surface temperature; PP: primary productivity; K490: diffuse attenuation coefficient; Chl: chlorophyll concentration; PIC: particulate inorganic carbon concentration; POC: particulate organic carbon concentration; MEI: Multivariate El Niño/Southern Oscillation Index. 

PP K490 Chl PIC POC MEI
SST 0.00, -0.69 0.00, -0.70 -1.00, -0.66 -4.00, 0.38 0.00, -0.74 -4.00, 0.41
PP 0.00, 0.98 0.00, 0.75 -2.00, 0.61 0.00, 0.94 4.00, 0.54
K490 0.00, 0.77 -2.00, 0.56 0.00, 0.96 3.00, -0.57
Chl -4.00, 0.37 0.00, 0.76 5.00, -0.48
PIC 2.00, 0.54 7.00, -0.40
POC 4.00, -0.58

Spectral analysis

The spectral analysis revealed that the main frequency influencing the variables in this study was 0.0833333. The inverse of this frequency (1/0.0833333) represents a period of 12, which in this study implies an annual cycle represented by 12 months (Fig. 6). In this context, temperature and PIC had the highest and lowest values of the spectrum, respectively.

Figure 6.  Spectral analysis of the standardized environmental variables used in this study for the central Mexican Pacific. SST: Sea surface temperature; PP: primary productivity; K490: diffuse attenuation coefficient; Chl: chlorophyll concentration; PIC: particulate inorganic carbon; POC: particulate organic carbon. 

Discussion

SST, POC, K490, and Chl are essential variables to the dynamics of coastal and pelagic ecosystems (Palacios et al. 2006, Escalante et al. 2013, Audzijonyte et al. 2020). Unlike SST, which shows little variability based on the distance from the coast, PP, Chl, K490, and POC values increase toward the coast (Zeitzschel 1969, Mercado-Santana et al. 2017). SST is not the only variable that affects environmental conditions, local processes, such as land runoff during the rainy season, coastal upwelling, and human activities, can also modify the environment (Zamudio et al. 2001, García-Oliva et al. 2002). Therefore, other processes beyond SST also contribute to oceanographic environmental variability, especially in coastal areas such as the CMP. In this study, the anomalies were consistent over the months and years in that the greatest anomalies were present during the same period (February-March) each year. During this period in 2011, the CMP was affected by the intense effects of a cold ENSO phase.

Given that the correlations between the variables depend on the characteristics of the water masses that converge in the area during cold (December to May) and warm (June to November) seasons, the CC has a more notable influence on the study area during winter (Palacios et al. 2006), whereas the MCC has a greater influence during summer (Collins et al. 1997). The convergence of these currents in the CMP affects the oceanographic conditions present throughout the year. Our results showed that in months with low-temperature values, the mean values of PP, POC, Chl, and K490 increased, and the greatest Chl concentration was present during the cold season in 2011, which agrees with previous observations of the area (Cepeda-Morales et al. 2017). The presence of the CC, which is characterized by high oxygen concentrations (Cepeda-Morales et al. 2013), coincided with an increase in PP, POC, Chl, and K490 during the cold season from December to May (Collins et al. 1997). In addition, intensified northwesterly winds (López-Sandoval et al. 2009a, b) promote upwelling events in the coastal area during this season. Cold and oxygen-rich water coupled with upwelling promotes microalgae blooms (Anderson et al. 2002, Cepeda-Morales et al. 2009, Cortés-Lara et al. 2012), which increase PP and Chl values in the coastal area during the cold ENSO phase. The increase in the periodicity of ENSO may further extend the cold season with increased nutrient availability and turbidity in the CMP.

During the years included in this study (2010-2017), an El Niño event occurred (2015-2016) during which the influence of the MCC increased in the CMP. This current was characterized by higher temperatures than those of the CC and severely hypoxic waters (Cepeda-Morales et al. 2013). These characteristics, combined with the diminished intensity of the northwesterly winds in the area during the same months (López-Sandoval et al. 2009b), suppressed upwelling events. As oxygen and nutrients are essential for the growth of photosynthetic microorganisms that generate up to 90% of the Chl and biomass in the CMP (Cepeda-Morales et al. 2009), the biological processes in the region may have been severely compromised during this intense and extended warm ENSO phase. The time lags between the variables considered in this study and the effects of El Niño events have been extensively studied, and the results of this study agree with the lag of 1-3 months that was described by Kumar in 2003, which has been observed since 1950. However, this phenomenon has increased in periodicity and intensity as of late, which may be affecting the annual cycle of oceanographic variation in the CMP and the variables in this study. The lag observed in this study was not only found in SST but also in the cross-correlation analysis, in which we observed that with respect to the other variables, the effect of the 2015-2016 El Niño event was also detected months after the highest values of this MEI were present.

Unlike the other physicochemical variables, PIC showed no differences between the cold and warm seasons. However, PIC exhibited 2 periods of increase based on monthly measurements. The first period occurred from April to May and was related to the influence of the CC and the resuspension of seafloor material by upwelling events in the coastal zone. The second period, which was stronger than the first, was observed from August to October and was a consequence of the intense rainy season that resulted in the runoff of terrigenous materials via different rivers (the Ameca River, Grande de Santiago River, and Armeria River) and brackish water lagoons (Marismas Nacionales). During this period, PIC increased primarily during September (Figs. 1e, 2e).

Overall, the cold seasons observed in 2010 and 2011 in the CMP were related to a La Niña event. As such, this cold season coincided with the most abrupt changes in the other variables considered in this study. The interannual variation in SST during this period reflected a change of 6 °C in only 6 months, which was greater than the temperature difference observed over the same period in the other years considered in this study (i.e., 4 °C). While low La Niña MEI values were recorded in September 2010, the effects of this phenomena were observed 4 months later (January-February) in the CMP. This apparent lack of relationship between the MEI and the more extreme values can be attributed to the time it takes severe La Niña and El Niño events to influence the region. During 2015-2016, an increase in temperature was observed due to the strong and extended El Niño event, which further intensified the characteristics of the warm-water period. The effects of this event were observed after 6 months and were characterized by the lowest PP, Chl, K490, POC, and PIC values.

The intraannual variability in this study was consistent with that of previous studies of the GC. However, these studies also acknowledge that the variability of the CMP is particularly high compared to that of the rest of the gulf (Mercado-Santana et al. 2017). The PP values in Mexican seas indicate that a highly productive period occurs within the area covered by this study during the cold season, which is only surpassed by the high productivity of the northern GC (Sosa-Ávalos et al. 2017). Additionally, the CMP shows great environmental variability between the cold and warm periods throughout the year (Mercado-Santana et al. 2017, Sosa-Ávalos et al. 2017). This study shows that the main cycle of variation for the 6 oceanographic variables evaluated corresponds to 12 months, which further supports the importance and influence of currents on the CMP and its biological resources. Indeed, the CMP is considered a temporary convergence zone for 2 marine currents. The natural oceanographic variability of the CMP combined with the increasingly intense and long-lasting effects of ENSO events augment the existing oceanographic variation within the CMP region.

Acknowledgments

CVPS received a PhD Scholarship from the Consejo Nacional de Ciencia y Tecnología (CONACYT, 388102/252832). The present research was supported by the projects PROSNI2018/P3E/243453 and 249708 supervised by ALCM.

References

Álvarez-Borrego S, Lara-Lara JR. 1991. The physical environment and primary productivity of the Gulf of California. In: Dauphin JP, Simonet BRT (eds.), The Gulf and Peninsular Province of the Californias. Tulsa (OK): AAPG. p. 555-567. http://doi.org/10.1306/M47542C26 [ Links ]

Alvarez-Borrego SA, Schwartzlose RA. 1979. Masas de agua del Golfo de California = Water masses of the Gulf of California. Cienc Mar. 6(1-2):43-63. https://doi.org/10.7773/cm.v6i1.350 [ Links ]

Anderson DM, Glibert PM, Burkholder JM. 2002. Harmful algal blooms and eutrophication: Nutrient sources, composition, and consequences. Estuaries. 25:704-726. https://doi.org/10.1007/BF02804901 [ Links ]

Audzijonyte A, Richards SA, Stuart-Smith RD, Pecl G, Edgar GJ, Barrett NS, Payne N, Blanchard JL. 2020. Fish body sizes change with temperature but not all species shrink with warming. Nat Ecol Evol. 4:809-814. https://doi.org/10.1038/s41559-020-1171-0 [ Links ]

Badan-Dangon A. 1998. Coastal circulation from the Galápagos to the Gulf of California. In: Robinson AR, Brink, KH (eds.), The Sea, The Global Coastal Ocean Regional Studies and Syntheses, vol. 11. Cambridge (MA): John Wiley and Sons. p. 315-334. [ Links ]

Ballance LT, Pitman RL, Fiedler PC. 2006. Oceanographic influences on seabirds and cetaceans of the eastern tropical Pacific: A review. Prog Oceanogr. 69(2-4):360-390. https://doi.org/10.1016/j.pocean.2006.03.013 [ Links ]

Brander K. 2010. Impacts of climate change on fisheries. J Mar Syst. 79(3-4):389-402. https://doi.org/10.1016/j.jmarsys.2008.12.015 [ Links ]

Brodeur RD, Pearcy WG, Ralston S. 2003. Abundance and distribution patterns of nekton and micronekton in the Northern California Current transition zone. J Oceanogr. 59:515-535.https://doi.org/10.1023/A:1025548801541 [ Links ]

Carpenter KE, Abrar M, Aeby G, Aronson RB, Banks S, Bruckner A, Chiriboga A, Cortés J, Delbeek JC, DeVantier L, et al. 2008. One-third of reef-building corals face elevated extinction risk from climate change and local impacts. Science. 321(5888):560-563. https://doi.org/10.1126/science.1159196 [ Links ]

Cepeda-Morales J, Beier E, Gaxiola-Castro G, Lavín MF, Godínez VM. 2009. Effect of the oxygen minimum zone on the second chlorophyll maximum in the Eastern Tropical Pacific off Mexico = Efecto de la zona del mínimo de oxígeno en el segundo máximo de clorofila en el Pacífico Oriental Tropical Mexicano. Cienc Mar. 35(4):389-403. https://doi.org/10.7773/cm.v35i4.1622 [ Links ]

Cepeda-Morales J, Gaxiola-Castro G, Beier E, Godínez VM. 2013. The mechanisms involved in defining the northern boundary of the shallow oxygen minimum zone in the eastern tropical Pacific Ocean off Mexico. Deep Sea Res Part I. 76:1-12. https://doi.org/10.1016/j.dsr.2013.02.004 [ Links ]

Cepeda-Morales J, Hernández-Vásquez F, Rivera-Caicedo J, Romero-Bañuelos C, Inda-Díaz E, Hernández-Almeida O. 2017. Seasonal variability of satellite derived chlorophyll and sea surface temperatura on the continental shelf of Nayarit, Mexico = Variabilidad estacional de clorofila y temperatura superficial del mar satelital en la plataforma continental de Nayarit, México. BioCien. 4(6):1-17. https://doi.org/10.15741/revbio.04.06.07 [ Links ]

Collins CA, Garfield N, Mascarenhas ASJ, Spearman MG, Rago TA. 1997. Ocean current across the entrance to the Gulf of California. J Geophys Res. 102(C9):20927-20936. https://doi.org/10.1029/97JC01302 [ Links ]

Cortés-Lara M, Cortés-Altamirano R, Cupul-Magaña AL, Rodríguez-Nava LV, Vega-Villasante F. 2012. Guía de Florecimientos Microalgales (2000-2011) Causantes de Mareas Rojas en la Bahía de Banderas Jalisco-Nayarit. Guadalajara (Mexico): Universidad de Guadalajara. 104 p. [ Links ]

Di Lorenzo E, Cobb KM, Furtado JC, Schneider N, Anderson BT, Bracco A, Alexander MA, Vimont DJ. 2010. Central Pacific El Niño and decadal climate change in the North Pacific Ocean. Nat Geosci. 3:762-765. https://doi.org/10.1038/ngeo984 [ Links ]

Donelson JM, Munday PL, McCormick MI, Pankhurst NW, Pankhurst PM. 2010. Effects of elevated water temperature and food availability on the reproductive performance of a coral reef fish. Mar Ecol Prog Ser. 401:233-243. https://doi.org/10.3354/meps08366 [ Links ]

Espinosa-Carreón LE, Valdez-Holguín JE. 2007. Variabilidad inter-anual de clorofila en el golfo de California = Gulf of California interannual chlorophyll variability. Ecol Appl. 6:83-92. [ Links ]

Escalante F, Valdez-Holguín JE, Álvarez-Borrego S, Lara-Lara JR. 2013. Temporal and spatial variation of sea surface temperature, chlorophyll a, and primary productivity in the Gulf of California = Variación temporal y espacial de temperatura superficial del mar, clorofila a y productividad primaria en el golfo de California. Cienc Mar . 39(2):203-215. https://doi.org/10.7773/cm.v39i2.2233 [ Links ]

Fiedler PC, Talley LD. 2006. Hydrography of the eastern tropical Pacific: A review. Prog Oceanogr . 69(2-4):143-180. https://doi.org/10.1016/j.pocean.2006.03.008 [ Links ]

Filonov AE, Tereshchenko IE, Monzón CO, González-Ruelas ME, Godínez-Domínguez E. 2000. Variabilidad estacional de los campos de temperatura y salinidad de la zona costera de los estados de Jalisco y Colima, México = Seasonal variability of the temperature and salinity fields in the coastal zone of the states of Jalisco and Colima, Mexico. Cienc Mar . 26(2):303-321. http://doi.org/10.7773/cm.v26i2.577 [ Links ]

Gafar NA, Eyre BD, Schulz KG. 2019. Particulate inorganic to organic carbon production as a predictor for coccolithophorid sensitivity to ongoing ocean acidification. Limnol Ocean Lett. 4(3):62-70. https://doi.org/10.1002/lol2.10105 [ Links ]

García-Oliva F, Camou A, Maass M. 2002. El clima de la región central de la costa del Pacífico mexicano. In: Noguera FA, Vega JH, García-Aldrete AN, Quesada M (eds.), Historia Natural del Bosque Caducifolio de Chamela. Mexico: Instituto de Biología de la UNAM. p. 3-10. [ Links ]

García-Oliva F, Ezcurra E, Galicia L. 1991. Pattern of rainfall distribution in the central Pacific coast of Mexico. Geogr Ann Ser A. 73(3-4):179-186. https://doi.org/10.1080/04353676.1991.11880343 [ Links ]

[IPCC] Intergovernmental Panel on Climate Change. 2002. Climate change and biodiversity, Technical Paper V. [place unknown]: IPCC. 86 p; accessed 2020 Jul 17. https://wedocs.unep.org/20.500.11822/32705. [ Links ]

Green BS, Fisher R. 2004. Temperature influences swimming speed, growth and larval duration in coral reef fish larvae. J Exp Mar Biol Ecol. 299:115-132. https://doi.org/10.1016/j.jembe.2003.09.001 [ Links ]

Harley CDG, Randall Hughes A, Hultgren KM, Miner BG, Sorte CJB, Thornber CS, Rodriguez LF, Tomanek L, Williams SL. 2006. The impacts of climate change in coastal marine systems. Ecol Lett. 9(2):228-241. https://doi.org/10.1111/j.1461-0248.2005.00871.x [ Links ]

Hopkins J, Henson SA, Poulton AJ, Balch WM. 2019. Regional characteristics of the temporal variability in the global particulate inorganic carbon inventory. Glob Biogeo Cycles. 33(11):1328-1338. https://doi.org/10.1029/2019GB006300 [ Links ]

Kumar A, Hoerling MP. 2003. The nature and causes for the delayed atmospheric response to El Niño. J Clim. 16(9):1391-1403. https://doi.org/10.1175/1520-0442-16.9.1391 [ Links ]

Lavín MF, Beier E, Gómez-Valdés J, Godínez VM, García J. 2006. On the summer poleward coastal current off SW México. Geophys Res Lett. 33(2):2-5. https://doi.org/10.1029/2005GL024686 [ Links ]

Lavín MF, Marinone SG. 2003. An overview of the physical oceanography of the Gulf of California. In: Velasco Fuentes OU, Sheinbaum J, Ochoa J (eds.), Nonlinear Processes in Geophysical Fluid Dynamics. Dordrecht (Netherlands): Springer. 378 p. https://doi.org/10.1007/978-94-010-0074-1_11 [ Links ]

Ledesma-Vázquez J, Johnson ME, González- Yajimovich O, Santamaría-del-Angel E. 2009. Gulf of California Geography, Geological Origins, Oceanography, and Sedimentation Patterns. In: Johnson ME, Ledesma-Vázquez J (eds.), Atlas of Coastal Ecosystems in the Gulf of California: Tracking Limestone Deposits on The Margin of a Young Sea. Tucson (AZ): The University of Arizona Press. p. 1-10. [ Links ]

Lee T, McPhaden MJ. 2010. Increasing intensity of El Niño in the central-equatorial Pacific. Geophys Res Lett . 37(14):1-5. https://doi.org/10.1029/2010GL044007 [ Links ]

López-Sandoval DC, Lara-Lara JR, Álvarez-Borrego S. 2009a. Phytoplankton production by remote sensing in the region off Cabo Corrientes, Mexico = Producción fitoplanctónica por sensores remotos en la región frente a Cabo Corrientes, México. Hidrobiol. 19:185-192. [ Links ]

López-Sandoval DC, Lara-Lara JR, Lavín MF, Álvarez-Borrego S, Gaxiola-Castro G. 2009b. Primary productivity in the eastern tropical Pacific off Cabo Corrientes, Mexico = Productividad primaria en el Pacífico oriental tropical adyacente a Cabo Corrientes, México. Cienc Mar . 35(2):169-182. https://doi.org/10.7773/cm.v35i2.1530 [ Links ]

Magaña-Rueda V, Pérez JL, Conde C. 1997. El fenómeno de El Niño y la Oscilación del Sur y sus impactos en México. Ciencias. 51:14-18. [ Links ]

Maul GA, Williams F, Roffer M, Sousa FM. 1984. Remotely sensed oceanographic patterns and variability of bluefin tuna catch in the Gulf of Mexico. Oceanol Acta. 7(4):469-479. [ Links ]

McPhaden MJ, Zhang D. 2002. Slowdown of the meridional overturning circulation in the upper Pacific Ocean. Nature. 415:603-608. https://doi.org/10.1038/415603a [ Links ]

Mercado-Santana JA, Santamaría-del-Ángel E, González-Silvera A, Sánchez-Velasco L, Gracia-Escobar MF, Millán-Núñez R, Torres-Navarrete C. 2017. Productivity in the Gulf of California large marine ecosystem. Environ Dev. 22:18-29. https://doi.org/10.1016/j.envdev.2017.01.003 [ Links ]

Moncayo-Estrada R, Castro-Aguirre JL, De La-Cruz Agüero J. 2006. Lista sistemática de la ictiofauna de Bahía de Banderas. México. Rev Mex Biodivers. 77:67-80. [ Links ]

Páez-Osuna F, Sanchez-Cabeza JA, Ruiz-Fernández AC, Alonso-Rodríguez R, Piñón-Gimate A, Cardoso-Mohedano JG, Flores-Verdugo FJ, Carballo JL, Cisneros-Mata MA, Álvarez-Borrego S. 2016. Environmental status of the Gulf of California: A review of responses to climate change and climate variability. Earth-Sci Rev. 162:253-268. https://doi.org/10.1016/j.earscirev.2016.09.015 [ Links ]

Palacios DM, Bograd SJ, Foley DG, Schwing FB. 2006. Oceanographic characteristics of biological hot spots in the North Pacific: A remote sensing perspective. Deep Sea Res Part II. 53(3-4):250-269. https://doi.org/10.1016/j.dsr2.2006.03.004 [ Links ]

Pantoja DA, Marinone SG, Parés-Sierra A, Gómez-Valdivia F. 2012. Numerical modeling of seasonal and mesoscale hydrography and circulation in the Mexican Central Pacific = Modelación numérica de la hidrografía y circulación estacional y de mesoescala en el Pacífico central mexicano. Cienc Mar . 38(2):363-379. https://doi.org/10.7773/cm.v38i2.2007 [ Links ]

Reyes-Bonilla H. 2003. Coral reefs of the Pacific coast of México. In: Cortés J (ed.), Latin American Coral Reefs. Amsterdam (Netherlands): Elsevier Science B.V. p. 331-349. https://doi.org/10.1016/B978-044451388-5/50015-1 [ Links ]

Reyes-Coca S, Troncoso-Gaytán R. 2001. “El Niño Oscilación del Sur” y los fenómenos hidrometeorológicos en Baja California: el evento de 1997-1998. INP. SAGARPA. México. Ciencia Pesquera. 15:89-96. [ Links ]

Ribas-Ribas M, Hernández-Ayón JM, Camacho-Ibar VF, Cabello-Pasini A, Mejia-Trejo A, Durazo R, Galindo-Bect S, Souza AJ, Forja JM, Siqueiros-Valencia A. 2011. Effects of upwelling, tides and biological processes on the inorganic carbon system of a coastal lagoon in Baja California. Estuar Coast Shelf Sci. 95(4):367-376. https://doi.org/10.1016/j.ecss.2011.09.017 [ Links ]

R Core Team. 2021. R: A language and environment for statistical computing. Vienna (Austria): R Foundation for Statistical Computing; accessed 08 December 2021. https://www.R-project.org/. [ Links ]

Robertson DR, Cramer KL. 2009. Shore fishes and biogeographic subdivisions of the Tropical Eastern Pacific. Mar Ecol Prog Ser . 380:1-17. https://doi.org/10.3354/meps07925 [ Links ]

Rodríguez-Troncoso AP, Carpizo-Ituarte E, Cupul-Magaña AL. 2016. Physiological response to high temperature in the Tropical Eastern Pacific coralPocillopora verrucosa. Mar Ecol. 37:1168-1175. https://doi.org/10.1111/maec.12392 [ Links ]

Simons RA, Mendelssohn R. 2012. ERDDAP-A Brokering Data Server for Gridded and Tabular Datasets. In: [AGU] American Geophysical Union (ed.), Fall Meeting Abstracts Vol. 2012. San Francisco (CA): American Geophysical Union. p. IN21B-1473; accessed 2018 May. https://coastwatch.pfeg.noaa.gov/erddap/index. [ Links ]

Sosa-Ávalos R, Espinosa-Carreón TL, Santamaría-del-Ángel E, Cepeda-Morales J, Mercado-Santana A, De la Cruz-Orozco ME, Lara-Lara JR, Álvarez-Borrego S, Gaxiola-Castro G. 2017. Diagnóstico y análisis de la producción primaria neta en los mares de México. Elem Pol Púb. 2(3):164-185. [ Links ]

Spalding MD, Fox HE, Allen GR, Davidson N, Ferdaña ZA, Finlayson M, Halpern BS, Jorge MA, Lombana A, Lourie SA, et al. 2007. Marine ecoregions of the world: a bioregionalization of coastal and shelf areas. BioScience. 57(7):573-583. https://doi.org/10.1641/B570707 [ Links ]

Speers AE, Besedin EY, Palardy JE, Moore C. 2016. Impacts of climate change and ocean acidification on coral reef fisheries: An integrated ecological-economic model. Ecol Econ. 128:33-43. https://doi.org/10.1016/j.ecolecon.2016.04.012 [ Links ]

Wild C, Hoegh-Guldberg O, Naumann MS, Colombo-Pallotta MF, Ateweberhan M, Fitt WK, Iglesias-Prieto R, Palmer C, Bythell JC, Ortiz J, et al. 2011. Climate change impedes reef ecosystem engineers. Mar Freshw Res. 62(2):205-215. https://doi.org/10.1071/MF10254 [ Links ]

Wolter K, Timlin MS. 2011. El Niño/Southern Oscillation behaviour since 1871 as diagnosed in an extended multivariate ENSO index. Int J Clim atol. 31(7):1074-1087. https://doi.org/10.1002/joc.2336 [ Links ]

Zamudio L, Leonardi AP, Meyers SD, O’Brien JJ. 2001. ENSO and eddies on the southwest coast of Mexico. Geophys Res Lett . 28(1):13-16. https://doi.org/10.1029/2000GL011814 [ Links ]

Zeitzschel B. 1969. Primary productivity in the Gulf of California. Mar Biol. 3:201-207. https://doi.org/10.1007/BF00360952 [ Links ]

Received: March 04, 2021; Accepted: May 27, 2022; Published: March 29, 2023

*Corresponding author: E-mail: amilcar.cupul@gmail.com

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License