Rome Didn't Fall in A Day.









Objective Truth Exists and is Accessible to Everyone.

All Human Problems can be Solved with Enough Knowledge, Wealth, Social Cooperation and Time.


Photo: Rusty Peak, Anchorage, Alaska


Translate

Thursday, May 3, 2012

The Envelope, Probability Theory, and the Black Swan


The Envelope, Probability Theory, and the Black Swan
The Envelope
Test pilots have a term that describes the known operating characteristics of an airplane: “The Envelope”.   The Envelope is the range of known, safe, operating parameters, such as top speed, stall speed, angle of attack, maximum bank, etc.  When the test pilot tries to extend the range of operating parameters, he is operating “outside of the Envelope”, and the airplane may perform in unexpected or dangerous ways.  There are abrupt discontinuities in the physics of flight, and test pilots’ reputation for courage is well-deserved. 

The concept of “The Envelope” applies to our total experience, as individuals, as businesses, as nations, or as a species.     When we are within the envelope of our experience, events unfold more or less as we encountered them before.  When we are outside the envelope of our experience, our theories about probability break down.   We are in a place where established rules do not apply, and where our experience is irrelevant, or worse, misleading.  We are in the realm of the Black Swan.   We will return to the Black Swan, but first we need to step back and consider the history of probability.
---------
Probability Theory
Our lives are ruled by probability.  Will the stock market rise or fall?  Will tomorrow be sunny or rainy?  Will I get a raise at work?   Will the Dodgers win the pennant?  Will I have an automobile accident?  Will I win the lottery?  Will I be late to the airport?  Will she get pregnant?  Will l die of cancer?   From the mundane to the life-shattering, every day is filled with uncertainty about events, and throughout our lives, we develop ways to understand those probabilities and to anticipate the results.

A Priori Knowledge and Discrete Outcomes
Gamblers were the first to take an interest in probability theory, for obvious reasons.  In the 16th and 17th centuries, gamblers and mathematicians (often being one in the same) began to develop the theory of probability in the simplest cases.  They learned how to calculate the odds for any event with a set of known a priori parameters and discrete outcomes.  A priori (“from earlier”) means that we inspected the dice before the throw, or counted the cards before the deal.   Discrete outcomes are critical; when we throw dice we do not expect the outcome to be eight and three-quarters, or the king of spades.  If we consider all of the possible discrete outcomes, we can calculate precisely the odds of winning or losing (assuming that the dealer isn’t cheating). 

Non-A Priori Cases
But what about situations without a priori knowledge? 
In the eighteenth century, mathematicians considered the problem of a bag filled with an unknown number of white balls and black balls.  We are not given the opportunity to inspect the bag, or count the balls.   The probability of drawing either white or black can only be discovered through experience.  As our experience grows, we gain confidence about the ratio of black to white balls in the bag, although we can never be completely certain of the total probability until we have drawn the final ball from the bag.  Modern sampling theory can quantify the degree of confidence in the probability as a function of how many balls we have drawn (i.e. the extent of our experience).
--------
Bayes’ Theorem
Bayes’ Theorem, developed in the eighteenth century, applies to probability problems without  knowledge of a priori conditions.   Bayes’ theorem is used to combine a subjective estimate of probability with previous experience about the actual rate of occurrence.

Suppose a bird watcher claims to have seen a rare species of duck.  Previous experience shows that in nature, the common bird is observed 99% of the time; the rare bird 1% of the time.  And suppose our experience also shows that amateur bird watchers accurately identify birds 90% of the time, and mis-identify birds 10% of the time. 
Bayes’ theorem gives the probability that the bird watcher actually saw the rare species as only about 8%.   (Bayesian calculator:  http://psych.fullerton.edu/mbirnbaum/bayes/BayesCalc.htm)

Consider the case of the Ivory-Billed Woodpecker, which is generally regarded as extinct since the last confirmed sighting in 1944.   In 2004, scientists from Cornell University claimed a sighting of the bird, creating a swarm of media interest.   Here is a blurry 2-second video which documents their sighting:
Is this the long-missing Ivory-Billed Woodpecker, or the similar, but relatively common Pileated Woodpecker?  What does Bayes’ Theorem say?


The actual population of the Ivory-Billed (assuming it exists) must be very low, compared to the Pileated Woodpecker.  Let’s assume there may be a total population of 10 Ivory-Billed Woodpeckers, and about 10,000 Pileated Woodpeckers.   Further, let’s assume that our scientists, paddling a canoe through a swamp, can identify a flying bird correctly 85% of the time.  Then we can do the math:  The probability that the observed bird was actually the rare Ivory-Billed Woodpecker is only 5 out of 1000.

This result brings to mind a rule of scientific analysis from Cornell’s most famous scientist, Carl Sagan:  Extraordinary claims require extraordinary evidence.”   Bayes’ Theorem is the mathematical expression of that principle.

Bayesian Theory in Geological Estimates of Success for Petroleum Prospects
One application of Bayes theorem is to calculate the true probability of an event, using subjective guesses about an event, and the actual rate of occurrence.   The method is essentially a “force-fit” of subjective guesses to the actual rate of occurrence, using the historical accuracy of previous guesses.   

We can extend the Bayesian method to subjective estimates of success used by petroleum geologists in prospect appraisal.  

In oil exploration, geologists assign a probability of commercial success to each prospect.  Over time, the cumulative experience of success and failure provides a means to review the accuracy of the predictions, provided that the predictions were made using the same methodology.

The following chart shows 74 exploration prospects drilled between 1996 and 2000.  The probability of success is the vertical axis, and prospects are shown in rank order, according to geological chance of success.  Successes are color coded in red, and failures in blue.  The chart demonstrates that geologists' estimates have merit; successful prospects generally occur on the left side of the chart.  But let's take a closer look.   Are the estimates quantitatively correct?  Can the estimates be improved with the Bayesian technique?


2  The second chart shows the prospect portfolio, roughly in thirds according to risk.  The upper third (highest chance of success) slightly outperformed the estimates, with 57% actual success, compared to 48% estimated success.  The middle group moderately underperformed the estimates, with 17% actual success compared to 28% in estimates.  And the bottom third substantially underperformed the estimates, with only 7% success, compared to 18% predicted success.



We can look at the performance of the entire portfolio by summing the estimated probability of success, to create a curve showing the cumulative predicted number of discoveries across the portfolio (blue curve).  We can compare the actual cumulative discoveries (red curve), which rises by integer steps over the successful prospects.  Actual results closely parallel the predictions to the midpoint (about 30% chance of success).  Actual results then trail the predictions to the bottom third (about 20% chance of success), where the actual results go flat, showing no success corresponding to prospects estimated at less than 20% chance of success.

Using Excel, we can run a regression on the curve representing the cumulative actual discoveries, relating the relationship between the rate of predicted to actual discoveries.  This function adjusts the estimated probabilities to actual results, and provides a predictive means to forecast future probabilities.  











The same function can be applied in a predictive fashion to a new portfolio of prospects.   The second group shares characteristics of the first group of prospects.   Success is concentrated in the lower-risk part of the portfolio.  Actual success is greater than predicted in the low-risk part of the spectrum, and success is almost absent in the higher-risk part of the spectrum.  The function derived from the first prospect group is not a perfect fit, but improves the fit of pre-drill estimates to actual results.  The cumulative success for the program using the Baysian-adjusted probability of success  very close to the actual success of the program.

Tracking predictions and results, combined with Bayesian methods allow the calculation of true probabilites from subjective estimates.  The method can be used iteratively to adapt to changes in exploration technology or improved estimates by the geological staff.  The process provides a way to obtain quantitatively better risk-adjusted investment decisions, and to concentrate attention on prospects most likely to yield success.
-------
Estimation of Low-Probability Events
People are very good at estimating probabilities in the middle range.  We fairly accurately assess the flip of a coin, the outcome of a college football game, or even the chance of a full house in five-card poker (14.4%). But people are terrible at estimating chances on the ends of the probability spectrum.  We are usually unable to discern the difference between the probability of events at 1:100, and 1:1000, or even 1:10,000 chances.   The same phenomenon occurs at the other end of the spectrum, for events of very high likelihood, but less than certainty.  We are simply unable to sense or quantify the difference.   

An example of this problem is shown in the risk assessments for the Space Shuttle program.  Richard Feynman wrote a stinging critique of the NASA risk estimates following the Challenger disaster.
Feynman's report includes this quote:
(Engineers at Rocketdyne, the manufacturer, estimate the total probability 
[of catastrophic failure]as 1/10,000. Engineers at marshal estimate it as 1/300, while NASA management, to whom these engineers report, claims it is 1/100,000. An independent engineer consulting for NASA thought 1 or 2 per 100 a reasonable estimate.)
The actual rate of failure was 2 disasters out of 135 missions, or about 1/67.

How is it possible that such wildly differing estimates existed, regarding the safety of such an important project?

Part of the problem is sampling.  For rare events, we must make a large number of observations to detect and quantify a possibility.  If we walk across a lake on thin ice ten times and do not fall through the ice, we can conclude that the chance of falling through the ice is probably less than 1:10.  It does not mean that walking on thin ice is safe, or that we can safely cross the lake 100 times.   For rare events, we simply cannot gain enough experience to adequate grasp the true probability.  This is particularly problematic for risky and dangerous events.

There are many other facets to the problem, including self-interest of management, and various other sources of bias which Feynman discusses in his report.

But the simple lesson that I take away is that people simply cannot understand risk in the range of low probability events.
------------
Excuse me, I will finish this post soon!
Black Swan Events
Author Nicholas Taleb introduced the term Black Swan into the modern vocabulary of risk in his book "Fooled by Randomness", in 2004.       



Wednesday, April 18, 2012

Modeling Global CO2 Cycles

This article is the fifth post in a series about Global CO2 trends and seasonal cycles.

1)  The Keeling Curve
      http://dougrobbins.blogspot.com/2011/05/keeling-curve.html
2)  The Keeling Curve and Seasonal Carbon Cycles
3)   Seasonal Carbon Isotope Cycles
4)   Long-Term Trends in Atmospheric CO2
5)   Modeling Global CO2 Cycles

In this post, I generate a model for the global CO2 record from 1971 - 2009.  Inputs to the model include agricultural biomass, fossil fuel emissions, absorption of excess CO2 by carbon sinks, and atmospheric mixing between the Northern and Southern Hemispheres.
The final model begins in the year 1971, and yields a set of CO2 curves, by latitude, that closely matches the actual record.

The important thing is that this model is made entirely by data resulting from human influences.   Natural factors also exist and clearly influence global CO2.  But the quantitative influence of agriculture and fossil fuel use is more than enough to model the annual cycles and long-term rise of CO2 in the atmosphere.

The ease with which the model was created, and the lack of any reasonable, quantifiable alternatives, indicate that changes in atmospheric CO2 are primarily the result of human activity.

Modeling
Modeling performs an important function in science.   After gathering data and making observations about pertinent parameters of a problem, it is necessary to put the pieces together in a quantitative model, to see if the parameters are acting as we expect.

The model “talks back” to us, in a way.   The model will show a fit to real-world data when our assumptions are reasonable, and a mis-fit to the data when the model is built on incorrect assumptions.   In this case, the model shows that the seasonal cyclicity of the Keeling curve is not the result of seasonal fossil fuel use, but instead is the result of seasonal photosynthesis and oxidation.  There is a good quantitative and geographic fit between the observed carbon cycles and the volume of carbon in agricultural biomass.   The increase of amplitude in CO2 cycles in proportion to population growth also supports the idea that the observed annual cycles are largely the result of agriculture.

This model considers separately the CO2 flux in the Northern Hemisphere and the Southern Hemisphere, and matches observations for the rate of mixing between the hemispheres.  A more complex model could be built, perhaps at increments of 5 or 10 degrees of latitude, and more closely identify the location of agriculture and fossil-fuel emissions, and that might be useful to addressing deeper questions.  But I believe a model should be simple in essence; sufficiently complex to answer the question at hand, and not any more complex.  This model is intended to answer the question of human influences on global CO2, and the division of the globe into two hemispheres is sufficient to answer that question.

Fossil Fuel Use and Annual Cycles
Fossil-fuel use has a cyclicity with the appropriate seasonal peaks for the Northern Hemisphere.  I spent quite a bit of time finding data on global fossil-fuel usage, assuming that it was a major factor in the Keeling curve annual cycles.  I extrapolated the seasonal patterns of use for coal, natural gas, and oil in the United States to global figures from the International Energy Agency (IEA) for the three commodities.  (I made an adjustment for the summer "air-conditioning"  bump in coal and natural gas, when extrapolating to the entire world.)
I converted the volumes (in gigatonnes) of CO2 emitted from fossil fuels to atmospheric CO2 concentrations in parts per million for the northern hemisphere and southern hemisphere.


In my last post,
http://dougrobbins.blogspot.com/2012/03/long-term-trends-in-atmospheric-co2.html, we saw that fossil fuel emissions are responsible for the long-term increase of atmospheric CO2.  But the peak of fossil fuel usage is in the months of December and January, but the steepest gains of CO2 occur in September through November.   In the fall months in high northern latitudes, CO2 concentrations rise by about 15 ppm.  The sum of fossil fuel emissions in those months total only about 1.5 ppm, an order of magnitude lower than the rise seen in the seasonal data.  Factors other than fossil fuels are causing most of the seasonal fluctuation of atmospheric CO2.

The Model

We can build a simple model to help identify the significant parameters driving global CO2 cycles, and quantify parameters wherever possible.   From our earlier observations, we can construct the model using the following parameters:
  • CO2 taken up by Plants during the growing season
  • Oxidation  of carbon in plants following the growing season
  • CO2 emissions from Fossil Fuels
  • Absorption of CO2 by carbon sinks (e.g. oceans)
  • Exchange of CO2 between Northern and Southern Hemispheres.
 We already observed that the behavior of CO2 cycles differs greatly by latitude.  The Northern Hemisphere, with 68% of the world’s landmass, 88% of the world’s population, and 83% of the worlds GDP, has a cycle showing very large seasonal fluctuations in CO2.  The Southern Hemisphere shows much less annual CO2 fluctuation.  The sharpest and largest fluctuations in CO2 occur in the summer and fall months of the Northern Hemisphere.  In the summer, plants are taking up carbon through photosynthesis, and atmospheric CO2 declines.  Immediately following the growing season, CO2 concentration rebounds sharply, as plants give CO2 back to the atmosphere through oxidation.



I constructed a model for annual CO2 uptake through photosynthesis, beginning with the volume of biomass generated through agriculture.  Agriculture generates about 140 gigatonnes of biomass every year.
http://www.unep.or.jp/Ietc/Publications/spc/WasteAgriculturalBiomassEST_Compendium.pdf
Adjustments for moisture content (50%), carbon content (45%), and conversion to CO2 (3.67x) results in about 96 gigatonnes of CO2 removed from the Northern Hemisphere atmosphere annually.   Keep in mind that this is only half of the air on the planet.  Thus, during the growing season, CO2 in the Northern Hemisphere falls sharply.

In the model, I distributed agricultural carbon and fossil fuel use according to economic output by hemisphere.  The Northern Hemisphere represents 83% of global economic output, and the Southern Hemisphere represents 17% of global economic output.

I assigned the 96 gigatonnes of agricultural CO2 intake in the summer growing months, as shown in the following graph.  For the oxidation part of the cycle, we can observe a very sharp rebound in CO2 in the data during the fall months.  It is possible that some of the rebound is from CO2 sinks, seeking equilibrium after the change during the growing season.  However, isotope data shows an equally sharp rebound. (http://dougrobbins.blogspot.com/2012/03/seasonal-carbon-isotope-cycles.html)   It appears to me that vegetation is giving back to the atmosphere the very same CO2 that was absorbed during the summer.
I adopted an oxidation/respiration model to return the CO2 to the atmosphere as a zero-sum annual exchange.  I tried an exponential decline for the oxidation part of the cycle, then tweaked it to match the annual CO2 cycles of the Northern Hemisphere high latitudes (with the long-term trend removed).

This model produced a surprisingly easy fit to the high latitude data of the Northern Hemisphere (see below).
Note that the long-term rising CO2 trend has been removed from the real world data, and there are no CO2 emissions from fossil fuels in the model at this point.

The total volume of vegetation includes both natural as well as agricultural biomass.  I found a single estimate for the net global annual uptake by plants, of about 60 gigatonne of carbon, or about 220 gigatonnes of CO2.   According to these estimates, agriculture represents 52% of the total annual carbon uptake by plants.  Of the total annual carbon cycle, some carbon is exchanged with carbon sinks (soil) and some carbon exchanged with the atmosphere.  The fit of the model to observed data shows that the NET amount of carbon taken out of the atmosphere by plants, and returned by oxidation, is very close to the volume of carbon taken up by agricultural activity.

I followed a similar procedure to model the Southern Hemisphere.  I found that 17% of global agricultural biomass produced CO2 fluctuations that were far too large to match the data in the Southern Hemisphere. I found a good match by using only 5% of global agricultural biomass.   The chart below shows the model parameters.
The following chart shows the match of the model to the data from the Southern Hemisphere.

Annual cycles from intermediate latitudes have lower amplitude than cycles from the high northern latitudes.  This was the topic of an earlier post: http://dougrobbins.blogspot.com/2012/03/keeling-curve-and-seasonal-carbon.html.  The Northern Hemisphere, with its large CO2 fluctuations, dominates global CO2 cycles.  CO2 cycles from low latitudes in the Southern Hemisphere (pink) follow the seasonal pattern of the Northern Hemisphere, showing the range and influence of atmospheric mixing between Northern and Southern Hemispheres.
I tried a simple mixing model to represent the cycles observed in intermediate latitudes.   The chart below shows a 50%-50% mixture at the equator, and 70% - 30% mixtures at intermediate latitudes.   As shown in the data above, the cycles of intermediate latitude (pink line) in the Southern Hemisphere follow the seasonal pattern of the Northern Hemisphere.
A more sophisticated model could be created, using greater detail in the location of agriculture by latitude, but I think this model demonstrates that atmospheric mixing between the Northern and Southern Hemispheres can reasonably explain the range of amplitude in CO2 cycles in intermediate latitudes.






Global CO2 data show distinctive characteristics of annual cyclicity and a long-term rising trend ("the Keeling Curve").   Subtle aspects of the curve include a rising rate of increase, and an increase in the amplitude of the cycles.

The final model runs from the year 1971 to 2009.  As a starting point, the model used values for the average CO2 concentration of the Northern and Southern Hemispheres in 1971, of 327 and 325 parts per million CO2, respectively.

The photosynthetic model, which was developed for the year 2009, was adjusted for earlier years as a function of global population. This resulted in cycles with increasing amplitude through the range of the model.
Agricultural production was assumed to vary directly as a function of population, but incremental agriculture was assumed to displace natural vegetation.  Growth of CO2 intake through photosynthesis was increased at a rate of 50% of incremental agricultural output (back-calculated from the 2009 model).

Carbon dioxide from fossil fuel emissions was added, according to estimates from IEA and the BP statistical review of world energy.  Annual figures given in these reports were scheduled on a monthly basis, by analogy to US monthly consumption of coal, natural gas, and oil, as described above.  As noted in a previous post, about 40% of fossil fuel CO2 emissions are absorbed by carbon sinks, including the ocean.  This fraction of new carbon emissions was removed from the model on a monthly basis.

The Northern Hemisphere receives the bulk of fossil fuel CO2 emissions, and modeled CO2 rises rapidly in Northern Hemisphere, unless a transfer to the Southern Hemisphere is allowed.  In an earlier post, we saw that rising CO2 in the Southern Hemisphere lags CO2 in the Northern Hemisphere, by a period of about 22 months.  In the model, I transfered half of the excess CO2 of the Northern Hemisphere to the Southern Hemisphere, using a lag of 22 months to represent the necessary mixing time.
Despite the general simplicity of the model, the resulting CO2 curve shows a reasonable correlation to actual data recorded across the global range of latitudes, and after 38 years of CO2 addition and subtraction, the model concludes at the appropriate concentrations of CO2 across the globe.

Conclusions:
1)  A model can be generated which provides a very good match to the long-term global CO2 record.  The model includes estimated fossil fuel use, absorption of CO2 by carbon sinks, carbon accumulation in agricultural biomass, and oxidation of agricultural biomass.  The volume of agricultural biomass was varied in the model according to world population growth.
2)  Surprisingly, fossil fuel use does not have a significant effect on seasonal CO2 cycles.   Known volumes and timing of fossil-fuel emissions do not match the cyclicity in CO2 observations.
3)  Photosynthesis in the Northern Hemisphere, dominates the seasonal cycles.  The volume of CO2 absorbed through agriculture closely matches the net volume of CO2 taken up by both natural and agricultural photosynthesis.
 4)  Oxidation of vegetation occurs quickly.  Three quarters of the net plant biomass is oxidized in the first three months following the growing season.  It seems likely to me that burning of agricultural waste accounts for some of the rapid oxidation following the growing season.
5)  The Northern Hemisphere dominates both seasonal and long-term trends in atmospheric CO2.
6)  Mixing between the hemispheres accounts very well for the gradation of cyclicity observed at intermediate latitudes.

The model shows that the long-term trend of rising CO2 is attributable to fossil-fuel emissions.  Fossil fuel emissions account quantitatively for the rise in CO2 over the last 38 years, and fit the data with regard to differences in concentration in the Northern and Southern Hemispheres.
The model also shows that the annual cyclicity of the biologic cycle is strongly influenced by agriculture.  Agricultural biomass alone can be used to model  and match observed data for seasonal CO2 cyclicity.

And finally, the uptake of CO2 through agriculture clearly outpaces emissions of CO2 from fossil fuels, at least on a seasonal basis.  As a tool for the management of CO2 concentrations, policy-makers should consider banning the burning of agricultural waste, and consider options for disposal of agricultural waste as a means of sequestering significant volumes of carbon.
--------------
Global CO2 concentration data in this report is credited to C. Keeling and others at the Scripps Institute of Oceanography, also Gaudry et al, Ciattaglia et al, Columbo and Santaguida, and Manning et al.  The data can be found on the Carbon Dioxide Information Analysis Center.
 http://cdiac.ornl.gov/trends/co2/
Data for CO2 released by fossil fuels is available from EIA CO2 Emissions from Fuel Consumption,
http://www.iea.org/co2highlights/co2highlights.pdf
And the BP Statistical Review of World Energy:
http://www.bp.com/sectionbodycopy.do?categoryId=7500&contentId=7068481

Monthly data for US fossil fuel consumption were taken from the EIA website:

Global population figures from 1970 - 2010 were taken from Wikipedia.
The estimate for annual global biomass, circa 2009 was taken from a UN report:
http://www.unep.or.jp/Ietc/Publications/spc/WasteAgriculturalBiomassEST_Compendium.pdf

Thursday, March 29, 2012

Long-term Trends in Atmospheric CO2


This article is the fourth post in a series about global CO2 trends and seasonal cycles.

1)  The Keeling Curve
2)  The Keeling Curve and Seasonal Carbon Cycles
3)   Seasonal Carbon Isotope Cycles
4)   Long-Term Trends in Atmospheric CO2
5)   Modeling Global CO2 Cycles

Long-Term Trends in Atmospheric CO2
Burning fossil fuels consumes oxygen and emits CO2.  This is simple chemistry that we learned in middle school, using a candle and a glass jar.  It is possible to calculate the quantity of CO2 released from the candle by the chemistry of the candle, and calculate the concentration of CO2 in the jar from the quantity of CO2 released by the candle.
We can also easily find the quantity of CO2 released into the atmosphere by the global consumption of coal, oil and natural gas.  Some respected organizations (the International Energy Agency, and the British Petroleum Statistical Review of World Energy) have collected the data on fossil fuels, and done the math on the quantity of CO2 released.
We can add the annual fossil fuel CO2 emissions to a baseline CO2 concentration, such as the 1970 global average, and compare the cumulative fossil fuel emissions to the observed change in world average CO2 concentration.
Only about 60% of the known CO2 emissions from burning fossil fuels winds up in the atmosphere.   We can infer that the other 40% is absorbed by various earth systems which act as carbon reservoirs or sinks.   Examples might include vegetation, the ocean, or precipitation of limestone.

In an earlier post, we explored the annual cycles of CO2 fluctuation, and how those cycles vary by latitude.
http://dougrobbins.blogspot.com/2012/03/keeling-curve-and-seasonal-carbon.html
We can also take the annual average for CO2 concentration in the Northern Hemisphere, and compare those readings to the annual average CO2 concentration in the Southern Hemisphere.  We can see a clear difference.  The Southern Hemisphere lags behind the Northern Hemisphere in terms of increasing CO2 by about 2.6 ppm.


Fossil fuel emissions are concentrated in the Northern Hemisphere; as we noted in the previous post, 90% of the world's population lives in the Northern Hemisphere, including the most industrialized economies.
It seems reasonable to conclude that the Northern Hemisphere leads the Southern Hemisphere in CO2 concentration gains, because the Northern Hemisphere is the source of the emissions.  The difference in CO2 concentration is very close to the quantity of global annual CO2 emissions.





The excess CO2 concentration of the Northern Hemisphere is closely matched, and easily explained by annual fossil fuel emissions of the Northern Hemisphere.  I allocated global CO2 emissions by hemisphere according to figures for national GDP.   In earlier posts, we have seen the strong correlation between GDP and energy use.  About 83% of global GDP, and by inference, CO2 emissions, occur in the Northern Hemisphere.  The excess CO2 delivered to the atmosphere in the Northern Hemisphere accounts very well for the difference in CO2 between the Northern and Southern Hemispheres.















The use of fossil fuels has grown exponentially since the industrial revolution ( M.K. Hubbert, 1956; D.H. Meadows, et al 2004).   Annual consumption has grown from essentially zero before the year 1800, to over 10.4 billion tonnes of oil equivalent in 2010.   

Air bubbles in ice cores from the Antarctic ice cap provide a record of historic and prehistoric concentrations of CO2 in the atmosphere.   The pre-industrial concentration of CO2 is generally reported around 280 ppm.   The record shows that CO2 levels were essentially constant from the year 1000 until the industrial revolution began, about 1800.

The increase in atmospheric CO2 fits well to an exponential function, reflecting the exponential growth of human population, and the exponential growth in the use of fossil fuels.  Trial and error produced this function, which ties pre-industrial CO2 concentration, and provides a good fit to the modern Keeling Curve.
CO2 concentration, ppm = e(n*0.001854) + 280,
where n = the number of months since Jan. 1800



Extrapolating the function forward, CO2 concentration would be expected to reach 450 ppm in the year 2031, and 500 ppm in the year 2042.

A carbon dioxide concentration of 450 ppm is sometimes cited as a theoretical "tipping point", beyond which climate change becomes irreversible, due to positive feedback mechanisms (i.e., release of greenhouse gasses from permafrost, dissolution of carbonate sediments due to ocean acidification, release of methane from gas hydrates, etc.).

Based on the consistency of the exponential increase in atmospheric CO2, the growth of world population, and industrialization of the world economy, it seems likely that atmospheric CO2 levels will continue to rise and likely exceed the 450 ppm and 500 ppm thresholds within the next 35 years.
----------------
This is fourth of five posts in a series about global atmospheric CO2.
The Keeling Curve
http://dougrobbins.blogspot.com/2011/05/keeling-curve.html

The Keeling Curve and Seasonal Carbon Cycles
Seasonal Carbon Isotope Cycles
Long-Term Trends in Atmospheric CO2
Modeling Global CO2 Cycles

-----------------
Global CO2 data is available from Keeling et. al., on the Carbon Dioxide Information Analysis Center website.  http://cdiac.ornl.gov/trends/co2/

The population chart was prepared by "Radical Cartographer" Bill Rankin.  http://www.radicalcartography.net/
I used the version of the map posted here: http://www.geekosystem.com/world-population-latitude-longitude/

Data for CO2 released by fossil fuels is available from EIA CO2 Emissions from Fuel Consumption,
http://www.iea.org/co2highlights/co2highlights.pdf

And the BP Statistical Review of World Energy:
http://www.bp.com/sectionbodycopy.do?categoryId=7500&contentId=7068481

Historic CO2 levels:


Exponential Growth in Fossil Fuels



BP Statistical Review of World Energy, Global Fossil Fuel Consumption: 

Thursday, March 22, 2012

Seasonal Carbon Isotope Cycles

The carbon isotope composition of CO2 in the atmosphere fluctuates in annual cycles, much like CO2 in the atmosphere itself.  Carbon 12 is the most common isotope, representing about 99% of atmospheric carbon. Carbon 13 represents most of the remaining percent.  Carbon 14 is an extremely small component (about one part per trillion), mostly generated by nuclear tests during the 1960s.  The amount of C14 in the air has been rapidly declining since the elimination of above-ground nuclear testing. 

The ratio of C13 to C12 is expressed as a standard measure: dC13/C12  (usually pronounced "del-thirteen").  The measure represents the ratio of C13 to C12, as compared to a standard ratio.  The equation for dC13/C12 is:

dC13/C12 =  ((C13/C12 ratio of sample / C13/C12 ratio of standard) -1) * 1000. 
The equation simply expresses the difference between the sample and the standard, expressed in tenths of percent.

The light isotope, C12, is more easily taken up in plants during photosynthesis.  Thus, plants, and plant-derived carbon (including you and me, and other things which eat plants) have negative dC13/C12 ratios.  Coal, which is derived from wood, and oil, which is derived from algae,  also have negative dC13/C12 ratios (in the range of -20 to -35).  Natural gas, which can be formed by several mechanisms, may have a very negative dC13/C12 ratio (from -20 to -50).

So when dC13/C12 is rising in the atmosphere, as happens in the Northern Hemisphere summer,  it is because plants are taking up the light C12, and the ratio of C13 remaining in the atmosphere is increasing.  When dC13/C12 is falling, as happens in the fall through spring, it is because animals and bacteria in the biosphere are respirating, giving back the C12 taken up by plants.  In addition, coal, oil and gas are being burned, giving back C12 taken up by plants long ago, and causing the dC13/C12 ratio to fall.

Let's look at the data.  
As seen in previous posts, atmospheric monitoring stations have collected data on CO2 concentrations and isotope rations across a wide range of latitudes.
Isotope trends show a similar seasonal cyclicity and amplitude dependence on latitude as seen in global CO2 concentrations.
This is a chart of C13/C12 observations in carbon dioxide, at latitudes ranging from the arctic circle to the south pole.  There is an annual, asymmetric cyclicity to the measurements, and a gradual downward trend, indicating progressively lighter isotopic composition in the atmosphere.

Let's take a more detailed look at the cycles.   This chart shows the C13/C12 readings from 2003 to 2006.  Observation stations are color-coded by latitude, with warm colors indicating the southern hemisphere, and cool colors showing the northern hemisphere.  The asymmetrical Northern Hemisphere cycles are exactly the inverse of what we saw in the previous CO2 charts.
Let's take a close look at the seasonal nature of the C13/C12 cycles.   The isotope signature rises sharply in the Northern Hemisphere summer, and falls through the remainder of the year, sharply at first, and then more gradually.  Southern hemisphere observations show a very weak opposite polarity to the northern hemisphere.  If we take the southern hemisphere readings as a baseline for global C13/C12, it is clear that there are forces producing both positive and negative seasonal deflections from that baseline.


















We can remove the long-term trend from the chart, to see the cyclicity better.   I subtracted a one-year moving average of measurements at each station from the data, to produce the chart of relative fluctuation.  Differences between northern and southern hemispheres are apparent, as is the asymmetry of the cycles.
Here is a closer look at the cycles with the long-term trend removed.  Observations from the northern hemisphere have very high amplitude cycles.  Southern hemisphere cycles are relatively flat, exactly as we saw previously in the previous CO2 composition charts.

We can make a plot of peak-to-trough amplitude by latitude, as we did for the CO2 cycles.  Annual cycles grow larger at higher latitudes in the Northern Hemisphere, falling slightly near the North Pole, exactly as seen in the CO2 chart.


And the amplitude of C13/C12 cycles by latitude can be compared to the distribution of population, as we did with cycles of atmospheric CO2.
Conclusions:
  • Northern Hemisphere C13/C12 observations show a seasonal cyclicity, rising in the summer and falling in the winter.  The cycles are consistent with light C12 being absorbed by plants during the growing season, and with light C12 being released to the atmosphere by plant oxidation and more fossil-fuel use during the winter.  
  • The cycles show both positive and negative deflections in the northern hemisphere, relative to the southern hemisphere baseline.  The absorption of C12 by plants in the growing season is the strongest and sharpest part of the cycle.  
  • The long term trend is toward more negative C13/C12, consistent with the accumulation of atmospheric CO2 from fossil fuel use.

I will leave a few loose ends and questions to address, but will place this on the blog today.   I will try to tie up some of the loose ends as soon as I can.
Loose Ends:  
>  The observed amplitudes of the C13/C12 isotope cycles should be compared to the estimated isotope changes produced by fossil fuel emissions.   These can be modeled, using known annual consumption volumes  and C13/C12 ratios.  Are cycles produced by annual fossil-fuel consumption of the same size as the observed cycles?
>  The magnitude of the long-term C13/C12 trend should be compared to estimates from fossil fuel use, to establish a reasonable origin for increasingly negative isotope signature in the atmosphere over time.
Questions:
>  As with the CO2 cycles, the sharpest annual changes in the C13/C12 cycle occur in the summer months, when light isotopes are being taken up by plant growth.  We can easily see that agriculture is concentrated in the Northern Hemisphere.  What is the estimated impact of agriculture on CO2 and isotope cycles?
>  What is the reason for the sharply negative rebound in C13/C12, following the summer growing season? >  We should be able to explain the full shape of the CO2 and isotope annual cycles.  The steepest part of the fall occurs before the main heating season for fossil fuel consumption.  The answer might be a rapid oxidation of plant matter by bacterial action, burning of agricultural waste, or some interaction of CO2 sources and sinks.  Demonstration of the actual mechanism, supported by data, would be helpful.
> The amplitude of C13/C12 cycles, like the amplitude of CO2 cycles, increases northward to a point beyond the Arctic Circle, before slightly diminishing at the highest latitudes observed.   It is not clear why this is so.  Some factors to consider are: 1)  the volume of air available for dispersion decreases per degree of latitude  northward, simply due to the curvature of the earth.  2)  Mixing with southern latitude air diminishes at higher northern latitudes.  3)  Winter fossil fuel usage per capita probably increases at highest northern latitudes.  For example, the gas utility serving Anchorage, Alaska, has about a 12-fold swing in fuel delivery in winter, as compared to summer months.
------
This article is the third of a series of articles about global CO2.  The final article consolidates and summarizes the previous posts.
2)  The Keeling Curve and Seasonal Carbon Cycles


3)   Seasonal Carbon Isotope Cycles


4)   Long-Term Trends in Atmospheric CO2
5)   Modeling Global CO2 Cycles

6)   The Keeling Curve Summary:  Seasonal CO2 cycles and Global CO2 Distribution
       http://dougrobbins.blogspot.com/2013/05/the-keeling-curve-seasonal-co2-cycles.html
-------
References:
References:
Atmospheric CO2 Carbon Isotope Data:
Keeling, R.F. S.C. Piper, A.F. Bollenbacher, and S.J. Walker. 2010. Monthly atmospheric 13C/12C isotopic ratios for 11 SIO stations. In Trends: A Compendium of Data on Global Change. Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, U.S. Department of Energy, Oak Ridge, Tenn., U.S.A.

Global Emissions average isotope data
Boden, T.A., G. Marland, and R.J. Andres. 2013. Global, Regional, and National Fossil-Fuel CO2 Emissions. Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, U.S. Department of Energy, Oak Ridge, Tenn., U.S.A. doi 10.3334/CDIAC/00001_V2013




All CO2 data in this article is credited to C. Keeling and other at the Scripps Institute of Oceanography, also Gaudry et al, Ciattaglia et al, Columbo and Santaguida, and Manning et al.  The data can be found on the Carbon Dioxide Information Analysis Center; http://cdiac.ornl.gov/trends/co2/

World Background Map for charts courtesy ESRI.

The population chart was prepared by "Radical Cartographer" Bill Rankin.  http://www.radicalcartography.net/
I used the version of the map posted here: http://www.geekosystem.com/world-population-latitude-longitude/