# Graph of July insolation at 65 N, CO_2 and temperature anomalies # of the last 800 ka; gnuplot script by Jan Hollan # Data sources: # Insolation values -800-800ka1361.txt have been computed online # at https://www.imcce.fr/Equipes/ASD/insola/earth/online/index.php # based on # A&A 428, 261-285 (2004), DOI: 10.1051/0004-6361:20041335 # Laskar, J., Robutel, P., Joutel, F., Gastineau, # M., Correia, A.C.M., Levrard, B. : 2004, # A long term numerical solution for the insolation quantities # of the Earth. # t.txt is a symlink to edc3deuttemp2007.txt, as available at # http://www.ncdc.noaa.gov/paleo/metadata/noaa-icecore-6080.html # # SUGGESTED DATA CITATION: Jouzel, J., et al. 2007. # EPICA Dome C Ice Core 800KYr Deuterium Data and Temperature Estimates. # IGBP PAGES/World Data Center for Paleoclimatology # Data Contribution Series # 2007-091. # NOAA/NCDC Paleoclimatology Program, Boulder CO, USA. # # ORIGINAL REFERENCE: # Jouzel, J., V. Masson-Delmotte, O. Cattani, G. Dreyfus, S. Falourd, # G. Hoffmann, B. Minster, J. Nouet, J.M. Barnola, J. Chappellaz, H. Fischer, # J.C. Gallet, S. Johnsen, M. Leuenberger, L. Loulergue, D. Luethi, H. Oerter, # F. Parrenin, G. Raisbeck, D. Raynaud, A. Schilt, J. Schwander, E. Selmo, # R. Souchez, R. Spahni, B. Stauffer, J.P. Steffensen, B. Stenni, T.F. Stocker, # J.L. Tison, M. Werner, and E.W. Wolff. 2007. # Orbital and Millennial Antarctic Climate Variability over the Past 800,000 Years. # Science, Vol. 317, No. 5839, pp.793-797, 10 August 2007. # comp800.txt is the bottom part of the # ftp://ftp.ncdc.noaa.gov/pub/data/paleo/icecore/antarctica/epica_domec/edc-co2-2008.txt # ORIGINAL REFERENCE: # Lüthi, D., M. Le Floch, B. Bereiter, T. Blunier, J.-M. Barnola, # U. Siegenthaler, D. Raynaud, J. Jouzel, H. Fischer, K. Kawamura, # and T.F. Stocker. 2008. # High-resolution carbon dioxide concentration record 650,000-800,000 # years before present. # Nature, Vol. 453, pp. 379-382, 15 May 2008. doi:10.1038/nature06949 # The very script to plot them, with Czech titles: set encoding iso_8859_1 set term pdf color enhanced font "Arial-Medium,5" set out "800-800ka-en65.pdf" set lmargin 13 set rmargin 5.8 set bmargin 3.2 set tmargin 3 set multiplot set ytics nomirror set grid set title \ "{/*1.2 Irradiance of atmosphere, CO_2, temperature anomaly: past 800 thousand years} \n \ {/*1.2 ^ and irradiance in next 800 thousand years, mid July at 65° latitude}" offset 0,-0.2 set ylabel \ "Future irradiance / (2 W/m^2) + 90\n\ ^ Past irradiance / (2 W/m^2) + 50,\n\ ^ Concentration / ppm,\n\ ^ Antarctic temperature anomaly · (8/K) + 225,\n\ " offset 2.5,0 set xlabel "Age or Time in future / ka" offset 0,0.2 set ytics 20 plot [800:0] [140:420] \ "-800-800ka1361.txt" using 1:(90 + $2/2) \ with lines lt 1 lw 2 title 'Future irradiance', \ "-800-800ka1361.txt" using (-$1):(50 + $2/2) \ with lines lt 1 lw 3 title 'Past irradiance',\ "comp800.txt" using ($1/1000):2 with lines lt 3 \ title 'CO_2', \ "t.txt" using ($3/1000):(225+$5*8) with lines lt 2 \ title 'Temperature anomaly' set title '' set xlabel '' set ylabel '' unset xtics unset ytics set y2tics (" 0" 0+225, " 2" 2*16+225, "-2" (-2*16)+225, "-4" (-4*16)+225) set y2label "Global temp. anomaly / 1 K {/Symbol »} Antarctic one / 2 K " \ offset graph -0.02,-0.03 plot [800:0] [140:420] '-' w p ls 9 lt 3 notitle 0 400 e unset multiplot quit