climate_indices icon indicating copy to clipboard operation
climate_indices copied to clipboard

Add a PET function using Hargreaves' equation

Open monocongo opened this issue 6 years ago • 5 comments

Begueria et al (2014) as well as Allen et al (1998) suggest Hargreaves as a better PET approximation method compared to Thornthwaite's when solar radiation, wind speed, and relative humidity data are unavailable.

The Hargreaves equation is spelled out in Droogers and Allen (2002), Allen et al (1998) and Hargreaves (1994).

There is R code for this, look into that and see what can be re-used.

Thanks/credit to Rebecca Ward for suggesting this.

Related to #101

monocongo avatar Apr 11 '18 20:04 monocongo

Another R implementation

monocongo avatar Apr 13 '18 20:04 monocongo

My colleague @erin-telfer has a notebook that takes netcdf as input (converts to xarray) to calculate PET using Hargreaves, this is how I calculated PET to use with the climate_indices module.

Could be a bit more helpful than the above R implementations. https://github.com/GeoscienceAustralia/dea-notebooks/blob/ErinT/ewater/NB07_calculate_PET_grids.ipynb

s-m-t-c avatar Jul 10 '18 04:07 s-m-t-c

Thanks for your help on this. Unfortunately, I can't reach that URL. ANy chance you can copy/paste the code someplace else where I can get to it? I'm not actively maintaining this code lately but I'd like to get to this eventually, and if this is as you describe then it should be relatively painless.

monocongo avatar Jul 10 '18 16:07 monocongo

This implementation looks solid: https://pyeto.readthedocs.io/en/latest/api.html#pyeto.hargreaves

monocongo avatar Sep 02 '20 17:09 monocongo

My colleague @erin-telfer has a notebook that takes netcdf as input (converts to xarray) to calculate PET using Hargreaves, this is how I calculated PET to use with the climate_indices module.

Could be a bit more helpful than the above R implementations. https://github.com/GeoscienceAustralia/dea-notebooks/blob/ErinT/ewater/NB07_calculate_PET_grids.ipynb

Could you please share another link to the notebook. Unable to reach the URL

SouhailAB avatar Apr 20 '21 11:04 SouhailAB