1. Introduction
Extratropical transient eddies largely regulate daily weather variations in midlatitudes, while interacting with the climatological-mean westerlies and low-frequency quasi-stationary circulation anomalies. These eddies play an important role in maintaining the time mean temperature field and hemispheric energy budget. Improved mechanistic understanding of these eddies is thus important for climate dynamics and regional weather forecasts. These transient eddies successively move eastward through a given midlatitude location, accounting locally for a major fraction of high-frequency fluctuations of a given atmospheric variable. Using a 10-yr atmospheric analysis data over the extratropical Northern Hemisphere, Blackmon (1976) decomposed the temporal variance of 500-hPa geopotential height locally into the contributions from high-pass-, bandpass-, and low-pass-filtered components, whose periods are less than 2 days, 2.5–6 days, and longer than 10 days, respectively. He suggested that maritime regions of large variance of the bandpass-filtered fluctuations closely correspond to those of frequent cyclone passage over the North Atlantic and Pacific basins, which were thus referred to as “storm tracks.” Blackmon et al. (1977) indicated that those storm tracks are characterized by strong lower-tropospheric poleward eddy heat flux based on the bandpass-filtered fluctuations, which is now commonly used as a measure of baroclinic eddy activity. In fact, poleward heat transport associated with high-frequency (sub-weekly) eddies are pronounced along the major storm tracks over the North Pacific (NP), North Atlantic, and south Indian Ocean. They tend to coincide closely with lower-tropospheric eddy-driven westerly jets and regions of strong near-surface baroclinicity anchored along major oceanic frontal zones (Nakamura et al. 2004).
Many studies have investigated climatological-mean, interannual or longer-scale variability, and seasonal evolution of storm-track activity based on similar techniques based on Eulerian statistics (e.g., Chang et al. 2002). Taking advantage of the Eulerian statistics, Lee and Kim (2003) and Nakamura et al. (2004) clarified the maintenance mechanisms for the westerly jets and storm tracks by separating eddy-driven, subpolar jets from the thermally driven subtropical jets. As argued by Nakamura et al. (2004), this clear separation holds over the North Atlantic throughout the cold season, while it holds over the NP only in the shoulder seasons (i.e., autumn and spring). Over the midwinter NP, a stronger subtropical jet is merged with a subpolar jet to form a hybrid jet.
According to linear theory of baroclinic instability, the stronger the vertical shear of zonal flow (equivalent to meridional temperature gradient through thermal wind balance) is, the larger the growth rate for most unstable eddies [
This counterintuitive phenomenon, referred to as the “midwinter minimum (MWM)” or “midwinter suppression” of the NP storm-track activity, has long been investigated from various viewpoints. Nakamura (1992) and Nakamura et al. (2002) argued that an excessive propagation speed of baroclinic eddies inhibits their development of eddies by reducing their residence time within the baroclinic zone over the western Pacific. It was suggested that the excessively strong lateral shear of the extremely strong westerly jet in midwinter may also lead to suppressed growth rate of eddies (Harnik and Chang 2004; Deng and Mak 2005). This is based on the “barotropic governor” mechanism proposed by James (1987), who showed that enhanced lateral shear of a westerly jet acts to suppress baroclinic eddy growth. The importance of the excessive strength of the westerlies has been substantiated by Afargan and Kaspi (2017), who observed a clear suppression of high-frequency eddy activity also over the North Atlantic during years of a particularly strong westerly jet. A similar result was obtained for a southward-shifted jet regime over the North Atlantic by Madonna et al. (2019).
Nakamura and Sampe (2002) pointed out that upper-tropospheric eddies propagating into the midwinter NP tend to be trapped into the southward and upward shifted westerly jet core and thus separated from the surface baroclinic zone, which is unfavorable for effective baroclinic growth of eddies (Nakamura et al. 2004). Vertical structure of the westerly Pacific jet varies between the stronger and more subtropical “merged” regime (Lachmy and Harnik 2014, 2016) in midwinter and the weaker eddy-driven regime in autumn and spring, which may be responsible for the MWM (Yuval et al. 2018; Yuval and Kaspi 2018). Following on this, by tracking upper and lower level eddies in the storm track region, Hadas and Kaspi (2021) showed that the baroclinic coupling of storms diminishes when the jet is enhanced, leading to reduced baroclinic eddy growth.
Furthermore, the effect of diabatic heating associated with low-level clouds in the cold sector of a cyclone was argued by Chang (2001) and Chang and Song (2006). Through a Lagrangian tracking method, Penny et al. (2010) focused on the “seeding effect” of upper-level cyclonic eddies propagating from the Asian Continent upstream of the Pacific storm track. They argued that fewer cyclonic eddies propagating from the upstream is responsible for the MWM, while no such clear relationship was found in an independent analysis by Chang and Guo (2012). Penny et al. (2013) argued that the discrepancies are partly due to different perspectives concerning whether cold seasons are viewed individually or as a single large dataset. In addition, Park et al. (2010) and Lee et al. (2013) suggested the potential importance of the upstream orography on the MWM.
Despite these efforts, the governing mechanisms for the MWM of the NP storm-track activity have yet to be uncovered and are still under debate. The phenomenon becomes one of the most difficult issues remaining in extratropical atmospheric and climate dynamics. The MWM was successfully simulated in AGCM experiments (Christoph et al. 1997; Zhang and Held 1999), which has not led to full understanding of the detailed mechanisms though. Zhang and Held (1999) successfully simulated the MWM in a stochastic linear storm-track model, implying that it is mainly caused by linear dry dynamics. However, a similar effort by Whitaker and Sardeshmukh (1998) was unsuccessful. Using an atmospheric general circulation model (AGCM) with an idealized zonally symmetric setting, which reproduced the climatological conditions in the Pacific over the different seasons, Yuval et al. (2018) were able to reproduce the MWM demonstrating that zonal asymmetries are not essential for reproducing the MWM. Using a similar model with a seasonal dependent forcing, Novak et al. (2020) suggested the importance of the equatorward shift of the westerly jet in midwinter for the suppressed eddy activity.
Chang (2001) and Zhao and Liang (2019) examined the energetics of migratory eddies along the Pacific storm track as an attempt to understand dynamical processes underlying the MWM. Their arguments are, however, based largely on energy conversion/generation rates, which is, by definition, dependent on eddy amplitude and thus insufficient for fully explaining the mechanisms for the suppressed eddy activity. As a measure independent of eddy amplitude, Chang (2001) evaluated a local “growth rate” of a given process, defined as the related energy conversion divided by local eddy kinetic energy (EKE). Schemm and Rivière (2019) evaluated efficiency of the baroclinic conversion term, defined as the conversion normalized by the product of eddy amplitude and magnitude of background baroclinicity (Rivière and Joly 2006; Rivière et al. 2018). They suggested that the reduced efficiency of baroclinic conversion owing to an anomalous poleward tilt of eddies over the western NP poleward of ∼40°N contributes to MWM. Schemm et al. (2021) argued that more rapid development and decay of cyclones over the Kuroshio in midwinter is consistent with the reduction. Their evaluation is, however, conducted only for baroclinic conversion, and thus insufficient for full understanding of the mechanisms for MWM. Additionally, the NP storm track is zonally elongated rather than confined to the vicinity of the Kuroshio, and its MWM is therefore a basin-scale, nonlocal phenomenon. Therefore, a more comprehensive, three-dimensionally integrated energetics of eddies propagating along the NP storm track, which is independent of eddy amplitude, is required to clarify the mechanisms for the MWM, which is the main purpose of the present study.
This paper is organized as follows. Section 2 describes the data and analysis methods used in this study. Section 3 explores the climatological seasonal evolution of various Eulerian eddy statistics for the NP storm track. In section 4, the full energetics for the NP storm track are investigated, with particular focus on their seasonal evolution. Section 5 addresses long-term modulations of the MWM. Section 6 is for summary and discussion.
2. Data and analysis methods
a. Observational data
In this study, 6-hourly atmospheric variables, including geopotential height, temperature, wind velocities, and diabatic heating rate in pressure coordinates as well as sea level pressure (SLP), are obtained from the Japanese 55-year reanalysis (JRA-55) by the Japan Meteorological Agency (JMA) (Kobayashi et al. 2015; Harada et al. 2016) for the period 1958–2017. The JRA-55 has been constructed with a four-dimensional variational data assimilation (4D-Var) system with TL319 horizontal resolution (equivalent to 55-km resolution) and 60 vertical levels up to the 0.1-hPa level. Those variables on a given pressure level are available on a 1.25° × 1.25° grid. Diabatic heating rate has been decomposed into five terms: parameterized convective precipitation, large-scale gridscale precipitation, shortwave radiation, longwave radiation, and vertical diffusion. When diabatic heating associated only with precipitation is considered, we take the sum of the first two terms.
At each grid point, sub-weekly fluctuations of a given variable associated with synoptic-scale transient eddies have been extracted from the 6-hourly atmospheric reanalysis as its local deviations from their low-pass-filtered fields with an 8-day cutoff Lanczos filter. Local activity of those transient eddies is evaluated, for example, as the temporal variance based on sub-weekly fluctuations of meridional velocity or the covariance representing poleward eddy heat flux. A region of particularly large variance or covariance corresponds to a “storm track,” along which transient eddies recurrently develop. Climatological-mean fields are calculated with a 31-day running mean.
b. Energetics
In (1), primes denote high-pass-filtered fields, subscripts “C” the climatological-mean fields, and
Hereafter, we signify the normalized rate simply as “λ,” and a subscript as listed in Table 1 is added in referring to a specific conversion/generation process. The unit of this normalized energy conversion (or generation) rate is day−1, and its reciprocal thus gives a time scale of how long it would take to replenish the total eddy energy within the storm-track domain solely by a given conversion/generation process. Because the energy conversion/generation rates, by definition, depend on eddy amplitude, it is essential to evaluate the conversion/generation as its normalized rate, which is independent of eddy amplitude, for the eddy energetics for storm tracks. In the following, a negative value of λ means that eddies are transferring energy to the mean flow or low-frequency variability.
List of the symbols for the normalized rates used in this study. All quantities are expressed in the unit of day−1.
3. Climatology of the North Pacific storm-track activity
We begin our analysis with exploring the seasonal evolution of the climatological activity of the NP storm track based on the atmospheric reanalysis that covers a longer period than in previous studies. Figure 1 shows the climatological-mean seasonal evolution of meridionally averaged eddy statistics around the axes of their local maxima at individual longitudes.
Climatological seasonality in storm-track activity measured as RMS of sub-weekly fluctuations in meridional wind velocity (m s−1) at the (a) 300-hPa level (
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
The variance of sub-weekly fluctuations in 300-hPa meridional wind (
Poleward eddy heat flux (V′T′) is a measure of baroclinic development of eddies and related to the baroclinic energy conversion (CP in section 2). We focus on the lower troposphere, where V′T′ is prominent throughout the cold season. At the 850-hPa,
To investigate the typical structure of transient eddies, the local correlation in high-pass-filtered 850-hPa fluctuations between meridional wind and temperature is calculated as
(a) As in Fig. 1c, but for the correlation (%) between V′ and T′ at 850 hPa within 150°E–180° around
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
The distinct MWM in
The poleward eddy heat flux is closely related to the upward wave-activity flux and thus the vertical structure of the MWM signature. Figure 3 shows zonal sections of the wave-activity flux defined by Takaya and Nakamura (2001) on 4 December, 24 January, and 19 March, which climatologically correspond to the first peak, MWM, and second peak, respectively, of the NP storm-track activity. The flux plotted does not include the contribution of the term proportional to eddy phase speed. The eastward component of the flux in the upper troposphere therefore represents the downstream development of transient eddies pointed out by Chang (1993), and the flux is indeed dominantly eastward (Fig. 3). In addition, the flux is strongly upward from the lower- to midtroposphere over the western and central NP. In this region, eddy baroclinic growth is prominent as marked by high positive V′–T′ correlation (Fig. 2a), and abundant heat release from the ocean occurs in the cold season around the oceanic frontal zone to maintain near-surface baroclinicity (bottom panels in Fig. 3; Hotta and Nakamura 2011). In the lower half of the troposphere, the upward wave-activity flux is strongest in midwinter (Fig. 3b), owing partly to the reduction of the stability parameter under the enhanced heat supply from the ocean (not shown). However, the upward flux at ∼400-hPa is weaker in midwinter than in the shoulder seasons, when the flux exhibits a deeper structure and the strong flux reaches above the 400-hPa level (in color). In early winter, the strongest upward flux at ∼400-hPa is consistent with a peak of
(top) Zonal sections of climatological-mean wave-activity flux associated with sub-weekly transient eddies (vectors) defined by Takaya and Nakamura (2001) without the term proportional to eddy phase speed for (a) 4 Dec, (b) 24 Jan, and (c) 19 Mar. Colors indicate its vertical component (Pa m s−2). Black contours denote climatological-mean EKE (m2 s−2). Purple dotted lines denote levels at which the climatological-mean westerly wind speed is 10 m s−1 as an approximate measure of transient eddy phase speed. Quantities shown are meridionally averaged over 30°–55°N. (bottom) Horizontal maps of climatological-mean sensible heat flux from the ocean (W m−2).
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
The covariance between high-pass-filtered diabatic heating and temperature (Q′T′) contributes to diabatic energy generation of eddies (CQ in section 2). We focus on the midtroposphere, where Q′T′ is large mainly due to condensational heating. As evident in Fig. 4a, meridionally averaged
(a) As in Fig. 1c, but for the covariance (K2 day−1) between high-pass-filtered diabatic heating associated with precipitation (Q′) and temperature (T′) at 600 hPa. Contours denote climatological-mean U300 (m s−1). (b) As in Fig. 2b, but for
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
The covariance between high-pass-filtered (sign-reversed) pressure velocity and temperature (−ω′T′) contributes to the transfer from EAPE to EKE (ET in section 2). The covariance −ω′T′ is prominent in the midtroposphere. Similar to
4. Energetics for the North Pacific storm track
With the Eulerian statistics for sub-weekly transient eddies discussed above, we examine the climatological energetics for the NP storm track in the framework described in section 2b. As a typical midwinter case, distributions of vertically integrated climatological-mean energy and conversion/generation terms are shown in Fig. 5 for the 24 January, when the storm-track activity shows its minimum climatologically. Vertically integrated climatological EKE and EAPE both exhibit their maxima along the NP storm track (Figs. 5a,b). The EAPE maximum is located upstream of the EKE maximum with smaller magnitude.
Climatological energetics of sub-weekly eddies for 24 Jan. Vertically integrated climatological-mean (a) EKE and (b) EAPE (105 J m−2). (c) Vertically integrated climatological-mean CK term (color; W m−2), extended E–P flux at 300 hPa (vectors; m2 s−2), and U300 (purple contours; m s−1). (d) Vertically integrated climatological-mean CP (color),
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
The barotropic energy conversion (CK) is negative around the jet exit region (Fig. 5c), where eddies are giving up EKE to the background westerlies. This is consistent with the westerly acceleration through vorticity flux by transient eddies represented as the diverging extended E–P flux (Trenberth 1986). In agreement with previous studies (e.g., Schemm and Rivière 2019), the baroclinic energy conversion (CP) is prominent in the storm-track core just east of Japan (Fig. 5d). This is due to the collocation of strongest poleward eddy heat flux and distinct meridional gradient of climatological-mean temperature (or large EGR) in the presence of the NP oceanic frontal zone (Hotta and Nakamura 2011) and the prominent Pacific jet. The CP maximum is located slightly upstream of the EAPE maximum. The amplitude of CP is much greater than that of CK. The APE generation through diabatic heating (CQ) is overall positive along the NP storm track (Fig. 5e), reflecting the positive Q′–T′ correlation in the free troposphere. The region of the positive CQ roughly collocates with the region of positive diabatic heating in the midtroposphere. The energy transfer (ET) from EAPE to EKE is largest in the western NP between the CP and EAPE maxima. The maximum value of ET is comparable to that of CP (Fig. 5f), indicating that most of the converted APE (CP) is transferred into EKE. The eastward energy flux (EF) is strongest along the storm track and deflected northward over the eastern NP (Fig. 5g), consistent with the distributions of EKE and EAPE (Figs. 5a,b).
Those energy conversion/generation terms are integrated horizontally over the NP domain (20°–65°N, 130°W–130°E) and vertically from the surface to 100 hPa, to evaluate the comprehensive energy budget for the NP storm track. Figure 6a shows the seasonal evolution of the three-dimensionally integrated EKE and EAPE. Both the integrated EKE and EAPE are reduced slightly in midwinter between their peaks in early winter and spring, in a manner consistent with the seasonality of
(a) Climatological-mean seasonal evolution of EKE (blue), EAPE (red), and EKE + EAPE (black) integrated three-dimensionally over the NP (1018 J). A tick mark on the abscissa represents the first day of a given calendar month. (b) As in (a), but for three-dimensionally integrated power of CK (solid blue), CKLF (dashed blue), CP (solid red), CPLF (dashed red), CQ (solid green), CQ associated only with precipitation (dashed green), EF (light blue), and ET (purple). Units: 1013 W. (c) As in (b), but for λCK (solid blue), λCKLF (dashed blue), λCP (solid red), λCPLF (dashed red), λCQ (solid green), λCQ associated only with precipitation (dashed green), λEF (light blue), and λET (purple). Black line denotes λTot relevant to the budget of EKE + EAPE (viz., CK + CP + CQ + EF + CKLF + CPLF). Units: day−1.
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
Among the energy conversion/generation terms, CP and ET are much larger than the others (Fig. 6b), reflecting the baroclinic nature of sub-weekly eddies. Those terms peak in early winter before suppressed slightly in midwinter. The larger CP in winter than in autumn is consistent with both Chang (2001) and Zhao and Liang (2019). The most enhanced downstream extension of
The climatological-mean barotropic conversion CK is systematically negative (Fig. 6b), indicating that EKE is converted into the mean westerlies, which is consistent with the acceleration of the midlatitude westerlies by transient eddies. More destructive (negative) CK in midwinter indicates stronger upper-tropospheric jet acceleration by transient eddies. The stronger negative feedback by eddies acting on a stronger westerly jet is consistent with results obtained through a dry quasigeostrophic model experiment by Robert et al. (2017).
The barotropic and baroclinic energy conversions from low-frequency variability to sub-weekly transients (CKLF and CPLF) are systematically negative and positive, respectively (Fig. 6b), yielding the slightly negative net energy conversion. Though relatively small, the magnitude of negative CKLF weakens around midwinter, acting to attenuate the midwinter suppression of the storm-track activity. The offsetting contributions between CKLF and CPLF are consistent with Lau and Nath (1991), who found that transient eddies along storm tracks act to reinforce monthly-mean circulation anomalies mainly through eddy vorticity fluxes while acting to weaken monthly-mean thermal anomalies. Differences in the relative amplitude of CKLF compared to CPLF between the present study and theirs can be attributable to different definitions of low-frequency variabilities.
The flux term EF is also negative throughout the cold season, indicating that eddy energy as a net is transported out of the NP storm track. It should be noted that there is still a relatively large positive remainder (∼1014 W in winter) in the energy budget of transient eddies, which is consistent with the previous studies (e.g., Chang et al. 2002). The remainder may be primarily due to dissipation associated with surface friction and gravity waves induced by orography, which was found to be quite effective (with time scale of less than 1 day) in the planetary boundary layer (Klinker and Sardeshmukh 1992).2
By definition, the energy conversion/generation terms are proportional to the squared amplitude of transient eddies. To avoid the dependence of those terms on eddy amplitude, it is vitally important to evaluate their normalized rates (λ) based on (4). Figure 6c indicates that λCP (red line) is much higher than λCK, λCQ, λEF, λCKLF, and λCPLF. We thus confirm that baroclinic eddy growth, especially through the poleward eddy heat flux (Fig. 7a), is of primary importance for maintaining the NP storm-track activity. The λCP is highest in early winter and then reduced somewhat in midwinter and further into early spring, which is consistent with the early-winter peak of the V′–T′ correlation in the mid and upper troposphere (Fig. 2b). Interestingly, λET (purple line) maximizes in January, indicating that the internal conversion from EAPE to EKE is most effective in midwinter. It is mainly contributed to by the lower and midtropospheric ET, as seen in the ω′–T′ correlation (Fig. 4d). By contrast, λET (light blue line) and λCK (blue line) are more destructive in early and midwinter than in spring, while λCQ (green line) is nearly constant throughout the cold season. The midtropospheric Q′–T′ correlation exhibits its midwinter maximum (green line in Fig. 4b), but its effect is offset by the MWM of the corresponding correlation in the upper troposphere (blue line) and autumn peak of the lower-tropospheric counterpart (black line) as well as the effect by vertical diffusion and radiative processes, as shown below.
(a) As in Fig. 6c, but for “zonal” (dashed lines), “meridional” (dotted lines), and total (solid lines) components of λCK (blue) and λCP (red). See the text for details. (b) As in Fig. 6b, but for the total λCQ (green), λCQ associated solely with precipitation (red), vertical diffusion (purple), and radiative processes (blue). (c) As in Fig. 6c, but for the total energy fluxes through the western (red), eastern (blue), southern (purple), and northern (green) boundaries. Black line denotes the net total energy flux. Positive (negative) values mean inward (outward) fluxes. (d),(e) As in (c), but for energy flux term associated solely with energy advected by climatological-mean flow in (d) and eddy geopotential flux in (e).
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
When CP is combined with the other processes (CK, CQ, EF, CKLF, and CPLF), the resultant λTot (black line in Fig. 6c) for the storm-track energetics shows a fairly apparent midwinter suppression. This result suggests that processes other than CP, especially CK and EF, are likely to contribute to the formation of the MWM, particularly to the early-spring recovery. The λCQ associated only with precipitation is kept nearly constant from early to midwinter (green dashed line) before slightly reduced into spring.
In the following, some of the conversion/generation terms are investigated in more detail. Decomposition of λCP into the two components related to the zonal heat flux
(a)–(c) Vertically integrated climatological-mean “zonal” component of the CK term (color; W m−2), extended E–P flux at 300 hPa (vectors; m2 s−2), and U300 (black contours; m s−1) for (a) 4 Dec, (b) 24 Jan, and (c) 19 Mar, respectively. (d)–(f) As in (a)–(c), respectively, but for “meridional” component of the CK term.
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
The total λCQ is modestly positive throughout the cold season (green line in Fig. 7b). Consistent with the climatological-mean midtropospheric Q′T′ (Fig. 4a), λCQ associated solely with precipitation and the corresponding latent heat release peaks in early winter (red line in Fig. 7b). It is systematically greater than the total λCQ, because the latter includes dominant thermal damping through heat exchanges with the underlying ocean (purple line in Fig. 7b). This damping effect is strongest in early winter over the relatively warm ocean, acting to weaken the early-winter peak of λCQ associated only with precipitation. Finally, λCQ associated solely with radiative processes is almost negligible if integrated three-dimensionally (blue line in Fig. 7b).
The seasonality of CQ is overall consistent with Chang (2001), who utilized a rather coarse-resolution AGCM experiment, and with Zhao and Liang (2019), who evaluated it as the residual of the EAPE equation. Compared with our result, a slightly stronger early-winter peak of CQ found by Zhao and Liang (2019) may be due to the different latitude of the northern boundary of the analysis domain and to their vertical integration above the 1000-hPa level. The larger total CQ in October than in midwinter especially over the western NP shown by Chang (2001) has been confirmed by the integration only over the western NP domain (not shown). The present study utilizes the individual diabatic heating components obtained in a more recent reanalysis data (JRA-55). Although the results depend on a particular forecast model and data assimilation system, the evaluation of the individual heating components has become possible by taking advantage of the reanalysis data.
The net total λEF can be decomposed into the contributions from energy fluxes passing through the four (western, eastern, northern, and southern) lateral boundaries and further into the advective and eddy flux components (Figs. 7c–e). In the net total λEF (Fig. 7c), the energy fluxes through the western and eastern boundaries (red and blue lines, respectively) are largely canceled out. Nevertheless, the flux out of the eastern boundary is systematically larger, especially in early winter, when the connection into the Atlantic storm track is most obvious (not shown). Both the energy inflow through the western boundary and outflow through the eastern boundary are mostly due to the advective effect by the mean westerlies (Fig. 7d). Still, the net zonal energy influx contributes to the MWM of the net λEF, especially for the spring recovery (light blue line in Fig. 6c). The corresponding energy fluxes through the southern and northern boundaries (purple and green lines, respectively, in Fig. 7c) are both relatively weak and outward, and they tend to peak in midwinter. Unlike the zonal energy fluxes, the eddy geopotential flux contributes substantially to the meridional energy outflow primarily through the northern boundary (Figs. 7d,e). Overall, the net contributions from eddies and the mean flow are comparable to the meridional energy outflow.
Recently, Zhao and Liang (2019) investigated individual energy conversion/generation terms in the energetics associated with the NP storm track, but their evaluations included direct contributions from eddy amplitude itself. Chang (2001) evaluated some energy conversion rates and their “growth rate” independent of eddy amplitude. Their evaluation was, however, conducted for local conversion/generation terms and only with EKE, which may not necessarily lead to full delineation of the maintenance mechanisms for the NP storm-track activity. This is because total eddy energy is better represented as the sum of EKE and EAPE, the latter of which is replenished mainly through the baroclinic conversion. In addition, the MWM is a basin-scale phenomenon, and therefore the three-dimensionally integrated energy budget analysis is required. Recently, Schemm and Rivière (2019) evaluated the efficiency of CP, in which baroclinic conversion rate is normalized locally by the product of total eddy energy and the magnitude of background baroclinicity. It is certainly advantageous for measuring the ability of eddies themselves in converting energy, though without representing an actual conversion rate of a given process for eddy energy. We point out that their evaluation was performed only for CP over a specific domain in the western NP. The present study performs a comprehensive evaluation of λ for each of the energy conversion/generation terms based on the three-dimensionally integrated statistics for the entire NP storm track, which can provide straightforward information about the ability of a given process to act against damping and dissipating processes with a dimensional time scale. We have revealed that λTot, in which all the processes are combined, is indeed reduced in midwinter in a manner consistent with the MWM of the storm-track activity, although its spring recovery is marked only as subtle peak of λTot compared to the eddy energy. These results are consistent with Hadas and Kaspi (2021), who showed, via tracking of upper- and lower-level eddy activity, a faster barotropization of the westward tilt of baroclinic eddies when the jet is stronger. This provides a mechanism leading to the energy conversions demonstrated in this section.
5. Long-term modulations of the MWM
Nakamura et al. (2002) found that the MWM of the NP storm-track activity became less distinct in a recent period after 1986 than in the earlier period. Figure 9 shows interannual variability of 5-yr running-mean storm-track activity including its seasonality. In agreement with Nakamura et al. (2002), the MWM of
(a) RMS of 5-yr mean of 31-day running mean
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
In recognition of the distinct transition from 1985/86 to 1986/87, climatological-mean fields for 1958/59–1985/86 and 1986/87–2016/17 are compared in the following to highlight the long-term modulations of the MWM of the NP storm-track activity. Figure 10 compares the climatological-mean seasonal evolutions of RMS
(a),(b) As in Fig. 1a, but for the periods of 1958/59–1985/86 and 1986/87–2016/17, respectively. (c),(d) As in Fig. 1d, but for the periods of 1958/59–1985/86 and 1986/87–2016/17, respectively.
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
Similar to
(a)–(d) As in Fig. 10c, but for (a)
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
Applying the same procedures as in section 4, we evaluate the energetics of sub-weekly transient eddies along the NP storm track separately for the earlier and later periods. In the earlier period, there are more apparent MWM signals in both spatially integrated EKE and EAPE (Fig. 12a) compared to those for the entire period (Fig. 6a), while the corresponding signals become ambiguous in the later period owing mainly to the enhanced midwinter eddy activity (Fig. 12b). This is consistent with the abovementioned results based on the Eulerian statistics.
(a),(b) As in Fig. 6a, but for the periods of 1958/59–1985/86 and 1986/87–2016/17, respectively. Units are 1018 J. (c),(d) As in Fig. 6c, but for the periods of 1958/59–1985/86 and 1986/87–2016/17, respectively. Units: day−1.
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
As illustrated in Figs. 12c and 12d, the normalized rates of energy conversion/generation terms evaluated with total eddy energy (EKE + EAPE) separately for the two periods exhibit their long-term modulations. In the earlier period, λCP (red line in Fig. 12c) undergoes a clearer midwinter reduction, leading to a noticeable midwinter suppression (with a slight minimum) of λTot (black line). In the later period, by contrast, the midwinter reduction of λCP after its early-winter peak is weaker, yielding no MWM of λTot (Fig. 12d). The weaker midwinter reduction of λCP for the later period is consistent with the more apparent midwinter maximum of the V′–T′ correlation than in the earlier period (Figs. 11d,h). Interestingly, the negative λCK (blue lines in Figs. 12c,d) from early to midwinter as an indication of eddy acceleration of the mean westerlies is enhanced in the later period, to which both the “zonal” and “meridional” components contribute (not shown). This indicates that the interdecadal enhancement of the midwinter Pacific storm-track activity is unlikely due to this barotropic effect. In fact, this interdecadal change in λCK cannot offset its counterpart in λCP. The midwinter λEF (light blue lines in Figs. 12c,d) is negative and roughly comparable between the two periods. From late winter into spring, the seasonal weakening of the negative λEF is more gradual in the later period, contributing to the weaker spring recovery of eddy energy. A comparable contribution is made by λCK in early spring (Fig. B1). The net energy conversion from low-frequency variabilities (dashed lines in Figs. 12c,d) is negative, and its contribution to the midwinter suppression of λTot is slightly enhanced in the earlier period. When all the processes are combined, λTot for the midwinter NP storm track is thus slightly higher in the later period than in the earlier period (black lines in Figs. 12c,d). We have confirmed that the long-term modulations discussed above are statistically significant (see appendix B). The results show that the long-term modulations of the MWM can be interpreted, at least in part, from the viewpoint of eddy energetics.
6. Summary and discussion
In this study, we have delineated the detailed seasonal evolution of climatological-mean Eulerian statistics and the energetics for sub-weekly transient eddies along the NP storm track. As its distinct characteristic, the MWM of eddy activity is prominent particularly in upper-tropospheric V′V′ and lower-tropospheric V′T′ and T′T′. The distinct midwinter reduction in the midtropospheric upward wave-activity flux is consistent with the vertical structure of the MWM in the upper-level eddy activity V′V′ and EKE.
To investigate the maintenance mechanisms for the MWM, we have comprehensively evaluated the energetics for sub-weekly transient eddies along the NP storm track. In recognition of the dependence of any conversion/generation terms to eddy amplitude, we have calculated their energy conversion/generation rate (λ) normalized by the eddy total energy (EKE + EAPE), which is independent of eddy amplitude. We have revealed that the net normalized energy conversion/generation rate (λTot) is indeed suppressed in midwinter. The reduction in the normalized EAPE conversion is found to play a substantial role especially in the reduction from its early-winter peak, whereas that of the energy outflux from the NP is particularly important for the spring recovery, among other processes contributing to the formation of the MWM in λTot. Our results suggest that interaction between transient eddies and climatological-mean background state, especially the westerly jet and associated baroclinic zone, is central to the MWM over the NP and multiple processes are responsible. The midwinter reduction of λCP indicates that baroclinic eddy growth is suppressed in the real atmosphere despite the midwinter maximum of EGR (Hadas and Kaspi 2021). It is consistent with Schemm and Rivière (2019), whose results are based on the further normalization by the magnitude of background baroclinicity and only for the baroclinic energy conversion. We have verified that our results are quite similar when the ERA-Interim reanalysis (Dee et al. 2011) is used for the period of 1979/80–2016/17 in place of the JRA-55 (not shown).
In addition, we have applied the framework of energetics to the long-term modulations of the MWM of the NP storm-track activity. As noted by Nakamura et al. (2002), the transition from 1985/86 to 1986/87 was striking, after which the MWM signal has become much less distinct in both
Considering the modulation of the total eddy energy, the present study has shown that the long-term modulation in the MWM of the NP storm-track activity is consistent with the corresponding modulation in λTot in which all the energy conversion/generation terms are combined. In the earlier period (1958/59–1985/86), when the MWM was more prominent, λTot for the NP storm track exhibited a distinct MWM. The higher midwinter λCP after the transition is consistent with the stronger positive V′–T′ correlation over the western NP observed in winters of weak midwinter suppression found by Nakamura et al. (2002). Deng and Mak (2006) suggested that the efficiency of the barotropic energy conversion (CK) was comparably important to the energy transfer (from EAPE to EKE) associated with vertical motion (ET) in the difference between winters of the stronger and weaker MWM. Unlike in the present study, however, they did not focus on the EAPE conversion from the background state (CP). In addition, the present study also reveals that the most effective (negative) CK in midwinter is mainly due to the eddy momentum flux acting on the enhanced jet diffluence in its exit but not due to the enhanced lateral shear of the jet as assumed in the “barotropic governor” mechanism.
To assess whether the more distinct MWM in λTot is primarily due to eddy Eulerian statistics or climatological-mean background state, Fig. 13 compares the energetics for migratory eddies along the NP storm track between two different combinations of the eddy statistics and the background state: one is based on the eddy statistics in the earlier period and the background state in the later period, and vice versa for the other. By definition, EKE is independent of the background state. The eddy contributions to EAPE in the earlier and later periods are similar to the actual EAPE in the respective periods (Fig. 12). The λTot shows its more distinct midwinter suppression if evaluated for the eddy component in the earlier period under the background state in the later period (Figs. 13b,d). This result is quite similar to our result of long-term modulations shown in Fig. 12. These results suggest that the long-term modulations in the midwinter suppression of λTot is contributed to predominantly by the modulated eddy Eulerian statistics, especially in the poleward eddy heat flux, rather than directly by interdecadal modulation in the background westerlies.
(a) Climatological-mean seasonal evolution of EKE (blue), EAPE (red), and EKE + EAPE (black) integrated three-dimensionally over the NP (1018 J), based on eddy statistics in 1958/59–1985/86 and the background state for 1986/87–2016/17. (b) As in (a), but for λCK (blue), λCP (red), λCQ (solid green), λCQ associated solely with precipitation (green dotted), λEF (light blue), and λEF (purple). Black line denotes λTot relevant to the budget of EKE + EAPE associated with the background state (viz., CK + CP + CQ + EF). Units: day−1. (c),(d) As in (a) and (b), respectively, but based on the eddy statistics in 1986/87–2016/17 and the background state in 1958/59–1985/86.
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
The present study is the first to demonstrate that the net normalized energy conversion/generation rate λTot for sub-weekly transients along the NP storm track is indeed reduced in midwinter, largely owing to the baroclinic energy conversion (CP) via poleward eddy heat flux with weakened temperature fluctuations especially for the reduction from early winter, among other mechanisms that also contribute positively to the MWM, including the energy flux term (EF) for the spring recovery. We have verified that the framework of energetics is useful for investigating the mechanisms for the MWM of the NP storm-track activity. It gives us a perspective of the MWM that encompasses various mechanisms proposed by previous studies. The barotropic conversion CK is presumably related to the excessively strong westerly jet, especially its diffluence in its exit region. The contribution of diabatic heating CQ includes the effect of low-level clouds in the cold sector of each cyclone. The net energy input through the lateral boundaries EF includes the “seeding effect” of cyclones from upstream. Specifically, the increasing λEF from the energy influx through the western boundary from midwinter into spring may be related to the upstream “seeding effect” (Penny et al. 2010; Zhao and Liang 2019), though its increase from autumn into midwinter is not consistent. In addition, the minimum or suppression of λCP is compatible with the suppressed baroclinic growth through the trapping of upper-tropospheric eddies into the subtropical jet core (Nakamura and Sampe 2002) over the central NP. The λCP is also compatible with the decreased Lagrangian coupling between the upper and lower part of the baroclinic wave over the NP (Hadas and Kaspi 2021).
Though its contribution to λTot is rather small, the positive zonal component of CP in midwinter (Fig. 7a) does not seem consistent with Schemm and Rivière (2019). It may be attributable to the difference between their analysis method and ours; they focused mainly on the energy conversion efficiency over the western NP north of ∼40°N, while the λCP is calculated for the entire troposphere and over the whole NP in this study. The sign of zonal eddy heat flux is opposite between the mid- to upper troposphere and lower troposphere (not shown), which ensures the advantage of the three-dimensional integration in the present study. Nonetheless, the discrepancy should be further investigated in future studies, as well as more detailed characteristics of the mechanisms for MWM. Particularly, the reduction in the upper-tropospheric V′–T′ correlation in midwinter and early spring, which probably contributes to the λCP reduction, will be addressed. Our results suggest that those processes previously discussed can be all operative in the formation of the MWM over the NP storm-track activity.
Recently, Schemm and Rivière (2019) and Schemm et al. (2021) have shown the usefulness of the Lagrangian perspective for studying MWM. The number of cyclones coming from the East China Sea (Schemm et al. 2021) may contribute also to the long-term modulation of MWM. Further investigation of the contribution from three-dimensional structure and characteristics of migratory cyclones to the MWM of the NP storm-track activity, focusing on Lagrangian perspective as in Schemm and Schneider (2018), will certainly be informative.
Acknowledgments.
The authors are grateful to the three anonymous reviewers for their sound criticism and constructive comments on the earlier versions of this paper. This study is supported in part by the Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT) through the Arctic Challenge for Sustainability (ArCS-II), by the Japan Science and Technology Agency through COI-NEXT JPMJPF2013, by the Japanese Ministry of Environment through Environment Research and Technology Development Fund JPMEERF20192004, and by the Japan Society for the Promotion of Science (JSPS) through Grants-in-Aid for Scientific Research JP18H01278, JP19H05702 (on Innovative Areas 6102), and 20H01970. It is also supported by JSPS-ISF Joint Research Project (JPJSBP120218403). Y.K. acknowledges support from the JSPS Invitational Fellowship for Research in Japan that supported a sabbatical at the University of Tokyo and ignited this collaboration, for support from the Research Center for Advanced Technology and Science at the University of Tokyo and the Israeli Science Foundation (Grant 996/20).
Data availability statement.
The JRA-55 atmospheric reanalysis is available online from the Japan Meteorological Agency at https://jra.kishou.go.jp/JRA-55/index_en.html as cited in Kobayashi et al. (2015) and Harada et al. (2016).
APPENDIX A
Derivation of CKLF
APPENDIX B
Statistical Significance of Long-Term Modulation
Statistical significance of the modulation of normalized conversion/generation rates is assessed here. To evaluate statistical significance of differences in the rates between the earlier and later periods by t test, standard deviations of the rates are calculated based on a random permutation of cold seasons by 20 times, the number of which is the same as that of the respective periods (i.e., 28 and 31), within the entire period (i.e., 59). Both climatological-mean eddy statistics and background states are calculated based on a given set of randomly selected cold seasons. Figure B1 shows that most of the distinct features in the differences of the normalized rates are statistically significant, including the midwinter reduction and early-winter/spring enhancement of λTot.
Differences in climatological-mean λCK (blue), λCP (red), λCQ (green), λEF (light blue), λEF (purple), and λTot (black) over the North Pacific between the earlier (1958–85) and the later (1986–2016) periods. Units: day−1. Closed circles indicate statistically significant differences at the 99% confidence level by Student’s t test.
Citation: Journal of Climate 35, 4; 10.1175/JCLI-D-21-0123.1
REFERENCES
Afargan, H., and Y. Kaspi, 2017: A midwinter minimum in North Atlantic storm track intensity in years of a strong jet. Geophys. Res. Lett., 44, 12–511, https://doi.org/10.1002/2017GL075136.
Blackmon, M. L., 1976: A climatological spectral study of the 500 mb geopotential height of the Northern Hemisphere. J. Atmos. Sci., 33, 1607–1623, https://doi.org/10.1175/1520-0469(1976)033<1607:ACSSOT>2.0.CO;2.
Blackmon, M. L., J. M. Wallace, N. C. Lau, and S. L. Mullen, 1977: An observational study of the Northern Hemisphere wintertime circulation. J. Atmos. Sci., 34, 1040–1053, https://doi.org/10.1175/1520-0469(1977)034<1040:AOSOTN>2.0.CO;2.
Chang, E. K., 1993: Downstream development of baroclinic waves as inferred from regression analysis. J. Atmos. Sci., 50, 2038–2053, https://doi.org/10.1175/1520-0469(1993)050<2038:DDOBWA>2.0.CO;2.
Chang, E. K., 2001: GCM and observational diagnoses of the seasonal and interannual variations of the Pacific storm track during the cool season. J. Atmos. Sci., 58, 1784–1800, https://doi.org/10.1175/1520-0469(2001)058<1784:GAODOT>2.0.CO;2.
Chang, E. K., and S. Song, 2006: The seasonal cycles in the distribution of precipitation around cyclones in the western North Pacific and Atlantic. J. Atmos. Sci., 63, 815–839, https://doi.org/10.1175/JAS3661.1.
Chang, E. K., and Y. Guo, 2012: Is Pacific storm-track activity correlated with the strength of upstream wave seeding? J. Climate, 25, 5768–5776, https://doi.org/10.1175/JCLI-D-11-00555.1.
Chang, E. K., S. Lee, and K. L. Swanson, 2002: Storm track dynamics. J. Climate, 15, 2163–2183, https://doi.org/10.1175/1520-0442(2002)015<02163:STD>2.0.CO;2.
Christoph, M., U. Ulbrich, and P. Speth, 1997: Midwinter suppression of Northern Hemisphere storm track activity in the real atmosphere and in GCM experiments. J. Atmos. Sci., 54, 1589–1599, https://doi.org/10.1175/1520-0469(1997)054<1589:MSONHS>2.0.CO;2.
Dee, D. P., and Coauthors, 2011: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828.
Deng, Y., and M. Mak, 2005: An idealized model study relevant to the dynamics of the midwinter minimum of the Pacific storm track. J. Atmos. Sci., 62, 1209–1225, https://doi.org/10.1175/JAS3400.1.
Deng, Y., and M. Mak, 2006: Nature of the differences in the intraseasonal variability of the Pacific and Atlantic storm tracks: A diagnostic study. J. Atmos. Sci., 63, 2602–2615, https://doi.org/10.1175/JAS3749.1.
Eady, E. T., 1949: Long waves and cyclone waves. Tellus, 1, 33–52, https://doi.org/10.3402/tellusa.v1i3.8507.
Hadas, O., and Y. Kaspi, 2021: Suppression of baroclinic eddies by strong jets. J. Atmos. Sci., 78, 2445–2457, https://doi.org/10.1175/JAS-D-20-0289.1.
Harada, Y., and Coauthors, 2016: The JRA-55 Reanalysis: Representation of atmospheric circulation and climate variability. J. Meteor. Soc. Japan, 94, 269–302, https://doi.org/10.2151/jmsj.2016-015.
Harnik, N., and E. K. Chang, 2004: The effects of variations in jet width on the growth of baroclinic waves: Implications for midwinter Pacific storm track variability. J. Atmos. Sci., 61, 23–40, https://doi.org/10.1175/1520-0469(2004)061<0023:TEOVIJ>2.0.CO;2.
Hotta, D., and H. Nakamura, 2011: On the significance of the sensible heat supply from the ocean in the maintenance of the mean baroclinicity along storm tracks. J. Climate, 24, 3377–3401, https://doi.org/10.1175/2010JCLI3910.1.
James, I. N., 1987: Suppression of baroclinic instability in horizontally sheared flows. J. Atmos. Sci., 44, 3710–3720, https://doi.org/10.1175/1520-0469(1987)044<3710:SOBIIH>2.0.CO;2.
Jiang, T., Y. Deng, and W. Li, 2013: Local kinetic energy budget of high-frequency and intermediate-frequency eddies: Winter climatology and interannual variability. Climate Dyn., 41, 961–976, https://doi.org/10.1007/s00382-013-1684-1.
Klinker, E., and P. D. Sardeshmukh, 1992: The diagnosis of mechanical dissipation in the atmosphere from large-scale balance requirements. J. Atmos. Sci., 49, 608–627, https://doi.org/10.1175/1520-0469(1992)049<0608:TDOMDI>2.0.CO;2.
Kobayashi, S., and Coauthors, 2015: The JRA-55 Reanalysis: General specifications and basic characteristics. J. Meteor. Soc. Japan, 93, 5–48, https://doi.org/10.2151/jmsj.2015-001.
Kosaka, Y., and H. Nakamura, 2010: Mechanisms of meridional teleconnection observed between a summer monsoon system and a subtropical anticyclone. Part I: The Pacific–Japan pattern. J. Climate, 23, 5085–5108, https://doi.org/10.1175/2010JCLI3413.1.
Lachmy, O., and N. Harnik, 2014: The transition to a subtropical jet regime and its maintenance. J. Atmos. Sci., 71, 1389–1409, https://doi.org/10.1175/JAS-D-13-0125.1.
Lachmy, O., and N. Harnik, 2016: Wave and jet maintenance in different flow regimes. J. Atmos. Sci., 73, 2465–2484, https://doi.org/10.1175/JAS-D-15-0321.1.
Lau, N. C., and M. J. Nath, 1991: Variability of the baroclinic and barotropic transient eddy forcing associated with monthly changes in the midlatitude storm tracks. J. Atmos. Sci., 48, 2589–2613, https://doi.org/10.1175/1520-0469(1991)048<2589:VOTBAB>2.0.CO;2.
Lee, S., and H. K. Kim, 2003: The dynamical relationship between subtropical and eddy-driven jets. J. Atmos. Sci., 60, 1490–1503, https://doi.org/10.1175/1520-0469(2003)060<1490:TDRBSA>2.0.CO;2.
Lee, S. S., J. Y. Lee, K. J. Ha, B. Wang, A. Kitoh, Y. Kajikawa, and M. Abe, 2013: Role of the Tibetan Plateau on the annual variation of mean atmospheric circulation and storm-track activity. J. Climate, 26, 5270–5286, https://doi.org/10.1175/JCLI-D-12-00213.1.
Madonna, E., C. Li, and J. J. Wettstein, 2019: Suppressed eddy driving during southward excursions of the North Atlantic jet on synoptic to seasonal time scales. Atmos. Sci. Lett., 20, e937, https://doi.org/10.1002/asl.937.
Martineau, P., H. Nakamura, Y. Kosaka, and A. Yamamoto, 2020: Importance of a vertically tilting structure for energizing the North Atlantic Oscillation. Sci. Rep., 10, 12671, https://doi.org/10.1038/s41598-020-69551-5.
Nakamura, H., 1992: Midwinter suppression of baroclinic wave activity in the Pacific. J. Atmos. Sci., 49, 1629–1642, https://doi.org/10.1175/1520-0469(1992)049<1629:MSOBWA>2.0.CO;2.
Nakamura, H., and T. Sampe, 2002: Trapping of synoptic-scale disturbances into the North-Pacific subtropical jet core in midwinter. Geophys. Res. Lett., 29, 2002, https://doi.org/10.1029/2002GL015535.
Nakamura, H., M. Tanaka, and J. M. Wallace, 1987: Horizontal structure and energetics of Northern Hemisphere wintertime teleconnection patterns. J. Atmos. Sci., 44, 3377–3391, https://doi.org/10.1175/1520-0469(1987)044<3377:HSAEON>2.0.CO;2.
Nakamura, H., T. Izumi, and T. Sampe, 2002: Interannual and decadal modulations recently observed in the Pacific storm track activity and East Asian winter monsoon. J. Climate, 15, 1855–1874, https://doi.org/10.1175/1520-0442(2002)015<1855:IADMRO>2.0.CO;2.
Nakamura, H., T. Sampe, Y. Tanimoto, and A. Shimpo, 2004: Observed associations among storm tracks, jet streams, and midlatitude oceanic fronts. Earth Climate: The Ocean–Atmosphere Interaction, Geophys. Monogr., Vol. 147, Amer. Geophys. Union, 329–345.
Novak, L., T. Schneider, and F. Ait-Chaalal, 2020: Midwinter suppression of storm tracks in an idealized zonally symmetric setting. J. Atmos. Sci., 77, 297–313, https://doi.org/10.1175/JAS-D-18-0353.1.
Park, H. S., J. C. Chiang, and S. W. Son, 2010: The role of the central Asian mountains on the midwinter suppression of North Pacific storminess. J. Atmos. Sci., 67, 3706–3720, https://doi.org/10.1175/2010JAS3349.1.
Penny, S., G. H. Roe, and D. S. Battisti, 2010: The source of the midwinter suppression in storminess over the North Pacific. J. Climate, 23, 634–648, https://doi.org/10.1175/2009JCLI2904.1.
Penny, S., D. S. Battisti, and G. H. Roe, 2013: Examining mechanisms of variability within the Pacific storm track: Upstream seeding and jet-core strength. J. Climate, 26, 5242–5259, https://doi.org/10.1175/JCLI-D-12-00017.1.
Rivière, G., and A. Joly, 2006: Role of the low-frequency deformation field on the explosive growth of extratropical cyclones at the jet exit. Part I: Barotropic critical region. J. Atmos. Sci., 63, 1965–1981, https://doi.org/10.1175/JAS3728.1.
Rivière, G., S. Berthou, G. Lapeyre, and M. Kageyama, 2018: On the reduced North Atlantic storminess during the last glacial period: The role of topography in shaping synoptic eddies. J. Climate, 31, 1637–1652, https://doi.org/10.1175/JCLI-D-17-0247.1.
Robert, L., G. Rivière, and F. Codron, 2017: Positive and negative eddy feedbacks acting on midlatitude jet variability in a three-level quasigeostrophic model. J. Atmos. Sci., 74, 1635–1649, https://doi.org/10.1175/JAS-D-16-0217.1.
Schemm, S., and T. Schneider, 2018: Eddy lifetime, number, and diffusivity and the suppression of eddy kinetic energy in midwinter. J. Climate, 31, 5649–5665, https://doi.org/10.1175/JCLI-D-17-0644.1.
Schemm, S., and G. Rivière, 2019: On the efficiency of baroclinic eddy growth and how it reduces the North Pacific storm-track intensity in midwinter. J. Climate, 32, 8373–8398, https://doi.org/10.1175/JCLI-D-19-0115.1.
Schemm, S., H. Wernli, and H. Binder, 2021: The storm-track suppression over the western North Pacific from a cyclone life-cycle perspective. Wea. Climate Dyn., 2, 55–69, https://doi.org/10.5194/wcd-2-55-2021.
Stephenson, D. B., 1995: The impact of changing the horizontal diffusion scheme on the northern winter climatology of a general circulation model. Quart. J. Roy. Meteor. Soc., 121, 211–226, https://doi.org/10.1002/qj.49712152110.
Sung, M. K., H. Y. Jang, B. M. Kim, S. W. Yeh, Y. S. Choi, and C. Yoo, 2019: Tropical influence on the North Pacific Oscillation drives winter extremes in North America. Nat. Climate Change, 9, 413–418, https://doi.org/10.1038/s41558-019-0461-5.
Takaya, K., and H. Nakamura, 2001: A formulation of a phase-independent wave-activity flux for stationary and migratory quasigeostrophic eddies on a zonally varying basic flow. J. Atmos. Sci., 58, 608–627, https://doi.org/10.1175/1520-0469(2001)058<0608:AFOAPI>2.0.CO;2.
Tanaka, S., K. Nishii, and H. Nakamura, 2016: Vertical structure and energetics of the western Pacific teleconnection pattern. J. Climate, 29, 6597–6616, https://doi.org/10.1175/JCLI-D-15-0549.1.
Trenberth, K. E., 1986: An assessment of the impact of transient eddies on the zonal flow during a blocking episode using localized Eliassen–Palm flux diagnostics. J. Atmos. Sci., 43, 2070–2087, https://doi.org/10.1175/1520-0469(1986)043<2070:AAOTIO>2.0.CO;2.
Whitaker, J. S., and P. D. Sardeshmukh, 1998: A linear theory of extratropical synoptic eddy statistics. J. Atmos. Sci., 55, 237–258, https://doi.org/10.1175/1520-0469(1998)055<0237:ALTOES>2.0.CO;2.
Yuval, J., and Y. Kaspi, 2018: Eddy sensitivity to jet characteristics. J. Atmos. Sci., 75, 1371–1383, https://doi.org/10.1175/JAS-D-17-0139.1.
Yuval, J., H. Afargan, and Y. Kaspi, 2018: The relation between the seasonal changes in jet characteristics and the Pacific midwinter minimum in eddy activity. Geophys. Res. Lett., 45, 9995–10 002, https://doi.org/10.1029/2018GL078678.
Zhang, Y., and I. M. Held, 1999: A linear stochastic model of a GCM’s midlatitude storm tracks. J. Atmos. Sci., 56, 3416–3435, https://doi.org/10.1175/1520-0469(1999)056<3416:ALSMOA>2.0.CO;2.
Zhao, Y., and X. S. Liang, 2019: Causes and underlying dynamic processes of the mid-winter suppression in the North Pacific storm track. Sci. China Earth Sci., 62, 872–890, https://doi.org/10.1007/s11430-018-9310-5.
The definition of the term is the same as that of efficiency used in some previous studies (e.g., Nakamura et al. 1987; Stephenson 1995; Kosaka and Nakamura 2010; Tanaka et al. 2016; Sung et al. 2019; Martineau et al. 2020). We nevertheless use the term normalized rate (λ) in this study, to avoid any confusion with the definition of efficiency as the ratio of energy conversion/generation divided by the product of eddy energy and the magnitude of background baroclinicity (e.g., Schemm and Rivière 2019), which is analogous to thermal efficiency in thermodynamics. Additionally, the energy conversion/generation terms are normalized by eddy total energy (EKE + EAPE) in this study, rather than solely by EKE as in Chang (2001). Normalized by both EKE and EAPE is advantageous, because a particular conversion/generation term is directly related to either EKE or EAPE, which is effectively transferred into the other by the ET term. The results are essentially the same when normalized solely by EKE.
The contributions from surface friction and gravity waves can be approximately estimated through parameterized surface momentum fluxes provided by JRA-55. For example, a horizontally averaged monthly-mean value of the sum of those fluxes multiplied by 10-m winds is ∼1.2 × 1014 W within the domain (130°E–130°W, 20°–65°N) in January 2017.