Tutorial 3: Reconstructing Past Changes in Terrestrial Climate#
Week 1, Day 4, Paleoclimate
Content creators: Sloane Garelick
Content reviewers: Yosmely Bermúdez, Dionessa Biton, Katrina Dobson, Maria Gonzalez, Will Gregory, Nahid Hasan, Paul Heubel, Sherry Mi, Beatriz Cosenza Muralles, Brodie Pearson, Jenna Pearson, Chi Zhang, Ohad Zivan
Content editors: Yosmely Bermúdez, Paul Heubel, Zahra Khodakaramimaghsoud, Jenna Pearson, Agustina Pesce, Chi Zhang, Ohad Zivan
Production editors: Wesley Banfield, Paul Heubel, Jenna Pearson, Konstantine Tsafatinos, Chi Zhang, Ohad Zivan
Our 2024 Sponsors: CMIP, NFDI4Earth
Tutorial Objectives#
Estimated timing of tutorial: 20 minutes
In this tutorial, we’ll explore the Euro2K proxy network, which is a subset of PAGES2K, the database we explored in the first tutorial. We will specifically focus on interpreting temperature change over the past 2,000 years as recorded by proxy records from tree rings, speleothems, and lake sediments. To analyze these datasets, we will group them by archive and create time series plots to assess temperature variations.
During this tutorial, you will:
Plot temperature records based on three different terrestrial proxies
Assess similarities and differences between the temperature records
Setup#
# installations ( uncomment and run this cell ONLY when using google colab or kaggle )
# !pip install LiPD
# !pip install cartopy
# !pip install pyleoclim==0.14.0
# imports
import pyleoclim as pyleo
import pandas as pd
import numpy as np
import os
import pooch
import tempfile
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from io import StringIO
import sys
Install and import feedback gadget#
Show code cell source
# @title Install and import feedback gadget
!pip3 install vibecheck datatops --quiet
from vibecheck import DatatopsContentReviewContainer
def content_review(notebook_section: str):
return DatatopsContentReviewContainer(
"", # No text prompt
notebook_section,
{
"url": "https://pmyvdlilci.execute-api.us-east-1.amazonaws.com/klab",
"name": "comptools_4clim",
"user_key": "l5jpxuee",
},
).render()
feedback_prefix = "W1D4_T3"
[notice] A new release of pip is available: 24.2 -> 24.3.1
[notice] To update, run: pip install --upgrade pip
Figure Settings#
Show code cell source
# @title Figure Settings
import ipywidgets as widgets # interactive display
%config InlineBackend.figure_format = 'retina'
plt.style.use(
"https://raw.githubusercontent.com/neuromatch/climate-course-content/main/cma.mplstyle"
)
Helper functions#
Show code cell source
# @title Helper functions
def pooch_load(filelocation=None, filename=None, processor=None):
# this is different for each day
shared_location = "/home/jovyan/shared/Data/tutorials/W1D4_Paleoclimate"
user_temp_cache = tempfile.gettempdir()
if os.path.exists(os.path.join(shared_location, filename)):
file = os.path.join(shared_location, filename)
else:
file = pooch.retrieve(
filelocation,
known_hash=None,
fname=os.path.join(user_temp_cache, filename),
processor=processor,
)
return file
class SupressOutputs(list):
def __enter__(self):
self._stdout = sys.stdout
sys.stdout = self._stringio = StringIO()
return self
def __exit__(self, *args):
self.extend(self._stringio.getvalue().splitlines())
del self._stringio # free up some memory
sys.stdout = self._stdout
Video 1: Terrestrial Climate Proxies#
Submit your feedback#
Show code cell source
# @title Submit your feedback
content_review(f"{feedback_prefix}_Terrestrial_Climate_Proxies_Video")
If you want to download the slides: https://osf.io/download/6gcsh/
Submit your feedback#
Show code cell source
# @title Submit your feedback
content_review(f"{feedback_prefix}_Terrestrial_Climate_Proxies_Slides")
Section 1: Loading Terrestrial Paleoclimate Records#
First, we need to download the data. Similar to Tutorial 1, the data is stored as a LiPD file, and we will be using Pyleoclim to format and interpret the data.
# set the name to save the Euro2k data
fname = "euro2k_data"
# download the data
lipd_file_path = pooch.retrieve(
url="https://osf.io/7ezp3/download/",
known_hash=None,
path="./",
fname=fname,
processor=pooch.Unzip(),
)
# the LiPD object can be used to load datasets stored in the LiPD format.
# in this first case study, we will load an entire library of LiPD files:
with SupressOutputs():
d_euro = pyleo.Lipd(os.path.join(".", f"{fname}.unzip", "Euro2k"))
/tmp/ipykernel_118657/2313300970.py:4: DeprecationWarning: The Lipd class is being deprecated and will be removed in Pyleoclim v1.0.0. Functionalities will instead be handled by the pyLipd package.
d_euro = pyleo.Lipd(os.path.join(".", f"{fname}.unzip", "Euro2k"))
Section 2: Temperature Reconstructions#
Before plotting, let’s narrow the data down a bit. We can filter all of the data so that we only keep reconstructions of temperature from terrestrial archives (e.g. tree rings, speleothems and lake sediments). This is accomplished with the function below.
def filter_data(dataset, archive_type, variable_name):
"""
Return a MultipleSeries object with the variable record (variable_name) for a given archive_type and coordinates.
"""
# Create a list of dictionary that can be iterated upon using Lipd.to_tso method
ts_list = dataset.to_tso()
# Append the correct indices for a given value of archive_type and variable_name
indices = []
lat = []
lon = []
for idx, item in enumerate(ts_list):
# Check that it is available to avoid errors on the loop
if "archiveType" in item.keys():
# If it's a archive_type, then proceed to the next step
if item["archiveType"] == archive_type:
if item["paleoData_variableName"] == variable_name:
indices.append(idx)
print(indices)
# Create a list of LipdSeries for the given indices
ts_list_archive_type = []
for indice in indices:
ts_list_archive_type.append(pyleo.LipdSeries(ts_list[indice]))
# save lat and lons of proxies
lat.append(ts_list[indice]["geo_meanLat"])
lon.append(ts_list[indice]["geo_meanLon"])
return pyleo.MultipleSeries(ts_list_archive_type), lat, lon
In the function above, the Lipd.to_tso() method is used to obtain a list of dictionaries that can be iterated upon.
ts_list = d_euro.to_tso()
extracting paleoData...
extracting: Ocn-RedSea.Felis.2000
extracting: Eur-CoastofPortugal.Abrantes.2011
extracting: Eur-MaritimeFrenchAlps.Büntgen.2012
extracting: Eur-LakeSilvaplana.Larocque-Tobler.2010
extracting: Eur-EuropeanAlps.Büntgen.2011
extracting: Eur-CentralandEasternPyrenees.Pla.2004
extracting: Eur-TatraMountains.Büntgen.2013
extracting: Arc-Jamtland.Wilson.2016
extracting: Ocn-AqabaJordanAQ18.Heiss.1999
extracting: Eur-EasternCarpathianMountains.Popa.2008
extracting: Eur-LakeSilvaplana.Trachsel.2010
extracting: Eur-Lötschental.Büntgen.2006
extracting: Eur-Tallinn.Tarand.2001
extracting: Arc-Forfjorddalen.McCarroll.2013
extracting: Ocn-AqabaJordanAQ19.Heiss.1999
extracting: Arc-GulfofAlaska.Wilson.2014
extracting: Eur-NorthernSpain.Martín-Chivelet.2011
extracting: Arc-Tjeggelvas.Bjorklund.2012
extracting: Arc-AkademiiNaukIceCap.Opel.2013
extracting: Arc-PolarUrals.Wilson.2015
extracting: Eur-Stockholm.Leijonhufvud.2009
extracting: Eur-NorthernScandinavia.Esper.2012
extracting: Eur-SpannagelCave.Mangini.2005
extracting: Arc-Kittelfjall.Bjorklund.2012
extracting: Eur-RAPiD-17-5P.Moffa-Sanchez.2014
extracting: Eur-Seebergsee.Larocque-Tobler.2012
extracting: Arc-Indigirka.Hughes.1999
extracting: Eur-CentralEurope.Dobrovolný.2009
extracting: Eur-SpanishPyrenees.Dorado-Linan.2012
extracting: Eur-FinnishLakelands.Helama.2014
extracting: Arc-Tornetrask.Melvin.2012
Created time series: 73 entries
Dictionaries are native to Python and can be explored as shown below.
# look at available entries for just one time-series
ts_list[0].keys()
dict_keys(['mode', 'time_id', '@context', 'archiveType', 'createdBy', 'dataContributor', 'dataSetName', 'googleDataURL', 'googleMetadataWorksheet', 'googleSpreadSheetKey', 'originalDataURL', 'studyName', 'tagMD5', 'pub1_author', 'pub1_citeKey', 'pub1_dataUrl', 'pub1_issue', 'pub1_journal', 'pub1_pages', 'pub1_publisher', 'pub1_title', 'pub1_type', 'pub1_volume', 'pub1_year', 'pub1_doi', 'pub2_author', 'pub2_citeKey', 'pub2_dataUrl', 'pub2_issue', 'pub2_journal', 'pub2_pages', 'pub2_publisher', 'pub2_title', 'pub2_type', 'pub2_volume', 'pub2_year', 'pub2_doi', 'pub3_author', 'pub3_Urldate', 'pub3_citeKey', 'pub3_institution', 'pub3_title', 'pub3_type', 'pub3_url', 'geo_type', 'geo_meanLon', 'geo_meanLat', 'geo_meanElev', 'geo_ocean', 'geo_pages2kRegion', 'geo_siteName', 'lipdVersion', 'tableType', 'paleoData_paleoDataTableName', 'paleoData_paleoDataMD5', 'paleoData_googleWorkSheetKey', 'paleoData_measurementTableName', 'paleoData_measurementTableMD5', 'paleoData_filename', 'paleoData_tableName', 'paleoData_missingValue', 'year', 'yearUnits', 'paleoData_QCCertification', 'paleoData_TSid', 'paleoData_WDSPaleoUrl', 'paleoData_archiveType', 'paleoData_hasMaxValue', 'paleoData_hasMeanValue', 'paleoData_hasMedianValue', 'paleoData_hasMinValue', 'paleoData_inCompilation', 'paleoData_iso2kUI', 'paleoData_notes', 'paleoData_ocean2kID', 'paleoData_pages2kID', 'paleoData_paleoMeasurementTableMD5', 'paleoData_proxy', 'paleoData_proxyObservationType', 'paleoData_sensorGenus', 'paleoData_sensorSpecies', 'paleoData_units', 'paleoData_useInGlobalTemperatureAnalysis', 'paleoData_variableName', 'paleoData_variableType', 'paleoData_hasResolution_hasMinValue', 'paleoData_hasResolution_hasMaxValue', 'paleoData_hasResolution_hasMeanValue', 'paleoData_hasResolution_hasMedianValue', 'paleoData_interpretation', 'paleoData_number', 'paleoData_values'])
# print relevant information for all entries
for idx, item in enumerate(ts_list):
print(str(idx) + ": " + item["dataSetName"] +
": " + item["paleoData_variableName"])
0: Ocn-RedSea.Felis.2000: d18O
1: Ocn-RedSea.Felis.2000: year
2: Eur-CoastofPortugal.Abrantes.2011: temperature
3: Eur-CoastofPortugal.Abrantes.2011: year
4: Eur-MaritimeFrenchAlps.Büntgen.2012: trsgi
5: Eur-MaritimeFrenchAlps.Büntgen.2012: year
6: Eur-LakeSilvaplana.Larocque-Tobler.2010: temperature
7: Eur-LakeSilvaplana.Larocque-Tobler.2010: year
8: Eur-EuropeanAlps.Büntgen.2011: trsgi
9: Eur-EuropeanAlps.Büntgen.2011: year
10: Eur-CentralandEasternPyrenees.Pla.2004: sampleID
11: Eur-CentralandEasternPyrenees.Pla.2004: year
12: Eur-CentralandEasternPyrenees.Pla.2004: age
13: Eur-CentralandEasternPyrenees.Pla.2004: temperature
14: Eur-CentralandEasternPyrenees.Pla.2004: uncertainty_temperature
15: Eur-TatraMountains.Büntgen.2013: trsgi
16: Eur-TatraMountains.Büntgen.2013: year
17: Arc-Jamtland.Wilson.2016: MXD
18: Arc-Jamtland.Wilson.2016: year
19: Ocn-AqabaJordanAQ18.Heiss.1999: d18O
20: Ocn-AqabaJordanAQ18.Heiss.1999: year
21: Ocn-AqabaJordanAQ18.Heiss.1999: d13C
22: Eur-EasternCarpathianMountains.Popa.2008: trsgi
23: Eur-EasternCarpathianMountains.Popa.2008: year
24: Eur-LakeSilvaplana.Trachsel.2010: temperature
25: Eur-LakeSilvaplana.Trachsel.2010: year
26: Eur-Lötschental.Büntgen.2006: MXD
27: Eur-Lötschental.Büntgen.2006: year
28: Eur-Tallinn.Tarand.2001: temperature
29: Eur-Tallinn.Tarand.2001: year
30: Eur-Tallinn.Tarand.2001: JulianDay
31: Arc-Forfjorddalen.McCarroll.2013: MXD
32: Arc-Forfjorddalen.McCarroll.2013: year
33: Ocn-AqabaJordanAQ19.Heiss.1999: d18O
34: Ocn-AqabaJordanAQ19.Heiss.1999: year
35: Ocn-AqabaJordanAQ19.Heiss.1999: d13C
36: Arc-GulfofAlaska.Wilson.2014: temperature
37: Arc-GulfofAlaska.Wilson.2014: year
38: Eur-NorthernSpain.Martín-Chivelet.2011: d18O
39: Eur-NorthernSpain.Martín-Chivelet.2011: year
40: Arc-Tjeggelvas.Bjorklund.2012: density
41: Arc-Tjeggelvas.Bjorklund.2012: year
42: Arc-AkademiiNaukIceCap.Opel.2013: thickness
43: Arc-AkademiiNaukIceCap.Opel.2013: year
44: Arc-AkademiiNaukIceCap.Opel.2013: age
45: Arc-AkademiiNaukIceCap.Opel.2013: d18O
46: Arc-AkademiiNaukIceCap.Opel.2013: Na
47: Arc-PolarUrals.Wilson.2015: density
48: Arc-PolarUrals.Wilson.2015: year
49: Eur-Stockholm.Leijonhufvud.2009: temperature
50: Eur-Stockholm.Leijonhufvud.2009: year
51: Eur-NorthernScandinavia.Esper.2012: MXD
52: Eur-NorthernScandinavia.Esper.2012: year
53: Eur-SpannagelCave.Mangini.2005: d18O
54: Eur-SpannagelCave.Mangini.2005: year
55: Arc-Kittelfjall.Bjorklund.2012: density
56: Arc-Kittelfjall.Bjorklund.2012: year
57: Eur-RAPiD-17-5P.Moffa-Sanchez.2014: d18O
58: Eur-RAPiD-17-5P.Moffa-Sanchez.2014: year
59: Eur-Seebergsee.Larocque-Tobler.2012: temperature
60: Eur-Seebergsee.Larocque-Tobler.2012: year
61: Arc-Indigirka.Hughes.1999: trsgi
62: Arc-Indigirka.Hughes.1999: year
63: Eur-CentralEurope.Dobrovolný.2009: temperature
64: Eur-CentralEurope.Dobrovolný.2009: year
65: Eur-SpanishPyrenees.Dorado-Linan.2012: trsgi
66: Eur-SpanishPyrenees.Dorado-Linan.2012: year
67: Eur-FinnishLakelands.Helama.2014: temperature
68: Eur-FinnishLakelands.Helama.2014: year
69: Arc-Tornetrask.Melvin.2012: temperature
70: Arc-Tornetrask.Melvin.2012: year
71: Arc-Tornetrask.Melvin.2012: temperature
72: Arc-Tornetrask.Melvin.2012: sampleDensity
Now let’s use our pre-defined function to create a new list that only has temperature reconstructions based on proxies from lake sediments:
ms_euro_lake, euro_lake_lat, euro_lake_lon = filter_data(
d_euro, "lake sediment", "temperature"
)
extracting paleoData...
extracting: Ocn-RedSea.Felis.2000
extracting: Eur-CoastofPortugal.Abrantes.2011
extracting: Eur-MaritimeFrenchAlps.Büntgen.2012
extracting: Eur-LakeSilvaplana.Larocque-Tobler.2010
extracting: Eur-EuropeanAlps.Büntgen.2011
extracting: Eur-CentralandEasternPyrenees.Pla.2004
extracting: Eur-TatraMountains.Büntgen.2013
extracting: Arc-Jamtland.Wilson.2016
extracting: Ocn-AqabaJordanAQ18.Heiss.1999
extracting: Eur-EasternCarpathianMountains.Popa.2008
extracting: Eur-LakeSilvaplana.Trachsel.2010
extracting: Eur-Lötschental.Büntgen.2006
extracting: Eur-Tallinn.Tarand.2001
extracting: Arc-Forfjorddalen.McCarroll.2013
extracting: Ocn-AqabaJordanAQ19.Heiss.1999
extracting: Arc-GulfofAlaska.Wilson.2014
extracting: Eur-NorthernSpain.Martín-Chivelet.2011
extracting: Arc-Tjeggelvas.Bjorklund.2012
extracting: Arc-AkademiiNaukIceCap.Opel.2013
extracting: Arc-PolarUrals.Wilson.2015
extracting: Eur-Stockholm.Leijonhufvud.2009
extracting: Eur-NorthernScandinavia.Esper.2012
extracting: Eur-SpannagelCave.Mangini.2005
extracting: Arc-Kittelfjall.Bjorklund.2012
extracting: Eur-RAPiD-17-5P.Moffa-Sanchez.2014
extracting: Eur-Seebergsee.Larocque-Tobler.2012
extracting: Arc-Indigirka.Hughes.1999
extracting: Eur-CentralEurope.Dobrovolný.2009
extracting: Eur-SpanishPyrenees.Dorado-Linan.2012
extracting: Eur-FinnishLakelands.Helama.2014
extracting: Arc-Tornetrask.Melvin.2012
Created time series: 73 entries
[6, 13, 24, 59]
Both age and year information are available, using age
/tmp/ipykernel_118657/2522667653.py:22: DeprecationWarning: The LipdSeries class is being deprecated and will be removed in Pyleoclim v1.0.0. It will be replaced by the geoSeries class (currently in development).
ts_list_archive_type.append(pyleo.LipdSeries(ts_list[indice]))
And a new list that only has temperature reconstructions based on proxies from tree rings:
ms_euro_tree, euro_tree_lat, euro_tree_lon = filter_data(
d_euro, "tree", "temperature")
extracting paleoData...
extracting: Ocn-RedSea.Felis.2000
extracting: Eur-CoastofPortugal.Abrantes.2011
extracting: Eur-MaritimeFrenchAlps.Büntgen.2012
extracting: Eur-LakeSilvaplana.Larocque-Tobler.2010
extracting: Eur-EuropeanAlps.Büntgen.2011
extracting: Eur-CentralandEasternPyrenees.Pla.2004
extracting: Eur-TatraMountains.Büntgen.2013
extracting: Arc-Jamtland.Wilson.2016
extracting: Ocn-AqabaJordanAQ18.Heiss.1999
extracting: Eur-EasternCarpathianMountains.Popa.2008
extracting: Eur-LakeSilvaplana.Trachsel.2010
extracting: Eur-Lötschental.Büntgen.2006
extracting: Eur-Tallinn.Tarand.2001
extracting: Arc-Forfjorddalen.McCarroll.2013
extracting: Ocn-AqabaJordanAQ19.Heiss.1999
extracting: Arc-GulfofAlaska.Wilson.2014
extracting: Eur-NorthernSpain.Martín-Chivelet.2011
extracting: Arc-Tjeggelvas.Bjorklund.2012
extracting: Arc-AkademiiNaukIceCap.Opel.2013
extracting: Arc-PolarUrals.Wilson.2015
extracting: Eur-Stockholm.Leijonhufvud.2009
extracting: Eur-NorthernScandinavia.Esper.2012
extracting: Eur-SpannagelCave.Mangini.2005
extracting: Arc-Kittelfjall.Bjorklund.2012
extracting: Eur-RAPiD-17-5P.Moffa-Sanchez.2014
extracting: Eur-Seebergsee.Larocque-Tobler.2012
extracting: Arc-Indigirka.Hughes.1999
extracting: Eur-CentralEurope.Dobrovolný.2009
extracting: Eur-SpanishPyrenees.Dorado-Linan.2012
extracting: Eur-FinnishLakelands.Helama.2014
extracting: Arc-Tornetrask.Melvin.2012
Created time series: 73 entries
[36, 67, 69, 71]
/tmp/ipykernel_118657/2522667653.py:22: DeprecationWarning: The LipdSeries class is being deprecated and will be removed in Pyleoclim v1.0.0. It will be replaced by the geoSeries class (currently in development).
ts_list_archive_type.append(pyleo.LipdSeries(ts_list[indice]))
And finally, a new list that only has temperature information based on proxies from speleothems:
ms_euro_spel, euro_spel_lat, euro_spel_lon = filter_data(
d_euro, "speleothem", "d18O")
extracting paleoData...
extracting: Ocn-RedSea.Felis.2000
extracting: Eur-CoastofPortugal.Abrantes.2011
extracting: Eur-MaritimeFrenchAlps.Büntgen.2012
extracting: Eur-LakeSilvaplana.Larocque-Tobler.2010
extracting: Eur-EuropeanAlps.Büntgen.2011
extracting: Eur-CentralandEasternPyrenees.Pla.2004
extracting: Eur-TatraMountains.Büntgen.2013
extracting: Arc-Jamtland.Wilson.2016
extracting: Ocn-AqabaJordanAQ18.Heiss.1999
extracting: Eur-EasternCarpathianMountains.Popa.2008
extracting: Eur-LakeSilvaplana.Trachsel.2010
extracting: Eur-Lötschental.Büntgen.2006
extracting: Eur-Tallinn.Tarand.2001
extracting: Arc-Forfjorddalen.McCarroll.2013
extracting: Ocn-AqabaJordanAQ19.Heiss.1999
extracting: Arc-GulfofAlaska.Wilson.2014
extracting: Eur-NorthernSpain.Martín-Chivelet.2011
extracting: Arc-Tjeggelvas.Bjorklund.2012
extracting: Arc-AkademiiNaukIceCap.Opel.2013
extracting: Arc-PolarUrals.Wilson.2015
extracting: Eur-Stockholm.Leijonhufvud.2009
extracting: Eur-NorthernScandinavia.Esper.2012
extracting: Eur-SpannagelCave.Mangini.2005
extracting: Arc-Kittelfjall.Bjorklund.2012
extracting: Eur-RAPiD-17-5P.Moffa-Sanchez.2014
extracting: Eur-Seebergsee.Larocque-Tobler.2012
extracting: Arc-Indigirka.Hughes.1999
extracting: Eur-CentralEurope.Dobrovolný.2009
extracting: Eur-SpanishPyrenees.Dorado-Linan.2012
extracting: Eur-FinnishLakelands.Helama.2014
extracting: Arc-Tornetrask.Melvin.2012
Created time series: 73 entries
[38, 53]
/tmp/ipykernel_118657/2522667653.py:22: DeprecationWarning: The LipdSeries class is being deprecated and will be removed in Pyleoclim v1.0.0. It will be replaced by the geoSeries class (currently in development).
ts_list_archive_type.append(pyleo.LipdSeries(ts_list[indice]))
Coding Exercises 2#
Using the coordinate information output from the filter_data() function, make a plot of the locations of the proxies using the markers and colors from Tutorial 1. Note that data close together may not be very visible.
# initiate plot
fig = plt.figure()
# set base map projection
ax = plt.axes(projection=ccrs.Robinson())
ax.set_global()
# add land fratures using gray color
ax.add_feature(cfeature.LAND, facecolor="k")
# add coastlines
ax.add_feature(cfeature.COASTLINE)
#################################################
## TODO for students: add the proxy locations ##
# Remove or comment the following line of code once you have completed the exercise:
raise NotImplementedError("Student exercise: Add the proxy locations via a scatter plot method.")
#################################################
# add the proxy locations
_ = ax.scatter(
...,
...,
transform=ccrs.Geodetic(),
label=...,
s=50,
marker="d",
color=[0.52734375, 0.8046875, 0.97916667],
edgecolor="k",
zorder=2,
)
_ = ax.scatter(
...,
...,
transform=ccrs.Geodetic(),
label=...,
s=50,
marker="p",
color=[0.73828125, 0.71484375, 0.41796875],
edgecolor="k",
zorder=2,
)
_ = ax.scatter(
...,
...,
transform=ccrs.Geodetic(),
label=...,
s=50,
marker="8",
color=[1, 0, 0],
edgecolor="k",
zorder=2,
)
# change the map view to zoom in on central Pacific
ax.set_extent((0, 360, 0, 90), crs=ccrs.PlateCarree())
ax.legend(
scatterpoints=1,
bbox_to_anchor=(0, -0.4),
loc="lower left",
ncol=3,
fontsize=15,
)
Example output:
Submit your feedback#
Show code cell source
# @title Submit your feedback
content_review(f"{feedback_prefix}_Coding_Exercises_2")
Since we are going to compare temperature datasets based on different terrestrial climate archives (lake sediments, tree rings and speleothems), the quantitative values of the measurements in each record will differ (i.e., the lake sediment and tree ring data are temperature in degrees C, but the speleothem data is oxygen isotopes in per mille). Therefore, to more easily and accurately compare temperature between the records, it’s helpful to standardize the data as we did in Tutorial 2. The .standardize() function removes the estimated mean of the time series and divides by its estimated standard deviation.
# standardize the data
spel_stnd = ms_euro_spel.standardize()
lake_stnd = ms_euro_lake.standardize()
tree_stnd = ms_euro_tree.standardize()
Now we can use Pyleoclim functions to create three stacked plots of this data with lake sediment records on top, tree ring reconstructions in the middle and speleothem records on the bottom.
Note that the colors used for the time series in each plot are the default colors generated by the function, so the corresponding colors in each of the three plots are not relevant.
# note the x-axis is years before present,
# so reading from left to right corresponds to moving back in time
ax = lake_stnd.stackplot(
label_x_loc=1.7,
xlim=[0, 2000],
v_shift_factor=1,
figsize=[9, 5],
time_unit="yrs BP",
)
ax[0].suptitle("Lake Cores", y=1.2)
ax = tree_stnd.stackplot(
label_x_loc=1.7,
xlim=[0, 2000],
v_shift_factor=1,
figsize=[9, 5],
time_unit="yrs BP",
)
ax[0].suptitle("Tree Rings", y=1.2)
# recall d18O is a proxy for SST, and that more positive d18O means colder SST
ax = spel_stnd.stackplot(
label_x_loc=1.7,
xlim=[0, 2000],
v_shift_factor=1,
figsize=[9, 5],
time_unit="yrs BP",
)
ax[0].suptitle("Speleothems", y=1.2)
Text(0.5, 1.2, 'Speleothems')
/opt/hostedtoolcache/Python/3.9.20/x64/lib/python3.9/site-packages/IPython/core/events.py:82: UserWarning: This figure includes Axes that are not compatible with tight_layout, so results might be incorrect.
func(*args, **kwargs)
/opt/hostedtoolcache/Python/3.9.20/x64/lib/python3.9/site-packages/IPython/core/pylabtools.py:152: UserWarning: This figure includes Axes that are not compatible with tight_layout, so results might be incorrect.
fig.canvas.print_figure(bytes_io, **kw)
Questions 2#
Using the plots we just made (and recalling that all of these records are from Europe), let’s make some inferences about the temperature data over the past 2,000 years:
Recall that δ18O is a proxy for SST, and that more positive δ18O means colder SST. Do the temperature records based on a single proxy type record similar patterns?
Do the three proxy types collectively record similar patterns?
What might be causing the more frequent variations in temperature?
Submit your feedback#
Show code cell source
# @title Submit your feedback
content_review(f"{feedback_prefix}_Questions_2")
Summary#
In this tutorial, we explored how to use the Euro2k proxy network to investigate changes in temperature over the past 2,000 years from tree rings, speleothems, and lake sediments. To analyze these diverse datasets, we categorized them based on their archive type and constructed time series plots.
Resources#
Code for this tutorial is based on an existing notebook from LinkedEarth that provides instruction on working with LiPD files.
Data from the following sources are used in this tutorial:
Euro2k database: PAGES2k Consortium., Emile-Geay, J., McKay, N. et al. A global multiproxy database for temperature reconstructions of the Common Era. Sci Data 4, 170088 (2017). https://doi.org/10.1038/sdata.2017.88