the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
The storm-track suppression over the western North Pacific from a cyclone life-cycle perspective
Sebastian Schemm
Heini Wernli
Hanin Binder
Surface cyclones that feed the western part of the North Pacific storm track and experience a midwinter suppression originate from three regions: the East China Sea (∼30∘ N), the Kuroshio extension (∼35∘ N), and downstream of Kamchatka (∼53∘ N). In midwinter, in terms of cyclone numbers, Kuroshio (45 %) and Kamchatka (40 %) cyclones dominate in the region where eddy kinetic energy is suppressed, while the relevance of East China Sea cyclones increases from winter (15 %) to spring (20 %). The equatorward movement of the baroclinicity and the associated upper-level jet toward midwinter influence cyclones from the three genesis regions in different ways. In January, Kamchatka cyclones are less numerous and less intense, and their lifetime shortens, broadly consistent with the reduced baroclinicity in which they grow. The opposite is found for East China Sea cyclones, which in winter live longer, are more intense, and experience more frequently explosive deepening. The fraction of explosive East China Sea cyclones is particularly high in January when they benefit from the increased baroclinicity in their environment. Again, a different and more complex behavior is found for Kuroshio cyclones. In midwinter, their number increases, but their lifetime decreases; on average they reach higher intensity in terms of minimum sea level pressure, but the fraction of explosively deepening cyclones decreases and the latitude where maximum growth occurs shifts equatorward. Therefore, the life cycle of Kuroshio cyclones seems to be accelerated in midwinter with a stronger and earlier but also shorter deepening phase followed by an earlier decay. Once they reach the latitude where eddy kinetic energy is suppressed in midwinter, their baroclinic conversion efficiency is strongly reduced. Together, this detailed cyclone life-cycle analysis reveals that the North Pacific storm-track suppression in midwinter is related to fewer and weaker Kamchatka cyclones and to more equatorward intensifying and then more rapidly decaying Kuroshio cyclones. The less numerous cyclone branch from the East China Sea partially opposes the midwinter suppression. The cyclones passing through the suppressed region over the western North Pacific do not propagate far downstream and decay in the central North Pacific. The behavior of cyclones in the eastern North Pacific requires further analysis.
- Article
(5865 KB) - Full-text XML
-
Supplement
(1034 KB) - BibTeX
- EndNote
Nakamura (1992) identified the contrasting intraseasonal cycles of the North Atlantic and North Pacific storm tracks in winter. While different measures of storm-track activity over the Atlantic experience a single peak in midwinter, the Pacific storm-track activity has two peaks, one in late autumn and another one in early spring (Nakamura, 1992). Over the North Atlantic, the seasonal cycle of the storm-track activity is broadly consistent with the seasonal cycle of mean baroclinicity, but this is not the case over the North Pacific where the midwinter suppression of the storm-track activity occurs at the time of maximum surface baroclinicity and jet strength. The midwinter suppression affects a large number of eddy measures, such as the eddy heat and momentum fluxes and the eddy kinetic energy, as well as the baroclinic and barotropic conversion rates, whereas measures of the mean background flow, such as the monthly mean jet strength and the Eady growth rate, are not affected by the suppression (Schemm and Schneider, 2018). The midwinter suppression is most pronounced in the upper troposphere and almost absent in the lower troposphere. Also, the total number of surface cyclones does not experience a suppression (Schemm and Schneider, 2018). The atypical intraseasonal cycle of baroclinic waves over the North Pacific has triggered considerable research during the last 2 decades. In the following, we give an overview of the current understanding of the midwinter suppression.
Factors that have been suggested to contribute to the midwinter suppression can be categorized into contributions from barotropic, baroclinic, and upstream seeding processes. Among the barotropic contributions are the increase in the horizontal wind shear near the jet and the deformation acting on baroclinic wave packets in a stronger and more narrow jet stream during midwinter (James, 1987; Nakamura, 1993; Harnik and Chang, 2004; Deng and Mak, 2005). However, in idealized simulations it was shown that such a barotropic governing mechanism is not symmetric in time around the suppression, and the increase in the horizontal shear lags the onset of the suppression (Novak et al., 2020). With regard to upstream seeding processes, a reduction in the amplitude and frequency of upper-level eddies propagating from upstream into the North Pacific has been suggested to play a crucial role in the formation of the midwinter suppression (Penny et al., 2010, 2011, 2013), but this is strongly debated because baroclinic growth is decorrelated with the strength of upstream seeding (Chang and Guo, 2011, 2012). Moreover, a reduced midwinter suppression also occurs in simulations with the upstream Asian mountains removed (Park et al., 2010). The increase in the velocity of eddy propagation along the baroclinic zone and a reduced lifetime have also been shown to be insufficient to explain the suppression (Chang, 2001; Nakamura and Sampe, 2002). The reduction in the lifetime also occurs over the North Atlantic, which does not exhibit a suppression in most winters (Schemm and Schneider, 2018).
There is mounting evidence that the shift to a subtropical jet regime in the western North Pacific is essential for the formation of the midwinter suppression. During the shift to a subtropical jet regime, the subtropical thermally driven jet increases its strength, shifts equatorward, and extends into the zonal direction. Idealized studies have shown that the transition to a subtropical jet regime is able to reproduce a realistic midwinter suppression and thus ruled out the absolute necessity of zonal asymmetries, such as upstream mountains, for its formation (Yuval et al., 2018; Novak et al., 2020). The importance of a subtropical jet regime for the formation of a storm-track suppression is reinforced by the fact that a mild suppression is also observed over the North Atlantic in years with strong subtropical jets (Penny et al., 2013; Afargan and Kaspi, 2017). Chang (2001) and Nakamura and Sampe (2002) already hinted at the potential key role of the subtropical jet and its meridional displacement relative to the low-level zone of highest baroclinicity, but the exact mechanism that reduces baroclinic growth has remained unclear. Such a mechanism was suggested by Schemm and Rivière (2019), who showed that in the subtropical jet regime the ability of eddies to extract eddy energy from the mean baroclinicity is reduced because of a reduction in the baroclinic conversion efficiency. Eddies from the northern seeding branch (Chang, 2005) propagate more equatorward and towards the subtropical jet in midwinter, and during this equatorward propagation they acquire a stronger than usual poleward tilt with height, which reduces the eddy efficiency because of a weaker alignment between the mean baroclinicity and the eddy heat flux (Schemm and Rivière, 2019). Schemm and Rivière (2019) quantified baroclinic conversion for all upper-level eddies and surface cyclones in the western North Pacific. However, as shown in Sect. 4 of this study, the surface-cyclone tracks in this sector of the Pacific emerge from three different regions: (i) downstream of Kamchatka, (ii) over the Kuroshio extension, and (iii) over the East China Sea. So far, it has remained unclear whether the suppression affects the cyclones from these genesis regions in a similar way.
In this study, we investigate midwinter changes in surface-cyclone life cycles over the western North Pacific according to their genesis region. Surface cyclones are an important subcategory of the wide distribution of flow features collectively termed “eddies”. Upper-level cyclonic eddies, some of them shallow, correspond to troughs. Once a trough interacts with a surface eddy, they mutually amplify (Hoskins et al., 1985, Fig. 21) and propagate poleward in tandem (Gilet et al., 2009; Rivière et al., 2012; Oruba et al., 2013). The combined system develops into a mature low-pressure system corresponding to a deep cyclonic eddy. With surface-cyclone tracks we thus identify particularly strong cyclonic eddies that play an essential role in the overall storm-track climatology. Over the North Pacific, we expect to find different life-cycle characteristics in midwinter compared to November and March because in midwinter the cyclones typically form on the poleward flank of a strong subtropical jet, whereas in November and March they usually develop on the equatorward flank of a more poleward-located jet. Schemm and Schneider (2018) have already shown that for the entire North Pacific the lifetime of surface cyclones is shortest in midwinter. Here we study in detail all surface cyclones that affect the region of the midwinter suppression in the western North Pacific between October and April and quantify their frequency, lifetime, intensity, baroclinic conversion rates, and other characteristics according to their genesis region. This approach will serve to address the following questions.
-
What is the relative contribution of different genesis regions to the surface-cyclone frequency in the region affected by the midwinter suppression?
-
Are there any differences in the character of the surface cyclones of different origin between midwinter and the shoulder months? For example, how do their number, lifetime, and time to maximum intensity vary during the cold season?
-
Is the suppression of the baroclinic conversion during midwinter equally strong for cyclones of different origin?
To answer these questions, we use an object-based surface-cyclone tracking algorithm and evaluate baroclinic conversion rates obtained from band-pass-filtered data along individual cyclone tracks. With this approach, we combine two complementary perspectives on storm-track dynamics.
Our study is organized as follows. In Sect. 2 we introduce the data and methods used. In Sect. 3 we describe the midwinter evolution of eddy kinetic energy (EKE) over the North Pacific and define specific target regions characterized by an increase and decrease in EKE during winter, respectively. The surface-cyclone tracks and the genesis regions of cyclones that propagate through the target regions are presented in Sect. 4. Section 5 presents a detailed analysis of changes in different life-cycle characteristics. Baroclinic conversion along cyclone tracks of different origin is studied in Sect. 6. We conclude our study in Sect. 7.
The analysis period is October to April 1979–2018. All diagnostics rely on 6-hourly ERA-Interim data that are interpolated to a 1∘ grid. ERA-Interim is publicly available for download via the ECMWF archive at https://apps.ecmwf.int/datasets/ (last access: January 2021).
2.1 Surface-cyclone tracks and surface cyclogenesis
For the identification and tracking of surface cyclones, we make use of the algorithm introduced by Wernli and Schwierz (2006) and refined by Sprenger et al. (2017). The detection of surface cyclones is based on a contour search in the mean sea level pressure (SLP) field at intervals of 0.5 hPa. To obtain a cyclone mask at each time step, all grid points inside the outermost closed contour, which must not exceed 7500 km in length, are labeled with 1 and all others with 0. The obtained binary cyclone fields are used to compute cyclone frequencies. Cyclone centers are defined as the grid point with minimum SLP inside the outermost closed contour. The cyclone centers are tracked using 6-hourly cyclone center positions, and a track is accepted if it exists for a period of at least 1 d. The first time step along each track is defined as the genesis time step, and the SLP minimum defines the genesis location. The algorithm contributed to the cyclone identification and tracking intercomparison project of Neu et al. (2013).
2.2 Baroclinicity and baroclinic conversion
The background baroclinicity and the corresponding baroclinic conversion are defined based on tendency equations for eddy kinetic and available potential energy (Lorenz, 1955; Orlanski and Katzfey, 1991; Chang, 2001). A detailed derivation of both tendency equations using 10 d high-pass-filtered input data is given in Schemm and Rivière (2019). The baroclinic conversion to eddy energy, which is the sum of eddy kinetic and eddy available potential energy, is the scalar product between the eddy heat flux and the background baroclinicity :
where v′ denotes the high-pass-filtered horizontal wind, the low-pass-filtered potential temperature, and S the static stability in pressure coordinates . The reference potential temperature θR is computed from monthly mean data, and h denotes the scale height. Background baroclinicity is defined as the horizontal gradient of the low-pass-filtered potential temperature divided by the static stability, . The background baroclinicity is closely related to the Eady growth rate (Lindzen and Farrell, 1980). The scalar product that defines the baroclinic conversion can further be decomposed into contributions from the background baroclinicity and the baroclinic conversion efficiency, which is a measure for how efficiently cyclones convert the baroclinicity into eddy energy. The baroclinic conversion efficiency is connected to the vertical tilt of the growing cyclone and is maximized if the vertical tilt is such that the eddy heat flux aligns with the mean baroclinicity – for more details see Schemm and Rivière (2019). We base our analysis on the conversion rates at the 500 hPa level. This is a pragmatic choice because the suppression increases in amplitude with altitude, while baroclinic conversion is largest in the lower troposphere (Schemm and Schneider, 2018). At the 500 hPa level, the suppression of baroclinic conversion is a well-marked feature.
This section recapitulates the seasonal cycle and transition of the North Pacific storm track as seen in EKE at 500 hPa. In November, EKE increases relative to October across the entire North Pacific (red shading in Fig. 1a), with the maximum increase over the eastern North Pacific near 150∘ W. In January, EKE is reduced compared to December in an elongated band across the Pacific north of ∼43∘ (blue shading in Fig. 1b). Equatorward of this latitude, however, EKE increases (red shading in Fig. 1b). This dipole pattern in EKE tendency results from the equatorward shift of the North Pacific jet during midwinter. Absolute values of EKE (black contours in Fig. 1b) are reduced compared to those in November, which is due to the stronger EKE reduction poleward of ∼43∘ N compared to the simultaneous increase equatorward (see Fig. 1 in Schemm and Schneider, 2018, for a more detailed discussion of month-to-month EKE variations). From February to March, EKE increases again in a meridionally confined band between ∼40∘ and 50∘ N and decreases equatorward of ∼40∘. Based on theses patterns of intraseasonal changes in EKE, we select in the following surface cyclones that propagate through one of the regions with a midwinter (December to January) EKE decline or increase, respectively. The two target regions are indicated as gray boxes in Fig. 1 and are located at the entrance of the storm track (EKE maximum). One of our aims is to assess if and how the EKE tendency dipole seen in Fig. 1b is linked to characteristics of the surface-cyclone tracks. We focus in particular on cyclones that propagate through the northern target region, which exhibits a decline in EKE during midwinter and is therefore essential for understanding the midwinter suppression phenomenon.
4.1 Surface cyclogenesis
Consideration is first given to cyclogenesis associated with surface-cyclone tracks that propagate through the northern target region, i.e., the region where EKE decreases during midwinter (Fig. 1b). To this end, we extract all tracks for which the cyclone center is inside the target region for at least one time step. In November, these surface cyclones originate from two preferred regions (Fig. 2a). The first region is located over the Kuroshio extension (near 35∘ N) and the second region downstream of the Kamchatka Peninsula (near 53∘ N). This pattern and the number of events are fairly similar in January (Fig. 2c). In March (Fig. 2e), a third cyclogenesis region emerges southwest of Japan over the East China Sea (near 30∘ N), while the other two cyclogenesis regions retain fairly similar frequencies. For cyclones propagating through the target region of the midwinter suppression, there is no signal of a suppression in the genesis and therefore the number of these cyclones (Table 1).
The southern target region, i.e., the region where EKE increases during midwinter (Fig. 1b), is fed exclusively by surface cyclones with genesis over the Kuroshio extension during November (Fig. 2b). In January, a second but weaker cyclogenesis region emerges southwest of Japan over the East China Sea. In March (Fig. 2f), the two genesis regions exhibit similar cyclogenesis frequencies and contribute equally to the cyclone tracks in the southern target region. Notably, for this region there is also no midwinter suppression in the cyclone frequency. This result is in agreement with the findings of Schemm and Schneider (2018) that the suppression is connected to a reduction in the cyclones' intensity rather than frequency.
Surface cyclones over the North Pacific are known to be triggered by two upper-level seeding branches: a northern branch over Siberia and a southern branch along the subtropical jet across southern Asia (Chang, 2005). We briefly report on the upper-level seeding associated with the three preferred regions of surface cyclogenesis (Fig. 2) by means of lagged 300 hPa geopotential anomalies. We find that Kamchatka and Kuroshio cyclogenesis is triggered by waves entering the North Pacific from the northern seeding branch over Siberia (Fig. S1 in the Supplement), which in January have a more equatorward propagation direction, in agreement with Schemm and Rivière (2019). After surface cyclogenesis, the upper-level wave packet retains its overall more zonal orientation (Fig. S1). Cyclogenesis over the East China Sea is associated with the southern seeding branch and an upper-level trough downstream of the genesis location (not shown). This behavior was already recognized by Chang (2005), who noted that “cyclogenesis for these cases is probably not triggered by the (upper-level) wave packet” (Chang, 2005, p. 1998). The genesis of these cyclones seems to be connected to a bottom-up development, as is the case for diabatic Rossby waves (e.g., Boettcher and Wernli, 2013).
4.2 Relative surface-cyclone frequencies
In the previous section, we showed that surface cyclogenesis downstream of Kamchatka, over the Kuroshio extension, and in late winter over the East China Sea contributes to the surface-cyclone tracks in the northern target region, where EKE exhibits a midwinter suppression. To study the relative importance of the tracks generated in the different genesis regions for the total cyclone frequency in the target region and elsewhere, we group the cyclone tracks into two categories. The first contains tracks with genesis over the Kuroshio extension or the East China Sea, which enter the target region from the south. The second category contains tracks with genesis near Kamchatka that propagate through the target region. Next, we compute cyclone frequency fields for the two categories at every grid point and divide them by the total cyclone frequency field obtained from all cyclone tracks. The obtained relative contributions are shown for November, January, and March in Fig. 3. During all months, Kamchatka cyclones contribute up to 40 %–50 % along the poleward side of the target region and only 10 %–20 % along the equatorward side. The relative contribution of cyclone tracks entering the target region from the south is 80 %–90 % along the equatorward side of the box and decreases towards higher latitudes to 50 %–60 % along the poleward side of the target region. The relative contributions indicate the tendency of the selected cyclone tracks to propagate poleward. For example, less than 20 % of the cyclone tracks propagated from the Kuroshio region across the Pacific and into the Gulf of Alaska. Figure 3 further corroborates the fact that the midwinter suppression is related to a change in the characteristics of these cyclones and not in their frequencies. Also, we cannot focus on either Kuroshio or Kamchatka cyclones only because they both substantially contribute to the total cyclone frequency in the target region. In the following, we study statistics of cyclone characteristics for the different genesis regions in greater detail.
In the previous section, we showed that the northwestern Pacific surface storm track is fed by three preferred cyclogenesis regions: (i) the East China Sea, (ii) Kuroshio, and (iii) Kamchatka; neither of the three exhibits a midwinter suppression in terms of cyclogenesis frequency (Fig. 2). Next, we investigate several life-cycle characteristics. As in the previous sections, our focus is on the northern target region (gray box in Fig. 1), where EKE decreases during midwinter.
The lifetime from genesis to lysis of Kamchatka and Kuroshio cyclones is shortest in January and larger in November and March (Table 2). This is what one might expect from the strong midwinter jet and the fact that over the North Atlantic the shortest cyclone lifetimes are also observed during midwinter (Schemm and Schneider, 2018). However, it could also result from the fact that these poleward-propagating cyclones leave the more equatorward-located baroclinic zone earlier. Finally, the lifetime of East China Sea cyclones increases from November to March.
The time to maximum deepening since cyclogenesis is shortest during midwinter, independent of the cyclogenesis region (Table 3). Maximum deepening is measured in Bergeron, which is the 24 h change in sea level pressure along a cyclone track normalized to 60∘ N (Sanders and Gyakum, 1980). Because the cyclogenesis regions exhibit almost no variations in terms of their exact location (Fig. 2), the reduced time to maximum deepening for Kuroshio and East China Sea cyclones could result from the more equatorward location of the zone of highest baroclinicity. These poleward-propagating systems eventually leave the baroclinic zone earlier in midwinter, which explains the reduced lifetime and the shorter time to maximum deepening. For Kuroshio and East China Sea cyclones, the mean latitude where maximum deepening occurs is therefore also shifted by around 2∘ equatorward in January compared to November (not shown), which is in agreement with the earlier deepening.
With regard to the minimum sea level pressure as a measure of the storm intensity, Kamchatka cyclones become less intense from November to March (Table 4). East China Sea cyclones are most intense during midwinter, but they contribute only 22 % to the total cyclone number (Table 1). Kuroshio cyclones are also most intense in January, but the change in minimum SLP between November and January is small. The equatorward movement of the baroclinic zone in midwinter seems to be beneficial for the intensification of East China Sea cyclones in January. Kamchatka cyclones, however, become less intense, a result that can, at least qualitatively, be expected from the equatorward retreat of the baroclinic zone. For Kuroshio cyclones the situation is complex. While there is a weak reduction in the minimum SLP from November to January, the fraction of life cycles that satisfy the Sanders and Gyakum (1980) criterion for explosive deepening, known as “bomb cyclogenesis”1, is reduced in midwinter. In contrast, for the East China Sea cyclones the bomb fraction peaks in January, which is in agreement with the increase in the background baroclinicity. Kuroshio cyclones thus appear to deepen rapidly in a short time period (Table 3), in agreement with the midwinter peak in baroclinicity, but after reaching their strongest intensity they also decay rapidly, as indicated by the shortest lifetime during midwinter (Table 2). This suggests that for Kuroshio cyclones the peak in baroclinic conversion occurs earlier during their life cycle and the short deepening is also more intense, but thereafter they move out of the zone of high baroclinicity relatively soon, resulting in less intense cyclones at higher latitudes. This could also explain the dipole pattern in EKE shown in Fig. 1b because in January approximately 50 % of all cyclone tracks that propagate through the northern target region also propagate through the southern target region. Thus, to better understand the intensification, in the next section we investigate baroclinic conversion first over the two target regions from an Eulerian viewpoint and afterward along the different tracks using feature-based cyclone tracking.
6.1 Baroclinic conversion over target regions (Eulerian perspective)
EKE has a baroclinic and barotropic source, and both are known to be affected by a midwinter suppression (Schemm and Schneider, 2018). In general, however, the dominant source of EKE is baroclinic conversion. In the following, we first diagnose variations in baroclinic conversion and its link to surface cyclones in the northern and southern target regions on synoptic timescales. The daily mean values of baroclinic conversion, averaged over the northern and southern target regions (see Fig. 1), are shown in Fig. 4 for November, January, and March. Baroclinic conversion often peaks at regular intervals of 6–10 d (see the zoom-in for November 2009 as well as January and March 2010 in the right panels of Fig. 4). These baroclinic conversion bursts result mostly, but not exclusively, from the propagation of deep synoptic systems through the target region. In the northern target region in November 2009 (black contours in the top right panel of Fig. 4), the first and third bursts (labeled 1 and 3 in Fig. 4) are associated with Kuroshio cyclones, while the fourth one is associated with a Kamchatka cyclone that propagates north of the target region but still affects a broader region around it. The second burst (13–15 November) is not associated with a surface cyclone but with a jet streak development at the edge of an upper-level trough. In November and March (black and blue contours in the top right panel of Fig. 4), the amplitude of the bursts exceeds those in January by a factor of 2–3 (orange contour). The opposite is found in the southern target region, where the baroclinic conversion bursts in January exceed those in November and March. However, the monthly differences in the southern target region are smaller than in the northern target region. Furthermore, the difference between January and March is less clear in the southern target region.
The above findings suggest that the characteristics of the synoptic systems that propagate through the two target regions shown in Fig. 1 clearly differ between the three months. In January, the associated baroclinic conversion is lower than in November and March in the northern target region and vice versa for the southern region. There are several possible explanations for this behavior. For instance, (i) in midwinter, baroclinic conversion is reduced along the entire life cycle of cyclones in the northern region, or (ii) the life cycles of those cyclones that propagate through both target regions, which are about ∼50 % of all cyclones that enter the northern target region from the south, have an earlier baroclinic conversion peak (in the southern target region) and reduced baroclinic conversion later in the northern target region. As we show below, the first scenario applies to Kamchatka cyclones and the second one to Kuroshio cyclones. But first, we explore how baroclinic conversion changes during days when a cyclone propagates through the northern target region.
6.2 Baroclinic conversion in the northern target region associated with surface cyclones
In order to quantify the contribution of surface cyclones to the climatological monthly mean baroclinic conversion in the northern region, we split all days into cyclone and non-cyclone days using the surface-cyclone tracks. Thereby, we essentially separate deep eddies that extend throughout the troposphere, like mature cyclones, and shallow diabatically maintained low-level cyclonic eddies, like diabatic Rossby waves, from upper-level shallow eddies, like troughs or ridges. Anticyclonic eddies are also excluded. Technically, a surface-cyclone track may propagate outside the northern target region but nevertheless affect the baroclinic conversion inside the target region. We therefore define a cyclone day as a time step when 25 % of the northern target region is covered by a cyclone mask (see Sect. 2 for details). This results in about 50 % cyclone and 50 % non-cyclone days.
Baroclinic conversion is, as expected, larger during cyclone days compared to non-cyclone days (Fig. 5), but baroclinic conversion is not zero during non-cyclone days. Only in January is the median of the baroclinic conversion distribution near zero for non-cyclone days. The fact that baroclinic conversion is not zero during non-cyclone days can be explained by baroclinic conversion related to an upper-level trough propagating over the northern target region that is not accompanied by a surface cyclone, as is the case for the second burst in Fig. 4 (upper right panel). As clearly shown in Fig. 5, the midwinter suppression affects baroclinic conversion during cyclone and non-cyclone days. Yet, the two distributions differ significantly from each other, in particular in January during midwinter suppression. To test this statistically, we compute 10 000 distributions, each of which consists of randomly selected cyclone and non-cyclone days with replacement. Each randomized distribution is of equal size as the original cyclone day distribution. For each randomized distribution, we compute the mean baroclinic conversion, and from the 10 000 mean values we compute percentiles 97.5 and 2.5, which are shown as red confidence intervals in Fig. 5. The mean baroclinic conversion values of the cyclone (non-cyclone) day distribution is above (below) the 97.5 % (2.5 %) confidence intervals in each month. We therefore conclude that the two distributions significantly differ from each other and from a randomized selection. Based on Fig. 5, we conclude that baroclinic conversion in the target region is reduced in midwinter both during cyclone days and non-cyclone days. However, since the baroclinic conversion during cyclone days is higher than during non-cyclone days, the cyclone days contribute in absolute terms more to the total baroclinic conversion in the Pacific storm track. Nevertheless, the relative contribution to the suppression is fairly similar. The average baroclinic conversion on a cyclone day is reduced from 17 to J kg s−1 and on a non-cyclone day from 13 to J kg s−1, and the number of days in each category is close to 50 % (Fig. 5). The separation into cyclone days and non-cyclone days is not clear-cut. We defined non-cyclone days as those days during which the target region is less than 25 % covered with a cyclone mask. Baroclinic conversion on non-cyclone days might thus still be associated with surface cyclones in close proximity to the target region. In a next step, we investigate in more detail those cyclone tracks that enter the northern target region from the south and compare the baroclinic conversion along these tracks before and after entrance.
6.3 Baroclinic conversion along cyclone tracks outside and inside the northern target region (Lagrangian perspective)
Cyclones that feed the northern target region are Kamchatka cyclones, with genesis inside the target region, and Kuroshio and East China Sea cyclones, which enter the target region from the south. We therefore group all time steps along Kuroshio and East China Sea cyclone tracks into two periods, before and after entering the northern target region. The idea is to see whether the maximum baroclinic conversion occurs earlier during the life cycle in January, as suggested in Sect. 5 and based on Table 3, and therefore outside the northern target region. In the following, we discuss box-and-whisker diagrams of baroclinic conversion, the background baroclinicity (Fig. 6a), and the baroclinic conversion efficiency (Fig. 6b) separately for cyclones that enter the target region from the south (Kuroshio and East China Sea cyclones) for time steps before and after entering the target region and for Kamchatka cyclones, which reside inside the target region. For all cyclones, baroclinic conversion and its efficiency are averaged within a 1000 km radius around the cyclone center.
Before entering the northern target region from the south, baroclinic conversion along cyclone tracks is larger in January than in November and March (black boxes in Fig. 6). The distribution of baroclinic conversion outside the target region exhibits a seasonal cycle that is qualitatively in agreement with the seasonal cycle of the mean baroclinicity equatorward of the target region. The difference between January and March is small, which is a result of the increasing influence of East China Sea cyclone tracks towards late winter and early spring. East China Sea cyclones deepen on average more rapidly than Kuroshio cyclones (Table 4), and they are particularly frequent in March (Fig. 2 and Table 1).
After entering the northern target region, baroclinic conversion associated with Kuroshio and East China Sea cyclones is reduced in January compared to November and March, which reflects the midwinter suppression (gray boxes in Fig. 6). Kuroshio cyclones spend most of their life cycle in the northern target region (the percentage of time steps in each category is shown below each box in Fig. 6). In January, the fraction of time steps outside the northern target region is lower than in November and March, suggesting that in January cyclones propagate poleward faster and hence out of the zone of high baroclinicity as originally hypothesized by Nakamura (1992). On their way poleward, they also become less efficient in converting the mean baroclinicity (Fig. 6b). The reduction of baroclinic conversion in the northern target region in January occurs despite the fact that the mean baroclinicity along the tracks of Kuroshio cyclones is only marginally reduced compared to other months (horizontal red bars on top of the gray box-and-whisker diagrams in Fig. 6). Equatorward of the target region, the baroclinic conversion efficiency in January is higher compared to November, but once Kuroshio cyclones have entered the target region the conversion efficiency declines and is lower in January. In November, the mean efficiency even increases when Kuroshio cyclones enter the northern target region. The conversion budget discussed in Schemm and Rivière (2019) showed that the reduced baroclinic conversion in this region indeed results from a reduction in both the mean baroclinicity and the conversion efficiency, with the conversion efficiency making the larger contribution to the reduction. Overall, the results presented here indicate that Kuroshio cyclones in January deepen rapidly equatorward of the target region and their growth is even stronger compared with November and March, but on their way poleward their conversion efficiency decreases, in agreement with the dipole anomaly seen in Fig. 4 in Schemm and Rivière (2019). The stronger growth in January equatorward of the target region appears to accelerate the life cycle, and the cyclones seem to reach a stage when they become less efficient at converting the mean baroclinicity into eddy energy earlier during the life cycle.
For Kamchatka cyclones, baroclinic conversion is reduced in midwinter compared to the shoulder months (blue boxes in Fig. 6), in agreement with the reduced baroclinicity and reduced efficiency (Fig. 6b). The re-intensification of Kamchatka cyclones during March occurs despite no notable change in the mean baroclinicity (horizontal red bars on top of the blue box-and-whisker diagrams in Fig. 6), a finding that again points towards an increase in the baroclinic conversion efficiency as an important moderating process. Kamchatka cyclones contribute with 40 % to the storm track over the northern target region, and their weakening is therefore an important contribution to the suppression.
In summary, maximum baroclinic conversion along the surface-cyclone tracks with genesis over the Kuroshio extension is climatologically largest in January but occurs equatorward of the northern target region and therefore earlier during the cyclone life cycle. The larger conversion is thus in agreement with an overall higher mean baroclinicity over the North Pacific, but the equatorward shift of the strengthened baroclinic zone causes earlier intensification because the zone shifts towards the preferred region of Kuroshio cyclogenesis. Kuroshio cyclones not only leave the zone of highest baroclinicity faster in midwinter, but on their way poleward they also become less efficient at converting the mean baroclinicity into eddy energy. For Kamchatka cyclones, the change throughout the cold season is in fairly good agreement with what must be expected from a reduction in the background baroclinicity during midwinter and an increase in early spring.
6.4 Baroclinic conversion along cyclone tracks in a target region that shifts with the maximum monthly mean baroclinic conversion
In the previous sections, attention was given to baroclinic conversion in a target region that is centered on the location of the maximum reduction in EKE during midwinter. In this section, we explore baroclinic conversion associated with cyclone tracks that propagate through a target region that shifts with the location of the maximum monthly mean baroclinic conversion (black target box in Fig. 7). The maximum is located over the western North Pacific throughout the winter and shifts equatorward in January. The target box, which has the same latitudinal and longitudinal range in each month, is shifted latitudinally only. Baroclinic conversion decreases north of approximately 41∘ N in January and increases south of it (blue shading in Fig. 7a). Next, in line with our previous analysis, we identify all surface-cyclone tracks that propagate into this target region from upstream. Figure 7 shows the mean position of the selected tracks (red lines in Fig. 7). The black dots along this mean track indicate the mean location of cyclogenesis, the location of maximum deepening (defined as the largest 6-hourly reduction in mean SLP), the location of maximum intensity (defined as the minimum in SLP), and the mean location of cyclolysis. The first result is that the mean location at maximum deepening (second black dot along the red line) is either exactly at (November, March) or very close to the location of the maximum monthly mean baroclinic conversion in all months (black contours in Fig. 7). This underlines the fact that the selected tracks play an important role in shaping the monthly mean baroclinic conversion. The mean location of the maximum deepening is located farthest equatorward during January (Fig. 7b). We split these tracks into time steps before and after the cyclones have crossed the latitude of the maximum deepening (dashed lines in Fig. 7) and compute box-and-whisker plots of the corresponding baroclinic conversion rates. The results show that in every month, baroclinic conversion rates are higher before the maximum deepening is reached. When comparing the three months, we find that baroclinic conversion is highest during January (black boxes in Fig. 8). An important feature of the suppression is the drop in baroclinic conversion in January after maximum deepening (difference between the black and gray boxes in Fig. 8), which is larger in January compared to November and March. Consequently, we arrive at similar conclusions as in the previous section. During January, the cyclones benefit from the increased baroclinicity early during the life cycle, but due to the equatorward shift of the baroclinic zone the maximum in baroclinic conversion occurs at lower latitudes (second black dot along the red mean track in Fig. 7). The cyclones leave the baroclinic zone early on their way poleward, which can be deduced from the mean cyclone tracks in Fig. 7, and baroclinic conversion is reduced to levels below that found during November and March after the cyclones cross the latitude where maximum deepening is observed.
This study presents a systematic analysis of the characteristics of cyclone life cycles over the North Pacific with a particular focus on surface cyclones that propagate through the western North Pacific, where EKE decreases during midwinter (referred to as the northern target region shown in Fig. 1), and on surface cyclones that propagate through a target region that shifts latitudinally with the maximum in baroclinic conversion (Fig. 7). The goal of this study is to enrich the existing literature on the midwinter suppression of the North Pacific storm track with a systematic surface-cyclone life-cycle perspective to understand how cyclone life cycles change during midwinter in the western North Pacific.
The surface-cyclone tracks feeding the storm track in the western North Pacific originate from three preferred regions: (i) downstream of Kamchatka, (ii) over the Kuroshio extension, and (iii) over the East China Sea (Fig. 2). Kuroshio and Kamchatka cyclones dominate the total cyclone number over the western North Pacific throughout the cold season, while East China Sea cyclones become relevant during spring. Kamchatka and Kuroshio cyclones are preferentially triggered by upper-level waves entering the Pacific through the northern seeding branch, while East China Sea cyclones are low-level features at genesis (Chang, 2005). The analyzed tracks have their lysis mostly poleward of their genesis location in the mid-Pacific. Cyclones in the eastern North Pacific thus require further analysis.
Our key findings can be summarized as follows. The equatorward movement of the baroclinic zone in midwinter affects the life cycles of cyclones from all three genesis regions, but in a different way.
-
Kamchatka cyclones develop in midwinter in a region of reduced baroclinicity. Compared to November, their lifetime decreases, the time to maximum deepening since genesis is reduced, and they become less intense. They contribute about 40 % to the total cyclone number in winter over the western North Pacific, where EKE is suppressed. The weakening of Kamchatka cyclones is thus a crucial contribution to the suppression. Interestingly, despite the reduced baroclinicity in January, the number of Kamchatka cyclones is not reduced in midwinter. Kamchatka cyclones do not re-intensify during March, and thus they do not benefit from the poleward movement of the baroclinic zone in spring.
-
East China Sea cyclones benefit from the equatorward movement of the baroclinic zone in midwinter. Compared to November, they become more intense, the fraction of explosively deepening cyclones increases, and their lifetime increases. They become weaker in March, but the fraction of explosively deepening life cycles remains higher than for Kuroshio and Kamchatka cyclones. In addition, their lifetime is longer in March compared to January. In March, East China Sea cyclones contribute nearly 22 % to the total cyclone number over the northwestern Pacific, while in fall and winter their contribution is approximately 15 %. Thus, they seem to play a role in the re-intensification of the storm track during spring.
-
The changes in the life cycles of Kuroshio cyclones are the most complex, but understanding these changes is crucial because Kuroshio cyclones contribute the most to the total cyclone number in the northern target region in midwinter (45 %). Compared to the shoulder months, in January the lifetime of Kuroshio cyclones and the time to maximum deepening are shortest. The fraction of cyclones satisfying the “bomb cyclogenesis” criterion (Sanders and Gyakum, 1980) first decreases from November to January but then remains at similar levels until March. However, the highest values in 6-hourly baroclinic conversion rates are found in January, but these occur at lower latitudes and they are sustained for a reduced number of time steps relative to the shoulder months. In terms of minimum sea level pressure, Kuroshio cyclones are, however, most intense in January.
Overall, during midwinter, the life cycle of a Kuroshio cyclone is best characterized by a short and intense early deepening, in agreement with the higher baroclinicity, followed by a fast decay and poleward propagation away from the more equatorward-located baroclinic zone. According to this interpretation, we observe an acceleration of the Kuroshio life cycle during midwinter. This interpretation is in agreement with the idea of a reduced baroclinic conversion efficiency because the efficiency is dictated by the vertical tilt of a cyclone (Schemm and Rivière, 2019). Acceleration of the life cycle with intense early growth results in cyclones that acquire a rather inefficient vertical tilt earlier in the life cycle. Kuroshio cyclones are thus in different stages of their life cycle in different months at similar latitudes. The stronger but earlier deepening followed by an earlier decay is the cyclone life-cycle perspective on the midwinter suppression over the western North Pacific.
It is also important to mention a few caveats of this study. Our results are based on a single object-based cyclone detection scheme. It is known that cyclone tracks are sensitive to the identification and tracking scheme (Neu et al., 2013); this also holds true for the genesis location. While they typically agree on deep systems, a higher sensitivity must be expected for shallow systems, as is the case for Kamchatka cyclones in January. Further, we ignore short-lived systems with a lifetime of less than 24 h. Such systems might become more frequent in midwinter because the lifetime of all systems is reduced in January. Baroclinic conversion also occurs in the absence of surface cyclones, for example through the propagation of an upper-level trough, like the second baroclinic conversion peak in Fig. 4. Baroclinic conversion during non-cyclone days is also affected by the midwinter suppression, and this reduction is not explained by our study. Our study is focused on the western North Pacific, where climatological mean baroclinicity is highest and the reduction in baroclinic conversion and EKE is thus most surprising. Cyclone tracks that feed the eastern North Pacific are generated over the central eastern Pacific (Hoskins and Hodges, 2002; Wernli and Schwierz, 2006) where climatologically baroclinicity is much lower compared with the western Pacific. Cyclone tracks over the eastern North Pacific not surprisingly have a large fraction of secondary cyclones (Fig. 5b in Schemm et al., 2018). Mechanisms responsible for the suppression in the eastern North Pacific thus require further analysis.
ERA-Interim is publicly available for download via the ECMWF archive at https://apps.ecmwf.int/datasets/ (last access: January 2021) (ECMWF, 2021). The cyclone track data is made available on request.
The supplement related to this article is available online at: https://doi.org/10.5194/wcd-2-55-2021-supplement.
All three authors contributed to the discussion and final interpretation of the results. HB and SeS performed the analyses. SeS wrote the publication, supported by HB and HW.
The authors declare that they have no conflict of interest.
The authors acknowledge the discussions on midwinter suppression with various colleagues at the 2018 storm-track workshop in Utö, Sweden, and thank two anonymous reviewers for their constructive comments that helped to improve this paper. We thank MeteoSwiss and the ECMWF for access to the ERA-Interim reanalyses.
This research has been supported by the European Research Council (grant no. GLAD (848698)) and the Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (grant no. 185049).
This paper was edited by David Battisti and reviewed by two anonymous referees.
Afargan, H. and Kaspi, Y.: A midwinter minimum in North Atlantic storm track intensity in years of a strong jet, Geophys. Res. Lett., 44, 12511–12518, https://doi.org/10.1002/2017GL075136, 2017. a
Boettcher, M. and Wernli, H.: A 10-yr climatology of diabatic Rossby waves in the Northern Hemisphere, Mon. Weather Rev., 141, 1139–1154, https://doi.org/10.1175/MWR-D-12-00012.1, 2013. a
Chang, E. K. M.: 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, 2001. a, b, c
Chang, E. K. M.: The impact of wave packets propagating across Asia on Pacific cyclone development, Mon. Weather Rev., 133, 1998–2015, https://doi.org/10.1175/MWR2953.1, 2005. a, b, c, d, e
Chang, E. K. M. and Guo, Y.: Comments on the source of the midwinter suppression in storminess over the North Pacific, J. Climate, 24, 5187–5191, https://doi.org/10.1175/2011JCLI3987.1, 2011. a
Chang, E. K. M. and Guo, Y.: 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, 2012. a
Deng, Y. and Mak, M.: 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, 2005. a
ECMWF: Public Datasets, available at: https://apps.ecmwf.int/datasets/, last access: January 2021. a
Gilet, J.-B., Plu, M., and Riviẽre, G.: Nonlinear baroclinic dynamics of surface cyclones crossing a zonal jet, J. Atmos. Sci., 66, 3021–3041, https://doi.org/10.1175/2009JAS3086.1, 2009. a
Harnik, N. and Chang, E. K. M.: 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, 2004. a
Hoskins, B. J. and Hodges, K. I.: New Perspectives on the Northern Hemisphere winter storm tracks, J. Atmos. Sci., 59, 1041–1061, https://doi.org/10.1175/1520-0469(2002)059<1041:NPOTNH>2.0.CO;2, 2002. a
Hoskins, B. J., McIntyre, M. E., and Robertson, A. W.: On the use and significance of isentropic potential vorticity maps, Q. J. Roy. Meteorol. Soc., 111, 877–946, 1985. a
James, I. N.: 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, 1987. a
Lindzen, R. S. and Farrell, B.: A simple approximate result for the maximum growth rate of baroclinic instabilities, J. Atmos. Sci., 37, 1648–1654, https://doi.org/10.1175/1520-0469(1980)037<1648:ASARFT>2.0.CO;2, 1980. a
Lorenz, E. N.: Available potential energy and the maintenance of the general circulation, Tellus, 7, 157–167, https://doi.org/10.3402/tellusa.v7i2.8796, 1955. a
Nakamura, H.: 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, 1992. a, b, c
Nakamura, H. and Sampe, T.: Trapping of synoptic-scale disturbances into the North-Pacific subtropical jet core in midwinter, Geophys. Res. Lett., 29, 8-1–8-4, https://doi.org/10.1029/2002GL015535, 2002. a, b
Nakamura, N.: An illustrative mmdel of instabilities in meridionally and vertically sheared flows, J. Atmos. Sci., 50, 357–376, https://doi.org/10.1175/1520-0469(1993)050<0357:AIMOII>2.0.CO;2, 1993. a
Neu, U., Akperov, M. G., Bellenbaum, N., Benestad, R., Blender, R., Caballero, R., Cocozza, A., Dacre, H. F., Feng, Y., Fraedrich, K., Grieger, J., Gulev, S., Hanley, J., Hewson, T., Inatsu, M., Keay, K., Kew, S. F., Kindem, I., Leckebusch, G. C., Liberato, M. L. R., Lionello, P., Mokhov, I. I., Pinto, J. G., Raible, C. C., Reale, M., Rudeva, I., Schuster, M., Simmonds, I., Sinclair, M., Sprenger, M., Tilinina, N. D., Trigo, I. F., Ulbrich, S., Ulbrich, U., Wang, X. L., and Wernli, H.: IMILAST: A community effort to intercompare extratropical cyclone detection and tracking algorithms, B. Am. Meteorol. Soc., 94, 529–547, https://doi.org/10.1175/BAMS-D-11-00154.1, 2013. a, b
Novak, L., Schneider, T., and Ait-Chaalal, F.: 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, 2020. a, b
Orlanski, I. and Katzfey, J.: The life cycle of a cyclone wave in the Southern Hemisphere. Part I: Eddy energy budget, J. Atmos. Sci., 48, 1972–1998, https://doi.org/10.1175/1520-0469(1991)048<1972:TLCOAC>2.0.CO;2, 1991. a
Oruba, L., Lapeyre, G., and Rivière, G.: On the poleward motion of midlatitude cyclones in a baroclinic meandering jet, J. Atmos. Sci., 70, 2629–2649, https://doi.org/10.1175/JAS-D-12-0341.1, 2013. a
Park, H.-S., Chiang, J. C. H., and Son, S.-W.: 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, 2010. a
Penny, S., Roe, G. H., and Battisti, D. S.: The source of the midwinter suppression in storminess over the North Pacific, J. Climate, 23, 634–648, https://doi.org/10.1175/2009JCLI2904.1, 2010. a
Penny, S. M., Roe, G. H., and Battisti, D. S.: Reply, J. Climate, 24, 5192–5194, https://doi.org/10.1175/2011JCLI4187.1, 2011. a
Penny, S. M., Battisti, D. S., and Roe, G. H.: 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, 2013. a, b
Rivière, G., Arbogast, P., Lapeyre, G., and Maynard, K.: A potential vorticity perspective on the motion of a mid-latitude winter storm, Geophys. Res. Lett., 39, L12808, https://doi.org/10.1029/2012GL052440, 2012. a
Sanders, F. and Gyakum, J. R.: Synoptic-dynamic climatology of the “bomb”, Mon. Weather Rev., 108, 1589–1606, https://doi.org/10.1175/1520-0493(1980)108<1589:SDCOT>2.0.CO;2, 1980. a, b, c, d
Schemm, S. and Rivière, G.: 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, 2019. a, b, c, d, e, f, g, h, i
Schemm, S. and Schneider, T.: 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, 2018. a, b, c, d, e, f, g, h, i
Schemm, S., Sprenger, M., and Wernli, H.: When during their life cycle are extratropical cyclones attended by fronts?, B. Am. Meteorol. Soc., 99, 149–165, https://doi.org/10.1175/BAMS-D-16-0261.1, 2018. a
Sprenger, M., Fragkoulidis, G., Binder, H., Croci-Maspoli, M., Graf, P., Grams, C. M., Knippertz, P., Madonna, E., Schemm, S., Škerlak, B., and Wernli, H.: Global climatologies of Eulerian and Lagrangian flow features based on ERA-Interim, B. Am. Meteorol. Soc., 98, 1739–1748, https://doi.org/10.1175/BAMS-D-15-00299.1, 2017. a
Wernli, H. and Schwierz, C.: Surface Cyclones in the ERA-40 Dataset (1958–2001). Part I: Novel identification method and global climatology, J. Atmos. Sci., 63, 2486–2507, https://doi.org/10.1175/JAS3766.1, 2006. a, b
Yuval, J., Afargan, H., and Kaspi, Y.: The relation between the seasonal changes in jet characteristics and the Pacific midwinter minimum in eddy activity, Geophys. Res. Lett., 45, 9995–10002, https://doi.org/10.1029/2018GL078678, 2018. a
A change in SLP larger than 24 hPa within 24 h normalized to 60∘ N (Sanders and Gyakum, 1980).
- Abstract
- Introduction
- Data and methods
- EKE of the North Pacific storm track in midwinter
- Surface-cyclone view on the North Pacific storm track in midwinter
- Detailed characteristics of surface-cyclone life cycle
- Baroclinic conversion and its relationship with surface-cyclone tracks
- Conclusions
- Data availability
- Author contributions
- Competing interests
- Acknowledgements
- Financial support
- Review statement
- References
- Supplement
- Abstract
- Introduction
- Data and methods
- EKE of the North Pacific storm track in midwinter
- Surface-cyclone view on the North Pacific storm track in midwinter
- Detailed characteristics of surface-cyclone life cycle
- Baroclinic conversion and its relationship with surface-cyclone tracks
- Conclusions
- Data availability
- Author contributions
- Competing interests
- Acknowledgements
- Financial support
- Review statement
- References
- Supplement