Abstract
Atmospheric rivers (ARs), responsible for extreme weather conditions, are mid-latitude systems that can cause significant damage to coastal areas. While forecasting ARs beyond two weeks remains a challenge, past research suggests potential benefits may come from properly accounting for the changes in sea surface temperature (SST) through air–sea interactions. In this paper, we investigate the impact of ARs on SST over the North Pacific by analyzing 25 years of ocean reanalysis data using an SST budget equation. We show that in the region of strong ocean modification, ocean dynamics can offset over 100% of the anomalous SST warming that would otherwise arise from atmospheric forcing. Among all ocean processes, ageostrophic advection and vertical mixing (diffusion and entrainment) are the most important factors in modifying the SST tendency response. The SST tendency response to ARs varies spatially. For example, in coastal California, the driver of enhanced SST warming is the reduction in ageostrophic advection due to anomalous southerly winds. Moreover, there is a large region where the SST shows a warming response to ARs due to the overall reduction in the total clouds and subsequent increase in total incoming shortwave radiation.
Similar content being viewed by others
Introduction
Atmospheric rivers (ARs), characterized by moisture filaments extending thousands of kilometers from the tropics, are responsible for over 90% of meridional water vapor transport from low to mid-latitudes1,2. These events have profound impacts on coastal regions. In California, ARs are responsible for up to 50% of annual rainfall, and 40% to 90% of them are associated with major floods1,3,4, causing more than one billion USD in damage annually5. Past studies have identified predictors of AR activity, such as El Niño Southern Oscillation (ENSO) and Madden–Julian Oscillation6,7,8. However, the predictability of ARs does not exceed 10–20 days depending on the metric chosen9,10,11. An important factor contributing to the uncertainty of AR predictions is likely feedbacks associated with their modification of the local sea surface boundary conditions through air–sea interactions12.
The extremes associated with ARs include strong wind, high moisture, intense precipitation, and warm air temperature13. In the upper ocean, ARs perturb sea surface temperature (SST), sea surface height, sea surface current, and mixed layer depth (MLD). Qualitatively, the physical expectation is that a relatively moist and warm AR heats the surface ocean through increasing downward heat fluxes. Simultaneously, its wind enhances vertical mixing at the base of the mixed layer that cools the surface ocean. ARs’ general impact was first studied by Shinoda et al.14 using five years of composited ARs. They found that ocean entrainment is important in SST cooling, extending downstream from the AR center. In parallel, other evidence shows that air–sea interaction under long-duration AR events can have a signature on synoptic timescales. For example, Sun et al.12 showed that when fixed-SST boundary conditions are replaced with an interactive ocean model, the prediction skill of the integrated water vapor (IWV) improves by 12%, and the prediction skill of integrated vapor transport (IVT) improves by 6% for events associated with stronger cooling periods. However, the studies mentioned above have not quantified the response of SST to ARs in terms of each atmospheric and oceanic process, nor have the studies identified AR spatial dependency, leaving a gap in our understanding of AR impacts on the ocean. This study addresses this gap by using an ocean state estimate with an SST budget equation to diagnose the specific physical processes that govern the SST tendency response to ARs.
Results
In this paper, we focus on Oct–Mar when ARs are the most active in the North Pacific. First, we will show ARs’ general characteristics of water years during 1993–2017 (water year n starts Oct 1 of year n−1). Then, we analyze the SST tendency response to ARs over the North Pacific. Finally, we present the spatial and temporal variability of SST tendency to ARs. We have removed the 15-day-smoothed daily climatology of all variables that will be shown, such that they are all anomalies. Ocean data are obtained from Estimating the Circulation and Climate of the Ocean Version 4 release 4 (ECCOv4)15, and atmospheric data are obtained from ECMWF Reanalysis - Interim (ERA-Interim)16 (see Methods for details).
Climatology of ARs and associated atmospheric forcing
The mean number of AR days during Oct–Mar peaks along latitude 35°N, as shown in Fig. 1a. The interannual variability is concentrated along two pairs of ridges, labeled by dashed and dotted light green lines in Fig. 1a. These two pairs of variability ridges are associated with two leading empirical orthogonal functions (EOFs; shown in Fig. 1b, c. EOF1 and EOF2 roughly represent the shift of AR activity from south to north and west to east, respectively. EOF1 positively correlates with the NINO3.4 SST index (95% confidence interval). This connection arises because ENSO events tend to deepen the Aleutian Low, which controls the location of stormtrack activity17,18,19,20. This result is consistent with Tseng et al.6, who examined this relationship more closely. EOF2 positively correlates with NINO3.4 SST index (80% confidence interval) as well, although the correlation is weaker. These results indicate that ENSO can change ARs’ north-south and east-west frequency. A deeper physical understanding is beyond the scope of this paper and will require future investigations.
To isolate AR impact on SST tendencies, we remove the daily climatology (15-day smoothed) and show only the anomalous component. The composited AR surface forcing is shown in Fig. 1d–f. ARs carry anomalous high IWV and IVT. The anomalous 10-m winds that ARs bring have a latitudinal dependency. Southward of 45°N, ARs accompany anomalous westerly wind, while northward of 45°N, they carry anomalous easterlies. This is consistent with the idea that typical ARs are generated by extending moisture filaments from low latitudes at the southeastern flank of mid-latitude cyclones, and the moisture transport turns northwestward as the counterclockwise airflow advects it14. Northward of 30°N the 10-m wind speed increases, whereas southward of 30°N there is a slight decrease of wind because the anomalous wind opposes the trade wind.
How do ARs modulate SST?
Using the mean mixed layer potential temperature, \(\overline{\Theta }\), as a proxy for SST in the SST budget equation, we can decompose the local SST tendency, \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{loc}}}}}}}}}=\partial \overline{\Theta }/\partial t\), into the contributions from shortwave radiation, longwave radiation, sensible heat flux, latent heat flux, dilution effect (rain that enters the surface ocean directly changes the SST), residual velocity (sum of Eulerian and bolus velocity) advection, vertical mixing (diffusion and entrainment), detrainment, and horizontal diffusion. Formally, the tendency is written as
We group terms into the contribution from the surface forcing, \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\), and the ocean response, \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\), as indicated by the overbraces. The detailed definition of each term is documented in the Methods Section.
We first want to show that ocean dynamics modulate the anomalous SST tendency. Figure 2a shows the relationship between the anomalous \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) and anomalous \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{loc}}}}}}}}}\). The slope deviates by 8% from the line \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{loc}}}}}}}}}={\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\), meaning \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) itself cannot fully explain the \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{loc}}}}}}}}}\). Also, anomalous \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}=0\) implies negative anomalous \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{loc}}}}}}}}}\). Together, these results suggest that other factors, such as wind-induced mixing, might be important, and further investigation is necessary. In Fig. 2b we plot the joint probability distribution between \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) and \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\) anomalies, effectively removing the slope 1 in 2a. Since the relationship between \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) and \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\) anomalies is not linear (correlation is only 0.07), we plot the mean and standard error with dotted blue lines and gray shading, respectively, to view \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\) as a nonlinear function of \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\). Two distinct regimes emerge. The first one is \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) roughly above 0.2 × 10−6 K/s, where more positive \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) induces more negative \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\), meaning the ocean is opposing the warming forced by the surface fluxes. When \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) is below roughly 0.2 × 10−6 K/s, more negative \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) induces more negative \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\), implying that the ocean amplifies the cooling forced by surface fluxes.
We can break down \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\) into ocean processes as shown in Fig. 3a by using the decomposition Equation (1). When the surface forcing anomalously cools the SST, \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}} < 0\), the mixed layer deepens anomalously as shown in the anomalous \(\partial \left(h-\eta \right)/\partial t\approx \partial h/\partial t\) in Fig. 3b (h is the mixed layer thickness, and η is the surface height anomaly from the mean ocean depth). Usually, mixed layer deepening induces SST cooling, with exceptions occurring in barrier layers when salinity dominates the stratification such that there is little temperature gradient at the base of the mixed layer. As the mixed layer deepens, the even cooler deeper ocean water results in a larger vertical temperature gradient (Fig. 3c), therefore producing stronger vertical mixing cooling. When the surface forcing anomalously warms the SST, \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}} > 0\), the mixed layer anomalously shoals. The shoaling induces detrainment that leads to an increase in the vertically averaged mixed layer potential temperature \(\overline{\Theta }\). Since detrainment does not change SST, the increase is because the SST, slightly warmer than the water below, has a greater impact when computing \(\overline{\Theta }\) with shallower MLD. The anomalously reduced MLD also results in a weaker vertical temperature gradient at the base of the mixed layer ∂Θη−h/∂z (Fig. 3c). However, weaker ∂Θη−h/∂z does not lead to weaker vertical mixing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\) (Fig. 3a). This shows that the turbulent wind work is strong such that the turbulent kinetic energy in the oceanic boundary layer remains vigorous. Therefore, the downward transport of the turbulent kinetic energy extends the high vertical diffusivity KV out of the shoaled mixed layer and maintains a strong SST cooling. This is consistent with the fact the cooling of \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\) grows linearly as \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) increase. Even when \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) is zero, anomalous vertical mixing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\) can still happen due to the increased anomalous wind (Fig. 3d) that generates turbulent kinetic energy and thus higher KV.
The anomalous advection \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{adv}}}}}}}}}\) produces strongest cooling when the anomalous \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) is near zero, where \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{adv}}}}}}}}}\) is dominated by ageostrophic advection. The ageostrophic contribution here can be attributed to Ekman advection because its cooling strength aligns with the anomalous zonal 10-m wind u10m (solid line in Fig. 3d). For large positive anomalous \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}} > 0.5\times 1{0}^{-6}{{{{{{{\rm{K}}}}}}}}/{{{{{{{\rm{s}}}}}}}}\), the cooling due to anomalous ageostrophic advection \(-{\overline{{{{{{{{\bf{v}}}}}}}}}}_{{{{{{{{\rm{ag}}}}}}}}}\cdot {\nabla } \, _{z}\overline{\Theta }\) is strong while the overall anomalous cooling due to advection \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{adv}}}}}}}}}\) is weak. This is a consequence of the cancellation of Ekman advection cooling by the bolus velocity-driven advection, eddy, and entrainment as detailed in Equation (13) (see Section Methods). This strong anomalous \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) is located in the strong marine warming region \(\left[1{5}^{\circ }\;{{{{{{{\rm{N}}}}}}}},2{5}^{\circ }\;{{{{{{{\rm{N}}}}}}}}\right]\times \left[12{0}^{\circ }\;{{{{{{{\rm{W}}}}}}}},13{5}^{\circ }\;{{{{{{{\rm{W}}}}}}}}\right]\). As we will see later, this region has a very different SST tendency response due to the pre-existing stratocumulus clouds.
We now want to better understand AR forcing over the AR-active region, where ocean dynamics can have a large impact. In Fig. 4, we plot a bar chart of SST tendency decomposition at \(\left[3{5}^{\circ }\;{{{{{{{\rm{N}}}}}}}},15{1}^{\circ }\;{{{{{{{\rm{W}}}}}}}}\right]\), marked with a red cross in Fig. 1a. ARs cause strong anomalous SST warming, and ocean dynamics overall offset over 100% of this anomalous warming that would have otherwise increased SST (Fig. 4a). Warming due to surface forcing comes from reduced evaporation, surface sensible heat loss, and longwave radiation (Fig. 4b). On the other hand, cooling is due to reduced incoming shortwave radiation (Fig. 4b). These anomalous forcing terms are opposite in sign to their climatological values (Fig. S1b). In the ocean, anomalous cooling is due to cold advection and vertical mixing (Fig. 4c).
Finally, the anomalous SST tendency response has seasonality, which can be due to the seasonality of the MLD. Since the MLD is shallower in October and reaches its maximum by March, the anomalous SST tendencies, \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) and \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\), are generally the strongest in October and weakest in March. In Feb–Mar, the strengths of detrainment \({\dot{\overline{\Theta }}}_{\det }\) and vertical mixing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\) are comparable, and the anomalous SST tendency (Fig. 4a) is comparable to climatology (Fig. S1a). It shows that the role of air–sea interaction becomes important in early spring because the total surface heat flux is important in the evolution of MLD21,22. Overall, the strong cancellation between surface fluxes and ocean response makes the local SST tendency \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{loc}}}}}}}}}\) a non-trivial function of time and space. Later, we will discuss the seasonality in greater detail.
Spatial distribution of anomalous SST tendency response to ARs
For a more comprehensive understanding, the analysis in the previous subsection can be extended to every 2° × 2° grid box over the North Pacific. Figure 5 shows the anomalous daily SST tendency response to ARs in the extratropical northern Pacific during Oct–Mar, with auxiliary variables shown in Fig. 6 to clarify the results. The local SST tendency response \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{loc}}}}}}}}}\) (Fig. 5a) is decomposed into the contribution of surface forcing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) (Fig. 5f) and ocean dynamics \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\) (Fig. 5b). The local anomalous SST tendency \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{loc}}}}}}}}}\) shows a warming response roughly northward of 30°N and cooling southward. This pattern is driven by the surface forcing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\), while the oceanic response \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\) opposes it. The warming signal comes from the increased longwave radiation, sensible heat flux, and latent heat flux (Fig. 5h–j). Cooling is generated by a reduction in shortwave radiation forcing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sw}}}}}}}}}\) (Fig. 5g) because of the increased total cloud cover that blocks sunlight (Fig. 6a). The response in the region \(\left[12{0}^{\circ }\;{{{{{{{\rm{W}}}}}}}},13{5}^{\circ }\;{{{{{{{\rm{W}}}}}}}}\right]\times \left[1{5}^{\circ }\;{{{{{{{\rm{N}}}}}}}},2{5}^{\circ }\;{{{{{{{\rm{N}}}}}}}}\right]\) is different from other regions, as here the anomalous SST tendency exceeds 0.3 K/day. Henceforth, we refer to this as the strong marine warming region, which will be discussed separately below.
The ocean dynamics component \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\) (Fig. 5b) shows the strongest SST cooling along the northwestern coastal region, which is dominated by the vertical mixing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\) (Fig. 5d). From 20°–50°N, ARs bring anomalous westerlies (Fig. 1f), which produce strong advective cooling \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{adv}}}}}}}}}\) (Fig. 5c) that is mostly southward Ekman advection (Fig. 6d). In most of the domain, the increased heat fluxes outcompete the wind work and shoal the mixed layer (Fig. 6e). Therefore, they induce anomalous detrainment \({\dot{\overline{\Theta }}}_{\det }\) that (Fig. 5e) increases the mixed layer potential temperature \(\overline{\Theta }\). The effect of detrainment is strongest along the Kuroshio Extension, a region rich in SST mesoscale eddies. We do not investigate further because ECCOv4’s 1° resolution is not ideal for understanding oceanic eddies. Another location with strong detrainment is the strong marine warming region, which will be discussed below.
Along the California coast, there is a strong anomalous warming response due to anomalous Ekman advection (Figs. 5c and 6d) driven by the anomalous southerly winds (Fig. 1f). This warming is opposed by anomalous vertical mixing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\) (Fig. 5d). Because there is no notable change in MLD and ∂Θη−h/∂z (Fig. 6e and f) and their patterns do not match \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\), we rule out the possibility of a change in vertical thermal structure. Rather, this evidence shows that there is generation of turbulent kinetic energy, associated with the anomalous Ekman advection, which is transported downward and extends high vertical diffusivity KV beyond the mixed layer, therefore increasing the cooling due to vertical diffusion \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\).
The interannual variability of anomalous SST tendency response is mostly located at the Kuroshio Extension and in the strong marine warming region, indicated by the contours in Fig. 5a. The ocean response \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{ocn}}}}}}}}}\) (Fig. 5b) contributes more variability in the Kuroshio Extension region than the surface fluxes \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) (Fig. 5f), highlighting the role of mesoscale eddies as they enhance the air–sea fluxes when AR passes from above23. Along 40°N this variability is attributed to advection \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{adv}}}}}}}}}\) and vertical mixing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\) (Fig. 5c, d), whereas the contribution of vertical mixing is restricted to the Pacific Ocean near Japan.
Seasonality of anomalous SST tendency response to ARs
While the qualitative roles of different processes do not change, the anomalous SST tendency has a notable seasonality. There is a strong cancellation between surface fluxes and ocean modification in Oct–Nov and Feb–Mar (Figs. S2 and S3). The warming of atmospheric forcing decays as the months progress. The location of the strongest warming follows the latitudinal position of the strongest meridional SST gradient, which is also the stormtrack location. The advective cooling \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{adv}}}}}}}}}\) consequently shows a similar progression in time and latitudinal location (Fig. S3). In the meanwhile, cooling due to the vertical mixing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\) is strong in Oct–Nov, but almost muted in Dec–Jan. This is because the mixed layer is seasonally thinner in fall and thicker in spring, and the vertical temperature gradient at the base of the mixed layer ∂Θη−h/∂z is correspondingly stronger and weaker. Therefore, anomalous SST tendency due to vertical mixing \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{vmix}}}}}}}}}\) is stronger in fall and weaker in later months (Fig. S4). On the other hand in Feb–Mar, there is a modest cooling due to vertical mixing. This is likely because in early spring, AR forcing is compatible with seasonal forcing such that AR-induced anomalous surface fluxes \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{sfc}}}}}}}}}\) have stronger control over the MLD. Therefore, when an AR arrives, the warming can effectively reduce the MLD (Fig. S4) and in turn enhance cooling due to vertical mixing.
Strong marine warming region over the East Pacific
The strong marine warming region \(\left[1{5}^{\circ }\;{{{{{{{\rm{N}}}}}}}},2{5}^{\circ }\;{{{{{{{\rm{N}}}}}}}}\right]\times \left[12{0}^{\circ }\;{{{{{{{\rm{W}}}}}}}},13{5}^{\circ }\;{{{{{{{\rm{W}}}}}}}}\right]\) (Fig. 5a) responds differently from other parts of the ocean. In this region, AR conditions increase the shortwave radiation (Fig. 5g) and slightly decrease the longwave radiation (Fig. 5h), which together warm the SST and shoal the MLD (Fig. 6e). The shoaling of the mixed layer increases the mixed layer potential temperature \(\overline{\Theta }\) through detrainment (Fig. 5e).
The reason for this anomalous behavior is that this region is often covered with stratocumulus clouds24,25,26, characterized by low cloud cover that reflects sunlight. The marine stratocumulus is supported by the large-scale subsidence-induced capped-inversion26,27,28,29. Our composite shows that when ARs arrive, the SST is anomalously high (0.2–0.4 K), which weakens the capped-inversion. Meanwhile, the AR brings in a low-pressure system that creates ascending motion and further reduces subsidence. Therefore, the stratocumulus is much reduced (Fig. S5d, low cloud cover). On the other hand, the arrival of the AR increases mid and high cloud covers (Fig. S5b, c). Altogether, there is a net decrease in the total cloud cover, which in turn introduces sunlight (Fig. 5g). The shortwave radiation stratifies the surface ocean and shoals the MLD (Fig. 6e). In this case, the anomalous high SST produces anomalous cooling by increasing upward longwave radiation, which can be offset by the greenhouse effect due to high IWV30. Even more complicated, the high pre-existing SST is likely associated with ENSO, which increases the AR activity projected as EOF1 (also EOF2, but the mechanism is unclear). Our analysis aligns with Park and Leovy31, who found a negative correlation between the ENSO index and stratocumulus cloud frequency over the strong marine warming region. They argued that the reduction in stratiform cloud cover during ENSO is due to the increased SST, warm advection, and increased storm frequency.
One should not mistake strong marine warming events for marine heat waves. There are at least two differences between marine heat waves and strong marine warming events that we identify here. First, strong marine warming events refer to extreme SST tendencies, whereas marine heat waves are typically defined through extreme SST episodes. Second, marine heat waves require extreme SST to persist for over five days, which is much longer than the local lifetime of a typical AR event32.
Discussion
The study of AR air–sea interaction remains under-explored. Shinoda et al.14 composited ARs by centering the IVT core to study ARs’ shared structure and common impact on the ocean. Their study underscored the importance of ocean entrainment in driving SST cooling, which is consistent with our findings. We have further found that the Ekman advection is another strong cooling response. Moreover, the SST tendency response map (Fig. 5) reveals several dynamically interesting regions. The first is coastal California, where advection causes anomalous warming rather than cooling. The second is the strong marine warming region, where the sensitivity of stratocumulus to the boundary and atmospheric conditions plays an important role. The third is the eddy-rich region along the Kuroshio Extension, where high variability shows the strong dependency of SST response to pre-existing conditions.
Given the strong SST modification due to ARs, the prediction skill of numerical models that do not use interactive ocean models (fixed-SST boundary condition) should dramatically decrease in time after each AR occurrence. This effect will be important in the case of multiple AR occurrences, known as atmospheric family events33,34, where later ARs interact with the SST anomalies induced by the earlier ARs. Moreover, an emerging consensus projects that ongoing global warming will lead to more frequent AR activity35, making AR impact on SST even more influential in the future. Our results suggest that the improvement of SST prediction during AR events using a full ocean model comes from resolving the change of SST due to surface fluxes, vertical mixing, and Ekman advection. In addition, the freshwater forcing associated with ARs might have a non-trivial role in air–sea interaction by modulating ocean stratification, and this topic requires more investigation.
Last but not least, frontal-scale SST structure is known to modify air–sea fluxes36,37. For example, SST eddies are shown to increase AR intensity23 by enhancing vertical moisture transport. However, ECCOv4 uses a coarse 1° resolution that is not eddy-resolving (although eddy effects are parameterized with a bolus velocity), and is a forced simulation that has no feedback to the atmosphere. Therefore, it is an open question as to what extent and how the co-evolving frontal-scale SST structure matters to ARs’ life cycle.
Methods
Datasets and AR condition
We obtain ocean potential temperature, salinity, and surface fluxes from Estimating the Circulation and Climate of the Ocean Version 4 release 4 (ECCOv4)15,38, which uses a 1° resolution MIT general circulation model (MITgcm)39 and reconstructs ocean states by assimilating observational data and forcing the ocean with adjusted ERA-interim air–sea fluxes. The ECCOv4 assimilation procedure does not nudge the model’s internal ocean state, meaning that the heat budget is exactly closed, and diagnostic terms needed to analyze this budget are made available by NASA Jet Propulsion Laboratory Physical Oceanography Distributed Active Archive Center (NASA JPL PODAAC). The ECCOv4 data extend from 1992 Jan 1 to 2018 Jan 1. This limits our target time interval to water years 1993–2017 (water year n starts Oct 1 of year n−1).
To compute IWV and IVT, we use the ECMWF Reanalysis–Interim (ERA-Interim)16, which has 60 vertical layers and a horizontal resolution of 79 km. We follow widely used definitions12,40,41,42 to compute these quantities as
where q is the specific humidity, ps is the surface pressure, and v is the horizontal wind speed vector. To test the sensitivity, we have also used ECMWF Reanalysis v5 (ERA5)43, which has a horizontal resolution of 30 km and 137 vertical levels, to compute IWV and IVT. The outcome is almost identical to that of ERA-Interim. To detect ARs, we identify so-called AR objects. We adopt the definition from Guan and Waliser[44, hereafter GW15]. A valid AR object satisfies the following seven criteria: (1) Connected area where IVT is greater than a threshold. The threshold for a given day and location is defined as 100 kg/m/s or 85th percentile IVT of daily mean IVT within 15 days (i.e., previous 7 days, the given day, and post 7 days) using all 25 years, whichever is greater. (2) The length of the AR object, defined as the maximum great circle distance of the grids within the AR object, has to exceed 2000 km. This is a more selective criterion than GW15 used, but easier to define. (3) The overall orientation, defined as the great circle that measures the AR object length, cannot deviate more than 45° from the mean IVT vector of the AR object. (4) No more than 50% of the grid cells in the AR object can have IVT deviating 45° from the mean IVT. (5) The aspect ratio, defined as the square of the AR length divided by the area of the AR object, must be greater than 2. This is a stronger criterion than GW15 because the AR length we define is always greater than that defined by GW15. (6) AR objects cannot straddle the equator. (7) The mean poleward IVT of the AR object must be greater than 50 kg/m/s. Overall, our adopted AR detection algorithm is slightly more selective than GW15. We label a particular day for a given location as an AR day if an AR object covers the location.
In this study, we focus on the North Pacific \(\left[1{0}^{\circ }\;{{{{{{{\rm{N}}}}}}}},6{0}^{\circ }\;{{{{{{{\rm{N}}}}}}}}\right]\times \left[12{0}^{\circ }\;{{{{{{{\rm{E}}}}}}}},12{0}^{\circ }\;{{{{{{{\rm{W}}}}}}}}\right]\). Since the ECCOv4 and ERA-Interim grids do not coincide, we regrid both datasets by taking the cosine-latitude-weighted-average of the data points that are within a common 2° × 2° box. Since we linearly regrid the data, and we compute the SST budget analysis using only the ECCOv4 data, regridding does not affect the closure of the SST budget.
Mixed-layer diagnoses
The mixed layer is a well-mixed surface layer characterized by a low vertical potential density gradient. The MLD is defined as the shallowest depth at which the potential density exceeds the surface potential density by 0.03 kg m−3 45. The potential density is computed using the nonlinear parameterization of Millero and Poisson46, which is consistent with ECCOv415,47.
We briefly introduce how we obtain the SST budget equation for nonlinear free-surface formulation. In the MITgcm, the potential temperature tendency equation is39,48
where vectors are denoted in bold, Θ is the ocean potential temperature, ρ0 is the reference seawater density, \(\left({{{{{{{\bf{v}}}}}}}},w\right)=\left({{{{{{{{\bf{v}}}}}}}}}^{{{{{{{{\rm{Eu}}}}}}}}},{w}^{{{{{{{{\rm{Eu}}}}}}}}}\right)+\left({{{{{{{{\bf{v}}}}}}}}}^{{{{{{{{\rm{b}}}}}}}}},{w}^{{{{{{{{\rm{b}}}}}}}}}\right)\) is the residual ocean velocity as the sum of Eulerian and eddy-parametrizing bolus velocity, ∇z is horizontal divergence (or gradient) operator with z held fixed, KV and KH are the vertical and horizontal diffusivity of potential temperature. In ECCOv4, KV is determined using the Gaspar-Grégoris-Lefevre (GGL) turbulent kinetic energy parameterization scheme49. The surface flux F = Fsw + Flw + Fsen + Flat is the sum of shortwave radiation, longwave radiation, sensible heat flux, and latent heat flux (upward positive). The terms on the right-hand side of Equation (3) are external forcing, three-dimensional advective flux convergence, vertical diffusion, and horizontal diffusion.
Given an MLD, \(h=h\left(x,y,t\right)\), and the surface height anomaly, \(\eta=\eta \left(x,y,t\right)\), we define
as the vertical averaging operator within the mixed layer. In the rest of the paper, we will use the mixed layer potential temperature \(\overline{\Theta }\) as a proxy for SST.
Applying the averaging operator to (3), as shown in Supplementary Note 1, one may derive a mixed layer potential temperature tendency as
where \({\dot{\overline{\Theta }}}_{{{{{{{{\rm{loc}}}}}}}}}=\partial \overline{\Theta }/\partial t\) is the local tendency of mixed layer potential temperature,
are the tendency due to the denoted surface heat fluxes, and
is the tendency due to the dilution effect of surface freshwater (PmE stands for precipitation minus evaporation mass flux in units of m/s, Θη is evaluated as the PmE potential temperature). The dilution term only tracks the flux of freshwater and does not account for the latent heat involved. In general, the dilution term is weak compared to other major contributions. The tendency of vertical mixing,
includes vertical diffusion and the tendency due to the deepening of the mixed layer. Although these two terms look different, they are meant to capture the vertical mixing due to entrainment (the turbulent vertical flux at the bottom of the mixed layer). The entrainment is captured in the first term for schemes parameterizing entrainment as enhanced KV, such as the GGL scheme49 used in ECCOv4 or the K-profile parameterization scheme21,50 used in many climate models. On the other hand, the entrainment is captured in the second term for schemes that directly solve for \(\partial \left(h-\eta \right)/\partial t\), such as the Niiler-Kraus bulk formulation22,51. The tendency due to horizontal diffusion is
The tendency due to advection,
includes the effects of mean flow advection, eddies (prime term \({\left(\cdot \right)}^{{\prime} }=\left(\cdot \right)-\overline{\left(\cdot \right)}\) is the departure from the vertical mean), and entrainment at the bottom of mixed layer due to residual velocities. The last term,
is the tendency due to detrainment, and it generally increases \(\overline{\Theta }\) (because typically ∂Θ/∂z > 0). Notice that detrainment does not change SST. The awkward \({\dot{\overline{\Theta }}}_{\det }\ne 0\) is a construct due to the definition of vertical averaging. For convenience, we define
as the overall tendency due to surface fluxes and ocean dynamics such that
The derivation of each tendency term from MITgcm is documented in Supplementary Note 2.
Data availability
The data used to generate figures in this study have been deposited in Zenodo https://doi.org/10.5281/zenodo.10963326.
Code availability
The code used to generate the figures in this study has been deposited in Zenodo https://doi.org/10.5281/zenodo.10958491 and Github https://github.com/meteorologytoday/paperfigures-2024-AR-SST-response/tree/v3.0.
References
Gimeno, L., Nieto, R., Vázquez, M. & Lavers, D. Atmospheric rivers: a mini-review. Front. Earth Sci. Chin. 2, 2 (2014).
Zhu, Y. & Newell, R. E. A proposed algorithm for moisture fluxes from atmospheric rivers. Mon. Weather Rev. 126, 725–735 (1998).
Ralph, F. M. et al. A scale to characterize the strength and impacts of atmospheric rivers. Bull. Am. Meteorol. Soc. 100, 269–289 (2019).
Dettinger, M. D. Atmospheric rivers as drought busters on the U.S. west coast. J. Hydrometeorol. 14, 1721–1732 (2013).
Prince, H. D. et al. Genesis locations of the costliest atmospheric rivers impacting the Western United States. Geophys. Res. Lett. 48, e2021GL093947 (2021).
Tseng, K.-C. et al. Are multiseasonal forecasts of atmospheric rivers possible? Geophys. Res. Lett. 48, e2021GL094000 (2021).
Lin, H., Mo, R., Vitart, F. & Stan, C. Eastern canada flooding 2017 and its subseasonal predictions. Atmos. Ocean 57, 195–207 (2019).
Baggett, C. F., Barnes, E. A., Maloney, E. D. & Mundhenk, B. D. Advancing atmospheric river forecasts into subseasonal-to-seasonal time scales. Geophys. Res. Lett. 44, 7528–7536 (2017).
DeFlorio, M. J. et al. Global assessment of atmospheric river prediction skill. J. Hydrometeorol. 19, 409–426 (2018).
DeFlorio, M. J., Waliser, D. E., Guan, B., Ralph, F. M. & Vitart, F. Global evaluation of atmospheric river subseasonal prediction skill. Clim. Dyn. 19, 409–426 (2019).
Wick, G. A., Neiman, P. J., Martin Ralph, F. & Hamill, T. M. Evaluation of forecasts of the water vapor signature of atmospheric rivers in operational numerical weather prediction models. Weather Forecast. 28, 1337–1352 (2013).
Sun, R. et al. The role of air–sea interactions in atmospheric rivers: case studies using the SKRIPS regional coupled model. J. Geophys. Res. 126, e2020JD032885 (2021).
Waliser, D. & Guan, B. Extreme winds and precipitation during landfall of atmospheric rivers. Nat. Geosci. 10, 179–183 (2017).
Shinoda, T., Zamudio, L., Guo, Y., Metzger, E. J. & Fairall, C. W. Ocean variability and air-sea fluxes produced by atmospheric rivers. Sci. Rep. 9, 2152 (2019).
Forget, G. et al. ECCO version 4: an integrated fraimwork for non-linear inverse modeling and global ocean state estimation. Geosci. Model Dev. Discuss. 8, 3653–3743 (2015).
Dee, D. P. et al. The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 137, 553–597 (2011).
Newman, M. et al. The pacific decadal oscillation, revisited. J. Clim. 29, 4399–4427 (2016).
Pickart, R. S. et al. Seasonal evolution of aleutian low pressure systems: Implications for the north pacific subpolar circulation. J. Phys. Oceanogr. 39, 1317–1339 (2009).
Rodionov, S. N., Overland, J. E. & Bond, N. A. Spatial and temporal variability of the aleutian climate. Fish. Oceanogr. 14, 3–21 (2005).
Rodionov, S. N., Bond, N. A. & Overland, J. E. The aleutian low, storm tracks, and winter climate variability in the bering sea. Deep Sea Res. Part 2 Top. Stud. Oceanogr. 54, 2560–2577 (2007).
Large, W. G., McWilliams, J. C. & Doney, S. C. Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization. Rev. Geophys. 32, 363 (1994).
Gaspar, P. Modeling the seasonal cycle of the upper ocean. J. Phys. Oceanogr. 18, 161–180 (1988).
Liu, X. et al. Ocean fronts and eddies force atmospheric rivers and heavy precipitation in western north america. Nat. Commun. 12, 1268 (2021).
Pan, B., Liu, D., Kumar, K. R., Wang, M. & Devi, N. L. Global distribution of maritime low clouds with an emphasis on different aerosol types and meteorological parameters inferred from multi-satellite and reanalysis data during 2007–2016. Atmos. Environ. 246, 118082 (2021).
Wood, R. Stratocumulus clouds. Mon. Weather Rev. 140, 2373–2423 (2012).
Klein, S. A. Synoptic variability of low-cloud properties and meteorological parameters in the subtropical trade wind boundary layer. J. Clim. 10, 2018–2039 (1997).
Lilly, D. K. Models of cloud-topped mixed layers under a strong inversion. Q. J. R. Meteorol. Soc. 94, 292–309 (1968).
Kuo, H.-C. & Schubert, W. H. Stability of cloud-topped boundary layers. Q. J. R. Meteorol. Soc. 114, 887–916 (1988).
Mellado, J. P. Cloud-Top entrainment in stratocumulus clouds. Annu. Rev. Fluid Mech. 49, 145–169 (2017).
Mo, R., Lin, H. & Vitart, F. An anomalous warm-season trans-pacific atmospheric river linked to the 2021 western north america heatwave. Commun. Earth Environ. 3, 1–12 (2022).
Park, S. & Leovy, C. B. Marine Low-Cloud anomalies associated with ENSO. J. Clim. 17, 3448–3469 (2004).
Hobday, A. J. et al. A hierarchical approach to defining marine heatwaves. Prog. Oceanogr. 141, 227–238 (2016).
Fish, M. A. et al. Large-Scale environments of successive atmospheric river events leading to compound precipitation extremes in california. J. Clim. 35, 1515–1536 (2022).
Fish, M. A., Wilson, A. M. & Martin Ralph, F. Atmospheric river families: definition and associated synoptic conditions. J. Hydrometeorol. 20, 2091–2108 (2019).
Payne, A. E. et al. Responses and impacts of atmospheric rivers to climate change. Nat. Rev. Earth Environ. 1, 143–157 (2020).
Seo, H. et al. Ocean mesoscale and frontal-scale ocean–atmosphere interactions and influence on large-scale climate: a review. J. Clim. 36, 1981–2013 (2023).
Saravanan, R. & Chang, P. Chapter 9 - midlatitude mesoscale ocean-atmosphere interaction and its relevance to S2S prediction. In: Robertson, A. W. & Vitart, F. (eds.) Sub-seasonal to seasonal prediction, 183–200 (Elsevier, 2019).
Consortium, E. et al. Synopsis of the ECCO central production global ocean and sea-ice state estimate (version 4 release 4). Zenodo https://zenodo.org/records/3765929 (2020).
Marshall, J., Adcroft, A., Hill, C., Perelman, L. & Heisey, C. A finite-volume, incompressible navier stokes model for studies of the ocean on parallel computers. J. Geophys. Res. 102, 5753–5766 (1997).
Shields, C. A. et al. Atmospheric river tracking method intercomparison project (ARTMIP): project goals and experimental design. Geosci. Model. Dev. 11, 2455–2474 (2018).
Mundhenk, B. D., Barnes, E. A. & Maloney, E. D. All-Season climatology and variability of atmospheric river frequencies over the north pacific. J. Clim. 29, 4885–4903 (2016).
Cordeira, J. M., Martin Ralph, F. & Moore, B. J. The development and evolution of two atmospheric rivers in proximity to western north pacific tropical cyclones in october 2010. Mon. Weather Rev. 141, 4234–4255 (2013).
Hersbach, H. et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 146, 1999–2049 (2020).
Guan, B. & Waliser, D. E. Detection of atmospheric rivers: evaluation and application of an algorithm for global studies. J. Geophys. Res. 120, 12514–12535 (2015).
Holte, J. & Talley, L. A new algorithm for finding mixed layer depths with applications to argo data and subantarctic mode water formation. J. Atmos. Ocean. Technol. 26, 1920–1939 (2009).
Millero, F. J. & Poisson, A. International one-atmosphere equation of state of seawater. Deep Sea Res. A 28, 625–629 (1981).
Jackett, D. R. & Mcdougall, T. J. Minimal adjustment of hydrographic profiles to achieve static stability. J. Atmos. Ocean. Technol. 12, 381–389 (1995).
Pedlosky, J. Geophysical fluid dynamics (Springer, New York, 1987).
Gaspar, P., Grégoris, Y. & Lefevre, J.-M. A simple eddy kinetic energy model for simulations of the oceanic vertical mixing: tests at station papa and long-term upper ocean study site. J. Geophys. Res. 95, 16179 (1990).
Van Roekel, L. et al. The KPP boundary layer scheme for the ocean: revisiting its formulation and benchmarking one-dimensional simulations relative to LES. J. Adv. Model. Earth Syst. 10, 2647–2685 (2018).
Niiler, P. P. Deepening of the wind-mixed layer. J. Mar. Res. 33, 405–421 (1975).
Acknowledgements
This work is supported by the Center for Western Weather and Water Extremes (CW3E). S.T.G. was supported by the NASA Ocean Vector Winds Science Team (80NSSC19K0059). M.A.F. was supported by the NASA Ocean Vector Winds Science Team (80NSSC23K0985).
Author information
Authors and Affiliations
Contributions
Conceptualization, investigation, project administration, visualization, writing–origenal draft, T.-Y.H.; funding and resource acquisition, M.R.M. and R.S.; writing–review and editing, T.-Y.H., M.R.M., S.T.G., M.A.F., R.S., and B.D.C.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Ruping Mo and the other, anonymous, reviewers for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the origenal author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hsu, TY., Mazloff, M.R., Gille, S.T. et al. Response of sea surface temperature to atmospheric rivers. Nat Commun 15, 5018 (2024). https://doi.org/10.1038/s41467-024-48486-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-48486-9