Intensity: Intensity

In rjaneUCF/MultiHazard: Tools for modeling compound events

Intensity

Description

Calculates the "intensity" of extreme water levels, as defined in Wahl et al. (2011).

Usage

Intensity(Data, Cluster_Max, Base_Line = "Mean")

Arguments

Data	A data frame containing the water level time series. First column may be a "Date" object.
Cluster_Max	Numeric vector containing indexes of the water level peaks in Data. If analyzing a sample conditioned on water level derived using Con_Sample_2D() set equal to the $xcon output.
Base_Line	Vector of length one, specifying water level about which to calculate the intensity. Default is "Mean" where the mean of the entire time series is used as the baseline water level above which intensity is calculated.

Value

A data frame with the following columns:

• Pre.High Index of the water level column of Data containing the preceding high water level.

• Fol.High Index of the water level column of Data containing the following high water level.

• Pre.Low Index of the water level column of Data containing the preceding low water level.

• Fol.Low Index of the water level column of Data containing the following low water level.

• Intensity Intensity of the extreme water level event.

See Also

Decluster WL_Curve

Examples

#Decluster O-sWL series at S-13 using a runs method
S13.OsWL.Declust = Decluster(Data=S13.Detrend.df$OsWL,
                            SepCrit=24*7, u=0.99667)
#Calculate O-sWL of the identified cluster maximum
intensity = Intensity(Data=S13.Detrend.df,Cluster_Max=S13.OsWL.Declust$EventsMax)
#Plot O-sWL series identifying cluster maximum (in red) and print "intensity" above each maximum
plot(as.Date(S13.Detrend.df$Date_Time),
     S13.Detrend.df$OsWL)
points(as.Date(S13.Detrend.df$Date_Time[S13.OsWL.Declust$EventsMax]),
       S13.Detrend.df$OsWL[S13.OsWL.Declust$EventsMax],pch=16,col=2)
text(as.Date(S13.Detrend.df$Date_Time[S13.OsWL.Declust$EventsMax]),
     S13.Detrend.df$OsWL[S13.OsWL.Declust$EventsMax]+0.2,
     round(intensity$Intensity,0),cex=0.5)

rjaneUCF/MultiHazard documentation built on Sept. 25, 2024, 2:25 p.m.