Journal topic
Ann. Geophys., 38, 545–555, 2020
https://doi.org/10.5194/angeo-38-545-2020
Ann. Geophys., 38, 545–555, 2020
https://doi.org/10.5194/angeo-38-545-2020

Regular paper 21 Apr 2020

Regular paper | 21 Apr 2020

# A note on the statistical evidence for an influence of geomagnetic activity on Northern Hemisphere seasonal-mean stratospheric temperatures using the Japanese 55-year Reanalysis

A note on the statistical evidence for an influence of geomagnetic activity on Northern Hemisphere seasonal-mean stratospheric temperatures using the Japanese 55-year Reanalysis
Nazario Tartaglione1,2, Thomas Toniazzo1,2, Yvan Orsolini2,3,4, and Odd Helge Otterå1,2 Nazario Tartaglione et al.
• 1NORCE Climate, Bergen, Norway
• 2Bjerkness Centre for Climate Research, Bergen, Norway
• 3Birkeland Centre for Space Science, University of Bergen, Bergen, Norway
• 4Norwegian Institute for Air Research, Kjeller, Norway

Correspondence: Nazario Tartaglione (nazario.tartaglione@norceresearch.no)

Abstract

We employ JRA-55 (Japanese 55-year Reanalysis), a recent second-generation global reanalysis providing data of high quality in the stratosphere, to examine whether a distinguishable effect of geomagnetic activity on Northern Hemisphere stratospheric temperatures can be detected. We focus on how the statistical significance of stratospheric temperature differences may be robustly assessed during years with high and low geomagnetic activity. Two problems must be overcome. The first is the temporal autocorrelation of the data, which is addressed with a correction of the t statistics by means of the estimate of the number of independent values in the series of correlated values. The second is the problem of multiplicity due to strong spatial autocorrelations, which is addressed by means of a false discovery rate (FDR) procedure. We find that the statistical tests fail to formally reject the null hypothesis, i.e. no significant response to geomagnetic activity can be found in the seasonal-mean Northern Hemisphere stratospheric temperature record.

1 Introduction

There is a large interest in the potential climate impact of geomagnetic activity. One of the main mechanisms by which geomagnetic activity is thought to affect the middle atmosphere is through the production of nitrogen oxides (NOx), either by the continuous precipitation of auroral electrons penetrating into the lower thermosphere or by the more episodic precipitation of higher energy electrons into the mesosphere . Downward transport from the mesosphere to the stratosphere in winter results in the increased availability of NOx in the dark polar stratosphere, where it is long lived. NOx can catalytically reduce ozone concentrations as the Sun returns , and thus alter radiative heating rates, with potential observable impacts on stratospheric temperatures and possible implications also for surface air temperature (SAT). The amount of NOx in the middle atmosphere during late winter and spring depends on the cumulative effect of geomagnetic activity over the preceding months on the NOx reservoir (Jacob1999). Stratospheric NOx concentrations however also depend on the magnitude of the downward transport from this reservoir and are thereby affected by internal variability of the atmospheric circulation from year to year, especially in the Northern Hemisphere  .

The impact of energetic electron precipitation (EEP) driven by geomagnetic activity on NOx and ozone concentrations has been well documented after detailed satellite studies were carried out in the early 2000 . Several recent studies suggest a significant signal associated with geomagnetic activity in the observed climate. However, there remains considerable uncertainty regarding the precise attribution of such a signal, and the existence of a direct link between EEP and stratospheric and tropospheric temperatures has remained controversial. Among other things, the study of , henceforth S09, in particular, claims to find a significant, direct relationship between the SAT and geomagnetic activity based on reanalysis data from the European Centre for Medium Range Weather Forecasts (ECMWF). In essence, S09 finds that the hypothesis that geomagnetic activity influences the SAT is supported by reanalysis data, whereas the null hypothesis that the SAT is not influenced by the geomagnetic activity at all is rejected. S09 compares the seasonal SAT in years with high and low geomagnetic activity and also considered the separate effect of the variation in solar irradiance associated with the 11-year heliomagnetic cycle.

The selection of years in S09 was based on two indices, Ap and f10.7. Ap (Rostoker1972) provides a measure for the daily average level of geomagnetic activity. To account for the cumulative effect of NOx production, transport and diffusion processes, Ap was commonly averaged over 4 months from late autumn to winter . In particular, S09 used Ap averaged between October and January to define winters of high and low geomagnetic activity in the Northern Hemisphere. The second index, f10.7 (https://www.swpc.noaa.gov/phenomena/f107-cm-radio-emissions, last access: 13 March 2020), is an indicator of the phase and intensity of the solar cycle. By compositing separately on the basis of Ap and f10.7, S09 obtained different samples of seasonal-mean data for years with high geomagnetic activity and for years with low geomagnetic activity. They then computed the SAT differences of the seasonal means (December, January and February – DJF; March, April and May – MAM; June, July and August – JJA; and September, October and November – SON) between the two samples and employed a t test based on the set of daily means (Annika Seppälä, personal communication, 2018) used to compute the seasonal averages to discriminate against a null hypothesis of no effect. As a consequence of such a procedure, S09's claim of significance is marred by the presence of very strong temporal and spatial autocorrelation within the samples.

In this paper, we revisit the S09 hypothesis by adopting a rigorous methodology for significance testing on strongly autocorrelated data. We focus on wintertime stratospheric temperatures between 200 and 1 hPa, a prerequisite for possible surface impacts associated with EEP-related changes in ozone concentrations. Although our analysis is focused, in this paper on stratospheric temperature, we look at all the levels present in the dataset. We show that statistical testing appropriate to the data at hand is a crucial step in any analysis purporting to demonstrate an observed climate signal of geomagnetic activity.

Data and methods are described in Sect. 2, including a discussion on the problem of autocorrelation in time and space. In Sect. 3, the results obtained by applying the t test to the stratospheric temperatures are shown. The analysis is applied to four different cases: with no correction at all, with the temporal and the spatial autocorrelation correction applied separately, and with both the corrections applied. In Sect. 4 conclusions are drawn.

2 Data and methods

## 2.1 Data

To analyse the possible impact of geomagnetic activity in the stratosphere, we use the Japanese 55-year Reanalysis (JRA-55) covering more than 55 years, extending from 1958 to the present . Due to the selection of cases of high and low geomagnetic activity as in S09, only data up to 2006 are used here. In the JRA-55 reanalysis, ozone is used interactively in the radiation code, although it is treated differently in the pre- and post-1979 satellite era. This is an important asset for JRA-55, since the EEP will primarily affect NOx and ozone, and this feature is not commonly found in other reanalysis systems, such as the ECMWF reanalyses. Older-generation reanalyses tend to suffer from temporal inhomogeneities because of the sequential introduction of new satellite data during the assimilation period, especially in the Southern Hemisphere (SH) as shown recently by . For these various reasons, we restricted our analysis to the recent JRA-55 reanalysis. also used the JRA-55 reanalyses to investigate the signature of geomagnetic activity but focused exclusively on the SH. He found a temperature signal in the upper stratosphere, but only in July. The S09 selection shown in Table 1 is used to compute the significance of the seasonal differences. The criteria used to select the different years are based on the Ap and f10.7 values and are the same as used by S09. The definition of high and low geomagnetic activity is the same as S09. We hence investigate the potential signatures on stratospheric temperature during the same winters and in the following seasons of the same calendar year as S09 did for the SAT. The set of data is denominated N1 as in S09 (Table 1).

Table 1Years used to define the N1 set, following S09.

## 2.2 Data autocorrelation and statistical significance

S09 computed the SAT differences of the seasonal means (DJF, MAM, JJA and SON) between those selected high-Ap and low-Ap years and employed Welch's t test (hereafter only t test) to assess the likelihood of the differences given a null hypothesis of no effect. Such a test assumes a statistical model in which observations are normally distributed and statistically independent. In particular, the t test is sensitive to the temporal autocorrelation or serial correlation within the samples. When serial correlation is not taken into account in the data, statistically significant differences in two means, which may not be different at all, are found more frequently than expected . S09's analysis is affected by this problem, because the authors used daily-mean data in their t test, which are highly autocorrelated in time. As seasonal averages can still suffer from temporal autocorrelation, the serial dependence is checked by means of the Durbin–Watson test . While the serial correlation, in general, is reduced from seasonal averaging, it can still persist, especially in summer. To deal with such serial correlation, a correction is applied as suggested by . The temporal autocorrelation is not the only potential caveat that needs to be considered when testing a hypothesis. When performing a significance test simultaneously on many samples, one will at some point find statistically significant temperature differences simply by accident. Unfortunately, the dominant approach to the multiplicity problem is generally to test the single grid points and then to report them as “significant” when the null hypothesis is locally rejected (Wilks2016). Sometimes temporal and spatial autocorrelation is not addressed at all, but, there are some exceptions. , for instance, while trying to find a relationship between solar activity and surface air temperature, dealt with temporal and spatial autocorrelation using a Monte Carlo approach. To overcome this multiplicity problem in our analysis, we apply the false discovery rate controlling introduced by and proposed in the atmospheric sciences by .

## 2.3 Accounting for temporal autocorrelation

The t test is a widely used method for hypothesis testing within the climate community. It is however well known that the t test assumes a statistical model where observations are statistically independent, and it is widely, but incorrectly, believed that the t test is valid only for normally distributed outcomes. Several authors have shown that the t test is suitable under symmetric, not necessarily normal and asymmetric distributions. The t test is sensitive to time autocorrelation or serial correlation within the samples. The effect of serial correlation is, usually, to make comparisons of means which are too liberal. The null hypothesis, assuming equal means, is hence rejected more frequently than expected. Two separate reasons favour the use of seasonal-mean data instead of daily-mean data. The first reason is that any influence of EEP on temperature is expected to accumulate over seasonal timescales. The second reason is that daily temperatures are strongly serially correlated, whereas seasonal data have less correlation between 2 consecutive years for instance. In fact, one of the causes of the serial correlation is that the variable of interest varies seasonally. Nevertheless, even for seasonal means, it is important to account for serial correlations, as there may be other causes leading temporal autocorrelation, including persistence. Figure 1a shows the results of the Durbin–Watson test applied at the seasonal temperatures at 5 hPa. Similar pictures can be obtained by plotting the lag-1 autocorrelation (Fig. 1b), but the Durbin–Watson test, which is a classical test to check whether data are serially correlated, is better, compared to including the lagged response, as it tests for autocorrelation in the residuals, and it is suitable when in time series there are trends or seasonal patterns. When data are serially correlated, the test gives values close to zero, whereas when data are not correlated, the test statistic values, as a rule of thumb, are in the range of 1.5 to 2.5. There is also the possibility of serial anticorrelation: in such a case, the value would be above 2.5, but this situation was not found in our study.

Figure 1Results of Durbin–Watson test (a) and lag-1 autocorrelation (b) for JRA-55 stratospheric temperature at 5 hPa for the period between 1958–2006.

During the winter and spring seasons, the data generally do not have a very strong temporal autocorrelation, and the t test can be applied with a lower risk of obtaining false-positive outcomes. There are some regions where the temporal autocorrelation still persists, such as over North America. Local higher autocorrelation values during other seasons can also be a result due also to low-frequency variance caused by large-scale teleconnections (Madden1977). During the summer season – and to a large extent also in autumn – data are very autocorrelated, but they will be analysed in any case, as it is worthwhile as well to show how the procedure used to assess the possible impact of the geomagnetic activity responds to serially correlated data. In general, autocorrelation is mainly due to the persistence of temperature patterns year by year. For instance, this is the case for example of the large value of temperature autocorrelation found during the summer season. However, we cannot exclude other causes, including a possible impact of the solar activity.

Serial correlation can be corrected for by adopting, for example, the strategy suggested by . This procedure is valid under the assumption that the time series, from which the data are sampled, can be modelled as an autoregressive process of order 1 or AR(1). have shown that the AR(1) representation fits modelled stratospheric temperature data very well according to standard goodness-of-fit tests. found a similar result with temperature observed by radiosondes and satellites.

If EEP has a cumulative impact during the different seasons, it has to be shown that the means of two subsets with high-Ap (H) and low-Ap (L) values from the set N1 must be different.

To test the null hypothesis ${H}_{\mathrm{0}}:{\mathit{\mu }}_{\mathrm{H}}={\mathit{\mu }}_{\mathrm{L}}$ with the t statistics at the 5 % significance level one, let us apply the t test under the condition that the standard deviation is scaled by the equivalent sample sizes me and ne that can be computed by

$\begin{array}{}\text{(1)}& {n}_{\mathrm{e}}=n\left(\frac{\mathrm{1}-{\mathit{\rho }}_{\mathrm{1}}}{\mathrm{1}+{\mathit{\rho }}_{\mathrm{1}}}\right),\end{array}$

where n is the original size of one out of two samples and ρ1 is the parameter of the AR(1) process representing the autocorrelation at lag 1; this is similar for me. The t test is then corrected in the following way:

$\begin{array}{}\text{(2)}& t=\frac{\stackrel{\mathrm{‾}}{H}-\stackrel{\mathrm{‾}}{L}}{s\left(\frac{\mathrm{1}}{\sqrt{{m}_{\mathrm{e}}}}+\frac{\mathrm{1}}{\sqrt{{n}_{\mathrm{e}}}}\right)},\end{array}$

where $\stackrel{\mathrm{‾}}{H}$ and $\stackrel{\mathrm{‾}}{L}$ are the sample averages and s2 is the pooled variance,

$\begin{array}{}\text{(3)}& {s}^{\mathrm{2}}=\frac{\sum _{i=\mathrm{1}}^{m}{\left({H}_{i}-\stackrel{\mathrm{‾}}{H}\right)}^{\mathrm{2}}+\sum _{i=\mathrm{1}}^{n}{\left({L}_{i}-\stackrel{\mathrm{‾}}{L}\right)}^{\mathrm{2}}}{m+n-\mathrm{2}}.\end{array}$

Figure 2Northern Hemisphere seasonal differences in stratospheric temperature (high-Ap–low-Ap values) at 5 hPa without (a) and with (b) temporal (TCC) and spatial (FDR) autocorrelation correction. Grey areas represent statistically significant temperature differences at the 5 % confidence levels.

## 2.4 Accounting for spatial autocorrelation

Spatial autocorrelation produces the so-called multiplicity problem, which arises when testing a statistical hypothesis on many samples (the domain's grid points, in our case) simultaneously. A single hypothesis test allows for a null hypothesis and an alternative hypothesis. The alternative hypothesis will be favoured when an extreme value, usually with a probability (called value) that is less than 5 %, is found (Wilks2016). Making a statistical test on multiple points, for example within a spatial domain, means that more realizations will be available and that there will be many grid points where one is more likely to reject the null hypothesis. In an ideal situation, where the value is set to 0.05 and each point is statistically independent of the others, it is expected to be found that 5 % of the points will be statistically significant by accident. The situation is worse when the grid points are correlated, as is often the case when analysing meteorological and climate data. This problem, known in the literature as the multiplicity problem, has been encountered in several studies, although most of the studies in the atmospheric sciences have not properly addressed the issue yet (Wilks2016). Some solutions have been proposed, each having their own advantages and disadvantages. gives a brief historical outline and shows different solutions to this problem. One technique to address this issue is by using the false discovery rate . According to , the false discovery rate is the expectation of the fraction of true null hypothesis rejections among all the rejections, and it is the best available approach to analyse multiple hypothesis test results, even when those results are mutually correlated.

As stated by , the FDR (false discovery rate) procedure requires smaller values to reject the local null hypothesis, arising the standard of the test. For the sake of the reader, we will describe the FDR algorithm as described in . The algorithm operates on the collection of ${H}_{\mathrm{0}}:{\mathit{\mu }}_{\mathrm{H}}={\mathit{\mu }}_{\mathrm{L}}$ values from me (number of grid points) of local hypothesis tests pi, with i=1, …, N, which are sorted in ascending order. Rejection of the test happens when the pi values are not larger than a threshold level pFDR that is a function of the distribution of the sorted pi values. More specifically, to define which values pass the test, the following formula is used:

$\left[{p}_{i}:{p}_{i}\mathit{⩽}{\mathit{\alpha }}_{\mathrm{FDR}}\left(\frac{i}{N}\right)\right],$

where αFDR is the chosen FDR control level that here is taken to equal 0.05. For a given value of αFDR, the largest value of i such that ${p}_{i}\mathit{⩽}{\mathit{\alpha }}_{\mathrm{FDR}}\left(\frac{i}{N}\right)$ defines the threshold below which the local null hypotheses are rejected.

Figure 3Northern Hemisphere seasonal differences in stratospheric temperature (high-Ap–low-Ap values) at 5 hPa after applying the correction for serial dependence (TCC). Grey areas indicate statistically significant areas at the 5 % confidence level.

Figure 4Northern Hemisphere seasonal differences in stratospheric temperature (high-Ap–low-Ap values) at 5 hPa after applying the FDR correction. Grey areas indicate statistically significant areas at the 5 % confidence level – before FDR correction.

3 Results

## 3.1 Stratospheric levels

We start with the application of the t test on 5 hPa temperature (Fig. 2), which represents the level where the statistically significant area is the largest among all the examined pressure levels. There are large areas with a statistically significant temperature difference at the 5 % level, especially during winter and summer.

At 5 hPa, the area with significant differences covers most of the Northern Hemisphere in JJA, but, as can be seen from the analysis of the Durbin–Watson test, the summer season exhibits a large temporal autocorrelation. Hence, the statistically significant areas observed in JJA should originate from this autocorrelation. In winter, the area with significant points cover North America, another region where the Durbin–Watson test suggests serial correlation. It is clear from Fig. 2 that a possible impact, of the geomagnetic activity, if it exists, would be limited at higher latitudes, from 40 to 90.

Because of the strong temporal autocorrelation, it is expected that at least in summer these significant differences should be false-positive outcomes, and they should be reduced or completely removed when applying the serial correlation correction. In fact, by applying the correction of serial dependence to the 5 hPa temperature differences, the t test results change dramatically, as Fig. 3 shows. The statistically significant differences are removed everywhere in JJA. However, in DJF small areas with significant differences are still present at that level. The Durbin–Watson test somehow predicted that there will be not significant points after applying the Zweirs and von Storch algorithm in the areas where the Durbin–Watson test value was close to zero.

Figure 5As Fig. 2 but for the 1 hPa level.

On the other hand, the problem of multiplicity is solved here by means of the FDR procedure described in Sect. 2. When applying such a procedure without correcting the serial dependence, some significant temperature differences still persist at 5 hPa during summertime (Fig. 4). However, the only application of FDR is to remove all significant differences when it is applied to other pressure levels (e.g. 10 hPa). This result is important, as the FDR procedure is quite powerful in removing most of the false-positive differences, but, as Fig. 2 shows, it is not sufficient in presence of a strong temporal correlation that can still leave regions where the t test rejects the null hypothesis, when, in fact, it would be true. This result is particularly important, and it recommends the application of both the corrections strongly.

The application of such corrections dealing both with temporal and spatial autocorrelation removes all the statistically significant differences in the domain (Fig. 2b), and the t test with the combined correction fails to reject the null hypothesis.

A similar result is obtained for all the other levels in the dataset; temperature differences at 1 and 100 hPa temperatures are shown in Figs. 5 and 6 without and with both the corrections. The application of the false discovery rate on those fields eliminates all the significant temperature differences, showing that also at those levels, there is no detectable impact of geomagnetic activity on the atmospheric temperature.

Figure 6As Fig. 2 but for the 100 hPa level.

Figure 7Zonal mean temperature differences (high-Ap–low-Ap values) without (a) and with (b) temporal (TCC) and spatial (FDR) autocorrelation corrections. The grey areas indicate statistically significant temperature differences at the 5 % confidence level.

## 3.2 Zonally averaged temperature and 2 m temperatures

Several studies have shown the possible impact of EEP or energetic particle precipitation on the observations using zonal mean temperatures . Thus, we show how, without any correction, even the zonal mean temperature difference has areas that are statistically significant at the 5 % level (Fig. 7a). In particular, there are statistically significant areas in all the seasons but spring between 10 and 1 hPa. There are no statically significant areas (Fig. 7b) after applying the two corrections that account for spatial and temporal autocorrelations.

It is natural to think that EEP would influence upper- and mid-stratosphere temperatures through its impact on ozone. The results discussed in the previous sections suggest that the EEP influence on NH stratospheric temperatures is problematic to detect, as it is much weaker than other causes of variability, among which the internal dynamical variability is paramount. As this work is motivated by S09 that analysed the 2 m temperature, Fig. 8a shows the 2 m temperature difference (high-Ap–low-Ap values) without any correction. There are large areas where the seasonal temperature differences are statistically significant at the 5 % level.

Figure 8Temperature differences at 2 m (high-Ap–low-Ap values) without (a) and with (b) temporal (TCC) and spatial (FDR) autocorrelation corrections. The grey areas indicate statistically significant temperature differences at the 5 % level.

The application of both the spatial and temporal autocorrelation corrections remove almost all these areas. However, some small areas of statistically significant temperature differences are still present. They are in the polar region and over Russia during the winter season and over the Scandinavia during the spring (Fig. 8b). As it is not easy to explain these statistically significant surface temperature differences with a causal relationship with EEP, given the lack of signal aloft, there may be some other reasons that can justify this significance with other causes, among which a positive outcome is obtained by chance.

4 Conclusions

Climate data often exhibit temporal and spatial autocorrelations which should be taken into account when testing a hypothesis, a task that is often neglected (Wilks2016). The effect of temporal autocorrelation was addressed with an appropriate procedure described in . The problem of evaluating the results of multiple hypothesis tests in a spatial domain was further addressed by means of the false discovery rate procedure. In this paper, the possible impact of geomagnetic activity on the seasonal-mean stratospheric temperature in the JRA-55 reanalysis was evaluated by means of Welch's t test under four different cases: (1) with no correction of temporal and spatial autocorrelation, (2) with correction on temporal autocorrelation only, (3) with correction on spatial autocorrelation only, and finally (4) with both the corrections. Most of the cases examined show significant points when temporal and spatial autocorrelations are not corrected, while they do not show any significant point when including just one out of the two corrections. In other words, in most cases, there is not even a need to apply both corrections to infer that there is no impact of geomagnetic activity. However, the statistically significant temperature differences at 5 hPa show that it strongly recommended the application of both the corrections for the spatial and temporal autocorrelation. In some cases, like for the JJA temperature difference at 5 hPa, there are a few significant areas remaining when applying one out of the two corrections (Figs. 3 and 4), but those significant areas disappeared when both corrections were applied. Finally, the procedures to take into account these autocorrelations, the significance test typically fails to reject the null hypothesis. This result is found for all the pressure levels analysed and for zonally averaged temperature. The only temperature field that still has statistically significant differences after applying both the corrections is the 2 m temperature. There are two seasons, DJF and MAM, where small statistically significant areas are present in the polar region. In the absence of a signature aloft, we therefore conclude that, based on the JRA-55 reanalyses, not enough evidence is available at present to suggest that the null hypothesis of no impact of geomagnetic activity on NH stratospheric temperatures is false. A remaining caveat concerns the definition of seasons of high or low geomagnetic activity, which is here the same as in S09 and is based on a lagged 4-month-averaged Ap index (i.e. from October to January for wintertime geomagnetic activity). Some sensitivity studies to this definition, e.g. to treat more intense shorter episodes of EEP or to treat differently the seasonal lag or accumulation of EEP, are certainly warranted for future studies. It is clear that the absence or the presence of significance does not put an end to the research of a possible relationship between EEP and stratospheric temperature, which we suppose to be weak and consequently difficult to detect.

Data availability
Data availability.

Data can be downloaded from the Meteorological Research Institute, Japan Meteorological Agency, Japan, or from the Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory.

Author contributions
Author contributions.

NT performed the statistical tests; all the authors contributed to the interpretation of the data and wrote the paper.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

This work was funded by the SOLENA (Solar effects on natural climate variability in the North Atlantic and Arctic) project, The Research Council of Norway (grant no. 255276/E10). The authors acknowledge the Meteorological Research Institute, Japan Meteorological Agency, Japan, for the JRA-55 reanalysis. We are extremely grateful to two reviewers for their valuable comments, corrections and suggestions.

Financial support
Financial support.

This research has been supported by the Program for Space Research Project, The Research Council of Norway (grant no. 255276/E10).

Review statement
Review statement.

This paper was edited by Petr Pisoft and reviewed by two anonymous referees.

References

Andersson, M. E., Verronen, P. T., Rodger, C. J., Clilverd, M. A., and Seppälä, A.: Missing driver in the Sun Earth connection from energetic electron precipitation impacts mesospheric ozone, Nat. Commun., 5, 5197, https://doi.org/10.1038/ncomms6197, 2014. a

Baumgaertner, A. J. G., Jöckel, P., Dameris, M., and Crutzen, P. J.: Will climate change increase ozone depletion from low-energy-electron precipitation?, Atmos. Chem. Phys., 10, 9647–9656, https://doi.org/10.5194/acp-10-9647-2010, 2010. a

Benjamini, Y., and Hochberg,Y.: Controlling the false discovery rate: A practical and powerful approach to multiple testing, J. Roy. Stat. Soc. B, 57, 289–300, https://doi.org/10.1111/j.2517-6161.1995.tb02031.x, 1995. a, b

Brasseur, G. P. and Solomon, S.: Aeronomy of the Middle Atmosphere, 2nd revised Edn., D. Reidel Publishing Company, the Netherlands, 1986. a

Bucha, V.: Geomagnetic activity and North Atlantic Oscillation, Stud. Geophys. Geod., 58, 461–472, https://doi.org/10.1007/s11200-014-0508-z, 2014. a

Callis, L. B., Natarajan, M., Lambeth, J. D., and Baker, D. N.: Solar-atmospheric coupling by electrons (SOLACE). 2. Calculated stratospheric effects of precipitating electrons, 1979–1988, J. Geophys. Res., 103, 28421–28438, https://doi.org/10.1029/98JD02407, 1998. a

de Winter, J. C. F.: Using the Student's t-test with extremely small sample sizes, Pract. Assess. Res. Eval., 18, 10, https://doi.org/10.7275/e4r6-dj05, 2013. a

Durbin, J., and Watson, G. S.: Testing for Serial Correlation in Least Squares Regression, I., Biometrika, 37, 409–428, https://doi.org/10.1093/biomet/37.3-4.409, 1950. a, b

Efron, B.: Student's t-Test Under Symmetry Conditions, J. Am. Stat. Assoc., 64, 1278–1302, https://doi.org/10.1080/01621459.1969.10501056, 1969. a

Funke, B., López-Puertas, M., Gil-López, S., von Clarmann, T., Stiller, G. P., Fischer, H., and Kellmann, S.: Downward transport of upper atmospheric NOx into the polar strato sphere and lower mesosphere during the Antarctic 2003 and Arctic 2002/2003 winters, J. Geophys. Res., 110, D24308, https://doi.org/10.1029/2005JD006463, 2005. a, b

Funke, B., Puertas, M. L., Holt, L., Randall, C. E., Stiller, G. P., and von Clarmann, T.: Hemispheric distributions and interannual variability of NOy produced by energetic particle precipitation in 2002–2012, J. Geophys. Res., 119, 13565–13582, https://doi.org/10.1002/2014JD022423, 2014. a

Jacob, D. J.: Introduction to Atmospheric Chemistry, Princeton University Press, Princeton, New Jersey, 1999. a

Kobayashi, S., Ota, Y., Harada, Y., Ebita, A., Moriya, M., Onoda, H., Onogi, K., Kamahori, H., Kobayashi, C., Endo, H., Miyaoka, K., and Takahashi, K.: The JRA-55 Reanalysis: General specifications and basic characteristics, J. Meteorol. Soc. Jpn., 93, 5–48, https://doi.org/10.2151/jmsj.2015-001, 2015. a

Long, C. S., Fujiwara, M., Davis, S., Mitchell, D. M., and Wright, C. J.: Climatology and interannual variability of dynamic variables in multiple reanalyses evaluated by the SPARC Reanalysis Intercomparison Project (S-RIP), Atmos. Chem. Phys., 17, 14593–14629, https://doi.org/10.5194/acp-17-14593-2017, 2017. a

Lu, H., Jarvis, M. J., and Hibbins, R. E.: Possible solar wind effect on the northern annular mode and northern hemispheric circulation during winter and spring, J. Geophys. Res., 113, D23104, https://doi.org/10.1029/2008JD010848, 2008. a

Madden, R. A.: Estimates of the Autocorrelations and Spectra of Seasonal Mean Temperatures over North America, Mon. Weather Rev., 105, 9–18, https://doi.org/10.1175/1520-0493(1977)105<0009:EOTAAS>2.0.CO;2, 1977. a

Maliniemi, V., Asikainen, T., and Mursula, K.: Spatial distribution of Northern Hemisphere winter temperatures during different phases of the solar cycle, J. Geophys. Res.-Atmos., 119, 9752–9764, https://doi.org/10.1002/2013JD021343, 2014. a

Päivärinta, S.-M., Verronen, P. T., Funke, B., Gardini, A., Seppälä, A., and Andersson, M. E.: Transport versus energetic particle precipitation: Northern polar stratospheric NOx and ozone in January–March 2012, J. Geophys. Res.-Atmos., 121, 6085–6100, https://doi.org/10.1002/2015JD024217, 2016. a, b

Poncet, A., Courvoisier, D. S., Combescure, C., and Perneger, T. V.: Normality and sample size do not matter for the selection of an appropriate statistical test for two-group comparisons, Methodology, 12, 61–71, https://doi.org/10.1027/1614-2241/a000110, 2016. a

Randall, C. E., Harvey, V. L., Manney, G. L., Orsolini, Y., Codrescu, M., Sioris, C., Brohede, S., Haley, C. S., Gordley, L. L., Zawodny, J. M., and Russell, J. M.: Stratospheric effects of energetic particle precipitation in 2003–2004, Geophys. Res. Lett., 32, L05802, https://doi.org/10.1029/2004GL022003, 2005. a

Randall, C. E., Harvey, V. L., Singleton, C. S., Bernath, P. F., Boone, C. D., and Kozyra, J. U.: Enhanced NOx in 2006 linked to strong upper stratospheric Arctic vortex, Geophys. Res. Lett., 33, L18811, https://doi.org/10.1029/2006GL027160, 2006. a

Rostoker, G.: Geomagnetic indices, Rev. Geophys. Space Phys., 10, 935–950, https://doi.org/10.1029/RG010i004p00935, 1972. a

Seidel, D. J. and Lanzante, J. R.: An assessment of three alternatives to linear trends for characterizing global atmospheric temperature changes, J. Geophys. Res., 109, D14108, https://doi.org/10.1029/2003JD004414, 2004. a

Seppälä, A., Randall, C. E., Clilverd, M. A., Rozanov, E., and Rodger, C. J.: Geomagnetic activity and polar surface air temperature variability, J. Geophys. Res., 114, A10312, https://doi.org/10.1029/2008JA014029, 2009. a, b, c

Seppälä, A., Lu, H., Clilverd, M. A., and Rodger, C. J.: Geomagnetic activity signatures in wintertime stratosphere wind, temperature, and wave response, J. Geophys. Res.-Atmos., 118, 2169–2183, https://doi.org/10.1002/jgrd.50236, 2013. a, b

Sinnhuber, M., Nieder, H., and Wieters, N.: Energetic Particle Precipitation and the Chemistry of the Mesosphere/Lower Thermosphere, Surv. Geophys., 33, 1281–1334, https://doi.org/10.1007/s10712-012-9201-3, 2012.  a

Tomikawa, Y.: Response of the Middle Atmosphere in the Southern Hemisphere to Energetic Particle Precipitation in the Latest Reanalysis Data, SOLA, 13A, 1–7, https://doi.org/10.2151/sola.13A-001, 2017. a, b, c

Vyushin, D. I., Kushner, P. J., and Zwiers, F.: Modeling and understanding persistence of climate variability, J. Geophys. Res., 117, D21106, https://doi.org/10.1029/2012JD018240, 2012. a

Wilks, D. S.: Statistical Methods in the Atmospheric Sciences, Elsevier Science, Burlington, MA, 2006. a, b

Wilks, D. S.: The stippling shows statistically significant grid-points. How research results are routinely overstated and overinterpreted, and what to do about it, B. Am. Meteorol. Soc., 97, 2263–2273, https://doi.org/10.1175/bams-d-15-00267.1, 2016. a, b, c, d, e, f, g, h, i

Zwiers, F. W. and von Storch, H.: Taking Serial Correlation into Account in Tests of the Mean, J. Climate, 8, 336–351, https://doi.org/10.1175/1520-0442(1995)008<0336:TSCIAI>2.0.CO;2, 1995. a, b, c, d