pRsim.weather {PRSim}

Weather simulation (temperature and precipitation) for multiple stations

Description

Applies the wavelet-based weather simulation algorithm to multiple sites (single site possible as well)

Usage

prsim.weather(data_p, data_t, station_id_p="Precip",
station_id_t="Temp", number_sim=1, win_h_length=15,
n_wave=100,verbose=TRUE,t_margin='sep',p_margin='egpd',...)

Arguments

data_p

list of precipitation data frames. One list entry, i.e. data frame, corresponds to one station/grid cell. Each data frame contains the time indications and precipitation of one station. See ‘Details’.

data_t

list of temperature data frames. One list entry, i.e. data frame, corresponds to one station/grid cell. Each data frame contains the time indications and temperature of one station. See ‘Details’.

station_id_p

identifies the precipitation variable name in case several time series are present in data_p. See ‘Details’.

station_id_t

identifies the temperature variable name in case several time series are present in data_t. See ‘Details’.

number_sim

number of simulations to be carried out.

win_h_length

(half-)length of moving window size.

t_margin

marginal distribution to be used for the backtransformation of temperature. Can be either "sep" or any type of CDF (see ‘Details’). "sep" uses the four-parameter skewed-exponential power for backtransformation. CDF allows for specifying any distribution ‘Examples’.

p_margin

marginal distribution to be used for the backtransformation of precipitation. Can be either "egpd" or any type of CDF (see ‘Details’). "egpd" uses the extended GPD for backtransformation. CDF allows for specifying any distribution ‘Examples’.

verbose

logical. Should progress be reported?

n_wave

number of scales to be considered in the continuous wavelet transform.

...

any other argument passed to the sub-function specifying the cdf for fitting. See ‘Details’ and ‘Examples’.

Details

Time can be given with three columns named "YYYY", "MM", "DD", or as in POSIXct format YYYY-MM-DD. All leap days (Feb 29th) will be omitted from the analysis, but no missing observations are allowed.

Stations are identified by list index.

Alternative distributions can be specified by providing three functions: (1) a function fitting the parameters of a distributions and providing a vector of these parameters as output (CDF_fit), (2) a function simulating random numbers from this distribution (rCDF), and (3) a function specifying the distribution (pCDF). See ‘Examples’ for the generalized beta for the second kind and for the Generalized Extreme Values (GEV) distribution.

Value

A list with elements temperature and precipitation of

simulation

data frames with time information, observations, and number_sim columns containing the simulated data.

Author(s)

Manuela Brunner

References

Brunner, M. I., and E. Gilleland (2021). Spatial compound hot-dry events in the United States: assessment using a multi-site multi-variable weather generator, in preparation.