TSA_RD: Detect RD date using the theory and function by Schmid et al....
In geocryology/PermafrostTools: PermafrostTools

Description Usage Arguments Details Value Author(s) Examples

Give one data.frame of a location or a list of data.frames of multiple locations and the function returns the Ripening Date for the snow periods.

1 2	TSA_RD(obs_stat, v = 0.25, v1 = 0.1, v2 = 0.3, lenghtsnow = 5, MDr.sd = 0.3, snow.period = NULL)

`obs_stat`	One dataframe of one location or a list of dataframes of multiple locations. A dataframe needs at least the following columns: Location name as loc_name; Mean daily GST as agg_avg; Minimum daily GST as agg_min; Maximum daily GST as agg_max; Daily sd as agg_sd; Date as time (in POSIXct and the following format: YYYY-MM-DD)
`v`	Threshold zero curtain (if dailymax & dailymin is within +/- v, this day is a zero curtain day) Default: 0.25
`v1`	Max. daily standard deviation indicating snow for POSITIVE Ground Surface Temperature. Default is 0.1.
`v2`	Max. daily standard deviation indicating snow for NEGATIVE Ground Surface Temperature. Default is 0.3.
`MDr.sd`	Threshold for the mean daily standard deviation for a snow period.
`snow.period`	In case the snow period is already available for one location, it can be inserted here. The data frame of the now period needs to have the following columns: loc_name, Snow.Start, Snow.End, MDr
`lengthsnow`	Number of days that the snow cover should at least have to be selected as a valid snow cover. Default is 5 days.

It is also possible to adjust the parameters for the TSA_snow_cover function in this function, as the warming periods are calculated for each snow cover period. Therefore, the TSA_snow_cover is run as well in this function.

Dataframe with fields: loc_name; Snow.Start; Snow.End; MDr

Thomas Knecht <t.knecht@hotmail.com>

1
2
3

con <- dbpf_con()
BC_SO01_01 <- TSA_data_import(con, "BC16-SO01_01")
RD_table <- TSA_RD(BC_SO01_01, v= 0.25, m1 = 3, m2 = 0.5, v1 = 0.1, v2 = 0.3, lenghtsnow = 5, MDr.sd = 0.3, snow.period=NULL)