GCAM v6 Documentation: Earth System Module – Hector v2.5.0

Documentation for GCAM
The Global Change Analysis Model

View the Project on GitHub JGCRI/gcam-doc


Earth System Module – Hector v2.5.0

This section describes the carbon-cycle climate module - Hector - that is available for use in GCAM. MAGICC5.3 (Wigley, 2008) has traditionally been the only climate module available in GCAM. Hector v2.5.0 is the default climate model (Hartin et al., 2015) within GCAM. Users still have the option of running MAGICC5.3 in GCAM5.1, but we will not be supporting this option going forward.

Hector, an open-source, object-oriented, reduced-form global climate carbon-cycle model, is written in C++. This model runs essentially instantaneously while still representing the most critical global-scale earth system processes. Hector has a three-part main carbon cycle: a one-pool atmosphere, three-pool land, and 4-pool ocean. The model’s terrestrial carbon cycle includes primary production and respiration fluxes, accommodating arbitrary geographic divisions into, e.g., ecological biomes or political units. Hector actively solves the inorganic carbon system in the surface ocean, directly calculating air– sea fluxes of carbon and ocean pH. Hector reproduces the global historical trends of atmospheric [CO2], radiative forcing, and surface temperatures. The model simulates all four Representative Concentration Pathways (RCPs) with equivalent rates of change of key variables over time compared to current observations, MAGICC, and models from CMIP5 (Hartin et al., 2015). Hector’s flexibility, open-source nature, and modular design facilitates a broad range of research in various areas.

There most notable change between Hector v1.1 and Hector v2.5.0 is the inclusion of a one-dimensional ocean heat diffusion model - DOECLIM (Kriegler, 2005; Tanaka and Kriegler, 2007). With this addition, Hector v2.5.0 exhibits improved vertical ocean heat uptake, as well as surface response to radiative forcing. (https://github.com/JGCRI/hector/pull/206) https://github.com/JGCRI/hector/releases

Hector Carbon Cycle diagram
Figure 1: Representation of Hector’s carbon cycle, land, atmosphere, and ocean. The atmosphere consists of one well-mixed box. The ocean consists of four boxes, with advection and water mass exchange simulating thermohaline circulation. At steady state, the high-latitude surface ocean takes up carbon from the atmosphere, while the low-latitude surface ocean off-gases carbon to the atmosphere. The land consists of a user-defined number of biomes or regions for vegetation, detritus and soil. At steady state the vegetation takes up carbon from the atmosphere while the detritus and soil release carbon back into the atmosphere. The earth pool is continually debited with each time step to act as a mass balance check on the carbon system.

GCAM-Hector interactions

Currently the GCAM sectors interact with Hector via emissions. At every time step, emissions from GCAM are passed to Hector. Hector converts these emissions to concentrations when necessary, and calculates the associated radiative forcing, as well as the response of the climate system and earth system (e.g., temperature, carbon-fluxes, etc.). Hector’s climate information can be used as a climate constraint for in a GCAM policy run.

Table 1: Emissions and sources from each sector passed to Hector.

Emission Sector Notes
CO2* AgLU, Energy  
CH4 AgLU, Energy, Industrial Processes  
N2O AgLU, Energy  
NH3 AgLU, Energy  
SO2 AgLU, Energy, Industrial Processes  
CO AgLU, Energy, Industrial Processes  
BC AgLU, Energy  
OC AgLU, Energy  
NOx AgLU, Energy, Industrial Processes  
NMVOC Energy, Industrial Processes  
C2F6 Energy, Industrial Processes  
CF4 Industrial Processes, Urban Processes  
SF6 Energy, Industrial Processes  
HFC134a Energy  
HFC32 Energy  
HFC125 Urban Processes  
HFC227ea Urban Processes  
HFC23 Urban Processes  
HFC236fa Urban Processes not included in Hector
HFC134a Industrial Processes  
HFC245fa Industrial Processes  
HFC365mfc Industrial Processes not included in Hector

* CO2 emissions from the AgLU sector are separate from CO2 emissions from the Energy sector. Any change in atmospheric carbon, occurs as a function of anthropogenic fossil fuel and industrial emissions (FA), land-use change emissions (FLC), and the atmospheric-ocean (FO) and atmosphere-land (FL) carbon fluxes.

dCatm/dt = FA(t) + FLC(t) - FO(t) - FL(t)

Land carbon pools change as a result of NPP, RH and land-use change fluxes, whose effects are partitioned among the carbon pools (Hartin et al., 2015).

Hector Outputs

At every time step Hector calculates and outputs key climate variables.

Atmosphere
  • Global mean temperature change
  • Total radiative forcing & radiative forcing of individual emissions
  • Atmospheric CO2 concentrations
Land
  • Air-land carbon fluxes
  • NPP - net primary production
  • RH - heterotrophic respiration
  • Carbon pools (vegetation, detritus, soil)
Ocean
  • Air-sea carbon fluxes
  • Carbon pools (high and low latitude surface, intermediate and deep)
  • Carbonate system (DIC, pCO2 , CO32-, pH, aragonite and calcite saturations)
  • Surface ocean temperature
  • Oceanic heat flux

Getting and Installing Hector for Use with GCAM

For users who are running GCAM with the Mac or Windows Release Package, Hector support is already compiled in. For users compiling from source or interested in getting the Hector source, please see the Hector section in How to Set Up and Build GCAM.

Policy options

Hector is a flexible, simple climate model. Users can run Hector under various configurations, parameters, and constraints.

Hector Configurations

By default, Hector’s carbon cycle model treats the entire land surface as a single, homogeneous ecosystem. However, it is possible to introduce some land surface heterogeneity by splitting the land surface into several different biomes with distinct parameters. This is explained in detail here.

Hector Parameters

Hector has many parameters can be user adjusted; see the online Hector manual for more information.

Constraints

The Hector model can be run subject to constraints that force the model to have a certain behavior. Technically, this means that the model’s components output user-provided data as opposed to their own calculations, similar to the data mode of a CESM sub-model. Currently, the available constraints include CO~2~ concentrations, total radiative forcing, and temperature.

IAMC Reference Card

Climate indicators

References

  1. Hartin, C. A., Patel, P., Schwarber, A., Link, R. P., and Bond-Lamberty, B. P.: A simple object-oriented and open-source model for scientific and policy analyses of the global climate system – Hector v1.0, Geosci. Model Dev., 8, 939-955, doi:10.5194/gmd-8-939-2015, 2015. link
  2. Hartin, C. A., Bond-Lamberty, B., Patel, P., and Mundra, A.: Ocean acidification over the next three centuries using a simple global climate carbon-cycle model: projections and sensitivities, Biogeosciences, 13, 4329-4342, doi:10.5194/bg-13-4329-2016, 2016. link
  3. Wigley, T. M. (2008), MAGICC/SENGEN 5.3: User manual (version 2), edited, p. 80, NCAR, Boulder CO.
  4. Online Hector manual