scenario_simulation
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
scenario_simulation [2020/04/26 05:37] – matsz | scenario_simulation [2023/09/08 11:51] – [2.1. Multinomial logit function] massfeller | ||
---|---|---|---|
Line 25: | Line 25: | ||
**Figure 13: Link of modules in CAPRI** | **Figure 13: Link of modules in CAPRI** | ||
- | {{::figure13.png? | + | {{::figure_13.png? |
=====Module for agricultural supply at regional level===== | =====Module for agricultural supply at regional level===== | ||
Line 117: | Line 117: | ||
Where “asym” is the land asymptote, i.e. the maximal amount of economically usable agricultural area in a region when the agricultural land rent goes towards infinity. For an application where the land market is used see Renwick et al. (2013). | Where “asym” is the land asymptote, i.e. the maximal amount of economically usable agricultural area in a region when the agricultural land rent goes towards infinity. For an application where the land market is used see Renwick et al. (2013). | ||
- | Set aside policies have changed frequently during CAP reforms. The recent specification is covered in the context of the premium modelling in Section [[Premium module]]. The obligatory set-aside restriction introduced by the McSharry reform 1992 and valid until the implementation of the Luxembourg compromise of June 2003 has been explicitly modelled through this equation: | + | Set aside policies have changed frequently during CAP reforms. The recent specification is covered in the context of the premium modelling in Section [[scenario simulation#Premium module]]. The obligatory set-aside restriction introduced by the McSharry reform 1992 and valid until the implementation of the Luxembourg compromise of June 2003 has been explicitly modelled through this equation: |
\begin{align} | \begin{align} | ||
Line 337: | Line 337: | ||
\end{matrix} | \end{matrix} | ||
\right] | \right] | ||
- | |||
\end{split} | \end{split} | ||
\end{align} | \end{align} | ||
Line 351: | Line 350: | ||
\end{align} | \end{align} | ||
- | The scaling factor to map from the legal quota legalquotA (as the B quota has been eliminated in the sugar reform, it holds that \(q^A = q^{A+B}) \)to the behavioural quota qA depends on the projected sugar beet sales quantity in the calibration point \(NETTRD_{SUGB}^{cal})\ : For a country with a high over quota production (say 40%) we would obtain a scaling factor of 1.31, such that this producer will behave like a moderate C-sugar producer: responsive to both the C-beet prices as well as to the quota beet price (and the legal quotas). Without this scaling factor, producers with significant over quota p | + | The scaling factor to map from the legal quota legalquotA (as the B quota has been eliminated in the sugar reform, it holds that \(q^A = q^{A+B} \) )to the behavioural quota qA depends on the projected sugar beet sales quantity in the calibration point \( NETTRD_{SUGB}^{cal} \) : For a country with a high over quota production (say 40%) we would obtain a scaling factor of 1.31, such that this producer will behave like a moderate C-sugar producer: responsive to both the C-beet prices as well as to the quota beet price (and the legal quotas). Without this scaling factor, producers with significant over quota p |
===Update note=== | ===Update note=== | ||
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 |
- | -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:// | + | - 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, | + | |
====Calibration of the regional programming models==== | ====Calibration of the regional programming models==== | ||
Line 393: | Line 392: | ||
Y_{j, | Y_{j, | ||
\end{equation} | \end{equation} | ||
+ | |||
+ | |||
+ | ==== LULUCF in the supply model of CAPRI ==== | ||
+ | |||
+ | === 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 ((https:// | ||
+ | )) and SUPREMA ((https:// | ||
+ | |||
+ | 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), | ||
+ | )), 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. | ||
+ | |||
+ | |||
+ | === A simple theory of land supply === | ||
+ | |||
+ | Recall the dual methodological changes attempted in this paper: | ||
+ | |||
+ | - Extend land use modelling to the entire land area, and | ||
+ | - Explicitly model transitions between each pair of land uses | ||
+ | |||
+ | In order to keep things as simple as possible, we opted for a theory where the decision of how much land to allocate to each use is independent of the explicit transitions between classes. This separation of decisions is simplifying the theoretical derivations, | ||
+ | |||
+ | The land supply and transformation model developed here is a bilevel optimization model. At the higher level (sometimes termed the //outer problem//), the land owner decides how much land to allocate to each aggregate land use based on the rents earned in each use and a set of parameters capturing the costs required in order to ensure that the land is available to the intended use. At the lower level (sometimes termed the //inner problem//), the transitions between land classes are modelled, with the condition that the total land needs of the outer problem are satisfied. The inner problem is modelled as a stochastic process involving no explicit economic model. | ||
+ | |||
+ | For the outer problem, i.e. the land owner’s problem, we propose a quadratic objective function that maximizes the sum of land rents minus a dual cost function. The parameters of the dual cost function were specified in two steps: | ||
+ | |||
+ | - A matrix of land supply elasticities was estimated (by TRUSTEE partner Jean Saveur Ay, CESEAR, Dijon (JSA). This estimation might be updated in future work or replaced with other sources for elasticities. | ||
+ | - The parameters of the dual cost function are specified so that the supply behaviour replicates the estimated elasticities as closely as possible while exactly replicating observed/ | ||
+ | |||
+ | 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, | ||
+ | |||
+ | === 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, | ||
+ | |||
+ | == Gamma density == | ||
+ | |||
+ | Since each transition is non-negative, | ||
+ | |||
+ | {{: | ||
+ | |||
+ | Figure 1: Gamma density graph for mode=1 and various standard deviations. “acc”=" | ||
+ | |||
+ | Let $i$ denote land use classes in CAPRI definition, whereas //l// and //k// are land uses in UNFCCC classification. Let $\text{LU}_{k}$ be total land use after transitions and $\text{LU}_{l}^{\text{initial}}$ be land use before transitions. Furthermore, | ||
+ | |||
+ | $${\max_{T_{\text{lk}}}{\log{\prod_{\text{lk}}^{}{f\left( T_{\text{lk}}|\alpha_{\text{lk}}, | ||
+ | |||
+ | $$\Rightarrow \max_{T_{\text{lk}}}\sum_{\text{lk}}^{}\left\lbrack \left( \alpha_{\text{lk}} - 1 \right)\log T_{\text{lk}} - \beta_{\text{lk}}T_{\text{lk}} \right\rbrack$$ | ||
+ | |||
+ | subject to | ||
+ | |||
+ | $$\text{LU}_{k} - \sum_{l}^{}T_{\text{lk}} = 0 \; \left\lbrack \tau_{k} \right\rbrack$$ | ||
+ | |||
+ | $$\text{LU}_{l}^{\text{initial}} - \sum_{k}^{}T_{\text{lk}} = 0\; | ||
+ | |||
+ | $$\text{LU}_{k} - \sum_{i}^{}{\text{shar}e_{\text{ki}}\text{LEV}L_{i}} = 0$$ | ||
+ | |||
+ | The last equation is needed to convert land use in UNFCCC classification to land use in CAPRI classification, | ||
+ | |||
+ | $$\ \left( \alpha_{\text{lk}} - 1 \right)T_{\text{lk}}^{- 1} - \beta_{\text{lk}} + \tau_{k}^{} + \tau_{l}^{\text{initial}} = 0$$ | ||
+ | |||
+ | The parameters $\alpha$ and $\beta$ of the gamma density function were computed by assuming that (i) the observed transitions are the mode of the density, and (ii) the standard deviation equals the mode. Then the parameters are obtained by solving the following quadratic system: | ||
+ | |||
+ | $$\text{mode} = \frac{\alpha - 1}{\beta}$$ | ||
+ | |||
+ | $$\text{variance} = \frac{\alpha}{\beta^{2}}$$ | ||
+ | |||
+ | == 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. | ||
+ | |||
+ | Initially the problem to link total to annual transitions has been solved by assuming a linear time path from the initial to the final period, but this was criticised as being an inconsistent time path (by FW). Ultimately the time path has been computed therefore in the supply model in line with a static Markov chain with constant probabilities $P_{\text{lk}}$ such that both land use $\text{LU}_{l}^{t}$ as well as transitions $T_{\text{lk}}^{t}$ in absolute ha require a time index (e_luOverTime in supply_model.gms). | ||
+ | |||
+ | $$\text{LU}_{k}^{t} - \sum_{l}^{}{P_{\text{lk}}\text{LU}_{l}^{t - 1}} = 0\ ,\ t = \{ 1,\ldots s\}$$ | ||
+ | |||
+ | Where $\text{LU}_{k}^{s}$ is the final land use in the simulation year s and $\text{LU}_{k}^{0} = \text{LU}_{k}^{\text{iniital}}$ is the initial land use. The transitions in ha in any year may be recovered from previous years land use and the annual (and constant) transition probabilities (e_LUCfromMatrix in supply_model.gms). | ||
+ | |||
+ | $$T_{\text{lk}}^{t} = P_{\text{lk}}*\text{LU}_{l}^{t - 1}$$ | ||
+ | |||
+ | The absolute transitions may enter the carbon accounting (ignored here) and if we substitute the last period’s transitions we are back to the condition for consistent land balancing in the final period from above: | ||
+ | |||
+ | $$\text{LU}_{k}^{s} = \sum_{l}^{}{P_{\text{lk}}\text{LU}_{l}^{s - 1}} = \sum_{l}^{}T_{\text{lk}}^{s}$$ | ||
+ | |||
+ | When using the transition probabilities in the consistency condition for initial land use we obtain | ||
+ | |||
+ | $$\text{LU}_{l}^{\text{initial}} - \sum_{k}^{}T_{\text{lk}}^{1} = 0$$ | ||
+ | |||
+ | $$\Longleftrightarrow \text{LU}_{l}^{\text{initial}} = \sum_{k}^{}{P_{\text{lk}}^{}\text{LU}}_{l}^{\text{iniital}}$$ | ||
+ | |||
+ | $$\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). | ||
+ | |||
+ | == 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< | ||
+ | |||
+ | $$G = Γ*h = \sum_{i, | ||
+ | |||
+ | 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, | ||
+ | + \sum_{i}^{}{ε_{\text{ii}}^{}\text{L}}_{ii}^{} $$ | ||
+ | |||
+ | Where Λ< | ||
+ | |||
+ | Using these average annual transitions for true (off-diagonal) LUC we may compute the final classes as follows: | ||
+ | |||
+ | $$ \text{LU}_{k, | ||
+ | |||
+ | 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 === | ||
+ | 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: | ||
+ | |||
+ | $$\max{\sum_{i}^{}{\text{LEV}L_{i}r_{i}} - \sum_{i}^{}{\text{LEV}L_{i}c_{i}} - \frac{1}{2}\sum_{\text{ij}}^{}{\text{LEV}L_{i}D_{\text{ij}}\text{LEV}L_{j}}}$$ | ||
+ | |||
+ | subject to, | ||
+ | |||
+ | $$\text{LU}_{k} - \sum_{i}^{}{\text{shar}e_{\text{ki}}\text{LEV}L_{i}} = 0$$ | ||
+ | |||
+ | $$\text{LU}_{k} - \sum_{l}^{}T_{\text{lk}} = 0\; | ||
+ | |||
+ | $$\text{LU}_{l}^{\text{initial}} - \sum_{k}^{}T_{\text{lk}} = 0\; | ||
+ | |||
+ | $$\ \left( \alpha_{\text{lk}} - 1 \right)T_{\text{lk}}^{- 1} - \beta_{\text{lk}} + \tau_{k}^{} + \tau_{l}^{\text{initial}} = 0$$ | ||
+ | |||
+ | The parameters of the inner model **α** and **β// | ||
+ | |||
+ | There are a few methodological and numerical challenges to overcome. In particular, we need to (i) analytically derive $\mathbf{\eta}\left( \mathbf{c}, | ||
+ | |||
+ | $$\sum_{i}^{}{\text{LEV}L_{i}} - \sum_{l}^{}{LU_{l}^{\text{initial}}} = 0$$ | ||
+ | |||
+ | Note that the second sum is a constant. This simplification is based on the observation that the land transitions don’t appear in the objective function of the outer problem, so that all solutions to the inner problems are equivalent from the perspective of the outer problem, and that any land use vector that preserves the initial land endowment is a feasible solution to the inner problem. | ||
+ | |||
+ | Next, we formulate the first order condition (FOC) of the modified outer problem to obtain land use as an implicit function of the parameters, $F\left( LEVL, | ||
+ | |||
+ | The first order conditions, and the implicit function, become | ||
+ | |||
+ | $$F\left( LEVL, | ||
+ | \frac{\partial\mathcal{L}}{\partial LEVL_{i}} = & r_{i} - c_{i} - \sum_{j}^{}{D_{\text{ij}}\text{LEV}L_{j}} - \lambda & = 0 \\ | ||
+ | \frac{\partial\mathcal{L}}{\partial\lambda} = & \sum_{i}^{}{\text{LEV}L_{i}} - \sum_{l}^{}{LU_{l}^{\text{initial}}} & = 0 \\ | ||
+ | \end{bmatrix}$$ | ||
+ | |||
+ | In order to apply the implicit function theorem((Recall that the implicit function theorem states that if F(x,p) = 0, then dx/dp = -[dF/ | ||
+ | )) we need to differentiate the FOC once w.r.t. the variables $\text{LEV}L_{i}$ and $\lambda$ and once with respect to the parameter of interest, $r_{j}$, invert the former and take the negative of the matrix product. If (currently) irrelevant parameter are omitted, the following matrix of $(N + 1) \times (N + 1)$ is obtained (the “+1” is the uninteresting derivative of total land rent $\lambda$ with respect to individual land class rent $r_{i}$) | ||
+ | |||
+ | $$\left\lbrack \frac{\partial LEVL}{\partial r} \right\rbrack = - \left\lbrack D_{LEVL, | ||
+ | |||
+ | $$\begin{bmatrix} | ||
+ | \frac{\partial LEVL}{\partial r} \\ | ||
+ | \frac{\partial\lambda}{\partial r} \\ | ||
+ | \end{bmatrix} = - \begin{bmatrix} | ||
+ | \frac{\partial F}{\partial LEVL} & \frac{\partial F}{\partial\lambda} \\ | ||
+ | \end{bmatrix}\left\lbrack \frac{\partial F}{\partial r} \right\rbrack$$ | ||
+ | |||
+ | Carrying out the differentiation specifically for land rent // | ||
+ | |||
+ | $$\begin{bmatrix} | ||
+ | \frac{\partial LEVL_{i}}{\partial r_{j}} \\ | ||
+ | \frac{\partial\lambda}{\partial r_{j}} \\ | ||
+ | \end{bmatrix} = - \begin{bmatrix} | ||
+ | \left\lbrack {- D}_{\text{ij}} \right\rbrack & - 1 \\ | ||
+ | - 1' & 0 \\ | ||
+ | \end{bmatrix}^{- 1}\begin{bmatrix} | ||
+ | I \\ | ||
+ | 0 \\ | ||
+ | \end{bmatrix}$$ | ||
+ | |||
+ | Discarding the last row of the resulting $(N + 1) \times N$ matrix finally lets us compute the elasticity as | ||
+ | |||
+ | $$\left\lbrack \eta_{\text{ij}} \right\rbrack = \left\lbrack \frac{\partial LEVL_{i}}{\partial r_{j}} \right\rbrack\left\lbrack \frac{r_{j}}{\text{LEV}L_{i}} \right\rbrack$$ | ||
+ | |||
+ | In the estimation, we assumed that the prior elasticity matrix is the mode of a density where each entry were independently distributed. Furthermore, | ||
+ | |||
+ | $$\max_{\eta, | ||
+ | |||
+ | subject to | ||
+ | |||
+ | $$\left\lbrack \frac{\partial LEVL_{i}}{\partial r_{j}} \right\rbrack = - \begin{bmatrix} | ||
+ | \left\lbrack {- D}_{\text{ij}} \right\rbrack & - 1 \\ | ||
+ | - 1' & 0 \\ | ||
+ | \end{bmatrix}^{- 1}\begin{bmatrix} | ||
+ | I \\ | ||
+ | 0 \\ | ||
+ | \end{bmatrix}$$ | ||
+ | |||
+ | $$\left\lbrack \eta_{\text{ij}} \right\rbrack = \left\lbrack \frac{\partial LEVL_{i}}{\partial r_{j}} \right\rbrack\left\lbrack \frac{r_{j}}{\text{LEV}L_{i}} \right\rbrack$$ | ||
+ | |||
+ | $$\begin{matrix} | ||
+ | & r_{i} - c_{i} - \sum_{j}^{}{D_{\text{ij}}\text{LEV}L_{j}} - \lambda & = 0 \\ | ||
+ | & \sum_{i}^{}{\text{LEV}L_{i}} - \sum_{l}^{}{LU_{l}^{\text{initial}}} & = 0 \\ | ||
+ | \end{matrix}$$ | ||
+ | |||
+ | and the curvature constraint using a stricter variant of the Cholesky factorization | ||
+ | |||
+ | $$D_{\text{ij}}\left( 1 - \delta I_{\text{ij}} \right) = \sum_{k}^{}{U_{\text{ki}}U_{\text{kj}}}$$ | ||
+ | |||
+ | 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 // | ||
+ | |||
+ | ==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 introduction has mentioned already three systems of area categories that need to be distinguished. The first one is the set of area aggregates with good coverage in statistics that has been investigated recently by JS Ay (2016), in the following “JSA”: | ||
+ | |||
+ | $$\text{LEVL} = \left\{ \text{ARAC}, | ||
+ | |||
+ | Where | ||
+ | |||
+ | ARAC = arable crops | ||
+ | |||
+ | FRUN = perennial crops | ||
+ | |||
+ | GRAS = permanent grassland | ||
+ | |||
+ | FORE = forest | ||
+ | |||
+ | ARTIF = artificial surfaces (settlements, | ||
+ | |||
+ | OLND = other land | ||
+ | |||
+ | The above categories are matching reasonably well with the definitions in JSA. A mismatch exists in the classification of paddy (part of ARAC in CAPRI but in the perennial group in JSA) and terrestrial wetlands (part of OLND in CAPRI and a separate category in JSA). Inland waters are considered exogenous in CAPRI and hence not included in the above set LEVL. | ||
+ | |||
+ | For carbon accounting we need to identify the six LU classes from IPCC recommendations and official UNFCCC reporting: | ||
+ | |||
+ | $$LU = \left\{ \text{CROP}, | ||
+ | |||
+ | which is typically indexed below with “l” or “k” ∈ LU and where | ||
+ | |||
+ | CROP = crop land (= sum of arable crops and perennial crops) | ||
+ | |||
+ | GRSLND = grassland in IPCC definition (includes some shrub land and other “nature land”, hence GRSLND> | ||
+ | |||
+ | WETLND = wetland (includes inland waters but also terrestrial wetlands) | ||
+ | |||
+ | RESLND = residual land is that part of OLND not allocated to grassland or wetland, hence RESLND< | ||
+ | |||
+ | FORE = forest | ||
+ | |||
+ | ARTIF = artificial surfaces | ||
+ | |||
+ | In the CAPRI database, in particular for its technical base year, we have estimated an allocation of other land OLND into its components attributable to the UNFCCC classes GRSLND, | ||
+ | |||
+ | $$\text{OLND}^{0} = {\text{OLND}G}^{0} + {\text{OLND}W}^{0} + {\text{OLND}R}^{0}$$ | ||
+ | |||
+ | Lacking better options to make the link between sets LEVL (activity level aggregates) and LU (UNFCCC classes, technically in CAPRI code: set “LUclass”) we will assume that these shares are fixed and may estimate the “mixed” LU areas from activity level aggregates as follows | ||
+ | |||
+ | ^// | ||
+ | |WETLND | ||
+ | |RESLND | ||
+ | |||
+ | which means that the mapping from set LEVL to set LU only uses some fixed shares of LEVL areas that are mapped to a certain LU: | ||
+ | |||
+ | $$LU_k=\sum_i{\text{share}_{\text{i, | ||
+ | |||
+ | where 0 ≤ // | ||
+ | |||
+ | ===Technical implementation=== | ||
+ | |||
+ | The key equations corresponding to the approach explained above are collected in file supply_model.gms or the included files supply/ | ||
+ | |||
+ | // | ||
+ | |||
+ | 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”, | ||
+ | |||
+ | Furthermore, | ||
+ | |||
+ | 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, | ||
+ | |||
+ | === 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, | ||
+ | |||
+ | 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** | ||
+ | |GWPA |Agricultural emissions | ||
+ | |CH4ENT | ||
+ | |CH4MAN | ||
+ | |CH4RIC | ||
+ | |N2OMAN | ||
+ | |N2OAPP | ||
+ | |N2OGRA | ||
+ | |N2OSYN | ||
+ | |N2OCRO | ||
+ | |N2OAMM | ||
+ | |N2OLEA | ||
+ | |N2OHIS | ||
+ | |GLUC |Emissions related to indirect land use changes | ||
+ | |CO2BIO | ||
+ | |CO2SOI | ||
+ | |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 | ||
+ | |N2OBUR | ||
+ | |N2OSOI | ||
+ | |GPRD |Emissions related to the production of non-agricultural inputs to agriculture | ||
+ | |N2OPRD | ||
+ | |O2PRD | ||
=====Premium module===== | =====Premium module===== | ||
Line 452: | Line 767: | ||
**Figure 14: Example of technical implementation of a premium scheme** | **Figure 14: Example of technical implementation of a premium scheme** | ||
- | {{:figure14.png? | + | {{:figure_14.png? |
The sets of payments, exemplified by DPGRCU in the figure, and the activity groups, exemplified by PGGRCU and PGPROT are defined in the file policy/ | The sets of payments, exemplified by DPGRCU in the figure, and the activity groups, exemplified by PGGRCU and PGPROT are defined in the file policy/ | ||
Line 473: | Line 788: | ||
**Figure 15: General flow of logic of CAPRI model as regards premiums** | **Figure 15: General flow of logic of CAPRI model as regards premiums** | ||
- | {{::fiugre15.png? | + | {{::figure_15.png? |
Generally, all attributes for a premium scheme are mapped down in space, e.g. from EU27 to EU 27 member states, from countries to NUTS1 regions inside the country, from there to the NUTS2 regions inside the NUTS1, and from NUTS2 regions to the farm types in a NUTS2 region (see // | Generally, all attributes for a premium scheme are mapped down in space, e.g. from EU27 to EU 27 member states, from countries to NUTS1 regions inside the country, from there to the NUTS2 regions inside the NUTS1, and from NUTS2 regions to the farm types in a NUTS2 region (see // | ||
Line 562: | Line 877: | ||
**Figure 16: General way of SFP implementation in CAPRI** | **Figure 16: General way of SFP implementation in CAPRI** | ||
- | {{::figure16.png? | + | {{::figure_16.png? |
In opposite to the reforms until Agenda 2000, there are hence in most cases not longer premium rates or individual ceilings in hectares found in legal texts. Rather, these are calculated by the model itself from the decoupled part of the “old” Mac Sharry and Agenda 2000 premiums which introduces additional complexity in the model code. | In opposite to the reforms until Agenda 2000, there are hence in most cases not longer premium rates or individual ceilings in hectares found in legal texts. Rather, these are calculated by the model itself from the decoupled part of the “old” Mac Sharry and Agenda 2000 premiums which introduces additional complexity in the model code. | ||
Line 680: | Line 995: | ||
**Figure 17: The logic of the CAP 2014-2020 reference policy as implemented in the CAPRI policy module** | **Figure 17: The logic of the CAP 2014-2020 reference policy as implemented in the CAPRI policy module** | ||
- | {{:figure17.png? | + | {{:figure_17.png? |
===Tradable Single Premium Scheme entitlements=== | ===Tradable Single Premium Scheme entitlements=== | ||
Line 816: | Line 1131: | ||
Finally, an extensification effect to the AE payments is introduced using the possibility to make technological variants differently eligible. | Finally, an extensification effect to the AE payments is introduced using the possibility to make technological variants differently eligible. | ||
- | {{: | + | {{: |
{{: | {{: | ||
Line 843: | Line 1159: | ||
{{: | {{: | ||
- | Currently, the following budget categories are supported (see ‘// | + | Currently, the following budget categories are supported (see ‘// |
{{: | {{: | ||
Line 922: | Line 1238: | ||
**Figure 18: Land supply curve examples** | **Figure 18: Land supply curve examples** | ||
- | {{::figure18.png? | + | {{::figure_18.png? |
In order to parameterize the land demand function, information about yield and supply elasticities is used. The marginal reaction of land to a marginal change in one of the prices is defined as the total supply effect minus the yield effect: | In order to parameterize the land demand function, information about yield and supply elasticities is used. The marginal reaction of land to a marginal change in one of the prices is defined as the total supply effect minus the yield effect: | ||
Line 961: | 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, | 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, | ||
+ | |||
+ | ====Land use, land use change and forestry (LULUCF) ==== | ||
+ | |||
+ | === 1. 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< | ||
+ | |||
+ | $${LT}_{ag} = \sum_{i}^{}{{LV}_{i}(\mathbf{P}, | ||
+ | |||
+ | Land demand depends on a vector of prices and the agricultural and rent R< | ||
+ | |||
+ | 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, | ||
+ | |||
+ | 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}/ | ||
+ | |||
+ | 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, | ||
+ | |||
+ | == 1.1.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< | ||
+ | |||
+ | $${LU}_{k, | ||
+ | |||
+ | Where the shares (probabilities) have to add up to one: | ||
+ | |||
+ | $$1 = \sum_{k}^{}{{sh}_{jk}, | ||
+ | |||
+ | The total areas converted from initial land use j into final land use k over the horizon (t-s) are denoted L< | ||
+ | |||
+ | $$\ \left( \lambda_{jk} - 1 \right)L_{jk}^{- 1} - \mu_{jk} + \tau_{k}^{} + \tau_{j}^{initial} = 0$$ | ||
+ | |||
+ | Where λ< | ||
+ | |||
+ | 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. | ||
+ | |||
+ | === 2. 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. | ||
+ | |||
+ | == 2.1. 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< | ||
+ | |||
+ | $$\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 δ< | ||
+ | |||
+ | ==== 2.2. 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< | ||
+ | |||
+ | $${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 (ϕ< | ||
+ | |||
+ | ==== 2.3. 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< | ||
+ | |||
+ | $$G = \Gamma \bullet h = \sum_{i, | ||
+ | |||
+ | 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, | ||
+ | |||
+ | Where Λ< | ||
+ | |||
+ | Using these average annual transitions for true (off-diagonal) LUC we may compute the final classes as follows: | ||
+ | |||
+ | $${LU}_{k, | ||
+ | |||
+ | 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, | ||
+ | |||
+ | ===== 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 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_ | ||
+ | |GWPCO2SOI_ | ||
+ | |GWPN2OSOI_ | ||
+ | |GWPCO2BUR_ | ||
+ | |GWPCH4BUR_ | ||
+ | |GWPN2OBUR_ | ||
+ | |GWPCO2HIS_ | ||
+ | |GWPCH42HIS_ | ||
+ | |||
+ | 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==== | ||
Line 1094: | Line 1564: | ||
The following table shows the substitution elasticities used for the different product groups. Compared to most other studies, we opted for a rather elastic substitution between products from different origins, as agricultural products are generally more uniform then aggregated product groups, as they can be found e.g. in CGE models. | The following table shows the substitution elasticities used for the different product groups. Compared to most other studies, we opted for a rather elastic substitution between products from different origins, as agricultural products are generally more uniform then aggregated product groups, as they can be found e.g. in CGE models. | ||
- | **Table 28: Substitution elasticities for the Armington CES utility aggregators((A sensitivity analysis on those elasticities is given in section [[Sensitivity analysis]]))** | + | **Table 28: Substitution elasticities for the Armington CES utility aggregators((A sensitivity analysis on those elasticities is given in section [[scenario simulation#Sensitivity analysis]]))** |
^Product (group) ^Substitution elasticity between domestic sales and imports | ^Product (group) ^Substitution elasticity between domestic sales and imports | ||
Line 1101: | Line 1571: | ||
|Other fruits | | |Other fruits | | ||
|Sugar| 12 | 12 | | |Sugar| 12 | 12 | | ||
- | |All other products| 8 | | + | |All other products| 8 | |
+ | Source: own calculations | ||
There are some specific settings, such as a value of 2 for rice and the EU15, 2.5 respectively 5 for Japan to account for its specific tariff system, as well as some lower values for EU’s Mediterrean partner countries. | There are some specific settings, such as a value of 2 for rice and the EU15, 2.5 respectively 5 for Japan to account for its specific tariff system, as well as some lower values for EU’s Mediterrean partner countries. | ||
Line 1109: | Line 1579: | ||
{{:: | {{:: | ||
- | The above “primal” formulation of the Armington approach in terms of quantity aggregators turned out numerically less stable in the implementaiotn than the dual representation in terms of price aggregators. The Armington approach suffers from two important shortcomings. First of all, a calibration to a zero flow is impossible so that only observed import flows react to policy changes while all others are fixed at zero level. For most simulation runs, that shortcoming should not be serious. If it is relevant, it may be overcome using the modified Armington approach as explained in Section [[Market module for agricultural outputs#Price linkages]]. | + | The above “primal” formulation of the Armington approach in terms of quantity aggregators turned out numerically less stable in the implementaiotn than the dual representation in terms of price aggregators. The Armington approach suffers from two important shortcomings. First of all, a calibration to a zero flow is impossible so that only observed import flows react to policy changes while all others are fixed at zero level. For most simulation runs, that shortcoming should not be serious. If it is relevant, it may be overcome using the modified Armington approach as explained in Section [[scenario simulation#Price linkages]]. |
Secondly, the Armington aggregator defines a utility aggregate and not a physical quantity. That second problem is healed by re-correcting in the post model part to physical quantities. Little empirical work can be found regarding the estimation of the functional parameters of Armington systems. Hence, substitution elasticities were chosen as to reflect product properties as shown above. | Secondly, the Armington aggregator defines a utility aggregate and not a physical quantity. That second problem is healed by re-correcting in the post model part to physical quantities. Little empirical work can be found regarding the estimation of the functional parameters of Armington systems. Hence, substitution elasticities were chosen as to reflect product properties as shown above. | ||
Line 1234: | Line 1704: | ||
**Figure 20: Witzke et al. calibration, | **Figure 20: Witzke et al. calibration, | ||
- | {{::figure20.png? | + | {{::figure_20.png? |
The additional commitment parameter involves another degree of freedom that needs to be eliminated with additional information. During the calibration this is provided by the expected imports from region 2 at the second hypothetical set of relative prices. Following the dual approach, the lower Armington nest is represented with Armington share-equations and with equations for the composite price indexes: | The additional commitment parameter involves another degree of freedom that needs to be eliminated with additional information. During the calibration this is provided by the expected imports from region 2 at the second hypothetical set of relative prices. Following the dual approach, the lower Armington nest is represented with Armington share-equations and with equations for the composite price indexes: | ||
Line 1250: | Line 1720: | ||
**Figure 21: GUI Option for the non-homothetic Armington system** | **Figure 21: GUI Option for the non-homothetic Armington system** | ||
- | {{::figure21.png?600}} | + | {{::figure_21.png?600}} |
The calibration of the non-homothetic Armington demand system does not require a full re-calibration of the complete CAPRI modelling system; it can be found under the workstep “Run scenario”, | The calibration of the non-homothetic Armington demand system does not require a full re-calibration of the complete CAPRI modelling system; it can be found under the workstep “Run scenario”, | ||
Line 1488: | Line 1958: | ||
\begin{equation} | \begin{equation} | ||
- | intd_{i,r} = (intk_{i, | + | intd_{i,r} = (intk_{i, |
\end{equation} | \end{equation} | ||
scenario_simulation.txt · Last modified: 2023/09/08 12:07 by massfeller