Atmospheric CO2 concentration
M. Delorme
31/05/2015
Graph of monthly mean carbon dioxide globally averaged over marine surface sites and graph of the carbon dioxide data, measured as the mole fraction in dry air, on Mauna Loa. This second graph constitutes the longest record of direct measurements of CO2 in the atmosphere.
library(ggplot2)
library(plyr)
library(Hmisc)
library(scales)
library(gridExtra)
# to format the date in english names
Sys.setlocale("LC_TIME", "en_US.UTF-8")1 Globally averaged marine surface monthly mean data (since 1980)
Data from Ed Dlugokencky and Pieter Tans, NOAA/ESRL
co2 <- read.table("ftp://aftp.cmdl.noaa.gov/products/trends/co2/co2_mm_gl.txt")
colnames(co2) <- c("year", "month", "decimal", "average", "trend")
co2$bp <- 1950 - co2$year # add before present years
ggplot(co2, aes(decimal, average)) +
geom_line(color = "pink") +
geom_point(color = "red", size = 1) +
stat_smooth() +
xlab("Year") +
ylab(expression(paste(CO[2], " fraction in dry air (", mu, "mol/mol)"))) +
ggtitle(substitute(paste("Globally averaged marine surface monthly mean ", CO[2], " concentration ", min, "-", max),
list(min = min(co2$year), max = max(co2$year)))) +
theme_bw()## geom_smooth: method="auto" and size of largest group is <1000, so using loess. Use 'method = x' to change the smoothing method.
co2## year month decimal average trend bp
## 1 1980 1 1980.042 338.45 337.82 -30
## 2 1980 2 1980.125 339.18 338.12 -30
## 3 1980 3 1980.208 339.49 338.14 -30
## 4 1980 4 1980.292 339.86 338.25 -30
## 5 1980 5 1980.375 340.30 338.78 -30
## 6 1980 6 1980.458 339.86 339.08 -30
## 7 1980 7 1980.542 338.32 339.17 -30
## 8 1980 8 1980.625 337.11 339.37 -30
## 9 1980 9 1980.708 336.95 339.33 -30
## 10 1980 10 1980.792 337.71 339.05 -30
## 11 1980 11 1980.875 338.83 339.17 -30
## 12 1980 12 1980.958 339.55 339.32 -30
## 13 1981 1 1981.042 340.10 339.47 -31
## 14 1981 2 1981.125 340.66 339.61 -31
## 15 1981 3 1981.208 341.29 339.95 -31
## 16 1981 4 1981.292 341.57 339.96 -31
## 17 1981 5 1981.375 341.28 339.76 -31
## 18 1981 6 1981.458 340.47 339.70 -31
## 19 1981 7 1981.542 339.10 339.95 -31
## 20 1981 8 1981.625 337.98 340.24 -31
## 21 1981 9 1981.708 337.86 340.24 -31
## 22 1981 10 1981.792 338.95 340.30 -31
## 23 1981 11 1981.875 340.06 340.39 -31
## 24 1981 12 1981.958 340.64 340.42 -31
## 25 1982 1 1982.042 341.27 340.66 -32
## 26 1982 2 1982.125 341.86 340.80 -32
## 27 1982 3 1982.208 342.13 340.75 -32
## 28 1982 4 1982.292 342.42 340.81 -32
## 29 1982 5 1982.375 342.27 340.75 -32
## 30 1982 6 1982.458 341.39 340.63 -32
## 31 1982 7 1982.542 339.66 340.51 -32
## 32 1982 8 1982.625 338.02 340.30 -32
## 33 1982 9 1982.708 338.08 340.44 -32
## 34 1982 10 1982.792 339.53 340.88 -32
## 35 1982 11 1982.875 340.85 341.19 -32
## 36 1982 12 1982.958 341.68 341.45 -32
## 37 1983 1 1983.042 342.28 341.62 -33
## 38 1983 2 1983.125 342.65 341.61 -33
## 39 1983 3 1983.208 342.94 341.60 -33
## 40 1983 4 1983.292 343.41 341.83 -33
## 41 1983 5 1983.375 343.72 342.20 -33
## 42 1983 6 1983.458 343.41 342.62 -33
## 43 1983 7 1983.542 342.05 342.82 -33
## 44 1983 8 1983.625 340.54 342.75 -33
## 45 1983 9 1983.708 340.48 342.89 -33
## 46 1983 10 1983.792 341.67 343.08 -33
## 47 1983 11 1983.875 342.71 343.06 -33
## 48 1983 12 1983.958 343.37 343.15 -33
## 49 1984 1 1984.042 344.22 343.56 -34
## 50 1984 2 1984.125 344.73 343.73 -34
## 51 1984 3 1984.208 344.85 343.50 -34
## 52 1984 4 1984.292 345.09 343.46 -34
## 53 1984 5 1984.375 345.24 343.70 -34
## 54 1984 6 1984.458 344.46 343.71 -34
## 55 1984 7 1984.542 343.10 343.91 -34
## 56 1984 8 1984.625 342.13 344.40 -34
## 57 1984 9 1984.708 342.03 344.47 -34
## 58 1984 10 1984.792 342.89 344.28 -34
## 59 1984 11 1984.875 344.07 344.38 -34
## 60 1984 12 1984.958 344.91 344.63 -34
## 61 1985 1 1985.042 345.26 344.58 -35
## 62 1985 2 1985.125 345.71 344.70 -35
## 63 1985 3 1985.208 346.57 345.25 -35
## 64 1985 4 1985.292 346.85 345.24 -35
## 65 1985 5 1985.375 346.67 345.10 -35
## 66 1985 6 1985.458 346.17 345.39 -35
## 67 1985 7 1985.542 344.87 345.69 -35
## 68 1985 8 1985.625 343.39 345.69 -35
## 69 1985 9 1985.708 343.33 345.78 -35
## 70 1985 10 1985.792 344.58 345.97 -35
## 71 1985 11 1985.875 345.71 346.03 -35
## 72 1985 12 1985.958 346.45 346.16 -35
## 73 1986 1 1986.042 347.02 346.33 -36
## 74 1986 2 1986.125 347.24 346.25 -36
## 75 1986 3 1986.208 347.62 346.29 -36
## 76 1986 4 1986.292 348.15 346.53 -36
## 77 1986 5 1986.375 348.25 346.68 -36
## 78 1986 6 1986.458 347.68 346.91 -36
## 79 1986 7 1986.542 346.22 347.05 -36
## 80 1986 8 1986.625 344.86 347.14 -36
## 81 1986 9 1986.708 344.76 347.16 -36
## 82 1986 10 1986.792 345.93 347.31 -36
## 83 1986 11 1986.875 347.25 347.59 -36
## 84 1986 12 1986.958 347.68 347.41 -36
## 85 1987 1 1987.042 347.92 347.14 -37
## 86 1987 2 1987.125 348.48 347.39 -37
## 87 1987 3 1987.208 349.21 347.80 -37
## 88 1987 4 1987.292 349.91 348.22 -37
## 89 1987 5 1987.375 350.19 348.62 -37
## 90 1987 6 1987.458 349.42 348.75 -37
## 91 1987 7 1987.542 347.88 348.85 -37
## 92 1987 8 1987.625 346.63 349.00 -37
## 93 1987 9 1987.708 346.67 349.12 -37
## 94 1987 10 1987.792 347.93 349.32 -37
## 95 1987 11 1987.875 349.20 349.53 -37
## 96 1987 12 1987.958 350.17 349.86 -37
## 97 1988 1 1988.042 350.90 350.11 -38
## 98 1988 2 1988.125 351.47 350.39 -38
## 99 1988 3 1988.208 351.86 350.40 -38
## 100 1988 4 1988.292 352.32 350.54 -38
## 101 1988 5 1988.375 352.46 350.83 -38
## 102 1988 6 1988.458 351.77 351.06 -38
## 103 1988 7 1988.542 350.31 351.29 -38
## 104 1988 8 1988.625 349.11 351.55 -38
## 105 1988 9 1988.708 349.27 351.79 -38
## 106 1988 10 1988.792 350.45 351.89 -38
## 107 1988 11 1988.875 351.64 352.00 -38
## 108 1988 12 1988.958 352.51 352.23 -38
## 109 1989 1 1989.042 353.05 352.23 -39
## 110 1989 2 1989.125 353.48 352.38 -39
## 111 1989 3 1989.208 354.03 352.60 -39
## 112 1989 4 1989.292 354.41 352.63 -39
## 113 1989 5 1989.375 354.25 352.57 -39
## 114 1989 6 1989.458 353.31 352.60 -39
## 115 1989 7 1989.542 351.55 352.57 -39
## 116 1989 8 1989.625 350.20 352.68 -39
## 117 1989 9 1989.708 350.53 353.07 -39
## 118 1989 10 1989.792 351.80 353.24 -39
## 119 1989 11 1989.875 352.98 353.32 -39
## 120 1989 12 1989.958 353.84 353.52 -39
## 121 1990 1 1990.042 354.40 353.65 -40
## 122 1990 2 1990.125 354.80 353.74 -40
## 123 1990 3 1990.208 355.11 353.69 -40
## 124 1990 4 1990.292 355.44 353.66 -40
## 125 1990 5 1990.375 355.35 353.67 -40
## 126 1990 6 1990.458 354.26 353.59 -40
## 127 1990 7 1990.542 352.65 353.72 -40
## 128 1990 8 1990.625 351.53 354.04 -40
## 129 1990 9 1990.708 351.72 354.26 -40
## 130 1990 10 1990.792 353.13 354.53 -40
## 131 1990 11 1990.875 354.36 354.65 -40
## 132 1990 12 1990.958 355.09 354.63 -40
## 133 1991 1 1991.042 355.65 354.89 -41
## 134 1991 2 1991.125 356.05 354.95 -41
## 135 1991 3 1991.208 356.61 355.18 -41
## 136 1991 4 1991.292 357.13 355.38 -41
## 137 1991 5 1991.375 357.02 355.39 -41
## 138 1991 6 1991.458 356.13 355.51 -41
## 139 1991 7 1991.542 354.45 355.55 -41
## 140 1991 8 1991.625 352.94 355.46 -41
## 141 1991 9 1991.708 352.86 355.39 -41
## 142 1991 10 1991.792 353.89 355.29 -41
## 143 1991 11 1991.875 355.04 355.29 -41
## 144 1991 12 1991.958 355.88 355.38 -41
## 145 1992 1 1992.042 356.47 355.67 -42
## 146 1992 2 1992.125 356.83 355.70 -42
## 147 1992 3 1992.208 357.20 355.73 -42
## 148 1992 4 1992.292 357.67 355.87 -42
## 149 1992 5 1992.375 357.73 356.09 -42
## 150 1992 6 1992.458 356.82 356.22 -42
## 151 1992 7 1992.542 355.12 356.26 -42
## 152 1992 8 1992.625 353.53 356.10 -42
## 153 1992 9 1992.708 353.45 356.02 -42
## 154 1992 10 1992.792 354.71 356.11 -42
## 155 1992 11 1992.875 355.89 356.14 -42
## 156 1992 12 1992.958 356.65 356.16 -42
## 157 1993 1 1993.042 357.12 356.28 -43
## 158 1993 2 1993.125 357.46 356.31 -43
## 159 1993 3 1993.208 357.87 356.43 -43
## 160 1993 4 1993.292 358.33 356.57 -43
## 161 1993 5 1993.375 358.29 356.65 -43
## 162 1993 6 1993.458 357.27 356.60 -43
## 163 1993 7 1993.542 355.63 356.64 -43
## 164 1993 8 1993.625 354.32 356.79 -43
## 165 1993 9 1993.708 354.36 356.95 -43
## 166 1993 10 1993.792 355.61 357.09 -43
## 167 1993 11 1993.875 356.85 357.19 -43
## 168 1993 12 1993.958 357.76 357.38 -43
## 169 1994 1 1994.042 358.36 357.50 -44
## 170 1994 2 1994.125 358.90 357.75 -44
## 171 1994 3 1994.208 359.24 357.79 -44
## 172 1994 4 1994.292 359.58 357.80 -44
## 173 1994 5 1994.375 359.60 357.92 -44
## 174 1994 6 1994.458 358.66 357.94 -44
## 175 1994 7 1994.542 357.16 358.13 -44
## 176 1994 8 1994.625 355.95 358.47 -44
## 177 1994 9 1994.708 355.90 358.56 -44
## 178 1994 10 1994.792 357.13 358.67 -44
## 179 1994 11 1994.875 358.60 358.95 -44
## 180 1994 12 1994.958 359.48 359.07 -44
## 181 1995 1 1995.042 360.03 359.18 -45
## 182 1995 2 1995.125 360.50 359.39 -45
## 183 1995 3 1995.208 360.97 359.60 -45
## 184 1995 4 1995.292 361.42 359.72 -45
## 185 1995 5 1995.375 361.32 359.70 -45
## 186 1995 6 1995.458 360.49 359.79 -45
## 187 1995 7 1995.542 358.86 359.85 -45
## 188 1995 8 1995.625 357.57 360.06 -45
## 189 1995 9 1995.708 357.90 360.51 -45
## 190 1995 10 1995.792 359.28 360.75 -45
## 191 1995 11 1995.875 360.63 360.90 -45
## 192 1995 12 1995.958 361.50 361.03 -45
## 193 1996 1 1996.042 362.02 361.09 -46
## 194 1996 2 1996.125 362.37 361.19 -46
## 195 1996 3 1996.208 362.68 361.25 -46
## 196 1996 4 1996.292 363.01 361.29 -46
## 197 1996 5 1996.375 363.15 361.55 -46
## 198 1996 6 1996.458 362.76 362.09 -46
## 199 1996 7 1996.542 361.48 362.53 -46
## 200 1996 8 1996.625 359.96 362.48 -46
## 201 1996 9 1996.708 359.57 362.19 -46
## 202 1996 10 1996.792 360.55 362.03 -46
## 203 1996 11 1996.875 361.66 361.95 -46
## 204 1996 12 1996.958 362.52 362.07 -46
## 205 1997 1 1997.042 363.15 362.20 -47
## 206 1997 2 1997.125 363.50 362.32 -47
## 207 1997 3 1997.208 363.87 362.46 -47
## 208 1997 4 1997.292 364.37 362.68 -47
## 209 1997 5 1997.375 364.41 362.85 -47
## 210 1997 6 1997.458 363.52 362.88 -47
## 211 1997 7 1997.542 361.83 362.87 -47
## 212 1997 8 1997.625 360.28 362.78 -47
## 213 1997 9 1997.708 360.27 362.88 -47
## 214 1997 10 1997.792 361.80 363.27 -47
## 215 1997 11 1997.875 363.49 363.75 -47
## 216 1997 12 1997.958 364.55 364.09 -47
## 217 1998 1 1998.042 365.06 364.12 -48
## 218 1998 2 1998.125 365.37 364.22 -48
## 219 1998 3 1998.208 365.84 364.45 -48
## 220 1998 4 1998.292 366.46 364.79 -48
## 221 1998 5 1998.375 366.70 365.16 -48
## 222 1998 6 1998.458 366.12 365.51 -48
## 223 1998 7 1998.542 364.66 365.71 -48
## 224 1998 8 1998.625 363.60 366.09 -48
## 225 1998 9 1998.708 363.87 366.46 -48
## 226 1998 10 1998.792 365.22 366.65 -48
## 227 1998 11 1998.875 366.49 366.72 -48
## 228 1998 12 1998.958 367.33 366.83 -48
## 229 1999 1 1999.042 367.98 367.05 -49
## 230 1999 2 1999.125 368.36 367.23 -49
## 231 1999 3 1999.208 368.74 367.36 -49
## 232 1999 4 1999.292 369.13 367.49 -49
## 233 1999 5 1999.375 369.06 367.55 -49
## 234 1999 6 1999.458 368.19 367.61 -49
## 235 1999 7 1999.542 366.47 367.51 -49
## 236 1999 8 1999.625 365.13 367.60 -49
## 237 1999 9 1999.708 365.32 367.88 -49
## 238 1999 10 1999.792 366.63 368.05 -49
## 239 1999 11 1999.875 367.92 368.13 -49
## 240 1999 12 1999.958 368.75 368.21 -49
## 241 2000 1 2000.042 369.26 368.35 -50
## 242 2000 2 2000.125 369.54 368.40 -50
## 243 2000 3 2000.208 369.85 368.46 -50
## 244 2000 4 2000.292 370.21 368.55 -50
## 245 2000 5 2000.375 370.11 368.58 -50
## 246 2000 6 2000.458 369.20 368.66 -50
## 247 2000 7 2000.542 367.79 368.93 -50
## 248 2000 8 2000.625 366.60 369.13 -50
## 249 2000 9 2000.708 366.56 369.11 -50
## 250 2000 10 2000.792 367.80 369.18 -50
## 251 2000 11 2000.875 369.18 369.33 -50
## 252 2000 12 2000.958 369.98 369.40 -50
## 253 2001 1 2001.042 370.52 369.64 -51
## 254 2001 2 2001.125 370.98 369.83 -51
## 255 2001 3 2001.208 371.37 369.96 -51
## 256 2001 4 2001.292 371.71 370.04 -51
## 257 2001 5 2001.375 371.63 370.11 -51
## 258 2001 6 2001.458 370.71 370.17 -51
## 259 2001 7 2001.542 369.30 370.43 -51
## 260 2001 8 2001.625 368.15 370.65 -51
## 261 2001 9 2001.708 368.19 370.72 -51
## 262 2001 10 2001.792 369.54 370.92 -51
## 263 2001 11 2001.875 370.90 371.09 -51
## 264 2001 12 2001.958 371.82 371.29 -51
## 265 2002 1 2002.042 372.34 371.38 -52
## 266 2002 2 2002.125 372.72 371.49 -52
## 267 2002 3 2002.208 373.22 371.71 -52
## 268 2002 4 2002.292 373.56 371.81 -52
## 269 2002 5 2002.375 373.55 371.97 -52
## 270 2002 6 2002.458 372.69 372.14 -52
## 271 2002 7 2002.542 371.26 372.41 -52
## 272 2002 8 2002.625 370.19 372.73 -52
## 273 2002 9 2002.708 370.49 373.11 -52
## 274 2002 10 2002.792 371.75 373.23 -52
## 275 2002 11 2002.875 373.12 373.39 -52
## 276 2002 12 2002.958 374.13 373.66 -52
## 277 2003 1 2003.042 374.78 373.78 -53
## 278 2003 2 2003.125 375.31 374.02 -53
## 279 2003 3 2003.208 375.75 374.18 -53
## 280 2003 4 2003.292 376.22 374.43 -53
## 281 2003 5 2003.375 376.35 374.75 -53
## 282 2003 6 2003.458 375.51 374.95 -53
## 283 2003 7 2003.542 373.96 375.10 -53
## 284 2003 8 2003.625 372.73 375.31 -53
## 285 2003 9 2003.708 372.93 375.60 -53
## 286 2003 10 2003.792 374.22 375.76 -53
## 287 2003 11 2003.875 375.51 375.83 -53
## 288 2003 12 2003.958 376.34 375.90 -53
## 289 2004 1 2004.042 376.99 376.02 -54
## 290 2004 2 2004.125 377.52 376.22 -54
## 291 2004 3 2004.208 377.94 376.35 -54
## 292 2004 4 2004.292 378.30 376.51 -54
## 293 2004 5 2004.375 378.23 376.63 -54
## 294 2004 6 2004.458 377.35 376.80 -54
## 295 2004 7 2004.542 375.81 376.99 -54
## 296 2004 8 2004.625 374.36 377.00 -54
## 297 2004 9 2004.708 374.25 376.92 -54
## 298 2004 10 2004.792 375.56 377.07 -54
## 299 2004 11 2004.875 377.02 377.31 -54
## 300 2004 12 2004.958 378.00 377.51 -54
## 301 2005 1 2005.042 378.59 377.63 -55
## 302 2005 2 2005.125 379.03 377.73 -55
## 303 2005 3 2005.208 379.66 378.06 -55
## 304 2005 4 2005.292 380.17 378.35 -55
## 305 2005 5 2005.375 380.26 378.63 -55
## 306 2005 6 2005.458 379.44 378.84 -55
## 307 2005 7 2005.542 377.75 378.91 -55
## 308 2005 8 2005.625 376.52 379.15 -55
## 309 2005 9 2005.708 376.59 379.28 -55
## 310 2005 10 2005.792 377.91 379.48 -55
## 311 2005 11 2005.875 379.40 379.73 -55
## 312 2005 12 2005.958 380.38 379.90 -55
## 313 2006 1 2006.042 381.10 380.11 -56
## 314 2006 2 2006.125 381.75 380.40 -56
## 315 2006 3 2006.208 382.16 380.54 -56
## 316 2006 4 2006.292 382.48 380.65 -56
## 317 2006 5 2006.375 382.42 380.76 -56
## 318 2006 6 2006.458 381.56 380.94 -56
## 319 2006 7 2006.542 379.88 381.07 -56
## 320 2006 8 2006.625 378.31 381.01 -56
## 321 2006 9 2006.708 378.42 381.18 -56
## 322 2006 10 2006.792 379.82 381.40 -56
## 323 2006 11 2006.875 381.18 381.51 -56
## 324 2006 12 2006.958 382.15 381.67 -56
## 325 2007 1 2007.042 382.82 381.82 -57
## 326 2007 2 2007.125 383.35 381.97 -57
## 327 2007 3 2007.208 383.80 382.17 -57
## 328 2007 4 2007.292 384.09 382.27 -57
## 329 2007 5 2007.375 383.96 382.32 -57
## 330 2007 6 2007.458 383.09 382.51 -57
## 331 2007 7 2007.542 381.40 382.64 -57
## 332 2007 8 2007.625 380.13 382.86 -57
## 333 2007 9 2007.708 380.50 383.27 -57
## 334 2007 10 2007.792 381.89 383.44 -57
## 335 2007 11 2007.875 383.21 383.49 -57
## 336 2007 12 2007.958 384.21 383.70 -57
## 337 2008 1 2008.042 384.98 383.97 -58
## 338 2008 2 2008.125 385.49 384.11 -58
## 339 2008 3 2008.208 385.90 384.29 -58
## 340 2008 4 2008.292 386.31 384.50 -58
## 341 2008 5 2008.375 386.29 384.67 -58
## 342 2008 6 2008.458 385.35 384.79 -58
## 343 2008 7 2008.542 383.82 385.13 -58
## 344 2008 8 2008.625 382.50 385.28 -58
## 345 2008 9 2008.708 382.29 385.03 -58
## 346 2008 10 2008.792 383.43 384.91 -58
## 347 2008 11 2008.875 384.92 385.14 -58
## 348 2008 12 2008.958 386.01 385.46 -58
## 349 2009 1 2009.042 386.77 385.74 -59
## 350 2009 2 2009.125 387.21 385.79 -59
## 351 2009 3 2009.208 387.43 385.78 -59
## 352 2009 4 2009.292 387.73 385.92 -59
## 353 2009 5 2009.375 387.71 386.12 -59
## 354 2009 6 2009.458 386.74 386.24 -59
## 355 2009 7 2009.542 384.87 386.26 -59
## 356 2009 8 2009.625 383.50 386.34 -59
## 357 2009 9 2009.708 383.71 386.43 -59
## 358 2009 10 2009.792 385.24 386.67 -59
## 359 2009 11 2009.875 386.77 386.95 -59
## 360 2009 12 2009.958 387.70 387.14 -59
## 361 2010 1 2010.042 388.45 387.44 -60
## 362 2010 2 2010.125 389.14 387.75 -60
## 363 2010 3 2010.208 389.47 387.85 -60
## 364 2010 4 2010.292 389.76 387.97 -60
## 365 2010 5 2010.375 389.72 388.13 -60
## 366 2010 6 2010.458 388.80 388.30 -60
## 367 2010 7 2010.542 387.21 388.58 -60
## 368 2010 8 2010.625 386.09 388.87 -60
## 369 2010 9 2010.708 386.53 389.22 -60
## 370 2010 10 2010.792 388.09 389.52 -60
## 371 2010 11 2010.875 389.50 389.68 -60
## 372 2010 12 2010.958 390.22 389.68 -60
## 373 2011 1 2011.042 390.74 389.68 -61
## 374 2011 2 2011.125 391.19 389.77 -61
## 375 2011 3 2011.208 391.50 389.89 -61
## 376 2011 4 2011.292 391.88 390.09 -61
## 377 2011 5 2011.375 391.89 390.28 -61
## 378 2011 6 2011.458 390.97 390.47 -61
## 379 2011 7 2011.542 389.05 390.41 -61
## 380 2011 8 2011.625 387.75 390.51 -61
## 381 2011 9 2011.708 388.16 390.87 -61
## 382 2011 10 2011.792 389.73 391.20 -61
## 383 2011 11 2011.875 391.07 391.28 -61
## 384 2011 12 2011.958 391.89 391.38 -61
## 385 2012 1 2012.042 392.43 391.40 -62
## 386 2012 2 2012.125 393.02 391.62 -62
## 387 2012 3 2012.208 393.54 391.94 -62
## 388 2012 4 2012.292 393.81 392.01 -62
## 389 2012 5 2012.375 393.69 392.10 -62
## 390 2012 6 2012.458 392.66 392.15 -62
## 391 2012 7 2012.542 390.88 392.25 -62
## 392 2012 8 2012.625 389.77 392.53 -62
## 393 2012 9 2012.708 390.37 393.07 -62
## 394 2012 10 2012.792 391.96 393.41 -62
## 395 2012 11 2012.875 393.37 393.57 -62
## 396 2012 12 2012.958 394.19 393.64 -62
## 397 2013 1 2013.042 394.90 393.87 -63
## 398 2013 2 2013.125 395.51 394.11 -63
## 399 2013 3 2013.208 396.06 394.47 -63
## 400 2013 4 2013.292 396.53 394.73 -63
## 401 2013 5 2013.375 396.58 394.99 -63
## 402 2013 6 2013.458 395.79 395.28 -63
## 403 2013 7 2013.542 394.28 395.65 -63
## 404 2013 8 2013.625 393.08 395.83 -63
## 405 2013 9 2013.708 393.07 395.77 -63
## 406 2013 10 2013.792 394.34 395.78 -63
## 407 2013 11 2013.875 395.77 395.97 -63
## 408 2013 12 2013.958 396.72 396.18 -63
## 409 2014 1 2014.042 397.44 396.41 -64
## 410 2014 2 2014.125 397.83 396.43 -64
## 411 2014 3 2014.208 398.05 396.46 -64
## 412 2014 4 2014.292 398.41 396.61 -64
## 413 2014 5 2014.375 398.49 396.91 -64
## 414 2014 6 2014.458 397.52 397.00 -64
## 415 2014 7 2014.542 395.90 397.28 -64
## 416 2014 8 2014.625 394.80 397.55 -64
## 417 2014 9 2014.708 394.91 397.61 -64
## 418 2014 10 2014.792 396.19 397.63 -64
## 419 2014 11 2014.875 397.70 397.90 -64
## 420 2014 12 2014.958 398.62 398.07 -64
## 421 2015 1 2015.042 399.24 398.21 -65
## 422 2015 2 2015.125 399.84 398.43 -65
## 423 2015 3 2015.208 400.34 398.75 -65
## 424 2015 4 2015.292 400.79 398.99 -65
## 425 2015 5 2015.375 400.65 399.07 -65
## 426 2015 6 2015.458 400.11 399.59 -652 Mauna Loa CO2 monthly mean data (since 1958)
Data from Dr. Pieter Tans, NOAA/ESRL and Dr. Ralph Keeling, Scripps Institution of Oceanography.
Code used in Wikimedia commons.
# get data
co2ml <- read.table("ftp://aftp.cmdl.noaa.gov/products/trends/co2/co2_mm_mlo.txt", na.strings = c(-99.99, -1))
colnames(co2ml) <- c("year", "month", "decimal", "average", "interpolated", "trend", "ndays")
# main plot
p1 <- ggplot(co2ml, aes(decimal, interpolated)) +
geom_line(color = "pink") +
geom_point(color = "red", size=1) +
stat_smooth() +
xlab("Year") +
ylab(expression(paste(CO[2], " fraction in dry air (", mu, "mol/mol)"))) +
ggtitle(substitute(paste("Mauna Loa monthly mean ", CO[2], " concentration ", min, "-", max), list(min = min(co2ml$year), max = max(co2ml$year)))) +
theme_bw()
# calculate the yearly means
m <- aggregate(co2ml$interpolated, list(co2ml$year), mean)
colnames(m) <- c("year", "year_mean")
# join the yearly means to the dataset to calculate the monthly deviations to the yearly mean
mco2ml <- join(co2ml, m, by = "year")
mco2ml$delta <- mco2ml$interpolated - mco2ml$year_mean
# virtual date for a nice x scale
mco2ml$vdate <- as.Date(paste0("2015-", sprintf("%02d", mco2ml$month), "-01"))
# inset plot
p2 <- ggplot(mco2ml, aes(vdate, delta)) +
geom_hline(yintercept = 0) +
stat_summary(fun.data = "mean_cl_boot", colour = "red") +
stat_smooth(span = 0.4) +
xlab("Month") +
scale_x_date(labels = date_format("%b")) +
ylab(expression(atop(paste(CO[2], " fraction in dry air (", mu, "mol/mol)"), "Departure from yearly average"))) +
ggtitle(expression(paste("Monthly ", CO[2], " variation (yearly average deviation)"))) +
theme_bw()
# merge the plots
p2_grob <- ggplotGrob(p2)## geom_smooth: method="auto" and size of largest group is <1000, so using loess. Use 'method = x' to change the smoothing method.pcombined <- p1 + annotation_custom(grob = p2_grob, xmin = 1955, xmax = 1991, ymin =357, ymax = 409)
pcombined## geom_smooth: method="auto" and size of largest group is <1000, so using loess. Use 'method = x' to change the smoothing method.
#ggsave(pcombined, file = "co2_mauna_loa.svg", width = 20, height = 20, units = "cm")
co2ml## year month decimal average interpolated trend ndays
## 1 1958 3 1958.208 315.71 315.71 314.62 NA
## 2 1958 4 1958.292 317.45 317.45 315.29 NA
## 3 1958 5 1958.375 317.50 317.50 314.71 NA
## 4 1958 6 1958.458 NA 317.10 314.85 NA
## 5 1958 7 1958.542 315.86 315.86 314.98 NA
## 6 1958 8 1958.625 314.93 314.93 315.94 NA
## 7 1958 9 1958.708 313.20 313.20 315.91 NA
## 8 1958 10 1958.792 NA 312.66 315.61 NA
## 9 1958 11 1958.875 313.33 313.33 315.31 NA
## 10 1958 12 1958.958 314.67 314.67 315.61 NA
## 11 1959 1 1959.042 315.62 315.62 315.70 NA
## 12 1959 2 1959.125 316.38 316.38 315.88 NA
## 13 1959 3 1959.208 316.71 316.71 315.62 NA
## 14 1959 4 1959.292 317.72 317.72 315.56 NA
## 15 1959 5 1959.375 318.29 318.29 315.50 NA
## 16 1959 6 1959.458 318.15 318.15 315.92 NA
## 17 1959 7 1959.542 316.54 316.54 315.66 NA
## 18 1959 8 1959.625 314.80 314.80 315.81 NA
## 19 1959 9 1959.708 313.84 313.84 316.55 NA
## 20 1959 10 1959.792 313.26 313.26 316.19 NA
## 21 1959 11 1959.875 314.80 314.80 316.78 NA
## 22 1959 12 1959.958 315.58 315.58 316.52 NA
## 23 1960 1 1960.042 316.43 316.43 316.51 NA
## 24 1960 2 1960.125 316.97 316.97 316.47 NA
## 25 1960 3 1960.208 317.58 317.58 316.49 NA
## 26 1960 4 1960.292 319.02 319.02 316.86 NA
## 27 1960 5 1960.375 320.03 320.03 317.24 NA
## 28 1960 6 1960.458 319.59 319.59 317.36 NA
## 29 1960 7 1960.542 318.18 318.18 317.30 NA
## 30 1960 8 1960.625 315.91 315.91 316.92 NA
## 31 1960 9 1960.708 314.16 314.16 316.87 NA
## 32 1960 10 1960.792 313.83 313.83 316.76 NA
## 33 1960 11 1960.875 315.00 315.00 316.98 NA
## 34 1960 12 1960.958 316.19 316.19 317.13 NA
## 35 1961 1 1961.042 316.93 316.93 317.03 NA
## 36 1961 2 1961.125 317.70 317.70 317.28 NA
## 37 1961 3 1961.208 318.54 318.54 317.47 NA
## 38 1961 4 1961.292 319.48 319.48 317.27 NA
## 39 1961 5 1961.375 320.58 320.58 317.70 NA
## 40 1961 6 1961.458 319.77 319.77 317.48 NA
## 41 1961 7 1961.542 318.57 318.57 317.70 NA
## 42 1961 8 1961.625 316.79 316.79 317.80 NA
## 43 1961 9 1961.708 314.80 314.80 317.49 NA
## 44 1961 10 1961.792 315.38 315.38 318.35 NA
## 45 1961 11 1961.875 316.10 316.10 318.13 NA
## 46 1961 12 1961.958 317.01 317.01 317.94 NA
## 47 1962 1 1962.042 317.94 317.94 318.06 NA
## 48 1962 2 1962.125 318.56 318.56 318.11 NA
## 49 1962 3 1962.208 319.68 319.68 318.57 NA
## 50 1962 4 1962.292 320.63 320.63 318.45 NA
## 51 1962 5 1962.375 321.01 321.01 318.20 NA
## 52 1962 6 1962.458 320.55 320.55 318.27 NA
## 53 1962 7 1962.542 319.58 319.58 318.67 NA
## 54 1962 8 1962.625 317.40 317.40 318.48 NA
## 55 1962 9 1962.708 316.26 316.26 319.03 NA
## 56 1962 10 1962.792 315.42 315.42 318.33 NA
## 57 1962 11 1962.875 316.69 316.69 318.62 NA
## 58 1962 12 1962.958 317.69 317.69 318.61 NA
## 59 1963 1 1963.042 318.74 318.74 318.91 NA
## 60 1963 2 1963.125 319.08 319.08 318.68 NA
## 61 1963 3 1963.208 319.86 319.86 318.69 NA
## 62 1963 4 1963.292 321.39 321.39 319.09 NA
## 63 1963 5 1963.375 322.25 322.25 319.39 NA
## 64 1963 6 1963.458 321.47 321.47 319.16 NA
## 65 1963 7 1963.542 319.74 319.74 318.77 NA
## 66 1963 8 1963.625 317.77 317.77 318.83 NA
## 67 1963 9 1963.708 316.21 316.21 319.06 NA
## 68 1963 10 1963.792 315.99 315.99 319.00 NA
## 69 1963 11 1963.875 317.12 317.12 319.10 NA
## 70 1963 12 1963.958 318.31 318.31 319.25 NA
## 71 1964 1 1964.042 319.57 319.57 319.67 NA
## 72 1964 2 1964.125 NA 320.07 319.61 NA
## 73 1964 3 1964.208 NA 320.73 319.55 NA
## 74 1964 4 1964.292 NA 321.77 319.48 NA
## 75 1964 5 1964.375 322.25 322.25 319.42 NA
## 76 1964 6 1964.458 321.89 321.89 319.69 NA
## 77 1964 7 1964.542 320.44 320.44 319.58 NA
## 78 1964 8 1964.625 318.70 318.70 319.81 NA
## 79 1964 9 1964.708 316.70 316.70 319.56 NA
## 80 1964 10 1964.792 316.79 316.79 319.78 NA
## 81 1964 11 1964.875 317.79 317.79 319.72 NA
## 82 1964 12 1964.958 318.71 318.71 319.59 NA
## 83 1965 1 1965.042 319.44 319.44 319.48 NA
## 84 1965 2 1965.125 320.44 320.44 319.97 NA
## 85 1965 3 1965.208 320.89 320.89 319.65 NA
## 86 1965 4 1965.292 322.13 322.13 319.80 NA
## 87 1965 5 1965.375 322.16 322.16 319.36 NA
## 88 1965 6 1965.458 321.87 321.87 319.65 NA
## 89 1965 7 1965.542 321.39 321.39 320.51 NA
## 90 1965 8 1965.625 318.81 318.81 319.93 NA
## 91 1965 9 1965.708 317.81 317.81 320.68 NA
## 92 1965 10 1965.792 317.30 317.30 320.36 NA
## 93 1965 11 1965.875 318.87 318.87 320.87 NA
## 94 1965 12 1965.958 319.42 319.42 320.26 NA
## 95 1966 1 1966.042 320.62 320.62 320.63 NA
## 96 1966 2 1966.125 321.59 321.59 321.10 NA
## 97 1966 3 1966.208 322.39 322.39 321.16 NA
## 98 1966 4 1966.292 323.87 323.87 321.51 NA
## 99 1966 5 1966.375 324.01 324.01 321.18 NA
## 100 1966 6 1966.458 323.75 323.75 321.52 NA
## 101 1966 7 1966.542 322.39 322.39 321.49 NA
## 102 1966 8 1966.625 320.37 320.37 321.50 NA
## 103 1966 9 1966.708 318.64 318.64 321.54 NA
## 104 1966 10 1966.792 318.10 318.10 321.18 NA
## 105 1966 11 1966.875 319.79 319.79 321.84 NA
## 106 1966 12 1966.958 321.08 321.08 321.95 NA
## 107 1967 1 1967.042 322.07 322.07 322.07 NA
## 108 1967 2 1967.125 322.50 322.50 321.94 NA
## 109 1967 3 1967.208 323.04 323.04 321.72 NA
## 110 1967 4 1967.292 324.42 324.42 322.05 NA
## 111 1967 5 1967.375 325.00 325.00 322.27 NA
## 112 1967 6 1967.458 324.09 324.09 321.94 NA
## 113 1967 7 1967.542 322.55 322.55 321.66 NA
## 114 1967 8 1967.625 320.92 320.92 322.04 NA
## 115 1967 9 1967.708 319.31 319.31 322.19 NA
## 116 1967 10 1967.792 319.31 319.31 322.36 NA
## 117 1967 11 1967.875 320.72 320.72 322.78 NA
## 118 1967 12 1967.958 321.96 321.96 322.86 NA
## 119 1968 1 1968.042 322.57 322.57 322.55 NA
## 120 1968 2 1968.125 323.15 323.15 322.56 NA
## 121 1968 3 1968.208 323.89 323.89 322.59 NA
## 122 1968 4 1968.292 325.02 325.02 322.73 NA
## 123 1968 5 1968.375 325.57 325.57 322.87 NA
## 124 1968 6 1968.458 325.36 325.36 323.20 NA
## 125 1968 7 1968.542 324.14 324.14 323.25 NA
## 126 1968 8 1968.625 322.03 322.03 323.15 NA
## 127 1968 9 1968.708 320.41 320.41 323.31 NA
## 128 1968 10 1968.792 320.25 320.25 323.32 NA
## 129 1968 11 1968.875 321.31 321.31 323.32 NA
## 130 1968 12 1968.958 322.84 322.84 323.69 NA
## 131 1969 1 1969.042 324.00 324.00 323.98 NA
## 132 1969 2 1969.125 324.42 324.42 323.89 NA
## 133 1969 3 1969.208 325.64 325.64 324.41 NA
## 134 1969 4 1969.292 326.66 326.66 324.35 NA
## 135 1969 5 1969.375 327.34 327.34 324.57 NA
## 136 1969 6 1969.458 326.76 326.76 324.63 NA
## 137 1969 7 1969.542 325.88 325.88 325.08 NA
## 138 1969 8 1969.625 323.67 323.67 324.80 NA
## 139 1969 9 1969.708 322.38 322.38 325.28 NA
## 140 1969 10 1969.792 321.78 321.78 324.84 NA
## 141 1969 11 1969.875 322.85 322.85 324.78 NA
## 142 1969 12 1969.958 324.11 324.11 324.88 NA
## 143 1970 1 1970.042 325.03 325.03 325.04 NA
## 144 1970 2 1970.125 325.99 325.99 325.42 NA
## 145 1970 3 1970.208 326.87 326.87 325.69 NA
## 146 1970 4 1970.292 328.13 328.13 325.86 NA
## 147 1970 5 1970.375 328.07 328.07 325.27 NA
## 148 1970 6 1970.458 327.66 327.66 325.52 NA
## 149 1970 7 1970.542 326.35 326.35 325.51 NA
## 150 1970 8 1970.625 324.69 324.69 325.76 NA
## 151 1970 9 1970.708 323.10 323.10 325.93 NA
## 152 1970 10 1970.792 323.16 323.16 326.15 NA
## 153 1970 11 1970.875 323.98 323.98 325.96 NA
## 154 1970 12 1970.958 325.13 325.13 326.06 NA
## 155 1971 1 1971.042 326.17 326.17 326.26 NA
## 156 1971 2 1971.125 326.68 326.68 326.09 NA
## 157 1971 3 1971.208 327.18 327.18 325.94 NA
## 158 1971 4 1971.292 327.78 327.78 325.47 NA
## 159 1971 5 1971.375 328.92 328.92 326.11 NA
## 160 1971 6 1971.458 328.57 328.57 326.41 NA
## 161 1971 7 1971.542 327.34 327.34 326.46 NA
## 162 1971 8 1971.625 325.46 325.46 326.49 NA
## 163 1971 9 1971.708 323.36 323.36 326.18 NA
## 164 1971 10 1971.792 323.57 323.57 326.58 NA
## 165 1971 11 1971.875 324.80 324.80 326.82 NA
## 166 1971 12 1971.958 326.01 326.01 327.02 NA
## 167 1972 1 1972.042 326.77 326.77 326.86 NA
## 168 1972 2 1972.125 327.63 327.63 327.05 NA
## 169 1972 3 1972.208 327.75 327.75 326.52 NA
## 170 1972 4 1972.292 329.72 329.72 327.41 NA
## 171 1972 5 1972.375 330.07 330.07 327.21 NA
## 172 1972 6 1972.458 329.09 329.09 326.94 NA
## 173 1972 7 1972.542 328.05 328.05 327.24 NA
## 174 1972 8 1972.625 326.32 326.32 327.37 NA
## 175 1972 9 1972.708 324.93 324.93 327.73 NA
## 176 1972 10 1972.792 325.06 325.06 328.05 NA
## 177 1972 11 1972.875 326.50 326.50 328.50 NA
## 178 1972 12 1972.958 327.55 327.55 328.56 NA
## 179 1973 1 1973.042 328.54 328.54 328.58 NA
## 180 1973 2 1973.125 329.56 329.56 328.87 NA
## 181 1973 3 1973.208 330.30 330.30 328.99 NA
## 182 1973 4 1973.292 331.50 331.50 329.12 NA
## 183 1973 5 1973.375 332.48 332.48 329.60 NA
## 184 1973 6 1973.458 332.07 332.07 329.96 NA
## 185 1973 7 1973.542 330.87 330.87 330.09 NA
## 186 1973 8 1973.625 329.31 329.31 330.43 NA
## 187 1973 9 1973.708 327.51 327.51 330.43 NA
## 188 1973 10 1973.792 327.18 327.18 330.23 NA
## 189 1973 11 1973.875 328.16 328.16 330.17 NA
## 190 1973 12 1973.958 328.64 328.64 329.66 NA
## 191 1974 1 1974.042 329.35 329.35 329.44 NA
## 192 1974 2 1974.125 330.71 330.71 330.11 NA
## 193 1974 3 1974.208 331.48 331.48 330.19 NA
## 194 1974 4 1974.292 332.65 332.65 330.25 NA
## 195 1974 5 1974.375 333.19 333.19 330.25 15
## 196 1974 6 1974.458 332.12 332.12 329.93 27
## 197 1974 7 1974.542 330.99 330.99 330.20 31
## 198 1974 8 1974.625 329.17 329.17 330.31 31
## 199 1974 9 1974.708 327.41 327.41 330.35 29
## 200 1974 10 1974.792 327.21 327.21 330.29 29
## 201 1974 11 1974.875 328.34 328.34 330.33 29
## 202 1974 12 1974.958 329.50 329.50 330.47 31
## 203 1975 1 1975.042 330.68 330.68 330.81 31
## 204 1975 2 1975.125 331.41 331.41 330.86 28
## 205 1975 3 1975.208 331.85 331.85 330.49 26
## 206 1975 4 1975.292 333.29 333.29 330.77 30
## 207 1975 5 1975.375 333.91 333.91 330.94 30
## 208 1975 6 1975.458 333.40 333.40 331.18 29
## 209 1975 7 1975.542 331.74 331.74 330.92 27
## 210 1975 8 1975.625 329.88 329.88 331.06 31
## 211 1975 9 1975.708 328.57 328.57 331.54 30
## 212 1975 10 1975.792 328.36 328.36 331.46 17
## 213 1975 11 1975.875 329.33 329.33 331.39 25
## 214 1975 12 1975.958 NA 330.59 331.60 1
## 215 1976 1 1976.042 331.66 331.66 331.82 25
## 216 1976 2 1976.125 332.75 332.75 332.19 29
## 217 1976 3 1976.208 333.46 333.46 331.99 29
## 218 1976 4 1976.292 334.78 334.78 332.24 25
## 219 1976 5 1976.375 334.78 334.78 331.86 31
## 220 1976 6 1976.458 334.06 334.06 331.74 25
## 221 1976 7 1976.542 332.95 332.95 332.06 21
## 222 1976 8 1976.625 330.64 330.64 331.87 31
## 223 1976 9 1976.708 328.96 328.96 331.94 28
## 224 1976 10 1976.792 328.77 328.77 331.92 30
## 225 1976 11 1976.875 330.18 330.18 332.33 26
## 226 1976 12 1976.958 331.65 331.65 332.69 24
## 227 1977 1 1977.042 332.69 332.69 332.75 29
## 228 1977 2 1977.125 333.23 333.23 332.60 26
## 229 1977 3 1977.208 334.97 334.97 333.36 29
## 230 1977 4 1977.292 336.03 336.03 333.39 29
## 231 1977 5 1977.375 336.82 336.82 333.87 30
## 232 1977 6 1977.458 336.10 336.10 333.76 25
## 233 1977 7 1977.542 334.79 334.79 333.96 28
## 234 1977 8 1977.625 332.53 332.53 333.92 30
## 235 1977 9 1977.708 331.19 331.19 334.27 27
## 236 1977 10 1977.792 331.21 331.21 334.43 30
## 237 1977 11 1977.875 332.35 332.35 334.59 28
## 238 1977 12 1977.958 333.47 333.47 334.48 31
## 239 1978 1 1978.042 335.10 335.10 335.11 26
## 240 1978 2 1978.125 335.26 335.26 334.58 28
## 241 1978 3 1978.208 336.61 336.61 334.95 30
## 242 1978 4 1978.292 337.77 337.77 335.11 28
## 243 1978 5 1978.375 338.01 338.01 335.05 30
## 244 1978 6 1978.458 337.98 337.98 335.60 27
## 245 1978 7 1978.542 336.48 336.48 335.68 31
## 246 1978 8 1978.625 334.37 334.37 335.88 31
## 247 1978 9 1978.708 332.33 332.33 335.49 24
## 248 1978 10 1978.792 332.41 332.41 335.66 31
## 249 1978 11 1978.875 333.76 333.76 336.01 28
## 250 1978 12 1978.958 334.83 334.83 335.79 31
## 251 1979 1 1979.042 336.21 336.21 336.23 31
## 252 1979 2 1979.125 336.65 336.65 335.94 27
## 253 1979 3 1979.208 338.13 338.13 336.40 29
## 254 1979 4 1979.292 338.94 338.94 336.21 24
## 255 1979 5 1979.375 339.00 339.00 336.01 20
## 256 1979 6 1979.458 339.20 339.20 336.79 30
## 257 1979 7 1979.542 337.60 337.60 336.78 29
## 258 1979 8 1979.625 335.56 335.56 337.10 29
## 259 1979 9 1979.708 333.93 333.93 337.17 26
## 260 1979 10 1979.792 334.12 334.12 337.44 28
## 261 1979 11 1979.875 335.26 335.26 337.54 30
## 262 1979 12 1979.958 336.78 336.78 337.76 24
## 263 1980 1 1980.042 337.80 337.80 337.95 30
## 264 1980 2 1980.125 338.28 338.28 337.73 28
## 265 1980 3 1980.208 340.04 340.04 338.43 31
## 266 1980 4 1980.292 340.86 340.86 338.20 30
## 267 1980 5 1980.375 341.47 341.47 338.45 31
## 268 1980 6 1980.458 341.26 341.26 338.79 30
## 269 1980 7 1980.542 339.34 339.34 338.48 26
## 270 1980 8 1980.625 337.45 337.45 338.92 24
## 271 1980 9 1980.708 336.10 336.10 339.27 26
## 272 1980 10 1980.792 336.05 336.05 339.32 29
## 273 1980 11 1980.875 337.21 337.21 339.38 27
## 274 1980 12 1980.958 338.29 338.29 339.23 27
## 275 1981 1 1981.042 339.36 339.36 339.41 31
## 276 1981 2 1981.125 340.51 340.51 339.89 28
## 277 1981 3 1981.208 341.57 341.57 339.98 31
## 278 1981 4 1981.292 342.56 342.56 339.93 30
## 279 1981 5 1981.375 343.01 343.01 339.98 31
## 280 1981 6 1981.458 342.49 342.49 340.02 28
## 281 1981 7 1981.542 340.68 340.68 339.81 31
## 282 1981 8 1981.625 338.49 338.49 339.97 28
## 283 1981 9 1981.708 336.92 336.92 340.17 30
## 284 1981 10 1981.792 337.12 337.12 340.44 31
## 285 1981 11 1981.875 338.59 338.59 340.76 30
## 286 1981 12 1981.958 339.90 339.90 340.81 28
## 287 1982 1 1982.042 340.92 340.92 341.03 30
## 288 1982 2 1982.125 341.69 341.69 341.05 27
## 289 1982 3 1982.208 342.86 342.86 341.21 21
## 290 1982 4 1982.292 343.92 343.92 341.31 9
## 291 1982 5 1982.375 344.67 344.67 341.63 31
## 292 1982 6 1982.458 343.78 343.78 341.33 29
## 293 1982 7 1982.542 342.23 342.23 341.45 28
## 294 1982 8 1982.625 340.11 340.11 341.61 27
## 295 1982 9 1982.708 338.32 338.32 341.53 26
## 296 1982 10 1982.792 338.39 338.39 341.72 30
## 297 1982 11 1982.875 339.48 339.48 341.62 28
## 298 1982 12 1982.958 340.88 340.88 341.77 31
## 299 1983 1 1983.042 341.64 341.64 341.79 31
## 300 1983 2 1983.125 342.87 342.87 342.30 27
## 301 1983 3 1983.208 343.59 343.59 342.04 30
## 302 1983 4 1983.292 345.25 345.25 342.63 26
## 303 1983 5 1983.375 345.96 345.96 342.80 29
## 304 1983 6 1983.458 345.52 345.52 343.03 25
## 305 1983 7 1983.542 344.15 344.15 343.42 30
## 306 1983 8 1983.625 342.25 342.25 343.72 28
## 307 1983 9 1983.708 340.17 340.17 343.34 27
## 308 1983 10 1983.792 340.30 340.30 343.63 31
## 309 1983 11 1983.875 341.53 341.53 343.66 29
## 310 1983 12 1983.958 343.07 343.07 343.97 26
## 311 1984 1 1984.042 344.05 344.05 344.15 31
## 312 1984 2 1984.125 344.77 344.77 344.21 29
## 313 1984 3 1984.208 345.46 345.46 344.00 23
## 314 1984 4 1984.292 NA 346.77 344.17 2
## 315 1984 5 1984.375 347.55 347.55 344.35 27
## 316 1984 6 1984.458 346.98 346.98 344.50 25
## 317 1984 7 1984.542 345.55 345.55 344.82 29
## 318 1984 8 1984.625 343.20 343.20 344.67 24
## 319 1984 9 1984.708 341.35 341.35 344.54 27
## 320 1984 10 1984.792 341.68 341.68 344.99 28
## 321 1984 11 1984.875 343.06 343.06 345.16 28
## 322 1984 12 1984.958 344.54 344.54 345.42 28
## 323 1985 1 1985.042 345.25 345.25 345.35 29
## 324 1985 2 1985.125 346.06 346.06 345.51 26
## 325 1985 3 1985.208 347.66 347.66 346.27 28
## 326 1985 4 1985.292 348.20 348.20 345.64 27
## 327 1985 5 1985.375 348.92 348.92 345.75 31
## 328 1985 6 1985.458 348.40 348.40 345.94 30
## 329 1985 7 1985.542 346.65 346.65 345.86 28
## 330 1985 8 1985.625 344.85 344.85 346.22 27
## 331 1985 9 1985.708 343.20 343.20 346.32 27
## 332 1985 10 1985.792 343.08 343.08 346.36 27
## 333 1985 11 1985.875 344.40 344.40 346.53 29
## 334 1985 12 1985.958 345.82 345.82 346.75 29
## 335 1986 1 1986.042 346.54 346.54 346.56 26
## 336 1986 2 1986.125 347.13 347.13 346.59 26
## 337 1986 3 1986.208 348.05 348.05 346.71 31
## 338 1986 4 1986.292 349.77 349.77 347.21 30
## 339 1986 5 1986.375 350.53 350.53 347.39 30
## 340 1986 6 1986.458 349.90 349.90 347.47 24
## 341 1986 7 1986.542 348.11 348.11 347.30 26
## 342 1986 8 1986.625 346.09 346.09 347.49 29
## 343 1986 9 1986.708 345.01 345.01 348.16 25
## 344 1986 10 1986.792 344.47 344.47 347.71 30
## 345 1986 11 1986.875 345.86 345.86 347.99 30
## 346 1986 12 1986.958 347.15 347.15 348.06 30
## 347 1987 1 1987.042 348.38 348.38 348.31 30
## 348 1987 2 1987.125 348.70 348.70 348.01 28
## 349 1987 3 1987.208 349.72 349.72 348.31 27
## 350 1987 4 1987.292 351.32 351.32 348.77 30
## 351 1987 5 1987.375 352.14 352.14 349.03 31
## 352 1987 6 1987.458 351.61 351.61 349.23 29
## 353 1987 7 1987.542 349.91 349.91 349.20 28
## 354 1987 8 1987.625 347.84 347.84 349.34 27
## 355 1987 9 1987.708 346.52 346.52 349.71 29
## 356 1987 10 1987.792 346.65 346.65 349.88 30
## 357 1987 11 1987.875 347.96 347.96 350.07 27
## 358 1987 12 1987.958 349.18 349.18 350.07 30
## 359 1988 1 1988.042 350.38 350.38 350.36 29
## 360 1988 2 1988.125 351.68 351.68 351.01 27
## 361 1988 3 1988.208 352.24 352.24 350.76 30
## 362 1988 4 1988.292 353.66 353.66 351.00 30
## 363 1988 5 1988.375 354.18 354.18 350.99 31
## 364 1988 6 1988.458 353.68 353.68 351.30 29
## 365 1988 7 1988.542 352.58 352.58 351.94 30
## 366 1988 8 1988.625 350.66 350.66 352.19 27
## 367 1988 9 1988.708 349.03 349.03 352.22 29
## 368 1988 10 1988.792 349.08 349.08 352.34 31
## 369 1988 11 1988.875 350.15 350.15 352.29 27
## 370 1988 12 1988.958 351.44 351.44 352.36 30
## 371 1989 1 1989.042 352.89 352.89 352.83 31
## 372 1989 2 1989.125 353.24 353.24 352.55 28
## 373 1989 3 1989.208 353.80 353.80 352.36 30
## 374 1989 4 1989.292 355.59 355.59 352.93 30
## 375 1989 5 1989.375 355.89 355.89 352.65 31
## 376 1989 6 1989.458 355.30 355.30 352.88 30
## 377 1989 7 1989.542 353.98 353.98 353.33 31
## 378 1989 8 1989.625 351.53 351.53 353.11 29
## 379 1989 9 1989.708 350.02 350.02 353.27 28
## 380 1989 10 1989.792 350.29 350.29 353.53 30
## 381 1989 11 1989.875 351.44 351.44 353.60 30
## 382 1989 12 1989.958 352.84 352.84 353.78 28
## 383 1990 1 1990.042 353.79 353.79 353.74 30
## 384 1990 2 1990.125 354.88 354.88 354.20 28
## 385 1990 3 1990.208 355.65 355.65 354.21 31
## 386 1990 4 1990.292 356.28 356.28 353.66 29
## 387 1990 5 1990.375 357.29 357.29 354.10 30
## 388 1990 6 1990.458 356.32 356.32 353.92 30
## 389 1990 7 1990.542 354.88 354.88 354.27 29
## 390 1990 8 1990.625 352.89 352.89 354.51 30
## 391 1990 9 1990.708 351.28 351.28 354.53 27
## 392 1990 10 1990.792 351.59 351.59 354.81 31
## 393 1990 11 1990.875 353.05 353.05 355.10 29
## 394 1990 12 1990.958 354.27 354.27 355.12 30
## 395 1991 1 1991.042 354.87 354.87 354.85 29
## 396 1991 2 1991.125 355.68 355.68 355.01 28
## 397 1991 3 1991.208 357.06 357.06 355.57 27
## 398 1991 4 1991.292 358.51 358.51 355.88 30
## 399 1991 5 1991.375 359.09 359.09 355.94 31
## 400 1991 6 1991.458 358.10 358.10 355.75 30
## 401 1991 7 1991.542 356.12 356.12 355.51 29
## 402 1991 8 1991.625 353.89 353.89 355.50 31
## 403 1991 9 1991.708 352.30 352.30 355.58 30
## 404 1991 10 1991.792 352.32 352.32 355.55 31
## 405 1991 11 1991.875 353.79 353.79 355.78 30
## 406 1991 12 1991.958 355.07 355.07 355.89 31
## 407 1992 1 1992.042 356.17 356.17 356.12 30
## 408 1992 2 1992.125 356.93 356.93 356.23 29
## 409 1992 3 1992.208 357.82 357.82 356.28 29
## 410 1992 4 1992.292 359.00 359.00 356.30 30
## 411 1992 5 1992.375 359.55 359.55 356.33 31
## 412 1992 6 1992.458 359.32 359.32 356.92 30
## 413 1992 7 1992.542 356.85 356.85 356.27 30
## 414 1992 8 1992.625 354.91 354.91 356.61 30
## 415 1992 9 1992.708 352.93 352.93 356.27 29
## 416 1992 10 1992.792 353.31 353.31 356.57 30
## 417 1992 11 1992.875 354.27 354.27 356.32 30
## 418 1992 12 1992.958 355.53 355.53 356.37 31
## 419 1993 1 1993.042 356.86 356.86 356.84 28
## 420 1993 2 1993.125 357.27 357.27 356.43 28
## 421 1993 3 1993.208 358.36 358.36 356.70 31
## 422 1993 4 1993.292 359.27 359.27 356.60 30
## 423 1993 5 1993.375 360.19 360.19 357.00 31
## 424 1993 6 1993.458 359.52 359.52 357.09 30
## 425 1993 7 1993.542 357.42 357.42 356.83 27
## 426 1993 8 1993.625 355.46 355.46 357.14 28
## 427 1993 9 1993.708 354.10 354.10 357.45 26
## 428 1993 10 1993.792 354.12 354.12 357.46 30
## 429 1993 11 1993.875 355.40 355.40 357.53 30
## 430 1993 12 1993.958 356.84 356.84 357.76 31
## 431 1994 1 1994.042 358.22 358.22 358.13 27
## 432 1994 2 1994.125 358.98 358.98 358.15 27
## 433 1994 3 1994.208 359.91 359.91 358.25 30
## 434 1994 4 1994.292 361.32 361.32 358.57 30
## 435 1994 5 1994.375 361.68 361.68 358.44 30
## 436 1994 6 1994.458 360.80 360.80 358.39 30
## 437 1994 7 1994.542 359.39 359.39 358.79 31
## 438 1994 8 1994.625 357.42 357.42 359.12 31
## 439 1994 9 1994.708 355.63 355.63 359.05 30
## 440 1994 10 1994.792 356.09 356.09 359.50 31
## 441 1994 11 1994.875 357.56 357.56 359.72 30
## 442 1994 12 1994.958 358.87 358.87 359.75 31
## 443 1995 1 1995.042 359.87 359.87 359.76 31
## 444 1995 2 1995.125 360.79 360.79 359.98 28
## 445 1995 3 1995.208 361.77 361.77 360.19 30
## 446 1995 4 1995.292 363.23 363.23 360.56 30
## 447 1995 5 1995.375 363.77 363.77 360.61 29
## 448 1995 6 1995.458 363.22 363.22 360.82 30
## 449 1995 7 1995.542 361.70 361.70 361.03 31
## 450 1995 8 1995.625 359.11 359.11 360.74 30
## 451 1995 9 1995.708 358.11 358.11 361.48 30
## 452 1995 10 1995.792 357.97 357.97 361.34 31
## 453 1995 11 1995.875 359.40 359.40 361.54 30
## 454 1995 12 1995.958 360.61 360.61 361.49 30
## 455 1996 1 1996.042 362.04 362.04 361.84 31
## 456 1996 2 1996.125 363.17 363.17 362.29 29
## 457 1996 3 1996.208 364.17 364.17 362.55 31
## 458 1996 4 1996.292 364.51 364.51 361.81 29
## 459 1996 5 1996.375 365.16 365.16 362.10 31
## 460 1996 6 1996.458 364.93 364.93 362.64 30
## 461 1996 7 1996.542 363.53 363.53 362.83 31
## 462 1996 8 1996.625 361.38 361.38 362.99 30
## 463 1996 9 1996.708 359.60 359.60 363.03 28
## 464 1996 10 1996.792 359.54 359.54 362.91 31
## 465 1996 11 1996.875 360.84 360.84 362.99 30
## 466 1996 12 1996.958 362.18 362.18 363.05 31
## 467 1997 1 1997.042 363.04 363.04 362.85 31
## 468 1997 2 1997.125 364.09 364.09 363.28 28
## 469 1997 3 1997.208 364.47 364.47 362.86 31
## 470 1997 4 1997.292 366.25 366.25 363.59 29
## 471 1997 5 1997.375 366.69 366.69 363.76 30
## 472 1997 6 1997.458 365.59 365.59 363.39 28
## 473 1997 7 1997.542 364.34 364.34 363.63 26
## 474 1997 8 1997.625 362.20 362.20 363.78 27
## 475 1997 9 1997.708 360.31 360.31 363.69 29
## 476 1997 10 1997.792 360.71 360.71 364.01 31
## 477 1997 11 1997.875 362.44 362.44 364.51 29
## 478 1997 12 1997.958 364.33 364.33 365.12 31
## 479 1998 1 1998.042 365.18 365.18 365.01 30
## 480 1998 2 1998.125 365.98 365.98 365.18 27
## 481 1998 3 1998.208 367.13 367.13 365.50 31
## 482 1998 4 1998.292 368.61 368.61 366.00 30
## 483 1998 5 1998.375 369.49 369.49 366.60 30
## 484 1998 6 1998.458 368.95 368.95 366.73 30
## 485 1998 7 1998.542 367.74 367.74 367.05 30
## 486 1998 8 1998.625 365.79 365.79 367.37 31
## 487 1998 9 1998.708 364.01 364.01 367.34 29
## 488 1998 10 1998.792 364.35 364.35 367.63 31
## 489 1998 11 1998.875 365.52 365.52 367.58 23
## 490 1998 12 1998.958 367.08 367.08 367.85 26
## 491 1999 1 1999.042 368.12 368.12 367.92 27
## 492 1999 2 1999.125 368.98 368.98 368.18 23
## 493 1999 3 1999.208 369.60 369.60 368.00 25
## 494 1999 4 1999.292 370.96 370.96 368.44 30
## 495 1999 5 1999.375 370.77 370.77 367.95 30
## 496 1999 6 1999.458 370.33 370.33 368.14 30
## 497 1999 7 1999.542 369.28 369.28 368.60 30
## 498 1999 8 1999.625 366.86 366.86 368.42 30
## 499 1999 9 1999.708 364.94 364.94 368.24 28
## 500 1999 10 1999.792 365.35 365.35 368.62 31
## 501 1999 11 1999.875 366.68 366.68 368.68 29
## 502 1999 12 1999.958 368.04 368.04 368.72 29
## 503 2000 1 2000.042 369.25 369.25 369.06 27
## 504 2000 2 2000.125 369.50 369.50 368.79 23
## 505 2000 3 2000.208 370.56 370.56 369.04 31
## 506 2000 4 2000.292 371.82 371.82 369.29 29
## 507 2000 5 2000.375 371.51 371.51 368.64 29
## 508 2000 6 2000.458 371.71 371.71 369.49 29
## 509 2000 7 2000.542 369.84 369.84 369.18 28
## 510 2000 8 2000.625 368.20 368.20 369.78 28
## 511 2000 9 2000.708 366.91 366.91 370.19 30
## 512 2000 10 2000.792 366.99 366.99 370.22 31
## 513 2000 11 2000.875 368.33 368.33 370.29 25
## 514 2000 12 2000.958 369.67 369.67 370.31 30
## 515 2001 1 2001.042 370.52 370.52 370.36 29
## 516 2001 2 2001.125 371.49 371.49 370.77 28
## 517 2001 3 2001.208 372.53 372.53 370.97 29
## 518 2001 4 2001.292 373.37 373.37 370.79 30
## 519 2001 5 2001.375 373.82 373.82 370.92 25
## 520 2001 6 2001.458 373.18 373.18 370.97 28
## 521 2001 7 2001.542 371.57 371.57 370.94 30
## 522 2001 8 2001.625 369.63 369.63 371.21 30
## 523 2001 9 2001.708 368.16 368.16 371.41 30
## 524 2001 10 2001.792 368.42 368.42 371.65 31
## 525 2001 11 2001.875 369.69 369.69 371.70 27
## 526 2001 12 2001.958 371.18 371.18 371.89 31
## 527 2002 1 2002.042 372.45 372.45 372.23 29
## 528 2002 2 2002.125 373.14 373.14 372.26 27
## 529 2002 3 2002.208 373.93 373.93 372.22 27
## 530 2002 4 2002.292 375.00 375.00 372.33 30
## 531 2002 5 2002.375 375.65 375.65 372.74 31
## 532 2002 6 2002.458 375.50 375.50 373.29 30
## 533 2002 7 2002.542 374.00 374.00 373.36 31
## 534 2002 8 2002.625 371.83 371.83 373.47 31
## 535 2002 9 2002.708 370.66 370.66 374.01 29
## 536 2002 10 2002.792 370.51 370.51 373.85 31
## 537 2002 11 2002.875 372.20 372.20 374.33 30
## 538 2002 12 2002.958 373.71 373.71 374.50 30
## 539 2003 1 2003.042 374.87 374.87 374.63 31
## 540 2003 2 2003.125 375.62 375.62 374.75 27
## 541 2003 3 2003.208 376.48 376.48 374.79 31
## 542 2003 4 2003.292 377.74 377.74 375.02 27
## 543 2003 5 2003.375 378.50 378.50 375.47 31
## 544 2003 6 2003.458 378.18 378.18 375.95 29
## 545 2003 7 2003.542 376.72 376.72 376.16 28
## 546 2003 8 2003.625 374.31 374.31 375.98 26
## 547 2003 9 2003.708 373.20 373.20 376.51 28
## 548 2003 10 2003.792 373.10 373.10 376.45 30
## 549 2003 11 2003.875 374.64 374.64 376.83 25
## 550 2003 12 2003.958 375.93 375.93 376.74 30
## 551 2004 1 2004.042 377.00 377.00 376.79 31
## 552 2004 2 2004.125 377.87 377.87 377.01 27
## 553 2004 3 2004.208 378.73 378.73 377.06 30
## 554 2004 4 2004.292 380.41 380.41 377.65 29
## 555 2004 5 2004.375 380.63 380.63 377.51 28
## 556 2004 6 2004.458 379.56 379.56 377.28 26
## 557 2004 7 2004.542 377.61 377.61 377.03 28
## 558 2004 8 2004.625 376.15 376.15 377.85 19
## 559 2004 9 2004.708 374.11 374.11 377.49 16
## 560 2004 10 2004.792 374.44 374.44 377.80 31
## 561 2004 11 2004.875 375.93 375.93 378.13 27
## 562 2004 12 2004.958 377.45 377.45 378.30 31
## 563 2005 1 2005.042 378.47 378.47 378.21 31
## 564 2005 2 2005.125 379.76 379.76 378.95 27
## 565 2005 3 2005.208 381.14 381.14 379.63 27
## 566 2005 4 2005.292 382.20 382.20 379.54 30
## 567 2005 5 2005.375 382.47 382.47 379.33 31
## 568 2005 6 2005.458 382.20 382.20 379.87 30
## 569 2005 7 2005.542 380.78 380.78 380.17 30
## 570 2005 8 2005.625 378.73 378.73 380.38 30
## 571 2005 9 2005.708 376.66 376.66 379.98 30
## 572 2005 10 2005.792 376.98 376.98 380.29 14
## 573 2005 11 2005.875 378.29 378.29 380.47 24
## 574 2005 12 2005.958 379.92 379.92 380.75 29
## 575 2006 1 2006.042 381.35 381.35 381.11 26
## 576 2006 2 2006.125 382.16 382.16 381.35 25
## 577 2006 3 2006.208 382.66 382.66 381.11 30
## 578 2006 4 2006.292 384.73 384.73 382.02 27
## 579 2006 5 2006.375 384.98 384.98 381.82 29
## 580 2006 6 2006.458 384.09 384.09 381.79 29
## 581 2006 7 2006.542 382.38 382.38 381.82 27
## 582 2006 8 2006.625 380.45 380.45 382.10 28
## 583 2006 9 2006.708 378.92 378.92 382.22 29
## 584 2006 10 2006.792 379.16 379.16 382.48 25
## 585 2006 11 2006.875 380.18 380.18 382.38 30
## 586 2006 12 2006.958 381.79 381.79 382.65 31
## 587 2007 1 2007.042 382.93 382.93 382.74 23
## 588 2007 2 2007.125 383.81 383.81 382.95 21
## 589 2007 3 2007.208 384.56 384.56 382.92 31
## 590 2007 4 2007.292 386.40 386.40 383.58 28
## 591 2007 5 2007.375 386.58 386.58 383.34 31
## 592 2007 6 2007.458 386.05 386.05 383.77 30
## 593 2007 7 2007.542 384.49 384.49 383.99 31
## 594 2007 8 2007.625 382.00 382.00 383.69 30
## 595 2007 9 2007.708 380.90 380.90 384.25 25
## 596 2007 10 2007.792 381.14 381.14 384.48 30
## 597 2007 11 2007.875 382.42 382.42 384.64 30
## 598 2007 12 2007.958 383.89 383.89 384.81 28
## 599 2008 1 2008.042 385.44 385.44 385.26 31
## 600 2008 2 2008.125 385.73 385.73 384.91 28
## 601 2008 3 2008.208 385.97 385.97 384.40 30
## 602 2008 4 2008.292 387.16 387.16 384.50 28
## 603 2008 5 2008.375 388.50 388.50 385.38 26
## 604 2008 6 2008.458 387.88 387.88 385.56 24
## 605 2008 7 2008.542 386.43 386.43 385.87 10
## 606 2008 8 2008.625 384.15 384.15 385.82 28
## 607 2008 9 2008.708 383.09 383.09 386.34 30
## 608 2008 10 2008.792 382.99 382.99 386.26 24
## 609 2008 11 2008.875 384.13 384.13 386.32 29
## 610 2008 12 2008.958 385.56 385.56 386.43 30
## 611 2009 1 2009.042 386.94 386.94 386.66 31
## 612 2009 2 2009.125 387.42 387.42 386.61 28
## 613 2009 3 2009.208 388.77 388.77 387.30 30
## 614 2009 4 2009.292 389.44 389.44 386.78 28
## 615 2009 5 2009.375 390.19 390.19 387.02 31
## 616 2009 6 2009.458 389.45 389.45 387.17 29
## 617 2009 7 2009.542 387.78 387.78 387.31 30
## 618 2009 8 2009.625 385.92 385.92 387.62 30
## 619 2009 9 2009.708 384.79 384.79 388.03 27
## 620 2009 10 2009.792 384.39 384.39 387.63 31
## 621 2009 11 2009.875 386.00 386.00 388.16 30
## 622 2009 12 2009.958 387.31 387.31 388.13 21
## 623 2010 1 2010.042 388.50 388.50 388.23 28
## 624 2010 2 2010.125 389.94 389.94 389.06 25
## 625 2010 3 2010.208 391.09 391.09 389.60 28
## 626 2010 4 2010.292 392.52 392.52 389.94 28
## 627 2010 5 2010.375 393.04 393.04 389.86 30
## 628 2010 6 2010.458 392.15 392.15 389.84 29
## 629 2010 7 2010.542 390.22 390.22 389.75 31
## 630 2010 8 2010.625 388.26 388.26 389.94 31
## 631 2010 9 2010.708 386.83 386.83 390.09 30
## 632 2010 10 2010.792 387.20 387.20 390.49 31
## 633 2010 11 2010.875 388.65 388.65 390.79 30
## 634 2010 12 2010.958 389.73 389.73 390.54 31
## 635 2011 1 2011.042 391.25 391.25 390.97 29
## 636 2011 2 2011.125 391.82 391.82 390.96 28
## 637 2011 3 2011.208 392.49 392.49 390.99 30
## 638 2011 4 2011.292 393.34 393.34 390.69 28
## 639 2011 5 2011.375 394.21 394.21 390.97 31
## 640 2011 6 2011.458 393.72 393.72 391.38 28
## 641 2011 7 2011.542 392.42 392.42 391.96 31
## 642 2011 8 2011.625 390.19 390.19 391.88 29
## 643 2011 9 2011.708 389.04 389.04 392.34 28
## 644 2011 10 2011.792 388.96 388.96 392.30 31
## 645 2011 11 2011.875 390.24 390.24 392.42 29
## 646 2011 12 2011.958 391.83 391.83 392.64 30
## 647 2012 1 2012.042 393.12 393.12 392.87 30
## 648 2012 2 2012.125 393.60 393.60 392.82 27
## 649 2012 3 2012.208 394.45 394.45 392.98 31
## 650 2012 4 2012.292 396.18 396.18 393.50 30
## 651 2012 5 2012.375 396.78 396.78 393.55 31
## 652 2012 6 2012.458 395.83 395.83 393.52 28
## 653 2012 7 2012.542 394.30 394.30 393.83 29
## 654 2012 8 2012.625 392.41 392.41 394.12 31
## 655 2012 9 2012.708 391.06 391.06 394.33 29
## 656 2012 10 2012.792 391.01 391.01 394.30 30
## 657 2012 11 2012.875 392.81 392.81 394.94 29
## 658 2012 12 2012.958 394.28 394.28 395.06 30
## 659 2013 1 2013.042 395.54 395.54 395.29 31
## 660 2013 2 2013.125 396.80 396.80 396.02 28
## 661 2013 3 2013.208 397.31 397.31 395.85 30
## 662 2013 4 2013.292 398.35 398.35 395.67 26
## 663 2013 5 2013.375 399.76 399.76 396.53 31
## 664 2013 6 2013.458 398.58 398.58 396.27 27
## 665 2013 7 2013.542 397.20 397.20 396.73 28
## 666 2013 8 2013.625 395.15 395.15 396.85 31
## 667 2013 9 2013.708 393.51 393.51 396.79 28
## 668 2013 10 2013.792 393.66 393.66 396.95 30
## 669 2013 11 2013.875 395.11 395.11 397.24 30
## 670 2013 12 2013.958 396.81 396.81 397.59 31
## 671 2014 1 2014.042 397.80 397.80 397.54 31
## 672 2014 2 2014.125 397.91 397.91 397.13 27
## 673 2014 3 2014.208 399.58 399.58 398.11 28
## 674 2014 4 2014.292 401.29 401.29 398.61 30
## 675 2014 5 2014.375 401.78 401.78 398.54 24
## 676 2014 6 2014.458 401.15 401.15 398.84 30
## 677 2014 7 2014.542 399.00 399.00 398.53 31
## 678 2014 8 2014.625 397.01 397.01 398.71 25
## 679 2014 9 2014.708 395.26 395.26 398.54 29
## 680 2014 10 2014.792 395.93 395.93 399.22 26
## 681 2014 11 2014.875 397.13 397.13 399.26 28
## 682 2014 12 2014.958 398.78 398.78 399.56 30
## 683 2015 1 2015.042 399.96 399.96 399.71 30
## 684 2015 2 2015.125 400.26 400.26 399.48 27
## 685 2015 3 2015.208 401.52 401.52 400.06 24
## 686 2015 4 2015.292 403.26 403.26 400.58 27
## 687 2015 5 2015.375 403.94 403.94 400.71 30
## 688 2015 6 2015.458 402.80 402.80 400.49 28
## 689 2015 7 2015.542 401.30 401.30 400.84 233 From Antarctic ice cores (last 800 000 years)
fancy_scientific <- function(l) {
# turn in to character string in scientific notation
l <- format(l, scientific = TRUE)
# quote the part before the exponent to keep all the digits
l <- gsub("^(.*)e", "'\\1'e", l)
# turn the 'e+' into plotmath format
l <- gsub("e", "%*%10^", l)
# return this as an expression
parse(text = l)
}
groupstr <- function(x) {
#x <- sprintf("%d", x)
prettyNum(sprintf("%d", x), big.mark = ",")
}
#vostok99 <- read.table("ftp://ftp.ncdc.noaa.gov/pub/data/paleo/icecore/antarctica/vostok/co2nat.txt", skip = 1)
# metadata ftp://ftp.ncdc.noaa.gov/pub/data/paleo/icecore/antarctica/vostok/readme_petit1999.txt
#colnames(vostok99) <- c("year.bp", "co2.ppm")
# epica2010 <- read.table("ftp://ftp.ncdc.noaa.gov/pub/data/paleo/icecore/antarctica/epica_domec/edc2010d13co2.txt", sep = "", skip = 117, na.strings = "", fill = TRUE, header = TRUE, nrows = 63, colClasses = "real")
#epica <- read.table("ftp://ftp.ncdc.noaa.gov/pub/data/paleo/icecore/antarctica/epica_domec/dc_co2_hol_fl02.txt", skip = 3)
#colnames(epica) <- c("depth", "age", "CO2", "uncertainty")
# Composite CO2 record (0-800 kyr BP)
#
# 0-22 kyr Dome C (Monnin et al. 2001) measured at University of Bern
# 22-393 kyr Vostok (Petit et al. 1999; Pepin et al. 2001; Raynaud et al. 2005) measured at LGGE in Grenoble
# 393-416 kyr Dome C (Siegenthaler et al. 2005) measured at LGGE in Grenoble
# 416-664 kyr Dome C (Siegenthaler et al. 2005) measured at University of Bern
# 664-800 kyr Dome C (Luethi et al. (sub)) measured at University of Bern
# Timescale EDC3_gas_a
icecore.1 <- read.table("ftp://ftp.ncdc.noaa.gov/pub/data/paleo/icecore/antarctica/epica_domec/edc-co2-2008.txt", header = TRUE, skip = 773)
# Historical CO2 Record from the Siple Station Ice Core
# Metadata http://cdiac.ornl.gov/trends/co2/siple.html
icecore.2 <- read.table("http://cdiac.ornl.gov/ftp/trends/co2/siple2.013", skip = 36, nrows = 22)
colnames(icecore.2) <- c("depth", "airyearAD", "co2")
icecore.2$bp <- 1950 - icecore.2$airyearAD # add before present years
# Use data from Mauna Loa (cf. above)
m$bp <- 1950 - m$year # add before present years
i <- ggplot(icecore.1, aes(Age.yrBP., CO2.ppmv.)) +
geom_line(colour = "blue") +
geom_line(data = icecore.2, aes(bp, co2), colour = "blue") +
geom_line(data = m, aes(bp, year_mean), colour = "red") +
#geom_line(data=co2, aes(bp, average)) +
#geom_line(data=epica2010, aes(EDC3gasage, CO2)) +
#geom_line(data=epica, aes(age, CO2)) +
#geom_point(color="red", size=1) +
#stat_smooth(method = "loess", span = 0.2) +
xlab("Gas age (year BP)") +
scale_x_reverse(labels = fancy_scientific) + #labels = comma) +
ylab(expression(paste(CO[2], " (ppmv)"))) +
ggtitle(expression(atop(paste(CO[2], " concentration "), "Antarctic ice cores & current data from Mauna Loa"))) +
theme_bw()
i
#ggsave(i, file = "co2_800_ky.svg", width = 20, height = 10, units = "cm")
icecore.1## Age.yrBP. CO2.ppmv.
## 1 137 280.4
## 2 268 274.9
## 3 279 277.9
## 4 395 279.1
## 5 404 281.9
## 6 485 277.7
## 7 559 281.1
## 8 672 282.2
## 9 754 280.1
## 10 877 278.4
## 11 950 276.6
## 12 1060 279.1
## 13 1153 277.7
## 14 1233 278.7
## 15 1350 277.4
## 16 1453 279.2
## 17 1552 280.0
## 18 1638 278.9
## 19 1733 278.7
## 20 1812 278.0
## 21 1931 276.9
## 22 2057 276.7
## 23 2128 276.7
## 24 2212 277.6
## 25 2334 277.9
## 26 2433 273.9
## 27 2536 278.9
## 28 2604 275.3
## 29 2728 274.7
## 30 2806 276.3
## 31 2902 274.6
## 32 3053 276.3
## 33 3116 273.1
## 34 3215 274.0
## 35 3336 275.0
## 36 3453 273.4
## 37 3523 273.0
## 38 3622 271.5
## 39 3721 275.4
## 40 3790 274.9
## 41 3910 271.7
## 42 4004 271.6
## 43 4096 272.8
## 44 4161 271.5
## 45 4324 271.1
## 46 4374 269.1
## 47 4480 269.8
## 48 4573 271.5
## 49 4703 270.7
## 50 4766 269.3
## 51 4874 268.6
## 52 5004 269.8
## 53 5094 267.6
## 54 5160 265.3
## 55 5274 265.2
## 56 5370 267.6
## 57 5476 265.9
## 58 5562 265.5
## 59 5657 260.7
## 60 5716 266.7
## 61 5855 265.5
## 62 5998 263.2
## 63 6039 262.7
## 64 6131 261.2
## 65 6263 261.1
## 66 6354 259.4
## 67 6434 262.1
## 68 6545 262.9
## 69 6617 258.1
## 70 6713 257.6
## 71 6838 262.3
## 72 6941 263.0
## 73 7028 260.7
## 74 7112 258.4
## 75 7234 260.1
## 76 7320 260.4
## 77 7413 259.7
## 78 7507 259.2
## 79 7590 260.8
## 80 7691 259.6
## 81 7781 259.3
## 82 7876 258.3
## 83 7990 261.3
## 84 8050 260.7
## 85 8181 261.8
## 86 8281 259.0
## 87 8387 260.9
## 88 8477 260.4
## 89 8579 259.3
## 90 8653 262.0
## 91 8784 263.7
## 92 8869 263.8
## 93 8973 265.2
## 94 9092 260.6
## 95 9140 260.9
## 96 9232 263.0
## 97 9317 263.8
## 98 9536 264.4
## 99 9597 264.2
## 100 9721 264.0
## 101 9807 263.4
## 102 9909 265.7
## 103 9983 264.9
## 104 10088 267.5
## 105 10209 266.9
## 106 10294 266.0
## 107 10417 265.1
## 108 10527 267.6
## 109 10621 264.8
## 110 10744 264.8
## 111 10805 265.0
## 112 10827 265.3
## 113 10895 264.4
## 114 10933 264.1
## 115 11014 264.2
## 116 11087 264.5
## 117 11136 264.0
## 118 11201 263.0
## 119 11236 265.2
## 120 11278 258.8
## 121 11338 260.8
## 122 11392 255.4
## 123 11436 253.9
## 124 11469 253.8
## 125 11580 250.7
## 126 11635 249.7
## 127 11676 251.1
## 128 11727 250.7
## 129 11819 245.3
## 130 11896 245.3
## 131 11958 246.6
## 132 12050 243.2
## 133 12122 240.3
## 134 12371 237.5
## 135 12496 237.6
## 136 12642 234.2
## 137 12760 238.3
## 138 12942 237.3
## 139 13090 237.9
## 140 13241 237.6
## 141 13440 236.4
## 142 13542 239.2
## 143 13653 238.6
## 144 13804 238.6
## 145 13948 239.1
## 146 14303 228.5
## 147 14550 228.4
## 148 14725 226.1
## 149 14890 225.2
## 150 15012 224.5
## 151 15233 222.0
## 152 15438 221.0
## 153 15570 220.9
## 154 15742 219.4
## 155 15886 214.0
## 156 16073 207.5
## 157 16260 207.7
## 158 16452 202.9
## 159 16659 200.8
## 160 16870 195.2
## 161 17111 193.9
## 162 17375 191.0
## 163 17565 188.5
## 164 17809 188.5
## 165 17943 189.2
## 166 18285 187.0
## 167 18541 188.6
## 168 18828 189.4
## 169 18868 192.3
## 170 18921 188.3
## 171 19347 188.7
## 172 19509 188.8
## 173 19748 190.0
## 174 19988 188.0
## 175 20168 188.2
## 176 20197 195.0
## 177 20502 187.8
## 178 20748 186.9
## 179 21011 186.5
## 180 21257 184.7
## 181 21507 186.1
## 182 21854 185.7
## 183 22015 184.4
## 184 22827 189.2
## 185 25994 191.6
## 186 29063 188.5
## 187 30020 191.7
## 188 35009 205.3
## 189 37471 209.1
## 190 43500 209.1
## 191 47336 189.3
## 192 48854 188.4
## 193 49690 210.1
## 194 50663 215.7
## 195 52382 190.4
## 196 57088 221.7
## 197 57657 210.4
## 198 62859 195.3
## 199 64939 191.4
## 200 65939 194.9
## 201 71049 227.3
## 202 73227 229.1
## 203 77150 217.1
## 204 78183 221.7
## 205 80614 230.9
## 206 82417 241.1
## 207 83333 236.4
## 208 84016 228.0
## 209 85020 214.2
## 210 86181 217.0
## 211 87917 208.0
## 212 90357 224.2
## 213 91249 228.3
## 214 94353 232.0
## 215 99849 225.8
## 216 100837 230.8
## 217 101749 236.9
## 218 103465 228.1
## 219 104704 236.9
## 220 105636 230.6
## 221 107579 238.2
## 222 108153 245.6
## 223 109804 251.2
## 224 111862 256.7
## 225 113262 266.3
## 226 114096 261.4
## 227 114601 274.5
## 228 115118 273.2
## 229 116501 262.5
## 230 117750 267.6
## 231 118649 273.7
## 232 119672 271.9
## 233 120382 265.2
## 234 121017 277.6
## 235 122344 272.1
## 236 123070 276.4
## 237 124213 268.7
## 238 124257 266.6
## 239 124789 266.3
## 240 125081 279.7
## 241 125262 273.0
## 242 125434 277.1
## 243 126347 273.7
## 244 126598 267.1
## 245 126886 262.5
## 246 127132 262.6
## 247 127622 275.3
## 248 127907 275.6
## 249 128344 274.0
## 250 128372 287.1
## 251 128609 286.8
## 252 128866 282.6
## 253 129146 264.1
## 254 129340 263.4
## 255 129652 257.9
## 256 129736 259.0
## 257 131329 245.0
## 258 131455 240.4
## 259 131728 228.9
## 260 132492 223.5
## 261 133069 223.9
## 262 133105 220.3
## 263 133427 210.6
## 264 134123 208.9
## 265 134287 203.7
## 266 134960 204.5
## 267 135114 200.4
## 268 135207 198.0
## 269 135603 198.0
## 270 135883 201.7
## 271 136011 200.7
## 272 136251 202.4
## 273 136567 195.8
## 274 136655 201.1
## 275 137293 194.3
## 276 137633 193.4
## 277 137982 194.2
## 278 138185 190.2
## 279 139275 192.3
## 280 139617 196.5
## 281 140899 195.6
## 282 140960 196.4
## 283 142058 190.4
## 284 145363 196.9
## 285 148831 203.0
## 286 149803 191.9
## 287 149921 188.9
## 288 151423 200.6
## 289 154480 189.0
## 290 155395 185.5
## 291 155813 187.5
## 292 159562 204.3
## 293 159943 196.5
## 294 161679 191.6
## 295 162228 190.1
## 296 163024 186.7
## 297 163698 183.8
## 298 164439 196.6
## 299 167183 197.8
## 300 169492 197.7
## 301 171351 196.0
## 302 172434 190.3
## 303 173135 189.4
## 304 173394 190.1
## 305 175306 207.7
## 306 177139 213.2
## 307 178179 217.7
## 308 179117 198.1
## 309 180068 199.7
## 310 182046 203.4
## 311 184685 210.7
## 312 186697 231.3
## 313 188436 231.4
## 314 189076 220.3
## 315 190352 218.0
## 316 192910 226.5
## 317 193481 220.0
## 318 195017 226.4
## 319 195673 241.2
## 320 196179 242.6
## 321 198972 250.9
## 322 199918 239.1
## 323 201163 247.6
## 324 202413 244.4
## 325 202874 231.9
## 326 203317 232.2
## 327 203837 228.6
## 328 204480 226.3
## 329 204908 229.4
## 330 205362 231.4
## 331 205952 238.1
## 332 206144 237.2
## 333 206810 230.0
## 334 207120 240.5
## 335 207544 242.2
## 336 208064 244.6
## 337 208254 243.9
## 338 208803 247.2
## 339 208995 252.0
## 340 209432 246.9
## 341 209817 239.5
## 342 210154 257.4
## 343 210813 243.4
## 344 211858 251.2
## 345 212716 241.4
## 346 213270 240.3
## 347 213536 242.6
## 348 213984 247.5
## 349 214459 251.7
## 350 214794 251.1
## 351 215382 245.3
## 352 216352 240.5
## 353 218361 214.1
## 354 219116 216.1
## 355 220058 207.1
## 356 220739 208.8
## 357 222030 205.6
## 358 224269 203.3
## 359 224923 215.7
## 360 225909 235.5
## 361 226711 234.5
## 362 227027 233.1
## 363 227776 224.5
## 364 229423 232.4
## 365 230422 233.9
## 366 231066 241.6
## 367 234817 245.2
## 368 235480 252.1
## 369 236114 241.4
## 370 236734 247.4
## 371 237294 243.1
## 372 237868 239.1
## 373 238558 245.6
## 374 239013 245.8
## 375 239477 247.4
## 376 239973 252.8
## 377 240945 259.7
## 378 241366 263.2
## 379 242007 279.0
## 380 242346 280.2
## 381 243071 263.7
## 382 243429 252.3
## 383 243856 249.9
## 384 244347 236.7
## 385 244864 230.4
## 386 245441 219.4
## 387 247681 214.7
## 388 250133 200.2
## 389 251005 213.9
## 390 251864 195.4
## 391 252739 196.7
## 392 253636 195.4
## 393 255498 199.0
## 394 256309 201.9
## 395 257100 204.0
## 396 258661 203.9
## 397 260353 209.6
## 398 260916 205.7
## 399 262092 208.9
## 400 262930 214.6
## 401 264509 228.1
## 402 265653 199.9
## 403 266326 211.7
## 404 267442 188.7
## 405 268181 187.2
## 406 269154 194.2
## 407 270222 198.8
## 408 271256 184.7
## 409 272311 190.4
## 410 273310 193.9
## 411 274321 194.1
## 412 275350 198.4
## 413 276326 193.2
## 414 277361 202.2
## 415 278350 204.5
## 416 279320 211.0
## 417 280269 215.3
## 418 281758 223.7
## 419 283785 231.3
## 420 285380 228.0
## 421 286128 226.4
## 422 287036 231.4
## 423 287538 230.4
## 424 288233 231.0
## 425 290153 234.9
## 426 290730 220.4
## 427 291367 217.1
## 428 292405 207.6
## 429 293057 206.0
## 430 293779 206.7
## 431 294568 212.7
## 432 295930 213.1
## 433 297921 217.1
## 434 298741 224.4
## 435 299621 231.0
## 436 300472 236.1
## 437 301402 239.0
## 438 302280 236.0
## 439 303226 240.2
## 440 304232 240.7
## 441 304901 250.2
## 442 305655 248.6
## 443 306547 244.8
## 444 307495 225.8
## 445 308744 227.8
## 446 309389 226.2
## 447 310168 233.2
## 448 311013 237.8
## 449 311868 239.0
## 450 312676 241.9
## 451 313414 251.6
## 452 314144 256.7
## 453 314867 257.1
## 454 315572 246.8
## 455 316200 272.6
## 456 316897 251.6
## 457 317734 245.2
## 458 318559 233.4
## 459 319480 255.8
## 460 320358 249.2
## 461 322015 257.2
## 462 322767 260.4
## 463 323526 260.3
## 464 324280 260.5
## 465 324971 266.2
## 466 325720 264.0
## 467 326321 266.1
## 468 326972 270.1
## 469 327590 271.9
## 470 328221 275.1
## 471 329475 265.0
## 472 330208 271.7
## 473 330740 272.6
## 474 331438 273.1
## 475 331944 282.4
## 476 332285 289.1
## 477 332462 288.4
## 478 332919 298.6
## 479 333380 278.1
## 480 333890 285.8
## 481 334261 278.6
## 482 334748 270.5
## 483 335287 255.7
## 484 335918 241.9
## 485 336725 239.6
## 486 337391 234.2
## 487 339298 250.1
## 488 340456 200.7
## 489 341802 205.2
## 490 343282 204.8
## 491 344446 211.9
## 492 345980 220.3
## 493 348298 221.1
## 494 349688 216.2
## 495 350925 209.4
## 496 352275 209.2
## 497 355302 193.0
## 498 356898 186.1
## 499 358712 185.8
## 500 360957 201.2
## 501 363385 206.3
## 502 366235 201.9
## 503 367856 199.9
## 504 369446 214.7
## 505 371090 224.6
## 506 372646 229.6
## 507 376568 227.0
## 508 378096 240.0
## 509 379696 239.1
## 510 381132 246.8
## 511 382670 245.8
## 512 384534 258.1
## 513 385398 264.6
## 514 386144 259.2
## 515 388062 255.2
## 516 389946 250.1
## 517 391896 266.3
## 518 392544 259.5
## 519 393579 273.6
## 520 394560 260.7
## 521 398086 276.3
## 522 399722 277.1
## 523 400504 283.2
## 524 402000 283.1
## 525 402731 275.7
## 526 404181 276.5
## 527 404927 280.5
## 528 406368 279.6
## 529 407093 285.6
## 530 408600 284.5
## 531 409383 275.2
## 532 409383 274.2
## 533 410206 282.6
## 534 411071 283.5
## 535 412962 274.9
## 536 413948 264.9
## 537 414963 271.6
## 538 415717 276.4
## 539 416193 271.7
## 540 417191 273.4
## 541 417698 271.8
## 542 418245 274.6
## 543 419260 273.7
## 544 419808 271.2
## 545 420350 273.8
## 546 421484 268.6
## 547 422074 266.4
## 548 422649 270.6
## 549 423764 267.7
## 550 424332 268.3
## 551 424840 270.8
## 552 425242 270.0
## 553 425569 265.4
## 554 425975 255.3
## 555 426598 252.1
## 556 427285 248.2
## 557 427566 242.5
## 558 429006 219.7
## 559 429876 227.2
## 560 431063 211.5
## 561 431445 207.6
## 562 432599 209.8
## 563 433674 207.5
## 564 434804 200.3
## 565 435989 201.7
## 566 437152 201.3
## 567 438356 202.0
## 568 439565 199.1
## 569 441220 201.1
## 570 442411 203.5
## 571 443562 208.1
## 572 444650 201.7
## 573 445829 201.2
## 574 446984 204.9
## 575 448103 201.9
## 576 449244 198.4
## 577 450455 193.3
## 578 451593 192.5
## 579 452283 199.1
## 580 452795 204.3
## 581 453514 203.3
## 582 454023 208.3
## 583 455279 202.4
## 584 456591 195.5
## 585 458049 190.7
## 586 459430 194.4
## 587 460792 199.9
## 588 461687 205.2
## 589 462133 210.0
## 590 463456 208.1
## 591 464866 204.4
## 592 466265 203.4
## 593 467602 205.5
## 594 468323 206.5
## 595 468810 215.5
## 596 469470 218.7
## 597 469941 229.2
## 598 470597 232.7
## 599 471046 243.7
## 600 471763 243.9
## 601 472095 245.6
## 602 473102 241.2
## 603 474162 233.1
## 604 475218 232.8
## 605 476240 236.6
## 606 477162 239.1
## 607 477719 236.5
## 608 478124 231.3
## 609 478721 220.8
## 610 479144 218.8
## 611 480108 223.8
## 612 481035 223.4
## 613 481932 227.3
## 614 482803 231.4
## 615 483649 231.4
## 616 484465 231.2
## 617 485276 233.0
## 618 486043 232.9
## 619 486816 231.3
## 620 487541 237.1
## 621 488268 242.0
## 622 488689 240.9
## 623 488996 252.8
## 624 489475 249.3
## 625 489722 242.6
## 626 490445 243.3
## 627 491191 246.7
## 628 491959 243.7
## 629 492738 239.7
## 630 493507 239.1
## 631 494311 238.3
## 632 495085 238.5
## 633 495853 237.6
## 634 496608 236.4
## 635 497366 236.5
## 636 498104 232.8
## 637 498865 230.6
## 638 499610 230.8
## 639 500363 228.2
## 640 501130 230.8
## 641 501866 231.2
## 642 502625 232.0
## 643 503375 232.2
## 644 504083 234.1
## 645 504777 233.9
## 646 505516 235.5
## 647 506287 236.0
## 648 507011 238.2
## 649 507714 237.0
## 650 508476 240.3
## 651 509247 233.4
## 652 510008 236.4
## 653 510765 235.2
## 654 511515 242.0
## 655 512269 238.2
## 656 512997 239.2
## 657 513724 241.5
## 658 514429 243.0
## 659 515113 247.2
## 660 515771 246.2
## 661 516417 245.5
## 662 517069 245.3
## 663 517706 247.7
## 664 518334 245.5
## 665 518965 243.5
## 666 520892 241.9
## 667 521552 236.8
## 668 522216 235.8
## 669 522923 233.7
## 670 523666 230.3
## 671 524409 227.8
## 672 525164 224.7
## 673 525908 222.7
## 674 526661 221.2
## 675 527454 220.3
## 676 528262 220.7
## 677 528747 214.5
## 678 529101 211.4
## 679 529640 204.4
## 680 530016 200.0
## 681 530685 195.1
## 682 531035 193.8
## 683 531662 190.5
## 684 532119 199.8
## 685 533246 199.4
## 686 534364 205.0
## 687 535483 206.4
## 688 536622 212.9
## 689 537770 211.2
## 690 538959 205.0
## 691 540154 203.6
## 692 541345 208.1
## 693 542502 210.3
## 694 543676 209.6
## 695 544804 209.7
## 696 545956 204.6
## 697 547013 202.5
## 698 547573 208.8
## 699 547972 213.6
## 700 548933 215.2
## 701 549830 219.9
## 702 550637 221.7
## 703 551124 222.3
## 704 551452 226.3
## 705 551979 222.1
## 706 552250 220.3
## 707 552961 223.6
## 708 553626 224.3
## 709 554254 226.7
## 710 554831 233.9
## 711 555345 241.0
## 712 555637 243.0
## 713 555839 249.1
## 714 556143 244.0
## 715 556364 250.5
## 716 556724 245.6
## 717 556889 240.4
## 718 557385 238.1
## 719 557887 234.5
## 720 558183 229.7
## 721 558393 230.2
## 722 558929 230.3
## 723 559461 228.3
## 724 559972 231.6
## 725 560459 231.9
## 726 560930 234.6
## 727 561309 234.1
## 728 561846 240.1
## 729 562272 242.3
## 730 562690 245.7
## 731 563115 245.8
## 732 563517 247.6
## 733 563918 251.4
## 734 564311 252.4
## 735 564699 252.6
## 736 565077 251.4
## 737 565466 253.7
## 738 565851 254.3
## 739 566233 253.9
## 740 566619 254.5
## 741 567007 253.2
## 742 567381 253.9
## 743 567777 252.8
## 744 568167 253.0
## 745 568556 250.2
## 746 568946 251.3
## 747 569331 250.0
## 748 569728 251.3
## 749 570139 251.8
## 750 570548 251.8
## 751 570955 249.6
## 752 571367 251.6
## 753 571775 250.3
## 754 572190 246.3
## 755 572619 247.7
## 756 573037 249.2
## 757 573462 248.7
## 758 573902 251.8
## 759 574353 251.9
## 760 574803 252.1
## 761 575265 248.9
## 762 575730 252.5
## 763 576222 252.9
## 764 576723 251.3
## 765 577242 251.3
## 766 577769 251.1
## 767 578305 249.1
## 768 578857 252.9
## 769 579435 251.5
## 770 579821 244.1
## 771 580019 243.8
## 772 580400 236.2
## 773 580695 230.3
## 774 581160 225.2
## 775 581507 219.4
## 776 582000 215.4
## 777 582391 209.9
## 778 583034 206.7
## 779 583413 210.6
## 780 586406 224.0
## 781 586893 226.0
## 782 587647 234.4
## 783 588034 238.8
## 784 588716 246.3
## 785 589224 250.2
## 786 589868 248.1
## 787 590335 243.6
## 788 590893 237.4
## 789 592279 225.7
## 790 592697 229.4
## 791 593415 233.2
## 792 593870 238.1
## 793 594639 238.0
## 794 595189 232.9
## 795 595950 221.2
## 796 596515 216.4
## 797 597443 216.3
## 798 597966 219.0
## 799 598833 226.0
## 800 599434 229.4
## 801 600882 232.5
## 802 602278 237.9
## 803 603673 239.1
## 804 605076 244.5
## 805 606374 248.5
## 806 607708 254.6
## 807 608929 259.2
## 808 610136 257.7
## 809 611319 259.6
## 810 612514 258.0
## 811 613709 256.0
## 812 614923 252.5
## 813 616227 252.4
## 814 617544 252.3
## 815 618955 247.9
## 816 620338 245.4
## 817 621760 243.0
## 818 623214 238.6
## 819 624690 237.3
## 820 625570 234.3
## 821 626215 227.8
## 822 627211 214.9
## 823 627999 205.0
## 824 629446 198.7
## 825 630295 199.9
## 826 632001 199.5
## 827 632856 195.6
## 828 634592 192.1
## 829 635489 193.9
## 830 636957 193.0
## 831 638127 194.3
## 832 639848 187.7
## 833 640760 189.1
## 834 642405 190.2
## 835 643204 190.4
## 836 644662 194.1
## 837 645470 194.6
## 838 646723 188.6
## 839 647588 190.5
## 840 648946 192.1
## 841 649598 194.8
## 842 650891 187.8
## 843 651556 192.0
## 844 652802 185.8
## 845 653448 185.3
## 846 654473 191.3
## 847 655202 193.2
## 848 656286 198.4
## 849 656882 198.7
## 850 657923 191.3
## 851 658474 189.2
## 852 659524 182.4
## 853 660084 183.9
## 854 660981 186.6
## 855 661650 189.8
## 856 662669 187.8
## 857 663206 189.0
## 858 664192 189.1
## 859 664690 185.8
## 860 665645 178.2
## 861 666174 178.5
## 862 666995 172.7
## 863 667569 171.6
## 864 668447 175.6
## 865 668934 178.5
## 866 669751 185.5
## 867 670176 189.3
## 868 670916 192.6
## 869 671286 194.6
## 870 671926 192.5
## 871 672386 194.8
## 872 673095 198.7
## 873 673475 202.3
## 874 674173 209.2
## 875 674521 213.7
## 876 675589 220.0
## 877 676170 217.4
## 878 676524 214.0
## 879 677068 215.6
## 880 677380 217.4
## 881 677984 226.5
## 882 678291 230.0
## 883 678914 230.8
## 884 679252 230.3
## 885 679732 223.0
## 886 680081 225.0
## 887 680608 218.4
## 888 680920 216.2
## 889 681531 217.3
## 890 681849 219.7
## 891 682446 222.6
## 892 682762 221.3
## 893 683262 218.7
## 894 683612 217.5
## 895 684142 217.5
## 896 684419 219.4
## 897 684949 219.5
## 898 685217 219.7
## 899 685716 223.5
## 900 685980 224.4
## 901 686425 224.3
## 902 686726 227.3
## 903 687227 228.7
## 904 687470 228.8
## 905 687949 232.3
## 906 688195 233.0
## 907 688652 235.8
## 908 688881 234.9
## 909 689304 235.8
## 910 689603 239.9
## 911 690057 235.6
## 912 690294 234.8
## 913 690771 240.0
## 914 690999 234.0
## 915 691432 235.6
## 916 691672 236.6
## 917 692057 235.5
## 918 692338 235.9
## 919 692765 232.9
## 920 692990 234.8
## 921 693433 234.1
## 922 693639 234.1
## 923 694068 235.4
## 924 694300 234.4
## 925 694672 235.8
## 926 694941 236.9
## 927 695383 236.7
## 928 695618 237.6
## 929 696062 234.6
## 930 696287 234.8
## 931 696717 235.4
## 932 696963 232.9
## 933 697331 232.4
## 934 697660 232.2
## 935 698120 231.0
## 936 698386 228.0
## 937 698880 226.4
## 938 699102 227.7
## 939 699633 227.6
## 940 699909 226.0
## 941 700389 228.7
## 942 700744 227.6
## 943 701242 228.2
## 944 701606 228.8
## 945 702158 229.3
## 946 702381 230.3
## 947 702988 232.4
## 948 703384 232.8
## 949 703934 231.2
## 950 704211 231.2
## 951 704973 233.9
## 952 705410 228.5
## 953 706141 228.6
## 954 706409 228.1
## 955 709680 219.7
## 956 710237 220.2
## 957 711104 218.9
## 958 711421 218.7
## 959 712833 218.5
## 960 713843 222.3
## 961 714288 221.7
## 962 715180 217.7
## 963 715750 222.1
## 964 716704 208.1
## 965 717065 205.0
## 966 718064 192.9
## 967 718779 184.3
## 968 719794 187.8
## 969 720337 195.2
## 970 721318 206.3
## 971 722039 211.3
## 972 723133 213.4
## 973 723464 213.2
## 974 724445 210.3
## 975 725059 206.3
## 976 725937 206.6
## 977 726322 208.1
## 978 727192 204.5
## 979 727732 205.5
## 980 728568 205.1
## 981 728921 205.4
## 982 730299 204.0
## 983 730949 204.1
## 984 731329 204.9
## 985 732141 203.8
## 986 732548 205.9
## 987 733261 203.9
## 988 733543 204.4
## 989 734223 205.3
## 990 734604 207.4
## 991 735167 201.7
## 992 735459 201.3
## 993 736079 204.1
## 994 736460 209.7
## 995 737034 208.1
## 996 737266 211.3
## 997 737889 210.6
## 998 738272 206.5
## 999 738844 200.3
## 1000 739137 197.8
## 1001 739842 185.0
## 1002 740277 178.9
## 1003 741021 178.6
## 1004 741307 176.6
## 1005 742075 177.9
## 1006 742558 182.3
## 1007 743233 181.7
## 1008 743589 180.4
## 1009 744354 181.2
## 1010 744843 182.0
## 1011 745836 180.2
## 1012 746643 183.1
## 1013 747110 180.8
## 1014 747759 180.8
## 1015 748101 180.4
## 1016 748913 181.4
## 1017 749401 182.9
## 1018 750194 184.8
## 1019 750476 184.0
## 1020 751301 186.2
## 1021 751830 190.1
## 1022 752605 188.0
## 1023 752978 196.5
## 1024 753984 191.5
## 1025 754394 195.8
## 1026 755241 200.5
## 1027 755495 201.9
## 1028 756357 206.0
## 1029 756871 210.7
## 1030 757561 215.0
## 1031 757962 219.2
## 1032 758259 215.3
## 1033 758475 217.6
## 1034 758883 217.9
## 1035 759306 214.9
## 1036 760202 208.8
## 1037 760500 206.9
## 1038 761453 206.1
## 1039 762064 212.0
## 1040 762871 216.0
## 1041 763319 215.0
## 1042 763686 225.7
## 1043 763934 221.4
## 1044 764249 224.1
## 1045 764833 222.6
## 1046 765735 216.8
## 1047 766110 215.7
## 1048 767096 213.5
## 1049 767737 217.7
## 1050 768655 222.4
## 1051 769050 221.4
## 1052 769514 228.1
## 1053 769787 226.0
## 1054 770133 231.5
## 1055 770447 235.6
## 1056 770724 238.2
## 1057 770994 235.6
## 1058 771260 237.7
## 1059 771708 229.5
## 1060 772076 230.9
## 1061 773043 230.0
## 1062 773680 233.3
## 1063 774548 236.9
## 1064 775051 240.5
## 1065 776037 238.8
## 1066 776636 243.9
## 1067 777556 246.3
## 1068 777900 246.1
## 1069 778798 245.8
## 1070 779364 243.9
## 1071 780268 248.5
## 1072 780612 242.8
## 1073 781535 248.7
## 1074 782059 246.2
## 1075 783263 252.1
## 1076 784147 246.9
## 1077 784692 250.7
## 1078 785413 253.3
## 1079 785789 255.4
## 1080 786587 260.3
## 1081 787076 256.9
## 1082 787805 248.1
## 1083 788033 247.9
## 1084 788887 229.5
## 1085 789381 226.3
## 1086 790153 218.2
## 1087 790538 221.3
## 1088 791491 215.4
## 1089 792081 205.1
## 1090 792943 204.0
## 1091 793260 209.0
## 1092 794608 199.4
## 1093 795202 195.2
## 1094 796467 189.3
## 1095 797099 188.4
## 1096 798512 191.0icecore.2## depth airyearAD co2 bp
## 1 187.70 1744 276.8 206
## 2 177.50 1764 276.7 186
## 3 168.30 1791 279.7 159
## 4 154.89 1816 283.8 134
## 5 142.75 1839 283.1 111
## 6 140.75 1843 287.4 107
## 7 138.20 1847 286.8 103
## 8 134.47 1854 288.2 96
## 9 126.80 1869 289.3 81
## 10 123.80 1874 289.5 76
## 11 121.80 1878 290.3 72
## 12 116.82 1887 292.3 63
## 13 110.20 1899 295.8 51
## 14 108.80 1903 294.8 47
## 15 107.20 1905 296.9 45
## 16 105.25 1909 299.2 41
## 17 101.80 1915 300.5 35
## 18 98.80 1921 301.6 29
## 19 95.17 1927 305.5 23
## 20 90.77 1935 306.6 15
## 21 86.80 1943 307.9 7
## 22 81.22 1953 312.7 -3m## year year_mean bp
## 1 1958 315.2410 -8
## 2 1959 315.9742 -9
## 3 1960 316.9075 -10
## 4 1961 317.6375 -11
## 5 1962 318.4508 -12
## 6 1963 318.9942 -13
## 7 1964 319.6175 -14
## 8 1965 320.0442 -15
## 9 1966 321.3833 -16
## 10 1967 322.1575 -17
## 11 1968 323.0450 -18
## 12 1969 324.6242 -19
## 13 1970 325.6800 -20
## 14 1971 326.3200 -21
## 15 1972 327.4533 -22
## 16 1973 329.6767 -23
## 17 1974 330.1767 -24
## 18 1975 331.0842 -25
## 19 1976 332.0533 -26
## 20 1977 333.7817 -27
## 21 1978 335.4092 -28
## 22 1979 336.7817 -29
## 23 1980 338.6792 -30
## 24 1981 340.1000 -31
## 25 1982 341.4375 -32
## 26 1983 343.0250 -33
## 27 1984 344.5800 -34
## 28 1985 346.0408 -35
## 29 1986 347.3842 -36
## 30 1987 349.1608 -37
## 31 1988 351.5633 -38
## 32 1989 353.0675 -39
## 33 1990 354.3475 -40
## 34 1991 355.5667 -41
## 35 1992 356.3825 -42
## 36 1993 357.0675 -43
## 37 1994 358.8225 -44
## 38 1995 360.7958 -45
## 39 1996 362.5875 -46
## 40 1997 363.7050 -47
## 41 1998 366.6525 -48
## 42 1999 368.3258 -49
## 43 2000 369.5242 -50
## 44 2001 371.1300 -51
## 45 2002 373.2150 -52
## 46 2003 375.7742 -53
## 47 2004 377.4908 -54
## 48 2005 379.8000 -55
## 49 2006 381.9042 -56
## 50 2007 383.7642 -57
## 51 2008 385.5858 -58
## 52 2009 387.3667 -59
## 53 2010 389.8442 -60
## 54 2011 391.6258 -61
## 55 2012 393.8192 -62
## 56 2013 396.4817 -63
## 57 2014 398.5517 -64
## 58 2015 401.8629 -65