Reconstructing Past Climate by Using Proxy Data and a Linear Climate Model

Thesis (Master's)--University of Washington, 2016-06

(Data D1) Annual climate proxy data from the NGRIP ice core, Greenland

The last two abrupt warmings at the onset of our present warm interglacial period, interrupted by the Younger Dryas cooling event, were investigated at high temporal resolution from the North Greenland Ice Core Project ice core. The deuterium excess, a proxy of Greenland precipitation moisture source, switched mode within 1 to 3 years over these transitions and initiated a more gradual change (over 50 years) of the Greenland air temperature, as recorded by stable water isotopes. The onsets of both abrupt Greenland warmings were slightly preceded by decreasing Greenland dust deposition, reflecting the wetting of Asian deserts. A northern shift of the Intertropical Convergence Zone could be the trigger of these abrupt shifts of Northern Hemisphere atmospheric circulation, resulting in changes of 2 to 4 kelvin in Greenland moisture source temperature from one year to the next.

Orbital-resolution Pliocene-Pleistocene simulations of CO2, sea level, global temperature and benthic d18O, and d11B-based proxy-CO2 data

During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.

Dynamical structure in paleoclimate data

Deterministic chaos in dynamical systems offers a new paradigm for understanding irregular fluctuations. The theory of chaotic dynamical systems includes methods that can test whether any given set of time series data, such as paleoclimate proxy data, are consistent with a deterministic interpretation. Paleoclimate data with annual resolution and absolute dating provide multiple channels of concurrent time series; these multiple time series can be treated as potential phase space coordinates to test whether interannual climate variability is deterministic. Dynamical structure tests which take advantage of such multichannel data are proposed and illustrated by application to a simple synthetic model of chaos, and to two paleoclimate proxy data series.

Testing the timing of radiocarbon-dated events between proxy archives

For interpreting past changes on a regional or global scale, the timings of proxy-inferred events are usually aligned with data from other locations. However, too often chronological uncertainties are ignored in proxy diagrams and multisite comparisons, making it possible for researchers to fall into the trap of sucking separate events into one illusionary event (or vice versa). Here we largely solve this "suck in and smear syndrome" for radiocarbon (14C) dated sequences. In a Bayesian framework, millions of plausible age-models are constructed to quantify the chronological uncertainties within and between proxy archives. We test the technique on replicated high-resolution 14C-dated peat cores deposited during the "Little Ice Age" (c. AD 1400-1900), a period characterized by abrupt climate changes and severe 14C calibration problems. Owing to internal variability in proxy data and uncertainties in age-models, these (and possibly many more) archives are not consistent in recording decadal climate change. Through explicit statistical tests of palaeoenvironmental hypotheses, we can move forward to systematic interpretations of proxy data. However, chronological uncertainties of non-annually resolved palaeoclimate records are too large for answering decadal timescale questions.

Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr – Part 1: A new geomagnetic data composite

We present a new composite of geomagnetic activity which is designed to be as homogeneous in its construction as possible. This is done by only combining data that, by virtue of the locations of the source observatories used, have similar responses to solar wind and IMF (interplanetary magnetic field) variations. This will enable us (in Part 2, Lockwood et al., 2013a) to use the new index to reconstruct the interplanetary magnetic field, B, back to 1846 with a full analysis of errors. Allowance is made for the effects of secular change in the geomagnetic field. The composite uses interdiurnal variation data from Helsinki for 1845–1890 (inclusive) and 1893–1896 and from Eskdalemuir from 1911 to the present. The gaps are filled using data from the Potsdam (1891–1892 and 1897–1907) and the nearby Seddin observatories (1908–1910) and intercalibration achieved using the Potsdam–Seddin sequence. The new index is termed IDV(1d) because it employs many of the principles of the IDV index derived by Svalgaard and Cliver (2010), inspired by the u index of Bartels (1932); however, we revert to using one-day (1d) means, as employed by Bartels, because the use of near-midnight values in IDV introduces contamination by the substorm current wedge auroral electrojet, giving noise and a dependence on solar wind speed that varies with latitude. The composite is compared with independent, early data from European-sector stations, Greenwich, St Petersburg, Parc St Maur, and Ekaterinburg, as well as the composite u index, compiled from 2–6 stations by Bartels, and the IDV index of Svalgaard and Cliver. Agreement is found to be extremely good in all cases, except two. Firstly, the Greenwich data are shown to have gradually degraded in quality until new instrumentation was installed in 1915. Secondly, we infer that the Bartels u index is increasingly unreliable before about 1886 and overestimates the solar cycle amplitude between 1872 and 1883 and this is amplified in the proxy data used before 1872. This is therefore also true of the IDV index which makes direct use of the u index values.

A multi-proxy approach for revealing recent climatic changes in the Russian Altai

For the first time we present a multi-proxy data set for the Russian Altai, consisting of Siberian larch tree-ring width (TRW), latewood density (MXD), δ13C and δ18O in cellulose chronologies obtained for the period 1779–2007 and cell wall thickness (CWT) for 1900–2008. All of these parameters agree well between each other in the high-frequency variability, while the low-frequency climate information shows systematic differences. The correlation analysis with temperature and precipitation data from the closest weather station and gridded data revealed that annual TRW, MXD, CWT, and δ13C data contain a strong summer temperature signal, while δ18O in cellulose represents a mixed summer and winter temperature and precipitation signal. The temperature and precipitation reconstructions from the Belukha ice core and Teletskoe lake sediments were used to investigate the correspondence of different independent proxies. Low frequency patterns in TRW and δ13C chronologies are consistent with temperature reconstructions from nearby Belukha ice core and Teletskoe lake sediments showing a pronounced warming trend in the last century. Their combination could be used for the regional temperature reconstruction. The long-term δ18O trend agrees with the precipitation reconstruction from the Teletskoe lake sediment indicating more humid conditions during the twentieth century. Therefore, these two proxies could be combined for the precipitation reconstruction.