GitLab

Created an R script for manipulating the GHG EFs and to adjust them to the new zones

Created a new folder with R script to calculate NP contents in R for forage crops from Nutrient tool database https://plantsorig.sc.egov.usda.gov/npk/main

external function script

Added the suitability crops to the excel > use the distinct crop list -> the other is not updated with suitability crops.

Creating a dataset with rangeland, arable land, arable pasture, natural pasture per zone (adm0, aez, soilagri). draft, not finished.

Final dataset for the globe and for EU28. Also cifos crop names are added. Issues are the not available disaggregation and the impication on nutrient intake (Rest, tnut, Veggie). Which nutrient content should be taken as a proxy?

all crops done. treenuts, veggies, EU subsetting. Function R script created. Suitability file created and transformed, linked to zones AEZ ADM0 Soilagri.

New zonation based on ADM0 - AEZ - DOMSOIL_IPCC(mineral and Organic) - Validation with EUROSTAT and FAO. Data transformation. Different types of agricultural land extractions

Function, defining the different surogate crop species for aggregated groups.

Peat done (except of merging the yield and area of grass >> still some issues). Single files are fine for crop and grassland area now with peat and nonpeat.
SPAM extraction from GEOTIFF file >> New zonal statistics function > much more efficient now. Replaced left_join with bind function. Now much less memory demanding.
Extracted factors from columnames.

Extraction of SPAM geotiff files based on ULrikes recommendation. Then I can use GAUL as a consistent dataset.

Deleted all old carbage from the script.

Started with the creation of the EU crop area validation data >> still needs finishing (see issues on git)

> No NA's > Only the **ratio of harvested area** for disaggregation. Like this, there is more consistency with the presence of nuts/veggies in adm1 zones. Otherwise, it happened that there is area but no production vice versa. Then you divide by 0 = inf.
> No inf values: The inf values came from the division of P/H=Y calculation. The P was in all cases very low so I decided to simply recode them to 0. Then Yields would be also 0.
> rounding fixed > 1 decimal like in the original dataset so we do not lose the level of detail.
> Yield conversion (production in mt) checked and correct.
> Aggregation to adm1 checked and correct.

> No NA's > Only the **ratio of harvested area** for disaggregation. Like this, there is more consistency with the presence of nuts/veggies in adm1 zones. Otherwise, it happened that there is area but no production vice versa. Then you divide by 0 = inf.
> No inf values: The inf values came from the division of P/H=Y calculation. The P was in all cases very low so I decided to simply recode them to 0. Then Yields would be also 0.
> rounding fixed > 1 decimal like in the original dataset so we do not lose the level of detail.
> Yield conversion (production in mt) checked and correct.
> Aggregation to adm1 checked and correct.

I fixed some computation problems. When I directly combined the SPAM datasets before running all the aggregation functions on them, my computer with 100 GB ram could not process it. So I changed it back to the old version with single files and then at last combine the datasets. Now it worked. Issue: No difference when disaggregating with harv and prod propotions on final values. Check why!

Veggies extraction from monfreda are now updated and code is much more concise with regex. Worked up to the adm1 aggregation of SPAM crop data line 652.

updated and validated the treenut script. Also the R.script now contains all the script needed to run the adm1 aggregation, treenut extraction. Before this was seperated in the treenuts script.