Many large-scale phenomena occur in the Pacific Ocean that act on varied time-scales, from the El Niño- Southern Oscillation (ENSO) with month-year cycles to decadal-to-multidecadal frequency events such as the Pacific Decadal Oscillation (PDO). Their effects are detectable in global to local ecosystems, and their aggregate contributions establish climatic shifts that have shown fluctuations in the atmospheric and oceanic conditions such as sea temperature [1] , wind fields associated with atmospheric pressure variations [2] , ocean currents [3] , coastal upwellings [4] , and mixed layer depths [5] . Many of these climatic variations occur via atmospheric- oceanic teleconnections [6] , which extend from the troposphere to the ocean surface, including fluctuations in sea level pressure that are closely linked to changes in the surface winds, sea surface temperature, heat content, mixed layer, and thermocline depth. Climate variability due to large-scale influences in the North Pacific receives considerable attention due to its impact on tropical and extra-tropical climate [7] and the weather over North America [8] .

The ENSO, perhaps the most studied large-scale phenomenon, is comprised of two phases. The warm El Niño phase, which has been studied extensively in the Pacific Ocean [6] [9] , is characterized by the weakening of trade winds, warming of the sea surface layer in the tropical-subtropical eastern Pacific, a switch from low to high atmospheric pressure near Darwin, Australia, and the opposite effect near the Tahiti Islands, which is termed the Southern Oscillation. The La Niña cold phase of the ENSO is much less studied in the Pacific. It is characterized by an intensification of the trade winds, the cooling of the sea surface layer in the tropical-sub- tropical Eastern Pacific, very low atmospheric pressure near Darwin, Australia, and very high pressure near the Tahiti Islands [10] [11] .

The Pacific Decadal Oscillation (PDO), another large-scale phenomenon, is a pattern of ocean variability over the entire Pacific that is similar to the ENSO in some respects, but with a much longer cycle [12] [13] . It is also defined by two phases; the positive phase in the North Pacific occurs when the sea surface temperature anomalies are cold in the central North Pacific and warm along the Pacific coast and when the sea level pressure is below average in the North Pacific [14] , while the converse occurs during the negative phase. Both phases are calculated by the standardized difference between sea surface temperatures in the north-central Pacific and Gulf of Alaska. The phases of the PDO may be important in enhancing or dampening the ENSO impacts [15] .

The Pacific/North American (PNA) pattern is one of the most prominent modes of low-frequency variability in the Northern hemisphere extra-tropics, and it describes the variation of atmospheric circulation patterns over the Pacific Ocean and North America that involve changes in the atmospheric pressure between the Aleutian Low and the high pressure over the Rocky Mountains [16] [17] . The PNA has the most impact on climate variables during the winter [18] [19] . As with the first two phenomena, the PNA also presents two phases, with the positive phase usually related to El Niño and the negative phase to La Niña. However, there is also a weak connection between the PNA and the ENSO [20] -[22] . During the positive phase of the PNA, the Aleutian Low in the intermountain region of North America shows anomalous low pressures, while the high pressure over the Rocky Mountains and Northeast Pacific Coast strengthens, increasing the pressure gradient between the two centers of circulation, in addition to increasing the wind speed, which steers more storms into the Northwest Pacific. The negative phase of the PNA is essentially the reverse pattern of the positive phase [16] .

The North Pacific pattern is quite similar to the PNA, but the variables involved in this atmospheric process are considered differently. For the PNA, [16] considered four centers of action in the mid-tropospheric height field, one each of one sign near Hawaii and along the west coast of North America, and one each of the opposite sign over the north Pacific and southeastern United States. According to [23] , these four centers of action and their resulting index have proven to be useful, but they also consider it inappropriate to weight the four centers of the PNA similarly, as the North Pacific has by far the most prominent height field. They proposed the North Pacific pattern and its index as a much more robust, but simpler measure of the atmospheric processes in North America.

The Pacific Ocean exerts a strong influence over the oceanographic conditions of the Gulf of California, the sea that shelters La Paz Bay, due its connection to the eastern tropical Pacific [24] . The dynamic forcing from the Pacific Ocean over the gulf is one of the most important oceanographic features, because it integrates relevant phenomena such as salt and heat global balances, thermohaline circulation, and barotropic ocean circulation [25] [26] .

Studies have demonstrated the effects of these large-scale phenomena on the Gulf of California region, and satellite analysis has shown the influence of the ENSO over the sea surface temperature variability. Based on the satellite records of the sea surface temperature anomalies, [27] [28] showed that during El Niño, the temperature in the gulf registered up to 3˚C above normal, while during La Niña, it was up to 3˚C below normal.

Similarly, studies on the PDO have shown its influence on pressure, wind, temperature, and the precipitation patterns of the North Pacific [12] [29] [30] . Its temporal modulations are linked to several important biological and ecosystem variables in the ocean [31] [32] . Nevertheless, other parameters such as the decadal fluctuations in salinity, nutrients, and chlorophyll a in the eastern North Pacific are often poorly correlated with the PDO [33] .

The relationship between the sea surface temperature over the Pacific Ocean and the PNA has also been studied by [34] - [36] , who determined that the PNA pattern is positively correlated with the sea surface temperature anomalies over the tropical Pacific [34] [37] and negatively correlated with those anomalies over the North Pacific [38] . The link between the PNA and the tropical sea surface temperature anomalies is mostly attributed to the ENSO phenomenon. During an El Niño event, atmosphere-ocean interactions generate a seesaw pattern with opposing anomalies of the surface pressure in the eastern and western tropical Pacific [36] .

On the other hand, a dependence of the western North Pacific (28˚ - 31.5˚N, 133˚ - 137˚E) mixed layer depth on decadal and inter-decadal time-scales (1950-1997) and large-scale, low-frequency North Pacific climatic variability has been reported [39] . These authors found an inverse relationship between the mixed layer depth and the NP index in 1962-1965 and 1973-1997 and suggested that its close linkage is caused by the alternating wind stress fields. Studies of the eastern North Pacific mixed layer depth variability (Baja California eastern boundary currents) are extremely important for evaluating the effects of local and basin-scale changes in the ocean surface. In particular, the low-frequency mixed layer depth changes in the eastern North Pacific are related to large-scale atmospheric forcing on the Pacific Ocean [40] .

Some large-scale phenomena effects on La Paz Bay have already been reported, primarily related to El Niño events, such as the sinking of the thermocline and isotherms and a warming of the entire water column (>4˚C) during the summer-autumn of 1997. Nevertheless, some inconsistencies trying to relate other physical events with the ENSO cycles are present, such as the unusual 100 m mixed layer depth during the winter of 2002 in a non-ENSO period [41] , and the absence of a mixed layer during the La Niña period in the winter of 1996 and 2009 [42] [43] , a period when theoretically it should have been dominated by strong and persistent winds.

It is important to mention that La Paz Bay is a coastal water body of economic interest due to its tourism industry and fishing activities that are maintained by the high biological productivity of the bay, as well as by its contiguity to La Paz City, the capital of Baja California Sur. Thus, a better understanding of the processes that affect that productivity such as the North Pacific large-scale phenomena is required, as most studies have only focused on the ENSO influence, particularly the temperature variations. Consequently, we analyzed the variability of two hydrographic characteristics, sea surface temperature from 2000 to 2010 and mixed layer depth from 1994 to 2009, to determine their relationships with the climatic phenomena of the North Pacific.

The sea surface temperature annual averages (2000-2010) showed that 2003, 2004, 2006, 2009, and 2010 were the warm years (Figure 2), and 2000, 2001, 2002, 2005, 2007, and 2008 were the cold years, with abrupt changes from 2006 to 2007 and from 2008 to 2009.

Two periods were clearly distinguished in the annual cycle of the sea surface temperature in La Paz Bay (Figure 3), a cold period with an average temperature range of 20.5˚C - 26.0˚C during December-May and a warm period of 26.0˚C - 31.0˚C during June-November. Figure 3 also shows the variations in 2007 and 2009, the coldest and warmest years of the study period. From January to April, the temperatures of both years were slightly above the mean, but from May to December the temperatures were above the mean in 2009 and below it in 2007, with the exception of October-November, when all of the temperatures were quite similar.

Table 1 shows the annual variability of the two sea surface temperature anomalies (SSTa) and the climate indices involved in this study, revealing that the SSTa variability is similar to that of the ONI and MEI climate indices in relation to the ENSO events. The Spearman correlation analysis corroborates this similarity, showing high correlations among the SSTa and the ONI and MEI indices.

Given the ONI had the highest correlation with the SSTa, we plotted both annual variabilities (Figure 4) and found a high similarity between them, reaffirming that the sea surface temperature changes occurring in the El Niño 3.4 zone have an influence on the changes in this parameter in La Paz Bay. However, a high discrepancy was observed in 2000 and 2002 when the sea surface temperature anomaly in the bay was negative and the ONI was positive and above 0.5.

In the case of the mixed layer depth variability, the larger depths were measured during December, February, and March, corresponding to the cold period of the bay, and during October and November, the months when temperature decreases begin in the bay. Standing out, the largest mixed layer depth of 100 m was during February 2002, while the common depth during this period was ~50 m (Table 2). This contrasts with the absence of the mixed layer during March of 1996 and 2009. Table 2 also shows the Spearman correlation analysis, which only found a statistically significant negative correlation between the mixed layer depth and the NP index.

Because the highest correlation was between the mixed layer and the NP index time series, both variables were plotted, revealing an inverse relationship between them (Figure 5). This negative and proportional relationship prevailed during most of the analysis period, with the exceptions of February 2002, when the 100 m mixed layer depth was registered and the relationship became very disproportional, and for March of 1996 and 2009, when the proportion of the negative relationship is unclear. Discarding these months the negative relationship becomes extremely high.

To analyze the anomalous conditions in the mixed layer depths described previously for 1996, 2002, and 2009, plots were made of the annual patterns for the NP, PNA, NOI, and SOI indices of these years (Figure 6). These plots show the contribution of the variability of each large-scale atmospheric related index that had an influence during the years of the anomalous mixed layer depths. In the case of the 100 m mixed layer depth in February 2002, the NP and SOI (Figure 6(a) and Figure 6(d)) values were lower than those in 2009, while the PNA and NOI (Figure 6(b) and Figure 6(c)) were higher than those in 2009. In the case of March 2009, the NP, NOI, and

SOI were higher than in 2002 (Figure 6(a), Figure 6(c) and Figure 6(d)). The differences among the indices for these winter months and the correlations shown in Table 2 demonstrated that, not only does the NP phenomenon influence the mixed layer depth, but other atmospheric phenomena such as the PNA, NOI, and SOI do as well. In the particular case of 1996, different patterns registered throughout the year from those in 2002 and 2009, but in March, the PNA and SOI (Figure 6(b) and Figure 6(d)) were higher than those of 2002 and 2009.

The annual variability of the La Paz Bay sea surface temperature showed that, although the bay is a semi-en- closed body of water, it has interannual variability and abrupt changes, such as occurred during the 2007 and 2009 years. This indicates La Paz Bay is highly dynamic and influenced by the North Pacific, which agrees with reports by [24] [41] [42] [45] [48] , who described the seasonal changes of this bay and emphasized the influence of the ocean-atmospheric processes that determine its hydrographic conditions.

The mean annual thermal cycle occurs in two phases (warm and cold) with their respective transition periods having similar temperatures. This differs somewhat from the findings of [41] [42] [48] , who consider the bay’s hydrographic variability based on four seasons. The differences in the annual thermal cycle between 2007 and 2009 showed that the changes in the temperature of La Paz Bay occurred not only during the El Niño events, when the warm waters from the tropical regions influence the marine environment, but also during La Niña, when the changes in La Paz Bay also occurred alongside decreases in the temperature during the warm phase of the annual cycle. This contrasts with reports by [41] [42] [48] [51] - [53] , who only focused their research efforts on the temperature increases in La Paz Bay during the El Niño events.

The anomalous surface temperature variations of La Paz Bay showed a high correlation among the ENSO indices, corroborating that the thermal fluctuations of the bay are caused by both phases of this phenomenon (El Niño/La Niña), primarily during the warm period of the annual thermal cycle of the bay through the entrance of surface waters from the tropical Pacific. This coincides with observations by [24] [41] [42] [45] [48] , who registered tropical waters during both the warm and cold periods of the La Paz Bay annual cycle. However, other studies exist showing that the La Niña effects are also reflected in sea surface temperature decreases in the Gulf of California, such as those by [27] [28] . Based on the sea surface temperature anomalies obtained from satellite records, they showed that during El Niño/La Niña events, the monthly mean temperatures of the Gulf of California vary by up to 3˚C from the normal conditions.

The variability of the mixed layer depth showed a high dependence on the wind conditions during the cold phase of the annual thermal cycle of the bay. [41] [45] showed that when the temperature decreases during the cold phase of the annual thermal cycle of the bay, the north winds velocities gradually rise, thereby increasing the mixed layer depth. Furthermore, [54] mentioned that the wind patterns in the Gulf of California correspond to a large-scale annual cycle of the Northeast Pacific. This was confirmed with the negative correlation obtained between the variability of the mixed layer depth and the NP index, showing that the wind patterns that regulate the depth’s dimensions are determined by the atmospheric pressure variability of the Northeast Pacific.

Although the NP fluctuations are the primary factor generating the mixed layer variations, other phenomena must also be taken into consideration, as their indices were not significantly correlated. Such was the case in 1996, 2002, and 2009 when the anomalous mixed layer depths were registered. These were caused by the pressure variations of two or more of the atmospheric phenomena we analyzed that resulted in the changes of the wind magnitude and direction. [40] attempted to relate the large-scale phenomena to the variability of the mixed layer depth at the tropical boundary of the California Current by using CTD observations from 1997-2007. They affirmed that the seasonal and interannual variability of the mixed layer depth was correlated with offshore Ekman transport and that the abrupt mixed layer depth change that occurred between January 1998 and January 2000 was associated with a strong El Niño-La Niña shift.

We thank the Instituto Politécnico Nacional (IPN), Consejo Nacional de Ciencia y Tecnología (CONACYT), and Comisión de Operación y Fomento de Actividades Académicas del IPN (COFAA) for financial help and scholarships, in addition to the Secretaría de Investigación y Posgrado-Becas de Estímulo Institucional de For- mación de Investigadores (SIP-BEIFI) and Becas de Estímulo al Desempeño de los Investigadores (SIP-EDI). We also thank American Manuscript Editors for editing this manuscript.