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_hist (t_ini10, j, land) |
Land area for historial time steps | \(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 | |
|---|---|
| ac | Age classes |
| ag_pools(c_pools) | Above ground carbon pools |
| bii_class_secd(bii_class44) | bii coefficent land cover classes secondary vegetation |
| bii_class44 | bii coefficent land cover classes |
| c_pools | Carbon pools |
| cell(i, j) | number of LPJ cells per region i |
| clcl | climate classification types |
| clcl58 | simple climate classes |
| consv_type | Type of land conservation |
| 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) |
| kcr(kve) | Cropping activities |
| land | Land pools |
| land58 | Peatland categories |
| landcover44 | land cover classes used in bii calculation |
| luh2_side_layers10 | side layers from LUH2 |
| manPeat58(land58) | Drained and managed peatland categories |
| marginal_land29 | Marginal land scenarios |
| 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 |
| relocation_target29 | Cropland requiring relocation based on different SNV targets |
| stockType | Carbon stock types |
| 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 |
| w | Water supply type |
Jan Philipp Dietrich, Florian Humpenoeder, Kristine Karstens
11_costs, 13_tc, 14_yields, 22_land_conservation, 29_cropland, 30_croparea, 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