Need help?
Call now 0207 118 0808

GET PRICE NOW

Writer's Profile
James Barton

Specialised Subjects

Civil Engineering, Engineering, Project Management

I have a Bachelor’s degree in Civil Engineering and recently completed a Master’s degree in Marine Civil Engineering from a UK University. In the course of studying for both my Bachelor’s and Master’s degrees, I have had the opportunity to write numerous academic papers on a wide range of subjects, ranging from science to management. For 12 months, I worked as a consultant engineer in the construction of roads, dams, irrigation and drainage schemes for developing countries around the world. I am planning to do a Ph.D. in Engineering but at the moment I am engaged in a graduate training scheme. I have always loved engineering activities.

The use of the harmonic analysis method for analysing tidal levels.

INTRODUCTION
Equilibrium Theory of the tides was the first attempt to describe tidal phenomena. Significant tide-producing forces on earth are caused by the relative configuration and motions of the earth-moon-sun system. Harmonic analysis is another way by which tides can be analysed. Since the tide is composed of various constituents that interact at the same time with the river flow, the resulting records of water levels in the river will show a rather complicated pattern. Harmonic analysis of the hourly water level of a station could be used to record individual interaction of each constituent of the tide in each month. The characteristics of the individual constituents of the tide can then be correlated with the river flow.  There are nine constituents of the tide but the first four are the main ones and are analysed, based on their characteristics which are:

Table 1: Major tidal constituents

No.

SYMBOL

CONSTITUTION

PERIOD, T (HOURS)

FREQUENCY, ω (RAD/HOUR)

1

M2

Lunar Semi-diurnal

12.42

0.5059

2

S2

Solar Semi-diurnal

12.00

0.5236

3

K1

Luni-Solar diurnal

23.93

0.2626

4

O1

Lunar declinational diurnal

25.82

0.2433

5

N2

Larger Lunar elliptic Semi-diurnal

12.66

0.4963

6

M4

First over-tide of M2

6.21

1.0118

7

S4

First over-tide of S2

6.00

1.0472

8

MS4

Compound Luni-solar

6.10

1.0300

9

M6

Second over-tide of M2

4.14

1.5177

The harmonic analysis of tides is developed from the following two parts: the resultant tide at any location is composed of a finite number of constituents, each with its own periodicity, phase angle and amplitude and also the constituents are each simple harmonics in time and are mutually independent. i.e.
ɳ=ao+1………………………..1
Where, ɳ = tidal level, ao= amplitude at mean water level, ai= amplitude at i-constituent, 2 = phase angle at i- constituent, t= time.

AIMS
The exercise is designed to illustrate how the measured water levels from an estuary can be analysed to produce harmonic tidal constants, how the tidal levels can be predicted using harmonic tidal constants and how to use dominate constants (M2, S2, K1 & O1) to identify a tidal type that can affect transport and engineering operations in an estuary.

HARMONIC ANALYSIS THEORY

Harmonic Analysis is the summation of the amplitude at water level and all the harmonic functions (M2, S2, K2, O1, N2, M4, S4, MS4, M6).
i.e.
ɳ=ao+

Where ωi = and φi= phase angle

n= number of samples
i= 1……9

Where ao can be calculated from
ɳ=ao+
We have
= dt +  [t-)dt] As
 [t-)dt] =0

So,
 = ……………………………………………………….2

While ai (amplitude of constituents) can be decided by the following derivation:

ɳ=ao+ (Multiply through by sin(t))
ɳ sin(t) =  sin(t) +  sin(t)

Integrate and take the limit of the above

dt=dt+=+]dt
=  dt + dt
Using the features gives,

 dt =   

Similarly we have:

 dt =   
Or for   we have,

=

 can be decided as
= ………………………..….3

Or simply,
= 2

Phase angle φi is derived from:


We have,
= (  )……………………………………………………4
=  ()

The tide is classified using the classification criteria given below:

F factor: F= ……………………………………………………………………5
where the symbols of the constituents indicate their respective amplitudes.

Table 2: Classification of tides

F

Category

< 0.25

Semidiurnal

0.25-3

Mixed (Mainly Semidiurnal 0.25 < F < 1.5)

Mixed (Mainly Diurnal 1.5 < F < 3)

>3

Diurnal

RESULTS AND DISCUSSION

SOLUTION TO QUESTION 1

Based on the given template, reconstructing the water level changes relative to Chart Datum over a period of at least half a year at quarter-hourly intervals using the harmonic tidal data from Liverpool (Gladstone Dock) was done by calculating for the harmonic functions of each constituent (M2, S2, K1, O1, N2, M4, S4, MS4 and M6) of the template. A sample of the results is shown in Table 2.

Theory
To obtain the values of A1, B1, A2, B2, A3, B3, A4, B4, A5, B5, A6, B6, A7, B7, A8, B8, A9, B9, manually, we would first need to calculate the individual component for each, using equation 6 and equation 7, after which, the sum of each column divided by the total number of samples (n) will give the corresponding value for each of the terms. To do this we need the parameters ωi for each of the tidal constituents which can be found in Table 1, above, over a range of time t

                          ηsin (ωit)..……………………………………….………………….6
ηcos (ωit)………………….…………………………………………7
Ai =)..……………………………………………..8
Bi =)..……………………………………………..9
Where, ɳ = tidal water level at a particular time (m)
ωi = angular frequency at a particular time (rad/hr)
t= time (hrs)
i= 1-9
n= number of samples considered
Substituting the parameters into the equations above we obtain the same values for A1…B9 as shown in Table 3 below.

SOLUTION TO QUESTION 2

Carrying out harmonic analysis on the data in Table 3, we can find out the harmonic constants based on the formulae given below (the summation of the observed data over a range of samples and then divided by the sum by the total number of samples being considered).
i.e.
From Table 3 to calculate the constants a0 ,  and:
For i =range of 1 – 9

 a0 =
where n= number of samples being considered (18801)

Therefore,

a0 =)/(18815-15+1)
a0 =
a0 =5.30525

while

 = =
or simply,
= 2
by substituting into the above equation for when i=1, 2, 3……….9, we have
for i=1
= 2
where A1=-1.130 and B1=0.977 (from Table 3)
= 2  =2.9880

Similarly, all other constants for each constituent are calculated using the same method as above.

Likewise for the phase angle from i=1, 2, 3,……..9 going by the formulae

For i=1
=  ()
where A1=-1.130 and B1=0.977

=  
=-0.85788
In the same way, the rest of the phase angle for each constituent is calculated. See Table 3 for the full results of all the constituents.

Table 3: The Harmonic Constants (a0, ai  and φi)

Port: P234

Consti

M2

S2

K1

O1

N2

M4

M6

S4

MS4

Site: Liverpool (Gladstone Dock)

i =

1

2

3

4

5

6

7

8

9

Latitude:

53.4497

ωi =

0.5059

0.5236

0.2626

0.2433

0.4963

1.0118

1.5177

1.0472

1.03

Longitude:

-3.0181

ai =

2.9882

1.0218

0.1221

0.1372

0.6168763

0.2341694

0.0505411

0.0221983

0.104331
Start Date: 01JAN2005-00.00.00

φi =

-0.8578

-0.7361

-0.8295

1.5594

-0.729238

0.0460049

-0.633921

0.1507001

1.106581
End Date: 31DEC2005-23.45.00

n =

18801

a0 =

5.3053

Contributor: Proudman Oceanographic Laboratory

F=

0.064664

(Semidiurnal)
Datum information: The data refer to Admiralty (ACD) Chart Datum
Parameter code: ASLVBG02 = Sea level Bubbler (second tide sensor) gauge

The calculated constants in summary

M2

S2

K1

O1

N2

M4

M6

S4

MS4

ai =

2.9882

1.0218

0.1221

0.1372

0.6168763

0.2341694

0.0505411

0.0221983

0.104331

φi =

-0.8578

-0.7361

-0.8295

1.5594

-0.729238

0.0460049

-0.633921

0.1507001

1.106581

n =

18801

a0 =

5.3053

SOLUTION TO QUESTION 3

Comparing the predicted and the observed tidal water level, Figure 1 shows the graph of the observed and predicted tidal water level against time.
The observed and predicted values are taken from Table 3, and t = 0, 0.25, 0.5…371.25 hours.

Figure 1: Plot of comparison between observed water level and predicted water level vs. time

Figure 1 above gives the graph of both observed and predicted tidal levels against time (every quarter of an hour). It can be seen that the predicted values obtained were all positive; there were no negative values. Also from the graph it can be seen that the predicted value has a very close relationship with the measured value. If more constituents are used in the calculation, accuracy is improved. (In this case, all nine constituents were used)

SOLUTION TO QUESTION 4

Plotting of the individual constituents (M2, S2, K1, & O1), as shown in Figure 2 below, is done by calculating the equations below. The summary is given in Table 4:
η = a0 +a1 cos(ω1t – φ1) + a2 cos(ω2t – φ2) + a3 cos(ω3t – φ3) + a4 cos(ω4t – φ4)……..10
M2 = a1 cos(ω1t – φ1)……………………………………………………………………………………….11
Ss = a2 cos(ω2t – φ2)……………………………………………………………………12
K1 = a3 cos(ω3t – φ3)…………………………………………………………………..13
O1 = a4 cos(ω4t – φ4)…………………………………………………………………..14
The coefficient tidal values are taken from Table 4, and t = 0, 0.25, 0.5…146.25 hours.
Table 4: The Values for the Four Main Constituents

Time(hours)

M2

S2

K1

O1

0

1.954719

0.75721

0.082444

0.001565

0.25

1.654006

0.661187

0.076358

0.009903

0.5

1.326871

0.55385

0.069944

0.018203

0.75

0.978539

0.437037

0.063229

0.026437

1

0.614576

0.312747

0.05624

0.034572

1.25

0.240795

0.183105

0.04901

0.04258

1.5

-0.13683

0.05033

0.041569

0.05043

1.75

-0.51227

-0.08331

0.033948

0.058094

2

-0.87953

-0.21552

0.026181

0.065543

2.25

-1.23274

-0.34404

0.018301

0.072749

2.5

-1.56626

-0.46668

0.010343

0.079687

2.75

-1.87475

-0.58133

0.00234

0.086329

3

-2.1533

-0.68603

-0.00567

0.092653

3.25

-2.39745

-0.779

-0.01366

0.098633

3.5

-2.6033

-0.85864

-0.02159

0.104249

3.75

-2.76756

-0.92358

-0.02943

0.10948

4

-2.88761

-0.97273

-0.03714

0.114305

4.25

-2.96154

-1.00523

-0.04469

0.118708

4.5

-2.98815

-1.02052

-0.05205

0.122671

4.75

-2.96703

-1.01836

-0.05918

0.126181

5

-2.89852

-0.99878

-0.06606

0.129224

5.25

-2.7837

-0.9621

-0.07265

0.131789

5.5

-2.62441

-0.90896

-0.07893

0.133867

5.75

-2.4232

-0.84027

-0.08487

0.13545

6

-2.18328

-0.7572

-0.09045

0.136532

6.25

-1.90848

-0.66118

-0.09563

0.137109

6.5

-1.60319

-0.55384

-0.10041

0.137178

6.75

-1.2723

-0.43703

-0.10475

0.136741

7

-0.92108

-0.31274

-0.10864

0.135797

138.5

-0.7076

-0.55371

0.107173

0.102999

138.75

-1.06816

-0.43688

0.110779

0.097298

139

-1.41166

-0.31258

0.113908

0.091237

139.25

-1.7326

-0.18293

0.116545

0.084839

139.5

-2.02587

-0.05016

0.118681

0.078127

139.75

-2.28677

0.083478

0.120305

0.071125

140

-2.51114

0.215686

0.121411

0.063861

140.25

-2.6954

0.344202

0.121994

0.056361

140.5

-2.8366

0.46683

0.122052

0.048652

140.75

-2.93249

0.58147

0.121583

0.040763

141

-2.98153

0.68616

0.120591

0.032724

141.25

-2.98294

0.779111

0.119079

0.024563

141.5

-2.93671

0.85873

0.117054

0.016311

141.75

-2.84355

0.923656

0.114525

0.008

142

-2.70498

0.972778

0.111502

-0.00034

142.25

-2.52319

1.005256

0.107999

-0.00868

142.5

-2.3011

1.020533

0.104031

-0.01699

142.75

-2.04224

1.018349

0.099614

-0.02523

143

-1.75077

0.99874

0.094769

-0.03339

143.25

-1.43132

0.962042

0.089515

-0.04142

143.5

-1.08901

0.908884

0.083875

-0.04929

143.75

-0.72931

0.840174

0.077874

-0.05698

144

-0.35795

0.757089

0.071537

-0.06447

144.25

0.019126

0.661049

0.064893

-0.07171

144.5

0.395895

0.553699

0.057968

-0.07869

144.75

0.76634

0.436875

0.050794

-0.08537

145

1.124543

0.312575

0.043401

-0.09175

145.25

1.464782

0.182928

0.035821

-0.09778

145.5

1.781622

0.05015

0.028087

-0.10345

145.75

2.070001

-0.08349

0.020232

-0.10874

146

2.325312

-0.21569

0.012289

-0.11362

146.25

2.543477

-0.34421

0.004294

-0.11809

The plotting produced is a cos function.

Figure 2: Plot of the individual constituents (M2, S2, K1 & O1) vs. time

SOLUTION TO QUESTION 5

In order to decide what type of tide it is at this location we use the formulae given below:

F factor: F=  =

Where K1= a3 = 0.1221
O1= a4 = 0.1372
M2= a1 = 2.9882
S2= a2 = 1.0218

Substituting into the equation above we have:

F=

F= 0.064663.

From Table 2 above, the value of F falls in the category of “SEMIDIURNAL” (0.064< F< 0.25).
Therefore the type of tide at this location is the SEMIDIURNAL TIDE.

SOLUTION TO QUESTION 6

  1. The harmonic analysis method is a very good method for predicting tidal level
  2. It can enable us to predict future values of the tidal level since we know that tidal levels in estuaries don’t change often. i.e. they have a similar pattern over a period of time
  3. The results that are obtained by this method will be more accurate if more constituents are considered
  4. Taking a regular and consistent reading/ measurement we can also reduce the error in the analysis, because if the reading is not done in a regular way then the result would be wrong and inaccurate.
  5. This method is more straightforward and easier to use.
  6. It allows a larger data sample in the analysis

CONCLUSION: The Harmonic analysis is a better method for analysing tidal levels. Further, it gives the option of using more than the four main constituents, which gives it an edge over other methods. In short, it is better, faster and easier to use.