hotsummerdays: Projection of hot and cold day statistics

In metno/esd: Climate analysis and empirical-statistical downscaling (ESD) package for monthly and daily data

Projection of hot and cold day statistics

Description

The functions hotsummerdays, heatwavespells, coldwinterdays, and coldspells estimate statistics for heatwaves/hot days or cold spells based on seasonal mean temperatures. The estimations are based on a regression analysis (GLM) between observed number of events or spell lengths and seasonal mean from station data. nwetdays estimates the number of days per year with precipitation amount exceeding a threshold values.

Usage

hotsummerdays(
  x,
  y = NULL,
  dse = NULL,
  it = "jja",
  threshold = 30,
  verbose = FALSE,
  plot = TRUE,
  nmin = 90,
  new = TRUE,
  ...
)

Arguments

x	station object, e.g. the temperature. Matches the element used in the dsensemble object 'dse'
y	station object which may be some statistics with dependency to x, e.g. snow depth.
dse	a dsensembel object. If NULL, then run DSensemble
it	Default season set for northern hemisphere. it here~~
threshold	Temperature threshold
verbose	TRUE for trouble shooting, debugging etc.
plot	TRUE - produce graphics
nmin	Minimum number of data points (e.g. days or months) with valid data accepted for annual estimate. NULL demands complete years.
new	if TRUE plot in new window

Details

The estimation of these statistics makes use of general linear models (GLMs) and take the counts to follow the 'Poisson family' whereas the spall lengths belong to the geometric distribution. The seasonal mean temperature or annual wet-mean precipitation are used as independent variable.

Author(s)

R.E. Benestad

Examples

data(ferder)
data(dse.ferder)
hw <- hotsummerdays(ferder,dse.ferder,threshold=20)

metno/esd documentation built on April 29, 2024, 3:34 p.m.