CAPRI Online Manual (update)

User Tools

Site Tools

Trace:
scenario_simulation

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
scenario_simulation [2022/01/17 09:50] – [Annex: Land supply and land transitions in the supply part of CAPRI] matszscenario_simulation [2023/09/08 12:07] (current) – [Land use, land use change and forestry (LULUCF)] massfeller
Line 357: Line 357:
 A number of recent developments are not covered in the previous exposition of supply model equations A number of recent developments are not covered in the previous exposition of supply model equations
  
-  -A series of projects have added a distinction of rainfed and irrigated varieties of most crop activities which is the core of the so-called “CAPRI-water” version of the system((A more complete presentation is given in [[https://ec.europa.eu/jrc/en/publication/eur-scientific-and-technical-research-reports/extension-capri-model-irrigation-sub-module]].)). +   - A series of projects have added a distinction of rainfed and irrigated varieties of most crop activities which is the core of the so-called “CAPRI-water” version of the system 
-  -Several projects have added endogenous GHG mitigation options((These are most completely included in the “trunk” version of the CAPRI system. For details, see, for example, [[http://publications.jrc.ec.europa.eu/repository/bitstream/JRC101396/jrc101396_ecampa2_final_report.pdf]].))  +  - Several projects have added endogenous GHG mitigation options 
-  -Several new equations serve to explicitly represent environmental constraints deriving from the Nitrates Directive and the NEC directive((These are most completely included in the “trunk” version of the CAPRI system but developments are still ongoing.)).  +  - Several new equations serve to explicitly represent environmental constraints deriving from the Nitrates Directive and the NEC directive
-  -A complete area balance monitoring the land use changes according to the six UNFCCC land use types (cropland, grassland, forest land, wetland, settlements, residual land) has been introduced for carbon accounting+
  
 ====Calibration of the regional programming models====   ====Calibration of the regional programming models====  
Line 395: Line 394:
  
  
-====Annex: Land supply and land transitions in the supply part of CAPRI====+==== LULUCF in the supply model of CAPRI ====
  
-**Introduction**+=== Introduction ===
  
-This technical paper explains how the most aggregate level of the CAPRI area allocation in the context of the supply models has been re-specified in the TRUSTEE project((https://www.trustee-project.eu/ +This technical paper explains how the most aggregate level of the CAPRI area allocation in the context of the supply models has been re-specified in the TRUSTEE ((https://www.trustee-project.eu/ 
-)) and subsequently adopted in the CAPRI trunk. The former specification for land supply and transformation functions focused on agricultural land use and the transformation of agricultural land between arable land and grass land((See https://svn1.agp.uni-bonn.de/svn/capri/trunk/doc/landSupplyCAPRI_v5.pdf (Torbjörn Jansson, Wolfgang Britz, Alan Renwick and Peter Verburg (2010) Modelling CAP reform and land abandonment in the European Union with a focus on Germany.) +)) and SUPREMA ((https://www.suprema-project.eu/)) projects and subsequently adopted in the CAPRI trunk. The former specification for land supply and transformation functions focused on agricultural land use and the transformation of agricultural land between arable land and grass land.
-)).+
  
 During the subsequent period, CAPRI was increasingly adapted to analyses of greenhouse gas (GHG) emission studies. Examples include CAPRI-ECC, GGELS, ECAMPA-X, AgCLim50-X, (European Commission, Joint Research Centre), ClipByFood (Swedish Energy Board), SUPREMA (H2020). This vein of research is very likely to gain in importance in the future. During the subsequent period, CAPRI was increasingly adapted to analyses of greenhouse gas (GHG) emission studies. Examples include CAPRI-ECC, GGELS, ECAMPA-X, AgCLim50-X, (European Commission, Joint Research Centre), ClipByFood (Swedish Energy Board), SUPREMA (H2020). This vein of research is very likely to gain in importance in the future.
  
 In order to improve land related climate gas modelling within CAPRI, it was deemed appropriate to (1) extend the land use modelled to //all// available land in the EU (i.e. not only agriculture), and (2) to explicitly model //transitions// between land use classes. The pioneering work was carried out within the TRUSTEE project((https://www.trustee-project.eu/ In order to improve land related climate gas modelling within CAPRI, it was deemed appropriate to (1) extend the land use modelled to //all// available land in the EU (i.e. not only agriculture), and (2) to explicitly model //transitions// between land use classes. The pioneering work was carried out within the TRUSTEE project((https://www.trustee-project.eu/
-)), but as always, an operational version emerged only after integrating efforts by researchers in several projects working at various institutions.+)), but as always, an operational version emerged only after integrating efforts by researchers in several projects working at various institutions. Within the SUPREMA project another important change in the depiction of land use change was made: the Markov chain approach was replaced by prespecifying the total land transitions as average transitions per year times the projection. This paper focusses on the theory applied while data and technical implementation are only briefly covered.
  
-This paper focusses on the theory applied while data and technical implementation are only briefly covered. 
  
-**A simple theory of land supply**+=== A simple theory of land supply ===
  
 Recall the dual methodological changes attempted in this paper: Recall the dual methodological changes attempted in this paper:
Line 428: Line 425:
 The model is somewhat complicated by the fact that land use classes in CAPRI are defined somewhat differently compared to the UNFCCC accounting and also in the land transition data set. Therefore, some of the land classes used in the land transitions are different from the ones used in the land supply model. In particular, “Other land”, “Wetlands” and “Pasture” are differently defined. To reconcile the differences, we assumed constant shares of the intersections of the different sets, as explained below. The model is somewhat complicated by the fact that land use classes in CAPRI are defined somewhat differently compared to the UNFCCC accounting and also in the land transition data set. Therefore, some of the land classes used in the land transitions are different from the ones used in the land supply model. In particular, “Other land”, “Wetlands” and “Pasture” are differently defined. To reconcile the differences, we assumed constant shares of the intersections of the different sets, as explained below.
  
-**Inner model – land transitions**+=== Inner model – transitions ===
  
 A vector of supply of land of various types could result from a wide range of different transitions. The inner model determines the matrix of land transitions that is “most likely”. The concept of “most likely” is formalized by assuming a joint density function for the land transitions, based on the historically observed transitions. The model then is to find the transition matrix that maximizes the joint density function. A vector of supply of land of various types could result from a wide range of different transitions. The inner model determines the matrix of land transitions that is “most likely”. The concept of “most likely” is formalized by assuming a joint density function for the land transitions, based on the historically observed transitions. The model then is to find the transition matrix that maximizes the joint density function.
 +
 +== Gamma density ==
  
 Since each transition is non-negative, but in principle unlimited upwards, we opted for a gamma density function, that has the support $\lbrack 0,\infty\rbrack$. For those that cannot immediately recall what the gamma density function looks like, and as entertainment for those that can, Figure 1 shows the graph of the density function for different parameters, all derived from an assumed mode of “1” and different assumed ratios “mode/standard deviations” (that we called “acc” for “accuracy” in the figure). Since each transition is non-negative, but in principle unlimited upwards, we opted for a gamma density function, that has the support $\lbrack 0,\infty\rbrack$. For those that cannot immediately recall what the gamma density function looks like, and as entertainment for those that can, Figure 1 shows the graph of the density function for different parameters, all derived from an assumed mode of “1” and different assumed ratios “mode/standard deviations” (that we called “acc” for “accuracy” in the figure).
Line 462: Line 461:
 $$\text{variance} = \frac{\alpha}{\beta^{2}}$$ $$\text{variance} = \frac{\alpha}{\beta^{2}}$$
  
-**Land use transitions as implemented in CAPRI**+== Annual transitions via Marcov chain in basic model ==
  
 The implementation in CAPRI differs from the above general framework in that it explicitly identifies the //annual// transitions in year t $T_{\text{lk}}^{t}$ from the initial $\text{LU}_{l}^{\text{initial}}$ land use to the final land use $\text{LU}_{k}$. This is necessary to identify the annual carbon effects occurring only in the final year in order to add them to the current GHG emissions, say from mineral fertiliser application in the final simulation year. If the initial year is the base year = 2008 and projection is for 2030, then the carbon effects related to the change from the 2008 $\text{LU}_{l}^{\text{initial}}$ to the final land use $\text{LU}_{k}$ (=$T_{\text{lk}}$in the above notation, without time index) refer to a period of 22 years that cannot reasonably be aggregated with the “running” non-CO2 effects from the final year 2030. Furthermore the historical time series used to determine the mode of the gamma density for the transitions also refer to annual transitions. The implementation in CAPRI differs from the above general framework in that it explicitly identifies the //annual// transitions in year t $T_{\text{lk}}^{t}$ from the initial $\text{LU}_{l}^{\text{initial}}$ land use to the final land use $\text{LU}_{k}$. This is necessary to identify the annual carbon effects occurring only in the final year in order to add them to the current GHG emissions, say from mineral fertiliser application in the final simulation year. If the initial year is the base year = 2008 and projection is for 2030, then the carbon effects related to the change from the 2008 $\text{LU}_{l}^{\text{initial}}$ to the final land use $\text{LU}_{k}$ (=$T_{\text{lk}}$in the above notation, without time index) refer to a period of 22 years that cannot reasonably be aggregated with the “running” non-CO2 effects from the final year 2030. Furthermore the historical time series used to determine the mode of the gamma density for the transitions also refer to annual transitions.
Line 486: Line 485:
 $$\Leftrightarrow 1 = \sum_{k}^{}P_{\text{lk}}$$ $$\Leftrightarrow 1 = \sum_{k}^{}P_{\text{lk}}$$
  
-So the simple condition is that probabilities have to add up to one (e_addUpTransMatrix in supply_model.gms). In this form the model is currently implemented in CAPRI.+So the simple condition is that probabilities have to add up to one (e_addUpTransMatrix in supply_model.gms).  
 + 
 +== Annual transitions if SUPREMA is active == 
 + 
 +As the use of the Marcov-chain approach allows the annual transitions to be explicit model variables that could be used to compute annual carbon effects but leads to computational limitations especially in the market model a new approach was developed under SUPREMA (i.e. if %supremaSup% == on) by re-specifying the total land transitions as average transitions per year times the projection horizon and by considering for the remaining class without land use change (on the diagonal of the land transition matrix) only the annual carbon effects per ha, relevant for the case of gains via forest management. 
 + 
 +The new accounting in the CAPRI global supply model may be explained as follows, starting from a calculation of the total GHG effects G over horizon h = t-s from total land transitions L<sub>lk</sub> and carbon effects per ha for the whole period e<sub>lk</sub>: 
 + 
 +$$G = Γ*h = \sum_{i,l}^{}{e_{\text{il}}^{}\text{L}}_{il}^{}$$ 
 + 
 +Where Γ collects the annual GHG effects that correspond to the total GHG effects divided by the time horizon G / h. These annual effects may be calculated as based on average annual transitions and annual effects for the remaining class as follows: 
 + 
 +$$Γ= \sum_{i,l}^{}{e_{\text{il}}^{}\text{L}}_{il}^{}/h = \sum_{i≠l}^{}{e_{\text{il}}^{}\text{Λ}}_{il}^{} 
 ++ \sum_{i}^{}{ε_{\text{ii}}^{}\text{L}}_{ii}^{} $$ 
 + 
 +Where Λ<sub>il</sub> = L<sub>il</sub> / h is the average land use change per year and ε<sub>ii</sub> is the annual carbon effect on a remaining class (relevant might be an annual increase due to growing forests while this will be zero for most effects based on IPCC default assumptions). 
 + 
 +Using these average annual transitions for true (off-diagonal) LUC we may compute the final classes as follows: 
 + 
 +$$ \text{LU}_{k,t} = \sum_{l}^{}{L_{\text{lk}}^{}\text{}}_{}^{} = \sum_{l≠k}^{}{Λ_{\text{lk}}^{}*h +\text{LU}_{kk}\text{}}_{}^{}$$ 
 + 
 +While adding up of shares (or probabilities) of LUC from class I to k over all receiving classes k continues to hold as stated above. It should be highlighted that the land use accounting implemented under SUPREMA avoids the need to explicitly trace the annual transitions in the form of a Markov chain and thereby economised on equations and variables. 
  
-**Outer model – land supply** 
  
 +=== Outer model – land supply ===
 The outer problem is defined as a maximization of the sum of land rents minus a quadratic cost term, subject to the first order optimality conditions of the inner problem: The outer problem is defined as a maximization of the sum of land rents minus a quadratic cost term, subject to the first order optimality conditions of the inner problem:
  
Line 577: Line 598:
 where $\delta$ is a small positive number and $I_{\text{ij}}$ entries of the identity matrix such that the factor $(1 - \delta I_{\text{ij}})$ shrinks the diagonal of the D-matrix, ensuring //strict// positive definiteness instead of //semi-//definiteness. We used $\delta = 0.05$. Furthermore, the Lagrange multiplier of the total land constraint, $\lambda$, was fixed at the weighted average of the rents $r_{i}$, i.e. $\lambda = \frac{\sum_{i}^{}{\text{LEV}L_{i}r_{i}}}{\sum_{i}^{}{\text{LEV}L_{i}}}$. Without the latter assumption, the parameters c and D are not uniquely identified. where $\delta$ is a small positive number and $I_{\text{ij}}$ entries of the identity matrix such that the factor $(1 - \delta I_{\text{ij}})$ shrinks the diagonal of the D-matrix, ensuring //strict// positive definiteness instead of //semi-//definiteness. We used $\delta = 0.05$. Furthermore, the Lagrange multiplier of the total land constraint, $\lambda$, was fixed at the weighted average of the rents $r_{i}$, i.e. $\lambda = \frac{\sum_{i}^{}{\text{LEV}L_{i}r_{i}}}{\sum_{i}^{}{\text{LEV}L_{i}}}$. Without the latter assumption, the parameters c and D are not uniquely identified.
  
-**Prior elasticities and area mappings**+==Prior elasticities and area mappings==
  
 The empirical evidence obtained in the TRUSTEE project applied to prior elasticities for land categories based on Corine Land Cover (CLC) data. These categories are also covered in the CAPRI database based on various sources (see the database section in the CAPRI documentation): The empirical evidence obtained in the TRUSTEE project applied to prior elasticities for land categories based on Corine Land Cover (CLC) data. These categories are also covered in the CAPRI database based on various sources (see the database section in the CAPRI documentation):
Line 635: Line 656:
 where 0 ≤ //share<sub>i,k</sub>// ≤ 1. where 0 ≤ //share<sub>i,k</sub>// ≤ 1.
  
-**Technical implementation**+===Technical implementation===
  
 The key equations corresponding to the approach explained above are collected in file supply_model.gms or the included files supply/declare_calibration_models_for_luc.gms and supply/declare_calibration_models_for_land_supply.gms. The declarations of parameters, variables, equations, models and even some sets only used in the calibration given in these files are included by the “supply_model.gms” only if “BASELINE==ON” or if it was a CAPREG base year task that was carried out. Loading of priors, initialisation of parameters and variables for the calibration as well as the organisation of solve attempts are handled in new sections of file “cal_land_nests.gms”, in turn called by the gams file “prep_cal.gms”. This implies that the land supply and land use change calibrations were inserted before the ordinary calibration of the supply models. The key equations corresponding to the approach explained above are collected in file supply_model.gms or the included files supply/declare_calibration_models_for_luc.gms and supply/declare_calibration_models_for_land_supply.gms. The declarations of parameters, variables, equations, models and even some sets only used in the calibration given in these files are included by the “supply_model.gms” only if “BASELINE==ON” or if it was a CAPREG base year task that was carried out. Loading of priors, initialisation of parameters and variables for the calibration as well as the organisation of solve attempts are handled in new sections of file “cal_land_nests.gms”, in turn called by the gams file “prep_cal.gms”. This implies that the land supply and land use change calibrations were inserted before the ordinary calibration of the supply models.
  
-The new land supply specification is only activated if the global variable %trustee_land%==on which may be set via the CAPRI GUI. In order to store the results of the calibration in a compact way that is compatible with the existing code, the existing parameter files “pmppar_XX.gdx” was used. The parameters of the land supply functions, called “c” and “D” above, were stored on two parameters “p_pmpCnstLandTypes” and “p_pmpQuadLandTypes”. As a new symbol (p_pmpCnstLandTypes) is introduced in an existing file, the first run of CAPRI after setting %trustee_land%==on may give errors if the file exists already but has been used with the previous land supply specification before. In this case it helps to delete or rename the old pmppar files.+//SupremaSup// should be active together with //trustee_land// to have smoother adjustments which may be set via the CAPRI GUI. In order to store the results of the calibration in a compact way that is compatible with the existing code, the existing parameter files “pmppar_XX.gdx” was used. The parameters of the land supply functions, called “c” and “D” above, were stored on two parameters “p_pmpCnstLandTypes” and “p_pmpQuadLandTypes”. As a new symbol (p_pmpCnstLandTypes) is introduced in an existing file, the first run of CAPRI after setting %trustee_land%==on may give errors if the file exists already but has been used with the previous land supply specification before. In this case it helps to delete or rename the old pmppar files.
  
 At this point, it should also be explained that rents for non-agricultural land types were entirely based on assumptions (a certain ratio to agricultural rents). As there were no plans to run scenarios with modified non-agricultural rents, these land rents //r// used in calibration for those land types were subtracted from the “c-paramter”, so that it is implicitly stored in p_pmpCnstLandTypes and enters the objective function through the PMP terms. This requires changes if the rents shall be modified or if non-agricultural production shall be included in some simplified form. At this point, it should also be explained that rents for non-agricultural land types were entirely based on assumptions (a certain ratio to agricultural rents). As there were no plans to run scenarios with modified non-agricultural rents, these land rents //r// used in calibration for those land types were subtracted from the “c-paramter”, so that it is implicitly stored in p_pmpCnstLandTypes and enters the objective function through the PMP terms. This requires changes if the rents shall be modified or if non-agricultural production shall be included in some simplified form.
Line 646: Line 667:
  
 More detailed explanations on the technical implementation are covered elsewhere, for example in the “Training material” included in the EcAMPA-4 deliverable D5. More detailed explanations on the technical implementation are covered elsewhere, for example in the “Training material” included in the EcAMPA-4 deliverable D5.
 +
 +Concerning the improvements made under SUPREMA from a technical perspective, the changes are merged to the trunk. The approach is controlled by globals in capmod\set_global_variables.gms. If the global variable %supremaSup% == on, the yearly transition rate p_lucAnnualFac_sup is calculated. If it is %supremaSup% == off, the old approach using the Marcov chain is used with the respective variable v_luYearly. The FOC-approach to calculate LUC as described above is standard and independent from if the global variable supremaSup is on or off.
 +
 +=== Emission Equations ===
 +
 +Under EcAMPA 3 and partly in earlier projects (inter alia EcAMPA 2) new modelling outputs have been developed for indicators without matching reporting infrastructure helping users to organise the additional information. This applied for example to
 +
 +1) Additional CAPRI results on land use results related to the complete area coverage, mappings to UNFCCC area categories and their transitions;
 +
 +2) The carbon effects linked to these land transitions.
 +
 +Furthermore, additional non-CO2 and CO2 related mitigation measures had been included under EcAMPA 3.
 +
 +The scenarios including the emission equations are only run if %ghgabatement% == on, otherwise emissions are only calculated and not simulated.
 +
 +The following emission equations have been implemented:
 +
 +^**Code**          ^**Description**                                                                                                ^
 +|GWPA              |Agricultural emissions                                                                                         |
 +|CH4ENT            |Methane emissions from enteric fermentation                                                                    |
 +|CH4MAN            |Methane emissions from manure management                                                                       |
 +|CH4RIC            |Methane emissions from rice production                                                                         |
 +|N2OMAN            |Direct nitrous oxide emissions stemming from manure management (only housing and storage)                      |
 +|N2OAPP            |Direct nitrous oxide emissions stemming from manure application on soils except grazings per animal activity   |
 +|N2OGRA            |Direct nitrous oxide emissions stemming from manure managment on grazings                                      |
 +|N2OSYN            |Direct nitrous oxide emissions from anorganic fertilizer application                                           |
 +|N2OCRO            |Direct nitrous oxide emissions from crop residues                                                              |
 +|N2OAMM            |Indirect nitrous oxide emissions from ammonia volatilisation                                                   |
 +|N2OLEA            |Indirect nitrous oxide emissions from leaching                                                                 |
 +|N2OHIS            |Direct nitrous oxide emissions from cultivation of histosols                                                   |
 +|GLUC              |Emissions related to indirect land use changes                                                                 |
 +|CO2BIO            |Carbon dioxide emissions from land use change due to losses of carbon in biomass and litter                    |
 +|CO2SOI            |Carbon dioxide emissions from land use change due to soil carbon losses                                        |
 +|CO2HIS\\ \\ CH4HIS|Carbon dioxide emissions from the cultivation of histosols\\ \\ Methane emissions from cultivation of histosols|
 +|CO2LIM\\ \\ CO2BUR|Carbon dioxide emissions from limestone and dolomit\\ \\ Carbon dioxide emissions from burning                 |
 +|CH4BUR            |Methane emissions from burning                                                                                 |
 +|N2OBUR            |Nitrous oxide emissions from burning                                                                           |
 +|N2OSOI            |N2O emissions from land use change due to soil carbon losses                                                   |
 +|GPRD              |Emissions related to the production of non-agricultural inputs to agriculture                                  |
 +|N2OPRD            |Nitrous oxide emissions during fertilizer production                                                           |
 +|O2PRD             |Carbon Dioxide emissions during fertilizer production                                                          |
  
 =====Premium module===== =====Premium module=====
Line 1215: Line 1277:
  
 The disadvantage of the behavioural functions above is the fact that they might generate non-positive values. That situation might be interpreted as a combination of prices where the marginal costs exceed marginal revenues. Accordingly, a fudging function is applied for supply, feed and (see below) processing demand which ensures strictly positive quantities. That fudging function is highly non-linear, and therefore only switched on on demand. The disadvantage of the behavioural functions above is the fact that they might generate non-positive values. That situation might be interpreted as a combination of prices where the marginal costs exceed marginal revenues. Accordingly, a fudging function is applied for supply, feed and (see below) processing demand which ensures strictly positive quantities. That fudging function is highly non-linear, and therefore only switched on on demand.
 +
 +====Land use, land use change and forestry (LULUCF) ====
 +
 +===LULUCF in the basic model ===
 +
 +Before SUPREMA LULUCF and area-based carbon accounting were not depicted in the global market model. Land demand was conceptually derived from maximising farmers profit. Land supply was represented with a function that links supply to agricultural land rents with an elasticity. Non-agricultural land use that complements farm land to give the total region area was disaggregated into forestry, built up areas (urban or “artificial” land) and a remaining “other land” category. There was neither a mapping of land use categories in the market model to the UNFCCC categories, nor a modelling of the transition matrix accompanied by a very limited product-based carbon accounting which was not in line with IPCC.
 +
 +The pre-SUPREMA specification may be described as follows.
 +
 +Agricultural outputs i (barley, wheat, beef ...) have land requirements LV<sub>i</sub> derived from production of these outputs (via yields that respond to prices according to yield elasticities). Adding up all land requirements gives total agricultural land (LT<sub>ag</sub>).
 +
 +$${LT}_{ag} = \sum_{i}^{}{{LV}_{i}(\mathbf{P},R_{ag})}$$
 +
 +Land demand depends on a vector of prices and the agricultural and rent R<sub>ag</sub> (treated separately). Total agricultural land is just one of several land types (l) that play a role:
 +
 +l = {ag, tc, pc, fd, no, fr, ur, ot, iw}, where
 +
 +ag = total agricultural land
 +
 +tc = temporary (non-fodder) crops
 +
 +pc = permanent crops
 +
 +fd = temporary fodder, permanent grassland and fallow land
 +
 +no = non-agricultural land
 +
 +fr = forest land
 +
 +ur = settlements, industrial, built up md any other artificial areas
 +
 +ot = other land
 +
 +iw = inland waters (exogenous)
 +
 +Matching with land demand there is a land supply function for total agricultural land
 +
 +$${LT}_{ag} = \alpha{R_{ag}}^{\beta}$$
 +
 +Given agricultural land and an exogenous region area as well as exogenous inland waters permits to compute total non-agricultural land residually:
 +
 +$${LT}_{no} = T - {LT}_{ag} - {LT}_{iw}$$
 +
 +This total non-agricultural land (beyond inland waters) is currently allocated to non-agricultural land types n = { fr, ur, ot} according to the shares of “intermediate” areas for non-agricultural land types:
 +
 +$${LT}_{n} = {LT}_{no}*{\widehat{LT}}_{n}/\sum_{n}^{}{\widehat{LT}}_{n}$$
 +
 +The “intermediate” areas in turn result from the change in the non-agricultural area against the baseline, considering elasticities that reflect the responsiveness of land types to imbalances:
 +
 +$${{\widehat{LT}}_{n} = LT}_{n}^{0}*\left( \frac{{LT}_{no}^{}}{{LT}_{no}^{0}} \right)^{\gamma_{n}}$$
 +
 +The concept of elasticities of land types to imbalances expresses the expectation that any disequilibrium in the land balance is very unlikely to be removed by changes in settlement area, and probably only to a small extent by changes in forest land and therefore most of all by changes in other land category.
 +
 +In spite of full consistency, we have changed this specification under SUPREMA, as 1) it is difficult to reconcile with welfare accounting, 2) it turned out that the scaling mechanism may dominate the planned responsiveness of land types, 3) the asymmetric specifications for supply of agricultural and non-agricultural land is ad-hoc and intransparent and 4) it does not link to standard empirical parameter estimation.
 +
 +== Land transitions via Gamma density and Marcov chain  ==
 +
 +The area-based carbon modelling and accounting requires the land transition matrix describing how an initial allocation of land uses (either from the base year or from an intermediate simulation year) is transformed into the currently simulated one. The transition matrix may be expressed in terms of absolute areas L<sub>jk</sub> changing from land use LU<sub>j</sub> in the initial year s into another land use LU<sub>k</sub> in the final year t or in terms of a transition matrix sh<sub>jk</sub> giving the share (probability in a Markov chain) of initial land use LU<sub>j</sub> converted into the final LU<sub>k</sub> over the whole horizon of (t-s):
 +
 +$${LU}_{k,t} = \sum_{j}^{}{{sh}_{jk}{LU}_{j,s} = \sum_{j}^{}L_{jk}}$$
 +
 +Where the shares (probabilities) have to add up to one:
 +
 +$$1 = \sum_{k}^{}{{sh}_{jk},\ \forall j}$$
 +
 +The total areas converted from initial land use j into final land use k over the horizon (t-s) are denoted L<sub>jk</sub> above. For those land transitions we would expect that the pattern of changes resembles that observed in the past, at least if total land uses LU<sub>j</sub> change similarly as in the past. This expectation corresponds to the most likely land transitions maximising a Gamma density, giving for each transition a corresponding FOC:
 +
 +$$\ \left( \lambda_{jk} - 1 \right)L_{jk}^{- 1} - \mu_{jk} + \tau_{k}^{} + \tau_{j}^{initial} = 0$$
 +
 +Where λ<sub>jk</sub> and μ<sub>jk</sub> are parameters related to the mode (determined from the database or baseline projection) and standard deviation (assumed = 1) of the Gamma density. The variables τ<sub>k</sub> and τ<sub>j</sub> are shadow values paired with the final year land use accounting from transition probabilities and the adding up condition for probabilities.
 +
 +The original specification for land transitions as used in the CAPRI supply models involve 6x(t-s) = 120 equations for a 20 years time horizon to represent a Markov chain of annual land transitions for each region. The advantage of this specification was that annual transitions were explicit model variables that could be used to compute annual carbon effects. These were hence comparable to annual non-CO2 effects related to agricultural production and could be added therefore to obtain total GHG effects from the LULUCF sector and non-CO2 GHGs. Having annual land transitions in each regional was also acceptable from a computational viewpoint for the relatively small regional supply models of CAPRI (about 1500 equations).
 +
 +However, in the global market model all regions (about 80 with agents like farmers, consumers or landowners) have to be solved simultaneously such that the additional equations and variables for the extended land use modelling and carbon accounting (addressed in the following section) could increase solution time beyond critical limits. Given that the standard market model already includes about 80000 equations the above framework was adjusted to give the land transitions in //one// step for the change from the initial years to the final year t, while still considering that we need annual carbon effects for comparability with the annual non-CO2 emissions. This has been achieved
 +
 +  * by re-specifying the total land transitions as average transitions per year times the projection horizon and
 +  * by considering for the remaining class without land use change (on the diagonal of the land transition matrix) only the annual carbon effects per ha, relevant for the case of gains via forest management.
 +
 +=== LULUCF and carbon accounting if SUPREMA is active ===
 +
 +Within the SUPREMA project two major changes were made:
 +
 +First, integration across spatial scales was improved: a) the land activity and LULUCF representation was extended to non-European countries and b) the product-based carbon accounting was replaced by an area-based carbon accounting.
 +
 +Second, the methodological approach was changed including a) statistical estimation of land use changes assuming a gamma density as in the supply model, b) re-specification of the total land transitions as average transitions per year times the projection horizon as in the supply model (replacement of the Markov chain approach) and c) representation of the disaggregated land supply in the market model through multinomial logit form to. These changes in the SUPREMA project allow for a more symmetric land use representation and carbon accounting between the supply models for European NUTS2 regions and in the global market model of CAPRI.
 +
 +== Multinomial logit function ==
 +
 +Under SUPREMA we have introduced a multinomial logit form for land supply of all major endogenous land types f = g = h = m = {ag, fr, ur, ot}. This approach is conceptually fully in line with land supply in the regional supply models. In this way we have integrated and replaced the above separate treatment of land supply for agricultural and non-agricultural land:
 +
 +$${LT}_{m} = {SH}_{m}*T$$
 +
 +where the area share of land type m is
 +
 +$${SH}_{g} = \frac{\exp\left( \delta_{g0} + \sum_{f}^{}{\delta_{gf}R_{f}} \right)}{\sum_{m}^{}{\exp\left( \delta_{m0} + \sum_{f}^{}{\delta_{mf}R_{f}} \right)}}$$
 +
 +and the elasticity of share g (and due to constant region area also land type LT<sub>g</sub>) with respect to rent R<sub>h</sub> may be derived as
 +
 +$$\varepsilon_{gh} = \frac{\partial{SH}_{g}}{\partial R_{h}}\frac{R_{h}}{{SH}_{g}} = R_{h}\left( \delta_{gh} - \sum_{m}^{}{\delta_{mh}{SH}_{m}} \right)$$
 +
 +which permits to make use of the same empirical information (on elasticities of agricultural land supply) and assumptions (on the ranking of responsiveness of non-agricultural areas) that have been used so far in the pre-SUPREMA version. For this purpose, a calibration problem has been set up that minimises weighted squared differences to the starting values by modifying parameters δ<sub>mh</sub>. Due to its symmetric treatment of all major land uses the system also includes supply elasticities for non-agricultural areas. In “standard” scenarios these are unlikely to play a relevant role. The key parameters are the “cross-rent” elasticities of non-agricultural areas with respect to agricultural rents as these are steering now which non-agricultural areas are increasing if agricultural area declines and vice versa. However, in the context of global carbon price scenarios also the supply elasticities of non-agricultural areas play an important role even though we do not introduce assumptions on changing prices of urban land, forest land or other land. This is because a global carbon price creates an endogenous mark-up for land rental prices of non-agricultural areas that reflect the value of the carbon effects from changing land use. In particular the rental price of forest land R<sub>fr</sub> is strongly reduced and might become negative in scenarios with high carbon prices.
 +
 +== Spatial extension ==
 +
 +Under SUPREMA the land use categories of the market model are mapped to the UNFCCC categories. The mapping of market model land types LT<sub>l</sub> to UNFCCC land use LU<sub>k</sub> will rely on the most recent historical shares ϕ<sub>kl</sub> of UNFCCC land use k in CAPRI land type l:
 +
 +$${LU}_{k} = \sum_{l}^{}{\varphi_{kl}{LT}_{l}}$$
 +
 +These shares are trivially zero or one in case that certain land types like “temporary non-fodder crops” (tc) and permanent crops (pc) are exclusively mapped to one UNFCCC category (cropland). The remainder to total cropland derives from temporary fodder and fallow land which is a fraction of total fodder area with the remainder being (productive) permanent grassland. The allocation of “other land” (ot) to grassland (ϕ<sub>glot</sub>), wetland (ϕ<sub>wlot</sub>) and residual land (ϕ<sub>rlot</sub>) may occur as in the European database but required that some CAPRI code in use for the supply models was transferred into the context of the extended global market model.
 +
 +== Land transitions as average transitions per year times the projection horizon  ==
 +
 +The new accounting in the CAPRI global market model may be explained as follows, starting from a calculation of the total GHG effects G over horizon h = t-s from total land transitions L<sub>jk</sub> and carbon effects per ha for the whole period e<sub>jk</sub>:
 +
 +$$G = \Gamma \bullet h = \sum_{i,j}^{}{e_{ij}L_{ij}}$$
 +
 +Where Γ collects the annual GHG effects that correspond to the total GHG effects divided by the time horizon G / h. These annual effects may be calculated as based on average annual transitions and annual effects for the remaining class as follows:
 +
 +$
 +$$$\Gamma = \sum_{i,j}^{}{e_{ij}L_{ij}/h} = \sum_{i \neq j}^{}{e_{ij}\Lambda_{ij}} + \sum_{i}^{}{\varepsilon_{ii}L_{ii}}$$
 +
 +Where Λ<sub>ij</sub> = L<sub>ij</sub> / h is the average land use change per year and ε<sub>ii</sub> is the annual carbon effect on a remaining class (relevant might be an annual increase due to growing forests while this will be zero for most effects based on IPCC default assumptions).
 +
 +Using these average annual transitions for true (off-diagonal) LUC we may compute the final classes as follows:
 +
 +$${LU}_{k,t} = \sum_{j}^{}L_{jk} = \sum_{j \neq k}^{}\Lambda_{jk} \bullet h + L_{kk}$$
 +
 +While adding up of shares (or probabilities) of LUC from class I to k over all receiving classes k continues to hold as stated above. It should be highlighted that the land use accounting implemented under SUPREMA avoids the need to explicitly trace the annual transitions in the form of a Markov chain and thereby economised on equations and variables. In this form LUC by CAPRI region and the associated accounting of carbon effects turned out computationally feasible even though the number of equations in the global market model increased from about 78000 to about 83000. Apart from feasibility the format above also permitted to retain the typical CAPRI accounting identity that some total “quantity” (“GROF”) should be computable as the effects “per activity” times activity levels. It was therefore also adopted in the CAPRI supply models.
 +
 +=== Technical aspects ===
 +
 +Concerning the improvements made under SUPREMA from a technical perspective, the changes are merged to the trunk. The approach is controlled by globals in capmod\set_global_variables.gms. If the global variable %supremaMrk% == on, the yearly transition rate p_lucAnnualFac_sup is calculated, land activity is expaned to non-european countries and the multinomial logit form approach is used to model land supply responsiveness. If it is %supremaMrk% == off, the old approach using the Marcov chain is used with the respective variable v_luYearly (for details on Marcov chain approach see [[module_for_agricultural_supply_at_regional_level|supply model description]]). The FOC-approach to calculate LUC as described above is standard and independent from if the global variable supremaMrk is on or off.
 +
 +=== Carbon accounting ===
 +
 +A last recent change concerns the transfer of the existing carbon accounting equations from the supply model to the global market model. These equations run if the global variable %supremaMrk% == on. More information on the equations can be found in the [[module_for_agricultural_supply_at_regional_level|description of the supply model]]. The equations concerned are, indicated with their present “CAPRI names” in the supply models and plus “Mrk” in the market model:
 +
 +Table 3. Equations concerning mitigation modelling in CAPRI
 +
 +^__Supply model__^__Market model__  ^
 +|GWPCO2BIO_      |GWPCO2BIOMrk_(RMS)|
 +|GWPCO2SOI_      |GWPCO2SOIMrk_(RMS)|
 +|GWPN2OSOI_      |GWPN2OSOIMrk_(RMS)|
 +|GWPCO2BUR_      |GWPCO2BURMrk_(RMS)|
 +|GWPCH4BUR_      |GWPCH4BURMrk_(RMS)|
 +|GWPN2OBUR_      |GWPN2OBURMrk_(RMS)|
 +|GWPCO2HIS_      |GWPCO2HISMrk_(RMS)|
 +|GWPCH42HIS_     |GWPCH4HISMrk_(RMS)|
 +
 +As most CAPRI regions combine only a small number of climate zones and we may assume for an illustrative calculation three (out of 9). In this case we would have 22 additional equations for carbon accounting per region and 1738 equations in total (on top of those for land use modelling mentioned above), confirming the order of magnitude for the additional equations that was mentioned above.
 +
 +For the technical coefficients we could rely on the FAO data compiled for the implementation of LULUCF accounting in the supply models. Here they served often only a fall-back solution in case that some European dataset was missing, but for the global market model the FAO data are often the only source of data readily available.
 +
 +There is one new element required for the planned implementation of a carbon tax deriving from land use and land use changes: In the supply models the tax is simply added as a cost element in the existing income accounting for the regional farms to make it effective. In the global market model there is no explicit income accounting. The carbon tax levied so far on non-CO2 emissions has been translated therefore into a tax on outputs, depending on the product based non-CO2 emission factors of outputs that may be changed implicitly. For LULUCF an explicit tax on land use has been introduced. This required to treat land demand and land supply for temporary and permanent crops and fodder as separate qualities with distinct rental prices for market clearing. 
  
 ====Behavioural equations for final demand==== ====Behavioural equations for final demand====
scenario_simulation.1642413054.txt.gz · Last modified: 2022/11/07 10:23 (external edit)
Except where otherwise noted, content on this wiki is licensed under the following license: CC0 1.0 Universal
CC0 1.0 Universal Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki