i.evapo.senay
Actual evapotranspiration, method after Senay (2007)
i.evapo.senay [-edb] temperature=name albedo=name elevation=name [latitude=name] [dayofyear=name] [transmissivitysingleway=name] [waterdensity=float] [slope=name] [aspect=name] [atmosphericemissivity=float] [airtemperature=name] [surfaceemissivity=name] ndvi=name [diurnaletpotential=name] output=name [evapfr=name] [--overwrite] [--verbose] [--quiet] [--qq] [--ui]
Example:
i.evapo.senay temperature=name albedo=name elevation=name ndvi=name output=name
grass.script.run_command("i.evapo.senay", temperature, albedo, elevation, latitude=None, dayofyear=None, transmissivitysingleway=None, waterdensity=1005.0, slope=None, aspect=None, atmosphericemissivity=None, airtemperature=None, surfaceemissivity=None, ndvi, diurnaletpotential=None, output, evapfr=None, flags=None, overwrite=False, verbose=False, quiet=False, superquiet=False)
Example:
gs.run_command("i.evapo.senay", temperature="name", albedo="name", elevation="name", ndvi="name", output="name")
Parameters
temperature=name [required]
Name of the temperature layer [Degree Kelvin]
albedo=name [required]
Name of the Albedo layer [0.0-1.0]
elevation=name [required]
Name of the elevation layer [m]
latitude=name
Name of the degree latitude layer [dd.ddd]
dayofyear=name
Name of the Day of Year layer [0.0-366.0]
transmissivitysingleway=name
Name of the single-way transmissivity layer [0.0-1.0]
waterdensity=float
Value of the density of fresh water ~[1000-1020]
Default: 1005.0
slope=name
Name of the Slope layer ~[0-90]
aspect=name
Name of the Aspect layer ~[0-360]
atmosphericemissivity=float
Value of the apparent atmospheric emissivity (Bandara, 1998 used 0.845 for Sri Lanka)
airtemperature=name
Name of the Air Temperature layer [Kelvin], use with -b
surfaceemissivity=name
Name of the Surface Emissivity layer [-], use with -b
ndvi=name [required]
Name of the NDVI layer [-]
diurnaletpotential=name
Name of the ET Potential layer [mm/day]
output=name [required]
Name of the output Actual ET layer
evapfr=name
Name of the output evaporative fraction layer
-e
ET Potential Map as input (By-Pass creation of one)
-d
Slope/Aspect correction
-b
Net Radiation Bandara (1998), generic Longwave calculation, need apparent atmospheric emissivity, Air temperature and surface emissivity inputs
--overwrite
Allow output files to overwrite existing files
--help
Print usage summary
--verbose
Verbose module output
--quiet
Quiet module output
--qq
Very quiet module output
--ui
Force launching GUI dialog
temperature : str, required
Name of the temperature layer [Degree Kelvin]
Used as: input, raster, name
albedo : str, required
Name of the Albedo layer [0.0-1.0]
Used as: input, raster, name
elevation : str, required
Name of the elevation layer [m]
Used as: input, raster, name
latitude : str, optional
Name of the degree latitude layer [dd.ddd]
Used as: input, raster, name
dayofyear : str, optional
Name of the Day of Year layer [0.0-366.0]
Used as: input, raster, name
transmissivitysingleway : str, optional
Name of the single-way transmissivity layer [0.0-1.0]
Used as: input, raster, name
waterdensity : float, optional
Value of the density of fresh water ~[1000-1020]
Default: 1005.0
slope : str, optional
Name of the Slope layer ~[0-90]
Used as: input, raster, name
aspect : str, optional
Name of the Aspect layer ~[0-360]
Used as: input, raster, name
atmosphericemissivity : float, optional
Value of the apparent atmospheric emissivity (Bandara, 1998 used 0.845 for Sri Lanka)
airtemperature : str, optional
Name of the Air Temperature layer [Kelvin], use with -b
Used as: input, raster, name
surfaceemissivity : str, optional
Name of the Surface Emissivity layer [-], use with -b
Used as: input, raster, name
ndvi : str, required
Name of the NDVI layer [-]
Used as: input, raster, name
diurnaletpotential : str, optional
Name of the ET Potential layer [mm/day]
Used as: input, raster, name
output : str, required
Name of the output Actual ET layer
Used as: output, raster, name
evapfr : str, optional
Name of the output evaporative fraction layer
Used as: output, raster, name
flags : str, optional
Allowed values: e, d, b
e
ET Potential Map as input (By-Pass creation of one)
d
Slope/Aspect correction
b
Net Radiation Bandara (1998), generic Longwave calculation, need apparent atmospheric emissivity, Air temperature and surface emissivity inputs
overwrite: bool, optional
Allow output files to overwrite existing files
Default: False
verbose: bool, optional
Verbose module output
Default: False
quiet: bool, optional
Quiet module output
Default: False
superquiet: bool, optional
Very quiet module output
Default: False
DESCRIPTION
i.evapo.senay Calculates the diurnal actual evapotranspiration after Senay (2007). This is converting all Net radiation from the diurnal period into ET, then uses Senay equation for evaporative fraction.
It takes input maps of Albedo, surface skin temperature, latitude, day of year, single-way transmissivity and takes input value of the density of fresh water.
DEM is used for calculating min and max temperature for Senay equation.
The "-s" flag permits output map of evaporative fraction from Senay.
NOTES
If you are trying to map irrigated crops, and you know there is at least one crop pixel that is evapotranspiring at maximum (ETa=ETpot), then read this.
i.evapo.senay is highly dependent on the wet pixel being the lowest temperature in the crop pixels to work for non water stressed crops, force it that way, even if it breaks non crop areas. I suggest you reduce your region to the irrigation system boundaries, checking that it includes a bit of dry area for the hot/dry pixel.
Since it is a direct relationship to LST, evaporative fraction can be very sensitive to the kind of pixel sample you feed it with.
TODO
SEE ALSO
r.sun, i.albedo, i.eb.eta, i.eb.evapfr, i.evapo.potrad
REFERENCES
- Senay 2007
Chemin, Y., 2012. A Distributed Benchmarking Framework for Actual ET Models, in: Irmak, A. (Ed.), Evapotranspiration - Remote Sensing and Modeling. InTech. (PDF, DOI: 10.5772/23571)
AUTHOR
Yann Chemin, International Rice Research Institute, The Philippines.
SOURCE CODE
Available at: i.evapo.senay source code
(history)
Latest change: Friday Feb 21 10:10:05 2025 in commit 7d78fe3