Statistics of extreme stormsurge levels at tidegauges
|
|
- Noreen Wilkerson
- 5 years ago
- Views:
Transcription
1 Statistics of extreme stormsurge levels at tidegauges along the Dutch coast Douwe Dillingh National Institute for Coastal and Marine Management / RIKZ P.O. Box 20907, 2500 EX The Hague The Netherlands d.dillingh@rikz.rws.minvenw.nl Abstract The basic standard for the design of a seadefence structure in the Dutch coastal region is the stormsurge level with a frequency of exceedance of 10"* per year, the so-called basic design level. A new extensive and thorough statistical research for these levels was done on the datasets of five main tidal stations. Lot of attention was given to the preprocessing of the data in order to obtain independent and identically distributed data. Different extreme value models were applied, such as the Generalised Extreme Value Distribution and the Generalised Pareto Distribution. Also a modern distribution free extreme value model was considered. New frequency curves for all relevant tide gauge stations were derived, not only based on statistics, but also on the physical coherence using hydrodynamic models (separate topic by Marc E. Philippart on this congress). 1 Introduction After the disastrous stormsurge during the night of february 1* 1953 the Delta Committee (1960) determined the water levels which served as the basis for improvement of the major seadefence structures in the Dutch coastal region: the stormsurge levels with a frequency of exceedance of 10~* per year, the basic design levels (Delta Committee*). In view of the uncertainties in these, particularly for the western Wadden Sea because of the short time series available as a consequence of the radical changing of the tidal regime by the closure of the Zuiderzee in 1932, the Delta
2 450 Coastal Engineering and Marina Developments Committee recommended to continue the studies with the use of new data. Figure 1: Position of the main tidal stations along the Dutch coast, the light area of the Netherlands is vulnerable for flooding without safety measures. In 1984 an extensive and thorough statistical research was started on the datasets of the five main tidal stations Vlissingen, Hoek van Holland, Den Helder, Harlingen and Delfzijl. Datasets with a time span of about hundred years were available, except for Den Helder and Harlingen, for which timeseries were usable ever since 1932 (western Wadden Sea). The frequencies of exceedance at interest were in the range of 10"* to 1(T*. The locations are marked in fig Approach to the problem The aim is to estimate the probability of exceedence of extreme high water levels that have not occurred yet. Application of extreme value theory demands independent and identically distributed (i.i.d.) data. Two kinds of data were considered: the high water level (HW) and the wind set up (HW-setup), defined as the difference between the HW and the corresponding astronomical high water, regardless a timeshift. The preprocessing of the data is described in the next section. Four statistical models were applied to the datasets of the five main stations, one is chosen. The statistical results were checked by considering the coherence of the stormsurge levels between the stations, using hydrodynamic models. Starting
3 Coastal Engineering and Marina Developments 451 from the statistically determined 1(T* quantile of Hoek van Holland other estimates were obtained for the 10~* quantile. The results of both estimates (statistics and hydrodynamic models) were integrated by averaging with the use of the standarddeviation for weighting factors. In this way 10"^ quantiles were determined for the stations Delfzijl, Harlingen, Den Helder, Vlissingen, West-Terschelling, Umuiden, Terneuzen en Hansweert (see fig. 1). Starting from the thus known 1CT* quantile and the empericai frequency distributions (from the i.i.d. data) frequency curves were established from which also other quantiles may be read. The hydrodynamic models were also used to interpolate between these stations. With the results frequency curves could be derived for the other tide gauges for which (relatively) short timeseries were available. 3 Preprocessing of the data The processing of the measured data (high water levels) to obtain i.i.d. data included the following steps. 1. Small gaps in the data were filled up by estimations, using relations with surrounding tidegauges. All datasets were checked on inaccuracies and consistency of filtering secundary fluctuations. All astronomical high water levels were recalculated, using one method for all years. 2. The wish to deal with real HW setups caused by wind led to the introduction of a threshold value for the HW setup. A minimum level of 30 cm was holded on basis of inaccuracies in the measured HW's and calculated astronomical HW's and of and airpressure variations. Empericai frequency curves (see fig. 4) of high waters show a deviating behaviour (stronger curvature) at the lower levels. Beside that low levels are here of little interest. The threshold is taken as the highest astronomical HW in the concerned dataset plus 30 cm. 3. The data should be samples from the same distribution function. The distribution of the HW setups, the stochastic part of the high water levels, caused by meteorological phenomena, of the stations Hoek van Holland and Delfzijl were visualised by means of boxplots. A threshold level for HW setup of 70 cm was (rather arbitrarily) applied. Boxplots were made for periods of a month, starting on the 1* and the 15* day. By putting the boxplots in line one gets a good impression of the behaviour of the distribution function over the year. The period from 1 October up to 15 march was chosen as the stormseason. This result was checked later with selected data for autocorreation. 4. Autocorrelation in the series of HW setups exists on different timescales: from within a single storm until correation between circulationpattems. Removal of all autocorrelation leaves little data left. It was tried to remove most autocorrelation by demanding that every selected HW setup must be larger than the i preceding and and i succeeding ones. The best way to estimate i turned out to be a clusteranalysis (Diggle^). In case of independancy of HW setups there is now clustering in time, no forming of
4 452 Coastal Engineering and Marina Developments groups. The analysis led to the choice i=4. Selection of HW levels was done through the selection of the corresponding HW setups. 5. Until sofar it was assumed that there was no change of the storminess (i.e. the climate of the HW setups during the considered observation period). This assumption had to be checked. In fig. 2 the number of selected HW setups per stormseason (with threshold level of 90 cm) for the station Delfzijl has been plotted against the year in which the month of January of the considered stormseason falls (more thresholds were applied). These figure clearly shows no trend. This was confirmed by a statistical test (Spearman test). average Figure 2: Delfzijl 1881/ /91, number of selected HW stormseason setups > 90 cm per 6. Plots of the yearly mean high water levels show a claerly rising trend, relative to the NAP, the Dutch ordnance datum, which lies at present about mean sea level. Because of human interventions and morfological changes the mean high water levels rise somewhat faster than the mean sea levels. Per station the dataset had to be corrected for the locally appeared changes of the mean high water levels, which were generally gradual but sometimes abrupt. The gradual changes were approximated by linear regression. All high waterlevels were corrected to the situation of 1985, the end of the used timeseries. 7. If the extreme value theory is applied to the HW setup instead of the HW levels, than the resulting distribution has to be combined with that of the astronomical high waters in order to obtain a useful final result. For this convolution independence of both parameters is demanded. This was checked by plotting and by testing (Spearman test with ties). Independance was concluded. 4 Statistical models For a good understanding this section gives a short description of the relevant general theory. More details (parameter estimation and confidence intervals) are to be found in the literature. The concern of this research is to estimate a high water level such that the probability of reaching or exceeding that level in a arbitrarily given year equals a chosen value po. Values between 10"^ and 10~* are of interest. For this N i.i.d. high water levels XI,...,XN are available. It is assumed that the probability of exceedance of extreme high water levels outside the stormseason is negligable to that within the stormseason, so that in the following we speak about
5 stormseason instead of year. Coastal Engineering and Marina Developments 453 The data XI,...,XN are samples from an unknown distribution F. If q is the average number of observations per stormseason (q=n/s), than for the stormseason maximum M* of a abitrary stormseason applies From this follows For the probability of exceedance p from F of XQ applies F(XQ) = 1-p, so that For small values of po follows from equations (2) and (3) and If F were known, then the demanded level XQ could be calculated from (4) ^=F-V--; «(5) F* is the inverse function of F. Distribution free method (DF) This method was developed during the present research (De Haan^). Be Xi, Xz,... i.i.d. data with distribution F, then for the probability distribution of the maximum Mn of the first n quantities Xi,...,Xn P( M»<X) = P( X, < X,..., Xn<x) = F"(X) ^ From this follows for every x for which F(x) < 1 limp(mn<x) = 0 /I-^«> It is said that this equation degenerates in its limit. By introducing properly chosen constants an>0 en b,, the above expression converges for the most known distribution functions F to one of the so-called extreme value distributions G?: an (?) for large n, in which
6 454 Coastal Engineering and Marina Developments for all x for wich 1+yx > 0. y is called the extreme vakue index. y=0 yields the well-known Gumbel distribution (Go(x) = exp(-e~*)). These distributions F are said to belong to the domain of attraction of Gy. For the application for the issue in question the quotient n/k is important, k is a properly chosen numbers largest observations of the series of n observations. It should hold that n/k ><» if n»<*>, which means that with increasing n, k increases slower than n. If n/k is larg enough than and so n Xo-bl i j5i = 7_fY;cJ= j.fg/ - ^J7^ q a± * (10) After substitution of (8) in (10) and with l-e~* «x for small x and some rewriting: k * ( ) -1 np (ID The parameters y, a^ en b^ have to be estimated from the data (De Haan^). GPD method From the i.i.d. data Xi,...,X% with distribution function F, those data are selected that exceed a proper chosen threshold u. For sufficient high u the distribution FU of the u exceeding parts of the data, YI,...,YM, is well approximated by the socalled GPD distribution (Generalised Pareto Distribution), here denoted as G(y;a(u),y) (Pickands* and SmiuY). Condition is that F belongs to the domain of attraction of the extreme value distributions Gy. So for large u -7 <? (12) for all y for which l+y(y/a) > 0. The parameter y is here the same as the extreme value index of the extreme value Gy y^o yields the exponential distribution function: G(y; a,0) = l-exp(-y/a). Following the procedure for the derivation of (5) and with yo=xo-u en qu=m/s (average number exceedances of u per stormseason), then and so <? <7,,,y) (13) W / i / * O \ v xo = y<>+u ~ - i^-( ) J+w 7 <*u (14) (7u and y are estimated (maximum likelyhood method) for properly chosen u.
7 Coastal Engineering and Marina Developments 455 GED method The GED model (Generalised Extreme value Distribution) uses maxima in periods of equal length, usually a year. Here stormseasonmaxima are considered. Formula (7) is applied to the maximum X^a of each stormseason i (i=l,...,s). It is assumed that the distribution of de Xma%,i's, if properly normalised, equals exactly the extreme value distribution G?. If the constants to be introduced are denoted by U, (place) and a (scale) then r,/ P( max,/ in which G-^x) is given by formula (8). Formule (15) can be rewriten as x-fl (16) With formula (1) an estimate of XQ, the po-quantile is obtained by solution of ' P.- G,<^> " (17) And after some conversion (15) The parameters are estimated with the maximum likelyhood method. (18) 5 Results DF model Estimates for the 1(T* quantile of the high water levels are calculated for an increasing number of k, the number of the highest observations to be considered, and plotted against k. Fig. 3 gives as an example the result for Delfzijl. For small k the estimates behave instable as a function of k. Further k should be smaller than the value that corresponds with the threshold level. The estimate is chosen with the largest value of k in the interval in which the quantile estimates behave stable (i.e. independant of k). GPD model In principle the same procedure is followed for this model as for the DF method. The GPD model works with thresholds, but for every threshold a number of highest obsevations is connected. GED model In this model k plays no role, which makes the results less subjective.. The biggest disadvantage of this method is that a lot of information is not used.
8 456 Coastal Engineering and Marina Developments Figure 3: Delfzijl 188 1/ /85, DF-estimate of the 1CT* quantile of HW Convolution method (CON) In this method the GPD model is applied to the HW setups. Joining of the distributions function of HW setup and atronomical high waters to a probability finding about the high water level is done by means of the so-called convolution integral. The further procedure is the same as for the DF method and the GPD model. The confidence intervals are calculated by a parametric bootstrap method, using simulations according to a Monte Carlo procedure. The results of the four methods are given in table 1. Table 1: Estimates of the 10"* quantile (cm above NAP) and the corresponding 95% confidence intervals method DF GPD CON GED Vlissingen 540 ± ± ± ±185 Hoek van Holland 500 ± ± ± ±180 Den Helder 370 ± ± ± ±110 Harlingen 405 ± ± ± ± 65 Delfzijl 600 ± ± ± ±195 The confidence intervals are large, on the one side because the extrapolation is large, but especially because y was not given a certain value in advance. They are also symmetric (except the CON method) due to the way of calculation (asymptically normal), but will not be in reality. Preference was given to the results of the distribution free method (DF), because no parametric distribution is imposed on the data. The underboundary shows a descending path (fig 3). A
9 Coastal Engineering and Marina Developments 457 bootstrap approach for the DF confidence intervals yielded a better image. 6 Combination ofstatistics with physical coherence Combination of statistical results (DF method for all stations with long timeseries) and the results of the numerical models for the coherence yielded the basic design levels (10"* quantiles) that are given in table 2. Table 2: Basic design levels (10"* quantiles): Delta Committee, statistical result and final result. Station Vlissingen Terneuzen Hansweert Hoek van Holland Umuiden Den Helder West Terschelling Harlingen Delfzijl Delta Committee (situation 1950) Statistical result (situation 1985) Final result (situation 1985) Specially for the Western Waddensea the new basic design level are considerably lower than those of the Delta Committee (the Delta Committee used in fact y=0). 7 Frequency curves A frequency curve gives the relation between the high water level and the average number per stormseason (or year) that level is reached or exceeded. It has to link smoothly to the selected observations. The DF model is not fit for constructing such lines. Instead they are constructed on the basis of the GPD model. The GPD model has three parameters: threshold and scale and shape parameter. So three points have to be known. For all stations the levels with frequency of exceedance of 5.10"% 10"^ and 10"* respectively were starting point. The first two were derived from the data (selected and corrected), taking into account the mutual coherence of the stations and the available, sometimes very limited, length of the timeseries. The threshold for all stations is the relatively high level with frequency of exceedance of 5.10"*. These three points of the frequency curve have to be converted to corresponding points of the GPD distribution. If F»(y) is the distribution function of y;, being the part of the observations X; that exceeds the threshold level u (yj = Xj - u), then one can argue that 4, Fu(y) ~ 1- for x > u (19) 4u with q% = the average number of exceedances per stormseason of level x
10 458 Coastal Engineering and Marina Developments qu = the average number of exceedances per stormseason of threshold u With formula (12) follows x-u for x > u (20) Here for the chosen thresholds of the frequency curves curve for the station Delfzijl as an example. Delfzijl, frequency curve, situation 1985 *** emperical frequency curve for all high water levels ooo emperical frequency curve for selected high water levels calculated frequency curve, = 0.5. Fig. 4 shows the ^^^^, o^^* ^ c. + - a < Z -1 ^. p#* «r" "^ ^ ^ ^ S>,.^> %.. * > / f a\ rera jcnui ' nbeirofexceedances] «ryear/s ormseason UU- Figure 4: Delfzijl, emperical and calculated frequency curve (situation 1985) The results are valid for the situation in Rise of mean high water did not stop since. On basis of the measured rise in the past few decades all frequency curves (so incuding the design levels) are shifted up 5 cm, valid until the year Reports of the complete study are available (only in Dutch) at the RIKZ. References 1. Rapport Deltacommissie, Eindverslag en Interimadviezen, Deel I, Staatsdrukkerij- en Uitgeverijbedrijf, 's-gravenhage, 1960 (in Dutch). 2. Diggle P.J.,Statistical Analysis of Spatial Point Patterns, Academic Press, De Haan L., Fighting the arch-enemy with mathematics, Statistica Neerlandica, 44, 2, p.45-68, Pickands J., Statistical inference using extreme order statistics, Annals of Statistics, Vol.3, no.l, p , Smith R.L., Estimating tails of probability distributions, The Annals of Statistics, 15, p , 1987.
Uncertainties in extreme surge level estimates from observational records
Uncertainties in extreme surge level estimates from observational records By H.W. van den Brink, G.P. Können & J.D. Opsteegh Royal Netherlands Meteorological Institute, P.O. Box 21, 373 AE De Bilt, The
More informationPost flood analysis. Demands of the Storm Surge Warning Service (SVSD) of the Netherlands. Annette Zijderveld/ Martin Verlaan
Post flood analysis Demands of the Storm Surge Warning Service (SVSD) of the Netherlands Annette Zijderveld/ Martin Verlaan Objective Bringing in coastal issues to the FEWS user days (second attempt) Demonstrating
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationBayesian Modelling of Extreme Rainfall Data
Bayesian Modelling of Extreme Rainfall Data Elizabeth Smith A thesis submitted for the degree of Doctor of Philosophy at the University of Newcastle upon Tyne September 2005 UNIVERSITY OF NEWCASTLE Bayesian
More informationEstimation techniques for inhomogeneous hydraulic boundary conditions along coast lines with a case study to the Dutch Petten Sea Dike
Estimation techniques for inhomogeneous hydraulic boundary conditions along coast lines with a case study to the Dutch Petten Sea Dike P.H.A.J.M. van Gelder, H.G. Voortman, and J.K. Vrijling TU Delft,
More informationOn the Application of the Generalized Pareto Distribution for Statistical Extrapolation in the Assessment of Dynamic Stability in Irregular Waves
On the Application of the Generalized Pareto Distribution for Statistical Extrapolation in the Assessment of Dynamic Stability in Irregular Waves Bradley Campbell 1, Vadim Belenky 1, Vladas Pipiras 2 1.
More informationOn the lack of information in hydraulic engineering models
On the lack of information in hydraulic engineering models Karen A.H. Slijkhuis Risk Analysis Department, Ministry of Public Works and Water Management, Utrecht, the Netherlands Pieter H.A.J.M. van Gelder
More informationHIERARCHICAL MODELS IN EXTREME VALUE THEORY
HIERARCHICAL MODELS IN EXTREME VALUE THEORY Richard L. Smith Department of Statistics and Operations Research, University of North Carolina, Chapel Hill and Statistical and Applied Mathematical Sciences
More informationStatistics for extreme & sparse data
Statistics for extreme & sparse data University of Bath December 6, 2018 Plan 1 2 3 4 5 6 The Problem Climate Change = Bad! 4 key problems Volcanic eruptions/catastrophic event prediction. Windstorms
More informationDATA ASSIMILATION STUDY USING THE DUTCH CONTINENTAL SHELF MODEL WITH FULL MEASUREMENT
DATA ASSIMILATION STUDY USING THE DUTCH CONTINENTAL SHELF MODEL WITH FULL MEASUREMENT JULIUS H. SUMIHAR AND MARTIN VERLAAN,2 Faculty of Electrical Engineering, Mathematics and Computer Science, Department
More informationIT S TIME FOR AN UPDATE EXTREME WAVES AND DIRECTIONAL DISTRIBUTIONS ALONG THE NEW SOUTH WALES COASTLINE
IT S TIME FOR AN UPDATE EXTREME WAVES AND DIRECTIONAL DISTRIBUTIONS ALONG THE NEW SOUTH WALES COASTLINE M Glatz 1, M Fitzhenry 2, M Kulmar 1 1 Manly Hydraulics Laboratory, Department of Finance, Services
More informationOverview of Extreme Value Theory. Dr. Sawsan Hilal space
Overview of Extreme Value Theory Dr. Sawsan Hilal space Maths Department - University of Bahrain space November 2010 Outline Part-1: Univariate Extremes Motivation Threshold Exceedances Part-2: Bivariate
More informationA Conditional Approach to Modeling Multivariate Extremes
A Approach to ing Multivariate Extremes By Heffernan & Tawn Department of Statistics Purdue University s April 30, 2014 Outline s s Multivariate Extremes s A central aim of multivariate extremes is trying
More informationOverview of Extreme Value Analysis (EVA)
Overview of Extreme Value Analysis (EVA) Brian Reich North Carolina State University July 26, 2016 Rossbypalooza Chicago, IL Brian Reich Overview of Extreme Value Analysis (EVA) 1 / 24 Importance of extremes
More informationEstimation of Quantiles
9 Estimation of Quantiles The notion of quantiles was introduced in Section 3.2: recall that a quantile x α for an r.v. X is a constant such that P(X x α )=1 α. (9.1) In this chapter we examine quantiles
More informationON THE TWO STEP THRESHOLD SELECTION FOR OVER-THRESHOLD MODELLING
ON THE TWO STEP THRESHOLD SELECTION FOR OVER-THRESHOLD MODELLING Pietro Bernardara (1,2), Franck Mazas (3), Jérôme Weiss (1,2), Marc Andreewsky (1), Xavier Kergadallan (4), Michel Benoît (1,2), Luc Hamm
More informationKNMI-HYDRA project. Phase report 7. Estimation of extreme return levels of wind speed: an analysis of storm maxima
KNMI-HYDRA project Phase report 7 Estimation of extreme return levels of wind speed: an analysis of storm maxima KNMI, May 03 Estimation of extreme return levels of wind speed: an analysis of storm maxima
More informationSome conditional extremes of a Markov chain
Some conditional extremes of a Markov chain Seminar at Edinburgh University, November 2005 Adam Butler, Biomathematics & Statistics Scotland Jonathan Tawn, Lancaster University Acknowledgements: Janet
More informationRichard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC
EXTREME VALUE THEORY Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu AMS Committee on Probability and Statistics
More informationMFM Practitioner Module: Quantitiative Risk Management. John Dodson. October 14, 2015
MFM Practitioner Module: Quantitiative Risk Management October 14, 2015 The n-block maxima 1 is a random variable defined as M n max (X 1,..., X n ) for i.i.d. random variables X i with distribution function
More informationModels and estimation.
Bivariate generalized Pareto distribution practice: Models and estimation. Eötvös Loránd University, Budapest, Hungary 7 June 2011, ASMDA Conference, Rome, Italy Problem How can we properly estimate the
More informationModeling the North West European Shelf using Delft3D Flexible Mesh
Modeling the North West European Shelf using Delft3D Flexible Mesh 2nd JCOMM Scientific and Technical Symposium on Storm Surges, 8-13 Nov. 2015, Key West, USA; Firmijn Zijl Outline of this presentation
More informationBivariate generalized Pareto distribution
Bivariate generalized Pareto distribution in practice Eötvös Loránd University, Budapest, Hungary Minisymposium on Uncertainty Modelling 27 September 2011, CSASC 2011, Krems, Austria Outline Short summary
More informationPERFORMANCE OF PARAMETER ESTIMATION TECHNIQUES WITH INHOMOGENEOUS DATASETS OF EXTREME WATER LEVELS ALONG THE DUTCH COAST.
PERFORMANCE OF PARAMETER ESTIMATION TECHNIQUES WITH INHOMOGENEOUS DATASETS OF EXTREME WATER LEVELS ALONG THE DUTCH COAST. P.H.A.J.M. VAN GELDER TU Delft, Faculty of Civil Engineering, Stevinweg 1, 2628CN
More informationSEDIMENT IMPORT BY TIDAL INLETS SEDBOX -model for tidal inlets Marsdiep and Vlie, Wadden sea, The Netherlands by L.C. van Rijn
SEDIMENT IMPORT BY TIDAL INLETS SEDBOX -model for tidal inlets Marsdiep and Vlie, Wadden sea, The Netherlands by L.C. van Rijn 1. Physical processes of sandy tidal inlet systems 1.1 Introduction 1.2 Hydrodynamics
More informationModelação de valores extremos e sua importância na
Modelação de valores extremos e sua importância na segurança e saúde Margarida Brito Departamento de Matemática FCUP (FCUP) Valores Extremos - DemSSO 1 / 12 Motivation Consider the following events Occurance
More informationVerification of WAQUA/DCSMv5 s operational water level probability forecasts
Ministry of Infrastructure and Water Management Verification of WAQUA/DCSMv5 s operational water level probability forecasts N. Wagenaar KNMI Internal report IR-18-1 UTRECHT UNIVSERSITY DEPARTMENT OF
More informationFinancial Econometrics and Volatility Models Extreme Value Theory
Financial Econometrics and Volatility Models Extreme Value Theory Eric Zivot May 3, 2010 1 Lecture Outline Modeling Maxima and Worst Cases The Generalized Extreme Value Distribution Modeling Extremes Over
More informationProbability forecasts for water levels at the coast of The Netherlands
from Newsletter Number 114 Winter 7/8 METEOROLOGY Probability forecasts for water levels at the coast of The Netherlands doi:1.1957/gpsn56sc This article appeared in the Meteorology section of ECMWF Newsletter
More informationWei-han Liu Department of Banking and Finance Tamkang University. R/Finance 2009 Conference 1
Detecting Structural Breaks in Tail Behavior -From the Perspective of Fitting the Generalized Pareto Distribution Wei-han Liu Department of Banking and Finance Tamkang University R/Finance 2009 Conference
More informationWIND TRENDS IN THE HIGHLANDS AND ISLANDS OF SCOTLAND AND THEIR RELATION TO THE NORTH ATLANTIC OSCILLATION. European Way, Southampton, SO14 3ZH, UK
J 4A.11A WIND TRENDS IN THE HIGHLANDS AND ISLANDS OF SCOTLAND AND THEIR RELATION TO THE NORTH ATLANTIC OSCILLATION Gwenna G. Corbel a, *, John T. Allen b, Stuart W. Gibb a and David Woolf a a Environmental
More informationMath 576: Quantitative Risk Management
Math 576: Quantitative Risk Management Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 11 Haijun Li Math 576: Quantitative Risk Management Week 11 1 / 21 Outline 1
More informationAPPLICATION OF EXTREMAL THEORY TO THE PRECIPITATION SERIES IN NORTHERN MORAVIA
APPLICATION OF EXTREMAL THEORY TO THE PRECIPITATION SERIES IN NORTHERN MORAVIA DANIELA JARUŠKOVÁ Department of Mathematics, Czech Technical University, Prague; jarus@mat.fsv.cvut.cz 1. Introduction The
More informationEstimating Bivariate Tail: a copula based approach
Estimating Bivariate Tail: a copula based approach Elena Di Bernardino, Université Lyon 1 - ISFA, Institut de Science Financiere et d'assurances - AST&Risk (ANR Project) Joint work with Véronique Maume-Deschamps
More informationNSW Ocean Water Levels
NSW Ocean Water Levels B Modra 1, S Hesse 1 1 Manly Hydraulics Laboratory, NSW Public Works, Sydney, NSW Manly Hydraulics Laboratory (MHL) has collected ocean water level and tide data on behalf of the
More informationMETOCEAN CRITERIA FOR VIRGINIA OFFSHORE WIND TECHNOLOGY ADVANCEMENT PROJECT (VOWTAP)
METOCEAN CRITERIA FOR VIRGINIA OFFSHORE WIND TECHNOLOGY ADVANCEMENT PROJECT (VOWTAP) Report Number: C56462/7907/R7 Issue Date: 06 December 2013 This report is not to be used for contractual or engineering
More informationTIME-DEPENDENT EFFECTS OF NOURISHMENTS ON SHOREFACE BAR BEHAVIOUR. Abstract
TIME-DEPENDENT EFFECTS OF NOURISHMENTS ON SHOREFACE BAR BEHAVIOUR Tommer Vermaas 1, Edwin Elias 2, Ad van der Spek 3 and Rena Hoogland 4 Abstract In 2011/2012 a shoreface nourishment was placed at Heemskerk,
More informationStudy Of The Wind Speed, Rainfall And Storm Surges For The Scheldt Estuary In Belgium
Study Of The Wind Speed, Rainfall And Storm Surges For The Scheldt Estuary In Belgium Mohammad Abul Hossen, Farjana Akhter Abstract: The Belgian coast and the Scheldt estuary are important for the Belgian
More informationBest Fit Probability Distributions for Monthly Radiosonde Weather Data
Best Fit Probability Distributions for Monthly Radiosonde Weather Data Athulya P. S 1 and K. C James 2 1 M.Tech III Semester, 2 Professor Department of statistics Cochin University of Science and Technology
More informationExploring extreme sea level events
PSMSL 80 th Anniversary Workshop on Sea Level Science Exploring extreme sea level events Melisa Menéndez Environmental Hydraulics Institute (IH-Cantabria), Universidad de Cantabria, Spain menendezm@unican.es
More informationClassical Extreme Value Theory - An Introduction
Chapter 1 Classical Extreme Value Theory - An Introduction 1.1 Introduction Asymptotic theory of functions of random variables plays a very important role in modern statistics. The objective of the asymptotic
More informationSharp statistical tools Statistics for extremes
Sharp statistical tools Statistics for extremes Georg Lindgren Lund University October 18, 2012 SARMA Background Motivation We want to predict outside the range of observations Sums, averages and proportions
More informationThe Spatial Variation of the Maximum Possible Pollutant Concentration from Steady Sources
International Environmental Modelling and Software Society (iemss) 2010 International Congress on Environmental Modelling and Software Modelling for Environment s Sake, Fifth Biennial Meeting, Ottawa,
More informationFlood Frequency Mapping using Multirun results from Infoworks RS applied to the river basin of the Yser, Belgium
Flood Frequency Mapping using Multirun results from Infoworks RS applied to the river basin of the Yser, Belgium Ir. Sven Verbeke Aminal, division Water, Flemish Government, Belgium Introduction Aminal
More informationOn the modelling of extreme droughts
Modelling and Management of Sustainable Basin-scale Water Resource Systems (Proceedings of a Boulder Symposium, July 1995). IAHS Publ. no. 231, 1995. 377 _ On the modelling of extreme droughts HENRIK MADSEN
More informationEvaluating Hydrodynamic Uncertainty in Oil Spill Modeling
Evaluating Hydrodynamic Uncertainty in Oil Spill Modeling GIS in Water Resources (CE 394K) Term Project Fall 2011 Written by Xianlong Hou December 1, 2011 Table of contents: Introduction Methods: Data
More informationAnalysis methods of heavy-tailed data
Institute of Control Sciences Russian Academy of Sciences, Moscow, Russia February, 13-18, 2006, Bamberg, Germany June, 19-23, 2006, Brest, France May, 14-19, 2007, Trondheim, Norway PhD course Chapter
More informationAnalysing River Discharge Data for Flood Analysis
Analysing River Discharge Data for Flood Analysis This suggested response provides answers to the eleven questions posed on the web-site. Your report should have covered all of this material. (A) Produce
More informationBayesian Inference for Clustered Extremes
Newcastle University, Newcastle-upon-Tyne, U.K. lee.fawcett@ncl.ac.uk 20th TIES Conference: Bologna, Italy, July 2009 Structure of this talk 1. Motivation and background 2. Review of existing methods Limitations/difficulties
More informationUsing statistical methods to analyse environmental extremes.
Using statistical methods to analyse environmental extremes. Emma Eastoe Department of Mathematics and Statistics Lancaster University December 16, 2008 Focus of talk Discuss statistical models used to
More informationEmma Simpson. 6 September 2013
6 September 2013 Test What is? Beijing during periods of low and high air pollution Air pollution is composed of sulphur oxides, nitrogen oxides, carbon monoxide and particulates. Particulates are small
More informationQuantitative Modeling of Operational Risk: Between g-and-h and EVT
: Between g-and-h and EVT Paul Embrechts Matthias Degen Dominik Lambrigger ETH Zurich (www.math.ethz.ch/ embrechts) Outline Basel II LDA g-and-h Aggregation Conclusion and References What is Basel II?
More informationMediterranean Sea Level Changes over the Period
Mediterranean Sea Level Changes over the Period 1961-2 Boris SHIRMAN and Yossi MELZER, Israel Key words: Sea level monitoring, Tide gauge operation, Tide gauge Bench Marks. ABSTRACT The mean value of the
More informationIntroduction to Algorithmic Trading Strategies Lecture 10
Introduction to Algorithmic Trading Strategies Lecture 10 Risk Management Haksun Li haksun.li@numericalmethod.com www.numericalmethod.com Outline Value at Risk (VaR) Extreme Value Theory (EVT) References
More informationZwiers FW and Kharin VV Changes in the extremes of the climate simulated by CCC GCM2 under CO 2 doubling. J. Climate 11:
Statistical Analysis of EXTREMES in GEOPHYSICS Zwiers FW and Kharin VV. 1998. Changes in the extremes of the climate simulated by CCC GCM2 under CO 2 doubling. J. Climate 11:2200 2222. http://www.ral.ucar.edu/staff/ericg/readinggroup.html
More informationQuantile-quantile plots and the method of peaksover-threshold
Problems in SF2980 2009-11-09 12 6 4 2 0 2 4 6 0.15 0.10 0.05 0.00 0.05 0.10 0.15 Figure 2: qqplot of log-returns (x-axis) against quantiles of a standard t-distribution with 4 degrees of freedom (y-axis).
More informationSouth San Francisco Bay Shoreline Studies for EIA 11 with Project Conditions
South San Francisco Bay Shoreline Studies for EIA 11 with Project Conditions U.S. Army Corps of Engineers San Francisco District Ms. Lisa Andes Mr. Craig Conner Dr. Frank Wu Dr. Jen-Men Lo Dr. Michael
More informationCOMPOUND FLOODING: EXAMPLES,
1 COMPOUND FLOODING: EXAMPLES, METHODS, AND CHALLENGES Thomas Wahl Civil, Environmental and Construction Engineering & National Center for Integrated Coastal Research University of Central Florida Definition(s)
More informationFrequency Estimation of Rare Events by Adaptive Thresholding
Frequency Estimation of Rare Events by Adaptive Thresholding J. R. M. Hosking IBM Research Division 2009 IBM Corporation Motivation IBM Research When managing IT systems, there is a need to identify transactions
More informationFOWPI Metocean Workshop Modelling, Design Parameters and Weather Windows
FOWPI Metocean Workshop Modelling, Design Parameters and Weather Windows Jesper Skourup, Chief Specialist, COWI 1 The Project is funded by The European Union Agenda 1. Metocean Data Requirements 2. Site
More informationThe Analysis of Uncertainty of Climate Change by Means of SDSM Model Case Study: Kermanshah
World Applied Sciences Journal 23 (1): 1392-1398, 213 ISSN 1818-4952 IDOSI Publications, 213 DOI: 1.5829/idosi.wasj.213.23.1.3152 The Analysis of Uncertainty of Climate Change by Means of SDSM Model Case
More informationConstruction of confidence intervals for extreme rainfall quantiles
Risk Analysis VIII 93 Construction of confidence intervals for extreme rainfall quantiles A. T. Silva 1, M. M. Portela 1, J. Baez & M. Naghettini 3 1 Instituto Superior Técnico, Portugal Universidad Católica
More informationSystem Reliability Analysis Using Tail Modeling
System Reliability Analysis Using Tail Modeling Palaniappn Ramu 1, Nam H. Kim 2 and Raphael T. Haftka 3 University of Florida, Gainesville, Florida, 32611 and Nestor V. Queipo 4 University of Zulia, Maracaibo,
More informationON THE APPROXIMATION OF THE TOTAL AMOUNT OF CLAIMS. T. PENTIKAINEN Helsinki
ON THE APPROXIMATION OF THE TOTAL AMOUNT OF CLAIMS T. PENTIKAINEN Helsinki Several "short cut" methods exist to approximate the total amount of claims (= Z) of an insurance collective. The classical one
More informationTropical Cyclone Oswald Coastal Monitoring
Tropical Cyclone Oswald Coastal Monitoring Kieran Harper and John Maher Coastal Impacts Unit, Queensland Department of Science, Information Technology, Innovation and the Arts (DSITIA). 27 Quinlan St Deagon,
More informationEuropean Geosciences Union General Assembly Vienna, Austria 27 April - 02 May 2014
European Geosciences Union General Assembly Vienna, Austria 27 April - 02 May 2014 Projecting the Current & Future Impact of Storm Surges on Coastal Flood Extent at Pigeon Point, South-West Tobago, through
More informationSelection of Best Fit Probability Distribution for Flood Frequency Analysis in South West Western Australia
Abstract Selection of Best Fit Probability Distribution for Flood Frequency Analysis in South West Western Australia Benjamin P 1 and Ataur Rahman 2 1 Student, Western Sydney University, NSW, Australia
More informationInvestigation of an Automated Approach to Threshold Selection for Generalized Pareto
Investigation of an Automated Approach to Threshold Selection for Generalized Pareto Kate R. Saunders Supervisors: Peter Taylor & David Karoly University of Melbourne April 8, 2015 Outline 1 Extreme Value
More informationProbabilistic Assessment of Coastal Storm Hazards
Resilience of Coastal Infrastructure Conference Hato Rey, PR March 8-9, 2017 Probabilistic Assessment of Coastal Storm Hazards Dr. Norberto C. Nadal-Caraballo Leader, Coastal Hazards Group Team: Victor
More informationUnderstanding Weather and Climate Risk. Matthew Perry Sharing an Uncertain World Conference The Geological Society, 13 July 2017
Understanding Weather and Climate Risk Matthew Perry Sharing an Uncertain World Conference The Geological Society, 13 July 2017 What is risk in a weather and climate context? Hazard: something with the
More informationA quick introduction to Markov chains and Markov chain Monte Carlo (revised version)
A quick introduction to Markov chains and Markov chain Monte Carlo (revised version) Rasmus Waagepetersen Institute of Mathematical Sciences Aalborg University 1 Introduction These notes are intended to
More informationTidal stream atlases Reprinted by PC Maritime with kind permission of Proudman Laboratory
The generation of UK tidal stream atlases from regularly gridded hydrodynamic modelled data by Colin Bell and Lisa Carlin Proudman Oceanographic Laboratory http://www.pol.ac.uk http://www.pcmaritime.co.uk
More informationSpatial and temporal extremes of wildfire sizes in Portugal ( )
International Journal of Wildland Fire 2009, 18, 983 991. doi:10.1071/wf07044_ac Accessory publication Spatial and temporal extremes of wildfire sizes in Portugal (1984 2004) P. de Zea Bermudez A, J. Mendes
More informationPeaks-Over-Threshold Modelling of Environmental Data
U.U.D.M. Project Report 2014:33 Peaks-Over-Threshold Modelling of Environmental Data Esther Bommier Examensarbete i matematik, 30 hp Handledare och examinator: Jesper Rydén September 2014 Department of
More informationAssessing Dependence in Extreme Values
02/09/2016 1 Motivation Motivation 2 Comparison 3 Asymptotic Independence Component-wise Maxima Measures Estimation Limitations 4 Idea Results Motivation Given historical flood levels, how high should
More informationDoes k-th Moment Exist?
Does k-th Moment Exist? Hitomi, K. 1 and Y. Nishiyama 2 1 Kyoto Institute of Technology, Japan 2 Institute of Economic Research, Kyoto University, Japan Email: hitomi@kit.ac.jp Keywords: Existence of moments,
More informationPhysically-based risk assessment of hurricane storm surge in a changing climate
Physically-based risk assessment of hurricane storm surge in a changing climate Ning Lin Princeton University Department of Civil and Environmental Engineering Hurricane Ike 5 Year Workshop Rice University,
More informationBayesian inference. Fredrik Ronquist and Peter Beerli. October 3, 2007
Bayesian inference Fredrik Ronquist and Peter Beerli October 3, 2007 1 Introduction The last few decades has seen a growing interest in Bayesian inference, an alternative approach to statistical inference.
More informationGamma processes for time-dependent reliability of structures
Gamma processes for time-dependent reliability of structures J.M. van Noortwijk & M.J. Kallen HKV Consultants, Lelystad, and Delft University of Technology, Delft, Netherlands M.D. Pandey Department of
More informationWave-height hazard analysis in Eastern Coast of Spain - Bayesian approach using generalized Pareto distribution
Wave-height hazard analysis in Eastern Coast of Spain - Bayesian approach using generalized Pareto distribution J. J. Egozcue, V. Pawlowsky-Glahn, M. I. Ortego To cite this version: J. J. Egozcue, V. Pawlowsky-Glahn,
More informationMozambique. General Climate. UNDP Climate Change Country Profiles. C. McSweeney 1, M. New 1,2 and G. Lizcano 1
UNDP Climate Change Country Profiles Mozambique C. McSweeney 1, M. New 1,2 and G. Lizcano 1 1. School of Geography and Environment, University of Oxford. 2.Tyndall Centre for Climate Change Research http://country-profiles.geog.ox.ac.uk
More informationNOAA Inundation Dashboard
NOAA Inundation Dashboard Audra Luscher & Paul Fanelli NOS Center for Operational Oceanographic Products & Services (CO-OPS) Hampton Roads Sea Level Rise/Flooding Adaptation Forum July 29, 2016 http://tidesandcurrents.noaa.gov/inundationdb/newyork.html
More informationRainfall Analysis in Mumbai using Gumbel s Extreme Value Distribution Model
Rainfall Analysis in Mumbai using Gumbel s Extreme Value Distribution Model S. Sangita Mishra 1, M. Solapure 2, R. Satpathy 3, K. Phadtare 4, S. Pawar 5 1 Assistant Professor, 2,3,4,5 Undergraduate students
More informationLong Term Monitoring of the Mediterranean and Red Sea Levels in Israel
Long Term Monitoring of the Mediterranean and Red Sea Levels in Israel Boris SHIRMAN and Yossi MELZER, Israel Key words: tide-gauge station, mean sea level, geoid, meteorological factors. SUMMARY The Survey
More informationModelling Operational Risk Using Bayesian Inference
Pavel V. Shevchenko Modelling Operational Risk Using Bayesian Inference 4y Springer 1 Operational Risk and Basel II 1 1.1 Introduction to Operational Risk 1 1.2 Defining Operational Risk 4 1.3 Basel II
More informationSea Level Rise and Hurricane Florence storm surge research methodology
Sea Level Rise and Hurricane Florence storm surge research methodology Hurricane Florence storm surge analysis was conducted using a variety of input sources. In order to determine the maximum storm surge
More informationMAS2317/3317: Introduction to Bayesian Statistics
MAS2317/3317: Introduction to Bayesian Statistics Case Study 2: Bayesian Modelling of Extreme Rainfall Data Semester 2, 2014 2015 Motivation Over the last 30 years or so, interest in the use of statistical
More information5.2 Annual maximum sea levels in Venice
Chapter 5 Non stationary extremes 5.1 Introduction In the context of environmental processes, it is common to observe non stationarity for example, due to different seasons having different climate patterns,
More informationShape of the return probability density function and extreme value statistics
Shape of the return probability density function and extreme value statistics 13/09/03 Int. Workshop on Risk and Regulation, Budapest Overview I aim to elucidate a relation between one field of research
More informationSea Level Space Watch: Service Offering
Space for Smarter Government Programme Sea Level Space Watch: Service Offering General Sea Level Space Watch : Operational service developed to support national flood defence planning. Systematically updated
More informationFlood risk management the Dutch approach Patrick Poelmann
Flood risk management the Dutch approach Patrick Poelmann Vice governor Province of Noord Holland Flood risk management Introduction Risk management in 20 th century The Dutch Delta Works Risk management
More informationNew Classes of Multivariate Survival Functions
Xiao Qin 2 Richard L. Smith 2 Ruoen Ren School of Economics and Management Beihang University Beijing, China 2 Department of Statistics and Operations Research University of North Carolina Chapel Hill,
More informationSubject Index. Block maxima, 3 Bootstrap, 45, 67, 149, 176 Box-Cox transformation, 71, 85 Brownian noise, 219
Subject Index Entries in this index are generally sorted with page number as they appear in the text. Page numbers that are marked in bold face indicate that the entry appears in a title or subheading.
More informationExtreme Precipitation: An Application Modeling N-Year Return Levels at the Station Level
Extreme Precipitation: An Application Modeling N-Year Return Levels at the Station Level Presented by: Elizabeth Shamseldin Joint work with: Richard Smith, Doug Nychka, Steve Sain, Dan Cooley Statistics
More informationAccommodating measurement scale uncertainty in extreme value analysis of. northern North Sea storm severity
Introduction Model Analysis Conclusions Accommodating measurement scale uncertainty in extreme value analysis of northern North Sea storm severity Yanyun Wu, David Randell, Daniel Reeve Philip Jonathan,
More informationScarborough Tide Gauge
Tide Gauge Location OS: 504898E 488622N WGS84: Latitude: 54 16' 56.990"N Longitude: 00 23' 25.0279"W Instrument Valeport 740 (Druck Pressure Transducer) Benchmarks Benchmark Description TGBM = 4.18m above
More informationFinal Report V2 November 13, PREPARED BY: Tetra Tech, Inc Powers Ferry Rd. SE, Suite 202 Atlanta, Georgia Phone: (770)
PREPARED BY: Tetra Tech, Inc. 2110 Powers Ferry Rd. SE, Suite 202 Atlanta, Georgia 30339 Phone: (770) 850-0949 Final Report V2 November 13, 2015 PREPARED FOR: Department of the Army Savannah District,
More informationof the 7 stations. In case the number of daily ozone maxima in a month is less than 15, the corresponding monthly mean was not computed, being treated
Spatial Trends and Spatial Extremes in South Korean Ozone Seokhoon Yun University of Suwon, Department of Applied Statistics Suwon, Kyonggi-do 445-74 South Korea syun@mail.suwon.ac.kr Richard L. Smith
More informationExtreme Value Theory and Applications
Extreme Value Theory and Deauville - 04/10/2013 Extreme Value Theory and Introduction Asymptotic behavior of the Sum Extreme (from Latin exter, exterus, being on the outside) : Exceeding the ordinary,
More informationCalculating the Expexted Time to Control the Recruitment When System Fails
Available online at www.scinzer.com Scinzer Journal of Engineering, Vol 2, Issue 1, (2016): 1-5 Calculating the Expexted Time to Control the Recruitment When System Fails Kannadasan K*, Pandiyan P, Vinoth
More information