Verification of DMI wave forecasts 3rd quarter of 2002

Transcription

1 DANISH METEOROLOGICAL INSTITUTE TECHNICAL REPORT - Verification of DMI wave forecasts rd quarter of Jacob Woge Nielsen dmi.dk Copenhagen

2 ISSN -X ISSN - (trykt) (online)

3 Verification of DMI wave forecasts rd quarter of Jacob Woge Nielsen The Danish Meteorological Institute, Copenhagen, Denmark forecast Scatter diagram h, all stations cc=. ervation

4 Contents Introduction with Key Numbers DMI-WAM. Physical model Model set-up Weather model Wave data. Forecasts Observations Error measures Results. Significant wave height Extreme wave height Mean wave period Dominant wave period Mean wave direction Conclusion Appendix. Wave recorder information Observed wave statistics tables Wave height vs. wave period plots Wave forecast statistics tables Significant wave height plots Mean wave period plots Dominant wave period plots

5 Introduction with Key Numbers We analyse the quality of wave forecasts valid for the rd quarter of, produced by DMI-WAM - DMI s operational set-up of the rd generation wave model WAM Cycle. The significant wave height H s is of primary interest, but other wave parameters (period, direction, swell) are examined where the data material is adequate. Standard error measures (bias, rms error,..) are calculated both as a function of forecast range and of wave height. Special statistics are done for the highest waves. Grand averages are calculated as mean values over all stations, over all ranges, and for separate geographical regions. All model results are forecasts. This means that errors in the parameters do not necessarily imply a wave model error, but may reflect errors in the wind forecast. DMI has produced short-range operational wave forecasts since. A pre-operational validation study was carried out in [], a combined wave-wind validation study in [], and a verification pilot study in []. Previous quarterly report(s) are [],[]. Outline: Ch. describes the DMI wave model set-up, ch. lines out the data material, and in ch. we define the statistical error measures used to describe the forecast quality. Ch. presents and discusses the results. Ch. concludes the work. Comprehensive results for each station are found in the Appendix. References and lists of figures/tables are found at the end of the report. For convenience, the Table below shows Key Numbers pertaining to the full model system. Please refer to the Results section (ch. ) for a detailed explanation and discussion. Parameter H s T T p w bias cm.s.s o relative bias -% % % rms error cm.s.s st.dev o scatter index... corr.coeff..... peak bias cm rel. peak bias % Table. Key numbers

6 DMI-WAM DMI has run an operational wave forecasting service for Danish waters since, using the rd generation wave model WAM Cycle (described in [], []) forced by DMI s numerical weather prediction model HIRLAM. In, the geographical model domain includes a large part of the North Atlantic, the North Sea and Baltic Sea, and the Mediterranean. The DMI-WAM model set-up is described in detail below.. Physical model WAM Cycle solves the spectral + ~c rf ~ = S in + S nl + S ds + S cu + S bf where F (f ~x t) is spectral wave energy density, depending on wave frequency, wave direction, position and time; c(f d) is the depth-dependent wave group speed; S in is wind energy input; S nl is non-linear wave-wave interaction; S ds is wave energy dissipation through wave breaking (white capping); S cu is wave-current interaction; S bf is interaction with the sea bed through friction and wave refraction. DMI-WAM still lacks current data (S cu =) and information about sea ice (acts to dampen waves).. Model set-up DMI provides wave forecasts in three geographical domains as shown below: Figure. DMI wave models. North Atlantic, North Sea-Baltic, and Mediterranean model The model open boundaries are chosen as follows. The coarse grid North Atlantic model uses the JONSWAP wind-sea spectrum (see [], []). The fine grid North Sea - Baltic model is nested into the North Atlantic model, and uses time/space interpolated boundary wave spectra calculated by that model. The Mediterranean is treated as a closed basin, assuming no wave energy exchange with the Atlantic or the Black Sea. Please refer to Table for a model set-up summary. The wave forecasting system was coldstarted using developed sea. Subsequent model runs are initialised using the sea state at analysis time, calculated by the previous run.

7 Model North Atlantic North Sea - Baltic Mediterranean Space res. / o / o / o Time step min min min Frequencies Direction resol. o o o Forcing model(s) Hirlam G Hirlam E Hirlam E+G - resolution of. o. o. o /. o Longitudes o W- o E o W- o E o W- o E Latitudes o N- o N o N- o N. o N- o N Open boundaries JONSWAP Nested Closed basin Forecast range h h h Output time step h h h Schedule x daily x daily x daily Table. DMI-WAM set-up. The wave model frequencies range from. Hz to. Hz in % steps. The Mediterranean model patches Hirlam E+G to get maximum resolution. Changes to the model set-up usually require a new coldstart.. Weather model The forcing models are the DMI limited area numerical weather prediction models Hirlam-E and Hirlam-G. Both are currently being used in the DMI weather forecasting service. The G model embeds the E model, but in coarser spatial resolution (km vs. km on a rotated latitude-longitude grid). The wind vector at m height is interpolated linearly in time and space to match the spherical wave model grids. Figure. DMI Hirlam. The outermost box is the G model, the box covering most of Europe is the E model.

8 Wave data The verification data consists of operational DMI-WAM wave forecasts, and wave ervations from a number of fixed positions (buoys and platforms).. Forecasts Wave parameters output from DMI-WAM are shown in Table. H s H sw T T p T sw w sw DMI-WAM output Significant wave height Height of swell Mean wave period Dominant wave period Swell mean period Mean wave direction Swell direction Table. DMI-WAM wave parameters, obtained by a suitable integral of the wave energy spectrum. Tp is the discretized model frequency (inverse) containing the highest energy, picked from a set of predefined values. The forecasts are stored as hourly maps in model resolution. Time series for each station are sampled using nearest model grid point. This is done for each analysis and each parameter. During the rd quarter of, out of scheduled forecasts were produced.. Observations The wave recorder positions are shown in Fig. below. A total of stations that record more or less regularly are selected for verification. The mean data coverage is % (see Fig. for missing data). Comprehensive station information is found in Appendix.. Observed wave statistics and wave height/period diagrams are shown in Appendices. and.. Wave data is obtained from a number of sources, as indicated in table. SMHI, KDI, BSH and NCMR data are kindly provided by each agency in question. NDBC data is retrieved via the GTS. Table shows the number of stations for each wave parameter, and for each of geographical domains. NDBC SMHI KDI NCMR BSH Wave Data providers National Data Buoy Center (UK) Swedish Meteorological Institute Danish Coastal Authority National Center for Marine Research (Greece) Bundesamt für Seeschiffahrt in Hamburg Table. Wave data providers.

9 Figure. Wave recorder locations. Parameter H s T T p w H m Atlantic Scotland-Faroe Irish Sea Br. Channel North Sea Danish West Coast Kattegat - Baltic Sea - Mediterranean - Total Table. Number of wave stations in each domain, and for each wave parameter. Maximum wave height H m is not output by the wave model, and none of the fixed stations record swell. Standard sampling rate is hour. stations (cf. Table ) sample only every hours. Standard sampling accuracy is H s :.m, T :.s, T p :s, w : o. A number of buoys use higher accuracy for some parameters. Two buoys (, ) use low.m H s accuracy.

10 Missing data (cov=.).... almag oland truba arkon athos Figure. Missing data. St.,, and sample every hours, but every hour in storm situations. St.,,, and Athos sample every hours. At these stations, the data coverage is based on hour sampling.

11 Error measures Model errors are calculated using a -d residual matrix, built from all available ervations and forecasts. With the general formula residual = forecast ervation the matrix reads (brackets indicate a dependency) residual(station analysis range) with the number of stations depending on the parameter in question (cf. Table ), analysis every hours, and forecasts ranging from - hours in hour steps. With stations, the size of the residual matrix is about.*. By averaging the residual over all analyses, we get the model bias or mean error: bias(station range) Further linear averaging gives the bias for each forecast range (averaged over all stations), for each station (averaged over the full forecast range), and as a grand average. bias(range) bias(station) BIAS In the same way, the root mean square error rms(station range) is calculated and averaged using the residual squared. For the wave height only, the bias and rms error are also calculated as a function of wave height. The residual is sorted into ervation bins.m wide and averaged for each bin. bias(station bin) rms(station bin) Averaging over all stations gives the model error dependency on wave height, calculated both as an absolute value and as a relative error in %. rms The scatter index si = is obtained by normalising rms with the erved mean value. si may <> be used to intercompare rms errors at stations with large differences in wave climate. Averaging is done as above. Correlation coefficients cc(station block) are calculated using forecast pseudo time series, established by concatenating forecasts in hour range blocks. This gives coefficients valid for each of the range blocks -, -,.., - hours. Range block and station-dependent values (cc cc), and a grand average (CC), are calculated. A special peak bias pbias(station block) is calculated using the most extreme events at each station, allowing for a forecast phase error of a few hours. Peak biases are calculated both as absolute and relative values.

12 Results This section describes wave verification results for the rd quarter of, for significant wave height (H s ), mean and dominant wave period (T T p ), and mean wave direction ( w ). In the sections below, we discuss grand averages and regional averages for each wave parameter in turn. Detailed results for each station are found in the Appendix at the end of the report.. Significant wave height Table shows bias and relative bias, rms error, scatter index and correlation coefficient, averaged over the full forecast range. Fig. shows the short range (-h) scatter diagram. The error estimates are sorted out on geographical regions. Parameter #st bias rms si cc Region cm % cm Atlantic.. Scotland Faroe.. Irish Sea.. Br.Channel.. North Sea.. Danish West Coast.. Kattegat Baltic.. Mediterranean.. All Waters.. Table. Significant wave height results forecast Scatter diagram h, all stations cc=. ervation Figure. Significant wave height: short range (-h) scatter diagram There is a small negative bias and an rms error of. m. Scatter index is low (.) and correlation coefficient high (.). There is some regional spread. Waves are underestimated on average (negative bias) at the Atlantic, Irish Sea, Shetland, North Sea and Mediterranean stations, and overestimated (positive bias) in the British Channel. RMS errors range roughly from a quarter to half a metre. The scatter index (see Fig. ) is well below an acceptance level of. in most regions, but

13 exceeds at some stations in the British Channel, in the Baltic and in the Mediterranean. Scatter index (avg=.), h.... almag oland truba arkon athos Figure. Significant wave height: scatter index. A scatter index of.-. may be considered as a lower limit due to stochastic processes. The error dependency on forecast range and on wave height is shown in Fig.. The bias is independent of forecast range, while rms error, scatter index and correlation coefficient get slightly worse as forecast range increases. The rms error is significant already at analysis time since the model is initialised without any use of the erved sea state. Errors depend strongly on wave height. Small waves have small bias, while higher waves (-m) are underestimated by.m on average. For waves higher than that the data material gets very sparse, please refer to section. on peak errors. The rms error increases with wave height. The relative rms error is -% except for very high or very small waves. Results on significant wave height for each single station is shown in Table, Appendix.. Data sheets are presented in Appendix..

14 bias=. rms=. Average residual error [m] Error distribution [%], all stations r.bias=. r.rms=. Observation distribution [m], all stations Average Expl.Var/Corr.Coeff. [%] cc=. ev=. forecast range [hrs/] bias=. rms=. Error distribution [m], all stations.... Average scatter index [%] Figure. Significant wave height

15 . Extreme wave height The error on the highest waves erved is calculated for the highest events at each station, using the forecast error allowing for a few hours phase displacement. Table shows peak biases for each of the geographical domains, averaged over all forecast ranges. Parameter #st peak bias Region cm % Atlantic Scotland Faroe Irish Sea Br.Channel North Sea Danish West Coast Kattegat Baltic Mediterranean All Waters Table. Extreme wave height results There is a negative peak bias in most regions. On average, the system has a negative peak bias of -cm or -%. The dependency of the peak bias on forecast range is shown in Fig.. The peak bias doubles from -.m at short range to -.m beyond day Peak bias. Average peak error [m] forecast range [hrs/] Figure. Peak wave height errors. Table in Appendix. shows peak biases for each single station.

16 . Mean wave period The mean wave period T is recorded at stations. Grand averages are shown in Table, and a short-range scatter diagram in Fig.. Parameter #st bias rms si cc Region sec. % sec. Danish West Coast.... Kattegat Baltic.... Mediterranean.... All Waters.... Table. Mean wave period results At the Danish West Coast stations T is overestimated by almost %. The reason for this is unresolved. At Baltic and Meditteranean stations T is slightly overestimated. The scatter index at these stations is well below the acceptancy level of.. forecast cc=. Scatter diagram h, all stations ervation Figure. Mean wave period, -h range. Results on mean wave period for each station are shown i Table, Appendix.. Data sheets are presented in Appendix..

17 . Dominant wave period The dominant (or peak) wave period Tp is recorded at stations. Grand averages are shown in Table, short-range scatter diagrams in Fig.. Parameter #st bias rms si cc Region sec. % sec. Atlantic.... Scotland Faroe.... Irish Sea.... Br.Channel.... North Sea.... Danish West Coast.... Mediterranean.... All Waters.... Table. Dominant wave period results Tp errors are large, due partly to low recording and forecasting accuracy, and partly to the nonsmoothness of the series, with Tp shifting abruptly between a high and a low period peak. Even when the wave spectrum is rather well predicted, a small error in the shape of the spectrum may lead to very large Tp errors in situations with a two-peaked spectrum (swell and wind sea). Two examples are found in Figs. and, when a low-frequent swell suddenly dominates a very calm sea. forecast Scatter diagram h, all stations cc=. ervation Figure. Dominant wave period, -h range. Results on dominant wave period for each stations is found in Table, Appendix.. Data sheets for Mediterranean stations only in Appendix..

18 . Mean wave direction The mean wave direction w is recorded at stations. The results are presented in Table and the scatter diagram in Fig.. Parameter #st bias std cc Region deg. deg. Danish West Coast. Kattegat Baltic. Mediterranean. All Waters. Table. Mean wave direction results The mean wave direction predictions fit the ervations well, with almost no bias and a high correlation coefficient. The Athos buoy is an exception in this season, with a rather high directional error. forecast Scatter diagram h, all stations cc=. ervation Figure. Mean wave direction, -h range. Results on the mean wave direction for each station is found in Table, Appendix.. Data sheets are not shown.

19 Conclusion DMI wave forecasts valid for the rd quarter of are verified, using wave data from buoys. A Key Number Table is found in Introduction section (Table ). Main conclusions are: Significant wave height H s and mean wave direction w are well predicted We have some problems predicting wave period T and T p There is a large regional spread in forecast quality The H s error depends on wave height Very high waves are usually underestimated on average The forecast quality decreases slightly with increasing forecast range The significant wave height is recorded at all stations. The error distribution is examined in terms of forecast range, as a function of erved wave height, and for separate geographical regions. The ervation-forecast correlation is high,. on average. This falls off a bit as the prediction range increases. The average bias is small, independent of forecast range, There is some geographical spread, and a pronounced dependency on wave height. Low waves are overpredicted, while the highest waves are underestimated by about %. The average rms error is.m, increasing gradually with forecast range. For medium-sized waves (-m) the rms error is roughly %; larger and smaller waves have larger relative error. The average scatter index SI=.. stations have si>. (sometimes used as an acceptance level), mostly due to low recording accuracy or small average wave height. The highest waves are usually underestimated, with an average % peak bias. The peak error doubles beyond day two. Two types of wave period are recorded; mean and dominant (peak) wave period. The mean wave period is recorded at locations, half of which have a data interpretation problem. At the remaining stations the model overestimates the mean wave period by roughly.s, with a low scatter index, and a correlation coefficient of about.. Dominant wave period predictions are not good. Bias and rms errors range up to several seconds. This could be a data problem; a well predicted wave spectrum does not guarantee a correct dominant wave period in situations with two spectral maxima. Also, most stations sample only with s accuracy and so does the model; this in itself leads to large error measures. Mean wave direction predictions have insignificant bias, a standard deviation of about o, and a high correlation coefficient. Swell parameters are not recorded at any of the fixed positions. A few record maximum wave height but this is not predicted by the wave model.

20 Appendix This Appendix contains a wave recorder station table (below), ervation statistics tables, forecast statistics tables, wave height/period plots, and a plot sheet for each station and each parameter (H s T T p w ), arranged sequentially according to the station table.. Wave recorder information Station ID Agency Region lat. lon. t parameters almag SMHI Baltic.N.E h H s H m T oland SMHI Baltic.N.E h H s H m T truba SMHI Baltic.N.E h H s H m T arkona BSH Baltic.N.E h H s H m T w KDI D. West Coast.N.E h H s H m T T p w KDI D. West Coast.N.E h H s H m T T p w KDI D. West Coast.N.E h H s H m T T p w KDI D. West Coast.N.E h H s H m T T p w NDBC Mediterranean.N.W h H s T p NDBC Mediterranean.N.W h H s T p NDBC Atlantic.N.W h H s T p NDBC North Sea.N.E h H s T p NDBC Atlantic.N.W h H s T p NDBC B.Channel.N.W h H s T p NDBC B.Channel.N.W h H s T p NDBC Atlantic.N.W h H s T p NDBC Atlantic.N.W h H s T p NDBC Atlantic.N.W h H s T p NDBC Atlantic.N.W h H s T p NDBC North Sea.N. h H s T p NDBC North Sea.N. h H s T p NDBC North Sea.N.E h H s T p NDBC Atlantic.N.W h H s T p NDBC Irish Sea.N.W h H s T p NDBC Irish Sea.N.W h H s T p NDBC B.Channel.N.E h H s T p NDBC B.Channel.N. h H s T p NDBC North Sea.N.E h H s NDBC Atlantic.N.W h H s T p NDBC Scotland.N.W h H s T p athos NCMR Mediterranean.N.E h H s T T p w Table. Wave stations. Station name/number, driving agency, position, and wave parameters. SMHI=Swedish Meteorological Institute, BSH=Bundesamt für Seeschiffahrt in Hamburg, KDI=Coastal Authorities (Denmark), NDBC=National Data Buoy Center (UK), NCMR=National Center for Marine Research (Greece). H s =significant wave height, H m =maximum wave height, T=mean wave period, Tp=peak or dominant wave period, w=mean wave direction. t is the sampling rate in hours.

21 . Observed wave statistics tables Station min mean max stdev almag.... oland.... truba.... arkon athos.... Table. Observed wave height. The fraction of missing data is shown in Fig.

22 Station min mean max stdev almag.... oland.... truba.... arkon athos.... Table. Observed mean wave period Station min mean max stdev athos.... Table. Observed dominant wave period

23 . Wave height vs. wave period plots The relation between significant wave height H s and mean wave period T is shown on the diagrams, below, for those stations that record both quantities. At each station, there is a fair linear correlation between H s and T, superimposed by a weak swell component. st. oland, cc=. st. almag, cc=. st. truba, cc=. Mean wave period [s] Mean wave period [s] Mean wave period [s] Significant wave height [m] st. arkon, cc=. Significant wave height [m] st., cc=. Significant wave height [m] st., cc=. Mean wave period [s] Mean wave period [s] Mean wave period [s] Significant wave height [m] st., cc=. Significant wave height [m] st., cc=. Significant wave height [m] st. athos, cc=. Mean wave period [s] Mean wave period [s] Mean wave period [s] Significant wave height [m] Significant wave height [m] Significant wave height [m] Figure. Significant wave height vs. mean wave period

24 . Wave forecast statistics tables Parameter bias rms si cc Station cm % cm almag.. oland.. truba.. arkon athos.. Table. Predicted significant wave height

25 Parameter Obs peak Station m m % almag.. oland.. truba.. arkon athos.. Table. Average of top wave events (peaks) and corresponding mean peak error (peak bias)

26 Parameter bias rms si cc Station sec % sec almag.... oland.... truba.... arkon athos.... Table. Predicted mean wave period Parameter bias rms si cc Station sec % sec athos.... Table. Predicted dominant wave period

27 Parameter bias std cc Station deg deg arkon..... athos. Table. Predicted mean wave direction

28 . Significant wave height plots The following pages show significant wave height error statistics for each station.

29 Obs/pred. wave height [m] at st. almag min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st. almag bias=. rms=. Residual error [m] at st. almag bias=. rms=. Error distribution [m] at st. almag. si=. Scatter index [%] at st. almag.. Error distribution [%] at st. almag r.bias=. r.rms= Figure. Significant wave height: Almagrundet

30 Obs/pred. wave height [m] at st. oland min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st. oland bias=. rms=. Residual error [m] at st. oland bias=. rms=. Error distribution [m] at st. oland. si=. Scatter index [%] at st. oland.. Error distribution [%] at st. oland r.bias=. r.rms= Figure. Significant wave height: Øland

31 Obs/pred. wave height [m] at st. truba min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st. truba bias=. rms=. Residual error [m] at st. truba bias=. rms=. Error distribution [m] at st. truba. si=. Scatter index [%] at st. truba.. Error distribution [%] at st. truba r.bias=. r.rms= Figure. Significant wave height: Trubaduren

32 Obs/pred. wave height [m] at st. arkon min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st. arkon bias=. rms=. Residual error [m] at st. arkon bias=. rms=. Error distribution [m] at st. arkon. si=. Scatter index [%] at st. arkon.. Error distribution [%] at st. arkon r.bias=. r.rms= Figure. Significant wave height: Arkona

33 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

34 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:. The zig-zag curves are caused by irregular h sampling.

35 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:.

36 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

37 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

38 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

39 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

40 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

41 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

42 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:. The sampling accuracy is.m

43 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

44 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:.

45 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:.

46 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

47 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

48 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

49 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:.

50 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

51 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

52 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:. The sampling accuracy is.m

53 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

54 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

55 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

56 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:. The zig-zag curves are caused by irregular h sampling.

57 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:.

58 Obs/pred. wave height [m] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st bias=. rms=. Residual error [m] at st bias=. rms=. Error distribution [m] at st.. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= Figure. Significant wave height:

59 Obs/pred. wave height [m] at st. athos min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [m] at st. athos bias=. rms=. Residual error [m] at st. athos bias=. rms=. Error distribution [m] at st. athos. si=. Scatter index [%] at st. athos.. Error distribution [%] at st. athos r.bias=. r.rms= Figure. Significant wave height: Athos

60 . Mean wave period plots The following pages show mean wave period error statistics for each station. Only stations with reasonable statistics are included.

61 Obs/pred. mean wave period [s] at st. almag min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [s] at st. almag ervation [s] bias=. rms=. Residual error [s] at st. almag... bias=. rms=. Error distribution [s] at st. almag ervation [s]. si=. Scatter index [%] at st. almag.. Error distribution [%] at st. almag r.bias=. r.rms= ervation [s] Figure. Mean wave period: Almagrundet

62 Obs/pred. mean wave period [s] at st. oland min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [s] at st. oland ervation [s] bias=. rms=. Residual error [s] at st. oland... bias=. rms=. Error distribution [s] at st. oland ervation [s]. si=. Scatter index [%] at st. oland.. Error distribution [%] at st. oland r.bias=. r.rms= ervation [s] Figure. Mean wave period: Øland

63 Obs/pred. mean wave period [s] at st. truba min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [s] at st. truba ervation [s] bias=. rms=. Residual error [s] at st. truba... bias=. rms=. Error distribution [s] at st. truba ervation [s]. si=. Scatter index [%] at st. truba.. Error distribution [%] at st. truba r.bias=. r.rms= ervation [s] Figure. Mean wave period: Trubaduren

64 Obs/pred. mean wave period [s] at st. arkon min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [s] at st. arkon ervation [s] bias=. rms=. Residual error [s] at st. arkon... bias=. rms=. Error distribution [s] at st. arkon ervation [s]. si=. Scatter index [%] at st. arkon.. Error distribution [%] at st. arkon r.bias=. r.rms= ervation [s] Figure. Mean wave period: Arkona

65 Obs/pred. mean wave period [s] at st. athos min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [s] at st. athos ervation [s] bias=. rms=. Residual error [s] at st. athos... bias=. rms=. Error distribution [s] at st. athos ervation [s]. si=. Scatter index [%] at st. athos.. Error distribution [%] at st. athos r.bias=. r.rms= ervation [s] Figure. Mean wave period: Athos

66 . Dominant wave period plots The following pages show dominant wave period error statistics for each station. Only stations with reasonable statistics are included.

67 Obs/pred. dominant wave period [s] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [s] at st. ervation [s].. bias=. rms=. Residual error [s] at st.. bias=. rms=. Error distribution [s] at st..... ervation [s]. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= ervation [s] Figure. Dominant wave period:

68 Obs/pred. dominant wave period [s] at st. min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [s] at st. ervation [s].. bias=. rms=. Residual error [s] at st.. bias=. rms=. Error distribution [s] at st..... ervation [s]. si=. Scatter index [%] at st... Error distribution [%] at st. r.bias=. r.rms= ervation [s] Figure. Dominant wave period:

69 Obs/pred. dominant wave period [s] at st. athos min=. prog h max=. mean=. stdev=. / / / / / / Observation distribution [s] at st. athos ervation [s].. bias=. rms=. Residual error [s] at st. athos. bias=. rms=. Error distribution [s] at st. athos.... ervation [s]. si=. Scatter index [%] at st. athos.. Error distribution [%] at st. athos r.bias=. r.rms= ervation [s] Figure. Dominant wave period: Athos

70 References [] G. Komen et al. Dynamics and Modelling of Ocean Waves. Cambridge University Press,. [] The SWAMP group. Ocean Wave Modelling. Plenum Press, New York,. [] H. Günther, S. Hasselmann, and P. Janssen. Wamodel cycle (revised version). Technical Report, Deutches Klimarechnenzentrum,. [] J. W. Nielsen. Verification of DMI wave forecasts: st quarter of. Technical Report -, Danish Meteorological Institute,. [] J. W. Nielsen. Verification of DMI wave forecasts: nd quarter of. Technical Report -, Danish Meteorological Institute,. [] J. W. Nielsen, J. B. Jørgensen, and J. She. Verification of wave forecasts: DMI-WAM nov-dec. Technical Report -, Danish Meteorological Institute,. [] J. She. HIRLAM-WAM quality assessment on winds and waves in the North Sea. Scientific Report -, Danish Meteorological Institute,. [] J. She and J.W. Nielsen. Operational wave forecasts in Baltic and North Sea. Scientific Report -, Danish Meteorological Institute,.

71 List of Tables Key numbers DMI-WAM operational set-up DMI-WAM wave parameters Wave data providers Wave stations in each domain Significant wave height results Extreme wave height results Mean wave period results Dominant wave period results Mean wave direction results Wave stations and locations Observed wave height Observed mean wave period Observed dominant wave period Predicted significant wave height Predicted peak waves and errors Predicted mean wave period Predicted dominant wave period Predicted mean wave direction List of Figures DMI wave models DMI Hirlam Wave recorders Missing data Significant wave height: short range (-h) scatter diagram Significant wave height: scatter index Significant wave height Peak wave height errors Mean wave period Dominant wave period Mean wave direction Significant wave height vs. mean wave period Significant wave height: Almagrundet Significant wave height: Øland Significant wave height: Trubaduren Significant wave height: Arkona Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height:

72 Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Significant wave height: Athos Mean wave period: Almagrundet Mean wave period: Øland Mean wave period: Trubaduren Mean wave period: Arkona Mean wave period: Athos Dominant wave period: Dominant wave period: Dominant wave period: Athos