The peatland module calculates GHG emissions from degrading/drained peatlands.
| Description | Unit | A | B | |
|---|---|---|---|---|
| pcm_land (j, land) |
Land area in previous time step | \(10^6 ha\) | x | |
| pm_climate_class (j, clcl) |
Koeppen-Geiger climate classification on the simulation cluster level | \(1\) | x | |
| pm_interest (t_all, i) |
Interest rate in each region and timestep | \(\%/yr\) | x | |
| vm_land (j, land) |
Land area of the different land types | \(10^6 ha\) | x |
| Description | Unit | |
|---|---|---|
| vm_peatland_cost (j) |
One-time and recurring cost of managed peatland | \(10^6 USD_{05MER}/yr\) |
| vm_peatland_emis (j) |
GHG emissions from managed peatland | \(t CO2eq/year\) |
In this realization, peatlands do not exist. Therefore, GHG emissions from degrading peatlands are assumed zero.
Limitations Peatland area and associated GHG emissions are fixed to zero.
In this realization, the state of peatlands is modelled as described in Humpenöder et al. (2020). The initial map of peatland area consists of intact and degraded peatland area for the year 2015. Future peatland dynamics depend on the ratio of total peatland area and total land area. GHG emissions from degraded and rewetted peatlands as are calculated based on GHG emission factors from the “2013 Supplement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories: Wetlands”.
Land transition matrix for peatland area
\[\begin{multline*} \sum_{from58,to58} v58\_lu\_transitions(j2,from58,to58) = \sum_{man58,land58} pc58\_peatland\_man(j2,man58,land58) + pc58\_peatland\_intact(j2) \end{multline*}\]
\[\begin{multline*} \sum_{from58} v58\_lu\_transitions(j2,from58,to58) = v58\_peatland\_man(j2,"degrad","crop")\$sameas(to58,"degrad\_crop") + v58\_peatland\_man(j2,"degrad","past")\$sameas(to58,"degrad\_past") + v58\_peatland\_man(j2,"degrad","forestry")\$sameas(to58,"degrad\_forestry") + v58\_peatland\_man(j2,"unused","crop")\$sameas(to58,"unused\_crop") + v58\_peatland\_man(j2,"unused","past")\$sameas(to58,"unused\_past") + v58\_peatland\_man(j2,"unused","forestry")\$sameas(to58,"unused\_forestry") + v58\_peatland\_man(j2,"rewet","crop")\$sameas(to58,"rewet\_crop") + v58\_peatland\_man(j2,"rewet","past")\$sameas(to58,"rewet\_past") + v58\_peatland\_man(j2,"rewet","forestry")\$sameas(to58,"rewet\_forestry") + v58\_peatland\_intact(j2)\$sameas(to58,"intact") \end{multline*}\]
\[\begin{multline*} \sum_{to58} v58\_lu\_transitions(j2,from58,to58) = pc58\_peatland\_man(j2,"degrad","crop")\$sameas(from58,"degrad\_crop") + pc58\_peatland\_man(j2,"degrad","past")\$sameas(from58,"degrad\_past") + pc58\_peatland\_man(j2,"degrad","forestry")\$sameas(from58,"degrad\_forestry") + pc58\_peatland\_man(j2,"unused","crop")\$sameas(from58,"unused\_crop") + pc58\_peatland\_man(j2,"unused","past")\$sameas(from58,"unused\_past") + pc58\_peatland\_man(j2,"unused","forestry")\$sameas(from58,"unused\_forestry") + pc58\_peatland\_man(j2,"rewet","crop")\$sameas(from58,"rewet\_crop") + pc58\_peatland\_man(j2,"rewet","past")\$sameas(from58,"rewet\_past") + pc58\_peatland\_man(j2,"rewet","forestry")\$sameas(from58,"rewet\_forestry") + pc58\_peatland\_intact(j2)\$sameas(from58,"intact") \end{multline*}\]
The following two equations calculate land expansion and land contraction based on the above land transition matrix.
\[\begin{multline*} v58\_expansion(j2,to58) = \sum_{from58\$\left(not sameas(from58,to58)\right)}\left( v58\_lu\_transitions(j2,from58,to58)\right) \end{multline*}\]
\[\begin{multline*} v58\_reduction(j2,from58) = \sum_{to58\$\left(not sameas(from58,to58)\right)}\left( v58\_lu\_transitions(j2,from58,to58)\right) \end{multline*}\]
Future peatland degradation (v58_peatland_man) depends on changes of managed land, scaled with the ratio of total peatland area and total land area (p58_scaling_factor). By multiplying changes in managed land with this scaling factor we implicitly assume that intact peatlands are distributed equally within a grid cell. The following example illustrates the mechanism used for projecting peatland dynamics: In a given grid cell, the total land area is 50 Mha and the total peatland area is 10 Mha. Therefore, the scaling factor is 0.2 (10 Mha divided by 50 Mha). If cropland expands by 5 Mha, 1 Mha of intact peatland is converted to degraded peatland (5 Mha0.2). If the total cell would become cropland, degraded peatland would equal to the total peatland area (50 Mha 0.2 = 10 Mha).
\[\begin{multline*} v58\_peatland\_man(j2,"degrad",land58) = pc58\_peatland\_man(j2,"degrad",land58) + \left(\left(vm\_land(j2,land58)-pcm\_land(j2,land58)\right) \cdot p58\_scaling\_factor(j2)\right)\$\left(\sum_{ct} m\_year(ct)>s58\_fix\_peatland\right) \end{multline*}\]
This constraint avoids the conversion of intact peatland into rewetted peatland.
\[\begin{multline*} \sum_{stat\_rewet58} v58\_expansion(j2,stat\_rewet58) \leq \sum_{stat\_degrad58}\left( v58\_reduction(j2,stat\_degrad58) + v58\_expansion(j2,stat\_degrad58)\right) - v58\_reduction(j2,"intact") \end{multline*}\]
Costs for peatland degradation and rewetting
\[\begin{multline*} vm\_peatland\_cost(j2) = v58\_peatland\_cost\_annuity(j2) + \sum_{land58} v58\_peatland\_man(j2,"rewet",land58) \cdot \sum_{ct} i58\_cost\_rewet\_recur(ct) + \sum_{degrad58,land58} v58\_peatland\_man(j2,degrad58,land58) \cdot \sum_{ct} i58\_cost\_degrad\_recur(ct) \end{multline*}\]
\[\begin{multline*} v58\_peatland\_cost\_annuity(j2) = \left(\sum_{stat\_rewet58} v58\_expansion(j2,stat\_rewet58) \cdot \sum_{ct} i58\_cost\_rewet\_onetime(ct) + \left(v58\_reduction(j2,"intact") + \sum_{stat\_rewet58} v58\_reduction(j2,stat\_rewet58)\right) \cdot \sum_{ct} i58\_cost\_degrad\_onetime(ct)\right) \cdot \sum_{cell(i2,j2),ct}\left(\frac{pm\_interest(ct,i2)}{\left(1+pm\_interest(ct,i2)\right)}\right) \end{multline*}\]
GHG emissions from managed peatlands (degraded and rewetted)
\[\begin{multline*} v58\_peatland\_emis(j2,emis58) = \sum_{man58,land58,clcl58}\left( v58\_peatland\_man(j2,man58,land58) \cdot p58\_mapping\_cell\_climate(j2,clcl58) \cdot p58\_ipcc\_wetland\_ef(clcl58,land58,emis58,man58)\right) \end{multline*}\]
\[\begin{multline*} vm\_peatland\_emis(j2) = \sum_{emis58} v58\_peatland\_emis(j2,emis58) \end{multline*}\]
Limitations Peatland area and GHG emissions are fixed to 2015 levels for the historic period, depending on
s58_fix_peatland. Organic carbon stocks in peatlands are not accounted for.
| Description | Unit | A | B | |
|---|---|---|---|---|
| f58_ipcc_wetland_ef (clcl58, land58, emis58, ef58) |
Wetland GWP100 emission factors | \(t CO2eq/ha\) | x | |
| f58_peatland_degrad (j) |
Degrading peatland area | \(10^6 ha\) | x | |
| f58_peatland_intact (j) |
Intact peatland area | \(10^6 ha\) | x | |
| i58_cost_degrad_onetime (t) |
One-time costs for peatland degradation | \(USD_{05MER}/ha\) | x | |
| i58_cost_degrad_recur (t) |
Recurring costs for degraded peatland | \(USD_{05MER}/ha\) | x | |
| i58_cost_rewet_onetime (t) |
One-time costs for peatland restoration | \(USD_{05MER}/ha\) | x | |
| i58_cost_rewet_recur (t) |
Recurring costs for rewetted peatland | \(USD_{05MER}/ha\) | x | |
| p58_intact_ratio (t, j) |
Ratio of intact and total peatland | \(1\) | x | |
| p58_ipcc_wetland_ef (clcl58, land58, emis58, man58) |
Wetland GWP100 emission factors | \(t CO2eq/ha\) | x | |
| p58_land_area (j) |
Total land area | \(10^6 ha\) | x | |
| p58_man_land_area (j) |
Total managed land | \(10^6 ha\) | x | |
| p58_mapping_cell_climate (j, clcl58) |
Mapping between cells and climate regions | \(binary\) | x | |
| p58_peatland_area (j) |
Total peatland area | \(10^6 ha\) | x | |
| p58_peatland_degrad_unused (j) |
Intermediate calculation in peatland initialization | \(10^6 ha\) | x | |
| p58_peatland_degrad_unused_weight (j, land58) |
Weight for intermediate calculation in peatland initialization | \(1\) | x | |
| p58_peatland_degrad_used (j) |
Intermediate calculation in peatland initialization | \(10^6 ha\) | x | |
| p58_peatland_degrad_weight (j, land58) |
Weight for peatland distribution to land58 | \(1\) | x | |
| p58_scaling_factor (j) |
Scaling factor for managed peatland | \(1\) | x | |
| pc58_man_land_shr (j, land58) |
Share of total managed land | \(1\) | x | |
| pc58_peatland_intact (j) |
Intact peatland | \(10^6 ha\) | x | |
| pc58_peatland_man (j, man58, land58) |
Managed peatland | \(10^6 ha\) | x | |
| pc58_peatland_man_inital (j, man58, land58) |
Managed peatland at initialization | \(10^6 ha\) | x | |
| q58_expansion (j, to58) |
Peatland expansion | \(10^6 ha\) | x | |
| q58_peatland_cost (j) |
One-time and recurring cost of peatland conversion and management | \(10^6 USD_{05MER}/yr\) | x | |
| q58_peatland_cost_annuity (j) |
Annuity costs of peatland conversion in the current timestep | \(10^6 USD_{05MER}/yr\) | x | |
| q58_peatland_degrad (j, land58) |
Constraint for peatland degradation | \(10^6 ha\) | x | |
| q58_peatland_emis (j) |
GHG emissions from managed peatland | \(t CO2eq/year\) | x | |
| q58_peatland_emis_detail (j, emis58) |
Detailed GHG emissions from managed peatland | \(t CO2eq/year\) | x | |
| q58_peatland_rewet (j) |
Constraint for peatland rewetting | \(10^6 ha\) | x | |
| q58_reduction (j, from58) |
Peatland reduction | \(10^6 ha\) | x | |
| q58_transition_from (j, from58) |
Peatland transitions from | \(10^6 ha\) | x | |
| q58_transition_matrix (j) |
Peatland transitions | \(10^6 ha\) | x | |
| q58_transition_to (j, to58) |
Peatland transitions to | \(10^6 ha\) | x | |
| s58_cost_degrad_onetime | One-time costs for peatland degradation | \(USD_{05MER}/ha\) | x | |
| s58_cost_degrad_recur | Recurring costs for degraded peatland | \(USD_{05MER}/ha\) | x | |
| s58_cost_rewet_onetime | One-time costs for peatland restoration | \(USD_{05MER}/ha\) | x | |
| s58_cost_rewet_recur | Recurring costs for rewetted peatland | \(USD_{05MER}/ha\) | x | |
| s58_fix_peatland | Year indicating until when peatland area should be fixed to 2015 levels | \(year\) | x | |
| s58_rewetting_switch | Peatland rewetting on (Inf) or off | \(0\) | x | |
| v58_expansion (j, stat58) |
Peatland expansion | \(10^6 ha\) | x | |
| v58_lu_transitions (j, from58, to58) |
Peatland transitions | \(10^6 ha\) | x | |
| v58_peatland_cost_annuity (j) |
Annuity costs of managed peatland expansion in the current timestep | \(10^6 USD_{05MER}/yr\) | x | |
| v58_peatland_emis (j, emis58) |
Detailed GHG emissions from managed peatland | \(t CO2eq/year\) | x | |
| v58_peatland_intact (j) |
Intact peatland | \(10^6 ha\) | x | |
| v58_peatland_man (j, man58, land58) |
Managed peatland | \(10^6 ha\) | x | |
| v58_reduction (j, stat58) |
Peatland reduction | \(10^6 ha\) | x |
| description | |
|---|---|
| cell(i, j) | number of LPJ cells per region i |
| clcl | climate classification types |
| clcl_mapping(clcl, clcl58) | Mapping between detailed and simple climate classes |
| clcl58 | simple climate classes |
| ct(t) | Current time period |
| degrad58(man58) | State of degraded peatland |
| ef58(man58) | Peatland emissions factors |
| emis58 | Wetland emission types |
| i | all economic regions |
| i2(i) | World regions (dynamic set) |
| j | number of LPJ cells |
| j2(j) | Spatial Clusters (dynamic set) |
| land | Land pools |
| land58(land) | Managed land types |
| man58 | State of managed peatland |
| stat_degrad58(stat58) | Peatland status degrad |
| stat_man58(stat58) | Peatland status managed land |
| stat_rewet58(stat58) | Peatland status rewet |
| stat58 | Peatland status |
| t_all(t_ext) | 5-year time periods |
| t(t_all) | Simulated time periods |
| type | GAMS variable attribute used for the output |
Florian Humpenöder
10_land, 11_costs, 12_interest_rate, 45_climate, 56_ghg_policy
Humpenöder, Florian, Kristine Karstens, Hermann Lotze-Campen, Jens Leifeld, Lorenzo Menichetti, Alexandra Barthelmes, and Alexander Popp. 2020. “Peatland Protection and Restoration Are Key for Climate Change Mitigation.” Environmental Research Letters 15 (10): 104093. https://doi.org/10.1088/1748-9326/abae2a.