Reconstructions of global temperature of the past 2000 years, using composite of different proxy methods

In the study of past climates ("paleoclimatology"), climate proxies are preserved physical characteristics of the past that stand in for direct meteorological measurements[1] and enable scientists to reconstruct the climatic conditions over a longer fraction of the Earth's history. Reliable global records of climate only began in the 1880s, and proxies provide the only means for scientists to determine climatic patterns before record-keeping began.

A large number of climate proxies have been studied from a variety of geologic contexts. Examples of proxies include stable isotope measurements from ice cores, growth rates in tree rings, species composition of sub-fossil pollen in lake sediment or foraminifera in ocean sediments, temperature profiles of boreholes, and stable isotopes and mineralogy of corals and carbonate speleothems. In each case, the proxy indicator has been influenced by a particular seasonal climate parameter (e.g., summer temperature or monsoon intensity) at the time in which they were laid down or grew. Interpretation of climate proxies requires a range of ancillary studies, including calibration of the sensitivity of the proxy to climate and cross-verification among proxy indicators.[2]

Proxies can be combined to produce temperature reconstructions longer than the instrumental temperature record and can inform discussions of global warming and climate history. The geographic distribution of proxy records, just like the instrumental record, is not at all uniform, with more records in the northern hemisphere.[3]


Main article: Proxy (statistics)

In science, it is sometimes necessary to study a variable which cannot be measured directly. This can be done by "proxy methods," in which a variable which correlates with the variable of interest is measured, and then used to infer the value of the variable of interest. Proxy methods are of particular use in the study of the past climate, beyond times when direct measurements of temperatures are available.

Most proxy records have to be calibrated against independent temperature measurements, or against a more directly calibrated proxy, during their period of overlap to estimate the relationship between temperature and the proxy. The longer history of the proxy is then used to reconstruct temperature from earlier periods.

Ice cores


Ice Core sample taken from drill. Photo by Lonnie Thompson, Byrd Polar Research Center.

Ice cores are cylindrical samples from within ice sheets in the Greenland, Antarctic, and North American regions.[4][5] First attempts of extraction occurred in 1956 as part of the International Geophysical Year. As original means of extraction, the U.S. Army's Cold Regions Research and Engineering Laboratory used an 80-foot (24 m)-long modified electrodrill in 1968 at Camp Century, Greenland, and Byrd Station, Antarctica. Their machinery could drill through 15–20 feet (4.6–6.1 m) of ice in 40–50 minutes. From 1300 to 3,000 feet (910 m) in depth, core samples were 4+14 inches (110 mm) in diameter and 10 to 20 feet (6.1 m) long. Deeper samples of 15 to 20 feet (6.1 m) long were not uncommon. Every subsequent drilling team improves their method with each new effort.[6]


δ18Oair and δDice for Vostok, Antarctica ice core.

The ratio between the 16O and 18O water molecule isotopologues in an ice core helps determine past temperatures and snow accumulations.[4] The heavier isotope (18O) condenses more readily as temperatures decrease and falls more easily as precipitation, while the lighter isotope (16O) needs colder conditions to precipitate. The farther north one needs to go to find elevated levels of the 18O isotopologue, the warmer the period.[further explanation needed][7]

In addition to oxygen isotopes, water contains hydrogen isotopes – 1H and 2H, usually referred to as H and D (for deuterium) – that are also used for temperature proxies. Normally, ice cores from Greenland are analyzed for δ18O and those from Antarctica for δ-deuterium.[why?] Those cores that analyze for both show a lack of agreement.[citation needed] (In the figure, δ18O is for the trapped air, not the ice. δD is for the ice.)

Air bubbles in the ice, which contain trapped greenhouse gases such as carbon dioxide and methane, are also helpful in determining past climate changes.[4]

From 1989 to 1992, the European Greenland Ice Core Drilling Project drilled in central Greenland at coordinates 72° 35' N, 37° 38' W. The ices in that core were 3840 years old at a depth of 770 m, 40,000 years old at 2521 m, and 200,000 years old or more at 3029 m bedrock.[8] Ice cores in Antarctica can reveal the climate records for the past 650,000 years.[4]

Location maps and a complete list of U.S. ice core drilling sites can be found on the website for the National Ice Core Laboratory.[5]

Tree rings

Tree rings seen in a cross section of a trunk of a tree.

Main article: Dendroclimatology

Dendroclimatology is the science of determining past climates from trees, primarily from properties of the annual tree rings. Tree rings are wider when conditions favor growth, narrower when times are difficult. Two primary factors are temperature and humidity / water availability. Other properties of the annual rings, such as maximum latewood density (MXD) have been shown to be better proxies than simple ring width. Using tree rings, scientists have estimated many local climates for hundreds to thousands of years previous. By combining multiple tree-ring studies (sometimes with other climate proxy records), scientists have estimated past regional and global climates (see Temperature record of the past 1000 years).

Fossil leaves

Paleoclimatologists often use leaf teeth to reconstruct mean annual temperature in past climates, and they use leaf size as a proxy for mean annual precipitation.[9] In the case of mean annual precipitation reconstructions, some researchers believe taphonomic processes cause smaller leaves to be overrepresented in the fossil record, which can bias reconstructions.[10] However, recent research suggests that the leaf fossil record may not be significantly biased toward small leaves.[11] New approaches retrieve data such as CO2 content of past atmospheres from fossil leaf stomata and isotope composition, measuring cellular CO2 concentrations. A 2014 study was able to use the carbon-13 isotope ratios to estimate the CO2 amounts of the past 400 million years, the findings hint at a higher climate sensitivity to CO2 concentrations.[12]


Borehole temperatures are used as temperature proxies. Since heat transfer through the ground is slow, temperature measurements at a series of different depths down the borehole, adjusted for the effect of rising heat from inside the Earth, can be "inverted" (a mathematical formula to solve matrix equations) to produce a non-unique series of surface temperature values. The solution is "non-unique" because there are multiple possible surface temperature reconstructions that can produce the same borehole temperature profile. In addition, due to physical limitations, the reconstructions are inevitably "smeared", and become more smeared further back in time. When reconstructing temperatures around 1500 AD, boreholes have a temporal resolution of a few centuries. At the start of the 20th century, their resolution is a few decades; hence they do not provide a useful check on the instrumental temperature record.[13][14] However, they are broadly comparable.[3] These confirmations have given paleoclimatologists the confidence that they can measure the temperature of 500 years ago. This is concluded by a depth scale of about 492 feet (150 meters) to measure the temperatures from 100 years ago and 1,640 feet (500 meters) to measure the temperatures from 1,000 years ago.[15]

Boreholes have a great advantage over many other proxies in that no calibration is required: they are actual temperatures. However, they record surface temperature not the near-surface temperature (1.5 meter) used for most "surface" weather observations. These can differ substantially under extreme conditions or when there is surface snow. In practice the effect on borehole temperature is believed to be generally small. A second source of error is contamination of the well by groundwater may affect the temperatures, since the water "carries" more modern temperatures with it. This effect is believed to be generally small, and more applicable at very humid sites.[13] It does not apply in ice cores where the site remains frozen all year.

More than 600 boreholes, on all continents, have been used as proxies for reconstructing surface temperatures.[14] The highest concentration of boreholes exist in North America and Europe. Their depths of drilling typically range from 200 to greater than 1,000 meters into the crust of the Earth or ice sheet.[15]

A small number of boreholes have been drilled in the ice sheets; the purity of the ice there permits longer reconstructions. Central Greenland borehole temperatures show "a warming over the last 150 years of approximately 1°C ± 0.2°C preceded by a few centuries of cool conditions. Preceding this was a warm period centered around A.D. 1000, which was warmer than the late 20th century by approximately 1°C." A borehole in the Antarctica icecap shows that the "temperature at A.D. 1 [was] approximately 1°C warmer than the late 20th century".[16]

Borehole temperatures in Greenland were responsible for an important revision to the isotopic temperature reconstruction, revealing that the former assumption that "spatial slope equals temporal slope" was incorrect.


Coral bleached due to changes in ocean water properties

Ocean coral skeletal rings, or bands, also share paleoclimatological information, similarly to tree rings. In 2002, a report was published on the findings of Drs. Lisa Greer and Peter Swart, associates of University of Miami at the time, in regard to stable oxygen isotopes in the calcium carbonate of coral. Cooler temperatures tend to cause coral to use heavier isotopes in its structure, while warmer temperatures result in more normal oxygen isotopes being built into the coral structure. Denser water salinity also tends to contain the heavier isotope. Greer's coral sample from the Atlantic Ocean was taken in 1994 and dated back to 1935. Greer recalls her conclusions, "When we look at the averaged annual data from 1935 to about 1994, we see it has the shape of a sine wave. It is periodic and has a significant pattern of oxygen isotope composition that has a peak at about every twelve to fifteen years." Surface water temperatures have coincided by also peaking every twelve and a half years. However, since recording this temperature has only been practiced for the last fifty years, correlation between recorded water temperature and coral structure can only be drawn so far back.[17]

Pollen grains

Pollen can be found in sediments. Plants produce pollen in large quantities and it is extremely resistant to decay. It is possible to identify a plant species from its pollen grain. The identified plant community of the area at the relative time from that sediment layer, will provide information about the climatic condition. The abundance of pollen of a given vegetation period or year depends partly on the weather conditions of the previous months, hence pollen density provides information on short-term climatic conditions.[18] The study of prehistoric pollen is palynology.

Dinoflagellate cysts

Cyst of a dinoflagellate Peridinium ovatum

Dinoflagellates occur in most aquatic environments and during their life cycle, some species produce highly resistant organic-walled cysts for a dormancy period when environmental conditions are not appropriate for growth. Their living depth is relatively shallow (dependent upon light penetration), and closely coupled to diatoms on which they feed. Their distribution patterns in surface waters are closely related to physical characteristics of the water bodies, and nearshore assemblages can also be distinguished from oceanic assemblages. The distribution of dinocysts in sediments has been relatively well documented and has contributed to understanding the average sea-surface conditions that determine the distribution pattern and abundances of the taxa ([19]). Several studies, including [20] and [21] have compiled box and gravity cores in the North Pacific analyzing them for palynological content to determine the distribution of dinocysts and their relationships with sea surface temperature, salinity, productivity and upwelling. Similarly,[22] and [23] use a box core at 576.5 m of water depth from 1992 in the central Santa Barbara Basin to determine oceanographic and climatic changes during the past 40 kyr in the area.

Lake and ocean sediments

Similar to their study on other proxies, paleoclimatologists examine oxygen isotopes in the contents of ocean sediments. Likewise, they measure the layers of varve (deposited fine and coarse silt or clay)[24] laminating lake sediments. Lake varves are primarily influenced by:

Diatoms, foraminifera, radiolarians, ostracods, and coccolithophores are examples of biotic proxies for lake and ocean conditions that are commonly used to reconstruct past climates. The distribution of the species of these and other aquatic creatures preserved in the sediments are useful proxies. The optimal conditions for species preserved in the sediment act as clues. Researchers use these clues to reveal what the climate and environment was like when the creatures died.[26] The oxygen isotope ratios in their shells can also be used as proxies for temperature.[27]

Water isotopes and temperature reconstruction

Ocean water is mostly H216O, with small amounts of HD16O and H218O, where D denotes deuterium, i.e. hydrogen with an extra neutron. In Vienna Standard Mean Ocean Water (VSMOW) the ratio of D to H is 155.76x10−6 and O-18 to O-16 is 2005.2x10−6. Isotope fractionation occurs during changes between condensed and vapour phases: the vapour pressure of heavier isotopes is lower, so vapour contains relatively more of the lighter isotopes and when the vapour condenses the precipitation preferentially contains heavier isotopes. The difference from VSMOW is expressed as δ18O = 1000‰ ; and a similar formula for δD. δ values for precipitation are always negative.[28] The major influence on δ is the difference between ocean temperatures where the moisture evaporated and the place where the final precipitation occurred; since ocean temperatures are relatively stable the δ value mostly reflects the temperature where precipitation occurs. Taking into account that the precipitation forms above the inversion layer, we are left with a linear relation:

δ 18O = aT + b

This is empirically calibrated from measurements of temperature and δ as a = 0.67 /°C for Greenland and 0.76 ‰/°C for East Antarctica. The calibration was initially done on the basis of spatial variations in temperature and it was assumed that this corresponded to temporal variations.[29] More recently, borehole thermometry has shown that for glacial-interglacial variations, a = 0.33 ‰/°C,[30] implying that glacial-interglacial temperature changes were twice as large as previously believed.

A study published in 2017 called the previous methodology to reconstruct paleo ocean temperatures 100 million years ago into question, suggesting it has been relatively stable during that time, much colder.[31]

Membrane lipids

A novel climate proxy obtained from peat (lignites, ancient peat) and soils, membrane lipids known as glycerol dialkyl glycerol tetraether (GDGT) is helping to study paleo environmental factors, which control relative distribution of differently branched GDGT isomers. The study authors note, "These branched membrane lipids are produced by an as yet unknown group of anaerobic soil bacteria."[32] As of 2018, there is a decade of research demonstrating that in mineral soils the degree of methylation of bacteria (brGDGTs), helps to calculate mean annual air temperatures. This proxy method was used to study the climate of the early Palaeogene, at the Cretaceous–Paleogene boundary, and researchers found that annual air temperatures, over land and at mid-latitude, averaged about 23–29 °C (± 4.7 °C), which is 5–10 °C higher than most previous findings.[33][34]


The skill of algorithms used to combine proxy records into an overall hemispheric temperature reconstruction may be tested using a technique known as "pseudoproxies". In this method, output from a climate model is sampled at locations corresponding to the known proxy network, and the temperature record produced is compared to the (known) overall temperature of the model.[35]

See also


  1. ^ "What Are "Proxy" Data? | National Centers for Environmental Information (NCEI) formerly known as National Climatic Data Center (NCDC)". Archived from the original on 2020-03-08. Retrieved 2017-10-12.
  2. ^ "Climate Change 2001: Palaeoclimate proxy indicators." Archived 2009-12-04 at the Wayback Machine
  3. ^ a b "Borehole Temperatures Confirm Global Warming Pattern."
  4. ^ a b c d Strom, Robert. Hot House. p. 255
  5. ^ a b "Core Location Maps." Archived 2009-11-10 at the Wayback Machine
  6. ^ Vardiman, Larry, Ph.D. Ice Cores and the Age of the Earth. p. 9-13
  7. ^ "Paleoclimatology: the Oxygen Balance."
  8. ^ "The GRIP Coring Effort."
  9. ^ Dana L. Royer; Peter Wilf; David A. Janesko; Elizabeth A. Kowalski; David L. Dilcher (1 July 2005). "Correlations of climate and plant ecology to leaf size and shape: potential proxies for the fossil record". American Journal of Botany. 92 (7): 1141–1151. doi:10.3732/ajb.92.7.1141. PMID 21646136.
  10. ^ David R. Greenwood (1994), "Palaeobotanical evidence for Tertiary climates", History of the Australian Vegetation: Cretaceous to Recent: 44–59
  11. ^ Eric R. Hagen; Dana Royer; Ryan A. Moye; Kirk R. Johnson (9 January 2019). "No large bias within species between the reconstructed areas of complete and fragmented fossil leaves". PALAIOS. 34 (1): 43–48. Bibcode:2019Palai..34...43H. doi:10.2110/palo.2018.091. S2CID 133599753.
  12. ^ Peter J. Franks; Dana Royer; David J. Beerling; Peter K. Van de Water; David J. Cantrill; Margaret M. Barbour; Joseph A. Berry (16 July 2014). "New constraints on atmospheric CO2 concentration for the Phanerozoic" (PDF). Geophysical Research Letters. 31 (13): 4685–4694. Bibcode:2014GeoRL..41.4685F. doi:10.1002/2014GL060457. hdl:10211.3/200431. S2CID 55701037. Archived from the original (PDF) on 12 August 2014. Retrieved 31 July 2014.
  13. ^ a b Council, National Research; Studies, Division on Earth Life; Climate, Board on Atmospheric Sciences and; Committee On Surface Temperature Reconstructions For The Last 2, 000 Years (2006). Surface Temperature Reconstructions for the Last 2,000 Years. CiteSeerX doi:10.17226/11676. ISBN 978-0-309-10225-4.((cite book)): CS1 maint: numeric names: authors list (link)
  14. ^ a b Pollack, H. N.; Huang, S.; Shen, P. Y. (2000). "Temperature trends over the past five centuries reconstructed from borehole temperatures" (PDF). Nature. 403 (6771): 756–758. Bibcode:2000Natur.403..756H. doi:10.1038/35001556. hdl:2027.42/62610. PMID 10693801. S2CID 4425128.
  15. ^ a b Environmental News Network staff. "Borehole temperatures confirm global warming." Archived 2009-10-29 at the Wayback Machine
  16. ^ BOREHOLES IN GLACIAL ICE Surface Temperature Reconstructions for the Last 2,000 Years (2006), pp 81,82 Board on Atmospheric Sciences and Climate (BASC), National Academy of Science, ISBN 978-0-309-10225-4
  17. ^ "Coral Layers Good Proxy for Atlantic Climate Cycles." Archived 2010-03-16 at the Wayback Machine
  18. ^ Bradley, R. S. and Jones, P. D. (eds) 1992: Climate since AD 1500. London: Routledge.
  19. ^ de Vernal, A.; Eynaud, F.; Henry, M.; Hillaire-Marcel, C.; Londeix, L.; Mangin, S.; Matthiessen, J.; Marret, F.; Radi, T.; Rochon, A.; Solignac, S.; Turon, J. -L. (1 April 2005). "Reconstruction of sea-surface conditions at middle to high latitudes of the Northern Hemisphere during the Last Glacial Maximum (LGM) based on dinoflagellate cyst assemblages". Quaternary Science Reviews. 24 (7–9): 897–924. Bibcode:2005QSRv...24..897D. doi:10.1016/j.quascirev.2004.06.014.
  20. ^ Radi, Taoufik; de Vernal, Anne (1 January 2004). "Dinocyst distribution in surface sediments from the northeastern Pacific margin (40–60°N) in relation to hydrographic conditions, productivity and upwelling". Review of Palaeobotany and Palynology. 128 (1–2): 169–193. Bibcode:2004RPaPa.128..169R. doi:10.1016/S0034-6667(03)00118-0.
  21. ^ Pospelova, Vera; de Vernal, Anne; Pedersen, Thomas F. (1 July 2008). "Distribution of dinoflagellate cysts in surface sediments from the northeastern Pacific Ocean (43–25°N) in relation to sea-surface temperature, salinity, productivity and coastal upwelling". Marine Micropaleontology. 68 (1–2): 21–48. Bibcode:2008MarMP..68...21P. doi:10.1016/j.marmicro.2008.01.008.
  22. ^ Pospelova, Vera; Pedersen, Thomas F.; de Vernal, Anne (1 June 2006). "Dinoflagellate cysts as indicators of climatic and oceanographic changes during the past 40 kyr in the Santa Barbara Basin, southern California". Paleoceanography. 21 (2): PA2010. Bibcode:2006PalOc..21.2010P. doi:10.1029/2005PA001251. ISSN 1944-9186.
  23. ^ Bringué, Manuel; Pospelova, Vera; Field, David B. (1 December 2014). "High resolution sedimentary record of dinoflagellate cysts reflects decadal variability and 20th century warming in the Santa Barbara Basin". Quaternary Science Reviews. 105: 86–101. Bibcode:2014QSRv..105...86B. doi:10.1016/j.quascirev.2014.09.022.
  24. ^ "Varve."
  25. ^ "Climate Change 2001: Palaeoclimate proxy indicators" Archived 2009-12-04 at the Wayback Machine
  26. ^ Bruckner, Monica. "Paleoclimatology: How Can We Infer Past Climates?". Montana State University.
  27. ^ Shemesh, A.; Charles, C. D.; Fairbanks, R. G. (1992-06-05). "Oxygen Isotopes in Biogenic Silica: Global Changes in Ocean Temperature and Isotopic Composition". Science. 256 (5062): 1434–1436. Bibcode:1992Sci...256.1434S. doi:10.1126/science.256.5062.1434. ISSN 0036-8075. PMID 17791613. S2CID 38840484.
  28. ^ Council, National Research; Studies, Division on Earth Life; Climate, Board on Atmospheric Sciences and; Committee On Surface Temperature Reconstructions For The Last 2, 000 Years (2006). Surface Temperature Reconstructions for the Last 2,000 Years. CiteSeerX doi:10.17226/11676. ISBN 978-0-309-10225-4.((cite book)): CS1 maint: numeric names: authors list (link)
  29. ^ Jouzel and Merlivat, 1984) Deuterium and oxygen 18 in precipitation: Modeling of the isotopic effects during snow formation, Journal of Geophysical Research: Atmospheres, Volume 89, Issue D7, Pages 11589–11829
  30. ^ Cuffey et al., 1995, Large Arctic temperature change at the Wisconsin– Holocene glacial transition, Science 270: 455–458
  31. ^ Bernard, S.; Daval, D.; Ackerer, P.; Pont, S.; Meibom, A. (2017-10-26). "Burial-induced oxygen-isotope re-equilibration of fossil foraminifera explains ocean paleotemperature paradoxes". Nature Communications. 8 (1): 1134. Bibcode:2017NatCo...8.1134B. doi:10.1038/s41467-017-01225-9. ISSN 2041-1723. PMC 5656689. PMID 29070888.
  32. ^ Johan W.H. Weijers; Stefan Schouten; Jurgen C. van den Donker; Ellen C. Hopmans; Jaap S. Sinninghe Damsté (2007). "Environmental controls on bacterial tetraether membrane lipid distribution in soils". Geochimica et Cosmochimica Acta. 71 (3): 703–713. Bibcode:2007GeCoA..71..703W. doi:10.1016/j.gca.2006.10.003. S2CID 84677185.
  33. ^ B. D. A. Naafs; M. Rohrssen; G. N. Inglis; O. Lähteenoja; S. J. Feakins; M. E. Collinson; E. M. Kennedy; P. K. Singh; M. P. Singh; D. J. Lunt; R. D. Pancost (2018). "High temperatures in the terrestrial mid-latitudes during the early Palaeogene" (PDF). Nature Geoscience. 11 (10): 766–771. Bibcode:2018NatGe..11..766N. doi:10.1038/s41561-018-0199-0. hdl:1983/82e93473-2a5d-4a6d-9ca1-da5ebf433d8b. S2CID 135045515.
  34. ^ University of Bristol (30 July 2018). "Ever-increasing CO2 levels could take us back to the tropical climate of Paleogene period". ScienceDaily.
  35. ^ Mann, M. E.; Rutherford, S. (31 May 2002), "Climate reconstruction using 'Pseudoproxies'", Geophysical Research Letters, 29 (10): 139–1–139–4, Bibcode:2002GeoRL..29.1501M, doi:10.1029/2001GL014554

Further reading