- 1Department of Oceanography and Coastal Sciences, College of the Coast and Environment, Louisiana State University, Baton Rouge, LA, United States
- 2Coastal Studies Institute, Louisiana State University, Baton Rouge, LA, United States
The tropospheric circumpolar vortex (CPV), an important signature of processes steering the general atmospheric circulation, surrounds each pole and is linked to the surface weather conditions. The CPV can be characterized by its area and circularity ratio (Rc), which both vary temporally. This research advances previous work identifying the daily 500-hPa Northern Hemispheric CPV (NHCPV) area, Rc, and temporal trends in its centroid by examining linear trends and periodic cycles in NHCPV area and Rc (1979–2017). Results suggest that NHCPV area has increased linearly over time. However, a more representative signal of the planetary warming may be the temporally weakening gradient which has blurred NHCPV distinctiveness—perhaps a new indicator of Arctic amplification. Rc displays opposing trends in subperiods and an insignificant overall trend. Distinct annual and semiannual cycles exist for area and Rc over all subperiods. These features of NHCPV change over time may impact surface weather/climate.
Introduction
The tropospheric circumpolar vortices (CPVs, Waugh et al., 2017) are the two hemispheric-scale, extratropical circulations that circumnavigate each polar region—one in the Northern Hemisphere (NHCPV) and the other in the Southern. The CPVs are positioned near the albeit discontinuous surface polar front at a given time, situated near the steepest geopotential height gradient, the sharpest gradient of air temperature, and the fastest flow of the quasi-westerly wind belt, including the polar front jet stream. Spatial features of the NHCPV area and circularity were examined by Bushra and Rohli (2019, 2021a), and the spatial and temporal features of NHCPV centroid position were described by Bushra and Rohli (2021b).
The amplitude and length of poleward (i.e., ridges in the Northern Hemisphere) and equatorward (i.e., troughs in the Northern Hemisphere) meanders in the Rossby waves in the CPV's westerly flow support surface storminess (Di Capua and Coumou, 2016) through upper-level divergence (ULD) and cyclonic vorticity advection (CVA), and suppress surface storminess by upper-level convergence (ULC) and anticyclonic vorticity advection (AVA). Likewise, ULC and AVA (ULD and CVA), driven by the Rossby waves, support (suppress) surface anticyclones. While the changes in the CPV's centroid location and dimensions (i.e., area and wave amplitudes in the form of ridges/troughs) are likely to affect the location and degree to which surface storms or anticyclones beneath the CPV are supported, surface features may also influence the CPV properties (Angell and Korshover, 1985; Kirby et al., 2002). Variabilities in the strength, position, wave amplitudes, and location of the CPV are well-known to be linked directly to variations in surface environmental features such as surface air temperature and wind (van den Broeke and van Lipzig, 2002; Moron et al., 2018), precipitation (Srinivas et al., 2018), sea surface temperature (Frauenfeld et al., 2005), ocean salinity (Chen et al., 2018), water vapor transport (Wang and Ding, 2009), air mass properties (Vanos and Cakmak, 2014), storm tracks (Kidston et al., 2015), ozone (Glovin et al., 2016), sea ice extent (Orme et al., 2017), aerosol distributions (Bartlett et al., 2018), and atmospheric blocking (Altenhoff et al., 2008; Tyrlis and Hoskins, 2008), just as Earth surface topography is associated directly to the steering circulation (Tang and Chan, 2016).
In an early study, Angell and Korshover (1977) noted that the 300-hPa annual-averaged NHCPV area increased from 1970 to 1975 despite simultaneously shrinking over time in winter. Angell (1998, 2006) later suggested that the 300-hPa-NHCPV area was decreasing over the 1963–2001 period, especially in the Western Hemisphere, at a rate of 1.5% per decade. Angell (1998) demonstrated that the most intense 300-hPa-NHCPV contraction began around 1990. By contrast, Davis and Benkovic (1992, 1994) noted an expansion of the 500-hPa January NHCPV from 1966 to 1990, primarily resulting from intensified regional troughing over the North Pacific Ocean and eastern North America and concurrent warm air mass intrusion into Alaska and western Canada. Burnett (1993) found similar results for the January NHCPV after the mid-1960s, with areal growth due to expansion along the Pacific and eastern North America/Atlantic sectors, but with contraction commencing in the late 1980s. Frauenfeld and Davis (2003) investigated temporal trends of NHCPV area at multiple constant-pressure levels over the 1949–2000 period; at all levels they found significant NHCPV expansion until 1970 and areal contraction afterward, with the primary areas of expansion/contraction over Asia, Europe, and North America, and less variability over the Northern Hemisphere oceans. Trends were stronger in the upper troposphere than near the surface (Frauenfeld and Davis, 2003), likely because of the complexity of the surface influence. Using threshold temperatures at the 850 hPa level as the indicator, Martin (2015) observed a shrinking of the Northern Hemisphere cold pool over time.
While all of the above research delineated the NHCPV area using a pre-determined isohypse, as recommended by Frauenfeld and Davis (2003), Thompson and Solomon (2002) applied vertical gradients of geopotential heights at the relatively coarse network of radiosonde-based stations to define the NHCPV. Bushra and Rohli (2019) delineated the CPV's leading edge at the sharpest latitudinal gradient of 500-hPa geopotential height.
Although the NHCPV's area has been studied for several decades, analysis of CPV shape had been largely restricted to reporting the relative area within various quadrants about the pole until Rohli et al. (2005), Wrona and Rohli (2007), and Bushra and Rohli (2019) borrowed from watershed hydrologists' drainage basin metrics (Chorley et al., 1984) by using a “circularity ratio” (Rc), ranging from 0 (i.e., infinitely amplified ridges and troughs) to 1.0 (i.e., perfectly circular, west-to-east-flowing CPV). CPV sinuosity (i.e., ratio of CPV length to the length of the corresponding parallel of latitude segment) has also been analyzed for the 500-hPa level in the North American sector (Vavrus et al., 2017) and at the hemispheric scale (Cattiaux et al., 2016; Di Capua and Coumou, 2016), using a fast Fourier transform (FFT) approach (Screen and Simmonds, 2013), and for potential-vorticity-based identification of Rossby wave initiation (Röthlisberger et al., 2016). However, despite these and a few other important studies (e.g., Zakinyan et al., 2016), there remains a need for examining and interpreting the CPV shape.
Analyzing temporal trends in the NHCPV and its Rc may reveal important features of hemispheric-scale response to surface warming. While air temperatures are increasing globally, the warming trend of the past several decades may be most pronounced in the high latitudes, especially in the Arctic (Holland and Bitz, 2003; Pithan and Mauritsen, 2014; Post et al., 2019). This “Arctic amplification” decreases the lateral temperature gradient and affects the tropics-to-poles transfer of energy in both the atmosphere (e.g., Cohen et al., 2014) and ocean (e.g., Rahmstorf et al., 2015; Praetorius, 2018). Arctic amplification has been proposed to weaken the zonal winds, which is connected to enhanced Rossby wave ridging and troughing and support for severe weather by increasing areas of ULD and CVA. Model-based research has posited that these meridional meanders will increase into the future (Francis and Vavrus, 2012, 2015; Peings and Magnusdottir, 2014), but Cattiaux et al. (2016) disagreed, with Vavrus (2018) summarizing the conflicting evidence and lamenting the lack of consensus.
Cyclical variation in the broad-scale steering circulation has included analysis of the energy distribution of the atmospheric motion and waves. This oscillation can be explained through their linearized/non-linearized energy transfer spectra with Fourier analysis over space and time. Understanding of the relationship between wave number and frequency, and variability of magnitude and phase could aid understanding of the dynamics of the energetics influencing surface weather. Hayashi (1974) performed space-time cross spectrum analysis to understand the effect of tropical disturbances on the general circulation. Chen et al. (1981) simulated winter wind fields over a 90-day period by applying modified Fourier analysis to characterize the wavenumber-frequency spectra. Hayashi (1982) used space-time spectral analysis to understand broad-scale atmospheric waves. Higuchi et al. (1999) employed multiresolution Fourier transform spectral analysis to resolve the temporal structure of the variation of the North Atlantic Oscillation (NAO; van Loon and Rogers, 1978; Wallace and Gutzler, 1981; Barnston and Livezey, 1987; Hurrell, 1995), a component of the NHCPV, in terms of various frequency components. Torrence and Compo (1998) explained the implications of wavelet analysis in the atmospheric sciences using the example of the El Niño–Southern Oscillation (ENSO; Philander, 1990) time series. Pelletier (1998) explored local variability of the high-frequency component of the temperature time series using a power spectrum. Ghil et al. (2002) applied time series analysis by treating the atmosphere as a non-linear dynamic system. Rodionov (2004) applied sequential algorithm in time series data to detect regime shifts and track magnitude changes. More recent research has included the examination of the daily 500-hPa geopotential height data of the Southern Hemisphere to analyze the space-time frequency characteristics of the expansion coefficients (Sun and Li, 2012).
A comprehensive analysis of the spatiotemporal cyclical nature of the NHCPV area and Rc remains lacking. Implications of the linear regression model represent the temporal trend over the time domain while the FFT extracts information on the variability of the magnitude and oscillations by analyzing the frequency variability.
Purpose
The purpose of this research is to evaluate temporal variations of the NHCPV's area and Rc at the 500-hPa geopotential height level. Time-series analysis will identify the periodicity/frequency characteristics and variability over time. Analysis of the intra- and inter-annual variabilities of the NHCPV area and Rc provides an index of the relationship between the atmospheric steering flow and surface warming, and it also contextualizes recent results (Bushra and Rohli, 2019). Analyzing the intra- and inter-annual variability of the NHCPV shape is important not only as a potential characterization of the response to Arctic amplification, but also as a possible indicator of baroclinic instability associated with severe extratropical weather in future work. Analysis of NHCPV area and Rc also informs studies on extratropical severe weather that rely on the juxtaposition of tropical and polar air masses (e.g., Overland et al., 2015).
Data and methods
To facilitate comparison of results to Bushra and Rohli (2019, 2021a,b), daily 500-hPa geopotential height data are obtained from the National Centers for Environmental Information (NCEI)/National Center for Atmospheric Research (NCAR) Reanalysis 2 project (Kanamitsu et al., 2002). The spatial resolution of this data set allows for delineation of longwave flow while filtering out shortwave flow, which assists in overcoming the computational challenges from daily CPV shapes perforated by cutoff lows, omega blocks, and Rex blocks (Lupo, 2021). At each 2.5° increment of longitude globally from 20°N to the North Pole, the daily 500-hPa NHCPV is delineated over the 1979–2017 period using the method described more fully in Bushra and Rohli (2019). Specifically, the NHCPV polygon is formed by connecting the dots representing the latitude of the steepest gradient of 500-hPa geopotential height for each meridian of longitude around the Northern Hemisphere. This approach circumvents the tenuous assumption that the leading edge of the CPV follows the same isohypse across space, despite differences in topographic, marine, and air mass influence. After Bézier smoothing (Jekeli, 2005), the NHCPV area is then measured, and Rc is calculated as
where A represents the area enclosed within the NHCPV and Ac denotes the area of the circle with the same perimeter as the NHCPV on that day. Because the calculation of segment length makes Rc resolution-dependent, it is important that the same data used in earlier calculations of NHCPV properties (Bushra and Rohli, 2019, 2021a,b) are used here.
Statistical regression is then used to analyze the linear trend of the NHCPV's area and Rc for the 500-hPa level for the entire time period and subperiods identified via breakpoint analysis (Bai and Perron, 2003; Muggeo, 2008). The number of breakpoints identified depends on the testing and deviations from stability in the linear regression models, where the coefficients shift from one stable regression relationship to a different one (Bai, 1994, 1997a,b; Bai and Perron, 1998). If there are m breakpoints, the number of segments in which the regression coefficients are constant will be m+1.
To remove high-frequency noise, a Butterworth (Butterworth, 1930) low-pass filter is applied at a cutoff point of 3.65 year−1 to the dataset, before applying FFT (Welch, 1967). Then, Fisher's Exact G Test for Multiple (Genetic) Time Series, or “fisher.g.test” (Ahdesmaki et al., 2021), is applied to the NHCPV's area and Rc data for the entire record and segmented time periods resulting from the breakpoint analysis to confirm that remaining noise is not problematic. Any p-values of < 0.001 support rejection of the null hypothesis of significant presence of white noise.
The presence of lag-1 autocorrelation between two consecutive observations is indicative of red noise (Schulz and Mudelsee, 2002) and has a power spectrum weighted toward low frequencies (Ghil et al., 2002). Application of the Durbin-Watson (D-W) test (Savin and White, 1977) with the test statistic of <2 (p-value < 0.001) for all data series confirms the presence of lag-1 autocorrelation in the data sets. To confirm that the resultant power spectra are not from red noise, the Cochrane–Orcutt (C-O) method (Montgomery et al., 2012) is used to remove any such autocorrelation to make the data sets independent in time (Vanhatalo and Kulahci, 2016). Again, application of D-W test with a test statistic of >2 (p-value < 0.001) suggests that lag-1 autocorrelation (i.e., red noise) is removed successfully from each data sets.
Finally, FFT is applied for spectral analysis. Spectral analysis identifies the space-time Fourier component of the NHCPV properties to evaluate the periodic characteristics of the area and Rc time series data for all time periods by decomposing their time series signal into the frequency domain (i.e., quantifying their magnitude of variations at different frequencies).
Results and discussion
Long-term trends
Linear regression analysis reveals that over the 1979–2017 time period, the area enclosed within the 500-hPa NHCPV shows a significant temporal increase (125 km2 day−1; p < 0.001; Figure 1A) despite coincident increases in global mean surface temperature (Folland et al., 2018). A warming world might have been expected to be characterized by a shrinking cold pool as represented by the NHCPV. Time series breakpoints are identified at a 95% confidence level in 1986, 1994, and 2001, with NHCPV area increasing significantly within all four segmented periods (Figure 1B), but with the 1979–1986 segment showing an increasing area only at a 90% confidence level (Table 1). The 95% confidence interval of each breakpoint location is shown in Figure 1B and Table 1. The greatest rate of NHCPV areal increase (1,305 km2 day−1) occurred over the 1994–2001 period, perhaps because of cooling at the end of the extended El Niño event in the early 1990s, which enlarges the Pacific sector of the NHCPV while shrinking the American sector (Frauenfeld and Davis, 2000). The eruption of Mount Pinatubo in the Philippines apparently was not an important factor, as it has been shown to be associated with shrinking, while strengthening, the NHCPV (Otterå, 2008). The lowest rate of 188 km2 day−1 occurred over the 2001–2017 period, possibly corresponding to shift in the Pacific Decadal Oscillation (PDO; Newman et al., 2003) to a cold phase (National Oceanic Atmospheric Administration, 2021), which might partially offset the tendency for NHCPV expansion, particularly in the Pacific basin. Both of these hypotheses invite further investigation, but the time series is not sufficiently long to identify multiple cycles of influence from either the Atlantic Multidecadal Oscillation (Kerr, 2000) or the PDO (Rohli et al., 2022). Note that because of the abrupt downward shift in area at the breakpoints, the minimum rate of areal growth in subperiods still exceeds the mean temporal increase over the entire time series.
Figure 1. (A) Statistically significant (p < 0.001) linearly increasing trend in daily 500-hPa NHCPV area over the 1979–2017 period. The smoothed black line, from Butterworth low-pass filtering, shows the annual cycle. (B) Linear regressions from breakpoint analysis of four subperiods of 1979–1986, 1986–1994, 1994–2001, and 2001–2017 for NHCPV area (p-values shown in Table 1) with the 95% confidence interval for the breakpoint positions shown along the x-axis and in Table 1. (C) As in (A) but showing statistically insignificant (p = 0.445) linearly decreasing trend in 500-hPa Rc. (D) As in (B), but for four subperiods of 1979–1988, 1988–1994, 1994–2001, and 2001–2017 for NHCPV Rc. (E) As in (A), but showing statistically significant (p = 0.001) linearly decreasing trend in 500-hPa mean maximum daily north-south 500-hPa geopotential height gradient. (F) As in (B), but for four subperiods of 1979–1986 (p = 0.016), 1986–1993 (p = 0.718), 1993–2004 (p < 0.001), and 2004–2017 (p < 0.001) for NHCPV gradients.
Table 1. Regression line with significance of slope, and 95% confidence interval of breakpoints, for NHCPV area and Rc.
Another possible explanation for the temporal increase in NHCPV area may be that in the warming world during these decades, the leading edge of the NHCPV may be becoming less distinct, leading to cases in which the subtropical jet (STJ) could possibly mistakenly be delineating the CPV, leading to an increasingly larger CPV area. This is a particular concern in the summer months, because the consideration of data points only poleward of 20°N would limit the possibilities for inclusion of the STJ in the cold season. The large areal extent of a STJ-based NHCPV of the warm season may overcompensate for the shrinking polar-front-based, cold-season NHCPV, leading to net NHCPV expansion.
Linear regression analysis of the 500-hPa NHCPV's Rc reveals no significant trend (−4.258 × 10−8 day−1; p = 0.445) over the entire time series (Figure 1C). Four time series segments are also identified for the 500-hPa NHCPV Rc, with breakpoints in 1988, 1994, and 2001 and the 95% confidence interval of each breakpoint location shown in Figure 1D and Table 1. These breakpoints are similar to those for the area analysis, probably suggesting that the trend of the undulation of the ridges and troughs is related to the NHCPV area, and that regime shifts for NHCPV area and circularity happen simultaneously.
In the 1979–1988 and 1994–2001 segments, Rc increases significantly (at 90 and 95% confidence levels, respectively) with the highest increasing rate of 0.02 × 10−8 day−1 over the latter period. However, these subperiods are interrupted by intervals of decreasing trends for 1988–1994 and 2001–2017 (at 95 and 90% confidence levels, respectively) with the greatest decreasing rate of −0.017 × 10−8 day−1 over the earlier period.
The apparently fluctuating but overall temporally increasing area of the CPV amid temporally warming conditions could leave less volume for the remaining cold air pool, leading to increasing circularity. It could also be related to the fact that, in contrast to the cold pool analysis of Martin (2015) which only considers DJF, the analysis here considers the entire year. But the trend identified here could be counteracted by a tendency for increased wave amplitudes amid a decreasing equator-to-pole temperature gradient, as observed by Martin (2021), as the forces for geostrophic flow weaken, in a manner analogous to the increased meandering of a stream amid a gentler gradient near its mouth. The net result would be a relationship between area and circularity but with no overall linear temporal trend in Rc despite a temporally decreasing NHCPV area.
This result suggests that the 500-hPa NHCPV's shape may be characterized by alternating increases and decreases in Rossby wave amplification over time. Decreasing circularity may be expected to occur in an environment of weakening equator-to-pole gradients of 500-hPa geopotential height and weakening zonality of winds, which may in turn amplify the Rossby waves (Francis and Vavrus, 2012), providing for increasing frequency and magnitude of weather extremes. But this pattern is “reset” eventually, as a trend for increasing circularity sets in for a period of time before the former pattern is re-established. The areal expansion contrasts with the notion that the 500-hPa NHCPV area would decrease over time as the cold pool retreats poleward in response to surface warming and instead suggests that other factors may be driving the result.
Therefore, in a follow-up analysis, it is confirmed that the daily-averaged 500-hPa geopotential height gradient that demarcates the leading edge of the NHCPV is indeed decreasing significantly, at a rate of 0.1661 m km−1 day−1 (p < 0.001) over the study period (Figure 1E). This value is calculated from the trendline (y = 89.44–1.661 × 10−4 (note here that the units of the slope of the equation are m m−1 day−1); p < 0.001) of the 14,245 daily mean values of the 144 daily longitudinal maximum 500-hPa daily geopotential height difference (in m) at a given longitude occurring across the horizontal distance between two adjacent grid points 2.5° of latitude (or 279.23 km) apart. This weakening gradient may be a stronger symptom of Arctic amplification than the NHCPV area. Decreasing trends within each of the four subperiods of the study period are also apparent but with different levels of statistical significance (Figure 1F).
Seasonal cycle
The power spectrum allows an assessment of the frequency and amplitudes of the 500-hPa NHCPV characteristics. The NHCPV area has a distinct annual cycle suggestive of the variation of solar radiation over each year in the Northern Hemisphere due to the tilt of Earth's axis of revolution and a semiannual cycle perhaps reflective of the cold-warm seasonal flow, over the complete study period (1979–2017). However, the magnitude of the annual cycle varies over the subperiods with an indistinct or weaker semiannual cycle (Figures 2A–E). The semiannual signal among these daily standardized values suggests that z-scores of the area might offer at least some degree of understanding for the z-scores of area 6 months into the future. FFT also reveals a high frequency intra-annual cycle of ~7.5 year−1 but only over the 1994–2001 period (Figure 2D). This frequency matches that of the 40–60-day oscillation, also known as the Madden-Julian Oscillation (MJO; Madden and Julian, 1994; Henderson et al., 2016; Gollan and Greatbatch, 2017), which is a tropical eastward-propagating sequence of convective activity that was identified via zonal wind anomalies in the tropical Pacific Ocean (Madden and Julian, 1971). Given the demonstrated link between tropical and extratropical circulation anomalies (e.g., Lau et al., 2012), a link between the MJO and NHCPV is a topic that warrants future research. The 1994–2001 subperiod showed the strongest increasing trend in NHCPV area.
Figure 2. Magnitudes of the power spectra of the daily NHCPV area over (A) 1979–2017, (B) 1979–1986, (C) 1986–1994, (D) 1994–2001, and (E) 2001–2017 periods. The peaks indicate the magnitude of the high-frequency variability. The highest peaks at 1 cycle year−1 indicate an annual cycle for all periods.
FFT analysis for the NHCPV Rc also shows a distinct annual cycle over the 1979–2017 period, but with less evidence of a semiannual cycle (Figures 3A–E). The annual signal has similar strength over the 1979–1988, 1988–1994, and 1994–2001 subperiods (Figures 3B–D), while the magnitude of 2001–2017 (Figure 3E) is greater. The smaller amplitudes in the earlier stage of the time series for the three subperiods imply that the shape of the CPV had less association to the annual cycle of solar radiation than in the 2001–2017 subperiod (i.e., the warmest decades). In addition, the 1979–1988 and 2001–2017 periods (Figures 3B, E) show a secondary cycle at 3 year−1.
Figure 3. As in Figure 2, but for Rc over (A) 1979–2017, (B) 1979–1988, (C) 1988–1994, (D) 1994–2001, and (E) 2001–2017 periods.
The reason for this cycle for only these two subperiods is unclear and deserves future research consideration. However, the NAO and the related Arctic Oscillation (Thompson and Wallace, 1998, 2000) are likely responses to the forcing mechanisms as they have been shown to be closely linked to the NHCPV, even at the daily time scale (Bushra and Rohli, 2021a).
Summary and conclusions
Analysis of daily variability of area and circularity of the NHCPV defined by the sharpest 500-hPa geopotential height gradient offers potential for improved understanding of the mid-latitude circulation. In this research, long-term changes in the area and Rc of the 500-hPa NHCPV are presented for the 1979–2017 period of record.
Linear regression analysis suggests that the NHCPV area has increased significantly over time, including temporal increases over each subperiod of the time series as revealed by breakpoint analysis. Concurrent with the areal change is a temporally significant trend toward a reduction in the 500-hPa geopotential height gradient that separates the polar cold pool from the temperate atmosphere and delineates the leading edge of the NHCPV. Thus, the distinctiveness of the NHCPV blurs over time, which may impact the trends and the NHCPV's covariation with surface weather. No linear change in Rc was observed over 1979–2017 periods, but the four subperiods are similar to those for the area analysis, suggesting that the regimes that govern area also govern circularity. Adjacent subperiods of Rc show reversing trends and different significance levels (i.e., 90–95%).
Application of the power spectrum to a low pass filtered time series reveals that both NHCPV area and Rc have distinct annual cycles but with varying amplitude. FFT also identifies intra-annual cycles for subperiods.
Collectively, these results are important indicators of the association of global temperature change at varying time scales with the upper-level steering circulation. Analyzing CPV characteristics reveals important information about how hemispheric circulation may be related to the surface environment and may have associations to extreme atmospheric events. Analysis of the intra-annual variability of the NHCPV area provides a tangible index of the link between atmospheric steering flow and the observed surface warming.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
NB developed the detailed methodology, collected and analyzed the data, and developed the initial text. RR helped develop the methodology and edited early and late drafts of the text. CL and PM conceptualized the methodologies and revised the text, provided feedback, and revised the text. RM helped to organize the paper, prepared the figures, and addressed reviewers comments. All authors contributed to the article and approved the submitted version.
Funding
The authors are grateful for funding from Louisiana State University's Economic Development Assistantship and Dissertation Year Fellowship programs. The publication of this article is subsidized by the LSU Libraries Open Access Author Fund.
Acknowledgments
The authors are grateful for funding from the Louisiana State University Economic Development Assistantship program for the completion of this project and NB doctoral degree program.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
Methods associated with the generation of the data used in this research to characterize the NHCPV properties are described in more detail in Bushra and Rohli (2019). The gridded data set was obtained from the National Centers for Environmental Prediction (NCEP)/Department of Energy (DOE) Reanalysis 2 project. The data are available at https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis2.pressure.html.
References
Ahdesmaki, M., Fokianos, K., and Strimmer, K. (2021). GeneCycle: Identification of Periodically Expressed Genes. (2012). R package version 1.1.5. Available online at: https://CRAN.R-project.org/package=GeneCycle (accessed May 03, 2021).
Altenhoff, A. M., Martius, O., Croci-Maspoli, M., Schwierz, C., and Davies, H. C. (2008). Linkage of atmospheric blocks and synoptic-scale Rossby waves: a climatological analysis. Tellus Dyn. Meteorol. Oceanogr. 60, 1053–1063. doi: 10.1111/j.1600-0870.2008.00354.x
Angell, J. K. (1998). Contraction of the 300 mbar north circumpolar vortex during 1963–1997 and its movement into the eastern hemisphere. J. Geophys. Res. Atmosph. 103, 25887–25893. doi: 10.1029/98JD02451
Angell, J. K. (2006). Changes in the 300-mb north circumpolar vortex, 1963–2001. J. Clim. 19, 2984–2994. doi: 10.1175/JCLI3778.1
Angell, J. K., and Korshover, J. (1977). Variation in size and location of 300 mb north circumpolar vortex between 1963 and 1975. Mon. Weather Rev. 105, 19–25.
Angell, J. K., and Korshover, J. (1985). Displacements of the north circumpolar vortex during El Niño, 1963–83. Mon. Weather Rev. 113, 1627–1630.
Bai, J. (1994). Least squares estimation of a shift in linear processes. J. Time Ser. Anal. 15, 453–472. doi: 10.1111/j.1467-9892.1994.tb00204.x
Bai, J. (1997a). Estimation of a change point in multiple regression models. Rev. Econ. Stat. 79, 551–563. doi: 10.1162/003465397557132
Bai, J. (1997b). Estimating multiple breaks one at a time. Econ. Theory 13, 315–352. doi: 10.1017/S0266466600005831
Bai, J., and Perron, P. (2003). Computation and analysis of multiple structural change models. J. Appl. Econ. 18, 1–22. doi: 10.1002/jae.659
Bai, J. S., and Perron, P. (1998). Estimating and testing linear models with multiple structural changes. Econometrica 66, 47–78. doi: 10.2307/2998540
Barnston, A. G., and Livezey, R. E. (1987). Classification, seasonality and persistence of low-frequency atmospheric circulation patterns. Mon. Weather Rev. 115, 1083–1126. doi: 10.1175/1520-0493(1987)115<1083:CSAPOL>2.0.CO;2
Bartlett, R. E., Bollasina, M. A., Booth, B. B., Dunstone, N. J., Marenco, F., Messori, G., et al. (2018). Do differences in future sulfate emission pathways matter for near-term climate? A case study for the Asian monsoon. Climate Dyn. 50, 1863–1880. doi: 10.1007/s00382-017-3726-6
Burnett, A. W. (1993). Size variations and long-wave circulation within the January northern-hemisphere circumpolar vortex-−1946–89. J. Clim. 6, 1914–1920.
Bushra, N., and Rohli, R. V. (2019). An objective procedure for delineating the circumpolar vortex. Earth Sp. Sci. 6, 774–783. doi: 10.1029/2019EA000590
Bushra, N., and Rohli, R. V. (2021a). Relationship between atmospheric teleconnections and the Northern Hemisphere's circumpolar vortex. Earth Sp. Sci. 8, e2021EA001802. doi: 10.1029/2021EA001802
Bushra, N., and Rohli, R. V. (2021b). Spatiotemporal trends and variability in the centroid of the Northern Hemisphere's circumpolar vortex. Earth Sp. Sci. 8, e2020EA001594. doi: 10.1029/2020EA001594
Cattiaux, J., Peings, Y., Saint-Martin, D., Trou-Kechout, N., and Vavrus, S. J. (2016). Sinuosity of midlatitude atmospheric flow in a warming world. Geophys. Res. Lett. 43, 8259–8268. doi: 10.1002/2016GL070309
Chen, G., Peng, L., and Ma, C. (2018). Climatology and seasonality of upper ocean salinity: a three-dimensional view from argo floats. Climate Dyn. 50, 2169–2182. doi: 10.1007/s00382-017-3742-6
Chen, T. C., Marshall, H. G., and Shukla, J. (1981). Spectral analysis and diagnosis of nonlinear interactions of large-scale moving waves at 200 mb in the GLAS general circulation model. Mon. Weather Rev. 109, 959–974.
Chorley, R. J., Schumm, S. A., and Sugden, D. E. (1984). Geomorphology. London, New York: Metheun and Co.
Cohen, J., Screen, J. A., Furtado, J. C., Barlow, M., Whittleston, D., Coumou, D., et al. (2014). Recent Arctic amplification and extreme mid-latitude weather. Nat. Geosci. 7, 627–637. doi: 10.1038/ngeo2234
Davis, R. E., and Benkovic, S. R. (1992). Climatological variations in the Northern-Hemisphere circumpolar vortex in January. Theor. Appl. Climatol. 46, 63–73. doi: 10.1007/BF00866086
Davis, R. E., and Benkovic, S. R. (1994). Spatial and temporal variations of the January circumpolar vortex over the Northern-Hemisphere. Int. J. Climatol. 14, 415–428. doi: 10.1002/joc.3370140406
Di Capua, G., and Coumou, D. (2016). Changes in meandering of the Northern Hemisphere circulation. Environ. Res. Lett. 11, 094028. doi: 10.1088/1748-9326/11/9/094028
Folland, C. K., Boucher, O., Colman, A., and Parker, D. E. (2018). Causes of irregularities in trends of global mean surface temperature since the late 19th century. Sci. Adv. 4, eaao5297. doi: 10.1126/sciadv.aao5297
Francis, J. A., and Vavrus, S. J. (2012). Evidence linking Arctic amplification to extreme weather in mid-latitudes. Geophys. Res. Lett. 39, L06801. doi: 10.1029/2012GL051000
Francis, J. A., and Vavrus, S. J. (2015). Evidence for a wavier jet stream in response to rapid Arctic warming. Environ. Res. Lett. 10, 014005. doi: 10.1088/1748-9326/10/1/014005
Frauenfeld, O. W., and Davis, R. E. (2000). The influence of El Nino-Southern Oscillation events on the Northern Hemisphere 500 hPa circumpolar vortex. Geophys. Res. Lett., 27, 537–540. doi: 10.1029/1999GL010996
Frauenfeld, O. W., and Davis, R. E. (2003). Northern Hemisphere circumpolar vortex trends and climate change implications. J. Geophys. Res. Atmosph. 108, 4423. doi: 10.1029/2002JD002958
Frauenfeld, O. W., Davis, R. E., and Mann, M. E. (2005). A distinctly interdecadal signal of Pacific Ocean–atmosphere interaction. J. Clim., 18, 1709–1718. doi: 10.1175/JCLI3367.1
Ghil, M., Allen, M. R., Dettinger, M. D., Ide, K., Kondrashov, D., Mann, M. E., et al. (2002). Advanced spectral methods for climatic time series. Rev. Geophys. 40, 3–1. doi: 10.1029/2000RG000092
Glovin, G. M., Arbetter, T. E., and Lynch, A. H. (2016). Wavelet analysis of polar vortex variability over the twentieth century. J. Geophys. Res. Atmosph. 121, 722–732. doi: 10.1002/2014JD022971
Gollan, G., and Greatbatch, R. J. (2017). The relationship between Northern Hemisphere winter blocking and tropical modes of variability. J. Clim. 30, 9321–9337. doi: 10.1175/JCLI-D-16-0742.1
Hayashi, Y. (1974). Spectral analysis of tropical disturbances appearing in a GFDL general circulation model. J. Atmosph. Sci. 31, 180–218.
Hayashi, Y. (1982). Space-time spectral analysis and its applications to atmospheric waves. J. Meteorol. Soc. Jpn. Ser. II 60, 156–171. doi: 10.2151/jmsj1965.60.1_156
Henderson, S. A., Maloney, E. D., and Barnes, E. A. (2016). The influence of the Madden–Julian oscillation on Northern Hemisphere winter blocking. J. Clim., 29, 4597–4616. doi: 10.1175/JCLI-D-15-0502.1
Higuchi, K., Huang, J., and Shabbar, A. (1999). A wavelet characterization of the North Atlantic oscillation variation and its relationship to the North Atlantic sea surface temperature. Int. J. Climatol. 19, 1119–1129.
Holland, M. M., and Bitz, C. M. (2003). Polar amplification of climate change in coupled models. Climate Dyn. 21, 221–232. doi: 10.1007/s00382-003-0332-6
Hurrell, J. W. (1995). Decadal trends in the North Atlantic Oscillation: Regional temperatures and precipitation. Science 269, 676–679. doi: 10.1126/science.269.5224.676
Jekeli, C. (2005). Spline Representations of Functions on a Sphere for Geopotential Modeling. Report No. 475. Columbus: Department of Civil and Environmental Engineering and Geodetic Science, The Ohio State University. p. 37.
Kanamitsu, M., Ebisuzaki, W., Woollen, J., Yang, S-K, Hnilo, J. J., Fiorino, M., and Potter, G.L. (2002). NCEP-DOE AMIP-II Reanalysis (R-2). Bull. Am. Meteorol. Soc. 83, 1631–1643. doi: 10.1175/BAMS-83-11-1631
Kerr, R. A. (2000). A North Atlantic climate pacemaker for the centuries. Science 288, 1984–1986 doi: 10.1126/science.288.5473.1984
Kidston, J., Scaife, A. A., Hardiman, S. C., Mitchell, D. M., Butchart, N., Baldwin, M. P., et al. (2015). Stratospheric influence on tropospheric jet streams, storm tracks and surface weather. Nat. Geosci. 8, 433–440. doi: 10.1038/ngeo2424
Kirby, M., Patterson, W., Mullins, H., and Burnett, A. (2002). Post-Younger Dryas climate interval linked to circumpolar vortex variability: Isotopic evidence from Fayetteville Green Lake, New York. Climate Dyn. 19, 321–330. doi: 10.1007/s00382-002-0227-y
Lau, W. K. M., Waliser, D. E., and Roundy, P. E. (2012). “Tropical–extratropical interactions,” in Intraseasonal Variability in the Atmosphere-Ocean Climate System (Springer Praxis Books. Springer, Berlin, Heidelberg). doi: 10.1007/978-3-642-13914-7_14
Lupo, A. R. (2021). Atmospheric blocking events: A review. Ann. N. Y. Acad. Sci. 1504, 5–24. doi: 10.1111/nyas.14557
Madden, R. A., and Julian, P. R. (1971). Detection of a 40–50 day oscillation in the zonal wind in the tropical Pacific. J. Atmospheric Sci. 28, 702–708.
Madden, R. A., and Julian, P. R. (1994). Observations of the 40–50-day tropical oscillation—a review. Mon. Weather Rev. 122, 814–837.
Martin, J. E. (2015). Contraction of the Northern Hemisphere, lower-tropospheric wintertime cold pool over the last 66 years. J. Clim. 28, 3764–3778. doi: 10.1175/JCLI-D-14-00496.1
Martin, J. E. (2021). Recent trends in the waviness of the Northern Hemisphere wintertime polar and subtropical jets. J. Geophys. Res. Atmosph. 126, e2020JD033668. doi: 10.1029/2020JD033668
Montgomery, D. C., Peck, E. A., and Vining, G. G. (2012). Introduction to Linear Regression Analysis, 5th Edition. John Wiley and Sons, New York.
Moron, V., Oueslati, B., Pohl, B., and Janicot, S. (2018). Daily weather types in February–June (1979–2016) and temperature variations in tropical North Africa. J. Appl. Meteorol. Climatol. 57, 1171–1195. doi: 10.1175/JAMC-D-17-0105.1
Muggeo, V. M. (2008). Segmented: an R package to fit regression models with broken-line relationships. R News 8, 20–25.
National Oceanic Atmospheric Administration (2021). Available online at: https://psl.noaa.gov/pdo/ (accessed July 2, 2022).
Newman, M., Compo, G. P., and Alexander, M. A. (2003). ENSO-forced variability of the Pacific decadal oscillation. J. Clim. 16, 3853–3857. doi: 10.1175/1520-0442(2003)016<3853:EVOTPD>2.0.CO;2
Orme, L. C., Charman, D. J., Reinhardt, L., Jones, R. T., Mitchell, F. J., Stefanini, B. S., et al. (2017). Past changes in the North Atlantic storm track driven by insolation and sea-ice forcing. Geology 45, 335–338. doi: 10.1130/G38521.1
Otterå, O. H. (2008). Simulating the effects of the 1991 Mount Pinatubo volcanic eruption using the ARPEGE atmosphere general circulation model. Adv. Atmosph. Sci. 25, 213–226. doi: 10.1007/s00376-008-0213-3
Overland, J., Francis, J. A., Hall, R., Hanna, E., Kim, S. J., and Vihma, T. (2015). The melting Arctic and midlatitude weather patterns: are they connected? J. Clim., 28, 7917–7932. doi: 10.1175/JCLI-D-14-00822.1
Peings, Y., and Magnusdottir, G. (2014). Forcing of the wintertime atmospheric circulation by the multidecadal fluctuations of the North Atlantic Ocean. Environ. Res. Lett. 9, 034018. doi: 10.1088/1748-9326/9/3/034018
Pelletier, J. D. (1998). The power spectral density of atmospheric temperature from time scales of 10–2 to 106 yr. Earth Planet. Sci. Lett. 158, 157–164. doi: 10.1016/S0012-821X(98)00051-X
Philander, S.G. (1990). El Niño, La Niña, and the Southern Oscillation, 293 pp. San Diego: Academic Press.
Pithan, F., and Mauritsen, T. (2014). Arctic amplification dominated by temperature feedbacks in contemporary climate models. Nat. Geosci. 7, 181–184, doi: 10.1038/ngeo2071
Post, E., Alley, R. B., Christensen, T. R., Macias-Fauria, M., Forbes, B. C., Gooseff, M. N., et al. (2019). The polar regions in a 2°C warmer world. Sci. Adv. 5, eaaw9883. doi: 10.1126/sciadv.aaw9883
Praetorius, S. K. (2018). North Atlantic circulation slows down. Nature 556, 180–181. doi: 10.1038/d41586-018-04086-4
Rahmstorf, S., Box, J. E., Feulner, G., Mann, M. E., Robinson, A., Rutherford, S., et al. (2015). Exceptional twentieth-century slowdown in Atlantic Ocean overturning circulation. Nat. Clim. Chang. 5, 475–480. doi: 10.1038/nclimate2554
Rodionov, S. N. (2004). A sequential algorithm for testing climate regime shifts. Geophys. Res. Lett. 31, L09204, doi: 10.1029/2004GL019448
Rohli, R. V., Snedden, G., Martin, E. R., and DeLong, K. L. (2022). Impacts of atmosphere-ocean teleconnections on the south-central United States. Front. Earth Sci. 10, 934654, doi: 10.3389/feart.2022.934654
Rohli, R. V., Wrona, K. M., and McHugh, M. J. (2005). January northern hemisphere circumpolar vortex variability and its relationship with hemispheric temperature and regional teleconnections. Int. J. Climatol. 25, 1421–1436. doi: 10.1002/joc.1204
Röthlisberger, M., Pfahl, S., and Martius, O. (2016). Regional-scale jet waviness modulates the occurrence of midlatitude weather extremes. Geophys. Res. Lett. 43, 10989–10997. doi: 10.1002/2016GL070944
Savin, N. E., and White, K. J. (1977). The Durbin-Watson test for serial correlation with extreme sample sizes or many regressors. Econ. J. Econ. Soc. 45, 1989–1996. doi: 10.2307/1914122
Schulz, M., and Mudelsee, M. (2002). REDFIT: Estimating red-noise spectra directly from unevenly spaced paleoclimatic time series. Comput. Geosci. 28, 421–426. doi: 10.1016/S0098-3004(01)00044-9
Screen, J. A., and Simmonds, I. (2013). Exploring links between Arctic amplification and mid-latitude weather. Geophys. Res. Lett. 40, 959–964. doi: 10.1002/grl.50174
Srinivas, G., Chowdary, J. S., Kosaka, Y., Gnanaseelan, C., Parekh, A., and Prasad, K. V. S. R. (2018). Influence of the Pacific–Japan Pattern on Indian summer monsoon rainfall. J. Clim. 31, 3943–3958. doi: 10.1175/JCLI-D-17-0408.1
Sun, C., and Li, J. (2012). Space–time spectral analysis of the Southern Hemisphere daily 500-hPa geopotential height. Mon. Weather Rev. 140, 3844–3856. doi: 10.1175/MWR-D-12-00019.1
Tang, C. K., and Chan, J. C. L. (2016). Idelized simulations of the effect of Taiwan topography on the tracks of tropical cyclones with different steering flow strengths. Q. J. R. Meteorol. Soc. 142, 3211–3221. doi: 10.1002/qj.2902
Thompson, D. W., and Wallace, J. M. (1998). The Arctic Oscillation signature in the wintertime geopotential height and temperature fields. Geophys. Res. Lett. 25, 1297–1300. doi: 10.1029/98GL00950
Thompson, D. W., and Wallace, J. M. (2000). Annular modes in the extratropical circulation. Part I: Month-to-month variability. J. Clim. 13, 1000–1016. doi: 10.1175/1520-0442(2000)013<1000:AMITEC>2.0.CO;2
Thompson, D. W. J., and Solomon, S. (2002). Interpretation of recent Southern Hemisphere climate change. Science 296, 895–899. doi: 10.1126/science.1069270
Torrence, C., and Compo, G. P. (1998). A ractical guide to wavelet analysis. Bull. Am. Meteorol. Soc. 79, 61–78.
Tyrlis, E., and Hoskins, B. J. (2008). Aspects of a Northern Hemisphere atmospheric blocking climatology. J. Atmosph. Sci. 65, 1638–1652. doi: 10.1175/2007JAS2337.1
van den Broeke, M. R., and van Lipzig, N. P. M. (2002). Impact of polar vortex variability on the wintertime low-level climate of East Antarctica: results of a regional climate model. Tellus Dynamic Meteorol. Oceanogr. 54, 485–496. doi: 10.3402/tellusa.v54i5.12161
van Loon, H., and Rogers, J. C. (1978). The seesaw in winter temperatures between Greenland and northern Europe. Part I: General description. Mon. Weather Rev. 106, 296–310.
Vanhatalo, E., and Kulahci, M. (2016). Impact of autocorrelation on principal components and their use in statistical process control. Quality Reliab. Eng. Int. 32, 1483–1500. doi: 10.1002/qre.1858
Vanos, J. K., and Cakmak, S. (2014). Changing air mass frequencies in Canada: Potential links and implications for human health. Int. J. Biometeorol. 58, 121–135. doi: 10.1007/s00484-013-0634-2
Vavrus, S. J. (2018). The influence of Arctic amplification on mid-latitude weather and climate. Current Climate Change Rep. 4, 238–249. doi: 10.1007/s40641-018-0105-2
Vavrus, S. J., Wang, F., Martin, J. E., Francis, J. A., Peings, Y., and Cattiaux, J. (2017). Changes in North American atmospheric circulation and extreme weather: Influence of Arctic amplification and Northern Hemisphere snow cover. J. Clim. 30, 4317–4333. doi: 10.1175/JCLI-D-16-0762.1
Wallace, J. M., and Gutzler, D. S. (1981). Teleconnections in the geopotential height field during the Northern Hemisphere winter. Mon. Weather Rev. 109, 784–812.
Wang, Z. Y., and Ding, Y. H. (2009). Impacts of the long-term change of the summer Asian polar vortex on the circulation system and the water vapor transport in East Asia. Chin. J. Geophys. 52, 20–29.
Waugh, D. W., Sobel, A. H., and Polvani, L. M. (2017). What is the polar vortex and how does it influence weather? Bull. Am. Meteorol. Soc. 98, 37–44. doi: 10.1175/BAMS-D-15-00212.1
Welch, P. (1967). The use of fast Fourier transform for the estimation of power spectra: a method based on time averaging over short, modified periodograms. IEEE Trans. Audio Electroacoust. 15, 70–73. doi: 10.1109/TAU.1967.1161901
Wrona, K. M., and Rohli, R. V. (2007). Seasonality of the northern hemisphere circumpolar vortex. Int. J. Climatol. 27, 697–713. doi: 10.1002/joc.1430
Keywords: tropospheric circumpolar vortex, 500-mb flow patterns, steering atmospheric circulation, fast Fourier transform, breakpoint analysis, climate change
Citation: Bushra N, Rohli RV, Li C, Miller PW and Mostafiz RB (2023) Changing features of the Northern Hemisphere 500-hPa circumpolar vortex. Front. Big Data 5:1009158. doi: 10.3389/fdata.2022.1009158
Received: 01 August 2022; Accepted: 07 December 2022;
Published: 09 January 2023.
Edited by:
Letizia Lusito, Inside Climate Service, ItalyReviewed by:
Antonio Speranza, University of Camerino, ItalyAnthony Lupo, University of Missouri, United States
Copyright © 2023 Bushra, Rohli, Li, Miller and Mostafiz. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nazla Bushra, bmJ1c2hyMSYjeDAwMDQwO2xzdS5lZHU=