• File
  • File history
  • File links

Co2-temperature-plot.svg‎ (SVG file, nominally 720 × 540 pixels, file size: 225 KB)


Summary

This figure shows historical CO2 (right axis) and reconstructed temperature (as a difference from the mean temperature for the last 100 years) records based on Antarctic ice cores, providing data for the last 800,000 years.

Licensing

This figure was produced by Leland McInnes using python and matplotlib and is licensed under the GFDL. All data is from publicly available sources.

Data Sources

  1. ( red ) EPICA Dome C temperature data : http://doi.pangaea.de/10.1594/PANGAEA.683655
  2. ( dark blue ) Vostok CO2 data : http://doi.pangaea.de/10.1594/PANGAEA.55501
  3. ( steel blue ) EPICA DomeC temperature data, 423-391 kybp : http://doi.pangaea.de/10.1594/PANGAEA.472482
  4. ( pale blue ) EPICA DomeC CO2 data, 650-413 kybp : http://doi.pangaea.de/10.1594/PANGAEA.472481
  5. ( cyan ) EPICA DomeC CO2 data, 800-650 kybp : http://doi.pangaea.de/10.1594/PANGAEA.710901
  6. Current CO2 level : http://en.wikipedia.org/wiki/Carbon_dioxide#Concentrations_of_CO2_in_atmosphere

Image Creation

The image was created with the following script, making use of the matplotlib library.

                        
                          # -*- coding: utf-8 -*-
                        
                        
                          import
                        
                        numpy
                        
                          as
                        
                        np
                        
                          import
                        
                        matplotlib.
                        
                          pylab
                        
                        
                          as
                        
                        plt
                        
                          from
                        
                        matplotlib.
                        
                          font_manager
                        
                        
                          import
                        
                        FontProperties
                        
                          def
                        
                        smooth
                        
                          (
                        
                        
                          signal
                        
                        , window_size
                        
                          )
                        
                        : extended_signal =
                        
                          signal
                        
                        
                        
                        +
                        
                          signal
                        
                        +
                        
                          signal
                        
                        
                        
                        s = np.
                        
                          array
                        
                        
                          (
                        
                        extended_signal
                        
                          )
                        
                        w = np.
                        
                          hamming
                        
                        
                          (
                        
                        window_size
                        
                          )
                        
                        y = np.
                        
                          convolve
                        
                        
                          (
                        
                        w/w.
                        
                          sum
                        
                        
                          (
                        
                        
                          )
                        
                        , s, mode=
                        
                          "same"
                        
                        
                          )
                        
                        
                          return
                        
                        y
                        
                        
                        temp_file =
                        
                          open
                        
                        
                          (
                        
                        
                          "EDC_dD_temp_estim.csv"
                        
                        
                          )
                        
                        data_rows =
                        
                        
                        temp_years =
                        
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          ]
                        
                        raw_temp_vals =
                        
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          ]
                        
                        temp_vals =
                        
                          list
                        
                        
                          (
                        
                        smooth
                        
                          (
                        
                        raw_temp_vals
                        
                        
                        , 50
                        
                          )
                        
                        
                          )
                        
                        +
                        
                          list
                        
                        
                          (
                        
                        smooth
                        
                          (
                        
                        raw_temp_vals
                        
                        
                        , 10
                        
                          )
                        
                        
                          )
                        
                        + raw_temp_vals
                        
                        
                        co2_file1 =
                        
                          open
                        
                        
                          (
                        
                        
                          "Vostok_CH4_CO2_age.csv"
                        
                        
                          )
                        
                        data_rows =
                        
                        
                        co2_years1 =
                        
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          if
                        
                        x
                        
                        
                        
                          !
                        
                        =
                        
                          "
                          
                            \n
                          
                          "
                        
                        
                          ]
                        
                        co2_vals1 =
                        
                        
                        .
                        
                          strip
                        
                        
                          (
                        
                        
                          )
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          if
                        
                        x
                        
                        
                        
                          !
                        
                        =
                        
                          "
                          
                            \n
                          
                          "
                        
                        
                          ]
                        
                        co2_file2 =
                        
                          open
                        
                        
                          (
                        
                        
                          "EDC_CO2_423-391_kyr_BP.csv"
                        
                        
                          )
                        
                        data_rows =
                        
                        
                        co2_years2 =
                        
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          ]
                        
                        co2_vals2 =
                        
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          ]
                        
                        co2_file3 =
                        
                          open
                        
                        
                          (
                        
                        
                          "EDC_CO2_650-413_kyr_BP.csv"
                        
                        
                          )
                        
                        data_rows =
                        
                        
                        co2_years3 =
                        
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          ]
                        
                        co2_vals3 =
                        
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          ]
                        
                        co2_file4 =
                        
                          open
                        
                        
                          (
                        
                        
                          "EDC99_CO2_bern.csv"
                        
                        
                          )
                        
                        data_rows =
                        
                        
                        co2_years4 =
                        
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          ]
                        
                        co2_vals4 =
                        
                        
                        
                          )
                        
                        
                          for
                        
                        x
                        
                          in
                        
                        data_rows
                        
                          ]
                        
                        plt.
                        
                          title
                        
                        
                          (
                        
                        
                          "Temperature and CO$_2$ Records"
                        
                        , size=
                        
                          28
                        
                        
                          )
                        
                        temp_ax = plt.
                        
                          axes
                        
                        
                          (
                        
                        
                          )
                        
                        temp_ax.
                        
                          yaxis
                        
                        .
                        
                          tick_left
                        
                        
                          (
                        
                        
                          )
                        
                        tline = plt.
                        
                          plot
                        
                        
                          (
                        
                        temp_years, temp_vals,
                        
                          "red"
                        
                        
                          )
                        
                        plt.
                        
                          ylabel
                        
                        
                          (
                        
                        u
                        
                          "Temperature anomaly (°C)"
                        
                        , size=
                        
                          16
                        
                        
                          )
                        
                        plt.
                        
                          xlabel
                        
                        
                          (
                        
                        
                          "Thousands of Years Ago"
                        
                        , size=16
                        
                          )
                        
                        plt.
                        
                          xlim
                        
                        
                          (
                        
                        800,0
                        
                          )
                        
                        plt.
                        
                          ylim
                        
                        
                          (
                        
                        -11,15
                        
                          )
                        
                        co2_ax = plt.
                        
                          axes
                        
                        
                          (
                        
                        
                        
                        , frameon=
                        
                          False
                        
                        
                          )
                        
                        co2_ax.
                        
                          yaxis
                        
                        .
                        
                          tick_right
                        
                        
                          (
                        
                        
                          )
                        
                        co2_ax.
                        
                          yaxis
                        
                        .
                        
                          set_label_position
                        
                        
                          (
                        
                        
                          "right"
                        
                        
                          )
                        
                        plt.
                        
                          xticks
                        
                        
                          (
                        
                        
                        
                        
                          )
                        
                        cline1 = plt.
                        
                          plot
                        
                        
                          (
                        
                        co2_years1, co2_vals1,
                        
                          "#0000AA"
                        
                        
                          )
                        
                        cline2 = plt.
                        
                          plot
                        
                        
                          (
                        
                        co2_years2, co2_vals2,
                        
                          "#2288AA"
                        
                        
                          )
                        
                        cline3 = plt.
                        
                          plot
                        
                        
                          (
                        
                        co2_years3, co2_vals3,
                        
                          "#44AAFF"
                        
                        
                          )
                        
                        cline4 = plt.
                        
                          plot
                        
                        
                          (
                        
                        co2_years4, co2_vals4,
                        
                          "#88DDFF"
                        
                        
                          )
                        
                        plt.
                        
                          ylabel
                        
                        
                          (
                        
                        
                          "Atmospheric CO$_2$ (ppm)"
                        
                        , size=
                        
                          16
                        
                        
                          )
                        
                        plt.
                        
                          xlim
                        
                        
                          (
                        
                        
                          800
                        
                        ,
                        
                          0
                        
                        
                          )
                        
                        plt.
                        
                          ylim
                        
                        
                          (
                        
                        
                          160
                        
                        ,
                        
                          420
                        
                        
                          )
                        
                        plt.
                        
                          annotate
                        
                        
                          (
                        
                        
                          "Current CO$_2$ level"
                        
                        ,
                        
                          (
                        
                        
                          2
                        
                        ,
                        
                          386.36
                        
                        
                          )
                        
                        ,
                        
                          (
                        
                        
                          300
                        
                        ,
                        
                          386.36
                        
                        
                          )
                        
                        , xycoords=
                        
                          "data"
                        
                        , arrowprops=
                        
                          dict
                        
                        
                          (
                        
                        arrowstyle=
                        
                          "->"
                        
                        
                          )
                        
                        
                          )
                        
                        plt.
                        
                          legend
                        
                        
                          (
                        
                        
                          (
                        
                        tline, cline1, cline2, cline3, cline4
                        
                          )
                        
                        ,
                        
                          (
                        
                        
                          "Temperature (EPICA Dome C)"
                        
                        ,
                        
                          "CO$_2$ (Vostok)"
                        
                        ,
                        
                          "CO$_2$ (EPICA Dome C)"
                        
                        ,
                        
                          "CO$_2$ (EPICA Dome C)"
                        
                        ,
                        
                          "CO$_2$ (EPICA Dome C)"
                        
                        
                          )
                        
                        ,
                        
                          "upper left"
                        
                        , prop=FontProperties
                        
                          (
                        
                        size=
                        
                          9
                        
                        
                          )
                        
                        , pad=
                        
                          0.1
                        
                        , handletextsep=
                        
                          0.005
                        
                        
                          )
                        
                        plt.
                        
                          savefig
                        
                        
                          (
                        
                        
                          "py-co2-temp.svg"
                        
                        
                          )
                        
                        plt.
                        
                          show
                        
                        
                          (
                        
                        
                          )
                        
                      

File history

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

(latest | earliest) View (newer 50) (older 50) (20 | 50 | 100 | 250 | 500)(latest | earliest) View (newer 50) (older 50) (20 | 50 | 100 | 250 | 500)
  • Edit this file using an external application

File links

The following pages on the English Wikipedia link to this file (pages on other projects are not listed):
  • User:Leland McInnes/Images
  • User:Isonomia/Global Fuel Depletion
  • The Great Global Warming Swindle

Hundred Chart and Skip Counting Patterns (Gr. 2) Printable (2nd Grade ...

Using a hundred chart and patterns, students explore multiplication.

...

Hundreds Chart

Free printable hundreds charts for early math practice. Children benefit from seeing the many patterns on these charts.

...

Printable 100s chart worksheets

blank 100's chart; fill in the hundred's chart; 2nd grade math ... Prev. Provide Your Own Instructions for the blank 100s ... Free Backstitch Motif Patterns - Free Printable Back Stitch ...

...

hundreds definition | Dictionary.com

Copy & paste this link to your blog or website to ... Printable blank hun... ... Math hundreds chart

...

Printable interactive 100 number chart worksheets- color number ...

Free interactive 100 chart worksheets. Children find number patterns and color code them with this interactive math game. Great way to reinforce skip counting - counting by 2's ...

...

Hundreds Chart Maths Worksheet, Australian Primary Mathematics

Hundreds Chart Maths Worksheet. Hundreds Chart - Print and keep a copy handy for grade year 4, Whole ... games, activities, educational, children, kids, online, free, fun, printable ...

...

Hundred Chart Printable (1st - 6th Grade) - TeacherVision.com

A printable hundred chart. ... Provided by Scott Foresman, an imprint of Pearson, the world's leading elementary educational publisher.

...

Printable Hundreds Chart: Customizable

A hundreds number chart that displays 1 to 100 in 10 rows and 10 columns and also gives the option for you to create your own version of the chart by picking a starting number and ...

...

Printable Hundreds Charts

Printable hundreds charts, ninety-nine charts, and place value charts for student reference.

...

Math Printables for Home

PRINTABLE Resources . Hundreds Chart - A 10 x 10 chart numbered to 100 to assist students with counting and number pattern recognition.

...