MAgPIE - An Open Source land-use modeling framework

4.3.5

created with goxygen 1.3.0

Yields (14_yields)

Description

The yields module simulates the temporal development of crop yields and pasture productivity. Spatially explicit information on pasture productivity and crop yields under both rainfed and irrigated conditions is provided by the global gridded crop model LPJmL (Lund-Potsdam-Jena with managed Land) (Bondeau et al. 2007). In the initial year of the simulation period, crop yields and pasture productivity are calibrated at the regional level to meet the observed cropland and pasture area as reported by FAO (FAOSTAT 2016). For the simulation of the temporal development of agricultural yields, the module receives information about the agricultural land use intensity represented by the \(\tau\) factor coming from the module 13_tc.

The module returns yields for all crops and for pasture, which is then used by the modules 30_crop and 31_past.

Interfaces

Input

Output

Realizations

(A) biocorrect

The biocorrect realization reads in the LPJmL data and performs several corrections. First, a bioenergy yield correction is performed. As there is currently no robust information on bioenergy yields available in (FAOSTAT 2016), it is assumed that the LPJmL yields for bioenergy correspond to the yields achieved under the highest currently observed value of the \(\tau\) factor representing agricultural land use intensity. Bioenergy yields are downscaled proportionally to the respective \(\tau\) factor in the given region. Second, yields for all other crops are corrected on the regional level by applying a calibration factor that does not differentiate between crops. Pasture yields have their own regional calibration factors. Applied calibration factors are derived via comparing cropland and pasture areas in the initial time step with values reported by FAO (FAOSTAT 2016).

\[\begin{multline*} vm\_yld(j2,kcr,w) = \sum_{ct}i14\_yields(ct,j2,kcr,w) \cdot \sum_{cell(i2,j2), supreg(h2,i2)}\left(\frac{vm\_tau(h2)}{fm\_tau1995(h2)}\right) \end{multline*}\]

For the current time step of the optimization, cellular yields of irrigated and rainfed crops are calculated by multiplying calibrated input yields from LPJmL with the intensification rate relative to the initial time step 1995.

\[\begin{multline*} vm\_yld(j2,"pasture",w) = \sum_{ct}i14\_yields(ct,j2,"pasture",w) \cdot \left(1 + s14\_yld\_past\_switch \cdot \left(\sum_{cell(i2,j2), supreg(h2,i2)}\left(\frac{vm\_tau(h2)}{fm\_tau1995(h2)}\right) - 1\right)\right) \end{multline*}\]

In the case of pasture yields, technological change cannot be fully translated into yield increases. To account for that, the parameter s14_yld_past_switch is defined to capture a certain magnitude of spillovers that can range from 0 (no spillover) to 1 (full spillover).

Limitations The magnitude of spillover effects from technological change in the crop sector towards improvements in pasture management is very uncertain.

(B) dynamic_aug18

In the dynamic_aug18 realization, the crop yield calculations are identical as in the above described realization (biocorrect). Additionally, this realization also calculates the growing stocks in commercial plantations and natural vegetation using LPJmL Carbon stocks.

\[\begin{multline*} vm\_yld(j2,kcr,w) = \sum_{ct}i14\_yields(ct,j2,kcr,w) \cdot \sum_{cell(i2,j2), supreg(h2,i2)}\left(\frac{vm\_tau(h2)}{fm\_tau1995(h2)}\right) \end{multline*}\]

Pasture yields are not linked to yield increases in the crop sector, but to an exogenous pasture management factor pm_past_mngmnt_factor:

\[\begin{multline*} vm\_yld(j2,"pasture",w) = \sum_{ct}\left(i14\_yields(ct,j2,"pasture",w) \cdot \sum\left(cell(i2,j2),pm\_past\_mngmnt\_factor(ct,i2)\right)\right) \end{multline*}\]

carbon density. We convert Carbon density in tC/ha to tDM/ha by using carbon fraction of s14_carbon_fraction in tC/tDM. For assessing wood harvesting we need only above ground biomass information, therefore we multiply with aboveground f14_aboveground_fraction. Additionally, we divide aboveground tree biomass by biomass conversion and expansion (BCE) factor to get stem biomass in tDM/ha.

Limitations The exogenous implementation of pasture intensification cannot capture feedbacks between land scarcity and efforts to improve pasture management.

(C) managementcalib_aug19

The managementcalib_aug19 realization reads in the LPJmL data and performs several corrections. First, a bioenergy yield correction is performed. As there is currently no robust information on bioenergy yields available in (FAOSTAT 2016), it is assumed that the LPJmL yields for bioenergy correspond to the yields achieved under the highest currently observed value of the \(\tau\) factor representing agricultural land-use intensity. Pasture yields are calculated based on pasture demand to account for in- and extensification of managed grasslands. Technological change spillover from crop yield increases of the time step before is included additionally. Second, yields for all other crops are corrected on the regional level by applying a calibration factor that does not differentiate between crops. Pasture yields have their own regional calibration factors. Applied calibration factors are derived via comparing cropland and pasture areas in the initial time step with values reported by FAO (FAOSTAT 2016). This realization also calculates the growing stocks in commercial plantations and natural vegetation using LPJmL Carbon stocks.

Technological change can increase the initial calibrated yields by:

\[\begin{multline*} vm\_yld(j2,kcr,w) = \sum_{ct}i14\_yields\_calib(ct,j2,kcr,w) \cdot \sum_{cell(i2,j2), supreg(h2,i2)}\left(\frac{ vm\_tau(h2) }{ fm\_tau1995(h2)}\right) \end{multline*}\]

For the current time step of the optimization, cellular yields of irrigated and rainfed crops are calculated by multiplying calibrated input yields from LPJmL with the intensification rate relative to the initial time step 1995. In the case of pasture yields, technological change cannot be fully translated into yield increases, to address that, an exogenous pasture management factor pm_past_mngmnt_factor is used to scale pasture yields based on the number of cattle reared to fulfill the domestic demand for ruminant livestock products in module 70.

Additionally, the parameter s14_yld_past_switch can be used to capture a certain magnitude of spillovers of the yield increase due to technological change from the time step before. It can range from 0 (no spillover) to 1 (full spillover).

\[\begin{multline*} vm\_yld(j2,"pasture",w) = \sum_{ct}\left(i14\_yields\_calib(ct,j2,"pasture",w) \cdot \sum_{cell(i2,j2)}pm\_past\_mngmnt\_factor(ct,i2)\right) \cdot \left(1 + s14\_yld\_past\_switch \cdot \left(\sum_{cell(i2,j2), supreg(h2,i2)}\left(\frac{ pcm\_tau(h2)}{fm\_tau1995(h2)}\right) - 1\right)\right) \end{multline*}\]

The following equations calibrate the cellular yield patterns (‘f14_yields’) to match FAO historical yields (‘f14_fao_yields_hist’) by calculating a calibration term called ‘i14_managementcalib’. For most cases, ‘i14_managementcalib’ is the ratio of the historical yields reported by FAO (‘f14_fao_yields_hist’) and regional mean yields (‘i14_modeled_yields_hist’) given historic crop area patterns (‘fm_croparea’) and cellular yields coming from crop models like LPJmL (‘f14_yields’). In these cases, ‘i14_managementcalib’ represents a purely relative calibration factor that depends only on the initial conditions of the starting year.

However, when FAO yields are significantly higher than given by the cellular yield inputs (underestimated baseline), the relative calibration terms can lead to unrealistically large yields in the case of future yield increases within the cellular yield patterns.

To address this issue, the factor ‘i14_lambda_yields’ determines the degree to which the baseline (FAO) is under- or overestimated and therefore controls whether the calibration factor is applied as an absolute or relative change. For overestimated FAO yields, ‘i14_lambda_yields’ is 1, which is equivalent to an entirely relative calibration. For underestimated yields, ‘i14_lambda_yields’ is calculated as the squared root of the ratio between LPJmL yields and FAO historical yields, and as ‘i14_lambda_yields’ approaches 0, it reduces the applied relative change resulting in a mean change increasingly similar to an additive term (Heinke et al. (2013)). This concept is referred to as limited calibration, as it limits the calibration to an additive term in case of a strongly underestimated baseline. The scalar

i14_croparea_total(t_all,w,j) = sum(kcr, fm_croparea(t_all,j,w,kcr));

Historic crop area patterns (‘fm_croprea’) are used to calculate regional yields (‘i14_modeled_yields_hist’) from the given cellular input pattern. In rare cases where a region has no crop area reported for a given crop type, the total crop area is used to calculate a proxy yield for the calibration, given by the following equation:

i14_modeled_yields_hist(t_past,i,knbe14)
   = (sum((cell(i,j),w), fm_croparea(t_past,j,w,knbe14) * f14_yields(t_past,j,knbe14,w)) /
      sum((cell(i,j),w), fm_croparea(t_past,j,w,knbe14)))$(sum((cell(i,j),w), fm_croparea(t_past,j,w,knbe14))>0)
   + (sum((cell(i,j),w), i14_croparea_total(t_past,w,j) * f14_yields(t_past,j,knbe14,w)) /
      sum((cell(i,j),w), i14_croparea_total(t_past,w,j)))$(sum((cell(i,j),w), fm_croparea(t_past,j,w,knbe14))=0);

The factor ‘i14_lambda_yields’ is calculated for the initial time step depending on the setting ‘s14_limit_calib’ and is then held constant for all other time steps. The regional FAO yield and regional yield of the crop model input of the initial time step is kept constant in the two parameters ‘i14_fao_yields_hist’ and ‘i14_modeled_yields_hist’:

loop(t,
     if(sum(sameas(t,"y1995"),1)=1,
          if    ((s14_limit_calib = 0),
               i14_lambda_yields(t,i,knbe14) = 1;
          Elseif (s14_limit_calib =1 ),
               i14_lambda_yields(t,i,knbe14) =
                    1$(f14_fao_yields_hist(t,i,knbe14) <= i14_modeled_yields_hist(t,i,knbe14))
                    + sqrt(i14_modeled_yields_hist(t,i,knbe14)/f14_fao_yields_hist(t,i,knbe14))$
                    (f14_fao_yields_hist(t,i,knbe14) > i14_modeled_yields_hist(t,i,knbe14));
          );
          i14_fao_yields_hist(t,i,knbe14) = f14_fao_yields_hist(t,i,knbe14);
     Else
          i14_modeled_yields_hist(t,i,knbe14) = i14_modeled_yields_hist(t-1,i,knbe14);
          i14_FAO_yields_hist(t,i,knbe14)  = i14_fao_yields_hist(t-1,i,knbe14);
          i14_lambda_yields(t,i,knbe14)   = i14_lambda_yields(t-1,i,knbe14);
     );
);

The calibrated cellular yield ‘i14_yields_calib’ is calculated for each time step depending on the constant values ‘i14_modeled_yields_hist’, ‘i14_fao_yields_hist’, ‘i14_lambda_yields’ and the uncalibrated, cellular yield ‘f14_yields’ following the idea of eq. (9) in Heinke et al. (2013):

i14_managementcalib(t,j,knbe14,w) =
  1 + (sum(cell(i,j), i14_fao_yields_hist(t,i,knbe14) - i14_modeled_yields_hist(t,i,knbe14)) /
                             f14_yields(t,j,knbe14,w) *
      (f14_yields(t,j,knbe14,w) / (sum(cell(i,j),i14_modeled_yields_hist(t,i,knbe14))+10**(-8))) **
                             sum(cell(i,j),i14_lambda_yields(t,i,knbe14)))$(f14_yields(t,j,knbe14,w)>0);
i14_yields_calib(t,j,knbe14,w)    = i14_managementcalib(t,j,knbe14,w) * f14_yields(t,j,knbe14,w);

Note that the calculation is split into two parts for better readability.

Calibrated yields are additionally adjusted by calibration factors ‘f14_yld_calib’ determined in a calibration run. As MAgPIE optimizes yield patterns and FAO regional yields are outlier corrected, historical production and croparea can only be reproduced with this additional step of correction:

i14_yields_calib(t,j,kcr,w)       = i14_yields_calib(t,j,kcr,w)      *sum(cell(i,j),f14_yld_calib(i,"crop"));
i14_yields_calib(t,j,"pasture",w) = i14_yields_calib(t,j,"pasture",w)*sum(cell(i,j),f14_yld_calib(i,"past"));

carbon density. We convert Carbon density in tC/ha to tDM/ha by using carbon fraction of s14_carbon_fraction in tC/tDM. For assessing wood harvesting we need only aboveground biomass information, therefore we multiply with aboveground f14_aboveground_fraction. Additionally, we divide aboveground tree biomass by biomass conversion and expansion (BCE) factor to get stem biomass in tDM/ha.

Limitations The exogenous implementation of pasture intensification cannot capture feedbacks between land scarcity and efforts to improve pasture management. And, the magnitude of spillover effects from technological change in the crop sector towards improvements in pasture management is very uncertain.

Definitions

Objects

Sets

Authors

Jan Philipp Dietrich, Isabelle Weindl, Florian Humpenöder, Anne Biewald, Kristine Karstens

See Also

09_drivers, 10_land, 13_tc, 30_crop, 31_past, 32_forestry, 35_natveg, 45_climate, 52_carbon, 70_livestock

References