Feature/simple crop model

[Bachibouzouk]

Summary of the discussion

The model runs through

Open questions:

• how many days should we simulate to get a valid answer from crop production ? --> whole year
• shouldn't we need to incorporate f(solar) ? --> yes
• is PAW equal to vwc in equation 12? --> not 100% sure

Type of change (CHANGELOG.md)

Added

• Facade SimpleCrop in wefe subpackage to simulate a crop growth without accounting for water in (following this reference https://doi.org/10.1016/j.eja.2019.01.009), a test model is provided in examples/scenarios/test_crop.

Changed

• Now validate_datapackage of facades must return the resource (modified or not)

[JulianBarinton]

Comments on open questions:

how many days should we simulate to get a valid answer from crop production ? -> we should at least simulate for a whole crop growth period... to be sure we can already simulate for the whole year.

shouldn't we need to incorporate f(solar) ? -> Yes let's include f(solar); @MaaJoo13 has created suitable code for that

is PAW equal to vwc in equation 12? -> see specific comment above regarding PAW

[JulianBarinton]

Open Question #2: Why is banana modelled as excess and has the marginal cost of -0.1? Which evapotranspiration time series to use for Et_0 (reference evapotransporation)? era5 provides actual evapotranspiration and potential evapotranspiration?

[MaaJoo13]

Hey you,

here are my thoughts on this:

1) f_solar/crop growth cycle

f_solar represents the crop growth cycle. In the APV input excel sheet (via mail before my holidays) I tried to explain how I incorporated that using sowing_date and harvest_date. That can be used for the crop model as well.

2) water stress factor

First, I did not use ERA5 evapotranspiration or vwc.

According to ARID, plant available water PAW is vwc rzd. awc is the average water holding capacity, it is basically an upper limit to vwc. So awc vwc does not really make sense.

So what we need for the water stress factor is et_o and PAW with PAW being computed from vwc and rzd.

Rzd is taken from SIMPLE, they have a table with data for every crop.

ET_o is calculated using the penman-monteith approach, making use of ghi, air temp, dew point temp and wind speed (all 4 taken from era5).

For vwc I then created a water balance according to ARID.

In general, I think we should have a quick meeting on how to handle the crop model.

[Bachibouzouk]

@MaaJoo13 TODOS:

• implement fSolar
• check vwc preprocessing (only if dew water point is available)
• run the model for several days to get the full plant life and check the production matches the data from https://doi.org/10.1016/j.eja.2019.01.009
• soil_dict parameters are not available for all locations anyway

Next step solar and water input, biomass production as output