Actions

Harmonic Analysis of Current Velocity: Difference between revisions

From BAWiki

imported>Lang Guenther
(first text version, incomplete)
 
(The LinkTitles extension automatically added links to existing pages (<a target="_blank" rel="nofollow noreferrer noopener" class="external free" href="https://github.com/bovender/LinkTitles">https://github.com/bovender/LinkTitles</a>).)
 
(8 intermediate revisions by 3 users not shown)
Line 3: Line 3:


=Motivation=
=Motivation=
For almost any location within an [[estuary]] the current velocity can be decomposed into a periodic part driven by the tides (astronomical and [[shallow water]] tides) as well as a non-periodic part driven by river [[runoff]], meteorological effects as well as local topographic features.
[[File:Tidal_current_ellipse_definition.png|200px|thumb|right|defintion of tidal ellipse]]
For each tidal constituent the tip of the current velocity vector describes an ellipse, the tidal ellipse (see accompanying figure). Within each period of the respective tidal constituent a full rotation is executed. Characteristic parameters for a tidal ellipse are semi major axis, semi minor axis, inclination of tidal ellipse as well as its eccentricity. In dependence on eccentricity the ellipse degenerates either to a circle (eccentricity +1.0 or -1.0) or a straight line (eccentricity 0.0). The beforementioned parameters are well suited to describe the tidal influence on current velocity.


=Basic Concepts=
=Basic Concepts=


Basic informations on  ''Tidal Harmonic Analysis'' are presented in [[Harmonic Analysis of Water Level]]. The methods used to compute various ellipse parameters using program [[NCANALYSE]] are described in more detail in publication [[File:zhigang_xu_2000_ellipse_parameters_conversion.pdf]].
Basic informations on  ''Tidal Harmonic Analysis'' are presented in [[Harmonic Analysis of Water Level]].  
The methods used to compute various ellipse parameters using program [[NCANALYSE]] are described in more detail in publication  
[[Media:zhigang_xu_2000_ellipse_parameters_conversion.pdf|zhigang_xu_2000_ellipse_parameters_conversion.pdf]].
 
The computational procedure is as follows (number of equation with reference to Zhigang Xu (2000) given in brackets):
# Compute amplitude and phase for u-component (x- or eastward) of current velocity (Eq. 1).
# Compute amplitude and phase for v-component (y- or northward) of current velocity  (Eq 2).
# Determine equivalent complex formulas for u and v (Eq. 9 und Eq. 10).
# Derive complex formulas for (mathematical) positive as well as negative rotating circular components (Eq. 7 and Eq. 8). Superposition of two components with different sense of orientation results in a tidal ellipse.
# From the lastly mentioned complex formulas magnitude an phase for positive and negative rotating components can be computed.
# With these results it is possible to compute semi major axis and semi minor axis (Eq. 14 and Eq. 19), axes inclination (Eq. 15 and Eq. 20) as well as eccentricity (Eq. 21).
The aforementioned calculations must be performed for any tidal constituent.


=Definition of Characteristic Numbers=
=Definition of Characteristic Numbers=


====Residual Tidal Constituent Velocity (Vres)====
The Least Square Fit Method (LSQF) is used to compute amplitude and phase for the u- and v-component of current velocity.
Residual value of current velocity (computed from harmonic analysis) within the period of data analysis.<br />
For this purpose program [[NCANALYSE]] makes use of methods available within the software package LAPACK.
The residual is computed for permanently wet locations only.<br />
 
<!-- [[Example figures: residual current velocity (computed from harmonic analysis)]]. <br /> -->
All quantities are computed for permanently wet areas only.
 
For each tidal constituent the following quantities are determined.
 
==Residual Current Velocity (Vres)==
 
Residual current velocity (computed from tidal constituent analysis) within the [[period of data analysis]].
This quantity describes the non-tidal (or non-periodc) contribution. It mainly depends on river flow,
meteorology and local topographic details, but is also dependent to some extent on the number of tidal
constituents taken into account.
<!-- [[Example figure: Residual current velocity (from tidal constituent analysis)]]. <br /> -->
 
==Amplitude and Phase of Tidal Constituents==
 
The subsequent quantities are essential to derive the various characteristic parameters for tidal ellipses
introduced below. Standard error for each of the four quantities is also computed. The quantities,
although important, are not well suited to be used for visualization.
 
===Amplitude of u-Component===
 
Amplitude of current velocity u-component (x- or eastward) of a tidal constituent.
<!-- visualization is not very meaningful for this type of quantity -->
 
===Amplitude of v-Component===
 
Amplitude of current velocity v-component (y- or northward) of a tidal constituent.
<!-- visualization is not very meaningful for this type of quantity -->
 
===Phase of u-Component===
 
Phase (in degree) of current velocity u-component (x- or eastward) of a tidal constituent.
<!-- visualization is not very meaningful for this type of quantity -->
 
===Phase of v-Component===
 
Phase (in degree) of current velocity v-component (y- or northward) of a tidal constituent.
<!-- visualization is not very meaningful for this type of quantity -->
 
==Characteristic Parameters for Tidal Ellipses==
 
From the above mentioned quantities (amplitudes and phases) all subsequent characteristic parameters for tidal ellipses can be derived.
 
===Semi Major Axis of Tidal Ellipse===
 
Semi major axis of a tidal constituent. Gives maximum current velocity within one period for the tidal constituent.
<!-- [[Example figure: Semi major axis for tidal constituent velocity]]. <br /> -->
 
===Semi Minor Axis of Tidal Ellipse===
 
Semi minor axis of a tidal constituent. Gives minimum current velocity within one period for the tidal constituent.
<!-- [[Example figure: Semi minor axis for tidal constituent velocity]]. <br /> -->
 
===Inclination (direction) of Major Axis===
 
Gives the (mathematically positive) direction (-180.0 to 180.0 degree) of the semi major axis for a tidal constituent velocity.
Is equivalent to the orientation of the maximum current induced by the tidal constituent.
<!-- visualization is not very meaningful for this type of quantity -->
 
===Inclination (direction) of Minor Axis===
 
Gives the (mathematically positive) direction (-180.0 to 180.0 degree) of the semi minor axis for a tidal constituent velocity.
Is equivalent to the orientation of the minimum current induced by the tidal constituent.
<!-- visualization is not very meaningful for this type of quantity -->
 
===Eccentricity of Tidal Ellipse===
 
Eccentricity is defined as the ratio between semi minor and semi major axis of a tidal ellipse. Special cases are:
* 0.0 : rectilinear current, e. g. due to topographic guidance;
* 1.0 : circular current velocity in mathematical positive direction (counter clockwise);
* -1.0 : circular current velocity in mathematical negative direction (clockwise).
<!-- [[Example figure: Eccentricity of tidal constituent velocity ellipse]]. <br /> -->
 
===Phase of maximum Current Velocity===
 
Phase (in dgree) of maximum current velocity for a tidal constituent. Restricted range -180 to 180 degree.
<!-- [[Example figure: Phase of maximum current velocity]]. <br /> -->
 
===Phase of minimum Current Velocity===
 
Phase (in dgree) of minimum current velocity for a tidal constituent. Restricted range -180 to 180 degree.
<!-- [[Example figure: Phase of maximum current velocity]]. <br /> -->
 
==Flow Duration Asymmetry (Skewness FDA)==
 
This parameter describes differences in the flow velocities around [[flood]] slack compared to those around ebb slack. Positive values indicate smaller velocities around [[flood]] slack. Please notice that a residual current is neglected in the formulas used.
 
===Skewness FDA===
Computation of skewness FDA (sum over all contributions of tidal constituents) according to [https://www.researchgate.net/publication/256682174_The_contribution_to_tidal_asymmetry_by_different_combinations_of_tidal_constituents Song (2011), Eq. 16] from amplitude, phase and frequency of tidal constituents.<br />
<!-- [[Example figure: Skewness FDA from all tidal constituents]]. <br /> -->


==Amplitude of Tidal Constituent Velocity==
===Skewness FDA for all pairs and triplets===
Computation of skewness FDA for all pairs and triplets of suitable tidal constituents according to [https://www.researchgate.net/publication/256682174_The_contribution_to_tidal_asymmetry_by_different_combinations_of_tidal_constituents Song (2011), Eq. 20 and Eq. 21] from amplitude, phase and frequency of contributing tidal constituents.<br />
<!-- [[Example figure: Skewness FDA for all pairs and triplets]]. <br /> -->


==Phase of Tidal Constituent Velocity==


==Error in Amplitude of Tidal Constituent Velocity==
==Flow Velocity Asymmetry (Skewness FVA)==


==Error in Phase of Tidal Constituent Velocity==
This parameter describes differences between maximum [[flood]] current and ebb current velocities. Positive values indicate larger [[flood]] current velocities. Please notice that a residual current is taken into account in the formulas used.


==Semi Major Axis of Tidal Constituent Ellipse==
===Skewness FVA===
Computation of skewness FVA (sum over all contributions of tidal constituents and a residual current) according to [https://www.researchgate.net/publication/299342942_Tidal_asymmetry_in_a_funnel-shaped_estuary_with_mixed_semidiurnal_tides Gong (2016), Eq. 24] from amplitude and phase of the tidal constituents as well as a residual current.<br />
<!-- [[Example figure: Skewenss FVA from all tidal constituents and a residual current]]. <br /> -->


==Semi Minor Axis of Tidal Constituent Ellipse==
===Skewness FVA for all pairs and triplets===
Computation of skewness FVA for all pairs and triplets formed out of suitable tidal constituents according to [https://www.researchgate.net/publication/299342942_Tidal_asymmetry_in_a_funnel-shaped_estuary_with_mixed_semidiurnal_tides Gong (2016), Eq. 24] from amplitude and phase of the tidal constituents as well as a residual current.<br />
<!-- [[Example figure: Skewness FVA for all pairs and triplets]]. <br /> -->


==Eccentricity of Tidal Constituent Ellipse==
===Skewness FVA from residual current===
Computation of skewness resulting from a residual current according to [https://www.researchgate.net/publication/299342942_Tidal_asymmetry_in_a_funnel-shaped_estuary_with_mixed_semidiurnal_tides Gong (2016), Eq. 24] from amplitude and phase of the tidal constituents as well as a residual current.<br />
<!-- [[Example figure: Skewness FVA for a residual current]]. <br /> -->


----
----

Latest revision as of 10:21, 21 October 2022


Motivation

For almost any location within an estuary the current velocity can be decomposed into a periodic part driven by the tides (astronomical and shallow water tides) as well as a non-periodic part driven by river runoff, meteorological effects as well as local topographic features.

defintion of tidal ellipse

For each tidal constituent the tip of the current velocity vector describes an ellipse, the tidal ellipse (see accompanying figure). Within each period of the respective tidal constituent a full rotation is executed. Characteristic parameters for a tidal ellipse are semi major axis, semi minor axis, inclination of tidal ellipse as well as its eccentricity. In dependence on eccentricity the ellipse degenerates either to a circle (eccentricity +1.0 or -1.0) or a straight line (eccentricity 0.0). The beforementioned parameters are well suited to describe the tidal influence on current velocity.

Basic Concepts

Basic informations on Tidal Harmonic Analysis are presented in Harmonic Analysis of Water Level. The methods used to compute various ellipse parameters using program NCANALYSE are described in more detail in publication zhigang_xu_2000_ellipse_parameters_conversion.pdf.

The computational procedure is as follows (number of equation with reference to Zhigang Xu (2000) given in brackets):

  1. Compute amplitude and phase for u-component (x- or eastward) of current velocity (Eq. 1).
  2. Compute amplitude and phase for v-component (y- or northward) of current velocity (Eq 2).
  3. Determine equivalent complex formulas for u and v (Eq. 9 und Eq. 10).
  4. Derive complex formulas for (mathematical) positive as well as negative rotating circular components (Eq. 7 and Eq. 8). Superposition of two components with different sense of orientation results in a tidal ellipse.
  5. From the lastly mentioned complex formulas magnitude an phase for positive and negative rotating components can be computed.
  6. With these results it is possible to compute semi major axis and semi minor axis (Eq. 14 and Eq. 19), axes inclination (Eq. 15 and Eq. 20) as well as eccentricity (Eq. 21).

The aforementioned calculations must be performed for any tidal constituent.

Definition of Characteristic Numbers

The Least Square Fit Method (LSQF) is used to compute amplitude and phase for the u- and v-component of current velocity. For this purpose program NCANALYSE makes use of methods available within the software package LAPACK.

All quantities are computed for permanently wet areas only.

For each tidal constituent the following quantities are determined.

Residual Current Velocity (Vres)

Residual current velocity (computed from tidal constituent analysis) within the period of data analysis. This quantity describes the non-tidal (or non-periodc) contribution. It mainly depends on river flow, meteorology and local topographic details, but is also dependent to some extent on the number of tidal constituents taken into account.

Amplitude and Phase of Tidal Constituents

The subsequent quantities are essential to derive the various characteristic parameters for tidal ellipses introduced below. Standard error for each of the four quantities is also computed. The quantities, although important, are not well suited to be used for visualization.

Amplitude of u-Component

Amplitude of current velocity u-component (x- or eastward) of a tidal constituent.

Amplitude of v-Component

Amplitude of current velocity v-component (y- or northward) of a tidal constituent.

Phase of u-Component

Phase (in degree) of current velocity u-component (x- or eastward) of a tidal constituent.

Phase of v-Component

Phase (in degree) of current velocity v-component (y- or northward) of a tidal constituent.

Characteristic Parameters for Tidal Ellipses

From the above mentioned quantities (amplitudes and phases) all subsequent characteristic parameters for tidal ellipses can be derived.

Semi Major Axis of Tidal Ellipse

Semi major axis of a tidal constituent. Gives maximum current velocity within one period for the tidal constituent.

Semi Minor Axis of Tidal Ellipse

Semi minor axis of a tidal constituent. Gives minimum current velocity within one period for the tidal constituent.

Inclination (direction) of Major Axis

Gives the (mathematically positive) direction (-180.0 to 180.0 degree) of the semi major axis for a tidal constituent velocity. Is equivalent to the orientation of the maximum current induced by the tidal constituent.

Inclination (direction) of Minor Axis

Gives the (mathematically positive) direction (-180.0 to 180.0 degree) of the semi minor axis for a tidal constituent velocity. Is equivalent to the orientation of the minimum current induced by the tidal constituent.

Eccentricity of Tidal Ellipse

Eccentricity is defined as the ratio between semi minor and semi major axis of a tidal ellipse. Special cases are:

  • 0.0 : rectilinear current, e. g. due to topographic guidance;
  • 1.0 : circular current velocity in mathematical positive direction (counter clockwise);
  • -1.0 : circular current velocity in mathematical negative direction (clockwise).

Phase of maximum Current Velocity

Phase (in dgree) of maximum current velocity for a tidal constituent. Restricted range -180 to 180 degree.

Phase of minimum Current Velocity

Phase (in dgree) of minimum current velocity for a tidal constituent. Restricted range -180 to 180 degree.

Flow Duration Asymmetry (Skewness FDA)

This parameter describes differences in the flow velocities around flood slack compared to those around ebb slack. Positive values indicate smaller velocities around flood slack. Please notice that a residual current is neglected in the formulas used.

Skewness FDA

Computation of skewness FDA (sum over all contributions of tidal constituents) according to Song (2011), Eq. 16 from amplitude, phase and frequency of tidal constituents.

Skewness FDA for all pairs and triplets

Computation of skewness FDA for all pairs and triplets of suitable tidal constituents according to Song (2011), Eq. 20 and Eq. 21 from amplitude, phase and frequency of contributing tidal constituents.


Flow Velocity Asymmetry (Skewness FVA)

This parameter describes differences between maximum flood current and ebb current velocities. Positive values indicate larger flood current velocities. Please notice that a residual current is taken into account in the formulas used.

Skewness FVA

Computation of skewness FVA (sum over all contributions of tidal constituents and a residual current) according to Gong (2016), Eq. 24 from amplitude and phase of the tidal constituents as well as a residual current.

Skewness FVA for all pairs and triplets

Computation of skewness FVA for all pairs and triplets formed out of suitable tidal constituents according to Gong (2016), Eq. 24 from amplitude and phase of the tidal constituents as well as a residual current.

Skewness FVA from residual current

Computation of skewness resulting from a residual current according to Gong (2016), Eq. 24 from amplitude and phase of the tidal constituents as well as a residual current.


back to Analysis of Calculated Results


Overview