ORIGINAL RESEARCH article

Front. Big Data, 09 January 2023

Sec. Data-driven Climate Sciences

Volume 5 - 2022 | https://doi.org/10.3389/fdata.2022.1009158

Changing features of the Northern Hemisphere 500-hPa circumpolar vortex

  • 1. Department of Oceanography and Coastal Sciences, College of the Coast and Environment, Louisiana State University, Baton Rouge, LA, United States

  • 2. Coastal Studies Institute, Louisiana State University, Baton Rouge, LA, United States

Article metrics

View details

2

Citations

3,1k

Views

1k

Downloads

Abstract

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

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

AreaRc
Sub-periodsRegression line (p-value)Breakpoints (confidence interval)Sub-periodsRegression line (p-value)Breakpoints (confidence interval)
01/01/1979–5/26/1986y = 4.662e+13 + 5.871e+08x (p = 0.086)01/01/1979–6/4/1988y = 8.577e-01 + 3.135e-05x (p = 0.061)
5/27/1986–10/25/1994y = 4.338e+13 + 3.769e+08x (p = 0.0413)5/27/1986 (6/2/1985 −5/19/1987)6/5/1988– 12/1/1994y = 8.720e-01 – 1.699e-06x (p = 0.031)6/5/1988 (11/22/1987–1/13/1989)
10/26/1994–3/24/2001y = 4.075e+13 + 1.305e+09x (p = 0.002)10/26/1994 (7/30/1994 −1/25/1995)12/2/1994–5/28/2001y = 8.429e-01 + 1.925e-06x (p = 0.017)12/2/1994 (3/26/1994–7/17/1995)
3/25/2001–12/31/2017y = 4.541e+13 + 1.884e+08x (p = 0.040)3/25/2001 (2/2/2000 −4/13/2002)5/28/2001– 12/31/2017y = 8.621e-01 – 1.307e-05x (p = 0.068)5/28/20013/12/1997–2/3/2003)

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

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

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.

Statements

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

  • 1

    AhdesmakiM.FokianosK.StrimmerK. (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).

  • 2

    AltenhoffA. M.MartiusO.Croci-MaspoliM.SchwierzC.DaviesH. C. (2008). Linkage of atmospheric blocks and synoptic-scale Rossby waves: a climatological analysis. Tellus Dyn. Meteorol. Oceanogr.60, 10531063. 10.1111/j.1600-0870.2008.00354.x

  • 3

    AngellJ. 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, 2588725893. 10.1029/98JD02451

  • 4

    AngellJ. K. (2006). Changes in the 300-mb north circumpolar vortex, 1963–2001. J. Clim.19, 29842994. 10.1175/JCLI3778.1

  • 5

    AngellJ. K.KorshoverJ. (1977). Variation in size and location of 300 mb north circumpolar vortex between 1963 and 1975. Mon. Weather Rev.105, 1925.

  • 6

    AngellJ. K.KorshoverJ. (1985). Displacements of the north circumpolar vortex during El Niño, 1963–83. Mon. Weather Rev.113, 16271630.

  • 7

    BaiJ. (1994). Least squares estimation of a shift in linear processes. J. Time Ser. Anal.15, 453472. 10.1111/j.1467-9892.1994.tb00204.x

  • 8

    BaiJ. (1997a). Estimation of a change point in multiple regression models. Rev. Econ. Stat.79, 551563. 10.1162/003465397557132

  • 9

    BaiJ. (1997b). Estimating multiple breaks one at a time. Econ. Theory13, 315352. 10.1017/S0266466600005831

  • 10

    BaiJ.PerronP. (2003). Computation and analysis of multiple structural change models. J. Appl. Econ.18, 122. 10.1002/jae.659

  • 11

    BaiJ. S.PerronP. (1998). Estimating and testing linear models with multiple structural changes. Econometrica66, 4778. 10.2307/2998540

  • 12

    BarnstonA. G.LivezeyR. E. (1987). Classification, seasonality and persistence of low-frequency atmospheric circulation patterns. Mon. Weather Rev.115, 10831126. 10.1175/1520-0493(1987)115<1083:CSAPOL>2.0.CO;2

  • 13

    BartlettR. E.BollasinaM. A.BoothB. B.DunstoneN. J.MarencoF.MessoriG.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, 18631880. 10.1007/s00382-017-3726-6

  • 14

    BurnettA. W. (1993). Size variations and long-wave circulation within the January northern-hemisphere circumpolar vortex-−1946–89. J. Clim.6, 19141920.

  • 15

    BushraN.RohliR. V. (2019). An objective procedure for delineating the circumpolar vortex. Earth Sp. Sci.6, 774783. 10.1029/2019EA000590

  • 16

    BushraN.RohliR. V. (2021a). Relationship between atmospheric teleconnections and the Northern Hemisphere's circumpolar vortex. Earth Sp. Sci.8, e2021EA001802. 10.1029/2021EA001802

  • 17

    BushraN.RohliR. V. (2021b). Spatiotemporal trends and variability in the centroid of the Northern Hemisphere's circumpolar vortex. Earth Sp. Sci.8, e2020EA001594. 10.1029/2020EA001594

  • 18

    ButterworthS. (1930). On the theory of filter amplifiers. Wireless Eng.7, 536541.

  • 19

    CattiauxJ.PeingsY.Saint-MartinD.Trou-KechoutN.VavrusS. J. (2016). Sinuosity of midlatitude atmospheric flow in a warming world. Geophys. Res. Lett.43, 82598268. 10.1002/2016GL070309

  • 20

    ChenG.PengL.MaC. (2018). Climatology and seasonality of upper ocean salinity: a three-dimensional view from argo floats. Climate Dyn.50, 21692182. 10.1007/s00382-017-3742-6

  • 21

    ChenT. C.MarshallH. G.ShuklaJ. (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, 959974.

  • 22

    ChorleyR. J.SchummS. A.SugdenD. E. (1984). Geomorphology. London, New York: Metheun and Co.

  • 23

    CohenJ.ScreenJ. A.FurtadoJ. C.BarlowM.WhittlestonD.CoumouD.et al. (2014). Recent Arctic amplification and extreme mid-latitude weather. Nat. Geosci.7, 627637. 10.1038/ngeo2234

  • 24

    DavisR. E.BenkovicS. R. (1992). Climatological variations in the Northern-Hemisphere circumpolar vortex in January. Theor. Appl. Climatol.46, 6373. 10.1007/BF00866086

  • 25

    DavisR. E.BenkovicS. R. (1994). Spatial and temporal variations of the January circumpolar vortex over the Northern-Hemisphere. Int. J. Climatol.14, 415428. 10.1002/joc.3370140406

  • 26

    Di CapuaG.CoumouD. (2016). Changes in meandering of the Northern Hemisphere circulation. Environ. Res. Lett.11, 094028. 10.1088/1748-9326/11/9/094028

  • 27

    FollandC. K.BoucherO.ColmanA.ParkerD. E. (2018). Causes of irregularities in trends of global mean surface temperature since the late 19th century. Sci. Adv. 4, eaao5297. 10.1126/sciadv.aao5297

  • 28

    FrancisJ. A.VavrusS. J. (2012). Evidence linking Arctic amplification to extreme weather in mid-latitudes. Geophys. Res. Lett.39, L06801. 10.1029/2012GL051000

  • 29

    FrancisJ. A.VavrusS. J. (2015). Evidence for a wavier jet stream in response to rapid Arctic warming. Environ. Res. Lett.10, 014005. 10.1088/1748-9326/10/1/014005

  • 30

    FrauenfeldO. W.DavisR. E. (2000). The influence of El Nino-Southern Oscillation events on the Northern Hemisphere 500 hPa circumpolar vortex. Geophys. Res. Lett., 27, 537540. 10.1029/1999GL010996

  • 31

    FrauenfeldO. W.DavisR. E. (2003). Northern Hemisphere circumpolar vortex trends and climate change implications. J. Geophys. Res. Atmosph.108, 4423. 10.1029/2002JD002958

  • 32

    FrauenfeldO. W.DavisR. E.MannM. E. (2005). A distinctly interdecadal signal of Pacific Ocean–atmosphere interaction. J. Clim., 18, 17091718. 10.1175/JCLI3367.1

  • 33

    GhilM.AllenM. R.DettingerM. D.IdeK.KondrashovD.MannM. E.et al. (2002). Advanced spectral methods for climatic time series. Rev. Geophys.40, 31. 10.1029/2000RG000092

  • 34

    GlovinG. M.ArbetterT. E.LynchA. H. (2016). Wavelet analysis of polar vortex variability over the twentieth century. J. Geophys. Res. Atmosph.121, 722732. 10.1002/2014JD022971

  • 35

    GollanG.GreatbatchR. J. (2017). The relationship between Northern Hemisphere winter blocking and tropical modes of variability. J. Clim.30, 93219337. 10.1175/JCLI-D-16-0742.1

  • 36

    HayashiY. (1974). Spectral analysis of tropical disturbances appearing in a GFDL general circulation model. J. Atmosph. Sci.31, 180218.

  • 37

    HayashiY. (1982). Space-time spectral analysis and its applications to atmospheric waves. J. Meteorol. Soc. Jpn. Ser. II60, 156171. 10.2151/jmsj1965.60.1_156

  • 38

    HendersonS. A.MaloneyE. D.BarnesE. A. (2016). The influence of the Madden–Julian oscillation on Northern Hemisphere winter blocking. J. Clim., 29, 45974616. 10.1175/JCLI-D-15-0502.1

  • 39

    HiguchiK.HuangJ.ShabbarA. (1999). A wavelet characterization of the North Atlantic oscillation variation and its relationship to the North Atlantic sea surface temperature. Int. J. Climatol.19, 11191129.

  • 40

    HollandM. M.BitzC. M. (2003). Polar amplification of climate change in coupled models. Climate Dyn.21, 221232. 10.1007/s00382-003-0332-6

  • 41

    HurrellJ. W. (1995). Decadal trends in the North Atlantic Oscillation: Regional temperatures and precipitation. Science269, 676679. 10.1126/science.269.5224.676

  • 42

    JekeliC. (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.

  • 43

    KanamitsuM.EbisuzakiW.WoollenJ.YangS-KHniloJ. J.FiorinoM.PotterG.L. (2002). NCEP-DOE AMIP-II Reanalysis (R-2). Bull. Am. Meteorol. Soc.83, 16311643. 10.1175/BAMS-83-11-1631

  • 44

    KerrR. A. (2000). A North Atlantic climate pacemaker for the centuries. Science288, 1984198610.1126/science.288.5473.1984

  • 45

    KidstonJ.ScaifeA. A.HardimanS. C.MitchellD. M.ButchartN.BaldwinM. P.et al. (2015). Stratospheric influence on tropospheric jet streams, storm tracks and surface weather. Nat. Geosci.8, 433440. 10.1038/ngeo2424

  • 46

    KirbyM.PattersonW.MullinsH.BurnettA. (2002). Post-Younger Dryas climate interval linked to circumpolar vortex variability: Isotopic evidence from Fayetteville Green Lake, New York. Climate Dyn.19, 321330. 10.1007/s00382-002-0227-y

  • 47

    LauW. K. M.WaliserD. E.RoundyP. E. (2012). “Tropical–extratropical interactions,” in Intraseasonal Variability in the Atmosphere-Ocean Climate System (Springer Praxis Books. Springer, Berlin, Heidelberg). 10.1007/978-3-642-13914-7_14

  • 48

    LupoA. R. (2021). Atmospheric blocking events: A review. Ann. N. Y. Acad. Sci.1504, 524. 10.1111/nyas.14557

  • 49

    MaddenR. A.JulianP. R. (1971). Detection of a 40–50 day oscillation in the zonal wind in the tropical Pacific. J. Atmospheric Sci.28, 702708.

  • 50

    MaddenR. A.JulianP. R. (1994). Observations of the 40–50-day tropical oscillation—a review. Mon. Weather Rev.122, 814837.

  • 51

    MartinJ. E. (2015). Contraction of the Northern Hemisphere, lower-tropospheric wintertime cold pool over the last 66 years. J. Clim.28, 37643778. 10.1175/JCLI-D-14-00496.1

  • 52

    MartinJ. E. (2021). Recent trends in the waviness of the Northern Hemisphere wintertime polar and subtropical jets. J. Geophys. Res. Atmosph.126, e2020JD033668. 10.1029/2020JD033668

  • 53

    MontgomeryD. C.PeckE. A.ViningG. G. (2012). Introduction to Linear Regression Analysis, 5th Edition. John Wiley and Sons, New York.

  • 54

    MoronV.OueslatiB.PohlB.JanicotS. (2018). Daily weather types in February–June (1979–2016) and temperature variations in tropical North Africa. J. Appl. Meteorol. Climatol.57, 11711195. 10.1175/JAMC-D-17-0105.1

  • 55

    MuggeoV. M. (2008). Segmented: an R package to fit regression models with broken-line relationships. R News8, 2025.

  • 56

    National Oceanic Atmospheric Administration (2021). Available online at: https://psl.noaa.gov/pdo/ (accessed July 2, 2022).

  • 57

    NewmanM.CompoG. P.AlexanderM. A. (2003). ENSO-forced variability of the Pacific decadal oscillation. J. Clim.16, 38533857. 10.1175/1520-0442(2003)016<3853:EVOTPD>2.0.CO;2

  • 58

    OrmeL. C.CharmanD. J.ReinhardtL.JonesR. T.MitchellF. J.StefaniniB. S.et al. (2017). Past changes in the North Atlantic storm track driven by insolation and sea-ice forcing. Geology45, 335338. 10.1130/G38521.1

  • 59

    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, 213226. 10.1007/s00376-008-0213-3

  • 60

    OverlandJ.FrancisJ. A.HallR.HannaE.KimS. J.VihmaT. (2015). The melting Arctic and midlatitude weather patterns: are they connected?J. Clim., 28, 79177932. 10.1175/JCLI-D-14-00822.1

  • 61

    PeingsY.MagnusdottirG. (2014). Forcing of the wintertime atmospheric circulation by the multidecadal fluctuations of the North Atlantic Ocean. Environ. Res. Lett.9, 034018. 10.1088/1748-9326/9/3/034018

  • 62

    PelletierJ. D. (1998). The power spectral density of atmospheric temperature from time scales of 10–2 to 106 yr. Earth Planet. Sci. Lett.158, 157164. 10.1016/S0012-821X(98)00051-X

  • 63

    PhilanderS.G. (1990). El Niño, La Niña, and the Southern Oscillation, 293 pp. San Diego: Academic Press.

  • 64

    PithanF.MauritsenT. (2014). Arctic amplification dominated by temperature feedbacks in contemporary climate models. Nat. Geosci.7, 181184, 10.1038/ngeo2071

  • 65

    PostE.AlleyR. B.ChristensenT. R.Macias-FauriaM.ForbesB. C.GooseffM. N.et al. (2019). The polar regions in a 2°C warmer world. Sci. Adv. 5, eaaw9883. 10.1126/sciadv.aaw9883

  • 66

    PraetoriusS. K. (2018). North Atlantic circulation slows down. Nature556, 180181. 10.1038/d41586-018-04086-4

  • 67

    RahmstorfS.BoxJ. E.FeulnerG.MannM. E.RobinsonA.RutherfordS.et al. (2015). Exceptional twentieth-century slowdown in Atlantic Ocean overturning circulation. Nat. Clim. Chang.5, 475480. 10.1038/nclimate2554

  • 68

    RodionovS. N. (2004). A sequential algorithm for testing climate regime shifts. Geophys. Res. Lett.31, L09204, 10.1029/2004GL019448

  • 69

    RohliR. V.SneddenG.MartinE. R.DeLongK. L. (2022). Impacts of atmosphere-ocean teleconnections on the south-central United States. Front. Earth Sci.10, 934654, 10.3389/feart.2022.934654

  • 70

    RohliR. V.WronaK. M.McHughM. J. (2005). January northern hemisphere circumpolar vortex variability and its relationship with hemispheric temperature and regional teleconnections. Int. J. Climatol.25, 14211436. 10.1002/joc.1204

  • 71

    RöthlisbergerM.PfahlS.MartiusO. (2016). Regional-scale jet waviness modulates the occurrence of midlatitude weather extremes. Geophys. Res. Lett.43, 1098910997. 10.1002/2016GL070944

  • 72

    SavinN. E.WhiteK. J. (1977). The Durbin-Watson test for serial correlation with extreme sample sizes or many regressors. Econ. J. Econ. Soc.45, 19891996. 10.2307/1914122

  • 73

    SchulzM.MudelseeM. (2002). REDFIT: Estimating red-noise spectra directly from unevenly spaced paleoclimatic time series. Comput. Geosci.28, 421426. 10.1016/S0098-3004(01)00044-9

  • 74

    ScreenJ. A.SimmondsI. (2013). Exploring links between Arctic amplification and mid-latitude weather. Geophys. Res. Lett.40, 959964. 10.1002/grl.50174

  • 75

    SrinivasG.ChowdaryJ. S.KosakaY.GnanaseelanC.ParekhA.PrasadK. V. S. R. (2018). Influence of the Pacific–Japan Pattern on Indian summer monsoon rainfall. J. Clim.31, 39433958. 10.1175/JCLI-D-17-0408.1

  • 76

    SunC.LiJ. (2012). Space–time spectral analysis of the Southern Hemisphere daily 500-hPa geopotential height. Mon. Weather Rev.140, 38443856. 10.1175/MWR-D-12-00019.1

  • 77

    TangC. K.ChanJ. 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, 32113221. 10.1002/qj.2902

  • 78

    ThompsonD. W.WallaceJ. M. (1998). The Arctic Oscillation signature in the wintertime geopotential height and temperature fields. Geophys. Res. Lett.25, 12971300. 10.1029/98GL00950

  • 79

    ThompsonD. W.WallaceJ. M. (2000). Annular modes in the extratropical circulation. Part I: Month-to-month variability. J. Clim.13, 10001016. 10.1175/1520-0442(2000)013<1000:AMITEC>2.0.CO;2

  • 80

    ThompsonD. W. J.SolomonS. (2002). Interpretation of recent Southern Hemisphere climate change. Science296, 895899. 10.1126/science.1069270

  • 81

    TorrenceC.CompoG. P. (1998). A ractical guide to wavelet analysis. Bull. Am. Meteorol. Soc.79, 6178.

  • 82

    TyrlisE.HoskinsB. J. (2008). Aspects of a Northern Hemisphere atmospheric blocking climatology. J. Atmosph. Sci.65, 16381652. 10.1175/2007JAS2337.1

  • 83

    van den BroekeM. R.van LipzigN. 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, 485496. 10.3402/tellusa.v54i5.12161

  • 84

    van LoonH.RogersJ. C. (1978). The seesaw in winter temperatures between Greenland and northern Europe. Part I: General description. Mon. Weather Rev.106, 296310.

  • 85

    VanhataloE.KulahciM. (2016). Impact of autocorrelation on principal components and their use in statistical process control. Quality Reliab. Eng. Int.32, 14831500. 10.1002/qre.1858

  • 86

    VanosJ. K.CakmakS. (2014). Changing air mass frequencies in Canada: Potential links and implications for human health. Int. J. Biometeorol.58, 121135. 10.1007/s00484-013-0634-2

  • 87

    VavrusS. J. (2018). The influence of Arctic amplification on mid-latitude weather and climate. Current Climate Change Rep.4, 238249. 10.1007/s40641-018-0105-2

  • 88

    VavrusS. J.WangF.MartinJ. E.FrancisJ. A.PeingsY.CattiauxJ. (2017). Changes in North American atmospheric circulation and extreme weather: Influence of Arctic amplification and Northern Hemisphere snow cover. J. Clim.30, 43174333. 10.1175/JCLI-D-16-0762.1

  • 89

    WallaceJ. M.GutzlerD. S. (1981). Teleconnections in the geopotential height field during the Northern Hemisphere winter. Mon. Weather Rev.109, 784812.

  • 90

    WangZ. Y.DingY. 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, 2029.

  • 91

    WaughD. W.SobelA. H.PolvaniL. M. (2017). What is the polar vortex and how does it influence weather?Bull. Am. Meteorol. Soc.98, 3744. 10.1175/BAMS-D-15-00212.1

  • 92

    WelchP. (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, 7073. 10.1109/TAU.1967.1161901

  • 93

    WronaK. M.RohliR. V. (2007). Seasonality of the northern hemisphere circumpolar vortex. Int. J. Climatol.27, 697713. 10.1002/joc.1430

  • 94

    ZakinyanR.ZakinyanA.RyzhkovR. (2016). Phases of the isobaric surface shapes in the geostrophic state of the atmosphere and connection to the polar vortices. Atmosphere7, 126. 10.3390/atmos7100126

Summary

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

Volume

5 - 2022

Edited by

Letizia Lusito, Inside Climate Service, Italy

Reviewed by

Antonio Speranza, University of Camerino, Italy; Anthony Lupo, University of Missouri, United States

Updates

Copyright

*Correspondence: Nazla Bushra ✉

This article was submitted to Data-driven Climate Sciences, a section of the journal Frontiers in Big Data

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics