California Niño/Niña

The present study shows the existence of intrinsic coastal air-sea coupled phenomenon in the coastal ocean off Baja California and California in boreal summer for the first time. It contributes significantly to the interannual sea surface temperature (SST) anomalies there. An initial decrease/increase in the equatorward alongshore surface winds weakens/strengthens the coastal upwelling and raises/lowers the coastal SSTs through oceanic mixed-layer processes. The resultant coastal warming/cooling, in turn, heats/cools the overlying atmosphere anomalously, decreases/increases the atmospheric pressure in the lower troposphere, generates an anomalous cross-shore pressure gradient, and thus reinforces or maintains the alongshore surface wind anomalies. The regional air-sea coupled phenomenon seems to be analogous to the well-known El Niño/Southern Oscillation (ENSO) in the tropical Pacific but with much smaller time and space scales, and may be referred to as California Niño/Niña in its intrinsic sense.

C limatological surface winds blow southward all year round along the western coast of Baja California and California and drive the surface water offshore owing to the Ekman transport 1 originated in the Coriolis force that deflects moving objects to the right of their paths in the Northern Hemisphere. To compensate for the surface water, cold and nutrient-rich subsurface water upwells and increases biological productivity in the coastal surface layer, which supports abundant marine species at different trophic levels. In some years, however, the upwelling is reduced (enhanced) and thus the coastal ocean becomes warmer (colder) than usual [2][3][4][5] . For example, the El Niño years of 1957/58, 1976/77, 1982/83 and 1997/98 witnessed pronounced warming in the northeast Pacific; the anomalous warming lasted for over one season with the maximum monthly SST anomalies up to 2-6uC [2][3][4][5] . The anomalous warming caused general decline in biomass of phytoplankton and zooplankton and intrusion of warm-water species into higher-than-usual latitudes of habitats [6][7][8][9] . Extreme climate was also observed in the nearby continent; increased winter storms hit California in 1982/83 and resulted in one of the wettest winters in record. In contrast, the coastal ocean was extremely cold in the 1999 La Niña year and California experienced cooler-than-normal summer 10 . Researches on the frequent co-occurrences of the tropical ENSO and the coastal warming/cooling have shown that ENSO can influence the coastal ocean via atmospheric teleconnections that modify strength of the alongshore surface winds 4,11,12 and/or via coastally trapped Kelvin waves that propagate poleward from the eastern tropical Pacific and influence the coastal thermocline depth, upwelling and SSTs 13 . However, not all of the coastal warming/cooling events are associated with ENSO. The basin-wide atmospheric circulation anomalies in mid-latitudes related to the Pacific Decadal Oscillation and the Aleutian Low are also regarded as contributors to the coastal SST variations 12,14,15 .
In general, subtropical oceans are passively forced by the overlying atmosphere and have weak influences on the atmosphere, compared to the atmospheric internal variability, in seasonal to interannual time scales 16 . However, the SST anomalies in the coastal upwelling ocean may affect the land-sea thermal contrast and cross-shore pressure gradient in the lower troposphere, thus influencing the strength of the alongshore surface winds. This hypothesis associated with the feedback between the coastal upwelling and alongshore surface winds is adopted to partly explain the seasonal intensification of subtropical highs in summer in the Northern Hemisphere 17 . Recently, the coastal ocean-atmosphere coupled process is named the coastal Bjerknes feedback and applied to explain the regionally intensified SST variability in the coastal ocean off western Australia 18 . Here, we demonstrate that the same coastal Bjerknes feedback contributes significantly to the coastal SST variations off Baja California and California in boreal summer and identify the existence of a regional coupled mode as the California Niño/Niña in relation to the pioneering work of Simpson (1983) among others that describe the warming events in the coastal ocean off western North America as California El Niño 3,19 . The important difference from previous works, however, is that we demonstrate the existence of an intrinsic regional coupled mode that is independent of ENSO events.
In the present study, we adopt Monthly National Oceanic and Atmospheric Administration (NOAA) Optimum Interpolation Sea Surface Temperature (OISST) version 2 20 , Pacific Fisheries Environmental Laboratory (PEFL) derived upwelling indices 21 , and National Centers for Environmental Prediction (NCEP)/ National Center for Atmospheric Research (NCAR) reanalysis 1 22 for the period of January 1982 to December 2011. Monthly averaged anomalies in SSTs, upwelling indices and atmospheric variables are calculated by removing the mean seasonal cycles. As shown in Fig. 1a, the coastal SSTs off Baja California show substantial interannual variations. The spatial distribution of monthly standard deviation of SST is similar to the first mode of empirical orthogonal function (EOF) analysis that explains ,47% of total variances (Fig. 1b). The California Niño/Niña index is defined as the SST anomalies averaged over 110uW-120uW and 20uN-30uN (the enclosed coastal ocean in Fig. 1a) where they show the largest interannual variability. In the past three decades, the California Niño/Niña index (Fig. 1c) is significantly correlated with the time series of the first EOF mode with the correlation coefficient up to 0.9, significant at a 99.99% confidence level by the two-tailed t test. Hence, both time series represent almost the same interannual SST anomalies off Baja California and California.
Seasonally stratified standard deviations of the California Niño/ Niña indices show strong interannual variations from boreal summer to winter. As mentioned above, the California Niño/Niña is strongly influenced by the tropical ENSO. However, the influences are much stronger in boreal winter than summer; the correlation coefficient between January-March (July-September, JAS) Niño3 and California Niño/Niña indices is ,0.65 (,0.32) in the past 30 years, significant at the 99.99% (90%) confidence level by the two tailed t test. If the interannual variations of the California Niño/Niña indices simultaneously related to ENSO are linearly regressed out, the California Niño/Niña shows the highest interannual variability in boreal summer. Since the mixed-layer depth is shallower and the isotherms at the surface and in the upper ocean off Baja California are more densely packed in boreal summer than winter (Supplementary Fig. 1), the highest interannual variability in boreal summer may be more closely related to the possible coastal Bjerknes feedback. The standard deviation and first EOF mode based on the JAS SST anomalies ( Supplementary Fig. 2) show spatial patterns almost identical to those based on the monthly SST anomalies for the whole calendar year (Fig. 1a-b) except that the first EOF mode of the former explains a little higher percentage of total variances. Hence, hereafter, we focus only on the California Niño/Niña in boreal summer (JAS) and investigate the above hypothesis in this report.

Results
As shown in Fig. 2a, the anomalous negative sea level pressure (SLP) and cyclonic atmospheric circulation appear off the coast of California three months before the peak phase of the California Niño. The associated alongshore component of the surface wind anomalies is poleward and thus reduces the year-round coastal upwelling and offshore Ekman transport, resulting in the significant www.nature.com/scientificreports positive SST anomalies in the coastal ocean. Over the anomalous warm coastal ocean, air parcels in the planetary boundary layer (PBL) are thus anomalously heated (Fig. 3a). This causes negative SLP anomalies that maintain/enhance the anomalous cross-shore pressure gradient in the lower troposphere and the poleward alongshore surface wind anomalies, which in turn raises the coastal SSTs. The close coupling between the lower troposphere and the coastal ocean during the development phase of California Niño can be clearly seen in Fig. 4a, where the lead-lag correlation coefficients between SST, alongshore surface wind and upwelling indices are plotted. It is shown that at the peak time, the correlation coefficient between the upwelling (also the alongshore wind) and California Niño/Niña indices amounts to 20.7, demonstrating clearly that the coastal Bjerknes feedback contributes to about a half of the total variances of the California Niño/Niña in boreal summer. The California Niño/Niña may be also influenced by the tropical ENSO in boreal summer. Hence, significant SST anomalies can be seen in the eastern tropical Pacific as well in the correlation figures (Fig. 2a). However, we stress that the close regional air-sea coupling off Baja California during the developing phase exists clearly even after the ENSO influences are linearly removed (Figs. 2-4b).
The climatological equatorward surface winds are weakened during the developing phase of the California Niño, but possible contributions from the latent heat flux anomalies to the positive SST anomalies averaged in the region of our interest are limited ( Supplementary Figs. 3-4). This is probably because of the spatial inconsistence of the latent heat flux anomalies. In the inshore ocean, the anomalously warm SSTs due to oceanic processes may offset or even surpass the effect of the reduced wind speed on evaporation through changes of the near-surface specific humidity; the latent heat flux anomalies there are insignificant until two months before the peak, and then become significantly increased from the inshore ocean to the overlying atmosphere. In contrast, the latent heat flux anomalies contribute significantly to the positive SST anomalies in the further offshore ocean particularly at the southern edge of the SLP anomalies where the anomalous cyclonic surface winds weaken the climatological anticyclonic surface winds and thus reduce the latent heat fluxes from the ocean to the atmosphere (Fig. 2,  Supplementary Figs. 1, 4). In addition, the subtropical ocean is known to be covered by low stratus clouds due to trade wind inversion related to subsidence of the subtropical high 23 . The SST-stratus cloud feedback, through which the warmer-than-normal SST reduces the amount of stratus clouds, increases the downward shortwave radiation and thus enhances the warm SST anomalies, may also contribute to the positive SST anomalies during the developing phase of the California Niño 24 (Supplementary Fig. 4). However, the surface thermal forcing as a whole may play a secondary role compared to the oceanic processes on the SST anomalies related to California Niño during the developing phase; the highest correlation coefficient between the net surface heat flux anomalies averaged over the region of our interest and the summer California Niño/Niña indices is ,0.4 at most when the California Niño/Niña indices lag two months for the case where the linear ENSO influences are excluded ( Supplementary Fig. 3b). In contrast, during the decaying phase, the accumulated positive SST anomalies off Baja California and California are dampened quickly by the surface heat fluxes, in particular the latent heat fluxes, to the overlying atmosphere as shown by much higher correlation coefficients with the California Niño/Niña indices ( Supplementary Figs. 3-4).
The anomalous oceanic advection may also contribute to the development of California Niño 3 . The poleward alongshore surface wind anomalies off Baja California during the developing phase of California Niño may reduce the southward transportation of cold water by the California Current, i.e., the equatorward wind-driven surface current 25 . The wind anomalies decrease the offshore Ekman transport and increase the coastal sea surface heights (SSHs)  Supplementary Fig. 5). The anomalous cross-shore gradient of SSHs thus formed is consistent with the poleward surface current anomalies. These processes are comparable to the zonal advective feedback in the tropical Pacific during the ENSO events 26 . Hence, in a broad sense, the coastal ocean is analogous to the tropical ocean as pointed out first by Japanese oceanographer Kozo Yoshida and Chinese oceanographer Han-Lee Mao in collaboration about a half century ago at the Scripps Institution of Ocenography 27,28 when they found a similarity between the equatorial upwelling and the coastal upwelling as the ocean responses to the surface winds. Here we have extended the similarity to ocean-atmosphere coupled processes.

Discussion
Although the generation mechanism of the California Niño/Niña in boreal summer is quite similar to that of the tropical ENSO, to recognize differences is important, too. The former has much smaller temporal and spatial scales; it is confined to the coastal ocean and lasts only for 1-2 seasons (Figs. 2-3). This is probably because, besides the coastal upwelling, the shallow mixed-layer and sharp vertical and horizontal temperature gradients in the upper ocean are required for the coastal Bjerknes feedback to operate effectively. Therefore, the California Niño/Niña generated by the coastal Bjerknes feedback is seasonally locked to boreal summer within the coastal ocean in the subtropical transition zone. Also, we note that the nature of California Niño/Niña in boreal summer on which we have focused in this study differs from that of California Niño/ Niña in boreal winter [2][3][4][5]19 . The latter is generally referred to the coastal ocean warming/cooling in the northeast Pacific along almost the entire western coast of North America. The broadly-extended coastal SST anomalies are attributed mainly to the external forcing of the tropical ENSO via both oceanic and atmospheric teleconnections and/or the ENSO-independent basin-scale atmospheric circulation anomalies in mid-latitudes related to the Aleutian Low. Even if the coastal Bjerknes feedback discussed here plays a role, it must be much less. Actually, the correlation coefficient between January-March California Niño/Niña and alongshore wind (upwelling) indices as defined in this study is only about 20.27 (20.16).
Since the coastal Bjerknes feedback can operate theoretically in the subtropical coastal oceans in summer, we expect that similar coastal Niño/Niña phenomena may exist along the eastern boundaries of the subtropical oceans under subtropical highs in both Northern and Southern Hemispheres. In the austral summer of 2010/11, extreme coastal warming took place off the western coast of Australia and devastated dramatically the marine ecosystem. This anomalous warming was called Ningaloo Niño 29 . Although it was attributed mainly to the oceanic and atmospheric teleconnections related to the co-occurring La Niña in the tropical Pacific, it may have been amplified by the coastal Bjerknes feedback 18,29,30 .
The California Niño/Niña lasts more than one season and may certainly impact the local marine ecosystem. For example, the spawning and meridional migration of Pacific sardine in the California Current Ecosystem may be influenced by the related SST anomalies 31 , and its abundance may be impacted by the upwelling-induced anomalies in the biological production 32 . Also, the California Niño/Niña may increase/decrease the summer precipitation in Baja California and surface temperature in the nearby coastal continent (Supplementary Fig. 6). Hence, accurate seasonal prediction of the anomalous events could be beneficial to the industrial management and associated economic activities. Since the California Niño/Niña in boreal summer is less influenced by ENSO than that in boreal winter, its seasonal prediction is challenging. So far, the seasonal predictability of subtropical SST anomalies is mostly originated from their robust relationship with the tropical ENSO [33][34][35] . For instance, the 2010/11 Ningaloo Niño can be successfully predicted nine months ahead in a coupled ocean-atmosphere model by virtue of the co-occurring La Niña 35 . We believe that efforts to predict regional intrinsic climate modes such as California Niño/Niña in boreal summer may contribute to enhancing skills of seasonal prediction for various societal applications as well as to deepening our knowledge of regional climate variations in the extratropics.

Methods
We use the National Oceanic and Atmospheric Administration (NOAA) Optimum Interpolation Sea Surface Temperature (OISST) version 2, Pacific Fisheries Environmental Laboratory (PEFL) derived upwelling indices, and National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) reanalysis 1 for the period of January 1982 to December 2011. The California Niño/Niña index is defined as the monthly SST anomalies averaged over 110uW-120uW and 20uN-30uN (the enclosed coastal ocean in Fig. 1a). The Niño3 index is the monthly SST anomalies averaged over 90uW-150uW and 5uS-5uN. The JAS California Niño/Niña index with variations related to ENSO linearly regressed out (K) is computed by the formula: Here, K 1 is the JAS California Niño/Niña index, K 2 the JAS Niño3 index, r 12 the correlation coefficient between them, D 1 the standard deviation of K 1 , and D 2 the standard deviation of K 2 . The alongshore surface wind index is the alongshore component of the 10-meter-height wind anomalies averaged over the same ocean as the California Niño/Niña index after projecting the surface wind along the coast of Baja California (321u to the north). The PEFL upwelling indices are computed from monthly-mean pressure fields and reflect the wind-induced coastal upwelling and off-shore Ekman transport at 15 standard locations along the west coast of North America 21 . Among them, three locations are within the coastal ocean of our interest (24uN 113uW, 27uN 115uW and 30uN 119uW) and thus their monthly anomalies are averaged to represent the upwelling index in this study. Three-month-running mean anomalies are adopted to represent the monthly anomalies in Figs. 2-4 and Supplementary Figs. 3-6, 8-10 to minimize the intraseasonal variations, but similar results can be obtained if monthly data are used ( Supplementary Fig. 7). The linear correlation and regression analyses are applied in this study to demonstrate the existence of coastal Bjerknes feedback for simplicity. However, there may be asymmetry between California Niño and Niña. For instance, California Niña may have larger anomalous amplitudes and spatial extension than Niño and last for a longer time ( Supplementary Figs. 8-10). Also, the JAS California Niño/Niña indices show a negative skewness of 20.11/20.20 with/without the simultaneous variations linearly related to ENSO. The asymmetry may be partly due to the co-occurrence of negative phase Pacific Decadal Oscillation with California Niña (Supplementary Fig. 8b) and partly due to non-linear processes; this needs to be studied further. We note that the results shown in this study do not change quantitatively even if the NOAA Extended Reconstructed Sea Surface Temperature (ERSST) version 3b, NCEP-National Energy Research Supercomputing Center of the Department of Energy (DOE) reanalysis 2 and European Centre for Medium-Range Weather Forecast Interim Reanalysis (ERA-Int) are used.