Differences

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

--- disaggregation_of_crop_areas [2020/03/28 08:20] – [Simulation model m_hpdCropSpat] matsz
+++ disaggregation_of_crop_areas [2020/03/28 09:01] – [Data sets] matsz
@@ Line 211: / Line 211: @@
 ====Data sets====
+Update pending
+FSS 2010 data at nested grid levels \\
+FSS 2010 data, gap-filled at 10km-NUTS3 overlay \\
+Forest map
+====Data preparation====
+===Re-mapping from FSS crops to CAPRI crops===
+  # Convert to CAPRI activities (posteact)
+  #As grids and NUTS are not always consistent - use the grid-HSUs to fill with crops
+  p_hsufraction(curgrid,%spatunit%_all)=p_hsu_grid10n23(%spatunit%_all,curgrid,"fracHSU");
+  grid10n23_hsu(curgrid,%spatunit%_all)$p_hsufraction(curgrid,%spatunit%_all)=yes;
+  cur%spatunit%(%spatunit%_all)=yes$sum(grid10n23_hsu(curgrid,%spatunit%_all),1);
+===Intersecting FSS 2010 10km x NUTS3 to FSU===
+The land use distribution model works at the spatial intersection of the prior and posterior spatial units. Historically (when CAPDIS was based on HSU) this intersection was an area determined by the fraction of the HSU that lies within different FSS-admin grid cells. This fraction \(f_{hg}\in[0,1]\).
+However after the update to the FSU, each spatial unit was fully in one FSS-admin grid cell only, thus \(f_{hg}\in\{0,1\}\). The following text is therefore relevant only if CAPDIS is run with ‘old’ HSU.
+  # Work on the intersection between grid and HSU (Line 585 ff)
+  p_levlunit(%region%,cur%spatunit%,"AREA")
+  =p_hsufraction(%region%,cur%spatunit%)*p_hsu(cur%spatunit%,"area");
+  " Work on intersection of HSU and grid for crops"' '" "'
+  # Distribute crops over intersected units –and scale to total area (should be already but not always is...)
+  p_levlunit(%region%,cur%spatunit%,%croptp%)
+  =p_hsufraction(%region%,cur%spatunit%)*p_levlmean(cur%spatunit%,%croptp%);
+  # Scale all areas such that the sum becomes the AREA of the unit (in case the HSU had to be split)
+  #    - if LAPM predictions are consistent with total area
+  p_temp3dim(%region%,cur%spatunit%,"allarea")
+  =sum(%croptp%$(not sameas(%croptp%,"FORE")),p_levlunit(%region%,cur%spatunit%,%croptp%));
+  p_levlunit(%region%,cur%spatunit%,%croptp%)
+	$(p_temp3dim(%region%,cur%spatunit%,"allarea") and not sameas(%croptp%,"FORE"))
+	=p_levlunit(%region%,cur%spatunit%,%croptp%)*(p_levlunit(%region%,cur%spatunit%,"AREA")-
+		p_levlunit(%region%,cur%spatunit%,"FORE"))
+		/p_temp3dim(%region%,cur%spatunit%,"allarea");
+In order to ensure consistency of total area between the prior and the posterior data sets, all data are re-mapped into their intersection. This is achieved with the fraction of the spatial units of one layer that is part of a unit of the second spatial layer:
+\begin{equation}
+a_{hg}=\sum_{h,g}\{a_h\cdot f_{hg} \}
+\end{equation}
+\(a_{hg}\) =	Area [parameter, km2] of unit //u// intersecting spatial unit //h// and grid cell //g//. \\
+\(a_h\)	= Area [parameter, km2] of spatial unit //h// \\
+\(f_{hg}\) =	Fraction of spatial unit// h//, which is covered by grid cell //g// \\
+Cultivated crop areas are re-mapped to the intersecting units proportionally to the area fraction, assuming homogeneous distribution of each crop within each spatial unit //h//.
+The LAPM predictions are not constrained to exhaust the total available area. However, this is a required characteristic in CAPRI. Therefore, the land use areas (crops, forest land and ‘other’ area) are scaled so that their sum matches the total available area. As forest areas are obtained from a data set, which is assumed to be of high precision, forest areas are excluded from scaling.
+\begin{equation}
+a_{c^o,hg}=\sum_{h,g} \left\{ a_{c,h} \cdot f_{hg} \cdot \frac{a_{hg} - a_{forest,hg}}{\sum_{c^{o^\prime}} a_{c^{o^\prime}}} \right\},
+\end{equation}
+\(c^*\)	Land use. c^*∈{c,other land} \\
+\(a_{c^o,hg}\)		Area [parameter, km<sup>2</sup>] cultivated with crop //c// or covered by ‘other land’ use excluding forest in unit //u// intersecting spatial unit //h// and grid cell //g//.  \\
+\(a_{c^o,h}\)		Area [parameter, km<sup>2</sup>] cultivated with crop //c// or covered by ‘other land’ use excluding forest in spatial unit //h// \\
+\(f_{hg}\)		Fraction of spatial unit //h//, which is covered by grid cell //g// \\