kunstler/climatenfi

climatenfi's Issues

Compare PET

TerraClimate provide PET at 4 x 4 km based on Penman Montieth
We can compute PET with Thornwaite based only on T and Lat
or with hargraeaves with tmin tmax and radiation (from WC mean over the period at 1x1km (kJ m-2 day-1) or time varying srad at 2 x 2 km from TerraClimate srad Downward surface shortwave radiation), units = W/m2).

We can compte these PET with this package https://rdrr.io/cran/SPEI/man/PET.html

Sum of growing degree days from monthly mean

In CHELSA
We interpolated mean monthly temperatures
(t) to a daily resolution using a B-cubic spline interpolation and used this as input data for the
growing degree days.

See also review of method by Julien

ERA5 land

Should we explore ERA5 - land as well (11 x 11 km but dowscaling possible )

https://www.erikkusch.com/project/krigr/

Unit Terraclimate

aet (Actual Evapotranspiration, monthly total), units = mm
def (Climate Water Deficit, monthly total), units = mm
pet (Potential evapotranspiration, monthly total), units = mm
ppt (Precipitation, monthly total), units = mm
q (Runoff, monthly total), units = mm
soil (Soil Moisture, total column - at end of month), units = mm
srad (Downward surface shortwave radiation), units = W/m2
swe (Snow water equivalent - at end of month), units = mm
tmax (Max Temperature, average for month), units = C
tmin (Min Temperature, average for month), units = C
vap (Vapor pressure, average for month), units = kPa
ws (Wind speed, average for month), units = m/s
vpd (Vapor Pressure Deficit, average for month), units = kpa
PDSI (Palmer Drought Severity Index, at end of month), units = unitless

Compute minimum soil water potential Soil PSI min

Benito Garzon et al. 2017 GEB computed Soil PSI min from

• ERA-INTERIM volumetric soil water content (spatial resolution (10 x 15 km)
  https://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/

• Soil water texture data from JRC Europe
  http://eusoils.jrc.ec.europa.eu/wrb/

https://esdac.jrc.ec.europa.eu/content/european-soil-database-v2-raster-library-1kmx1km
https://esdac.jrc.ec.europa.eu/content/sgdbe-attributes#TEXT_SRF_DOM dominant surface

The for each plot the max volumetric soil water content Wmax was computed as mean of monthly data for 5 wetter year and minimum for 5 driest year Wmin. For each soil category computed water potential at maximum hydratation PSImax

Then PSImin is computed PSImin = PSImax * (Wmin/Wmax)^b with b computed per soil type.

Extrac soil data

todo https://www.isric.org/explore/soilgrids

UPDATE download command

wget -P data/ --no-host-directories --force-directories --input-file=envidatS3paths_all.txt