Abstract
In the context of global warming, the question of why Antarctic sea ice extent (SIE) has increased is one of the most fundamental unsolved mysteries. Although many mechanisms have been proposed, it is still unclear whether the increasing trend is anthropogenically originated or only caused by internal natural variability. In this study, we employ a new method where the underlying natural persistence in the Antarctic SIE can be correctly accounted for. We find that the Antarctic SIE is not simply shortterm persistent as assumed in the standard significance analysis, but actually characterized by a combination of both short and longterm persistence. By generating surrogate data with the same persistence properties, the SIE trends over Antarctica (as well as five subregions) are evaluated using MonteCarlo simulations. It is found that the SIE trends over most subregions of Antarctica are not statistically significant. Only the SIE over Ross Sea has experienced a highly significant increasing trend (p = 0.008) which cannot be explained by natural variability. Influenced by the positive SIE trend over Ross Sea, the SIE over the entire Antarctica also increased over the past decades, but the trend is only at the edge of being significant (p = 0.034).
Introduction
Polar sea ice, as an important component of the climate system, can affect the earth system in many ways. It hinders exchanges between ocean and atmosphere, reflects solar radiation back to space, and even changes the habitat and species of plants and animals in the polar region^{1,2}. In the context of global warming, it is of great importance to monitor the variations of polar sea ice and further study the potential effects of polar sea ice changes^{3,4}. In the Arctic, a significant decreasing trend of the sea ice extent (SIE) has been well recognized, which is due to the Arctic warming^{5,6}. However, unlike the Arctic SIE, recent satellite observations surprisingly reveal a positive trend in the Antarctic SIE^{7,8}, which confounds the most trusted climate models in the world, including the Coupled Model Intercomparison Project Phase 5 (CMIP5) models^{9,10}. Therefore, understanding the unexpected increasing trend of the Antarctic SIE has become a research hotspot in the past decades.
Among the natural explanations are the effects of El Niño events, the Southern Annular Mode (SAM), wind patterns as well as geology, etc.^{11,12,13,14}. Among the anthropogenic explanations are, apart from global warming, the ozone depletion, which may both deepen the Amundsen Sea Low (ASL), strengthen the southernly winds over the Ross Sea, and finally lead to an increasing SIE trend over this area^{15,16,17,18}. The ozone hole also may induce a positive SAM trend, which in turn may lead to an increasing trend of the Antarctic SIE. However, coupled atmosphereoceansea ice models suggest that the ozone depletion may lead to a negative trend of the Antarctic SIE^{19,20,21}.
Up to now, the reasons why the SIE shows a positive trend are still in debate. But beyond the reasons for a positive trend, a more important question is, whether the trend is originated from anthropogenic forcings, or within the bounds of natural variability. All previous calculations of the trend significance concluded that the positive trend of the entire Antarctic SIE over the past 30 years is statistically highly significant^{8}, which indicates the positive trend is unexplainable using natural variability. However, from model simulations, previous studies argued that it may be difficult to attribute the positive SIE trend over Antarctica to anthropogenic forcings, since the distribution of the simulated natural trend may be wide enough to cover the observed Antarctic sea ice trend^{15,22,23,24}. Therefore, it is still an open question whether the Antarctic SIE trend is significant or not. To make a better estimation of the trend significance, it is crucial to know its accurate natural variability including its persistence properties.
Figure 1 shows the time evolution of the Antarctic SIE anomalies between November 1978 and December 2013 in the five Antarctic subregions^{25}. The SIE anomalies were calculated by subtracting the annual cycle (see the “Data” section). To show the temporal variation more clearly, we made a moving average over 30 months (see the green curves). In the figure, one can see clearly a pronounced mountainvalley structure which represents the variability caused by the natural persistence and anthropogenic forcings. In the previous calculations of the trend significance, it has been assumed that the natural part of the Antarctic SIE exhibits only shortterm persistence, which can be modeled by an autoregression process of first order (AR(1)). In this case, the autocorrelation function C(s) decays exponentially with time lag s as , and the statistical significance of the Antarctic SIE trend can be determined from the classic method developed by^{26} (see “Method” section). However, it has been well recognized that another type of persistence, longterm persistence, is more relevant in climate^{27} (see “Method” section). Different from shortterm persistence, longterm persistence indicates long lasting influences from the past histories and its autocorrelation function C(s) decays asymptotically as C(s) ~ (1 − γ)s^{−γ} (s > 0, see “Method” section). In general, both short and longterm persistence can contribute to the observed natural variability, and a correct estimation of the underlying persistence is a prerequisite for a proper estimation of the trend significance^{28,29,30,31,32}. So far, this kind of analysis has not been done for the Antarctic SIE.
In this study, we apply Detrended Fluctuation Analysis of second order^{33} (DFA2, see “Method” section) to the monthly SIE anomalies to estimate the persistence of the SIE. We show explicitely that the persistence can be modeled by a superposition of short and longterm persistent processes. By using extensive MonteCarlo simulations, we determine the desired distribution of natural trends, which finally allows us to determine how significant the SIE trends in the Antarctic regions are.
Figure 2 shows the DFA2 fluctuation function F(s) of the SIE for the 5 regions and the entire Antarctica (black dots). In the doublelogarithmic presentation one can see that the fluctuation functions F(s) show, for all SIE data, a crossover at s between 18 and 20 months. The steeper slopes below the crossovers indicate shortterm persistence. Above the crossover, F(s) increases by a power law, F(s) ~ s^{h}, with the Hurst exponent h between 0.6 and 0.8, indicating longterm persistence on large time scales. Accordingly, the case of pure shortterm persistence where h = 0.5 can be excluded. The Hurst exponent characterizes the strength of the longterm persistence: the larger h the higher the persistence. The figure suggests that the natural variability of the Antarctic SIE can be described by a superposition of short and longterm persistence processes, which can be modeled by (see “Method” section)
where a is the AR(1) parameter and η_{h}(i) is longterm persistent noise. We like to note that equation (1) reduces to the standard AR(1) model when η_{h}(i) is substituted by Gaussian white noise. On the other hand, for a = 0, it reduces to the standard model of longterm persistence.
It is not surprising that the Antarctic SIE shows, on large time scales, longterm persistence, since also the surrounding sea and air temperatures are known to be longterm persistent^{30,31,32}. A possible origin of the longterm persistence of the Antarctic SIE is the internal dynamics of sea ice. As discussed in refs 34 and 35, processes including the Southern Annular Mode (SAM), El NiñoSouthern Oscillation (ENSO), highfrequency weather forcings, as well as the South Pacific Intrinsic Mode (SPIM) can drive the multiscale sea ice variability, which may further induce the arise of longterm persistence.
Using equation (1), we have determined the appropriate parameters a and h for all Antarctic regions by best fits of the fluctuation functions F(s). The values of a and h are shown in Fig. 2. For all five regions, the estimated AR(1) parameters are similar ranging from a = 0.5 to 0.55, indicating the shortterm persistence of SIE over different Antarctic subregions may have similar dynamical sources. Regarding the longterm persistence, the parameter h varies from 0.6 in the “Ross Sea” and “Pacific Ocean”, to 0.8 in the “Indian Ocean”, which may imply different mechanisms in producing the longterm persistence. As discussed in refs 34 and 35, the different roles of SAM, ENSO, as well as the highfrequency weather forcings in driving seaice variability over different regions, may contribute to the different h values. Using these values of a and h, as the figure shows, the DFA2 results of the simulated data agree well with those of the observed SIE data.
For each SIE data set, we then used equation (1) with the proper a and h values, to generate more than 100,000 artificial data sets with the same length L = 422 which by definition mimic the natural variability of the considered SIE data set. By determining the relative trend x in each artificial record, we then arrive, for each of the 6 SIE data sets, at the desired probability density function P(x) and the related statistical significance S(x) of the observed trend, as well as the pvalues p(x) = 1 − S(x) (see “Method” section).
Figure 3 shows the pvalues for the artificial data sets (green line) simulating the natural variability of (a) the Weddell Sea, (b) Bellingshausen and Amundsen Seas, (c) Ross Sea, (d) Indian Ocean, (e) Pacific Ocean, and (f) entire Antarctica. In order to reveal how sensitively the statistical significance depends on the assumption of the underlying persistence of each record, we show also p(x) for (i) white noise and (ii) the corresponding AR(1) process. For both white noise and AR(1) processes, p(x) can be obtained straightforwardly from equation (7) in the Methods Section. In addition, (iii) we have determined p(x) for longterm persistent data with the same h values and the same length L = 422 as the original SIE data. For details of the calculation we refer to ref. 29.
The vertical dashed line in each panel represents the relative trend x of the SIE in each of the 6 Antarctic regions. The intersection of this vertical line with p(x) yields the desired p value (horizontal dashed line). The figure shows the dependence of p on the considered statistical model for the natural variability. Under the assumption of white noise (black line), the SIE trend in all regions is highly significant, with pvalues well below 0.01. Under the standard assumption of an AR(1) process (red line), the SIE trends in Ross Sea, Indian Ocean, and Antarctica are highly significant, with p values well below 0.01. For the other 3 regions, the trends are not significant since p is above 0.05. Under the assumption of a purely longterm persistent process (blue line), the trends in the Ross Sea and entire Antarctica are highly significant, with p values below 0.01. For the appropriate assumption of a combined short and longterm persistent process (equation (1)), we find that the trend in the SIE is not significant in 4 of the 5 Antarctic subregions where the p values are well above 0.2. Only in the Ross Sea the trend is highly significant (p = 0.008). This result suggests that only the positive trend of the SIE in the Ross Sea cannot be explained by its natural variability, whereas the trends in the other 4 subregions are well within the bounds of natural variability. In the entire Antarctica, the trend is only slightly significant, with p = 0.034 only little below the significance threshold 0.05.
In addition to the significance, we have also determined the maximum and minimum external trend (see equation (5) and (6) in the “Method” section) for the significance levels α = 0.05 and 0.01. Figure 4 summarizes all results for the SIE in the 6 Antarctic regions, (a) for the standard assumption of an AR(1) process, and (b) for the combined AR(1) and longterm persistent process. The numbers shown are the pvalues. The upper and lower bounds of the boxes represent the maximum and minimum external trend for the significance level α = 0.05, while the error bars stand for the maximum and minimum external trend for α = 0.01. The figure again shows that the assumption of an AR(1) process strongly underestimates the uncertainties in the Antarctic SIE and overestimates the significance. A remarkable difference between (a) and (b) occurs in the Indian Ocean. The conventional AR(1) treatment yields that the increase of the SIE is highly significant, with p = 0.008. On the other hand, when the considerable longterm persistence (h = 0.8) of the SIE in the Indian Ocean is taken into account, p becomes 0.26, which is well above the significance level 0.05.
Discussion and Conclusion
In this study, we presented evidence that the temporal evolution of the Antarctic SIE is characterized by a superposition of both, shortterm persistent and longterm persistent processes. By combining both processes (see equation (1)), we succeeded in obtaining a realistic statistical model for the natural variability of the Antarctic SIE. The model allowed us to determine, by MonteCarlo simulations, the natural distribution of trends and thus the statistical significance of the trends in the Antarctic SIE. We find that already in the considered small time window of 35 years the SIE trend over Ross Sea is highly statistically significant with p = 0.008, while the trends of the other four subregions are well within the bounds of natural variability, with pvalues well above 0.2. Accordingly, only the trend in the Ross Sea SIE cannot be explained by natural variability alone and external forcings must contribute to this trend. For the other four subregions, the SIE trends are not statistically significant, either because there is no anthropogenic trend or because the considered time period is not long enough to distinguish a deterministic anthropogenic trend from the natural fluctuations of the records. The SIE over the entire Antarctica increased over the past decades, but the trend is only at the edge of being significant (p = 0.034).
Our work, from a new perspective, resolved the discrepancy reported by the previous studies. On one hand, statistical models based on shortterm persistence only, suggest that the Antarctic SIE trends are highly significant and cannot be explained by natural variability^{8}. On the other hand, dynamical model simulations suggested that the observed SIE trends may be not strong enough to exceed the bounds of the natural variabilities of SIE^{22,24}. But it is not sure how trustworthy the simulations are, since they are unable to reproduce the observed increase in the Antarctic SIE. Here we have shown that the main reason for the great discrepancy between the previous statistical analysis and the dynamical model simulations is the inappropriate assumption of shortterm persistence in the statistical model. Possible defects in the ability of the dynamical models in simulating the Antarctic SIE^{9,10} do not play a major role here.
Based on our analysis, Ross Sea is the only subregion that has been affected strongly by the anthropogenic forcings. The mechanism may be as follows^{15,16,17,18,36}: Both global warming and ozone depletion may deepen the Amundsen Sea Low (ASL), weaken the westerlies near the Ross Sea and strengthen the southerly cold winds from the continent. This way, the sea ice in the Ross Sea may become more isolated and increase due to the southerly cold winds.
The fact that the CMIP5 models are not able to reproduce the increase of the SIE in Antarctica, points to a remarkable deficiency of the models, even though the models tried to implement the effects of the ozone depletion^{19,20,21}. Since the Ross Sea is the main driver of this increase, we suggest that the models concentrate on the Ross Sea in order to see what went wrong and how to improve the performance.
Data and Methods
Data
In this study, the trends of the SIE over Antarctica are investigated using the monthly sea ice data set from the US National Snow and Ice Data Center (NSIDC) (Data contributors: Stroeve, J. and Meier, W., http://nsidc.org/data/nsidc0192). Monthly mean sea ice concentrations were used to derive the extents, which is derived from passive microwave measurements at a grid cell size of 25 × 25 km from Nimbus7 SMMR; DMSPF8, F11, and F13 SSM/I; and from DMSPF17 SSMIS. NASA team Algorithms is applied in computing the total SIE. Total SIE is derived by summing the number of pixels with at least 15 percent ice concentration multiplied by the area per pixel. In this way, the entire area of any pixel with at least 15 percent ice concentration is considered to contribute to the total SIE. According to ref. 37, it is suggested that using monthly mean sea ice concentrations may provide a more accurate monthly ice extent since the effects of storms on the ice edge are less pronounced in the monthly mean sea ice data than in the daily sea ice data. Besides the SIE of the entire Antarctic, the SIE over the sector of Bellingshausen and Amundsen Sea, Southern Pacific Ocean, Southern Indian Ocean, Weddell Sea and Ross Sea is also calculated, with the same protocol of Zwall et al.^{25} and Turner et al.^{15}. SIE data employed in this study ranges from November 1978 to December 2013 (422 months in total). Before analysis, the SIE anomalies are calculated by subtracting the annual cycle, which is defined as the averages of each calendar month (see Fig. 1b).
Methods
Short and longterm persistence
In records with shortterm persistence, the autocorrelation function C(s) of the (detrended) anomalies decays exponentially as . In records with longterm persistence, C(s) decays algebraically as C(s) ≅ (1 − γ)s^{−γ}, s > 0, 0 < γ < 1. In some cases like the Antarctic SIE, the persistence of the record is characterized by the superposition of short and longterm persistence. All the above three cases can be modeled as following,
where a is the AR(1) parameter and η_{i} is longterm persistent noise. Basically, this is a generalization of an autoregression process of first order (AR(1)). By setting η_{i} as white noise, equation (2) will be degenerated as AR(1) process, which is useful for the modeling of shortterm persistent records. In this case, the persistent time . By setting the AR(1) parameter as a = 0, equation (2) will be able to model purely longterm persistent data. In this case, the persistent time s_{x} diverges. While by using appropriate AR(1) parameter a and longterm persistent data η_{i}, equation (2) can be used to simulate the natural variability of data with combined short and longterm persistence. Numerically, a longterm persistent data set {η_{1}, η_{2} … η_{N}} can be generated by the Fourier Filtering Technique^{38}. The data sets are characterized by the correlation exponent γ or the Hurst exponent h (see below). Then, by choosing y_{1} randomly from a Gaussian distribution, the values for y_{2}, y_{3} … y_{N} can be calculated recursively from (2).
Detrended Fluctuation Analysis
Since C(s) shows strong finite size effects, one usually considers the Detrended Fluctuation Analysis 2 (DFA2) to detect the underlying persistence in the record of interest. In DFA2, one first divides the record {y_{i}}, i = 1, …, L, into nonoverlapping windows μ of lengths s. Then one focuses, in each segment μ, on the cumulated sum Y_{i} of the data and determines the variance of the Y_{i} around the best polynomial fit of order 2. After averaging over all segments μ and taking the square root, one arrives at the desired fluctuation function F(s)^{33}. F(s) reflects the different kinds of persistence in the following way:
 i
Without any persistence, for simple white noise, F(s) increases as s^{1/2}.
 ii
In shortterm persistent records with a crossover time s_{x}, F(s) shows a crossover. For s well above s_{x}, F(s) increases as s^{1/2} as for white noise, while below s_{x}, F increases stronger than s^{1/2}. For very large s_{x}, F increases as s^{3/2} for s below s_{x}.
 iii
In purely longterm persistent records, F(s) increases as
where the exponent h can be associated with the Hurst exponent and is related to the correlation exponent γ by h = 1 − γ/2. For stationary records, the dimensionless Hurst exponent h is between 0.5 and 1; h characterizes the strength of the longterm persistence. For increasing h, the strength of the longterm persistence increases and the mountain valley structure caused by the persistence becomes more pronounced.
When both short and longterm persistence are present, F(s) shows a crossover similar to that one of shortterm persistent records. For s below s_{x}, F increases as in shortterm persistent records, while for s well above s_{x}, F increases as for purely longterm persistent records.
Significance of trends
We consider the record {y_{i}} with length L. From the regression line r_{i} = ai + b, we obtain the magnitude of the trend Δ = a(L − 1) and the fluctuations around the trend, characterized by the standard deviation . The relevant quantity we are interested in is the relative trend
For assessing if an observed trend in a data set may be due to its natural persistence or not, one needs to know the probability P(x, L)dx that in model records with the same persistence properties as the considered data set, a relative trend between x and x + dx occurs. From P one obtains the statistical significance S of the trend x, . The deviation of S from 1 is called pvalue.
If p(x, L) is below a certain significance level α (usually α is 0.05 or 0.01), one usually assumes that the considered trend cannot be fully explained by the natural variability of the record. The relation p(x_{α};L) = 1 − S(x_{α};L) = α defines the upper and lower limits ±x_{α} of the considered significance interval (also called confidence interval). By the above assumption, relative trends x between −x_{α}(L) and x_{α}(L) can be regarded as natural. If x is above x_{α}, the part x − x_{α} cannot be explained by the natural variability of the record and thus can be regarded as minimum external relative trend,
Similarly, the external relative trend cannot exceed
which thus represents the maximum external relative trend. According to (5) and (6), ±x_{α}(L) can be regarded as error bars for an external relative trend in a record of length L.
For Gaussian white noise, S is given by
with l(L) = L − 2 and a = ((L − 1)/)/ ≅ . _{2}F_{1} is the hypergeometric function and Γ denotes the Γfunction.
For shortterm persistent records described by an autoregressive process of first order (AR(1)), Eq. (7) remains the same, only l(L) has to be substituted by L(1 − C(1))/(1 + C(1)), where C(1) is the autocorrelation with lag of 1 (months).
For longterm persistent records characterized by a Hurst exponent h, S(x, L) has the same functional form as in Eq. (7). The parameters l and a depend on L and h and have been tabulated in ref. 39.
Additional Information
How to cite this article: Yuan, N. et al. Increase of the Antarctic Sea Ice Extent is highly significant only in the Ross Sea. Sci. Rep. 7, 41096; doi: 10.1038/srep41096 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
 1
Parkinson, C. L. Southern ocean sea ice and its wider linkages: insights revealed from models and observations. Antarctic Science. 16, 387–400 (2004).
 2
DeConto, R. M. & Pollard, D. Contribution of Antarctica to past and future sealevel rise. Nature. 531, 591–597 (2016).
 3
Yuan, X. & Martinson, D. G. Antarctic sea ice extent variability and its global connectivity. J. Climate. 13, 1697–1717 (2000).
 4
Gao, Y. et al. Arctic sea ice and Eurasian climate: A review. Adv. Atmos. Sci. 32, 92–114 (2015).
 5
Cavalieri, D. J. & Parkinson, C. L. 30Year satellite record reveals contrasting Arctic and Antarctic decadal sea ice variability. Geophys. Res. Lett. 30, 1970 (2003).
 6
Cavalieri, D. J. & Parkinson, C. L. Arctic sea ice variability and trends, 19792010. The Cryosphere. 6, 881–889 (2012).
 7
Cavalieri, D. J. & Parkinson, C. L. Antarctic sea ice variability and trends, 19792006. J. Geophys. Res. 113, C07004 (2008).
 8
Parkinson, C. L. & Cavalieri, D. J. Antarctic sea ice variability and trends, 19792010. The Cryosphere. 6, 871–880 (2012).
 9
Turner, J., Bracegirdle, T. J., Phillips, T., Marshall, G. J. & Hosking, J. S. An initial assessment of Antarctic sea ice extent in the CMIP5 models. J. Climate. 26, 1473–1484 (2013).
 10
Purich, A., Cai, W., England, M. H. & Cowan, T. Evidence for link between modelled trends in Antarctic sea ice and underestimated westerly wind changes. Nat. Commun. 7, 10409 (2016).
 11
Rind, D., Chandler, M., Lerner, J., Martinson, D. G. & Yuan, X. Climate response to basinspecific changes in latitudinal temperature gradients and implications for sea ice variability. J. Geophys. Res. 106, 20161–20173 (2001).
 12
Stammerjohn, S. E., Martinson, D. G., Smith, R. C., Yuan, X. & Rind, D. Trends in Antarctic annual sea ice retreat and advance and their relation to El NiñoSouthern Oscillation and Southern Annular Mode variability. J. Geophys. Res. 113, C03S90 (2008).
 13
Holland, P. R. & Kwok, R. Winddriven trends in Antarctic seaice drift. Nat. Geosci. 5, 872–875 (2012).
 14
Nghiem, S. V., Rigor, I. G., ClementeColon, P., Neumann, G. & Li, P. P. Geophysical constraints on the Antarctic sea ice cover. Remote Sens. Environ. 181, 281–292 (2016).
 15
Turner, J. et al. Nonannular atmospheric circulation change induced by stratospheric ozone depletion and its role in the recent increase of Antarctic sea ice extent. Geophys. Res. Lett. 36, L08502 (2009).
 16
Thompson, D. W. J. et al. Signatures of the Antarctic ozone hole in Southern Hemisphere surface climate change, Nat. Geosci., 4, 741–749 (2011).
 17
Previdi, M. & Polvani, L. M. Climate system response to stratospheric ozone depletion and recovery. Q. J. R. Meteorol. Soc. 140, 2401–2419 (2014).
 18
Raphael, M. N. et al. The Amundsen Sea Low variability, change, and impact on Antarctic climate. Bull. Amer. Meteor. Soc. 97, 111–121 (2016).
 19
Sigmond, M. & Fyfe, J. C. Has the ozone hole contributed to increased Antarctic sea ice extent? Geophys. Res. Lett. 37, L18502 (2010).
 20
Bitz, C. M. & Polvani, L. M. Antarctic climate response to stratospheric ozone depletion in a fine resolution ocean climate model. Geophys. Res. Lett. 39, L20705 (2012).
 21
Sigmond, M. & Fyfe, J. C. The Antarctic sea ice response to the ozone hole in climate models. J. Climate. 27, 1336–1342 (2014).
 22
Polvani, L. M. & Smith, K. L. Can natural variability explain observed Antarctic sea ice trends? New modeling evidence from CMIP5. Geophys. Res. Lett. 40, 3195–3199 (2013).
 23
Turner, J., Hosking, J. S., Bracegirdle, T. J., Marshall, G. J. & Phillips, T. Recent changes in Antarctic sea ice Phil. Trans. Roy. Soc. A 373, 20140163 (2015).
 24
Turner, J., Hosking, J. S., Marshall, G. J., Phillips, T. & Bracegirdle, T. J. Antarctic sea ice increase consistent with intrinsic variability of the Amundsen Sea Low. Clim. Dyn. 46, 2391–2402 (2016).
 25
Zwally, H. J., Comiso, J. C., Parkinson, C. L. & Cavalieri, D. J. Variability of Antarctic sea ice 1979–1998. J. Geophys. Res. 107, 3041 (2002).
 26
Santer, B. D. et al. Statistical significance of trends and trend differences in layeraverage atmospheric temperature time series. J. Geophys. Res. 105, 7337–7356 (2000).
 27
KoscielnyBunde, E. et al. Indication of a universal persistence law governing atmospheric variability. Phys. Rev. Lett. 81, 729–732 (1998).
 28
Cohn, T. A. & Lins, H. F. Nature’s style: Naturally trendy. J. Geophys. Res. 32, L23402 (2005).
 29
Lennartz, S. & Bunde, A. Trend evaluation in records with longterm memory: Application to global warming. Geophys. Res. Lett. 36, L16706 (2009).
 30
Bunde, A., Ludescher, J., Franzke, C. L. E. & Büntgen, U. How significant is West Antarctic warming. Nat. Geosci. 7, 246–247 (2014).
 31
Ludescher, J., Bunde, A., Franzke, C. E. & Schellnhuber, H. J. Longterm persistence enhances uncertainty about anthropogenic warming of Antarctica. Clim. Dyn. 457, 459–462 (2015).
 32
Yuan, N. et al. On the longterm climate memory in the surface air temperature records over Antarctica: A nonnegligible factor for trend evaluation. J. Climate. 28, 5922–5934 (2015).
 33
Kantelhardt, J. W., KoscielnyBunde, E., Rego, H. H. A., Havlin, S. & Bunde, A. Detecting longrange correlations with detrended fluctuation analysis. Physica A. 295, 441–454 (2001).
 34
O’Kane, T. J. et al. Decadal variability in an OGCM southern ocean: Intrinsic modes, forced modes and metastable states. Ocean Model. 69, 1–21 (2013).
 35
Matear, R. J., O’Kane, T. J., Risbey, J. S. & Chamberlain, M. Sources of heterogeneous variability and trends in Antarctic seaice. Nat. Commun. 6, 8656 (2015).
 36
Scarchilli, C., Frezzotti, M. & Ruti, P. M. Snow precipitation at four ice core sites in East Antarctica: provenance, seasonality and blocking factors. Clim. Dyn. 37, 2107–2125 (2011).
 37
Stroeve, J. Sea ice trends and climatologies from SMMR and SSM/ISSMIS, Version 1. [Indicate subset used]. Boulder. Colorado USA. NASA National Snow and Ice Data Center Distributed Active Archive Center. doi: http://dx.doi.org/10.5067/C2HMG48G83QQ. (2003).
 38
Turcotte, D. Fractals and Chaos in Geology and Geophysics. second edition, (Cambridge University Press, 1997).
 39
Tamazian, A., Ludescher, J. & Bunde, A. Significance of trends in longterm correlated records. Phys. Rev. E. 91, 032806 (2015).
 40
Schlitzer, R. Ocean Data View, http://odv.awi.de, (2016).
Acknowledgements
Many thanks are due to supports from National Key Research and Development Program (2016YFA0600404) and from National Key Basic Research Program of China (2013CBA01804). This work is also supported by National Natural Science Foundation of China (No. 41405074, No. 41675088, and No. 41425003) and the LOEWE excellence cluster FACE2FACE of the Hessen State Ministry of Higher Education, Research and the Arts. N. Y. thanks also the supports from CAS Pioneer Hundred Talents Program. We also thank Dr. Tingfeng Dou for helpful discussion.
Author information
Affiliations
Contributions
N.Y. and A.B. designed the study. N.Y., J.L. and M.D. performed the study. N.Y. and A.B. wrote the manuscript. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Yuan, N., Ding, M., Ludescher, J. et al. Increase of the Antarctic Sea Ice Extent is highly significant only in the Ross Sea. Sci Rep 7, 41096 (2017). https://doi.org/10.1038/srep41096
Received:
Accepted:
Published:
Further reading

How to determine the statistical significance of trends in seasonal records: application to Antarctic temperatures
Climate Dynamics (2021)

Role of seaice initialization in climate predictability over the Weddell Sea
Scientific Reports (2019)

Detecting the statistical significance of the trends in the Antarctic sea ice extent: an indication for a turning point
Climate Dynamics (2019)

On memory and nonmemory parts of surface air temperatures over China: can they be simulated by decadal hindcast experiments in CMIP5?
Climate Dynamics (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.