PotentialEvapotranspiration: Compute potential evapotranspiration#
- class PotentialEvapotranspiration(*args, **kwds)[source]#
Bases:
Component
Potential Evapotranspiration Component calculates spatially distributed potential evapotranspiration based on input radiation factor (spatial distribution of incoming radiation) using chosen method such as constant or Priestley Taylor. Ref: Xiaochi et. al. 2013 for ‘Cosine’ method and ASCE-EWRI Task Committee Report Jan 2005 for ‘PriestleyTaylor’ method. Note: Calling ‘PriestleyTaylor’ method would generate/overwrite shortwave & longwave radiation fields.
Code author: Sai Nudurupati and Erkan Istanbulluoglu
Examples
>>> from landlab import RasterModelGrid >>> from landlab.components.pet import PotentialEvapotranspiration
>>> grid = RasterModelGrid((5, 4), xy_spacing=(0.2, 0.2)) >>> grid['cell']['radiation__ratio_to_flat_surface'] = np.array([ ... 0.38488566, 0.38488566, ... 0.33309785, 0.33309785, ... 0.37381705, 0.37381705]) >>> PET = PotentialEvapotranspiration(grid) >>> PET.name 'PotentialEvapotranspiration' >>> PET.input_var_names ('radiation__ratio_to_flat_surface',) >>> sorted(PET.output_var_names) ['radiation__incoming_shortwave_flux', 'radiation__net_flux', 'radiation__net_longwave_flux', 'radiation__net_shortwave_flux', 'surface__potential_evapotranspiration_rate'] >>> sorted(PET.units) [('radiation__incoming_shortwave_flux', 'W/m^2'), ('radiation__net_flux', 'W/m^2'), ('radiation__net_longwave_flux', 'W/m^2'), ('radiation__net_shortwave_flux', 'W/m^2'), ('radiation__ratio_to_flat_surface', 'None'), ('surface__potential_evapotranspiration_rate', 'mm')] >>> PET.grid.number_of_cell_rows 3 >>> PET.grid.number_of_cell_columns 2 >>> PET.grid is grid True >>> pet_rate = grid.at_cell['surface__potential_evapotranspiration_rate'] >>> np.allclose(pet_rate, 0.) True >>> PET.current_time = 0.5 >>> PET.update() >>> np.allclose(pet_rate, 0.) False
References
Required Software Citation(s) Specific to this Component
None Listed
Additional References
ASCE-EWRI: The ASCE standardized reference evapotranspiration equation, in: Standardization of Reference Evapotranspiration Task Committee Final Report, edited by: Allen, R. G., Walter, I. A., Elliot, R. L., Howell, T. A., Itenfisu, D., Jensen, M. E., and Snyder, R. L., Technical Committee report to the Environmental and Water Resources Institute of the American Society of Civil Engineers from the Task Committee on Standardization of Reference Evapotranspiration, Reston, VA, USA, 2005.
Zhou, X., Istanbulluoglu, E., and Vivoni, E. R.: Modeling the ecohydrological role of aspect-controlled radiation on tree-grass-shrub coexistence in a semiarid climate, Water Resour. Res., 49, 2872– 2895, doi:10.1002/wrcr.20259, 2013.
- Parameters:
grid (RasterModelGrid) – A grid.
method (required for 'MeasuredRadiationPT') – Priestley Taylor method will spit out radiation outputs too.
priestley_taylor_constant (float, optional) – Alpha used in Priestley Taylor method.
albedo (float, optional) – Albedo.
latent_heat_of_vaporization (float, optional) – Latent heat of vaporization for water Pwhv (Wd/(m*mm^2)).
psychometric_const (float, optional) – Psychometric constant (kPa (deg C)^-1).
stefan_boltzmann_const (float, optional) – Stefan Boltzmann’s constant (W/(m^2K^-4)).
solar_const (float, optional) – Solar constant (W/m^2).
latitude (float, optional) – Latitude (radians).
elevation_of_measurement (float, optional) – Elevation at which measurement was taken (m).
adjustment_coeff (float, optional) – adjustment coeff to predict Rs from air temperature (deg C)^-0.5.
lt (float, optional) – lag between peak TmaxF and solar forcing (days).
nd (float, optional) – Number of days in year (days).
MeanTmaxF (float, optional) – Mean annual rate of TmaxF (mm/d).
delta_d (float, optional) – Calibrated difference between max & min daily TmaxF (mm/d).
current_time (float, required only for 'Cosine' method) – Current time (Years)
const_potential_evapotranspiration (float, optional for) – ‘Constant’ method Constant PET value to be spatially distributed.
Tmin (float, required for 'Priestley Taylor' method) – Minimum temperature of the day (deg C)
Tmax (float, required for 'Priestley Taylor' method) – Maximum temperature of the day (deg C)
Tavg (float, required for 'Priestley Taylor' and 'MeasuredRadiationPT') – methods Average temperature of the day (deg C)
float (obs_radiation) – Observed radiation (W/m^2)
method – Observed radiation (W/m^2)
- property Tavg#
Average temperature of the day (deg C)
Tavg: float, required for ‘Priestley Taylor’ and ‘MeasuredRadiationPT’ methods.
- property Tmax#
Maximum temperature of the day (deg C)
Tmax: float, required for ‘Priestley Taylor’ method.
- property Tmin#
Minimum temperature of the day (deg C)
Tmin: float, required for ‘Priestley Taylor’ method.
- __init__(grid, method='Cosine', priestley_taylor_const=1.26, albedo=0.6, latent_heat_of_vaporization=28.34, psychometric_const=0.066, stefan_boltzmann_const=5.67e-08, solar_const=1366.67, latitude=34.0, elevation_of_measurement=300, adjustment_coeff=0.18, lt=0.0, nd=365.0, MeanTmaxF=12.0, delta_d=5.0, current_time=None, const_potential_evapotranspiration=12.0, Tmin=0.0, Tmax=1.0, Tavg=0.5, obs_radiation=350.0)[source]#
- Parameters:
grid (RasterModelGrid) – A grid.
method (required for 'MeasuredRadiationPT') – Priestley Taylor method will spit out radiation outputs too.
priestley_taylor_constant (float, optional) – Alpha used in Priestley Taylor method.
albedo (float, optional) – Albedo.
latent_heat_of_vaporization (float, optional) – Latent heat of vaporization for water Pwhv (Wd/(m*mm^2)).
psychometric_const (float, optional) – Psychometric constant (kPa (deg C)^-1).
stefan_boltzmann_const (float, optional) – Stefan Boltzmann’s constant (W/(m^2K^-4)).
solar_const (float, optional) – Solar constant (W/m^2).
latitude (float, optional) – Latitude (radians).
elevation_of_measurement (float, optional) – Elevation at which measurement was taken (m).
adjustment_coeff (float, optional) – adjustment coeff to predict Rs from air temperature (deg C)^-0.5.
lt (float, optional) – lag between peak TmaxF and solar forcing (days).
nd (float, optional) – Number of days in year (days).
MeanTmaxF (float, optional) – Mean annual rate of TmaxF (mm/d).
delta_d (float, optional) – Calibrated difference between max & min daily TmaxF (mm/d).
current_time (float, required only for 'Cosine' method) – Current time (Years)
const_potential_evapotranspiration (float, optional for) – ‘Constant’ method Constant PET value to be spatially distributed.
Tmin (float, required for 'Priestley Taylor' method) – Minimum temperature of the day (deg C)
Tmax (float, required for 'Priestley Taylor' method) – Maximum temperature of the day (deg C)
Tavg (float, required for 'Priestley Taylor' and 'MeasuredRadiationPT') – methods Average temperature of the day (deg C)
float (obs_radiation) – Observed radiation (W/m^2)
method – Observed radiation (W/m^2)
- property const_potential_evapotranspiration#
Constant PET value to be spatially distributed.
Used by ‘Constant’ method.
- property obs_radiation#
Observed radiation (W/m^2)
obs_radiation float, required for ‘MeasuredRadiationPT’ method.
- update()[source]#
Update fields with current conditions.
If the ‘Constant’ method is used, this method looks to the value of the
const_potential_evapotranspiration
property.If the ‘PriestleyTaylor’ method is used, this method looks to the values of the
Tmin
,Tmax
, andTavg
properties.If the ‘MeasuredRadiationPT’ method is use this method looks to the values of the
Tavg
andobs_radiation
property.