Penn StateNASA

Lab 6: Carbon Cycle Modeling


The goals of this lab are to:

  1. Simulate the fate of CO2 added to the atmosphere
  2. Interpret what happens when we increase the rate of fossil fuel burning
  3. Evaluate how the rate of carbon addition affects the maximum temperature and minimum pH attained under fossil fuel burning
  4. Simulate the impact of the different emission scenarios on global temperature
  5. Evaluate how permafrost melting amplifies warming under the different emission scenarios

Please make sure that you read the Introduction to the lab.  Skipping it will result in a lot of confusion and a lower score on the assignment.


In the lab activity for this module, we will be working with a STELLA model of the global carbon cycle that is attached to the climate model we used in Module 3. This model incorporates the processes of carbon transfer in the terrestrial and oceanic realms discussed in the previous sections; it also includes the history (from 1880 to 2010) of human impacts on the carbon cycle in the form of emissions from burning fossil fuels, burning forests, and disrupting the soil. The model is initially set up to represent the carbon cycle in a steady state just before the industrial revolution, at which point human alterations to the carbon cycle began in earnest.  We will use this model to explore different carbon emissions scenarios for the future, to see how the climate and carbon cycle might respond.

Here is what the model looks like, in a very simple, stripped-down form:

Schematic to show a simplified version of the Pre-Industrial Carbon Cycle
Pre-Industrial Carbon Cycle
Credit: D.Bice

In this figure, you can see the initial amounts of carbon in each reservoir in GT and the annual flows of carbon between reservoirs in GT C/yr. A couple of things should be pointed out about this model. In some cases, flows have been combined into things called bi-flows that have arrows on each end; this means that carbon can flow either way. There is a bi-flow connecting the atmosphere and surface ocean that represents the two-way transfer of ocean-atmosphere exchange. There is another bi-flow connecting the surface and deep ocean that represents upwelling and downwelling combined into one.

The real model is quite a bit more complex-looking, as can be seen below:

Schematic to show a complex, more complete, version of the Pre-Industrial Carbon Cycle
Land use changes represent the combined effects of soil distribution, deforestation, etc. for the last 100 years, extrapolated into the future 200 years at the year 2000 level.
Credit: D. Bice

The complications here arise from the fact that most of the flows are expressed by rather complicated equations. Most of the flows have small red arrows attached to them; these show you what things affect the flow rate. Think of the circle in the middle of the flow pipe as being a valve on a water pipe that controls how much water moves through the pipe. In many cases, the red arrows come from the reservoir that is being drained by the flow; this means that the outflow is dependent on how much is in the reservoir — when you have more in a reservoir, the outflow is often greater, and in essence, the outflow is a percentage of how much is in the reservoir.

Note that the atmosphere reservoir is connected to a converter called pCO2 atm — this is the concentration of CO2 in the atmosphere and the units are in parts per million or ppm, the same units that CO2 concentrations are typically given in. In this model, the initial amount of carbon in the atmosphere gives a pCO2 value of 280 ppm (and by now, it is just over 400 ppm). The pCO2 atm converter is in turn connected to the same climate model we used in Module 3, where it determines the strength of the greenhouse effect.  The climate model calculates the temperature at each moment in time and then passes that information back to the carbon cycle model in the form of a converter called global temp change, which is the change in global temperature relative to the starting temperature — this is like a temperature anomaly. 

The global temp change converter is then attached to a couple of other converters that attach to the photosynthesis flow and the soil respiration flow. Both of these flows are sensitive to temperature and the temperature combines with something called a temperature sensitivity. You can see something above called the Tsens sr — this is the temperature sensitivity for soil respiration. Both photosynthesis and soil respiration are sensitive to the temperature; they increase with temperature. Global temp change is also connected to the surface temperature of the oceans (T surf) via a "ghosted" version of the converter — a dashed line version that helps eliminate so many long connecting arrows running all over the diagram.

The photosynthesis flow has lots of converters associated with it since it is dependent on numerous factors, but the main things are temperature and the atmospheric CO2 concentration.

The model also includes a whole set of connected converters in the upper right that do all of the calculations related to the carbon chemistry in the oceans; this is where the pH or acidity of the surface ocean is calculated.

At the very top right of the model, there is a converter called Observed Atm CO2 that contains the observed history of atmospheric CO2 concentration since 1880. This is in the model so that we can test how good the model is. If our carbon cycle model is good, then it should calculate a pCO2 that closely matches the observed record.

The model includes the history of carbon emissions from burning fossils fuels shown in the previous section; it also includes a history of land use changes that impact the carbon cycle. These land use changes are broken up into tree burning (that accompanies deforestation) and soil disruption (related to farming); they are then additions to the flow of carbon from the land biota and soil back into the atmosphere. These human alterations to the carbon cycle are shown in the model by clicking on the pink graph icons labeled ffb and land use changes on the right side of the model, as shown below.

As before, the model we will work runs on a browser; here is a link to the model. When you follow the link you should see a screen like this:

Screenshot of interactive model of the carbon cycle
Screenshot: Interactive model of the carbon cycle.
Credit: D. Bice

You can think of this as the control panel for the model, where you can run it, stop it, make changes, and look at the results in the form of different graphs. You can access the different graphs by clicking on the triangular tab at the lower left of the graph window. The controls here consist of a set of switches you can turn on or off — these will determine which of 3 different emissions scenarios are applied to the model. The three buttons on the right-hand side show what the 3 emissions scenarios look like.  The video below explains how to work the switches.  There is also a switch that can be used to either include or exclude the carbon emissions to the atmosphere from human-related land use changes such as deforestation and soil disruption, but we will not mess around with this switch — just leave it in the "on" or up position as it is in the diagram above.

This initial model is set up to run from 1880 to 2100. Later, we will work with a model that runs farther into the future.

Video: Carbon Cycle Intro (4:15)