Connection between Antarctic Ozone and Climate: Interannual Precipitation Changes in the Southern Hemisphere

: In this study, we explored the connection between anomalies in springtime Antarctic ozone and all-year precipitation in the Southern Hemisphere by using observations from 1960–2018 and coupled simulations for 1960–2050. The observations showed that this correlation was enhanced during the last several decades, when a simultaneously increased coupling between ozone and Southern Annular Mode (SAM) anomalies became broader, covering most of the following summer and part of the previous winter. For eastern Australia, the ozone–precipitation connection shows a greater persistence toward the following summer than for other regions. On the other hand, for South America, the ozone–precipitation correlation seems more robust, especially in the early summer. There, the correlation also covers part of the previous winter, suggesting that winter planetary waves could a ﬀ ect both parameters. Further, we estimated the sensitivity of precipitation to changes in Antarctic ozone. In both observations and simulations, we found comparable sensitivity values during the spring–summer period. Overall, our results indicate that ozone anomalies can be understood as a tracer of stratospheric circulation. However, simulations indicate that stratospheric ozone chemistry still contributes to strengthening the interannual relationship between ozone and surface climate. Because simulations reproduced most of the observed connections, we suggest that including ozone variability in seasonal forecasting systems can potentially improve predictions.


Introduction
Although mean global precipitation is likely to increase in subsequent decades, as a response to global warming, the geographical patterns of this change are expected to be diverse [1]. In addition to anthropogenic greenhouse gases (GHG), the Antarctic ozone (O 3 ) hole [2] is also thought to have driven the evolution of surface climate in the Southern Hemisphere (SH) during the last several decades [3,4]. Indeed, the ozone hole contributes to the temperature decrease in the Antarctic lower stratosphere [5] and, although ozone depletion occurs during the austral spring, the associated reduction in shortwave heating lasts long into the summer, accelerates the stratospheric zonal winds, and causes negative this study was exploring the possible use of springtime Antarctic ozone as a predictor of summertime precipitation, we further analyzed this connection for the other seasons.
The advantage of using ozone observations rather than other atmospheric dynamics proxies from reanalysis [18] is the longer extension of the ozone series compared to reanalysis products, which are mostly constrained by the duration of satellite observations. This longer time series allowed us to explore the relationship between ozone and surface climate even before the development of the ozone hole. On the other hand, the effect of ozone on tropospheric weather is thought to occur largely through stratospheric dynamics changes. Therefore, the use of stratospheric dynamics is valuable if the focus is on the impact of the stratospheric changes on tropospheric dynamics. From this perspective, the ozone time series at only one location would present some limitations if compared with the 3D fields in meteorological analyses.

Precipitation
Changes in precipitation were investigated with the Global Precipitation Climatology Centre (GPCC) Full Data Reanalysis V7 monthly precipitation dataset [31], which is based on quality-controlled data from 67,200 stations. To generate this dataset, precipitation anomalies recorded at the stations were interpolated and then superimposed on the GPCC climatology [31]. The spatial resolution of the selected product is 1 • × 1 • latitude by longitude. The data were analyzed for the period 1960-2018.

Total Column Ozone
Reliable and continuous satellite observations of the total column ozone (hereinafter total ozone) over Antarctica have been conducted since about 1979. Therefore, we exploited ground-based total ozone observations recorded in Antarctica to estimate the evolution of the ozone hole. Total ozone observations in Antarctica roughly started with The Third International Polar Year (1957Year ( -1958  Previous studies showed that shifts in the location of the polar vortex could bias observations of air masses recorded at stations near the vortex edge [32]. Therefore, according to their latitude, Faraday and Syowa stations are more affected by this shift than Halley and the South Pole. Although Halley is usually inside the polar vortex, here we focused on the South Pole because, potentially, it is the station more within the deep vortex core. In particular, we exploited total ozone observations recorded at the South Pole Station in November. Nevertheless, since the November time series of the South Pole and Halley are well correlated (e.g., R 2 = 0.89), the main results of this study do not substantially change if Halley instead of South Pole data are used. Besides the fact that the South Pole time series has almost no data gaps in November, we selected this month for three additional reasons: (i) the ozone radiative forcing peaks in the lower stratosphere in November (e.g., [24] and reference therein); (ii) ozone anomalies in November are much larger and easier identifiable than anomalies in September or October; (iii) during the investigated period, November ozone anomalies are well correlated with the SAM index.
Moreover, the total ozone dataset from the multi-sensor reanalysis version 2 (MSR2, [33]) has been used to explore the dependence of the connection between ozone and surface climate with the location of the ozone series.

Simulations
Moreover, when possible, observations were compared with a mean ensemble of transient simulations from a coupled climate model with interactive atmospheric chemistry and ocean and sea-ice  [34]). It has a well-resolved stratosphere and an internally generated quasi-biennial oscillation. The upper boundary is in the middle mesosphere (at 0.01 hPa) and the spatial resolution is approximately 2.8 • in latitude and longitude. Within the Coupled Model Intercomparison Project 5 (CMIP5) framework, Eyring et al. [12] evaluated MIROC ozone climatology, long-term changes, and their associated climate impacts against other models and observations and found a limited Antarctic ozone overestimation in spring. However, due to the inclusion of the interactive chemistry, the ozone response of the climate system in MIROC is fully consistent with the ozone forcing (i.e., the temporal evolution of the dynamics is coupled with the ozone-climate feedback). In contrast, the responses obtained from most of the CMIP5 models did not account for chemistry-climate feedbacks.
The MIROC simulations forced with all external drivers and labeled historical in CMIP5, covered the period 1960-2005, whereas the future projections following the Representative Concentration Pathway (RCP) 4.5, covered the period 2006-2050. Moreover, we ran additional MIROC simulations covering the same periods as above (i.e., 1960-2005 and 2006-2050) but forced with the same forcing values, except for the ozone-depleting substance (ODS) concentrations, which were kept at 1960 levels (i.e., ozone depletion could not develop during the following decades). An ensemble of nine transient historical simulations was used to compare with the observations, whereas the ensembles of fixed ODS and RCP 4.5 projections were based on three simulations.
Although ten of the CMIP5 models simulate atmospheric chemistry interactively, no fixed ODS simulations are freely available. We also revised the datasets of the Chemistry Climate Model Initiative (CCMI) [35] looking for comparable simulations performed with models including interactive atmospheric chemistry, fully interactive ocean, and multiple members available. However, as their future projections are based on RCP 6.0, we discarded their use in this study.
SAM indices from Marshall [36] were employed for the observational analysis, and the difference between the zonal mean sea level pressure at 40 • S and 65 • S [37] was used for the analysis of the simulations. On the other hand, the Niño 3.4 index was used for both datasets.
Observational and model-based datasets were investigated by correlation and composite analysis, and statistical significance was assessed by a two-sided Student's t-test. For simulations, we analyzed each member independently and then averaged the results. Although this procedure tends to smooth the magnitude of the final product, the results should be considered more robust. Figure 1a shows the temporal evolution of the total ozone at the South Pole in November as measured by ground-based observations (for 1960-2018) and reproduced by MIROC all drivers and MIROC fixed ODS simulations (for 1960-2050). The ozone depletion and the expected on-going recovery appear to be reproduced satisfactorily. Additionally, since about halfway through the 1970s, the ozone levels in MIROC all drivers are always well below the amounts of the fixed ODS simulations.

Results
Aside from the evident ODS-induced decreasing trend in ozone levels, Figure 1a also highlights large year-to-year variability. Therefore, Figure 1b was constructed to present the anomalies (i.e., residuals) calculated with respect to a third-degree polynomial fitted to the time series of ozone values of Figure 1a (note that the following results remain virtually the same if a different polynomial is employed). In general, the simulated anomalies are somewhat smaller than the observed anomalies (see the respective standard deviations in the upper corners). We will use these anomalies to build the correlations and composites presented later in the paper.
The bottom panels of Figure 1 show the correlation maps between the average precipitation anomalies in central-southern Chile (red box highlighted by the red arrow) and the precipitation anomalies in the given grid point for GPCC and MIROC (all drivers) in the austral spring (October-November-December, OND) for 1960-2018 and 1960-2005, respectively. As our objective was assessing the ozone variations in November as a potential predictor of precipitation in summer, here and Atmosphere 2020, 11, 579 5 of 18 in the rest of this paper, OND is assumed to be the reference three-month period. By removing the trends, we also removed at least part of the potential influence of GHGs on precipitation. Therefore, we found it appropriate to compare a somewhat longer observational period (Figure 1d) with the simulations (Figure 1c) to achieve a more robust statistic (however, the standard deviation of the 59-year long series is identical to that of the 46-year long series; see Figure 1b; the observed correlation map for 1960-2005 is virtually equivalent to Figure 1d and shown in Figure A1 of the Appendix A). Overall, the model does a good job of reproducing the main observed pattern, which shows an obvious positive correlation in southern South America but also a negative correlation in both Australia and South Africa. This pattern confirms previous findings [13] and points to the dependence of SH precipitation on large scale modes, specifically the SAM.  Figure A1 of the Appendix A). Overall, the model does a good job of reproducing the main observed pattern, which shows an obvious positive correlation in southern South America but also a negative correlation in both Australia and South Africa. This pattern confirms previous findings [13] and points to the dependence of SH precipitation on large scale modes, specifically the SAM.  Because the decreasing springtime ozone trend has been associated with the positive drift of the SAM in summer, to highlight the temporal evolution of their connection, Figure 2a shows the 31-year running correlation (with correlation values assigned to the central year of each 31-year-long period) between November total ozone (as in Figure 1a) and the SAM index in the given month (each computed as a mean of three months). To shed additional light on this connection, we also extended the analysis to the previous winter. In this way, we explored the opposite situation, i.e., the SAM in winter leads the ozone in November, likely via changes in winter planetary waves affecting both of them [26]. Overall, a negative correlation, peaking in spring and extending from the next summer to the previous winter, appears to be stronger during the last decades, characterized by the development of the ozone hole, than during the period before. Because the decreasing springtime ozone trend has been associated with the positive drift of the SAM in summer, to highlight the temporal evolution of their connection, Figure 2a shows the 31-year running correlation (with correlation values assigned to the central year of each 31-year-long period) between November total ozone (as in Figure 1a) and the SAM index in the given month (each computed as a mean of three months). To shed additional light on this connection, we also extended the analysis to the previous winter. In this way, we explored the opposite situation, i.e., the SAM in winter leads the ozone in November, likely via changes in winter planetary waves affecting both of them [26]. Overall, a negative correlation, peaking in spring and extending from the next summer to the previous winter, appears to be stronger during the last decades, characterized by the development of the ozone hole, than during the period before.   Here, the discontinuity between the last decades (e.g., 31-year running correlation centered around 1987-2003) and the previous period (e.g., around 1976-1987) is even more evident, and the most recent correlations (e.g., 1998-2003) extend well into the next summer. Figure 2c shows the running mean of the composite difference in the SAM index achieved by stratifying the data according to the ozone levels in November (i.e., low minus high O 3 years). Moreover, the running mean of the spread between low and high ozone anomalies (hereafter, ozone change) is shown as a percentage with respect to the mean on the opposite axis. The resulting SAM difference is positive during the entire time range investigated, but it peaks during the last decades in correspondence with larger ozone changes (red line).
Previous studies have highlighted the influence of the ozone hole on climate by using ground-based ozone observations recorded at various Antarctic stations [5,24,26,29], or ozone data averaged over the polar cap [25], as a proxy of the ozone hole. To evaluate the dependence of the connection between Antarctic ozone and precipitation with the location of the ozone series employed, Figure 2d-f shows the spatial correlation between anomalies of MSR2 total ozone in November for a given pixel and anomalies of GPCC precipitation in NDJ averaged over Australia, South Africa, and Chile (red boxes). We focused on the period of 1988-2018 when the correlation between SAM and O 3 anomaly was at the highest levels ( Figure 2b); nevertheless, the results are virtually equivalent to the case of using the full MSR2 dataset (i.e., 1979-2018; not shown). For Australia, the correlation is negative and statistically significant at the 95% level only within a limited Antarctic region (i.e., southern of 70 • S and roughly within the sector 40 • W-40 • E), but it is not even significant for South Africa. In contrast, a clear positive correlation extending over the whole of Antarctica is evident for Chile. Moreover, due to the influence of the Brewer-Dobson circulation on the ozone distribution, all regions show an evident correlation with ozone at tropical latitudes. The simulations (Figure 2g-i) reproduce both the spatial patterns and the sign of the observed correlations and also show a consistent correlation over the Indian Ocean that is not present in the observations. Overall, Figure 2d-i show that the precipitation averaged over the three investigated regions is always correlated with ozone over most of Antarctica. Nevertheless, it can be noted that, usually, ozone pixels at higher latitudes present greater correlation than pixels in other locations for both observations and simulations. This fact adds confidence in using the total ozone observations recorded at the South Pole station against other possible ground-based datasets when a single dataset has to be used as a proxy of the precipitation variability for the whole southern hemisphere. On the other hand, the availably of 3D ozone fields would allow using a proxy based on Antarctic zonal means averaged over high latitudes (e.g., greater than 65-70 • S) which is expected to be more representative of the ozone hole variability. Therefore, future studies focusing on precipitation prediction capabilities are likely to be further strengthened by incorporating observed or simulated 3D ozone fields.
To provide evidence of the possible use of the O 3 variability as a predictor of precipitation, the left column of Figure 3 shows the correlation maps between O 3 anomalies in November at the South Pole and the observed GPCC precipitation anomalies in (top to bottom) January-February-March (JFM), December-January-February (DJF), November-December-January (NDJ), and OND for 1988-2018 (i.e., during the period characterized by the highest correlation between SAM and O 3 anomalies).
Overall, the positive correlations along the Pacific coastline of South America and the negative correlations in South Africa and eastern Australia are evident. Over South America, correlations are stronger in spring/early summer (i.e., OND-NDJ), but they tend to persist longer over Australia. In the right column, we show the patterns simulated by MIROC for 1960-2005. Despite the expected influence of El Niño-Southern Oscillation (ENSO) on precipitation, MIROC confirms the main features of the left panels with a somewhat stronger persistence in summer.
Atmosphere 2020, 11, x FOR PEER REVIEW 8 of 19 influence of El Niño-Southern Oscillation (ENSO) on precipitation, MIROC confirms the main features of the left panels with a somewhat stronger persistence in summer.  Figure 4 shows the 31-year running correlations between O3 anomalies in November at the South Pole (top panels), as well as the SAM (central panels) and ENSO (bottom panels) anomalies in their given months, and the observed GPCC precipitation anomalies averaged over the three regions highlighted above. As a reference, the red vertical dotted lines (top panels) highlight the ozoneprecipitation correlation for the most recent 31-year-long period shown in Figure 3a-d. For Chile (left column), the correlation between the SAM and the precipitation is negative, especially during the last decades. The correlation attains its maximum in the late spring/early summer, is reduced in late summer, and then increases again during autumn. Except for autumn, a similar pattern, but with the opposite sign, emerges in the correlation between O3 and precipitation (Figure 4a), while a significant (positive) correlation between ENSO and precipitation (Figure 4g) is evident only in early spring.  Figure 4 shows the 31-year running correlations between O 3 anomalies in November at the South Pole (top panels), as well as the SAM (central panels) and ENSO (bottom panels) anomalies in their given months, and the observed GPCC precipitation anomalies averaged over the three regions highlighted above. As a reference, the red vertical dotted lines (top panels) highlight the ozone-precipitation correlation for the most recent 31-year-long period shown in Figure 3a-d. For Chile (left column), the correlation between the SAM and the precipitation is negative, especially during the last decades. The correlation attains its maximum in the late spring/early summer, is reduced in late summer, and then increases again during autumn. Except for autumn, a similar pattern, but with the opposite sign, emerges in the correlation between O 3 and precipitation (Figure 4a), while a significant (positive) correlation between ENSO and precipitation (Figure 4g) is evident only in early spring.
The correlations of O 3 and SAM with precipitation also have the opposite sign for Australia (central column), although they show a different temporal evolution. Indeed, the relationship between SAM (O 3 ) and precipitation maximizes in summer, mostly during the pre-ozone hole (ozone hole) period. Nevertheless, the O 3 -precipitation correlation persists longer toward summer than in the other regions.
The ozone variability in November also appears to be a potential predictor of precipitation for South Africa. However, it is worth noting that the impact of ENSO on precipitation maximizes in late spring/summer for both Australia and South Africa. It makes the interpretation of the results and their attribution more ambiguous for these regions.  Figure 3 (note the inverted color scale) for both observations and simulations. Especially in the observations, positive differences in the late spring are larger than those in summer and maximize around the Pacific coastline of South America and eastern Australia (note that the negative change over western Australia in late summer is not associated with any significant correlation in Figure 3a-b). Although statistically significant changes in MIROC are more widespread, their magnitude is usually smaller due to averaging.  The correlations of O3 and SAM with precipitation also have the opposite sign for Australia (central column), although they show a different temporal evolution. Indeed, the relationship between SAM (O3) and precipitation maximizes in summer, mostly during the pre-ozone hole (ozone hole) period. Nevertheless, the O3-precipitation correlation persists longer toward summer than in the other regions.
The ozone variability in November also appears to be a potential predictor of precipitation for South Africa. However, it is worth noting that the impact of ENSO on precipitation maximizes in late  are larger than those in summer and maximize around the Pacific coastline of South America and eastern Australia (note that the negative change over western Australia in late summer is not associated with any significant correlation in Figure 3a-b). Although statistically significant changes in MIROC are more widespread, their magnitude is usually smaller due to averaging. An issue in adopting composite difference analysis is the choice of a correct threshold, which should result in an optimal representation of the investigated phenomenon. For example, in Figure  5, we used a somewhat conservative ozone threshold (i.e., σ/2), so the magnitude and pattern could potentially be different with higher thresholds.
To assess this aspect, Figure 6 shows the composites (low minus high O3), similar to Figure 5, but for distinct seasons and computed by stratifying the time series of the SAM, Niño 3.4, and mean precipitation (averaged within the three investigated regions) as a function of the ozone change (i.e., the mean spread between high and low O3 anomalies resulting from an increasing O3 threshold with steps of 1 DU; the lowest value shows differences for the zero threshold, and the highest value shows differences for the threshold resultant from five samples in both high and low anomalies). As a reference, the red dotted lines highlight the composite difference for σ/2 shown in Figure 5. Overall, the composite difference of the observed SAM is positive and depicts a general increasing difference in SAM along with increasing ozone change (Figure 6b). Such differences have the same sign as the trends of the last decades and roughly span from the previous winter to the following summer. An issue in adopting composite difference analysis is the choice of a correct threshold, which should result in an optimal representation of the investigated phenomenon. For example, in Figure 5, we used a somewhat conservative ozone threshold (i.e., σ/2), so the magnitude and pattern could potentially be different with higher thresholds.
To assess this aspect, Figure 6 shows the composites (low minus high O 3 ), similar to Figure 5, but for distinct seasons and computed by stratifying the time series of the SAM, Niño 3.4, and mean precipitation (averaged within the three investigated regions) as a function of the ozone change (i.e., the mean spread between high and low O 3 anomalies resulting from an increasing O 3 threshold with steps of 1 DU; the lowest value shows differences for the zero threshold, and the highest value shows differences for the threshold resultant from five samples in both high and low anomalies). As a reference, the red dotted lines highlight the composite difference for σ/2 shown in Figure 5. Overall, the composite difference of the observed SAM is positive and depicts a general increasing difference in SAM along with increasing ozone change (Figure 6b). Such differences have the same sign as the trends of the last decades and roughly span from the previous winter to the following summer. Figure 6c-e shows the difference in precipitation within the various regions. While we observe a tendency toward increasing drying along with increasing ozone change for Chile (Figure 6c), Australia, and South Africa show the opposite trend (Figure 6d,e; however, for South Africa, the sign reverts for high ozone). For Chile and Australia, and for South Africa to a lesser extent, this difference in precipitation, which develops under a positive SAM, is characterized by a well-defined pattern that extends to the following summer. Nonetheless, the composites show some discontinuities. Precipitation differences associated with smaller ozone changes (i.e., O 3 < 40%) usually extend toward the following summer and match the simulations better (Figure 6h-j). In contrast, for Australia, the observed composite difference associated with larger changes also covers the previous winter. Finally, differences in precipitation associated with the highest ozone changes (i.e., O 3 > 50%) are coupled with a moderate positive difference in the Niño 3.4 index (Figure 6a), which is already developed in the previous winter and peaks in the following summer. For South Africa, this discontinuity reverses the sign of the precipitation difference; however, for Australia, it is associated with a sudden decrease in these differences. Further composite analysis of other modes of variability did not present any significant change (e.g., Figure A2 shows the composite differences for the Indian Ocean Dipole index [38] which is supposedly affecting Australian precipitation). Note that correlations between ENSO and precipitation for both regions are negative in summer (Figure 4h,i), and are likely to explain the change of pattern. Although the correlation between ENSO and precipitation is mostly low and not significant for Chile, ENSO tends to reinforce the already negative difference.   (Figure 6d,e; however, for South Africa, the sign reverts for high ozone). For Chile and Australia, and for South Africa to a lesser extent, this difference in precipitation, which develops under a positive SAM, is characterized by a well-defined pattern that extends to the following summer. Nonetheless, the composites show some discontinuities. Overall, MIROC (all drivers) historical simulations (Figure 6f-j) reproduced the observed patterns but with some differences. First, different from the observations, in MIROC, the peak in the SAM composite (Figure 6g), which occurs in NDJ, is mostly associated with the largest precipitation difference. In addition, in agreement with previous studies [26], simulated SAM differences persist well into the summer, more so than in the observations. Furthermore, although the largest change in the SAM occurs one month before the reference period (which is in OND) in the observations (Figure 6b), in the simulations (Figure 6g,l,q,v), it occurs one month later (i.e., SON for observations and NDJ for simulations). Finally, the observations show an increased positive SAM in winter, which is largely absent in the historical MIROC (all drivers) simulations.
MIROC (all drivers) does not show any evident difference in the ENSO composite (Figure 6f). Moreover, statistically significant changes in precipitation do not occur in winter; however, under comparable ozone change, the observed intensity of the precipitation difference is somewhat well reproduced for spring-summer, although MIROC precipitation peaks a month later (NDJ) than in the GPCC data.
In general, the future projections (Figures 6k-o and 6u-y, the latter with ODS fixed at 1960 levels) present less significant precipitation changes, and this is particularly evident for Australia (Figure 6n). Moreover, SAM differences extend more toward the previous winter (Figure 6l,v) than in the case of the respective historical simulations (Figure 6g,q). Nevertheless, although they are small, we note statistically significant negative differences in the ENSO composites (Figure 6k,u). Finally, simulations with ODS fixed at 1960 levels (Figure 6p-y) generally show composite differences characterized by patterns analogous to those of the corresponding historical simulations, although with lower significance (e.g., Figure 6q). Figure 7 shows the slope of the linear regression line and the Pearson correlation between the precipitation differences and the corresponding ozone changes shown in Figure 6. Here, the slope depicts the sensitivity of the composite difference in precipitation to the ozone change (i.e., ∆P/∆O 3 ). Because precipitation differences and ozone changes are usually scarcely correlated in winter (Figure 7d-f) and values are not statistically significant (except for observations in Australia in ASO; see Figure 7b), we focus the following discussion on the spring-summer period, although we also show the values for winter, for completeness. For Chile, the simulated (red lines) and observed (black lines) correlations are negative and peak at approximately −0.9 in both NDJ and DJF, with somewhat smaller values in JFM (Figure 7d). In NDJ, ∆P/∆O 3 is around −0.7 for both observations (black lines) and simulations (red lines; see Figure 7a); however, in DJF and JFM, the sensitivity of the observations becomes slightly greater than the sensitivity of the simulations. Nevertheless, limiting the observational dataset to ozone changes comparable to the simulations (e.g., O 3 < 40%, which corresponds to the discontinuity discussed above; see blue lines) makes its sensitivity much closer to the simulations in summer. The same occurs for Australia, with both the observed and simulated sensitivities peaking at around 0.4 in DJF. Here, it is worth noting that sensitivity change persists longer into the summer (FMA) for the simulations than for the observations (see [29]). Finally, no definitive conclusions can be drawn for South Africa, where the composite differences appear to be less robust than for the other regions. Indeed, Figure 7c

Discussion and Conclusions
In this study, we analyzed the interannual relationship between Antarctic ozone and precipitation in the SH with special care given to some regions where previous studies showed an ozone hole-induced influence on the trends of the surface climate. Although the main goal was evaluating late springtime (i.e., November) ozone anomalies as a predictor of summer precipitation, we also extended the investigation to other seasons.
First, we found that the correlation between O3 anomalies in November and SAM anomalies in late spring/summer was higher (lower) and more (less) persistent during the most recent (oldest) decades (Figure 2b). The increased correlation between O3 and SAM during the latter decades ( Figure  2a) has been ascribed to an overall ozone-induced strengthening of the vortex along with a consequent postponed final warming and deeper and long-lasting stratosphere-troposphere coupling, ultimately resulting in a more frequent positive phase of the SAM [39]. This mechanism also potentially worked at the interannual time scale.
However, the 31 years running correlation between November ozone and spring-summer SAM anomalies for MIROC (all drivers) simulations showed a robust correlation also in the earlier decades, e.g., during 1965-1995 (see Figure A3a; correlations values were assigned to the central years, i.e., 1980 for 1965-1995). Although observations showed less evident changes in this period (Figure 2b), we still note a corresponding peak in spring (SON) with a tendency for the signal to extend into the summer. Furthermore, the correlations between the observed ENSO and SAM anomalies ( Figure  A4a) indicate that they are not fully independent and that the sign of the correlation changes depending on the season [40]. It has been suggested that this connection can occur because larger sea

Discussion and Conclusions
In this study, we analyzed the interannual relationship between Antarctic ozone and precipitation in the SH with special care given to some regions where previous studies showed an ozone hole-induced influence on the trends of the surface climate. Although the main goal was evaluating late springtime (i.e., November) ozone anomalies as a predictor of summer precipitation, we also extended the investigation to other seasons.
First, we found that the correlation between O 3 anomalies in November and SAM anomalies in late spring/summer was higher (lower) and more (less) persistent during the most recent (oldest) decades (Figure 2b). The increased correlation between O 3 and SAM during the latter decades (Figure 2a) has been ascribed to an overall ozone-induced strengthening of the vortex along with a consequent postponed final warming and deeper and long-lasting stratosphere-troposphere coupling, ultimately resulting in a more frequent positive phase of the SAM [39]. This mechanism also potentially worked at the interannual time scale.
However, the 31 years running correlation between November ozone and spring-summer SAM anomalies for MIROC (all drivers) simulations showed a robust correlation also in the earlier decades, e.g., during 1965-1995 (see Figure A3a; correlations values were assigned to the central years, i.e., 1980 for 1965-1995). Although observations showed less evident changes in this period (Figure 2b), we still note a corresponding peak in spring (SON) with a tendency for the signal to extend into the summer. Furthermore, the correlations between the observed ENSO and SAM anomalies ( Figure A4a) indicate that they are not fully independent and that the sign of the correlation changes depending on the season [40]. It has been suggested that this connection can occur because larger sea surface temperature anomalies in the tropical Pacific during El Niño years force Rossby wavetrains, resulting in positive pressure anomalies at higher latitudes (i.e., Amundsen-Bellingshausen Sea) [41] and, therefore, in negative correlation. Since about 1999, the spatial structure of ENSO has moved from an eastern Pacific to a central Pacific type [42]. It has been suggested that this variation is also coupled with a change in the SAM [43]. Finally, although we did not find an evident correlation between November O 3 and ENSO anomalies ( Figure A4b), previous studies have suggested that more frequent vortex disturbances and positive ozone anomalies are associated with El Niño [44].
Stratifying the observed SAM index according to springtime O 3 levels (low minus high O 3 years) allowed the highlighting of positive SAM differences (see Figures 2c and 6b). This outcome occurred toward the next summer, when November O 3 was the driver of SAM changes in summer, as well as toward the previous winter, with late winter/spring SAM being the driver of the O 3 changes. The latter implies that both O 3 and SAM could be influenced by the winter planetary waves without a direct relationship between them. Overall, it appears that the largest changes in the SAM are actually associated with planetary wave disturbances.
We also found that the correlation between the observed O 3 and precipitation tended to be at a maximum during the last decades for all investigated regions (Figure 4a-c). Overall, Australia was the region with the longest persistence of correlation between O 3 and precipitation (Figure 3), although a generally more robust connection between precipitation and Antarctic ozone has been highlighted for Chile (Figures 2d and 4a).
A pronounced drying tendency in Chile with an increase of ozone change contrasts with the moistening tendency found in Australia and, to a lesser extent, South Africa. The magnitude and sign of the observed precipitation difference in spring-summer is similar to those of the simulations for comparable ozone change (e.g., <40%; see Figure 6). On the other hand, for more substantial ozone change, the observed precipitation differences become larger and cover the prior winter in Australia, and have the sign changed from positive to negative in South Africa, and keep increasing without a clear connection with the corresponding SAM difference in Chile ( Figure 6). Because the simulated ozone anomalies were usually smaller than the observed anomalies, we did not attempt to compare the composites associated with the largest ozone changes with the simulations. However, it is worth noting that these largest ozone changes were associated with positive differences in the Niño composite. This is consistent with previous findings suggesting that more frequent vortex disturbances and positive O 3 anomalies are associated with El Niño ( [44], and references therein).
We also evaluated the sensitivity of precipitation to ozone change. Under comparable ozone change, ∆P/∆O 3 for both observations and simulations peaked at around −0.7 in NDJ for Chile and at +0.4 in DJF for Australia, whereas no clear conclusions could be drawn for South Africa.
Overall, our results suggest that springtime ozone anomalies in Antarctica can be used as a predictor of subsequent summer precipitation, at least for Chile and Australia. Because MIROC-ESM-CHEM coupled experiments reproduced most of the observed connections, incorporating ozone variability into seasonal forecasting systems could potentially improve seasonal predictions.
Although our simulations with ODS, held fixed at 1960 levels, showed a less significant linkage between ozone and surface climate, the patterns remained similar to those of the historical simulations. This implies that ozone anomalies can be mostly perceived as a tracer of stratospheric circulation, while stratospheric ozone chemistry contributes to strengthening the relationship between ozone and surface climate. Indeed, under the current high levels of halogen, the Antarctic ozone is more sensitive to dynamical changes. Ozone depletion causes a reduction in shortwave heating, leading to stratospheric cooling and to more stable and persistent polar vortex conditions, which facilitate the formation of polar stratospheric clouds and lead in turn to enhanced ozone depletion. Then, it is worth noting that the future simulations, under the RCP 4.5 scenario, showed a diminishment of the connection between ozone and precipitation, as well as more consistent SAM differences extending to the prior winter, as in the observations.
The simulation results of this study are based solely on a single model (i.e., MIROC-ESM-CHEM). Therefore, potentially, they could be partially biased by some specific issue of this model, such as rainfall [45] and jet climatology [46], sea ice extension [29], and spatial resolution. Nevertheless, the main conclusions appear to be robust and consistent with our preliminary analysis (not shown) based on a subset of chemistry-climate models (CCMs) showing patterns in the ozone-surface climate comparable to MIROC simulations and observations in particular for South America. Future investigations based on a larger suite of CCMs will allow confirmation of such features and investigation of the mechanisms behind such changes. Figure A2. As in Figure 6a except for showing the composite difference (low minus high O3) for the observed Indian Ocean Dipole index (data from the Japan Agency for Marine-Earth Science and Technology). . Note that simulations are not seamless, i.e., historical and future runs cannot be joined, and this results in a gap around present years. The shaded lines indicate where the change was significant at the 95% level based on a two-tailed t-test. Figure A4. As in Figure 2b except for showing the correlation between ENSO and SAM (a) and the correlation between November total ozone and ENSO (b) for given month. The shaded lines indicate where the change was significant at the 95% level based on a two-tailed t-test. . Note that simulations are not seamless, i.e., historical and future runs cannot be joined, and this results in a gap around present years. The shaded lines indicate where the change was significant at the 95% level based on a two-tailed t-test. Figure A2. As in Figure 6a except for showing the composite difference (low minus high O3) for the observed Indian Ocean Dipole index (data from the Japan Agency for Marine-Earth Science and Technology). . Note that simulations are not seamless, i.e., historical and future runs cannot be joined, and this results in a gap around present years. The shaded lines indicate where the change was significant at the 95% level based on a two-tailed t-test. Figure A4. As in Figure 2b except for showing the correlation between ENSO and SAM (a) and the correlation between November total ozone and ENSO (b) for given month. The shaded lines indicate where the change was significant at the 95% level based on a two-tailed t-test.