Title: | Evapotranspiration R Recipes |
---|---|
Description: | An R-based application for exploratory data analysis of global EvapoTranspiration (ET) datasets. 'evapoRe' enables users to download, validate, visualize, and analyze multi-source ET data across various spatio-temporal scales. Also, the package offers calculation methods for estimating potential ET (PET), including temperature-based approaches described in : Oudin et al., (2005) <doi:10.1016/j.jhydrol.2004.08.026>. 'evapoRe' supports hydrological modeling, climate studies, agricultural research, and other data-driven fields by facilitating access to ET data and offering powerful analysis capabilities. Users can seamlessly integrate the package into their research applications and explore diverse ET data at different resolutions. |
Authors: | Akbar Rahmati Ziveh [aut, cre]
|
Maintainer: | Akbar Rahmati Ziveh <[email protected]> |
License: | GPL-3 |
Version: | 1.0.0 |
Built: | 2025-04-02 12:39:49 UTC |
Source: | https://github.com/akbarr1184/evapore |
The function detect_exeve
identifies extreme evaporation events based
on standardized evaporation and extreme thresholds.
detect_exeve(x, EXTREMES_THRES = 0.95, LOW_THRES = 0.8)
detect_exeve(x, EXTREMES_THRES = 0.95, LOW_THRES = 0.8)
x |
A |
EXTREMES_THRES |
Numeric. Quantile threshold used to define extreme evaporation events. Default is 0.95. |
LOW_THRES |
Numeric. Lower quantile threshold. Default is 0.80. |
A data.table
with original columns and added columns:
std_value
— standardized evaporation
pentad_std_q80
, pentad_std_q95
— pentad thresholds
value_above_low_thres
, extreme
, evap_event
— flags
above_low_thres_id
, extreme_id
, event_id
— event group IDs
evap <- readRDS("czechia_evap_gleam.rds") events <- detect_exeve(evap)
evap <- readRDS("czechia_evap_gleam.rds") events <- detect_exeve(evap)
The function download_data
downloads the selected data product.
download_data( data_name = "all", path = "", domain = "raw", time_res = "monthly", variable = "e" )
download_data( data_name = "all", path = "", domain = "raw", time_res = "monthly", variable = "e" )
data_name |
a character string with the name(s) of the desired data set. Suitable options are:
|
path |
a character string with the path where the database will be downloaded. |
domain |
a character string with the desired domain data set. Suitable options are:
|
time_res |
a character string with the desired time resolution. Suitable options are:
|
variable |
a character string specifying the variable to download (default = "e"). |
No return value, called to download the required data sets.
download_data("gldas-vic-v2-1", tempdir())
download_data("gldas-vic-v2-1", tempdir())
Downloading Temperature data from different datasets
download_t_data( data_name, path = "", domain = "raw", time_res = "monthly", variable = "all" )
download_t_data( data_name, path = "", domain = "raw", time_res = "monthly", variable = "all" )
data_name |
a character string indicating the dataset to download. Suitable options are:
|
path |
a character string with the path where the data will be downloaded. |
domain |
a character string with the desired domain data set. Suitable options are:
|
time_res |
a character string with the desired time resolution. Suitable options are:
|
variable |
a character string indicating the variable to download. Suitable options are: For TerraClimate dataset:
Use "all" to download all available variables for the dataset. |
No return value, called to download the required data sets.
download_t_data("cru", tempdir())
download_t_data("cru", tempdir())
A subset of GLDAS CLSM monthly Evapotranspiration data in mm over Spain. More detail about raw data can be found here.
gldas_clsm_esp_ts
gldas_clsm_esp_ts
A data.table with 120 obs. of 2 variables:
IDate format %Y-%m-%d
monthly average values
National Aeronautics and Space Administration (NASA)
Global GLDAS monthly Evapotranspiration data in mm. More details of the raw data can be found here.
gldas_clsm_global_ts
gldas_clsm_global_ts
A data.table with 120 obs. of 2 variables:
IDate format %Y-%m-%d
monthly average values
National Aeronautics and Space Administration (NASA)
A subset of GLDAS monthly Evapotranspiration data in mm over -10-40E, 30-45N. More details of the raw data can be found here.
gldas_clsm_subset_ts
gldas_clsm_subset_ts
A data.table with 120 obs. of 2 variables:
IDate format %Y-%m-%d
monthly average values
National Aeronautics and Space Administration (NASA)
The function muldpm
multiplies the value by days per month.
muldpm(x)
muldpm(x)
x |
a RasterBrick object with monthly data in [units/day] |
a RasterBrick object
tavg_brick <- raster::brick('terraclimate_tavg.nc') pet_od <- pet(method = "od", tavg = tavg_brick) pet_od <- muldpm(pet_od)
tavg_brick <- raster::brick('terraclimate_tavg.nc') pet_od <- pet(method = "od", tavg = tavg_brick) pet_od <- muldpm(pet_od)
The function pet
estimates PET by different methods
pet(x, y = NULL, z = NULL, method = "od")
pet(x, y = NULL, z = NULL, method = "od")
x |
a datatable with lon, lat, date, tavg, tmax, and tmin data, a Raster object containing tavg, or a file path to a .nc file containing tavg data. |
y |
a Raster object or file path to a .nc file containing tmax data (optional, required for certain methods). |
z |
a Raster object or file path to a .nc file containing tmin data (optional, required for certain methods). |
method |
a character string indicating the method to be used. Available options are:
|
A Raster object containing potential evapotranspiration in [mm/day], and a data.table object with potential evapotranspiration in [mm/month].
#Calculate PET by Oudin tavg <- raster::brick("terraclimate_tavg_land_19580101_20221231_025_monthly.nc") pet_oudin <- pet(tavg, method = "od") pet_oudin <- muldpm(pet_oudin)
#Calculate PET by Oudin tavg <- raster::brick("terraclimate_tavg_land_19580101_20221231_025_monthly.nc") pet_oudin <- pet(tavg, method = "od") pet_oudin <- muldpm(pet_oudin)
A subset of calculated monthly Potential Evapotranspiration data in mm over Spain. More details of the used method can be found ‘https://www.sciencedirect.com/science/article/pii/S0022169404004056’.
pet_oudin_esp_ts
pet_oudin_esp_ts
A data.table with 120 obs. of 2 variables:
IDate format %Y-%m-%d
monthly average values
Data was calculated using the Oudin method based on raw temperature data. More details of the raw data can be found ‘https://journals.ametsoc.org/view/journals/bams/103/3/BAMS-D-21-0145.1.xml’.
Monthly Potential Evapotranspiration data in mm calculated by Oudin method. More details of the used method can be found ‘https://www.sciencedirect.com/science/article/pii/S0022169404004056’.
pet_oudin_global_ts
pet_oudin_global_ts
A data.table with 120 obs. of 2 variables:
IDate format %Y-%m-%d
monthly average values
Data was calculated using the Oudin method based on raw temperature data. More details of the raw data can be found ‘https://journals.ametsoc.org/view/journals/bams/103/3/BAMS-D-21-0145.1.xml’.
A subset of Monthly Potential Evapotranspiration data in mm calculated by Oudin method over -10-40E, 30-45N. More details of the used method can be found ‘https://www.sciencedirect.com/science/article/pii/S0022169404004056’.
pet_oudin_subset_ts
pet_oudin_subset_ts
A data.table with 120 obs. of 2 variables:
IDate format %Y-%m-%d
monthly average values
Data was calculated using the Oudin method based on raw temperature data. More details of the raw data can be found ‘https://journals.ametsoc.org/view/journals/bams/103/3/BAMS-D-21-0145.1.xml’.