File:Ngrip20 yd pb 2.png

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

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

Captions

Captions

Add a one-line explanation of what this file represents

Summary

[edit]
Description
English: Younger dryas/Preboreal transiition in Greenland NGRIP oxygen isotope data. Oxygen isotope ratio changes with temperature.

Data for this file is from http://www.iceandclimate.nbi.ku.dk/data/NGRIP_d18O_and_dust_5cm.xls http://www.iceandclimate.nbi.ku.dk/data/ processed to csv

δ18O values and dust concentrations


The dataset provides NGRIP δ18O values, dust concentrations, and GICC05 ages in 5cm depth resolution for the period 0-60 ka (δ18O) and 10-60 ka (dust).

The dataset accompany the following papers:

NGRIP members, Nature, 431, 147-151, 2004. DOI: 10.1038/nature02805

Gkinis et al., Earth Planet. Sci. Lett., 405, 132-141, 2014. DOI: 10.1016/j.epsl.2014.08.022

Ruth et al., J. Geophys. Res., 108, 4098, 2003. DOI: 4010.1029/2002JD002376


    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

def running_mean(x, N): 

   cumsum = np.cumsum(np.insert(x, 0, 0)) 
   return (cumsum[N:] - cumsum[:-N]) / float(N)


datafilename="ngrip1.csv" captioni="YD/PB boundary in NGRIP ice core" savename="ngrip_dryas.svg"

figsizex=16 figsizey=8

  1. x = []
  2. y = []
  3. y2= []


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

lst1=['gicc05_age','delta_O18']

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

x0=dfin1['gicc05_age'] y0=dfin1['delta_O18']

  1. y20=dfin1['GISP_dO18']
  2. y30=dfin1['GISP2_dO18']


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


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


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

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


  1. y_running = running_mean(y, 31)


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


fig, ax1 = plt.subplots()


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

ax1.set_xlim(11000,12000) ax1.set_ylim(-33.0, -46.0)

  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="#b0b0ff", linewidth=1,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="#0000FF", linewidth=4, label="SavGol filter, 71 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.tick_params(axis='both', which='major', labelsize=size2)

ax1.xaxis.set_minor_locator(MultipleLocator(1000)) ax1.xaxis.set_minor_locator(MultipleLocator(200))

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 BP', color="darkgreen", fontsize=size2)


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()





Old stuff:


GICC05modelext time scale for the NGRIP ice core released 19 November 2010.

This image is gnuplotted from file

http://www.iceandclimate.nbi.ku.dk/data/2010-11-19_GICC05modelext_for_NGRIP.txt/

with this script


  1. !/gnuplot
  4. G N U P L O T
  5. Version 4.0 patchlevel 0
  6. last modified Thu Apr 15 14:44:22 CEST 2004
  7. System: MS-Windows 32 bit
  9. Copyright (C) 1986 - 1993, 1998, 2004
  10. Thomas Williams, Colin Kelley and many others
  12. This is gnuplot version 4.0. Please refer to the documentation
  13. for command syntax changes. The old syntax will be accepted
  14. throughout the 4.0 series, but all save files use the new syntax.
  16. Type `help` to access the on-line reference manual.
  17. The gnuplot FAQ is available from
  18. http://www.gnuplot.info/faq/
  20. Send comments and requests for help to
  21. <gnuplot-info@lists.sourceforge.net>
  22. Send bugs, suggestions and mods to
  23. <gnuplot-bugs@lists.sourceforge.net>
  25. set terminal png nocrop medium size 640,480
  26. set output 'ngrip20_yd_pb_2'

unset clip points set clip one unset clip two set bar 1.000000 set border 31 lt -1 lw 1.000 set xdata set ydata set zdata set x2data set y2data set timefmt x "%d/%m/%y,%H:%M" set timefmt y "%d/%m/%y,%H:%M" set timefmt z "%d/%m/%y,%H:%M" set timefmt x2 "%d/%m/%y,%H:%M" set timefmt y2 "%d/%m/%y,%H:%M" set timefmt cb "%d/%m/%y,%H:%M" set boxwidth set style fill empty border set dummy x,y set format x "% g" set format y "% g" set format x2 "% g" set format y2 "% g" set format z "% g" set format cb "% g" set angles radians unset grid set key title "" set key right top Right noreverse enhanced box linetype -2 linewidth 1.000 samplen 4 spacing 1 width 0 height 0 autotitles unset label unset arrow unset style line unset style arrow unset logscale set offsets 0, 0, 0, 0 set pointsize 1 set encoding default unset polar unset parametric unset decimalsign set view 60, 30, 1, 1 set samples 100, 100 set isosamples 10, 10 set surface unset contour set clabel '%8.3g' set mapping cartesian set datafile separator whitespace unset hidden3d set cntrparam order 4 set cntrparam linear set cntrparam levels auto 5 set cntrparam points 5 set size ratio 0 1,1 set origin 0,0 set style data points set style function lines set xzeroaxis lt -2 lw 1.000 set yzeroaxis lt -2 lw 1.000 set x2zeroaxis lt -2 lw 1.000 set y2zeroaxis lt -2 lw 1.000 set tics in set ticslevel 0.5 set ticscale 1 0.5 set mxtics default set mytics default set mztics default set mx2tics default set my2tics default set mcbtics default set xtics border mirror norotate 100 set ytics border mirror norotate autofreq set ztics border nomirror norotate autofreq set nox2tics set noy2tics set cbtics border mirror norotate autofreq set title "NGRIP 20yr YD/PB delta-O18" 0.000000,0.000000 font "" set timestamp "" bottom norotate 0.000000,0.000000 "" set rrange [ * : * ] noreverse nowriteback # (currently [0.000000:10.0000] ) set trange [ * : * ] noreverse nowriteback # (currently [-5.00000:5.00000] ) set urange [ * : * ] noreverse nowriteback # (currently [-5.00000:5.00000] ) set vrange [ * : * ] noreverse nowriteback # (currently [-5.00000:5.00000] ) set xlabel "Time yr BP" 0.000000,0.000000 font "" set x2label "" 0.000000,0.000000 font "" set xrange [ * : * ] noreverse nowriteback # (currently [-10.0000:10.0000] ) set x2range [ * : * ] noreverse nowriteback # (currently [-10.0000:10.0000] ) set ylabel "delta-O18" 0.000000,0.000000 font "" set y2label "" 0.000000,0.000000 font "" set yrange [ * : * ] noreverse nowriteback # (currently [-10.0000:10.0000] ) set y2range [ * : * ] noreverse nowriteback # (currently [-10.0000:10.0000] ) set zlabel "" 0.000000,0.000000 font "" set zrange [ * : * ] noreverse nowriteback # (currently [-10.0000:10.0000] ) set cblabel "" 0.000000,0.000000 font "" set cbrange [ * : * ] noreverse nowriteback # (currently [-10.0000:10.0000] ) set zero 1e-008 set lmargin -1 set bmargin -1 set rmargin -1 set tmargin -1 set locale "C" set pm3d scansautomatic flush begin noftriangles nohidden3d implicit corners2color mean unset pm3d set palette positive nops_allcF maxcolors 0 gamma 1.5 color model RGB set palette rgbformulae 7, 5, 15 set colorbox default set colorbox vertical origin 0.9,0.2 size 0.1,0.63 bdefault set loadpath set fontpath set fit noerrorvariables plot 'C:\gnuplot\bin\ng202.txt' using 1:3 with lines title "NGRIP delta-O8" lw 3

  1. EOF

Additional data form original datafile


% The splicing of the GICC05 and model time scale at 60.2 ka b2k (GICC05modelext) is described in the paper: % E.W. Wolff, J. Chappellaz, T. Blunier, S.O. Rasmussen, A. Svensson, % Millennial-scale variability during the last glacial: The ice core record % Quaternary Science Reviews 29, p. 2828-2838, 2010 % % Age NGRIP1/2 NGRIP Max. counting % (years depth d18O error

% b2k) (m) (permil) (years)
Date
Source Own work
Author Merikanto

Licensing

[edit]
I, the copyright holder of this work, hereby publish it under the following license:
w:en:Creative Commons
attribution share alike
This file is licensed under the Creative Commons Attribution-Share Alike 3.0 Unported license.
You are free:
  • to share – to copy, distribute and transmit the work
  • to remix – to adapt the work
Under the following conditions:
  • attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
  • share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.

File history

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

Date/TimeThumbnailDimensionsUserComment
current18:11, 12 September 2020Thumbnail for version as of 18:11, 12 September 20201,600 × 766 (207 KB)Merikanto (talk | contribs)New layout
15:37, 27 December 2013Thumbnail for version as of 15:37, 27 December 2013640 × 480 (6 KB)Merikanto~commonswiki (talk | contribs)User created page with UploadWizard

You cannot overwrite this file.

There are no pages that 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:Ngrip20_yd_pb_2.png&oldid=709485224"