The Dynamic Character of Northern Hemisphere Flow Regimes in a Near-Term Climate Change Projection

The dynamic character of an enstrophy-based diagnostic, previously used in the study of atmospheric blocking, is examined here, in near-term future simulations from the Institut Pierre Simon Laplace Climate Model version 4 (IPSL-CM4) and version 5 (IPSL-CM5) climate models of the Northern Hemisphere flow for moderate climate change scenarios. Previous research has shown that integrated regional enstrophy (IE) increases during blocking onset and decay, which is a reflection of planetary-scale instability. In addition, IE has been shown previously to increase during flow regime transitions in general, even those not associated with blocking events. Here, a 31-year IE diagnostic time series is examined for changes in short term (5–40 days) planetary-scale variability that may correspond flow regime changes in an increased carbon dioxide environment. The time-series analysis herein indicates that the IE diagnostic provides evidence for approximately 30–35 atmospheric flow regime transitions per year in a warmer climate, which is similar to that of the control run and the latest 30-year observed climate, as derived from re-analyses. This result has implications regarding the predictability of weather in a warmer world.


Introduction
The idea that Earth's atmosphere has preferred states or that certain states are persistent, at least in a qualitative sense, has been noted in the atmospheric science literature for more than 75 years [1].Later research found correlations between the large-scale flow pattern over the north Pacific and the eastern U.S. [2,3].At the same time, researchers identified persistent positive pressure anomalies that occurred in preferred regions of the Northern Hemisphere (e.g., [4]), called blocking anticyclones.Then, the work of [5] demonstrated, using a low-order primitive equation model for convection in a closed hydrodynamic system, that two flows which are initially similar may evolve very differently over time, or that predictability beyond a certain time frame for a system like the atmosphere is essentially impossible using the primitive equations.The solutions to the low-order model in [5] extreme (RCP 8.5) anthropogenic emissions scenarios.However, [26] notes increased inter-annual variability in their projected blocking climatologies.Additionally, a comparison of [30] and the model studies cited here demonstrates the well-known tendency for the underrepresentation of blocking in climate models.
In [31] and references therein, integrated enstrophy (IE) first published by [32] is demonstrated to be useful in identifying the onset and decay period for blocking.However, [33] states that this diagnostic may not be unambiguously associated only with blocking regime transitions, but is likely associated flow regime transitions in general.Also, [34] used IE to identify seasonal flow regime transitions in their study, and showed that the behavior of the IE diagnostic on the seasonal time scale is similar to that of the index in [6,7].
Then, the characteristics of the IE diagnostic in climate model simulations for the period 2020-2050 were obtained from simulations in the IPSL-CM4 and IPSL-CM5 models.In particular, hemispheric-scale flow regime transitions were identified in the 5-40 day range, which is consistent with the time period for these types of transitions (e.g., [6,7]).The IE diagnostic used to identify flow regime transition is related also to finite-time instability [32,35].While these methods have been applied previously to past blocking events, here an important objective is to apply the methodology (time-series analysis of IE) to a reanalysis data set and climate model simulations, in order to study flow regime transition.The authors are not aware of published studies that perform a similar analysis with IE.The results of time-series analysis indicate that the IE diagnostic supports approximately 30-35 large-scale flow regime transitions (not necessarily blocking as discussed above) per year in the current and in a warmer climate.The paper is organized as follows.The relevant theory is described in section two, where the methods developed in [32,36] are outlined.Section three contains our findings, and section four includes a discussion of the results and conclusions.

Model and Reanalysis Data
Simulations from the coupled ocean-atmosphere general circulation models, the IPSL-CM4 [27], with horizontal resolution 2.5 • × 3.75 • (lat/lon) under the SRES-A1B, and the IPSL-CM5-LR [37] with horizontal resolution 1.875 • × 3.75 • under a present-day control scenario (from IPSL-CM5) and the RCP 4.5 (intermediate emission) scenario, were the models and scenarios used here.The near-future scenarios covered the 31-year period from 2020 to 2050.The model was chosen because the data were available readily, and because this model is part of the Coupled Model Intercomparison Project (CMIP) suite [37].Also, the IPSL-CM5 model featured improvements, including those in the boundary layer physics, an improved ocean model, and air-sea ice interaction scheme, as well as improvements in the atmospheric chemistry.A complete review of improvements can be found in [37].Since the intermediate scenario was chosen as a mid-range scenario, it may be representative of the range of scenarios produced by [28].The variable used to calculate the IE diagnostic from each model and scenario were the Northern Hemisphere 500 hPa height (m).
In order to calculate the spectral characteristics of the daily IE field, obtained from the Northern Hemisphere 500 hPa heights for the years 2020-2050 from both versions of the IPSL model, a Fast Fourier Transform (FFT) was applied to the daily IE time series.A bandpass filter with was employed to isolate periods of approximately 5-40 days in the IE time series, which correspond to the general time-scale of planetary-scale flows.Additionally, [38] and references therein demonstrate that the planetary scale using spatial filtering techniques corresponds to low wave numbers, up to four or five, which have a time-scale of about six days as inferred from [39,40] or others.For the wavelet analysis, the Morlet wavelet was used as the wave basis, and this has been used in published climatological studies [41][42][43].The data were smoothed, using a moving average to remove seasonality, and the linear trend was subtracted before obtaining the power spectrum.The spectral density was calculated at the 95% confidence level.IE was also calculated from the set of Northern Hemisphere re-analyses, as well for comparison to the model output.The 500 hPa heights were obtained from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) re-analyses, which were archived on a 2.5 • × 2.5 • latitude/longitude grid from 1981 to 2010.

Methods
Here the IE diagnostic of [32] is introduced.Local Lyapunov exponents from the barotropic vorticity equation are defined in [27] as λ i (ξ 0 , T) = 1 2n log ν i for initial vorticity ξ 0 and time T = n ∆t.The ν i are the eigenvalues of M * n M n , where M n = ∏ k=n k=−n B(k∆t), and B(t) is the Jacobian operator, and n is a multiple of the time increment ∆t, and k an index in the product of M*M.The study of [32] also showed that the sum of the positive local Lyapunov exponents is highly correlated with the integral of enstrophy.This integral may be used as a stability indicator, and is defined as in [36]: In (1), A is the area and ξ is the relative vorticity; detailed explanations and derivations of this quantity can be found in [32,44].Briefly, the relative vorticity was calculated as the geostrophic value using the height field and second order finite differencing.Height fields were used, since they have less error than wind fields.This value can be integrated over the entire NH [32].
The IE diagnostic is based on the premise that changes in the stability in the Northern Hemisphere flow have been shown to coincide with large-scale flow regime transition as in [45].The time series of IE can be interpreted as follows: (a) peaks (local maxima or minima) in the IE time series are associated with a change in sign of the time derivative of IE as in [46].The authors of [46] also showed that the change in the sign of the time derivative of this quantity from positive to negative, associated with an IE maximum, indicated a local maximum in instability.Based on previous research ( [33,36]), instability is expected at or near block onset and decay ( [45,46]).The study of [45] argued that blocking existed within a large-scale flow regime.In their statistical analysis, they argued that the minimum in the probability density distribution between two large-scale flow regimes represented instability associated with regime change.This is consistent with the view of [46] and here.
Additionally, [46] provided an explanation for the reason large-scale flow and blocking regime changes are correlated with extrema in the IE time series.They showed that enstrophy can be written as: In (2) f, φ, β, σ, and u are the Coriolis parameter, geopotential height, the "Beta" effect (meridional change in Coriolis parameter), deformation (shearing and stretching), and the zonal wind, respectively.In other words, extrema in the geopotential height (the horizontal Laplacian term), as well as relative vorticity, zonal wind, and resultant deformation, all contribute to the enstrophy budget.For example, local maxima in the geopotential height field within a ridge or blocking event decrease the enstrophy.Then, with the changes in geopotential height as the large-scale flow regime breaks down, or when a blocking event decays, enstrophy is increased.The fields of westerlies and relative vorticity are known also to be different in higher-amplitude flows versus zonal flows.The study of [46] also showed that fluxes of enstrophy across a fluid boundary are balanced by the time rate of change of deformation, which leads to changes in the enstrophy as shown by: In Equation ( 3), C, v h , n, and ds are the fluid boundary, the horizontal wind (flow), the unit normal, and the change in position on C, respectively.

Results
The raw time series of the modelled RCP 4.5 Northern Hemisphere IE for 2020-2050 is shown in Figure 1 (IPSL-CM4 Figure 1a, IPSL-CM5 Figure 1b).Time series from both models exhibit an approximate annual cycle, which is a reflection of seasonality in the distribution of the 500 hPa height field.To determine if there are other significant spectral peaks, apart from the clear annual periodicity in IE, as indicated in Figure 1, a bandpass filter was used to filter out frequencies outside the 5-40 day range.The power spectra (both Fourier and Wavelet) obtained after filtering, smoothing, and detrending the time series from both model simulations are shown in Figures 2 and 3.In Figure 2, the power spectra for the IPSL-CM4 time series of IE have a clear spectral peak in the 10-12 day range.This implies that a hemispheric flow regime persists for about 10-12 days before a transition occurs, as identified by the IE diagnostic.The same result is evident for the IPSL-CM5 time series (Figure 3).This persistence is consistent with the results of [6,7] and others for large-scale flow regimes.
A comparison of the near future model RCP 4.5 time series (Figures 2 and 3) with the control model, using the IPSL-CM5 as an example, shows a smaller peak near 10-12 days (Figure 4a) consistent with the model simulations.However, there is a large peak in Figure 4a, near 20-30 days, that is apparent in the IPSL-CM5 control run (but not the IPSM-CM4).It is not immediately obvious why the peak at 20-30 days is strong in the control model simulation for the IPSL-CM5, and a possible explanation will be discussed in section four below.

Results
The raw time series of the modelled RCP 4.5 Northern Hemisphere IE for 2020-2050 is shown in Figure 1 (IPSL-CM4 Figure 1a, IPSL-CM5 Figure 1b).Time series from both models exhibit an approximate annual cycle, which is a reflection of seasonality in the distribution of the 500 hPa height field.To determine if there are other significant spectral peaks, apart from the clear annual periodicity in IE, as indicated in Figure 1, a bandpass filter was used to filter out frequencies outside the 5-40 day range.The power spectra (both Fourier and Wavelet) obtained after filtering, smoothing, and detrending the time series from both model simulations are shown in Figures 2 and 3.In Figure 2, the power spectra for the IPSL-CM4 time series of IE have a clear spectral peak in the 10-12 day range.This implies that a hemispheric flow regime persists for about 10-12 days before a transition occurs, as identified by the IE diagnostic.The same result is evident for the IPSL-CM5 time series (Figure 3).This persistence is consistent with the results of [6,7] and others for large-scale flow regimes.
A comparison of the near future model RCP 4.5 time series (Figures 2 and 3) with the control model, using the IPSL-CM5 as an example, shows a smaller peak near 10-12 days (Figure 4a) consistent with the model simulations.However, there is a large peak in Figure 4a, near 20-30 days, that is apparent in the IPSL-CM5 control run (but not the IPSM-CM4).It is not immediately obvious why the peak at 20-30 days is strong in the control model simulation for the IPSL-CM5, and a possible explanation will be discussed in section four below.The autocorrelation (see Figure 5a) of the IPSL-CM4 RCP 4.5 time series decays rapidly, but there is a "bump" at a lag correlation of approximately 10 days, as suggested by the spectral peak in Figure 2a.In Figure 3, the power spectrum for the IPSL-CM5 time series of the IE has a spectral peak at approximately 10-12 days, as shown above, but also has peaks beyond the 10-day period (around 20-30 days) that were not captured as strongly in the IPSL-CM4 time series of the IE.However, the autocorrelation in the time series for both versions of the IPSL model (Figure 5) are qualitatively similar.The autocorrelation of the control model simulation, as represented by IPSL-CM4 (Figure 4b), is also consistent with the RCP 4.5 IPSL-CM4 and CM5 model (Figure 5).
A comparison of the IE derived from the NCEP/NCAR re-analyses (Figure 6) for the period 1981-2010 demonstrates that the IPSL-CM5 represents the character of the large-scale dynamic flow regime properties more faithfully than that of the IPSL-CM4, and this is true for the control runs as well (Figure 5).These will be discussed further below.Additionally, there were no qualitative differences noted in the seasonal cycle between Figure 1 and the control run or NCEP re-analyses.The autocorrelation (see Figure 5a) of the IPSL-CM4 RCP 4.5 time series decays rapidly, but there is a "bump" at a lag correlation of approximately 10 days, as suggested by the spectral peak in Figure 2a.In Figure 3, the power spectrum for the IPSL-CM5 time series of the IE has a spectral peak at approximately 10-12 days, as shown above, but also has peaks beyond the 10-day period (around 20-30 days) that were not captured as strongly in the IPSL-CM4 time series of the IE.However, the autocorrelation in the time series for both versions of the IPSL model (Figure 5) are qualitatively similar.The autocorrelation of the control model simulation, as represented by IPSL-CM4 (Figure 4b), is also consistent with the RCP 4.5 IPSL-CM4 and CM5 model (Figure 5).
A comparison of the IE derived from the NCEP/NCAR re-analyses (Figure 6) for the period 1981-2010 demonstrates that the IPSL-CM5 represents the character of the large-scale dynamic flow regime properties more faithfully than that of the IPSL-CM4, and this is true for the control runs as well (Figure 5).These will be discussed further below.Additionally, there were no qualitative differences noted in the seasonal cycle between Figure 1 and the control run or NCEP re-analyses.As explained in section two, the IE achieves a relative maximum during flow regime transitions.The results of the spectral analysis of both models indicate that flow regimes persist for approximately 10 days, which suggests that the IE diagnostic may support, on average, up to approximately 35 flow regime transitions per year in a warmer climate.This number (35) is greater than the number of blocking events found in [26,29].In spite of this, both studies found that the average duration for blocking is close to 10 days in the IPSL-CM4 and IPSL-CM5 simulations.As discussed in sections one and two, not every flow regime can be associated with the occurrence of blocking, since it has been shown that occasionally blocking occurs simultaneously in two different locations of the Northern Hemisphere [30,48].This is an important reason for the difference in the annual average number of blocking events in the climatologies of [26,29] (approximately 22 events), and the approximately 30-35 flow regime transition events found here.As explained in section two, the IE achieves a relative maximum during flow regime transitions.The results of the spectral analysis of both models indicate that flow regimes persist for approximately 10 days, which suggests that the IE diagnostic may support, on average, up to approximately 35 flow regime transitions per year in a warmer climate.This number (35) is greater than the number of blocking events found in [26,29].In spite of this, both studies found that the average duration for blocking is close to 10 days in the IPSL-CM4 and IPSL-CM5 simulations.As discussed in sections one and two, not every flow regime can be associated with the occurrence of blocking, since it has been shown that occasionally blocking occurs simultaneously in two different locations of the Northern Hemisphere [30,48].This is an important reason for the difference in the annual average number of blocking events in the climatologies of [26,29] (approximately 22 events), and the approximately 30-35 flow regime transition events found here.

Discussion and Conclusions
This study examines the behavior of an enstrophy-based diagnostic (IE) and Northern Hemisphere flow regime dynamics described above, in two near-future moderate emission simulations from the IPSL-CM4 and IPSL-CM5 GCMs.The use of the IE diagnostic to study flow regime dynamics in the NCEP/NCAR re-analyses and climate simulations had not been performed previously.The time series of IE for the re-analyses, control simulations, and the 2020-2050 simulations does have an annual cycle.After filtering out low and high frequencies, there is a clear peak of IE at approximately 10-12 days in both the IPSL-CM4 and IPSL-CM5, and this corresponds well to the control run and re-analyses from the most recent 30-year period (1981-2010) derived from observations.Flow regime transition means that a particular atmospheric flow regime persists for 10-12 days, before transitioning to another flow regime.The spectral density peaks also in the 20-30 day range, especially in the NCEP/NCAR re-analyses and the IPSL-CM5 model as shown using FFT analysis.This spectral peak was not statistically significant in the wavelet analysis.Based on previous research ( [33,36,46] and others), the IE reaches a maximum during large-scale flow regime transitions, or at the onset and demise of blocking events.This work supported 30-35 flow regime transitions per year in the large-scale Northern Hemisphere flow for the IPSL-CM4 and CM5 simulations, while the annual mean number of blocking events in these same data sets was fewer [26][27][28][29].
Since these results imply that the IE indicates approximately 30-35 flow regime transitions per year on average in a near-future warmer climate, the implication is that the well-known dynamic predictability wall for large-scale flow (e.g., [49,50]) would still be similar to the present day.Further, the implication is that the dynamic character of the near-future large-scale flow regime in a moderate CO2 emissions scenario would be similar to that for the immediate past 30 years, or that the "predictability wall" for short-term weather forecasting does not change appreciably.In a future climate, weather and climate prediction beyond 10-12 days is still only possible using primarily statistical methodologies.While the dynamic character of the large-scale flow, as determined using the IE diagnostic, may not change significantly in the near-future, this does not preclude further changes in regional or global sensible climate that is expected to occur [28].
Support for this conclusion might be inferred from modern studies of jet stream behavior under increased CO2 conditions (e.g., [25,51,52]).All three of these studies [25,51,52] used the CMIP5 model ensembles, including the IPSL-CM4 and CM5 runs, in the RCP4.5 and RCP8.5 scenarios.In [25] they found a poleward shift for the Atlantic jet, and less variability in the latitudinal position of the jet in the increased CO2 runs (but there was also less variability in the control runs).This is supported by [51,52], which also found that the Atlantic basin jet stream included less positional variability, and what they described as more "pulsing" behavior (i.e., stronger variations in scalar wind speed).They [51,52] found the opposite result in the Pacific Region, and speculated that the connection between positional variability versus "pulsating type" behavior may be Rossby wave breaking.Thus, we speculate, based on these studies, that the hemispheric-scale "preferred mode" behavior can be

Discussion and Conclusions
This study examines the behavior of an enstrophy-based diagnostic (IE) and Northern Hemisphere flow regime dynamics described above, in two near-future moderate emission simulations from the IPSL-CM4 and IPSL-CM5 GCMs.The use of the IE diagnostic to study flow regime dynamics in the NCEP/NCAR re-analyses and climate simulations had not been performed previously.The time series of IE for the re-analyses, control simulations, and the 2020-2050 simulations does have an annual cycle.After filtering out low and high frequencies, there is a clear peak of IE at approximately 10-12 days in both the IPSL-CM4 and IPSL-CM5, and this corresponds well to the control run and re-analyses from the most recent 30-year period (1981-2010) derived from observations.Flow regime transition means that a particular atmospheric flow regime persists for 10-12 days, before transitioning to another flow regime.The spectral density peaks also in the 20-30 day range, especially in the NCEP/NCAR re-analyses and the IPSL-CM5 model as shown using FFT analysis.This spectral peak was not statistically significant in the wavelet analysis.Based on previous research ( [33,36,46] and others), the IE reaches a maximum during large-scale flow regime transitions, or at the onset and demise of blocking events.This work supported 30-35 flow regime transitions per year in the large-scale Northern Hemisphere flow for the IPSL-CM4 and CM5 simulations, while the annual mean number of blocking events in these same data sets was fewer [26][27][28][29].
Since these results imply that the IE indicates approximately 30-35 flow regime transitions per year on average in a near-future warmer climate, the implication is that the well-known dynamic predictability wall for large-scale flow (e.g., [49,50]) would still be similar to the present day.Further, the implication is that the dynamic character of the near-future large-scale flow regime in a moderate CO 2 emissions scenario would be similar to that for the immediate past 30 years, or that the "predictability wall" for short-term weather forecasting does not change appreciably.In a future climate, weather and climate prediction beyond 10-12 days is still only possible using primarily statistical methodologies.While the dynamic character of the large-scale flow, as determined using the IE diagnostic, may not change significantly in the near-future, this does not preclude further changes in regional or global sensible climate that is expected to occur [28].
Support for this conclusion might be inferred from modern studies of jet stream behavior under increased CO 2 conditions (e.g., [25,51,52]).All three of these studies [25,51,52] used the CMIP5 model ensembles, including the IPSL-CM4 and CM5 runs, in the RCP4.5 and RCP8.5 scenarios.In [25] they found a poleward shift for the Atlantic jet, and less variability in the latitudinal position of the jet in the increased CO 2 runs (but there was also less variability in the control runs).This is supported by [51,52], which also found that the Atlantic basin jet stream included less positional variability, and what they described as more "pulsing" behavior (i.e., stronger variations in scalar wind speed).They [51,52] found the opposite result in the Pacific Region, and speculated that the connection between positional variability versus "pulsating type" behavior may be Rossby wave breaking.Thus, we speculate, based on these studies, that the hemispheric-scale "preferred mode" behavior can be expressed as variability in the position or speed of the jet stream (meridional or zonal), or both.The IE diagnostic is integrated over the whole of the Northern Hemisphere, and thus will reflect the changes in the combined variability across both the Atlantic and Pacific Ocean basins.This may be because enstrophy is the square of the vorticity, and in natural coordinates vorticity can be expressed as the sum of the speed shear and curvature.Thus, if in an increased CO 2 atmosphere the preferred mode behavior is represented more by a "pulsating" jet, with less positional variability (at least in the Atlantic), the IE diagnostic should be able to capture this behavior as well.
However, it is conceded here that the discussion in relation to the studies above is within the context of projections for the end of the 21st century ( [51,52]).Nonetheless, [52] shows near-term (2020-2044) scenarios in relation to arctic amplification, and the trend toward a more poleward and less variable jet stream is evident, but less than that for the end-of-the-century projections.It is presumed the more "pulsating" type behavior described these studies is also occurring, though it is difficult to determine from the near-term statistics presented by [52].Additionally, it is noted here that [53] demonstrated that lower resolution models may not represent regional large-scale flow regimes and transitions as faithfully as much higher resolution simulations.
The spectral peak in the modelled IE time series, at approximately 10-12 days, compares favorably with the results of [26,29] for the projected duration of blocking in the mid-to-late 21st century, for the IPSL-CM4 and IPSL-CM5 simulations the peaks were close to 10 days.The model spectral peaks were also consistent with the re-analysis 500-hPa height data, as well as the control run.However, as discussed in section three, large-scale atmospheric flow regimes may occur with the absence of blocking, and at other times with simultaneously occurring blocking events.Simultaneous blocking does occur about nine percent of the time ( [30,48]) in the Northern Hemisphere.This is an important reason for the difference in the average number of blocking events in the climatology of [26,29] (approximately 22 events), and the result of approximately 30-35 flow regime transition events, as found here.Additionally, this work did not focus exclusively on atmospheric blocking, as there are several published studies on the future occurrence of these events in increased CO 2 climate scenarios, some of which show conflicting results e.g., [29,54].
Lastly, there was a strong spectral peak in the 20-30 day time period in the control model run that is also evident in the IPSL-CM5 simulation and the re-analyses.This peak is evident in the autocorrelation of the IE figures as well.While this time period is beyond the scope of our work here, predictability in the 20-30 day period was discussed in a recent study [55].That study found statistically significant peaks in this time frame, using the FFT and autocorrelation analysis applied to a 67-year time series of the Pacific North America Index [55].They attributed these to persistent and quasi-stationary Rossby wave trains occurring within the Pacific Ocean basin, and this information can be used to make weather predictions in the 2-4 week time period.

Figure 2 .
Figure 2. The Integrated Regional Enstrophy (IE) (× 10 −6 m 2 s −2 -ordinate) diagnostic derived from the Northern Hemisphere 500 hPa heights for the (a) band passed (5-40 day) Fast Fourier Transform (FFT) power spectra (ordinate) versus spectral period (days-abscissa) (top), and (b), and (c) the wavelet power spectra (time period is days) (bottom).In (b) and (c) the dashed lines represent statistical significance at 90 and 95%, and regions where edge effects become important [47], respectively.Only peaks with statistically significant power are shown in (c).

Figure 2 .
Figure 2. The Integrated Regional Enstrophy (IE) (× 10 −6 m 2 s −2 -ordinate) diagnostic derived from the Northern Hemisphere 500 hPa heights for the (a) band passed (5-40 day) Fast Fourier Transform (FFT) power spectra (ordinate) versus spectral period (days-abscissa) (top), and (b), and (c) the wavelet power spectra (time period is days) (bottom).In (b) and (c) the dashed lines represent statistical significance at 90 and 95%, and regions where edge effects become important [47], respectively.Only peaks with statistically significant power are shown in (c).

Figure 3 .
Figure 3. Same as in Figure 2, except with the IPSL-CM5 model power spectra.

Figure 4 .
Figure 4.The control model run (IPSL-CM5) for (a) IE band passed (5-40 day) FFT power spectra (ordinate) versus spectral period (days-abscissa) (top), and (b) the autocorrelation of the IE diagnostic (ordinate), where the lag is presented in days (abscissa).The blue dashed line is the zero correlation line, and the solid line is the 95% confidence level.

Figure 4 .
Figure 4.The control model run (IPSL-CM5) for (a) IE band passed (5-40 day) FFT power spectra (ordinate) versus spectral period (days-abscissa) (top), and (b) the autocorrelation of the IE diagnostic (ordinate), where the lag is presented in days (abscissa).The blue dashed line is the zero correlation line, and the solid line is the 95% confidence level.

Figure 4 .
Figure 4.The control model run (IPSL-CM5) for (a) IE band passed (5-40 day) FFT power spectra (ordinate) versus spectral period (days-abscissa) (top), and (b) the autocorrelation of the IE diagnostic (ordinate), where the lag is presented in days (abscissa).The blue dashed line is the zero correlation line, and the solid line is the 95% confidence level.

Figure 6 .
Figure 6.Same as in Figures 2 and 4, except using the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) re-analyses for the period 1981-2010.The top, middle, and bottom diagrams are the Fourier, Wavelet, and autocorrelation diagrams, respectively.

Figure 6 .
Figure 6.Same as in Figures 2 and 4, except using the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) re-analyses for the period 1981-2010.The top, middle, and bottom diagrams are the Fourier, Wavelet, and autocorrelation diagrams, respectively.