globalchange  > 气候变化事实与影响
DOI: doi:10.1038/nclimate2568
论文题名:
Causal feedbacks in climate change
作者: Egbert H. van Nes
刊名: Nature Climate Change
ISSN: 1758-960X
EISSN: 1758-7080
出版年: 2015-03-30
卷: Volume:5, 页码:Pages:445;448 (2015)
语种: 英语
英文关键词: Palaeoclimate
英文摘要:

The statistical association between temperature and greenhouse gases over glacial cycles is well documented1, but causality behind this correlation remains difficult to extract directly from the data. A time lag of CO2 behind Antarctic temperature—originally thought to hint at a driving role for temperature2, 3—is absent4, 5 at the last deglaciation, but recently confirmed at the last ice age inception6 and the end of the earlier termination II (ref. 7). We show that such variable time lags are typical for complex nonlinear systems such as the climate, prohibiting straightforward use of correlation lags to infer causation. However, an insight from dynamical systems theory8 now allows us to circumvent the classical challenges of unravelling causation from multivariate time series. We build on this insight to demonstrate directly from ice-core data that, over glacial–interglacial timescales, climate dynamics are largely driven by internal Earth system mechanisms, including a marked positive feedback effect from temperature variability on greenhouse-gas concentrations.

Earth system models9 have been an effective, albeit indirect, way to quantify causality in the climate system. The effects of CO2 and other greenhouse gases (GHGs) on Earth’s temperature are relatively well understood, but estimates of the effect of temperature variability on GHG dynamics remain uncertain10, 11, 12. Quantifying the actual strength of this effect is challenging, because it involves a plethora of mechanisms that are difficult to measure and sometimes oppose each other. For instance, increased photosynthesis at higher CO2 levels implies a negative feedback, whereas enhanced plant and soil respiration at higher temperatures leads to carbon release and a positive feedback13. A warmer climate may induce the release of CO2, CH4 and N2O from terrestrial ecosystems, especially in polar regions14. Furthermore, at higher temperatures, marine CaCO3 neutralization of anthropogenic CO2 decreases15, and methane is released from hydrate storages below the sea floor, which may amplify global warming16. Overall, higher global temperatures are believed to cause a net increase in atmospheric concentrations of GHGs, implying a positive feedback in warming10, 11, 17, 18, 19. However, given the complexity of the mechanisms and models, uncertainty over the feedback effect remains large.

This issue raises the question if there are more direct, model-independent estimates of the feedback effect based on the strikingly parallel dynamics of temperature and GHGs over the Pleistocene ice ages (Fig. 1a). Data-based approaches for unravelling the causation operating behind this correlation have hitherto largely focused on phase lags between past climate data sets3, but these lags vary over time. A slight lead of Antarctic temperature over CO2 variations has been argued to point to temperature as a driver of CO2 changes2. However, more recent studies cast doubt on the existence of a significant time lag of CO2 behind either Antarctic4 or global5 temperature at the last glacial termination, with variations in methane and temperature seeming nearly synchronous at the Bølling transition20. Meanwhile, the latest data on an earlier termination7 and inception6 show periods of significant time lags between CO2 and Antarctic temperature. A simple moving-window scan of optimal time displacement for correlation (Supplementary Fig. 1c) supports the emerging view that the time lag of CO2 behind temperature as recorded in the Vostok ice core1 has varied widely over the past 400 kyr. Although errors in dating may contribute to such variation, detailed recent studies6, 7, 21 confirm that these lags do vary substantially over time.

Figure 1: Causation inferred from time series of insolation, temperature and GHGs.
Causation inferred from time series of insolation, temperature and GHGs.

a, Fluctuations in orbitally driven insolation (65° N) and in temperature and GHGs inferred from the Vostok ice core1. b, Schematic of causality inferred from the patterns in this study. Temperature and GHGs exhibit strong feedbacks, whereas orbital forcing has a relatively small role in influencing temperature.

We used the Vostok ice core1 on local temperature, CO2, CH4, sodium, dust (http://www.ncdc.noaa.gov/paleo/icecore/antarctica/vostok/vostok_data.html) and July insolation28 at 65° N (http://www1.ncdc.noaa.gov/pub/data/paleo/climate_forcing/orbital_variations/berger_insolation). All time series were linearly interpolated to produce equidistant estimates spaced by 1 kyr. The CCM algorithm follows ref. 8, which is based on nonlinear state space reconstruction. This method reconstructs the manifold of the dynamical system based on one variable X only. E time-lagged values of X (time lags 0, τ, 2τ, …(E − 1)τ) are used as coordinate axes to reconstruct this ‘shadow attractor manifold’. The algorithm then finds points on this shadow manifold that are close together and tests whether paired observations of another variable Y are also close together. This is done by predicting each value of Y on the basis of the closest points in X using simplex projection (see details in ref. 8). After this each predicted value of Y is compared with the observed Y (using Pearson’s correlation ρ). This procedure is repeated using a subset of the time series of X with different lengths L. It is expected that the prediction improves with the length of the time series L until it converges to a maximum level. To measure the convergence we fitted an exponential function (ρ = ρmaxρ0ec(LL0)) to the relationship between predicted and observed values based on the different subsets of X. L = length of the used subset, L0 = first length used = 10, ρmax the maximum correlation coefficient, ρ0 the correlation at L = 10 and c a convergence speed. The cross-mapping variables are labelled following the convention of ref. 8, where ‘Y xmap X’ quantifies the causal effect of X on Y by predicting Xt from E lagged time-series fragments of Yt. We used E = 4, τ = 2 kyr by default as this combination gave a good unfolding of the attractor (visual inspection of Supplementary Fig. 14) and above E = 4 there was no clear improvement of the predictability of temperature and the GHGs (Supplementary Fig. 15) using simplex projection29. We also tested other values of embedding dimension E and embedding lag τ, which gave very similar results (Supplementary Fig. 3).

Robustness was tested by varying CCM parameters (Supplementary Fig. 3), exploring other ways of interpolation (Supplementary Fig. 4 and Supplementary Tables 1 and 2), and repeating the analysis on the oldest part of the EPICA ice-core data24, 25 and a high-resolution part of the ice cores of the past 22 kyr (ref. 21; Supplementary Fig. 2). In addition we examined the significance of the results using two conservative null models. For the default null model we generated 100 surrogate data sets which were randomly shifted in phase by choosing a random break point and swapping the order of both segments. This procedure destroys the dynamic interdependency between time series, but preserves nearly all short-term behaviour. We also used a second null model, generating 100 time series for each variable having the same frequency spectrum as observed, but with the frequencies randomly shifted in phase30 (see Supplementary Fig. 5). With each of these 4 × 100 time series all original variables were predicted using cross-mapping. An interaction is considered to be significant only if the CCM skill of the real time series is outside the range of the 5th and 95th percentiles computed from the randomly generated time series (that is, in Fig. 2 lines above the shaded areas indicate significant CCM).

For synchronous (that is, strongly correlated) variables, one-way causation is hard to distinguish from bi-directional causation. Therefore we also analysed the effect of time displacements on CCM skill. This was done by displacing the time series up to 10 kyr backwards and forwards before measuring the CCM skill using 500 bootstrapped library sets of length 100. The optimum CCM time lag was determined by finding the optimum in a Gaussian filtered relation (bandwidth = 5) of the CCM skill as function of the time lag. We did the same in a null model (see Supplementary Information).

  1. Petit, J. R. et al. Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica. Nature 399, 429436 (1999).
  2. Fischer, H., Wahlen, M., Smith, J., Mastroianni, D. & Deck, B. Ice core records of atmospheric CO2 around the last three glacial terminations. Science 283, 17121714 (1999).
  3. Shackleton, N. J. The 100,000-year ice-age cycle identified and found to lag temperature, carbon dioxide, and orbital eccentricity. Science 289, 18971902 (2000).
  4. Pedro, J. B., Rasmussen, S. O. & Van Ommen, T. D. Tightened constraints on the time-lag between Antarctic temperature and CO2 during the last deglaciation. Clim. Past 8, 12131221 (2012).
  5. Shakun, J. D. et al. Global warming preceded by increasing carbon dioxide concentrations during the last deglaciation. Nature 484, 4954 (2012).
  6. Landais, A. et al. Two-phase change in CO2, Antarctic temperature and global climate during Termination II. Nature Geosci. 6, 10621065 (2013).
URL: http://www.nature.com/nclimate/journal/v5/n5/full/nclimate2568.html
Citation statistics:
资源类型: 期刊论文
标识符: http://119.78.100.158/handle/2HF3EXSE/4788
Appears in Collections:气候变化事实与影响
科学计划与规划
气候变化与战略

Files in This Item:
File Name/ File Size Content Type Version Access License
nclimate2568.pdf(448KB)期刊论文作者接受稿开放获取View Download

Recommended Citation:
Egbert H. van Nes. Causal feedbacks in climate change[J]. Nature Climate Change,2015-03-30,Volume:5:Pages:445;448 (2015).
Service
Recommend this item
Sava as my favorate item
Show this item's statistics
Export Endnote File
Google Scholar
Similar articles in Google Scholar
[Egbert H. van Nes]'s Articles
百度学术
Similar articles in Baidu Scholar
[Egbert H. van Nes]'s Articles
CSDL cross search
Similar articles in CSDL Cross Search
[Egbert H. van Nes]‘s Articles
Related Copyright Policies
Null
收藏/分享
文件名: nclimate2568.pdf
格式: Adobe PDF
此文件暂不支持浏览
所有评论 (0)
暂无评论
 

Items in IR are protected by copyright, with all rights reserved, unless otherwise indicated.