File:Reuverian curve zachos 2001 raw data.PNG

From Wikimedia Commons, the free media repository
Jump to navigation Jump to search

Original file ‎(1,600 × 766 pixels, file size: 179 KB, MIME type: image/png)

Captions

Captions

Add a one-line explanation of what this file represents

Summary

[edit]
Description
English: Raw delta-O18 curve between 3,60-2,50 millios ago, that is Reuverian or Piacenzian stage. Based on James Zachos et al data file

2001 compliationdata.txt on web site

Zachos, J.C., Dickens, G.R., and Zeebe, R.E., 2008, An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature. |Vol 451|17 January 2008|doi:10.1038/nature

http://www.es.ucsc.edu/~jzachos/Publications.html. http://www.es.ucsc.edu/~jzachos/Data/2008CompilationData.xls

Y-axis is column 3 on file, that is is uncorrected delta_018. X axis is dating millions of years ago.

Python code to produce plot Note: You must preprocess data w/ spreadsheet and text editor

    3. drawing climate diagram in python 3
    4. version 2.11
    5. 11.9.2020


import matplotlib.pyplot as plt import numpy as np import pandas as pd from scipy import interpolate from matplotlib.ticker import (MultipleLocator, AutoMinorLocator) import scipy.signal

beginage=2.4 endage=3.8 beginy=1.8 endy=5.0

datafilename="compilation1.csv" captioni="Reuverian period - benthic D-O18" savename="lisiecki.svg"

figsizex=16 figsizey=8


size0=14 size1=16 size2=18 size3=24


dfin0=pd.read_csv(datafilename, sep=";")

  1. age_Ma;d18O_ajd;d13C

col1='age_Ma' col2='dO18_ajd'

lst1=[col1,col2]

dfin1 = dfin0[dfin0.columns.intersection(lst1)]

x0=dfin1[col1] y0=dfin1[col2]

x=np.array(x0) y=np.array(y0)


  1. y_savgol = scipy.signal.savgol_filter(y,31, 3)

y_savgol = scipy.signal.savgol_filter(y,31, 3)


  1. y_running = running_mean(y, 31)


  1. x_sm = np.array(x)
  2. y_sm = np.array(y)
  3. x_smooth = np.linspace(x_sm.min(), x_sm.max(), 20000)
  4. funk1 = interpolate.interp1d(x_sm, y_sm, kind="cubic")
  5. y_smooth = funk1(x_smooth)


fig, ax1 = plt.subplots()


  1. ax1.axis((11600,14000,0,ymax1))

ax1.set_xlim(beginage, endage) ax1.set_ylim(beginy, endy)

  1. ax1.set_ylim(-35.0, -42.0)

plt.gca().invert_xaxis() plt.gca().invert_yaxis()

ax1.set_ylabel('delta-O18', color='#0000ff', fontsize=size2+2)


ax1.plot(x,y, color="#005f00", linewidth=3,label="NGRIP delta-O18")

  1. ax1.plot(x_smooth,y_smooth, color="#0000ff", linewidth=3,label="NGRIP delta-O18")

ax1.plot(x,y_savgol, color="#a0FFa0", linewidth=1, label="SavGol filter, 21 and 3")

  1. ax1.plot(x,y_running, color="#FF0000", linewidth=3)


  1. ax1.plot(x,data_avg1, color="#ff0000", linewidth=2, linestyle=":", label="Average of NGRIP, GISP, GISP2 delta-O18")


ax1.text(3.44, 4.7, "Reuver A -->", fontsize=22, horizontalalignment='left') ax1.plot(3.44,4.9, marker="v", color="green",markersize=30) ax1.text(3.00, 4.7, "B1", fontsize=22, horizontalalignment='left') ax1.plot(3.00,4.9, marker="v", color="green",markersize=30) ax1.text(2.76, 4.7, "B2", fontsize=22, horizontalalignment='left') ax1.plot(2.76,4.9, marker="v", color="green",markersize=30) ax1.text(2.66, 4.7, "B3", fontsize=22, horizontalalignment='left') ax1.plot(2.66,4.9, marker="v", color="green",markersize=30) ax1.text(2.6, 4.7, "C1", fontsize=22, horizontalalignment='left') ax1.plot(2.6,4.9, marker="v", color="green",markersize=30) ax1.text(2.55, 4.7, "PT?", fontsize=22, horizontalalignment='left') ax1.plot(2.55,4.9, marker="v", color="green",markersize=30) ax1.text(2.44, 4.7, "PT?", fontsize=22, horizontalalignment='left') ax1.plot(2.44,4.9, marker="v", color="green",markersize=30)

ax1.plot(2.69,2.5, marker="o", color="blue",markersize=26) ax1.plot(2.52,2.5, marker="o", color="blue",markersize=26)

ax1.plot(2.85,2.5, marker="o", color="blue",markersize=26)

ax1.plot(3.125,2.5, marker="o", color="blue",markersize=26)

ax1.plot(3.3,2.5, marker="o", color="blue",markersize=26) ax1.plot(3.35,2.5, marker="o", color="blue",markersize=26)

    1. warm

ax1.plot(2.63,2.4, marker="o", color="red",markersize=26) ax1.plot(2.985,2.4, marker="o", color="red",markersize=26)

ax1.plot(3.4,2.4, marker="o", color="red",markersize=26)

ax1.tick_params(axis='both', which='major', labelsize=size2)

ax1.xaxis.set_minor_locator(MultipleLocator(1)) ax1.xaxis.set_minor_locator(MultipleLocator(0.05))

ax1.yaxis.set_minor_locator(MultipleLocator(1.0)) ax1.yaxis.set_minor_locator(MultipleLocator(0.1))

ax1.grid(which='major', linestyle='-', linewidth='0.1', color='black') ax1.grid(which='minor', linestyle=':', linewidth='0.1', color='black')


ax1.set_xlabel('Age ma BP', color="darkgreen", fontsize=size2+2)


ax1.set_title(captioni, fontsize=size3, color="#0000af")

plt.legend(fontsize=size0)

fig = plt.gcf() fig.set_size_inches(figsizex, figsizey, forward=True)


plt.savefig(savename, format="svg", dpi = 100)

plt.show()


Date 3 August 2008 (upload date)
Source Own work
Author Merikanto

Licensing

[edit]
Public domain I, the copyright holder of this work, release this work into the public domain. This applies worldwide.
In some countries this may not be legally possible; if so:
I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law.

File history

Click on a date/time to view the file as it appeared at that time.

Date/TimeThumbnailDimensionsUserComment
current18:03, 13 September 2020Thumbnail for version as of 18:03, 13 September 20201,600 × 766 (179 KB)Merikanto (talk | contribs)New layout
18:29, 3 August 2008Thumbnail for version as of 18:29, 3 August 2008840 × 350 (21 KB)Merikanto~commonswiki (talk | contribs)reuverian age paleoclimatology curve, basea zachos et al data 2001 Category:Paleoclimatology
17:31, 3 August 2008Thumbnail for version as of 17:31, 3 August 2008840 × 350 (21 KB)Merikanto~commonswiki (talk | contribs){{Information |Description={{en|1=Raw delta-O18 curve between 3,60-2,50 millios ago, that is Reuverian or Piacenzian stage. Based on James Zachos et al data file 2001 compliationdata.txt on web site http://www.es.ucsc.edu/~jzachos/Publications.html. Y-axi

You cannot overwrite this file.

There are no pages that use this file.

File usage on other wikis

The following other wikis use this file:

Metadata

This file contains additional information such as Exif metadata which may have been added by the digital camera, scanner, or software program used to create or digitize it. If the file has been modified from its original state, some details such as the timestamp may not fully reflect those of the original file. The timestamp is only as accurate as the clock in the camera, and it may be completely wrong.

Retrieved from "https://commons.wikimedia.org/w/index.php?title=File:Reuverian_curve_zachos_2001_raw_data.PNG&oldid=893515272"