Next Article in Journal
Research on Torque Performance of Marine Hybrid Excitation Synchronous Motors Based on PSO Optimization of Magnetic Permeability Structure
Previous Article in Journal
Aerodynamic Performance and Coupling Gain Effect of Archimedes Spiral Wind Turbine Array
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

El Niño–Southern Oscillation Diversity: Effect on Upwelling Center Intensity and Its Biological Response

by
Gabriel Santiago Gutiérrez-Cárdenas
1,
Enrique Morales-Acuña
2,*,
Leonardo Tenorio-Fernández
1,
Jaime Gómez-Gutiérrez
1,
Rafael Cervantes-Duarte
1,* and
Sergio Aguíñiga-García
1
1
Instituto Politécnico Nacional, Centro Interdisciplinario de Ciencias Marinas, La Paz 23096, Baja California Sur, Mexico
2
Instituto Politécnico Nacional, Centro Interdisciplinario de Investigación para el Desarrollo Integral Regional, Bulevar Juan de Dios Bátiz Paredes, Col. San Joachin, Guasave 81049, Sinaloa, Mexico
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(7), 1061; https://doi.org/10.3390/jmse12071061
Submission received: 11 May 2024 / Revised: 15 June 2024 / Accepted: 19 June 2024 / Published: 24 June 2024
(This article belongs to the Section Physical Oceanography)

Abstract

:
The influence of El Niño–Southern Oscillation (ENSO) on productivity at upwelling systems has been extensively studied. However, in recent decades, ENSO diversity has been documented; there are more frequent events with anomalies in the central Pacific (CP) causing ecological impacts that are different from the canonical events centered in the eastern Pacific (EP). We hypothesize that ENSO effects on upwelling intensity and its biological response are highly dependent on ENSO diversity. Wavelet coherence was computed for monthly standardized anomalies of Ekman transport and sea surface Chlorophyll-a (Chl-a) concentration from eight upwelling centers located along the west coast of the Baja California peninsula (1993–2022). Ekman transport (ET) showed a non-linear association with ENSO at interannual periodicities between 1.2 and 7 years, as well as semiannual scale variability. Coherence between ENSO and ET showed a reduction in upwelling intensity during CP ENSO events and an increased upwelling intensity after EP El Niño events with regional differences. Our results reveal an influence of northern Pacific circulation that subsequently influences ENSO events during its development and its peak. We detected an influence of SST and wind anomalies on the ENSO–Ekman transport connection associated with the northern Pacific Meridional Mode pattern. The CP shows an important role in sea surface Chl-a concentration anomalies (1998–2022). We highlight the conceptual relevance of distinguishing ENSO diversity, with influential ecological effects associated with CP ENSO events.

1. Introduction

In wind-driven upwelling systems, economic activities influenced by climate variability typically develop in the most productive coastal regions of the world’s ocean, contributing to over 20% of fish catch around the world and 7% of primary productivity [1,2]. In recent decades, the scientific community has focused its efforts on understanding the genesis and processes that give rise to these climatic-oceanographic events. Although the wind drag is the most relevant mechanism of physical force for upwelling genesis, coastline irregularities, stratification, thermocline depth, bathygraphy, and coastal topography, all have a great influence on the dynamics of the upwelling and influence on related physical, biogeochemical [3,4], hydrological [5], and biological processes [6]. As a result of the irregularities of the coastline and wind fluctuations, these factors are not uniform along the coast, causing the formation of upwelling centers [7], characterized by minimum surface temperatures [8,9,10], and high primary productivity [11,12].
The coastal upwelling system along the west of the Baja California peninsula, Mexico, located in the southern California Current System, is well known [13,14,15,16,17]. The influence of El Niño–Southern Oscillation (ENSO) is the main interannual timescale variability factor modulating the physical and biogeochemical environment. ENSO’s impacts on particular events have been inferred mostly from sea surface temperature (SST) [18], salinity [13,19], upwelling intensity [15,20], sea surface Chl-a concentration [21,22,23], primary productivity [24,25], and phytoplankton, mostly focusing on data from the El Niño 1997–1998 [13,20,21,24,26].
One of the main mechanisms that trigger environmental oceanographic changes along the west coast of the Baja California peninsula is atmospheric changes in the Aleutian Low (primary expansion) and northern Pacific High (contraction) during El Niño [27,28,29], and the opposite during La Niña events, both of which drive changes in the wind field [30]. In addition, environmental oceanographic changes in the region are influenced by the oceanographic modulation of thermocline anomalies that propagate eastward from the tropical Pacific and then poleward to the northern Pacific [27,31]. Generally, both mechanisms tend to inhibit (rise) the upwelling of nutrients to the euphotic zone during El Niño (La Niña) and trigger reduced (increased) primary productivity with an impact on various trophic levels. However, the relative influence of each mechanism (atmospheric changes and environmental oceanographic changes) is not fully understood yet, and often two seemingly similar events show differential biological impacts. For example, in the 1997–1998 El Niño event, the decrease in upwelling and productivity was attributed to changes in coastal SST [20]; in contrast, in the 2009–2010 El Niño event, southern water advection did not occur, and the atmospheric teleconnection associated with wind stress was widely recognized as a modulating factor [19].
Two types of ENSO events were recently well-identified according to the position of the maximum SST anomalies, distinguishing between eastern Pacific (EP) and central Pacific (CP) events [32,33,34,35]. This ENSO distinction is a key factor in ocean–atmospheric teleconnections that highlights the traditional use of SST anomalies that limit our understanding of the impact of ocean–atmospheric teleconnections [36,37,38,39]; thus, the distinction between EP and CP events is necessary to better comprehend ENSO teleconnections.
The CP and EP ENSO events have shown opposite and non-linear effects on the frequency and intensity of tropical cyclones [40], but also on South American and North American precipitation [41,42,43,44,45], Eurasian climate, North Atlantic Oscillation [38,39,46,47], western Pacific subtropical high [48], and on the biomass of the oceanic phytoplankton [49], with changes during their development or before their peak. The effects of CP and EP ENSO could be enhanced by global climatic change [50]. Additionally, both types of events do not show the same development and teleconnection patterns; CP variability has a strong connection to subtropics, and in some cases, CP El Niño events are preceded by northeastern Pacific wind and SST anomalies [51,52,53].
The extratropical ENSO–atmosphere–ocean teleconnection affects ENSO diversity [33,36,52]. CP El Niño events have been reported more frequent since the 1990s [54,55,56], with a distinct regional impact on nutrient, primary productivity, and phytoplankton biomass in the northern EP [49]. It is thus expected that ENSO diversity will have different effects on the intensity of upwelling centers and phytoplankton biomass along the west coast of the Baja California peninsula, a region located in the southern sector of the California Current System with profitable economic, biological, and environmental relevance.
We tested two hypotheses: (1) Ekman transport at upwelling centers will be reduced in EP El Niño and will be increased during both CP El Niño and La Niña events, with no substantial effect from EP La Niña; and (2) a substantial sea surface Chl-a reduction response is expected during the development of EP El Niño events, and similar effects during CP ENSO events. To address both hypotheses, we analyzed the diversity of ENSO effects on upwelling centers, using Ekman transport as a measure of upwelling intensity. We also analyzed the biological response of upwelling centers, using the concentration of sea surface Chl-a to distinguish between CP and EP ENSO events. The present study aims to establish the baseline of coastal upwelling centers and ENSO diversity on ocean and atmospheric regional scales along the west coast of the Baja California peninsula, Mexico.

2. Materials and Methods

2.1. Study Area

The west coast of the Baja California peninsula is in the southern California Current System (22–32° N, and 109.5–117° W), a transitional region in the northeastern Pacific characterized by the influence of three sea surface water masses: Subarctic Water (California Current), Subtropical Surface Water, and Tropical Surface Water. This region has a pronounced seasonal cycle that transports the northern coolest water along the west coast of the Baja California peninsula during the winter-spring seasons and the tropical and subtropical seawater masses during the summer–fall seasons [13,14,57]. Oceanographic variability separates the area into two regions divided by Punta Eugenia (~27.5° N); the northern Baja California region has subarctic characteristics, while the southern Baja California has warmer tropical/subtropical conditions. These differences result in a spatiotemporal complex transition region [13,18,21] with intense upwelling activity [20], and complex patterns of filaments, meanders, and mesoscale structures [58] that carry nutrients from coastal waters to offshore regions driven by a high variable seasonal and interannual variability in physical and pelagic biological characteristics [59].
The seasonal wind field along the west coast of the Baja California peninsula is controlled by the position and intensity of the northern Pacific High cell pressure [60,61], triggering prevailing northwest winds for the entire region with maximum intensity during spring and summer, and minimum intensity during autumn and winter [20]. Due to northwest winds, Coriolis force, and coastline orientation, upwelling is the primary nutrient input to the euphotic zone with a seasonal cycle that follows the changes in the wind stress field [20]. There are previously identified upwelling regions along the west coast of the Baja California peninsula, two in the north (Ensenada and northern Punta Eugenia) and two in the south (Punta Eugenia and Bahía Magdalena) [20]. Figure 1 shows the characteristics of the study area and the location of the eight upwelling centers recorded so far.

2.2. ENSO Diversity

ENSO events were identified using three distinct indices to represent ENSO diversity, and also the differences between CP and EP ENSO events for El Niño and La Niña phases. The indices were based on the SST of Niño 3 and Niño 4 regions (map details shown in: https://www.cpc.ncep.noaa.gov; accessed on 30 April 2024) to capture inherent differences in the maximum (minimum) longitudinal position anomaly of ENSO, as a key factor in the global effects of ENSO diversity [39,45]. Specifically, EP El Niño (EP La Niña) events are characterized by the warmest (coldest) anomalies located in the Niño 3 region; on the other hand, CP El Niño (CP La Niña) events are characterized by the maximum (minimum) anomalies located in the Niño 4 region [38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62].
We used the monthly Extended Reconstructed SST (ERSST V.5), which was calculated using 1991–2020 climatology data for each region and was normalized using the following equation [63]:
Z i = ( S S T A i X ¯ ) S D
where Z i and S S T A i are the standardized anomaly and the SST anomaly for each month i, respectively. X ¯ and SD are the mean and the standard deviation of the SST anomalies. The EP and CP SST anomalies indices were calculated following these equations [62]:
E P = Z N i ñ o 3 0.5 × Z N i ñ o 4
C P = Z N i ñ o 4 0.5 × Z N i ñ o 3
where Z N i ñ o 3 and Z N i ñ o 4 represent the normalized SST anomalies of Niño 3 (5° S–5° N, 150–90° W) and Niño 4 (5° S–5° N, 160–150° W) regions, respectively [62]. Figure 2 shows the EP ENSO and CP ENSO indices, comparing indices (black line) with Niño 3.4 normalized SST anomalies (colors). Table 1 shows EP and CP El Niño/La Niña events that are higher than one SD. When the EP or CP index is higher than one SD for December–February, the event is classified as an EP or CP ENSO, respectively. The ENSO composites for each type of ENSO are shown in Figure A1.

2.3. Wind Data and Ekman Transport

We used the Cross Calibrated Multiplatform 10 m wind monthly data from 1993 to 2022 from the NASA Physical Oceanography Distributed Archive Center. The data were calibrated using Variational Analysis Methods to produce a complete grid with high-resolution (0.25° × 0.25°) zonal wind (u) and vector wind (v) data sets. The data fit in a vector difference of 0.8 m s−1 to QuikSCAT with no wind direction bias and 0.6 m s−1 to high-quality TAO buoys. This data set is capable of revealing ocean surface wind field features that have been described in previous studies and provides high-spatial-resolution time averages of ocean surface winds [64]. This data set is available at https://www.remss.com/measurements/ccmp/ (accessed on 30 September 2023).
Eight upwelling centers on the west coast of the Baja California peninsula were selected to calculate alongshore wind stress and Ekman transport as a measure of upwelling intensity for each upwelling center. For calculations, we used a standardized method [65] that allows for (1) a systematic approach for Ekman transport estimation via the simplification of a high-resolution coastline without losing accurate determination of coastal angles, and (2) a higher Ekman transport spatial resolution calculation, allowing for the highest precision available for each coastal upwelling center.
The coastal angle estimation and Ekman transport method were applied separately for each upwelling center section. For this purpose, a Douglas–Peuker simplification method was used with ε = 0.25 by matching the wind spatial resolution. The coastal angle was estimated perpendicularly to the approximate coastline pointing offshore. Once the coastal angle was estimated, the wind stress was computed as follows:
τ x y = l a t l a t τ x c o s θ 90 ° + τ y s i n θ 90 ° ,
where τ x y is the alongshore wind stress, θ is the coastal angle, and l a t represents latitude in the [−90, 90] interval; all parameters were computed for each upwelling center. The Ekman transport (ET) is calculated as follows:
E T = τ x y ρ f ,
where ρ is the average density of seawater (1025 kg m−3) and f is the Coriolis parameter computed using the latitude of each upwelling center. These parameters were computed for each upwelling center shown in Figure 1.

2.4. Sea Surface Chlorophyll-a Concentration

The sea surface Chl-a concentration was obtained from the Copernicus Climate Change Service (C3S) v. 5.0.1 (doi:10.24381/cds.f85b319d), which contains daily composites of sensors including VIIRS, OLCI, MODIS Aqua, SeaWiFS, and MERIS. A more complete spatial coverage is thus obtained than when data from independent sensors are used. Data were taken for the period 1998–2022 with a 0.0416° × 0.0416° spatial resolution (data are available at https://cds.climate.copernicus.eu; accessed on 30 May 2024). Monthly time series of sea surface Chl-a concentration in mg m−3 were extracted for the location of each of the eight upwelling centers.

2.5. Wavelet Transform

The wavelet transformation analysis allows the estimation of the phase and amplitude of each component of a time series simultaneously [66], and follows the evolution of the distinct form of each signal. The modern wavelet analysis uses a convolution of the real-time series x t with multiple functions Ψ ( t : τ ,   a ) derived from a complex “mother” wavelet Ψ ( t ) by scaling by a , and converting time in τ :
T x = 1 a Ψ * a 1 t τ x t d t .
The complex conjugate is represented by the symbol *. For the Ψ t selection, it will have a similar pattern and x t characteristics. Higher T x values denote that x t follows the desired form. For the current study, the Morlet wavelet is selected as follows:
Ψ t = π 1 / 4 e t 2 / 2 e i ω t .
This mother wavelet is selected according to (1) improved localization and detection of scale and (2) better frequency resolution. Additionally, the Morlet wavelet splits the wavelet into the part that describes the characteristics of an oscillatory times series and another part that retains the phase information required for coherence with another time series [67], which is the main purpose of the present study.

2.6. Wavelet Coherence

The wavelet-based coherence is the method selected to compare ENSO indices with the Ekman transport and sea surface Chl-a concentration time series. The wavelet-based coherence evaluates the frequency and phase synchronization of signals [67,68], which in turn is based on cross-wavelet analysis concepts. The cross-wavelet transform was calculated according to the following equation [69]:
T x , y τ , a = a 1   T x τ , a ·   T y * τ , a ,
where Tx (τ, s) and Ty (τ, s) follow Equation (6)’s form. To compare both time series the energy density was calculated producing the cross-wavelet power spectrum that is useful to compare both time series:
P x , y τ , a = | T x , y τ , a | .
Therefore, the cross-wavelet power spectrum depends on the units of the time series. The wavelet coherency that avoids an erroneous interpretation of the results is defined as follows:
C x t , y t = a T x , y a a 1 T x τ , a 2 · a a 1 T y τ , a 2 .
The wavelet coherence is defined as follows:
C x , y 2 ( x , y ) = | a ( a 1 T x τ , a · T y * τ , a ) | 2 a a 1 T x τ , a 2 · a a 1 T y τ , a 2 .
In Equation (11), C x , y 2 ( x , y ) varies between 0 and 1, where 0 indicates no coherence between both time series, and 1 indicates that the coherence between both time series is the maximum [67].
To obtain information of phase and synchronization in terms of phase, the phase difference was calculated, indicating which signals leads the other [66,70]:
ϕ x τ , a ϕ y τ , a = A r g T x , y τ , a = tan 1 ( I m ( T x , y τ , a ) R e ( T x , y τ , a ) ) ,
where ϕ x ( τ , a ) and ϕ y ( τ , a ) are phases of each individual signal. The series moves in-phase (anti-phase) if the phase difference defined in Equation (12) has an absolute value below (greater than) π / 2 . The leading time series is defined by the sign of the difference between the phases. Both wavelet coherence and phase difference are displayed in the spectrum; C x , y 2 ( x , y ) is represented by colors and ϕ x τ , a ϕ y τ , a by arrows [68]. The phase difference (Equation (12)) varies from π to π . For a given period (P) (for example, P = 2 years), a difference of phase of π / 2 between x and y time series is equal to 0.5 years lag. Where the ENSO index shows a decrease and negative value (that represents La Niña), and the phase difference shows in-phase, a decrease in Ekman transport index occurs; whereas, when La Niña has an anti-phase, Ekman transport increase occurs. The sign represents which time series lead the other. Figure A2 helps to interpret phases and arrows.
The statistical significance of wavelet coherence was tested using the theoretical red noise wavelet power spectrum for the null hypothesis (“no joint periodicity”) through a Monte Carlo simulation using 1024 repeats. The cone of influence, contour lines, and arrows delimit the areas with >95% confidence interval (see [67,68,69,70,71] for details).

3. Results

The selected upwelling centers along the west coast of the Baja California peninsula show different Ekman transport intensities and seasonal cycles (Table 2 and Figure 3). The highest monthly mean upwelling intensity and the highest long-term variability were observed in the Bahía Magdalena region. The boxplot shows that the >75% monthly mean was <−0.4 × 10−8 m2 s−1 for the spring. Punta Eugenia and the north of Punta Eugenia showed similar interannual variability and a median between −0.2 × 10−8 and −0.5 × 10−8 m2 s−1. Ensenada upwelling centers showed the lowest Ekman transport of the study area, with lower interannual variability in comparison with other upwelling centers, mainly during winter. In December, the Ekman transport was below −0.5 × 10−9 m2 s−1 in 75% of the recorded years.
We successfully identified 20 ENSO events using the EP and CP indices (Table 1), achieving effective differentiation between EP and CP events for both El Niño and La Niña events. The composites are shown in Figure A1. The wavelet coherence between ENSO and monthly Ekman transport shows a similar coherence structure for the upwelling centers of the same region; however, there is a difference in coherence and phase between different regions and EP or CP indices (Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8).
In Ensenada, the periods with interannual significant coherence ranged between 1 and 2 years for both ENSO indices (Figure 4). For the EP ENSO, the interannual coherence of ~2 years indicates anti-phase coherence during 1997, 2007–2009, and 2014–2017 years, with a 2- to 5-month lag. The CP ENSO and Ekman transport wavelet coherence showed anti-phase coherence in the 1.2- to 2-year periods for La Niña events of 2000–2001 and 2010–2011, but also in-phase coherence for the 2016–2022 period, both with a lag between 0 and 4.5 months.
Differences in coherence between ENSO and Ekman transport indices in the north of Punta Eugenia are shown in Figure 5. Periods of 3 to 4 years show Ekman transport lead coherence between 1996 and 1998 with a 2-month lag. For 1.5- to 2-year periods, the 2015–2016 El Niño shows anti-phase lead Ekman transport. These observations imply an increase in Ekman transport during this El Niño event. In the case of CP ENSO, anti-phase coherence for 4 to 7 periods shows lead Ekman transport time series with a lag of 9 to 12 months.
Additionally, the coherence of CP and Ekman transport anomalies shows two periods of anti-phase lagging at 1.5 years. The annual periodicity of wavelet coherence, under multi-year CP La Niña conditions for the 1997–2002 period, shows 3 months lead Ekman transport. The 2015–2021 period shows ENSO events with predominantly positive anomalies at CP (Figure 2), and an ENSO lead with ~2–3 months lag upwelling reduction for both upwelling centers. In 2021–2022, CP La Niña events show in-phase coherence.
The coherence in Punta Eugenia shows differences between EP and CP ENSO with Ekman transport (Figure 6A–D). The EP ENSO wavelet coherence shows in-phase coherence in ~2- to 2.6-year periods between 1995 and 2008, and anti-phase coherence for a 2 year period between 2013 and 2016 for both upwelling centers (Figure 6A,B). The CP ENSO and Ekman transport wavelet coherence for Punta Eugenia shows a 4- to 7-year periodicity lead Ekman transport, with a lag around 8 to 14 months. A ~2.5-year periodicity shows a CP ENSO lead between 2009 and 2010 with a lag of around 5 months. Additionally, periods around 1.2 to 1.5 years show a CP index lead coherence between 2016 and 2020, a predominantly positive CP anomaly period with a lag of 3 to 5 months with upwelling reduction (Figure 2).
The coherence spectra of EP ENSO in Bahía Magdalena (located in the southern region of our study site) showed two main periods with interannual periodicity (Figure 7): ~2-year periodicity and the period of 1997–1998 showed an EP ENSO lead during a period with predominantly positive EP anomalies with upwelling intensity reduction and a lag response of around 4 months. The CP index did not show coherence during interannual periods.
The wavelet coherence in Ensenada at annual and semiannual timescales showed in-phase coherence with a lag of 2–3 months, predominantly under the conditions of CP La Niña in the year 2000. In 1995 the coherence showed an EP lead by ~2 months. In contrast, the CP showed an opposite direction of coherence. The same phenomenon was observed in 1998, with anti-phase coherence for EP and in-phase coherence for CP. Similarly, the periods 2007–2008, 2010–2011, 2013–2014, and 2022 showed both EP and CP indices lead. The lag response of the Ekman transport in relation to the ENSO index was between 0 and 2 months (Figure 4). In the north of Punta Eugenia, the coherences with a periodicity below 1 year showed anti-phase and in-phase coherence with a lag of 0 to 2 months. There was also an opposite coherence with EP and CP events for 1995, 1997, 1998, 2004–2006, 2012–2013, and 2022. The EP ENSO index showed a 3-month lag in Ekman transport response in 2003 (Figure 5).
Similarly, in Punta Eugenia, at annual and semiannual periodicities, the coherence between ENSO and Ekman transport showed an opposite direction for EP and CP indices in 1994–1995, 2019–2020, and 2022. The periods 2005–2006 and 2012–2013 showed an EP ENSO lead and Ekman transport leading above CP for the same period (Figure 6). In Bahía Magdalena, a 0.5-year periodicity was observed, with opposite coherence for EP and CP indices in 1998, 2009–2010, 2011–2013, 2015–2016, and 2019–2020. Additionally, the periods 2007–2008 and 2011–2012 showed coherence with the CP index. The lag response was between 0 and 3 months.

4. Discussion

4.1. ENSO and Upwelling Intensity: Interannual Coherence

The present study reveals an upwelling intensity response to ENSO diversity, with differences in ENSO and Ekman transport coherences for eight upwelling centers located along the west coast of the Baja California peninsula. We identified differences between EP and CP ENSO on the Ekman transport response in the time and periodicity domains for the 30-year period of study (Figure 4, Figure 5, Figure 6 and Figure 7). Our results show that not all ENSO events show coherence and that not all El Niño events with an Ekman transport anomaly response show a canonical El Niño and La Niña pattern in upwelling. These canonical patterns have been reported previously, and include a decrease (increase) in upwelling for El Niño (La Niña) [22,23], coupled with changes in primary productivity [24]. These observations imply that other factors that have not been analyzed in detail may affect the diversity of El Niño events and have not been properly conceptualized.
Our results show that the eight upwelling centers in the four main regions of the coast of the Baja California peninsula show similar responses in terms of the direction to ENSO at interannual wavelet coherence periodicity. The Ekman transport response was opposite for EP El Niño and CP El Niño in some cases, with differences in the periods where coherence between ENSO and Ekman transport was detected. In Ensenada, the 1997–1998 and 2015–2016 EP El Niño events were the strongest El Niño events of the last century [72] and the only two EP El Niño of last three decades [38,73]. Both EP El Niño events showed anti-phase coherence which led to an Ekman transport increase 1–4 months later, in the spring before the ENSO peak (~December). This is consistent with observations of upwelling favored by atmospheric forcing during the spring of 1998, attributed to the ending of EP El Niño [22].Our results thus reveal a consistent response of Ekman transport to EP El Niño. The same coherences were observed in upwelling centers in the north of Punta Eugenia (Figure 5A,B)
The EP ENSO index anti-phase coherence of 2007–2009 at 2-year periodicity coincides with the 2007–2008 La Niña event, which was an event with multiple location minima throughout the year, and minimum values of the EP index developing from May–March 2006 and extending into December 2007 [56]. According to our results, the 2007–2008 La Niña event led to a decrease in upwelling with a lag of 1–4 months, with Ekman transport standardized anomalies >1 during November (0 years) and January (1 year). Furthermore, in the spring of 2008, the Ekman transport showed a standardized anomaly of −2, with a higher increase in upwelling during this period. The CP ENSO index, which showed minimum values in February 2008 (Figure 2) led to an increase in upwelling with a lag of 3 months (Figure 4C,D), which suggests an opposite role of CP and EP negative anomalies on Ekman transport. Monthly anomalies in Ensenada in April 2007 reported upwelling with increased intensity, coinciding with changes in thermohaline conditions in surface and subsurface coastal sea waters [74].
In 2001–2002, there was an increase in upwelling during December–February, but it was not reported as an El Niño event in the eastern tropical Pacific. In this period, both CP and EP indices showed opposite coherence. Interestingly, this period was characterized by abrupt oceanic and atmospheric changes with westerly anomalies in December 2000; the SST became warmer in the CP which might have led to an El Niño event by early spring with a different development relative to a canonical El Niño pattern [75]. Both CP and EP indices not only show opposite coherence during this period but also opposite directions with positive CP anomalies and negative EP anomalies (Figure 2), which suggests non-El Niño upwelling anomalies, possibly related to opposite CP and EP indices and a prominent tropical Pacific gradient. Another significant period with a CP ENSO event that led to a 2–5-month lag coherence was in 2016–2021 (Figure 4C,D). In spite of the presence of multi-year La Niña events in 2017 and 2018, the CP index of the period 2016–2020 showed predominantly positive anomalies. During this period, upwelling centers in Ensenada showed multiple months with decreased Ekman transport, suggesting that the CP El Niño triggers a decrease in upwelling, contrasting with the 1997–1998 EP El Niño event.
In the north of Punta Eugenia, coherence showed that the Ekman transport led 4–5 months during 1997–2001 in a ~3.5-year periodicity. Initially, the Ekman transport on the west coast might not have led to an ENSO event. Furthermore, changes in Ekman transport before the peak period of the ENSO event suggest changes in the wind field during or prior to ENSO development, which in turn controls Ekman transport. It is well known that the seasonal and annual wind stress variability along the west coast of the Baja California peninsula is controlled by the Aleutian Low and the northern Pacific High cell pressure [13,60,61]; thus, ENSO events could trigger variations in both atmospheric systems [29,37]. Recent studies demonstrated that the strengthening and weakening of the spring Aleutian Low (developing ENSO year) could lead to Pacific warming in the following year [76]. This is in keeping with the lag of the ENSO–Ekman transport coherence. On the other hand, the North Pacific Oscillation, the second source of variability of the northern Pacific could trigger an ENSO event [36,77]. The same lag of ENSO–Ekman transport coherence was observed in a ~1.5 periodicity, with a 2–3-month lag in 1997, which is the period of EP El Niño development.
Our results highlighted the conceptual relevance of EP and CP indices distinction. A previous study reported a non-significant correlation between the coastal upwelling index and the Multivariate El Niño Index (MEI; r = −0.13) [18]. However, the coherence in our results shows a non-linear Ekman transport response at upwelling centers, which could result in lower linear metrics of correlation analysis. Our results also reveal that the use of a single index for ENSO limits our understanding of their impacts [36,37,38,39]. For instance, in the north of Punta Eugenia, there is an anti-phase lag of the La Niña multi-year event (1997–2003 and 2011–2013) (Figure 5C,D), which suggests an increase in upwelling in response to negative CP anomalies during both La Niña events.
The increase in upwelling was evident after 1997–1998. However, in 2000–2002 there were two CP La Niña events during which the decrease in upwelling was evident. After the El Niño event of 2010, there were multiple years with CP La Niña events, during which the anti-phase coherence was higher, suggesting a decrease in upwelling as a response to negative anomalies. Another period with CP ENSO coherence in the north of Punta Eugenia was 2016–2020 (Figure 5C,D); it was characterized by sustained positive anomalies in the CP (Figure 2). The coherence observed in 2016–2020 suggests an atmospheric northern Pacific response to tropical Pacific CP anomalies with a ~2-month lag and upwelling decreasing in the north of the Baja California peninsula. These differences in Ekman transport response validate the importance of addressing the impact of ENSO diversity on atmospheric impact on the Baja California peninsula.
Both Punta Eugenia and Bahía Magdalena showed an in-phase ~2–4-month lag response in Ekman transport. The maximum coherence was observed in 1996 (Figure 2) when the EP index showed predominantly negative anomalies (Figure 6A,B and Figure 7A,B). In spite of EP La Niña conditions, the conditions of the northern Pacific were unusual, showing positive SST anomalies that extended around the northeast Pacific and persisted well into the 1997–1998 El Niño event showing western tropical-northeastern Pacific positive anomalies and anomalous cyclonic wind anomalies [78]. On the other hand, both Punta Eugenia and Bahía Magdalena regions showed EP and Ekman transport anti-phase coherence from 2013 in the southern region to 2015 or 2016 (Figure 6A,B and Figure 7A,B). In 2013–2016, the oceanographic–atmospheric conditions in the northern Pacific were predominantly associated with “The Blob” (2013–2015) that originated in the north CP, a sustained multiple-year warming event in the northeast that was not driven by either ENSO or a positive North Pacific Meridional Mode (NPMM) pattern. Both the coherence of ENSO and Ekman transport suggest that the initial perturbations that trigger upwelling anomalies started in the north EP.
Both the 1996 and 2015 ENSO events had a similar winter warming of the northern Pacific [53] but with the difference that the predominant anomalies in the tropical Pacific were negative in 1996 and positive in 2015. After October 2013, there was a decrease in upwelling and it remained ET below normal until 2015 [15], similar to what we observed during the La Niña event of 1996 in the present study. The opposite tropical Pacific anomalies and the positive northern Pacific anomalies suggest an Ekman transport anomaly that is not driven by ENSO and is possibly associated with the NPMM and local atmospheric interannual changes in the northern Pacific. There is evidence that the northern Pacific anomalies could be triggered by a wind–evaporation–SST feedback [79], triggered in turn by westward wind anomalies in the western Pacific that are identified as ENSO precursors [53]; these anomalies are related to a third Pacific variability mode [80]. This event is not typically considered to be driven by atmospheric-oceanographic anomalies in local studies carried out along the west coast of the Baja California peninsula, Mexico.
In Punta Eugenia and Bahía Magdalena, the CP ENSO index and Ekman transport showed coherence in 2008–2010. In 2008–2010, in California [19], the importance of atmospheric teleconnection and the lack of tropical Pacific seawater advection to the region was high, which indicates that atmospheric teleconnection is a relevant mechanism for El Niño effects on the California Current System. Like in the northern regions, the CP positive anomaly period of 2016–2020 showed an interannual reduction response in Ekman transport upwelling centers with a 2–5-month lag. At this point, our results highlight the difference in coherence between EP and CP anomalies, and the non-linear teleconnection; we found quasi-opposite coherences between both EP and CP indices for the same anomaly direction at warm events in the tropical Pacific. These findings have great conceptual relevance to distinguish between EP and CP types of ENSO events, and also for tropical Pacific gradient of SSTs, in order to better address the atmospheric impact of ENSO events on the study area.

4.2. Annual and Semiannual Coherence: ENSO and Ekman Transport

The semiannual and annual coherence periods found in the present study could be related to (1) ENSO interaction with other quasi-oscillatory processes, (2) changes in northern Pacific pressure gyres triggered by tropical–extratropical Pacific anomalies, or (3) ENSO precursor processes driven by the main modes of variability (for example, northern Pacific SST [51], Aleutian Low [81] and North Pacific High/NPMM [53,82]). Our results indicate that neither the strongest ENSO events, nor the weaker ones, necessarily trigger SST anomalies in the California Current System that drive changes in atmospheric circulation, wind, and Ekman transport [19,31]. We showed ENSO and Ekman transport wavelet coherence for seven El Niño (La Niña) events along the west coast of the Baja California peninsula. Most of the semiannual coherence estimated between the ENSO and Ekman transport was observed in the spring and summer. In 1996 and 2000, there was a negative CP index and a positive EP index in Ensenada. In contrast, in the ENSO events of 1998, 2007, 2019, and 2022 (Figure 4C,D), coherence was the opposite. Moreover, 1998, 2005, 2008, 2013, and 2021 were opposite in the north of Punta Eugenia (Figure 5C,D), highlighting an opposite response at semiannual timescales to CP and EP events, apparently impacting the tropical Pacific gradient in the Ekman transport response in the study area. Several studies suggest a complex non-linear connection of the atmosphere response to ENSO diversity [27,37] that leads to opposite atmospheric and biological responses [49].
Punta Eugenia and Bahía Magdalena had the same opposite coherence as EP and CP in 1995, 1997, 2009, 2011, 2013, 2015, and 2019; and some coherences show Ekman transport leading to ENSO anomalies, suggesting the existence of changes in the northern Pacific circulation changes that precede the tropical SST variability, which could be captured by the time series variability of both EP and CP. This tropical–extratropical teleconnection was previously documented [36,83], and is in keeping with our results. It should be pointed out that the interaction between ENSO and atmospheric-oceanic circulation systems is not unidirectional; instead, ocean–atmospheric processes occurring in the northern Pacific can also lead to the development of ENSO events [76]. This explains why previous studies found opposite effects of ENSO on upwelling along the coast of California during El Niño events [31]. Our results show that this does not necessarily mean that the state of the tropical Pacific–northern Pacific teleconnection does not generate a local effect, even in the presence of an anomaly that is opposite to the canonical effect.
The atmospheric field in the northern Pacific, which is the key parameter driving upwelling, showed coherence in the spring–summer period in some time series. The spring–summer period is the key parameter driving upwelling, as the Aleutian Low and northern Pacific High pressure systems, as well as the northern Pacific circulation, play a critical role in ENSO development [76]. Our results thus suggest a complex one-way linear connection between ENSO and the upwelling response on the west coast of the Baja California peninsula. According to recent studies, this tropical–extratropical teleconnection was stronger in the late 1990s [76], since the CP ENSO became more frequent [33,34,54,56,84].

4.3. ENSO Development and Conditions in Pacific for Coherence Events

To understand why some ENSO events do not show typical canonical patterns on Ekman transport and how the development of previous tropical–extratropical Pacific conditions contributed to upwelling anomalies in the Baja California peninsula, we complemented our analysis with direct evaluation of the composite evolution of SST and 10-wind anomaly fields during ENSO (development year) events with significant wavelet coherence and events without it. Table A1 shows ENSO events and the lag of EP and CP ENSO events and Ekman transport for each region.
The evolution of El Niño events with significant wavelet coherence with Ekman transport in the Baja California peninsula shows positive February–April SST anomalies elongating from the CP to the subtropical eastern northern Pacific, as well as negative SST anomalies in the central subtropical Pacific region (Figure 8A). Additionally, western wind anomalies in CP, the meridional Pacific, and the south of Baja California peninsula are shown. This pattern is similar to the positive NPMM pattern [85], characterized by spring positive CP SST anomalies and negative EP SST anomalies and the weakening of the northern Pacific High [53]. Six El Niño events with significant coherence with upwelling intensity showed an NPMM pattern that was recently described as a key process to trigger northern Pacific variability through SST westerly anomalies and the southwestward extension of the northern Pacific positive SST in the following summer season (Figure 8C) (June–August), which is driven by wind–evaporation–SST feedback [79]. Previously, western anomalies near 170° E (Figure 8A) and northern positive SST anomalies were described as key factors of CP ENSO development [51,52], both of which are characteristics of NPMM. Thus, NPMM’s teleconnection with ENSO implies positive feedback, which excites mutually, especially with CP SST anomalies [51,86].
Additionally, the NPMM is related to the northern Pacific oscillation [27,87] through a seasonal foot-printing mechanism, in which SST anomalies are generated by midlatitude atmospheric variability in the previous winter force tropical atmosphere in the spring and summer [88,89]. This teleconnection possibly modulates wind stress anomalies along the Baja California peninsula via tropical–extratropical ocean–atmosphere interactions at interannual timescales. Additionally, it is possible that the coherence between tropical Pacific anomalies and wind field anomalies that drive Ekman transport is triggered by the development of NPMM associated with a preconditioned state of the northern Pacific and sustained by ENSO development via extratropical teleconnection [27,36]. Furthermore, the difference in coherence between EP and CP has a non-linear Ekman transport response to warm–cold anomalies, highlighting the conceptual relevance of taking into account ENSO diversity and the development of each ENSO event.
In El Niño events that do not show phase coherence with Ekman transport in the study area, the pattern of SST and 10 m wind anomalies are totally different with well-developed SST anomalies in the tropical Pacific in the boreal winter, but with no evident teleconnection in oceanic and atmospheric patterns with the northern Pacific, and an apparent absence of wind–evaporation–SST feedback (Figure 8B,D). These observations suggest the importance of the El Niño development and the NPMM-ENSO positive feedback to drive an atmospheric response in terms of upwelling intensity in Ekman transport.
La Niña events that exhibit coherence with Ekman transport were preceded by spring negative CP anomalies elongating from the CP to the subtropical eastern northern Pacific (the positive subtropical central Pacific region, Figure 9A), and also by wind anomalies towards the equator in the west of the Baja California peninsula, such as a negative NPMM [85].
La Niña events without coherence between ENSO and Ekman transport were preceded mainly by the warming of the tropical Pacific (Figure 9B) and reduced wind in the east and north of the Baja California peninsula. On the other hand, La Niña events with coherence showed negative anomalies in the tropical Pacific and equatorward 10-wind anomalies. These two patterns reveal that during the development of La Niña events, the connection of the tropical Pacific with the eastern northern Pacific is opposite in events with or without coherence. It is important to highlight that in the first scenario, the CP negative anomalies during the development of ENSO and peak event (Figure 9A,C,D) are shown. In contrast, in the no coherence La Niña scenario the predominant anomalies are positive during February–March (Figure 9B). During June–August and October–December, SST anomalies were located mainly in the EP (Figure 9D,F), highlighting the non-linearity of CP and EP events that have been documented in previous studies [39,47].

4.4. ENSO and Biological Response: Interannual Coherence

To analyze whether changes in ENSO diversity that led to changes in the upwelling centers could be predicted from sea surface Chl-a concentration, through an approach similar to that used for Ekman transport, we used satellite sea surface Chl-a concentration values extracted from each upwelling center. Wavelet coherence was calculated for the EP and CP indices, together with sea surface Chl-a monthly standardized anomalies at each upwelling center (Figure A3, Figure A4, Figure A5 and Figure A6). To perform this analysis, we only considered coherences for periodicities over 1.2 years. When EP La Niña events occurred that were followed by CP El Niño and mixed-type La Niña events, for the EP index, the coherence analysis was limited to 2006–2009 and a 2-year periodicity in the north of Punta Eugenia, in Punta Eugenia, and in Bahía Magdalena. In 2006–2009, anomalous cold water and high-sea-surface Chl-a concentrations were reported [15,74]. The coherence between Ekman transport and the EP index was only evident during this event in Ensenada, leading to an increase in upwelling. These observations suggest that EP ENSO events might be weakly connected to sea surface Chl-a concentration in upwelling centers. Previously, a low correlation between integrated Chl-a concentration (recorded in the column water) and the Multivariate ENSO Index (MEI.v2) was reported for the Punta Eugenia region [21]. However, our results show the highest coherence between the CP ENSO index and sea surface Chl-a concentration in upwelling centers on the west coast of the Baja California peninsula. The anti-phase coherence suggests that there is linear coherence with increased sea surface Chl-a concentration for CP La Niña events, and decreased sea surface Chl-a concentration for CP El Niño events, with 0- to 4-month lags (Figure A3C,D, Figure A4C,D, Figure A5C,D and Figure A6C,D) for the entire time series. Despite the differences in the non-linearities in ENSO–Ekman transport coherence, the Chl-a concentration shows a linear and consistent response to CP ENSO and a lower influence of EP ENSO on the sea surface Chl-a at upwelling centers in the Baja California peninsula. Differences in Ekman transport and sea surface Chl-a concentrations reflect the conceptual relevance of distinguishing ENSO events according to their development characteristics and diversity.
To complement our analysis, a composite analysis of surface Chl-a anomalies and wind anomalies was implemented for the same scenarios of Section 4.3 (Figure A7 and Figure A8). This analysis shows significant Chl-a anomalies in the coastal zone and in the regions where the upwelling centers are located. The coherence between ENSO and Ekman transport reflects changes in Chl-a only for the spring and summer seasons. Although the peak of ENSO events is usually in the boreal winter and wind anomalies are evident in this period, there was no response in Chl-a. This suggests that differences in Ekman coherence and Chl-a change between events: while ENSO can affect the wind field during other seasons, Chl-a will only respond to changes in wind during the spring and summer.
Furthermore, it is necessary to study not only the effects and impact of the NPMM on the biological response but also the effects on the oceanographic and biogeochemical characteristics that trigger EP ENSO and CP ENSO in the Baja California peninsula, which influence primary productivity and fisheries.

5. Conclusions

The upwelling centers located along the west coast of the Baja California peninsula, Mexico, show non-linearities in Ekman transport standardized anomalies and ENSO diversity in coherence over a time series. EP El Niño events show an increase in upwelling during boreal spring after the peak event, and CP El Niño and La Niña events trigger an Ekman transport decrease with latitudinal region differences. The non-linearities and Ekman transport response show a dependence on ENSO development and diversity, showing coherence between ENSO and Ekman transport when the CP and northeastern Pacific exhibit positive (negative) SST anomalies during the development of El Niño (La Niña) events in the spring and summer. The similarities with the northern NPMM pattern are possibly a key factor of atmospheric teleconnection between ENSO events and wind-forcing upwellings. The semiannual and annual coherences suggest an association with (1) other quasi-oscillatory processes, (2) changes in northern Pacific pressure gyres triggered by tropical–extratropical Pacific anomalies, or (3) ENSO precursor processes driven by the main variability modes, which are a matter for further study.
The sea surface Chl-a concentration shows coherence mainly with CP ENSO and shows a linear connection at interannual periodicities with anti-phase coherence: increased (decreased) sea surface Chl-a concentration during CP La Niña (CP El Niño) along the west coast of the Baja California peninsula upwelling centers. This suggests that changes in Ekman transport triggered by ENSO events are not the main factor of interannual changes in phytoplankton biomass.
Our results highlight the conceptual relevance of CP ENSO and ENSO diversity driving changes in sea surface Chl-a concentration. Further studies need to address ocean changes that drive ENSO diversity and its effects on upwelling centers in the Baja California peninsula to improve the identification of the impact of EP and CP ENSO on ecosystems and fisheries.

Author Contributions

Conceptualization, G.S.G.-C. and E.M.-A.; methodology, G.S.G.-C. and E.M.-A.; software, G.S.G.-C. and E.M.-A.; validation, G.S.G.-C., E.M.-A. and R.C.-D.; formal analysis, G.S.G.-C.; investigation, G.S.G.-C. and E.M.-A.; resources, G.S.G.-C. and E.M.-A.; data curation, G.S.G.-C. and E.M.-A.; writing—original draft preparation, G.S.G.-C.; writing—review and editing, G.S.G.-C., E.M.-A., R.C.-D., J.G.-G. and S.A.-G.; visualization, G.S.G.-C.; supervision, E.M.-A., R.C.-D., L.T.-F., J.G.-G. and S.A.-G.; project administration, R.C.-D. and S.A.-G.; funding acquisition, R.C.-D., J.G.-G. and S.A.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Consejo Nacional de Humanidades Ciencia y Tecnología (CONAHCyT), Instituto Politécnico Nacional SIP-IPN 20230854 institutional multidisciplinary research project “Perturbaciones ecosistémicas por impactos naturales y antropogénicos en ecosistemas costeros de México” and Instituto Politécnico Nacional SIP-IPN 20240636 institutional research project “Características oceanográficas y calidad del agua en mares y costas del noroeste de México”. Instituto Politécnico Nacional SIP-IPN full supported the Open Access cost of the present publication.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are freely available: wind speed and direction (https://www.remss.com/measurements/ccmp/, accessed on 18 November 2023), sea surface Chl-a concentration (https://cds.climate.copernicus.eu, accessed on 18 November 2023), and SST (https://coastwatch.pfeg.noaa.gov/erddap/griddap/nceiErsstv5.html, accessed on 18 December 2023). EM-A corresponding author can provide processing algorithms and routines under request.

Acknowledgments

GGC thanks CONAHCyT for an MSc scholarship (2022–2024), Instituto Politécnico Nacional and the Programa Institucional de Formación de Investigadores (PIFI) for the BEIFI-IPN scholarship as part of the SIP 20230854 institutional multidisciplinary research project: “Perturbaciones ecosistémicas por impactos naturales y antropogénicos en ecosistemas costeros de México” during August–December 2022, August–December 2023 and January–July 2024. EM-A thanks “Programa de estancias posdoctorales del CONAHCyT, investigadores e investigadoras por México convocatoria 2022(2)”. LT-F thanks CONAHCYT Catedras Patrimoniales program. JG-G, RC-D, SA-G thank the EDI-IPN and COFAA-IPN fellowships. EM-A, LT-F, JG-G, RC-D, SA-G thank the CONAHCYT-SNI fellowship. The authors thank NASA, ECMWF Earth observation component, the Copernicus Marine Service and Remote Sensing Systems for high-resolution data availability. All authors thank the anonymous reviewers for their helpful suggestions.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Figure A1. Composite anomalies of SST during December–February: (A) EP El Niño, (B) CP El Niño, (C) EP La Niña, and (D) CP La Niña from Table 1. Data: ERSSTV5 and 1991–2020 climatology were used.
Figure A1. Composite anomalies of SST during December–February: (A) EP El Niño, (B) CP El Niño, (C) EP La Niña, and (D) CP La Niña from Table 1. Data: ERSSTV5 and 1991–2020 climatology were used.
Jmse 12 01061 g0a1
Figure A2. Phase coherence and interpretation of wavelet coherence phase difference A r g T x , y τ , a arrows [68]. Shows the arrow (black) directions and their respective interpretation. When the phase shows in-phase (AC) or anti-phase (DF). When there is no lag between both series (A,D), when leading x (B,E) and when leading y (C,F).
Figure A2. Phase coherence and interpretation of wavelet coherence phase difference A r g T x , y τ , a arrows [68]. Shows the arrow (black) directions and their respective interpretation. When the phase shows in-phase (AC) or anti-phase (DF). When there is no lag between both series (A,D), when leading x (B,E) and when leading y (C,F).
Jmse 12 01061 g0a2
Figure A3. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly sea Surface Chl-a concentration in Ensenada upwelling centers: (A) EP San Antonio del Mar, (B) EP La Bocana, (C) CP San Antonio del Mar, and (D) CP La Bocana. Right arrows represent in-phase and left arrows represent anti-phase.
Figure A3. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly sea Surface Chl-a concentration in Ensenada upwelling centers: (A) EP San Antonio del Mar, (B) EP La Bocana, (C) CP San Antonio del Mar, and (D) CP La Bocana. Right arrows represent in-phase and left arrows represent anti-phase.
Jmse 12 01061 g0a3
Figure A4. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly sea surface Chl-a concentration in north of Punta Eugenia upwelling centers: (A) EP María, (B) EP Sta. Prieta, (C) CP María, and (D) CP Sta. Prieta. Right arrows represent in-phase and left arrows represent anti-phase.
Figure A4. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly sea surface Chl-a concentration in north of Punta Eugenia upwelling centers: (A) EP María, (B) EP Sta. Prieta, (C) CP María, and (D) CP Sta. Prieta. Right arrows represent in-phase and left arrows represent anti-phase.
Jmse 12 01061 g0a4
Figure A5. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly sea surface Chl-a concentration in Punta Eugenia upwelling centers: (A) EP San Pablo, (B) EP Clambey, (C) CP San Pablo, and (D) CP Clambey. Right arrows represent in-phase and left arrows represent anti-phase.
Figure A5. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly sea surface Chl-a concentration in Punta Eugenia upwelling centers: (A) EP San Pablo, (B) EP Clambey, (C) CP San Pablo, and (D) CP Clambey. Right arrows represent in-phase and left arrows represent anti-phase.
Jmse 12 01061 g0a5
Figure A6. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly sea surface Chl-a concentration in Bahía Magdalena upwelling centers: (A) EP Isla Sta. Margarita, (B) EP Pto. Magdalena, (C) CP Isla Sta. Margarita, and (D) CP Pto. Magdalena. Right arrows represent in-phase and left arrows represent anti-phase.
Figure A6. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly sea surface Chl-a concentration in Bahía Magdalena upwelling centers: (A) EP Isla Sta. Margarita, (B) EP Pto. Magdalena, (C) CP Isla Sta. Margarita, and (D) CP Pto. Magdalena. Right arrows represent in-phase and left arrows represent anti-phase.
Jmse 12 01061 g0a6
Table A1. A classification of ENSO events using the approach of [63] and phase differences of the eastern Pacific (EP) or central Pacific (CP) indices for the boreal winter December–January–February and Ekman transport for upwelling centers in the west coast of the Baja California peninsula. The ERSST v.5 database and climatology 1991–2020 were used to calculate both indices. The time month lag for coherence between EP Ekman transport and CP Ekman transport are shown for the study regions Ensenada (ENS), the north of Punta Eugenia (NPE), Punta Eugenia (PE) and Bahía Magdalena (BM). * denotes that the type of ENSO and the index were not the same (e.g., 1994/95 are CP ENSO events but showed coherence with the EP index). Bold numbers indicate that ENSO led coherence, and for the opposite, Ekman led in wavelet coherence.
Table A1. A classification of ENSO events using the approach of [63] and phase differences of the eastern Pacific (EP) or central Pacific (CP) indices for the boreal winter December–January–February and Ekman transport for upwelling centers in the west coast of the Baja California peninsula. The ERSST v.5 database and climatology 1991–2020 were used to calculate both indices. The time month lag for coherence between EP Ekman transport and CP Ekman transport are shown for the study regions Ensenada (ENS), the north of Punta Eugenia (NPE), Punta Eugenia (PE) and Bahía Magdalena (BM). * denotes that the type of ENSO and the index were not the same (e.g., 1994/95 are CP ENSO events but showed coherence with the EP index). Bold numbers indicate that ENSO led coherence, and for the opposite, Ekman led in wavelet coherence.
EventPhaseTypeCoherence Lead ET (Time Months)Coherence Chl-a
ENSNPEPEBMENSNPEPEBM
1994/95El NiñoCP1–2 * (1–2)0 *0–4 * ----
1996/97La NiñaEP ----
1997/98El NiñoEP3–40–24–51–5----
1998/99La NiñaCP 1
1999/00La NiñaCPc20–2 12 * (2)0 *0 *
2001/01La NiñaCP1–2 * (2–3)2 0 *0 *
2002/03El NiñoCP 0 *
2004/05El NiñoCP 1 * (2)
2007/08La NiñaCP 0 1–2 *
2008/09La NiñaCP0 *0 * 1–2 *1 *1–4
2009/10El NiñoCP 1–40–20–1 * (4–5)41–2 * (2)1 (3–4 *)
I confirm2010/11La NiñaCP0 1–2 *21 *1–2
2011/12La NiñaCP0 1–2 *0–2 * (4)1 *2–3
2014/15El NiñoCP 0 *4 * 4–5 1–2 *0 * (2–3)
2015/16El NiñoEP1–2 *4 *0–41–2 *
2017/18La NiñaEP4–51–2
2018/19El NiñoCP 0 * (0–1)
2019/20El NiñoCP 1 *5 0–1
2020/21La NiñaEP 0 *
2021/22La NiñaCP 0 0 *
Figure A7. Composite of sea surface Chl-a and 10 m wind anomalies (arrows) for El Niño events that show coherence with Ekman transport (A,C,E) and ENSO events that do not show significant coherence (B,D,F). SST anomalies are observed in February–April (A,B), June–August (C,D), and October–December (E,F). Light dots represent a 95% confidence level using two-tailed t-test. The 10 m wind anomalies below 80% of a two-tailed Student’s t-test were omitted.
Figure A7. Composite of sea surface Chl-a and 10 m wind anomalies (arrows) for El Niño events that show coherence with Ekman transport (A,C,E) and ENSO events that do not show significant coherence (B,D,F). SST anomalies are observed in February–April (A,B), June–August (C,D), and October–December (E,F). Light dots represent a 95% confidence level using two-tailed t-test. The 10 m wind anomalies below 80% of a two-tailed Student’s t-test were omitted.
Jmse 12 01061 g0a7
Figure A8. Composite sea surface Chl-a and 10 m wind anomalies (arrows) for La Niña development that show coherence with Ekman transport (A,C,E) and events that do not show significant coherence (B,D,F). There are SST anomalies shown for February–April (A,B), June–August (C,D), and October–December (E,F). Light dots represent a 95% confidence level using a two-tailed t-test. The 10 m wind anomalies below 80% are omitted.
Figure A8. Composite sea surface Chl-a and 10 m wind anomalies (arrows) for La Niña development that show coherence with Ekman transport (A,C,E) and events that do not show significant coherence (B,D,F). There are SST anomalies shown for February–April (A,B), June–August (C,D), and October–December (E,F). Light dots represent a 95% confidence level using a two-tailed t-test. The 10 m wind anomalies below 80% are omitted.
Jmse 12 01061 g0a8

References

  1. Bonino, G.; Di Lorenzo, E.; Masina, S.; Iovino, D. Interannual to decadal variability within and across the major Eastern Boundary Upwelling Systems. Sci. Rep. 2019, 9, 19949. [Google Scholar] [CrossRef] [PubMed]
  2. Pauly, D.; Christensen, V. Primary production sustain global fisheries. Nature 1995, 347, 255–257. [Google Scholar] [CrossRef]
  3. Messié, M.; Chavez, F.P. Seasonal regulation of primary production in eastern boundary upwelling systems. Prog. Oceanogr. 2015, 134, 1–18. [Google Scholar] [CrossRef]
  4. Chavez, F.P.; Messié, M. A comparison of Eastern Boundary Upwelling Ecosystems. Prog. Oceanogr. 2009, 83, 80–96. [Google Scholar] [CrossRef]
  5. Wu, W.; Zhai, F.; Gu, Y.; Liu, C.; Li, P. Weak local upwelling may elevate the risks of harmful algal blooms and hypoxia in shallow waters during the warm season. Environ. Res. Lett. 2023, 18, 114031. [Google Scholar] [CrossRef]
  6. Villegas, N.; Málikov, I.; Díaz, D. Variabilidad mensual de la velocidad de surgencia y clorofila a en la región del Panama Bight. Rev. Mutis 2016, 6, 82–94. [Google Scholar] [CrossRef]
  7. Capet, X.; Estrade, P.; Machu, E.; Ndoye, S.; Grelet, J.; Lazar, A.; Marié, L.; Dausse, D.; Brehmer, P. On the dynamics of the southern Senegal upwelling center: Observed variability from synoptic to superinertial scales. J. Phys. Oceanogr. 2017, 47, 155–180. [Google Scholar] [CrossRef]
  8. Kämpf, J.; Doubell, M.; Griffin, D.; Matthews, R.L.; Ward, T.M. Evidence of a large seasonal coastal upwelling system along the southern shelf of Australia. Geophys. Res. Lett. 2004, 31, 4–7. [Google Scholar] [CrossRef]
  9. Dang, X.; Bai, Y.; Gong, F.; Chen, X.; Zhu, Q.; Huang, H.; He, X. Different responses of phytoplankton to the ENSO in two upwelling systems of the South China Sea. Estuaries Coasts 2022, 45, 485–500. [Google Scholar] [CrossRef]
  10. Xiao, F.; Wu, Z.; Lyu, Y.; Zhang, Y. Abnormal strong upwelling off the coast of southeast Vietnam in the late summer of 2016: A comparison with the case in 1998. Atmosphere 2020, 11, 940. [Google Scholar] [CrossRef]
  11. Rueda-Roa, D.T.; Ezer, T.; Muller-Karger, F.E. Description and mechanisms of the mid-year upwelling in the southern Caribbean Sea from remote sensing and local data. J. Mar. Sci. Eng. 2018, 6, 36. [Google Scholar] [CrossRef]
  12. Rueda-Roa, D.T.; Muller-Karger, F.E. The southern Caribbean upwelling system: Sea surface temperature, wind forcing and chlorophyll concentration patterns. Deep. Res. Part I Oceanogr. Res. Pap. 2013, 78, 102–114. [Google Scholar] [CrossRef]
  13. Durazo, R. Climate and upper ocean variability off Baja California, Mexico: 1997–2008. Prog. Oceanogr. 2009, 83, 361–368. [Google Scholar] [CrossRef]
  14. Durazo, R. Seasonality of the transitional region of the California Current System off Baja California. J. Geophys. Res. Ocean. 2015, 120, 1173–1196. [Google Scholar] [CrossRef]
  15. González-Silveira, A.; Santamaría-del-Ángel, E.; Camacho-Ibar, V.; López-Calderón, J.; Santander-Cruz, J.; Mercado-Santana, A. The effect of cold and warm anomalies on the phytoplankton pigment composition in waters off northern Baja California Peninsula (México). J. Mar. Sci. Eng. 2020, 8, 533. [Google Scholar] [CrossRef]
  16. Kahru, M.; Di Lorenzo, E.; Manzano-Sarabia, M.; Mitchell, B.G. Spatial and temporal statistics of sea surface temperature and chlorophyll fronts in the California Current. J. Plankton Res. 2012, 34, 749–760. [Google Scholar] [CrossRef]
  17. Legaard, K.R.; Thomas, A.C. Spatial patterns in seasonal and interannual variability of chlorophyll and sea surface temperature in the California Current. J. Geophys. Res. Ocean. 2006, 111, 1–21. [Google Scholar] [CrossRef]
  18. Herrera-Cervantes, H.; Lluch-Cota, S.E.; Lluch-Cota, D.B.; Gutiérrez-de-Velasco, G. Interannual correlations between sea surface temperature and concentration of chlorophyll pigment off Punta Eugenia, Baja California, during different remote forcing conditions. Ocean Sci. 2014, 10, 345–355. [Google Scholar] [CrossRef]
  19. Todd, R.E.; Rudnick, D.L.; Davis, R.E.; Ohman, M.D. Underwater gliders reveal rapid arrival of El Niño effects off California’s coast. Geophys. Res. Lett. 2011, 38, 1–5. [Google Scholar] [CrossRef]
  20. Zaytsev, O.; Cervantes-Duarte, R.; Montante, O.; Gallegos-Garcia, A. Coastal upwelling activity on the Pacific shelf of the Baja California Peninsula. J. Oceanogr. 2003, 59, 489–502. [Google Scholar] [CrossRef]
  21. De la Cruz-Orozco, M.E.; Gómez-Ocampo, E.; Miranda-Bojórquez, L.E.; Cepeda-Morales, J.; Durazo, R.; Lavaniegos, B.E.; Espinosa-Carreón, T.L.; Sosa-Ávalos, R.; Aguirre-Hernández, E.; Gaxiola-Castro, G. Biomasa y producción del fitoplancton frente a la Península de Baja California: 1997–2016. Cienc. Mar. 2017, 43, 217–228. [Google Scholar] [CrossRef]
  22. Espinosa-Carreon, T.L.; Strub, P.T.; Beier, E.; Ocampo-Torres, F.; Gaxiola-Castro, G. Seasonal and interannual variability of satellite-derived chlorophyll pigment, surface height, and temperature off Baja California. J. Geophys. Res. Ocean. 2004, 109, 1–20. [Google Scholar] [CrossRef]
  23. Kahru, M.; Mitchell, B.G. Influence of the 1997-98 El Niño on the surface chlorophyll in the California Current. Geophys. Res. Lett. 2000, 27, 2937–2940. [Google Scholar] [CrossRef]
  24. Kahru, M.; Mitchell, B.G. Influence of the El Niño-La Niña cycle on satellite-derived primary production in the California Current. Geophys. Res. Lett. 2002, 29, 2–5. [Google Scholar] [CrossRef]
  25. Sosa-Ávalos, R.; Durazo, R.; Mitchell, B.G.; Cepeda-Morales, J.; Gaxiola-Castro, G. Parámetros fotosintéticos frente a Baja California: Una herramienta para estimar la producción primaria con datos de sensores remotos. Cienc. Mar. 2017, 43, 157–172. [Google Scholar] [CrossRef]
  26. Bograd, S.J.; Lynn, R.J. Physical-biological coupling in the California Current during the 1997-99 El Niño-La Niña cycle. Geophys. Res. Lett. 2001, 28, 275–278. [Google Scholar] [CrossRef]
  27. Capotondi, A.; Sardeshmukh, P.D.; Di Lorenzo, E.; Subramanian, A.C.; Miller, A.J. Predictability of US west coast ocean temperatures is not solely due to ENSO. Sci. Rep. 2019, 9, 10993. [Google Scholar] [CrossRef] [PubMed]
  28. Checkley, D.M.; Barth, J.A. Patterns and processes in the California Current System. Prog. Oceanogr. 2009, 83, 49–64. [Google Scholar] [CrossRef]
  29. Sung, M.K.; Kim, B.M.; An, S.I. Altered atmospheric responses to eastern Pacific and central Pacific El Niños over the North Atlantic region due to stratospheric interference. Clim. Dyn. 2014, 42, 159–170. [Google Scholar] [CrossRef]
  30. Patricola, C.M.; O’Brien, J.P.; Risser, M.D.; Rhoades, A.M.; O’Brien, T.A.; Ullrich, P.A.; Stone, D.A.; Collins, W.D. Maximizing ENSO as a source of western US hydroclimate predictability. Clim. Dyn. 2020, 54, 351–372. [Google Scholar] [CrossRef]
  31. Jacox, M.G.; Fiechter, J.; Moore, A.M.; Edwards, C.A. ENSO and the California Current coastal upwelling response. J. Geophys. Res. Ocean. 2015, 120, 1691–1702. [Google Scholar] [CrossRef]
  32. Ashok, K.; Yamagata, T. The El Niño with a difference. Nature 2009, 461, 481–484. [Google Scholar] [CrossRef] [PubMed]
  33. Capotondi, A.; Wittenberg, A.T.; Newman, M.; Di Lorenzo, E.; Yu, J.Y.; Braconnot, P.; Cole, J.; Dewitte, B.; Giese, B.; Guilyardi, E.; et al. Understanding ENSO diversity. Bull. Am. Meteorol. Soc. 2015, 96, 921–938. [Google Scholar] [CrossRef]
  34. Kao, H.Y.; Yu, J.Y. Contrasting Eastern-Pacific and Central-Pacific types of ENSO. J. Clim. 2009, 22, 615–632. [Google Scholar] [CrossRef]
  35. Takahashi, K.; Montecinos, A.; Goubanova, K.; Dewitte, B. ENSO regimes: Reinterpreting the canonical and Modoki El Niño. Geophys. Res. Lett. 2011, 38, 1–5. [Google Scholar] [CrossRef]
  36. Di Lorenzo, E.; Xu, T.; Zhao, Y.; Newman, M.; Capotondi, A.; Stevenson, S.; Amaya, D.J.; Anderson, B.T.; Ding, R.; Furtado, J.C.; et al. Modes and mechanisms of Pacific decadal-scale variability. Ann. Rev. Mar. Sci. 2023, 15, 249–275. [Google Scholar] [CrossRef] [PubMed]
  37. Okumura, Y.M. ENSO diversity from an atmospheric perspective. Curr. Clim. Change Rep. 2019, 5, 245–257. [Google Scholar] [CrossRef]
  38. Zhang, W.; Wang, Z.; Stuecker, M.F.; Turner, A.G.; Jin, F.F.; Geng, X. Impact of ENSO longitudinal position on teleconnections to the NAO. Clim. Dyn. 2018, 52, 257–274. [Google Scholar] [CrossRef]
  39. Gutiérrez-Cárdenas, G.S.; Díaz, D.C. El Niño-Southern Oscillation diversity and its relationship with the North Atlantic Oscillation—Atmospheric anomalies response over the North Atlantic and the Pacific. Atmosfera 2024, 38, 129–149. [Google Scholar] [CrossRef]
  40. Kim, H.-M.; Webster, P.J.; Curry, J.A. Ocean warming on North Atlantic tropical cyclones. Science 2009, 77, 77–80. [Google Scholar] [CrossRef]
  41. Beltrán, L.; Díaz, D.C. Oscilaciones macroclimáticas que afectan la oferta hídrica en la Cuenca del Río Gachaneca; Boyacá-Colombia. Rev. Bras. Meteorol. 2020, 35, 171–185. [Google Scholar] [CrossRef]
  42. Navarro-Monterroza, E.; Arias, P.A.; Vieira, S.C. El Niño/Southern Oscillation Modoki and its effects on the spatiotemporal variability of precipitation in Colombia. Rev. Acad. Colomb. Cienc. Exactas Fís. Nat. 2019, 43, 120–132. [Google Scholar] [CrossRef]
  43. Salas, H.D.; Poveda, G.; Mesa, Ó.J.; Marwan, N. Generalized synchronization between ENSO and hydrological variables in Colombia: A recurrence quantification approach. Front. Appl. Math. Stat. 2020, 6, 3. [Google Scholar] [CrossRef]
  44. Tedeschi, R.G.; Cavalcanti, I.F.A.; Grimm, A.M. Influences of two types of ENSO on South American precipitation. Int. J. Climatol. 2013, 33(6), 1382–1400. [Google Scholar] [CrossRef]
  45. Cai, W.; McPhaden, M.J.; Grimm, A.; Rodrigues, R.; Taschetto, A.; Garreaud, R.; Dewitte, B.; Germán, P.; Ham, Y.-G.; Santoso, A.; et al. Climate impacts of the El Niño—Southern Oscillation on South America. Nat. Rev. Earth Environ. 2020, 1, 215–231. [Google Scholar] [CrossRef]
  46. Li, Y.; Lau, N.C. Impact of ENSO on the atmospheric variability over the North Atlantic in late winter-role of transient eddies. J. Clim. 2012, 25, 320–342. [Google Scholar] [CrossRef]
  47. Zhang, W.; Wang, L.; Xiang, B.; Qi, L.; He, J. Impacts of two types of La Niña on the NAO during boreal winter. Clim. Dyn. 2015, 44, 1351–1366. [Google Scholar] [CrossRef]
  48. Chen, M.; Yu, J.-Y.; Wang, X.; Jiang, W. The Changing impact mechanisms of a diverse El Niño on the western Pacific Subtropical high. Geophys. Res. Lett. 2019, 46, 956–962. [Google Scholar] [CrossRef]
  49. Racault, M.F.; Sathyendranath, S.; Brewin, R.J.W.; Raitsos, D.E.; Jackson, T.; Platt, T. Impact of El Niño variability on oceanic phytoplankton. Front. Mar. Sci. 2017, 4, 133. [Google Scholar] [CrossRef]
  50. Lehodey, P.; Bertrand, A.; Hobday, A.J.; Kiyofuji, H.; McClatchie, S.; Menkès, C.E.; Pilling, G.; Polovina, J.; Tommasi, D. ENSO Impact on Marine Fisheries and Ecosystems. Chapter 19. In El Niño Southern Oscillation in a Changing Climate; McPhaden, M.J., Santoso, A., Cai, W., Eds.; AGU Publications: Washington, DC, USA, 2020; pp. 429–451. [Google Scholar] [CrossRef]
  51. Yu, J.-Y.; Kao, H.-Y.; Lee, T. Subtropics-related interannual sea surface temperature variability in the central Equatorial Pacific. Am. Meteorol. Soc. 2010, 23, 2869–2884. [Google Scholar] [CrossRef]
  52. Wang, X.; Guan, C.; Xin, R.; Wei, H.; Wang, L. The roles of tropical and subtropical wind stress anomalies in the El Niño Modoki onset. Clim. Dyn. 2019, 52, 6585–6597. [Google Scholar] [CrossRef]
  53. Amaya, D.J. The Pacific Meridional Mode and ENSO: A Review. Curr. Clim. Change Rep. 2019, 5, 296–307. [Google Scholar] [CrossRef]
  54. Kug, J.S.; Jin, F.F.; An, S. Il Two types of El Niño events: Cold tongue El Niño and warm pool El Niño. J. Clim. 2009, 22, 1499–1515. [Google Scholar] [CrossRef]
  55. Yeh, S.W.; Kug, J.S.; Dewitte, B.; Kwon, M.H.; Kirtman, B.P.; Jin, F.F. El Niño in a changing climate. Nature 2009, 461, 511–514. [Google Scholar] [CrossRef]
  56. Lee, T.; McPhaden, M.J. Increasing intensity of El Niño in the central-equatorial Pacific. Geophys. Res. Lett. 2010, 37, 1–5. [Google Scholar] [CrossRef]
  57. López-Aviles, B.; Beier, E.; Duran, R.; Gómez-Valdés, J.; Castro, R.; Sánchez-Velasco, L. The California Current System off Baja California Sur. Prog. Oceanogr. 2024, 222, 103225. [Google Scholar] [CrossRef]
  58. Kurczyn, J.A.; Beier, E.; Lavín, M.F.; Chaigneau, A. Mesoscale eddies in the northeastern Pacific tropical-subtropical transition zone: Statistical characterization from satellite altimetry. J. Geophys. Res. Ocean. 2012, 117, 1–17. [Google Scholar] [CrossRef]
  59. Espinosa-Carreón, T.L.; Gaxiola-Castro, G.; Beier, E.; Strub, P.T.; Kurczyn, J.A. Effects of mesoscale processes on phytoplankton chlorophyll off Baja California. J. Geophys. Res. Ocean. 2012, 117, C04005. [Google Scholar] [CrossRef]
  60. Castro, R.; Martínez, J.A. Variabilidad espacial y temporal del campo de viento. Dinámica del ecosistema pelágico frente a Baja California. In Dinámica del Ecosistema Pelágico Frente a Baja California 1997–2007, Diez años de investigaciones mexicanas de la Corriente de California; Gaxiola-Castro, G., Durazo, R., Eds.; Centro de Investigaciones Científicas y de Educación Superior de Ensenada: UABC; SEMARNAT: Instituto Nacional de Ecología: Ensenada, Mexico, 2010; pp. 129–147. ISBN 978-607-7908-30-2. [Google Scholar]
  61. Vazquez, H.J.; Gomez-Valdes, J. Wind events in a subtropical coastal upwelling region as detected by admittance analysis. Ocean Dyn. 2021, 71, 631. [Google Scholar] [CrossRef]
  62. Sullivan, A.; Luo, J.J.; Hirst, A.C.; Bi, D.; Cai, W.; He, J. Robust contribution of decadal anomalies to the frequency of central-Pacific El Niño. Sci. Rep. 2016, 6, 38540. [Google Scholar] [CrossRef]
  63. Santamaría-del-Ángel, E.; Sebastia-Frasquet, M.T.; González-Silveira, A.; Aguilar-Maldonado, J.; Mercado-Santana, A.; Herrera-Carmona, J.C. Uso potencial de las anomalías estandarizadas en la interpretación de fenómenos oceanográficos globales a escalas locales. In Tópicos de Agenda para la Sostenibilidad de Costas y Mares Mexicanos; Rivera-Arriaga, E., Sánchez-Gil, P., Gutiérrez, J., Eds.; Red RICOMAR; Universidad Autónoma de Campeche: Campeche, Mexico, 2019; pp. 193–211. ISBN 978-607-844-57-1. [Google Scholar]
  64. Atlas, R.; Hoffman, R.N.; Ardizzone, J.; Leidner, S.M.; Jusem, J.C.; Smith, D.K.; Gombos, D. A cross-calibrated, multiplatform ocean surface wind velocity product for meteorological and oceanographic applications. Bull. Am. Meteorol. Soc. 2011, 92, 157–174. [Google Scholar] [CrossRef]
  65. Ray, S.; Swain, D. A systematic method of estimation of alongshore windstress and Ekman transport associated with coastal upwelling. MethodsX 2023, 10, 102186. [Google Scholar] [CrossRef] [PubMed]
  66. Thomson, R.E.; Emery, W.J. Data Analysis in Physical Oceanography, 3rd ed.; Elsevier: Oxford, UK, 2014; ISBN 9780123877826. [Google Scholar]
  67. Díaz, D.; Villegas, N. Wavelet coherence between ENSO indices and two precipitation databases for the Andes region of Colombia. Atmosfera 2022, 35, 237–271. [Google Scholar] [CrossRef]
  68. Das, J.; Jha, S.; Goyal, M.K. On the relationship of climatic and monsoon teleconnections with monthly precipitation over meteorologically homogenous regions in India: Wavelet & global coherence approaches. Atmos. Res. 2020, 238, 104889. [Google Scholar] [CrossRef]
  69. Veleda, D.; Montagne, R.; Araujo, M. Cross-wavelet bias corrected by normalizing scales. Am. Meteorol. Soc. 2012, 29, 1401–1408. [Google Scholar] [CrossRef]
  70. Grinsted, A.; Moore, J.C.; Jevrejeva, S. Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear analysis of multivariate geoscientific data—Advanced methods, theory and application. Nonlinear Process. Geophys. 2004, 11, 561–566. [Google Scholar] [CrossRef]
  71. Torrence, C.; Compo, G.P. A Practical Guide to Wavelet Analysis. Bull. Am. Meteorol. Soc. 1998, 79(1), 61–78. [Google Scholar] [CrossRef]
  72. Barnard, P.L.; Hoover, D.; Hubbard, D.M.; Snyder, A.; Ludka, B.C.; Allan, J.; Kaminsky, G.M.; Ruggiero, P.; Gallien, T.W.; Gabel, L.; et al. Extreme oceanographic forcing and coastal response due to the 2015–2016 El Niño. Nat. Commun. 2017, 8, 6–13. [Google Scholar] [CrossRef]
  73. Cai, W.; Santoso, A.; Collins, M.; Dewitte, B.; Karamperidou, C.; Kug, J.-S.; Lengaigne, M.; McPhaden, M.J.; Stuecker, M.F.; Taschetto, A.S.; et al. Changing El Niño–Southern Oscillation in a warming climate. Nat. Rev. Earth Environ. 2021, 2, 628–644. [Google Scholar] [CrossRef]
  74. Linacre, L.; Durazo, R.; Hernández-Ayón, J.M.; Delgadillo-Hinojosa, F.; Cervantes-Díaz, G.; Lara-Lara, J.R.; Camacho-Ibar, V.; Siqueiros-Valencia, A.; Bazán-Guzmán, C. Temporal variability of the physical and chemical water characteristics at a coastal monitoring observatory: Station ENSENADA. Cont. Shelf Res. 2010, 30, 1730–1742. [Google Scholar] [CrossRef]
  75. McPhaden, M.J. Evolution of the 2002/03 El Niño. Am. Meteorol. Soc. 2004, 85(5), 677–695. [Google Scholar] [CrossRef]
  76. Chen, S.; Chen, W.; Yu, B.; Wu, R.; Graf, H.F.; Chen, L. Enhanced impact of the Aleutian Low on increasing the Central Pacific ENSO in recent decades. npj Clim. Atmos. Sci. 2023, 6, 29. [Google Scholar] [CrossRef]
  77. Anderson, B.T.; Perez, R.C.; Karspeck, A. Triggering of El Niño onset through trade wind—Induced charging of the equatorial Pacific. Geophys. Res. Lett. 2013, 40, 1212–1216. [Google Scholar] [CrossRef]
  78. Schwing, F.B.; Murphree, T.; Green, P.M. The evolution of oceanic and atmospheric anomalies in the northeast Pacific during El Niño and La Niña events of 1995-2001. Prog. Oceanogr. 2002, 54, 459–491. [Google Scholar] [CrossRef]
  79. Xie, S.P. A Dynamic ocean—atmosphere model of the Tropical Atlantic decadal variability. J. Clim. 1999, 12, 64–70. [Google Scholar] [CrossRef]
  80. Hartman, D. Pacific sea surface temperature and winter of 2014. Geophys. Res. Lett. 2015, 42, 1894–1902. [Google Scholar] [CrossRef]
  81. Chen, S.; Chen, W.; Wu, R.; Yu, B.; Graf, H.F. Potential impact of preceding Aleutian low variation on El Niño–Southern oscillation during the following winter. J. Clim. 2020, 33, 3061–3077. [Google Scholar] [CrossRef]
  82. Amaya, D.J.; Miller, A.J.; Xie, S. Kosaka Physical drivers of the summer 2019 North Pacific marine heathwave. Nat. Commun. 2019, 11, 1903. [Google Scholar] [CrossRef] [PubMed]
  83. Di Lorenzo, E.; Combes, V.; Keister, J.E.; Strub, P.T.; Thomas, A.C.; Franks, P.J.S.; Ohman, M.D.; Furtado, J.C.; Bracco, A.; Bograd, S.J.; et al. Synthesis of Pacific Ocean climate and ecosystem dynamics. Oceanography 2013, 26, 68–81. [Google Scholar] [CrossRef]
  84. Ashok, K.; Behera, S.K.; Rao, S.A.; Weng, H.; Yamagata, T. El Niño Modoki and its possible teleconnection. J. Geophys. Res. Ocean. 2007, 112, 1–27. [Google Scholar] [CrossRef]
  85. Zhao, J.; Kug, J.; Park, J.-H.; An, S.I. Diversity of North Pacific meridional mode and its distinct impacts on El Niño-Southern Oscillation. Geophys. Res. Lett. 2020, 47, e2020GL088993. [Google Scholar] [CrossRef]
  86. Yang, S.; Li, Z.; Yu, J.Y.; Hu, X.; Dong, W.; He, S. El Niño-Southern Oscillation and its impact in the changing climate. Natl. Sci. Rev. 2018, 5, 840–857. [Google Scholar] [CrossRef]
  87. Rogers, J.C. The North Pacific oscillation. J. Climatol. 1981, 1, 39–57. [Google Scholar] [CrossRef]
  88. Vimont, D.J.; Wallace, J.M.; Battisti, D.S. The seasonal footprinting mechanism in the Pacific: Implications for ENSO. J. Clim. 2003, 16, 2668–2675. [Google Scholar] [CrossRef]
  89. Vimont, D.J.; Battisti, D.S.; Hirts, A.C. Footprinting: A seasonal connection between the tropics and mid-latitudes. Geophys. Res. Lett. 2001, 28, 3923–3926. [Google Scholar] [CrossRef]
Figure 1. Main upwelling centers (red stars) identified along the west coast of the Baja California peninsula. The bathymetry (m) is shown. The dotted black line corresponds to the 200 m isobath.
Figure 1. Main upwelling centers (red stars) identified along the west coast of the Baja California peninsula. The bathymetry (m) is shown. The dotted black line corresponds to the 200 m isobath.
Jmse 12 01061 g001
Figure 2. El Niño–Southern Oscillation index in the 1993–2022 period for (A) eastern Pacific (EP, black line) and (B) central Pacific (CP, black line). Niño 3.4 standardized anomalies region denoting El Niño (red), La Niña (blue), and the neutral phase of ENSO (green).
Figure 2. El Niño–Southern Oscillation index in the 1993–2022 period for (A) eastern Pacific (EP, black line) and (B) central Pacific (CP, black line). Niño 3.4 standardized anomalies region denoting El Niño (red), La Niña (blue), and the neutral phase of ENSO (green).
Jmse 12 01061 g002
Figure 3. Boxplot showing the monthly mean and standard deviation of Ekman Transport, calculated between 1993 and 2022 for upwelling centers located along the west coast of the Baja California peninsula, Mexico. The median (blue line), quantile one (Q1) and three (Q3; box limits), interquartile range (IQR; length of the box), Q1 − 1.5 × IQR (lower boundary), Q3 + 1.5 × IQR (upper boundary), and outliers (circles) are shown.
Figure 3. Boxplot showing the monthly mean and standard deviation of Ekman Transport, calculated between 1993 and 2022 for upwelling centers located along the west coast of the Baja California peninsula, Mexico. The median (blue line), quantile one (Q1) and three (Q3; box limits), interquartile range (IQR; length of the box), Q1 − 1.5 × IQR (lower boundary), Q3 + 1.5 × IQR (upper boundary), and outliers (circles) are shown.
Jmse 12 01061 g003
Figure 4. Wavelet coherence C x , y 2 ( x , y ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly Ekman transport for Ensenada upwelling centers during EP and CP El Niño events: (A) EP San Antonio del Mar, (B) EP La Bocana, (C) CP San Antonio del Mar, and (D) CP La Bocana. Right arrow represents in-phase and the left arrow represents anti-phase.
Figure 4. Wavelet coherence C x , y 2 ( x , y ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly Ekman transport for Ensenada upwelling centers during EP and CP El Niño events: (A) EP San Antonio del Mar, (B) EP La Bocana, (C) CP San Antonio del Mar, and (D) CP La Bocana. Right arrow represents in-phase and the left arrow represents anti-phase.
Jmse 12 01061 g004
Figure 5. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for the ENSO index and monthly Ekman transport for north of Punta Eugenia upwelling centers during EP and CP El Niño events: (A) EP María, (B) EP Sta. Prieta, (C) CP María, and (D) CP Sta. Prieta. Right arrow represents in-phase and left arrow represents anti-phase.
Figure 5. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for the ENSO index and monthly Ekman transport for north of Punta Eugenia upwelling centers during EP and CP El Niño events: (A) EP María, (B) EP Sta. Prieta, (C) CP María, and (D) CP Sta. Prieta. Right arrow represents in-phase and left arrow represents anti-phase.
Jmse 12 01061 g005
Figure 6. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly Ekman transport for Punta Eugenia upwelling centers during EP and CP El Niño events: (A) EP San Pablo, (B) EP Clambey, (C) CP San Pablo, and (D) CP Clambey. Right arrows represent in-phase and left arrows represent anti-phase.
Figure 6. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for ENSO index and monthly Ekman transport for Punta Eugenia upwelling centers during EP and CP El Niño events: (A) EP San Pablo, (B) EP Clambey, (C) CP San Pablo, and (D) CP Clambey. Right arrows represent in-phase and left arrows represent anti-phase.
Jmse 12 01061 g006
Figure 7. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for the ENSO index and monthly Ekman transport for Bahía Magdalena upwelling centers during EP and CP El Niño events: (A) EP Isla Sta. Margarita, (B) EP Pto. Magdalena, (C) CP Isla Sta. Margarita, and (D) CP Pto. Magdalena. Right arrows represent in-phase and left arrows represent anti-phase.
Figure 7. Wavelet coherence ( C x , y 2 ( x , y ) ) (colors) and phase difference A r g T x , y τ , a (arrows) for the ENSO index and monthly Ekman transport for Bahía Magdalena upwelling centers during EP and CP El Niño events: (A) EP Isla Sta. Margarita, (B) EP Pto. Magdalena, (C) CP Isla Sta. Margarita, and (D) CP Pto. Magdalena. Right arrows represent in-phase and left arrows represent anti-phase.
Jmse 12 01061 g007
Figure 8. Composite of SST and 10 m wind anomalies (arrows) for El Niño development that show coherence with Ekman transport (A,C,E) and ENSO events that do not show significant coherence (B,D,F). SST anomalies are observed in February–April (A,B), June–August (C,D), and October–December (E,F). Light dots represent a 95% confidence level using a t-test for ERSSTV.5 NOAA and 1991–2020 climatology. The 10 m wind anomalies below an 80% significant two-tailed test were omitted.
Figure 8. Composite of SST and 10 m wind anomalies (arrows) for El Niño development that show coherence with Ekman transport (A,C,E) and ENSO events that do not show significant coherence (B,D,F). SST anomalies are observed in February–April (A,B), June–August (C,D), and October–December (E,F). Light dots represent a 95% confidence level using a t-test for ERSSTV.5 NOAA and 1991–2020 climatology. The 10 m wind anomalies below an 80% significant two-tailed test were omitted.
Jmse 12 01061 g008
Figure 9. Composite of SST and 10 m wind anomalies (arrows) for La Niña events that show coherence with Ekman transport (A,C,E) and ENSO events that do not show significant coherence (B,D,F). SST anomalies are observed in February–April (A,B), June–August (C,D), and October–December (E,F). Light dots represent a 95% confidence level using a t-test for ERSSTV.5 NOAA and 1991–2020 climatology. The 10 m wind anomalies below 80% of a two-tailed Student’s t-test were omitted.
Figure 9. Composite of SST and 10 m wind anomalies (arrows) for La Niña events that show coherence with Ekman transport (A,C,E) and ENSO events that do not show significant coherence (B,D,F). SST anomalies are observed in February–April (A,B), June–August (C,D), and October–December (E,F). Light dots represent a 95% confidence level using a t-test for ERSSTV.5 NOAA and 1991–2020 climatology. The 10 m wind anomalies below 80% of a two-tailed Student’s t-test were omitted.
Jmse 12 01061 g009
Table 1. ENSO events for EP and CP. The table shows the events in which the EP or CP index was higher than one SD for the ENSO peak in December–February.
Table 1. ENSO events for EP and CP. The table shows the events in which the EP or CP index was higher than one SD for the ENSO peak in December–February.
El NiñoLa Niña
EP1997/98, 2015/161996/97, 2017/18, 2021/22
CP1994/95, 2002/03, 2004/05, 2009/10, 2014/15, 2018/19, 2019/201998/99, 1999/00, 2000/01, 2007/08, 2008/09, 2010/11, 2011/12, 2020/21
Table 2. The mean and standard deviation of Ekman transport (m2 s−1) calculated at each coastal upwelling center located along the west coast of the Baja California peninsula, Mexico.
Table 2. The mean and standard deviation of Ekman transport (m2 s−1) calculated at each coastal upwelling center located along the west coast of the Baja California peninsula, Mexico.
RegionUpwelling CenterLong Term Mean
(m s−1)
Standard DeviationMonth
(Minimum Value)
Bahía MagdalenaIsla Santa Margarita−3.82 × 10−91.82 × 10−9June (−5.58 × 10−9)
Puerto Magdalena−4.18 × 10−91.96 × 10−9June (−5.94 × 10−9)
Punta EugeniaSan Pablo−2.34 × 10−90.97 × 10−9April (−3.43 × 10−9)
Clambey−3.12 × 10−91.24 × 10−9April (−4.71 × 10−9)
North of Punta EugeniaMaría−1.95 × 10−90.79 × 10−9April (−2.85 × 10−9)
Punta Prieta−3.38 × 10−91.35 × 10−9April (−4.85 × 10−9)
EnsenadaSan Antonio del Mar−1.85 × 10−104.13 × 10−10December (−1.98 × 10−9)
La Bocana−4.21 × 10−102.73 × 10−10December (−1.37 × 10−9)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gutiérrez-Cárdenas, G.S.; Morales-Acuña, E.; Tenorio-Fernández, L.; Gómez-Gutiérrez, J.; Cervantes-Duarte, R.; Aguíñiga-García, S. El Niño–Southern Oscillation Diversity: Effect on Upwelling Center Intensity and Its Biological Response. J. Mar. Sci. Eng. 2024, 12, 1061. https://doi.org/10.3390/jmse12071061

AMA Style

Gutiérrez-Cárdenas GS, Morales-Acuña E, Tenorio-Fernández L, Gómez-Gutiérrez J, Cervantes-Duarte R, Aguíñiga-García S. El Niño–Southern Oscillation Diversity: Effect on Upwelling Center Intensity and Its Biological Response. Journal of Marine Science and Engineering. 2024; 12(7):1061. https://doi.org/10.3390/jmse12071061

Chicago/Turabian Style

Gutiérrez-Cárdenas, Gabriel Santiago, Enrique Morales-Acuña, Leonardo Tenorio-Fernández, Jaime Gómez-Gutiérrez, Rafael Cervantes-Duarte, and Sergio Aguíñiga-García. 2024. "El Niño–Southern Oscillation Diversity: Effect on Upwelling Center Intensity and Its Biological Response" Journal of Marine Science and Engineering 12, no. 7: 1061. https://doi.org/10.3390/jmse12071061

APA Style

Gutiérrez-Cárdenas, G. S., Morales-Acuña, E., Tenorio-Fernández, L., Gómez-Gutiérrez, J., Cervantes-Duarte, R., & Aguíñiga-García, S. (2024). El Niño–Southern Oscillation Diversity: Effect on Upwelling Center Intensity and Its Biological Response. Journal of Marine Science and Engineering, 12(7), 1061. https://doi.org/10.3390/jmse12071061

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop