The land module coordinates and analyzes all land related activities by summing up all land types and calculating the gross changes in land use between two time steps of optimization given the recursive dynamic structure of MAgPIE model. The land module tracks land use transitions by directly counting sources and targets of conversions.
| Description | Unit | A | |
|---|---|---|---|
| vm_landdiff_forestry | Aggregated difference in forestry land compared to previous timestep | \(10^6 ha\) | x |
| vm_landdiff_natveg | Aggregated difference in natveg land compared to previous timestep | \(10^6 ha\) | x |
| Description | Unit | |
|---|---|---|
| fm_land_iso (t_ini10, iso, land) |
Land area for different land pools at ISO level | \(10^6 ha\) |
| fm_luh2_side_layers (j, luh2_side_layers10) |
luh2 side layers | \(grid cell share\) |
| pcm_land (j, land) |
Land area in previous time step including possible changes after optimization | \(10^6 ha\) |
| pm_land_start (j, land) |
Land initialization area | \(10^6 ha\) |
| vm_cost_land_transition (j) |
Costs for lu transitions | \(10^6 USD_{05MER}/yr\) |
| vm_land (j, land) |
Land area of the different land types | \(10^6 ha\) |
| vm_landdiff | Aggregated difference in land between current and previous time step | \(10^6 ha\) |
| vm_landexpansion (j, land) |
Land expansion | \(10^6 ha\) |
| vm_landreduction (j, land) |
Land reduction | \(10^6 ha\) |
| vm_lu_transitions (j, land_from, land_to) |
Land transitions between time steps | \(10^6 ha\) |
The landmatrix_dec18 realization tracks land use transitions by directly counting sources and targets of conversions.
The following three equations describe the general structure of the land transition matrix. The first equation defines the total amount of land to be constant over time. The two balancing variables v10_balance_positive and v10_balance_negative are needed to avoid technical infeasibilities due to small differences in accuracy between variables and parameters in GAMS. The use of v10_balance_positive and v10_balance_negative is minimized by putting a high cost factor on these variables (q10_cost). In practice, v10_balance_positive and v10_balance_negativeshould deviate from zero only in exceptional cases.
\[\begin{multline*} \sum_{land\_from,land\_to} vm\_lu\_transitions(j2,land\_from,land\_to) + v10\_balance\_positive(j2) - v10\_balance\_negative(j2) = \sum_{land} pcm\_land(j2,land) \end{multline*}\]
\[\begin{multline*} \sum_{land\_from} vm\_lu\_transitions(j2,land\_from,land\_to) = vm\_land(j2,land\_to) \end{multline*}\]
\[\begin{multline*} \sum_{land\_to} vm\_lu\_transitions(j2,land\_from,land\_to) = pcm\_land(j2,land\_from) \end{multline*}\]
The following two equations calculate land expansion and land contraction based on the above land transition matrix.
\[\begin{multline*} vm\_landexpansion(j2,land\_to) = \sum_{land\_from\$\left(not sameas(land\_from,land\_to)\right)}\left( vm\_lu\_transitions(j2,land\_from,land\_to)\right) \end{multline*}\]
\[\begin{multline*} vm\_landreduction(j2,land\_from) = \sum_{land\_to\$\left(not sameas(land\_from,land\_to)\right)}\left( vm\_lu\_transitions(j2,land\_from,land\_to)\right) \end{multline*}\]
Small costs of 1 $ per ha on gross land-use change avoid unrealistic patterns in the land transition matrix
\[\begin{multline*} vm\_cost\_land\_transition(j2) = \sum_{land}\left( vm\_landexpansion(j2,land) + vm\_landreduction(j2,land)\right) \cdot 1 + \left(v10\_balance\_positive(j2) + v10\_balance\_negative(j2)\right) \cdot s10\_cost\_balance \end{multline*}\]
The gross changes in land are calculated based on land expansion, land contraction and land changes from within the modules 35_natveg and 32_forestry:
\[\begin{multline*} vm\_landdiff = \sum_{j2,land}\left( vm\_landexpansion(j2,land) + vm\_landreduction(j2,land)\right) + vm\_landdiff\_natveg + vm\_landdiff\_forestry \end{multline*}\]
Some of the land use transitions are restricted: No planted forest on natveg areas
vm_lu_transitions.fx(j,"primforest","forestry") = 0;Conversions within natveg are not allowed
vm_lu_transitions.fx(j,"primforest","other") = 0;
vm_lu_transitions.fx(j,"secdforest","other") = 0;Primforest can only decrease
vm_lu_transitions.fx(j,land_from,"primforest") = 0;
vm_lu_transitions.up(j,"primforest","primforest") = Inf;Limitations This realization only accounts for net land use transitions.
| Description | Unit | A | |
|---|---|---|---|
| f10_land (t_ini10, j, land) |
Different land type areas | \(10^6 ha\) | x |
| q10_cost (j) |
Costs for lu transitions | \(10^6 USD_{05MER}/yr\) | x |
| q10_landdiff | Land difference constraint | \(10^6 ha\) | x |
| q10_landexpansion (j, land_to) |
Land expansion constraint | \(10^6 ha\) | x |
| q10_landreduction (j, land_from) |
Land reduction constraint | \(10^6 ha\) | x |
| q10_transition_from (j, land_from) |
Land transition constraint from | \(10^6 ha\) | x |
| q10_transition_matrix (j) |
Land transition constraint cell area | \(10^6 ha\) | x |
| q10_transition_to (j, land_to) |
Land transition constraint to | \(10^6 ha\) | x |
| s10_cost_balance | Artificial cost for balance variable | \(USD_{05MER}/ha\) | x |
| v10_balance_negative (j) |
Balance variable for land transitions | \(10^6 ha\) | x |
| v10_balance_positive (j) |
Balance variable for land transitions | \(10^6 ha\) | x |
| description | |
|---|---|
| cell(i, j) | number of LPJ cells per region i |
| clcl | climate classification types |
| clcl58 | simple climate classes |
| cost58 | annunity cost categories |
| emis_source | Emission sources |
| emis58 | Wetland emission types |
| i | all economic regions |
| iso | list of iso countries |
| j | number of LPJ cells |
| j2(j) | Spatial Clusters (dynamic set) |
| land | Land pools |
| land58 | Peatland categories |
| luh2_side_layers10 | side layers from LUH2 |
| manPeat58(land58) | Drained and managed peatland categories |
| poll58(pollutants) | Wetland emissions that can be taxed |
| pollutants(pollutants_all) | subset of pollutants_all that can be taxed |
| potnatveg(luh2_side_layers10) | potentially forested biomes |
| t_all(t_ext) | 5-year time periods |
| t_ini10 | Time periods with land initialization data |
| t(t_all) | Simulated time periods |
| type | GAMS variable attribute used for the output |
| type32 | plantation type |
Jan Philipp Dietrich, Florian Humpenoeder, Kristine Karstens
11_costs, 13_tc, 14_yields, 22_land_conservation, 30_crop, 31_past, 32_forestry, 34_urban, 35_natveg, 39_landconversion, 44_biodiversity, 50_nr_soil_budget, 56_ghg_policy, 58_peatland, 59_som, 71_disagg_lvst, 80_optimization