Model specifications for place-specific carbon sequestration and storage models developed for the Puget Sound, San Pedro, and Rocky Mountains case studies
Definitions
Model structure and assumptions
Case studies 1. Western Washington State 2. San Pedro River, Arizona and Sonora 3. Madagascar 4. Rocky Mountains
Carbon sequestration source models. Although carbon sequestration data are available globally at 1 km2 resolution, we developed simple Bayesian network models that include the influences on carbon sequestration (e.g., vegetation, soils, climate). Existing datasets can be used in ARIES to provide mean values to use in training finer-grained models, allowing estimation of carbon sequestration changes in scenarios or for up-scaled modeling of carbon sequestration where higher resolution input data are available. Based on the literature (e.g., Iverson et al. 1994, Eade and Moran 1996, Gaston et al. 1998, Chan et al. 2006, Naidoo and Ricketts 2006, Egoh et al. 2008, Wundscher et al. 2008, Auch 2010, Wendland et al. 2010, Tallis et al. 2013) and discussions with regional experts, we set carbon sequestration as a function of vegetation density and sequestration rate, two intermediate variables created to keep conditional probability tables tractable (Marcot et al. 2006). We set sequestration rate as a function of soil C:N ratio and the difference between mean summer high and winter low (in Madagascar and Western Washington), and as a function of land cover, vegetation type, and actual evapotranspiration (in Orange County). We set vegetation density as a function of hardwood:softwood ratio, percent tree canopy cover, and successional stage (in Western Washington), and percent tree canopy cover and forest degradation status (in Madagascar). For the San Pedro, Orange County, Rocky Mountain, and Vermont agricultural carbon models, we used a collapsed number of variables, removing the intermediate nodes for vegetation density and sequestration rate. For the San Pedro model, we estimate sequestration as a function of vegetation type, percent tree canopy cover, and mean annual precipitation. We added a node for bark beetle kill in the Rocky Mountain model, where outbreaks of pine and spruce beetles have killed large numbers of trees, lowering sequestration in forests with a substantial number of dead trees (Coops and Wulder 2010). For the Orange County model, we used the above noted variables as input nodes to sequestration rate, then combined sequestration rate with percent tree canopy cover to estimate annual vegetation and soil carbon sequestration. Actual evapotranspiration (AET) has been found to have a strong relationship with primary productivity, and therefore carbon sequestration (Lieth and Box 1972, Elegene et al. 1989, Metherell et al. 1993). This is especially true in water-limited regions such as semi-arid biomes, as with the Orange County case study (Claudio et al. 2006, Fuentes et al., 2006). Vegetation type can help to predict the quantity of vegetation sequestration and storage capacities from expected biomass for certain plant species (Kirby and Potvin 2007). In the Vermont model, we estimated sequestration as a function of vegetation carbon storage (itself a function of mean annual precipitation, vegetation type, and the difference between mean summer high and winter low) and soil C:N ratio (Liu et al. 2010). We used Jenks natural breaks to discretize summer high-winter low, soil C:N ratio, and actual evapotranspiration. We used equal intervals to discretize vegetation and soil carbon sequestration, hardwood:softwood ratio, and percent tree canopy cover. We based prior probabilities for the models on either the actual distribution of regional data (where we have these datasets), expert opinion (where consensus by experts was possible), or uninformed priors (where there was true uncertainty and a lack of consensus by experts). We filled out conditional probability tables by setting extremes set at both ends (i.e., “pegging the corners,” Marcot et al. 2006) and interpolating intermediate values. Where possible we used expert opinion about which variables are most influential, and which should have the greatest influence on the contingent probability tables, and what the general level of uncertainty was for that system (i.e., how wide to set the distribution of values across discrete states). All else being equal, we set vegetation density at its highest values at greater percent tree canopy cover, later successional stages, more softwoods, and no forest degradation (where applicable). We set sequestration rate with its highest values at higher C:N ratios, higher actual evapotranspiration, lower differences between mean summer high and winter low temperatures, and land cover and vegetation types with greater biomass (where applicable). We set sequestration to its greatest values at high levels of vegetation density and sequestration rate.
Potential stored carbon release sink models. Stored carbon release can be calculated probabilistically or deterministically, though deterministic calculations are likely more accurate assuming historical or forecasted land-use and fire data are available. Deterministic estimates of carbon loss due to land-use change can be calculated by modeling carbon storage under pre- and post-change conditions. Carbon loss due to fire can be estimated, for example, by overlaying fire boundary polgyons (e.g., GeoMAC 2013) with fuel consumption coefficients (Spracklen et al. 2009) and carbon pool data (Smith et al. 2006), or by applying more advanced models (e.g., Lutes 2013). This method was demonstrated for the U.S. Pacific Northwest by Bagstad et al. (in press). We have also estimated stored carbon release probabilistically, as a function of vegetation and soil carbon storage (the sum of vegetation carbon storage and soil carbon storage) and the risk of deforestation and/or fire, with greater stored carbon release at higher risk and carbon storage levels. Soil carbon storage is influenced by slope, soil pH, soil oxygen conditions (i.e., greater storage in wetlands where anaerobic conditions inhibit respiration), vegetation density (an intermediate variable incorporating tree canopy cover and degradation status in Madagascar, tree canopy cover, and vegetation type in the San Pedro, and successional stage, tree canopy cover, and hardwood:softwood ratio in Western Washington, noted as important determinants of carbon sequestration in the Pacific Northwest by Nelson et al. 2008), and soil carbon:nitrogen ratio. A simpler model was applied in the Rocky Mountains, incorporating only annual precipitation (Derner and Schuman 2007) and soil order (Buringh 1984). The importance of these variables in influencing soil carbon dynamics has been noted by previous authors. We set vegetation carbon storage as a function of the difference between mean summer high and winter low temperature (Auch 2010) and vegetation density, with population density added as an influence in Madagascar. For the San Pedro, we set vegetation carbon storage as a function of mean annual precipitation and vegetation density, and for the Rocky Mountains we set it as a function of vegetation type and tree canopy cover (Smith et al. 2006). For the Orange County model, deforestation was not considered as an influence on stored carbon release (though it would be included in non-urban areas within the same biome), slope was dropped as an influence on soil carbon storage (since slope/aspect influence AET and other water balance measurements in chaparral and scrub ecosystems, Miller 1947, Parsons 1973, Ng and Miller 1980), and actual evapotranspiration and percent tree canopy cover were added as influences on soil carbon storage. We set vegetation carbon storage as a function of land cover, vegetation type, percent tree canopy cover, and AET for the Orange County model. The Vermont model used soil tillage and biomass removal rate as influences on agricultural stored carbon release (Gollany et al. 2010, Gonzalez-Chavez et al. 2010). This model considered soil C:N ratio, biomass residue input (Hai et al. 2010), and vegetation type as influences on soil carbon storage and vegetation type, mean annual precipitation, and the difference between mean summer high and winter low temperature. Iverson et al. (1994) and Gaston et al. (1998) provide discretization of continuous variables for slope and population density. Bosworth and Tricou (1999) and Darby et al. (2009) provide discretization for vegetation carbon storage in the Vermont carbon model. We used Jenks natural breaks to discretize soil carbon storage, summer high-winter low, vegetation and soil carbon storage, soil C:N ratio, vegetation carbon storage, fire frequency, and actual evapotranspiration. We used equal intervals to discretize hardwood:softwood ratio and percent tree canopy cover. All else being equal, we set soil carbon storage at its highest values at low or high pH, high C:N ratio, level slopes, greater vegetation density and annual precipitation, and on anoxic (i.e., wetland) soils, and vice versa. We set vegetation carbon storage at its greatest values with low differences between mean summer high and winter low temperature, high vegetation density or tree canopy cover, and low population density (in Madagascar). We set stored carbon release at its highest with greater vegetation and soil carbon storage and greater deforestation and fire risk. The outputs of the carbon sink model are either the potential stored carbon release (probabilistic calculation of stored carbon release) or vegetation and soil carbon storage (deterministic calculation of stored carbon release).
Greenhouse gas emissions use models. The beneficiaries of carbon sequestration and storage are greenhouse gas emitters who release CO2 into the atmosphere. Spatially explicit data on greenhouse gas emissions exist for the United States. Globally, we use population density data multiplied by per capita emissions for the country or sub-national region of interest.
Carbon flow models. Since carbon dioxide is relatively quickly mixed in the atmosphere, the benefits of carbon sequestration and storage can be enjoyed by any human beneficiary on Earth, regardless of location. As such, no flow model is necessary for carbon sequestration and storage. However, for a given region, we can calculate the differential between carbon uptake by ecosystems (sequestration minus release of stored carbon) and anthropogenic carbon release. This information can be used in a flow model to show whether that region has a negative or positive carbon balance, i.e., whether its emissions are greater or less than the amount of carbon sequestered. Key outputs from the flow models include: 1. Carbon mitigation surplus: Calculated when local sequestration exceeds emissions plus atmospheric carbon sources. 2. Carbon mitigation deficit: Calculated when local emissions exceed net carbon uptake (sequestration minus stored carbon release).
Spatial data sources
Models | Data theme | Source | Spatial extent | Spatial resolution | Date |
Carbon sequestration | All models | NBII-Millennium Ecosystem Assessment | Global | 1 km | 2000 |
Carbon sequestration & Stored carbon release – Western WA | Forest successional stage | BLM/Interagency Vegetation Mapping Project | Western Washington & Oregon | 25 m | 1996 |
Hardwood: softwood ratio | BLM/Interagency Vegetation Mapping Project | Western Washington & Oregon | 25 m | 1996 | |
Carbon sequestration & Stored carbon release – San Pedro, Rocky Mountains | Mean annual precipitation | PRISM/Oregon State | United States | 800 m | 1971-2000 |
Carbon sequestration & Stored carbon release – San Pedro, Western WA, Rocky Mountains | Percent tree canopy cover | NLCD 2001 | United States | 30 m | 2001 |
Carbon sequestration & Stored carbon release – Madagascar | GLCF/Univ. of Maryland | Global | 1 km | 2000 | |
Carbon sequestration & Stored carbon release – Madagascar, Western WA | Soil C:N ratio | FAO soils | Global | 0.0833 min | 1970-1978 |
Carbon sequestration & Stored carbon release – Western WA | Summer high – winter low | PRISM/Oregon State | United States | 800 m | 1971-2000 |
Carbon sequestration & Stored carbon release – Madagascar | WorldClim | Global | 30 arc-seconds | 1950-2000 | |
Carbon sequestration & Stored carbon release – San Pedro, Rocky Mountains | Southwest Regional GAP Analysis (SWReGAP) | AZ, CO, NM, NV, UT | 30 m | 1999-2001 | |
Carbon sequestration & Stored carbon release – Rocky Mountains | Northwest GAP Analysis (NWGAP) | CA, ID, MT, OR, WA, WY | 30 m | 1999-2001 | |
Stored carbon release – Rocky Mountains | Bark beetle kill | USDA Forest Service-Region 2 | CO, KS, SD, WY | Vector polygon data | 2010-present |
Stored carbon release – Madagascar | Deforestation risk | GLCF/Univ. of Maryland | Global (processed only for Madagascar) | 250 m | 2001-2005 |
Stored carbon release – San Pedro | Southwest Regional Gap Analysis LULC & TNC fire data | AZ, CO, NM, NV, UT | 30 m | 2000 | |
Stored carbon release – Western WA | Washington DNR & Oregon Dept. of Forestry | Washington & Oregon | 1.5 derived from point data | 1970-2007 | |
Stored carbon release – Madagascar | Population density | LANDSCAN/Oak Ridge National Lab | Global | 30 arc-second | 2006 |
Stored carbon release – Madagascar, San Pedro, Western WA | Slope | Derived from global SRTM data | Global | 90 m | n/a |
Stored carbon release – Global models | Soil carbon storage | FAO soils | Global | 0.0833 min2 | 1970-1978 |
Stored carbon release – U.S. models | Soil carbon storage | SSURGO soils data | United States | 30 m | n/a |
Stored carbon release – San Pedro, Western WA | Soil oxygen conditions (e.g., wetlands) | NLCD 2001 | United States | 30 m | 2001 |
Stored carbon release – Madagascar | Kew Gardens Madagascar vegetation map | Madagascar | 30 m | 1999-2003 | |
Stored carbon release – Rocky Mountains | Soil order | STATSGO soils data | United States | Vector | n/a |
Stored carbon release – San Pedro, Western WA | Soil pH | SSURGO soils data | United States | 30 m | n/a |
Stored carbon release – Madagascar | FAO soils | Global | 0.0833 min | 1970-1978 | |
Stored carbon release – All U.S. case studies | Vegetation carbon storage | National Biomass and Carbon Dataset | United States | 30 m | 2000 |
Stored carbon release – Madagascar | CDIAC/Ruesch & Gibbs | Global | 1 km | 2000 | |
Use – All U.S. models | GHG emissions | VULCAN Project, Purdue Univ. | United States | 10 km | 2002 |
Use – Madagascar | Population density | LANDSCAN, Oak Ridge National Lab | Global | 30 arc-second | 2006 |
Per capita emissions | Energy Information Administration: International Energy Annual | Global | Aspatial, by country | 2006 |
References
Acknowledgements and additional contributors
Ted Auch and Serguei Krivov provided input on the initial ARIES carbon models. Mark Casias developed the case study for Orange County. Sam Gorton developed the agricultural carbon case study for Vermont. Dave Batker, Jim Pittman, and Paula Swedeen provided data and model review for the Western Washington case study. Miro Honzak provided data and model review for the Madagascar case study. An expert review panel including individuals from the U.S. Geological Survey, University of Arizona, Bureau of Land Management, and other organizations provided data and model review for the San Pedro case study. Initial ARIES data and models for the Rocky Mountains were developed by students participating in a graduate level ecosystem services modeling course taught in the University of Denver’s Department of Geography in the fall of 2011; James Reed, Darius Semmens, and Todd Hawbaker assisted with further data and model refinement.