Actions

EVENTFILTER: Difference between revisions

From BAWiki

imported>BAWiki 3
mNo edit summary
mNo edit summary
 
(11 intermediate revisions by 2 users not shown)
Line 10: Line 10:
filter considering fresh water flow<br />
filter considering fresh water flow<br />
desired and actual value<br />
desired and actual value<br />
|shortdescription=In order to analyse time series of high water or low water events you usually calculate their annual mean values. These values are influenced by meteorologic processes e.g. wind or fresh water flow in an estuary. This makes it difficult to compare the mean values of several years. The program eventFilter reduces the meteorologic influence by filtering the events and generates means you can better compare with each other.
|shortdescription=In order to analyse time series of high water or low water events you usually calculate their annual mean values. These values are influenced by meteorologic processes e.g. wind or fresh water flow in an estuary. This makes it difficult to compare the mean values of several years. The program EVENTFILTER reduces the meteorologic influence by filtering the events and generates means you can better compare with each other.
The input values are filtered by a [[#Coupled band-pass|coupled band-pass]] and / or under consideration of the [[#Filter considering fresh water flow|fresh water flow]].<br />  
The input values are filtered by a [[#Coupled band-pass|coupled band-pass]] and / or under consideration of the [[#Filter considering fresh water flow|fresh water flow]].<br />  
|inputfiles=
|inputfiles=
Line 23: Line 23:


|outputfiles=
|outputfiles=
# file(s) containing a table with the statistics of a single station (type [[EVENTSTATISTIC.DAT|eventstatistic.dat]]  (ASCII) or [[EVENTSTATISTIC.EXCEL|eventstatistic.excel]]
# file(s) containing a table with the statistics of a single station (type [http://www.baw.de/methoden_en/index.php5/EVENTSTATISTIC.DAT_and_EVENTSTATISTIC.EXCEL eventstatistic.dat (ASCII) or eventstatistic.excel (excel compatibly)])
#: (excel compatibly))
# informative printer file (eventFilter.sdr)
# informative printer file (eventFilter.sdr)
# (optional) trace of program execution (filetype eventFilter.trc)<br /><br />  
# (optional) trace of program execution (filetype eventFilter.trc)<br /><br />  
Line 31: Line 30:


====Preparations====
====Preparations====
The reference station should be situated in an area with few meteorological influence on high water and low water. In the German Bight Helgoland would be a suited place. At first eventFilter has to identify a reference event as high water or low water. You get the best results by analyzing the synthetic values. EventFilter compares date and time of a measured event with date and time of a neighboured synthetic event and identifies the actual high and low water events.
The reference station should be situated in an area with few meteorological influence on high water and low water. In the German Bight Helgoland would be a suited place. At first EVENTFILTER has to identify a reference event as high water or low water. You get the best results by analyzing the synthetic values. EVENTFILTER compares date and time of a measured event with date and time of a neighboured synthetic event and identifies the actual high and low water events.
(All events in the following pictures are linear connected and therefore not realistically represented.)
(All events in the following pictures are linear connected and therefore not realistically represented.)


[[File:Soll_ist_small.png‎]]
[http://www.baw.de/downloads/wasserbau/mathematische_verfahren/programmkennbl_de/pdf/helgoland_soll_ist.pdf Example: Helgoland desired and actual value]


 
The picture above shows that the identification works well even in the case of a storm flood. The water level differs between the synthetic red curve and the measured green curve up to 2 meters, whereas the time always differs less than an hour.
desired and actual value at the reference The picture above shows that the identification works well even in the case of a storm flood. The water level differs between the synthetic red curve and the measured green curve up to 2 meters, whereas the time always differs less than an hour.


The second way to identify high and low water is to use the mean flood period at the reference station. It should only be used if there are no synthetic values available.
The second way to identify high and low water is to use the mean flood period at the reference station. It should only be used if there are no synthetic values available.
Line 46: Line 44:
The algorithm for nipp tides works analogously.
The algorithm for nipp tides works analogously.


[http://www.baw.de/downloads/wasserbau/mathematische_verfahren/programmkennbl_de/pdf/helgoland_bandpass.pdf Example: Coupled band-pass Helgoland]


bandpass at the reference The diagram shows the described case. Thw1 has been marked by the filter though it has a high value. The tide calender confirms the supposition that it takes place during a spring tide period. Thw2 and Thw3 are marked because they are inside the limits. Thw4 is without astronomical cause to low and has not been marked.
The diagram shows the described case. Thw1 has been marked by the filter though it has a high value. The tide calender confirms the supposition that it takes place during a spring tide period. Thw2 and Thw3 are marked because they are inside the limits. Thw4 is without astronomical cause to low and has not been marked.
Tnw1 has been marked because Thw1 is to high. Tnw2 to Tnw4 are to low.
Tnw1 has been marked because Thw1 is to high. Tnw2 to Tnw4 are to low.


Line 59: Line 58:


The water level in the upstream region of an estuary is highly influenced by the amount of fresh water flow. But in other areas the water level is influenced by extremly high or low amounts, too. Therefore water flow measurements between a lower and upper limit are marked. The corresponding high and low water events are marked as well and their mean values get calculated.
The water level in the upstream region of an estuary is highly influenced by the amount of fresh water flow. But in other areas the water level is influenced by extremly high or low amounts, too. Therefore water flow measurements between a lower and upper limit are marked. The corresponding high and low water events are marked as well and their mean values get calculated.
fresh water flow and high / low water Example: The daily mean fresh water flow at Neu Darchau in the Elbe estuary is shown in the upper part of the picture above. The values 200m**3/s and 800m**3/s limit the filter. E2 is the first event that lays inside these boundaries and that gets marked.
 
The lower diagram shows the high and low water events at Hamburg St. Pauli 87km downstream from Neu Darchau. The events Tnw3, Thw3, Tnw4 and Thw4 are marked because they take place within 24 h after the marked event E2.
[http://www.baw.de/downloads/wasserbau/mathematische_verfahren/programmkennbl_de/pdf/stpauli_oberwasser.pdf Example: fresh water flow and high/low water]
 
 
Example: The daily mean fresh water flow at Neu Darchau in the Elbe estuary is shown in the upper part of the picture above. The values 200m**3/s and 800m**3/s limit the filter. E2 is the first event that lays inside these boundaries and that gets marked.
The lower diagram shows the high and low water events at Hamburg St. Pauli 87 km downstream from Neu Darchau. The events Tnw3, Thw3, Tnw4 and Thw4 are marked because they take place within 24 h after the marked event E2.


====Combination of both filters====
====Combination of both filters====
If you combine both filters the data is at first processed by the coupled band-pass and then by the fresh water flow filter. The combination works similar to a series connection. As a consequence the results in [[EVENTSTATISTIC.DAT|eventstatistic.dat]] and [[EVENTSTATISTIC.EXCEL|eventstatistic.excel]] named '*_Q' are based on events that have been marked in '''both''' filters.  
If you combine both filters the data is at first processed by the coupled band-pass and then by the fresh water flow filter. The combination works similar to a series connection. As a consequence the results in [[EVENTSTATISTIC.DAT and EVENTSTATISTIC.EXCEL|eventstatistic.dat/.excel]] named '*_Q' are based on events that have been marked in '''both''' filters.  


|preprocessor=-
|preprocessor=-
Line 69: Line 72:
|language=Fortran90
|language=Fortran90
|add_software= -
|add_software= -
|contact_original=[mailto:peter.schade@baw.de P. Schade]
|contact_original=P. Schade
|contact_maintenance=[mailto:peter.schade@baw.de P. Schade]
|contact_maintenance=[mailto:pre.proghome@baw.de working group PRE]
|documentation=$PROGHOME/examples/Frqwf/  
|documentation=$PROGHOME/examples/Frqwf/  
}}
}}

Latest revision as of 12:52, 10 October 2022

Basic Information

Name of Program

EVENTFILTER

Version-Date

March 2001

Description-Date

March 2001

Catchwords

time series of high water and low water events
statistical analysis
coupled band-pass filter
filter considering fresh water flow
desired and actual value

Short Description of Functionality

In order to analyse time series of high water or low water events you usually calculate their annual mean values. These values are influenced by meteorologic processes e.g. wind or fresh water flow in an estuary. This makes it difficult to compare the mean values of several years. The program EVENTFILTER reduces the meteorologic influence by filtering the events and generates means you can better compare with each other. The input values are filtered by a coupled band-pass and / or under consideration of the fresh water flow.

Input-Files

  1. general input data (filetype eventfilter.dat)
    The reference files 2. and 3. have to contain high water events as well as low water events.
  2. time series of measured actual high water and low water events at a reference station (file of type boewrt.dat)
  3. (optional) time series of synthetic, desired values of high water and low water at a reference, generated for example by the
    harmonic process (file of type boewrt.dat)
  4. measured time series of high water and low water events at further stations e.g. along an estuary (files of type boewrt.dat)
  5. (optional) time series of fresh water flow into an estuary (file of type boewrt.dat), needed for the filter considering fresh
    water flow. The timestep has to be 24 h.

Output-Files

  1. file(s) containing a table with the statistics of a single station (type eventstatistic.dat (ASCII) or eventstatistic.excel (excel compatibly))
  2. informative printer file (eventFilter.sdr)
  3. (optional) trace of program execution (filetype eventFilter.trc)

Methodology

Preparations

The reference station should be situated in an area with few meteorological influence on high water and low water. In the German Bight Helgoland would be a suited place. At first EVENTFILTER has to identify a reference event as high water or low water. You get the best results by analyzing the synthetic values. EVENTFILTER compares date and time of a measured event with date and time of a neighboured synthetic event and identifies the actual high and low water events. (All events in the following pictures are linear connected and therefore not realistically represented.)

Example: Helgoland desired and actual value

The picture above shows that the identification works well even in the case of a storm flood. The water level differs between the synthetic red curve and the measured green curve up to 2 meters, whereas the time always differs less than an hour.

The second way to identify high and low water is to use the mean flood period at the reference station. It should only be used if there are no synthetic values available.

Coupled band-pass

The coupled band-pass processes the measured high and low water at the reference station. At first all events inside a scope around the mean value are marked. Example high water: all high water events higher than mean high water - 0,1m and lower than mean high water + 0,1m are marked. The events that have not been marked are examined. In the example above a high water higher than the upper high water limit is marked when the low water before or after is lower than the lower low water limit. This is typical for a spring tide. This phenomenon is caused by astronomical forces and the events should be considered. The algorithm for nipp tides works analogously.

Example: Coupled band-pass Helgoland

The diagram shows the described case. Thw1 has been marked by the filter though it has a high value. The tide calender confirms the supposition that it takes place during a spring tide period. Thw2 and Thw3 are marked because they are inside the limits. Thw4 is without astronomical cause to low and has not been marked. Tnw1 has been marked because Thw1 is to high. Tnw2 to Tnw4 are to low.

Next mean values of the marked events are evaluated and the upper and lower limits are adapted to them.

The process of marking and evaluating is so often repeated until both mean values converge.

Then eventFilter identifies high and low water at the other stations by using the known event runtime between reference and station. The events corresponding the marked ones at the reference are marked as well and last not least their mean values get evaluated.

Filter considering fresh water flow

The water level in the upstream region of an estuary is highly influenced by the amount of fresh water flow. But in other areas the water level is influenced by extremly high or low amounts, too. Therefore water flow measurements between a lower and upper limit are marked. The corresponding high and low water events are marked as well and their mean values get calculated.

Example: fresh water flow and high/low water


Example: The daily mean fresh water flow at Neu Darchau in the Elbe estuary is shown in the upper part of the picture above. The values 200m**3/s and 800m**3/s limit the filter. E2 is the first event that lays inside these boundaries and that gets marked. The lower diagram shows the high and low water events at Hamburg St. Pauli 87 km downstream from Neu Darchau. The events Tnw3, Thw3, Tnw4 and Thw4 are marked because they take place within 24 h after the marked event E2.

Combination of both filters

If you combine both filters the data is at first processed by the coupled band-pass and then by the fresh water flow filter. The combination works similar to a series connection. As a consequence the results in eventstatistic.dat/.excel named '*_Q' are based on events that have been marked in both filters.

Program(s) to run before this Program

-

Program(s) to run after this Program

FFT, MESKOR

Additional Information

Language

Fortran90

Additional software

-

Original Version

P. Schade

Maintenance

working group PRE

Documentation/Literature

$PROGHOME/examples/Frqwf/


back to Program Descriptions


Overview