<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" dtd-version="3.0"><?xmltex \makeatother\@nolinetrue\makeatletter?>
  <front>
    <journal-meta>
<journal-id journal-id-type="publisher">OS</journal-id>
<journal-title-group>
<journal-title>Ocean Science</journal-title>
<abbrev-journal-title abbrev-type="publisher">OS</abbrev-journal-title>
<abbrev-journal-title abbrev-type="nlm-ta">Ocean Sci.</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">1812-0792</issn>
<publisher><publisher-name>Copernicus Publications</publisher-name>
<publisher-loc>Göttingen, Germany</publisher-loc>
</publisher>
</journal-meta>

    <article-meta>
      <article-id pub-id-type="doi">10.5194/os-13-599-2017</article-id><title-group><article-title>The double high tide at Port Ellen: Doodson's criterion revisited</article-title>
      </title-group><?xmltex \runningtitle{Double high tides}?><?xmltex \runningauthor{H. A. M. Byrne et al.}?>
      <contrib-group>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Byrne</surname><given-names>Hannah A. M.</given-names></name>
          
        <ext-link>https://orcid.org/0000-0001-6928-488X</ext-link></contrib>
        <contrib contrib-type="author" corresp="yes" rid="aff1">
          <name><surname>Green</surname><given-names>J. A. Mattias</given-names></name>
          <email>m.green@bangor.ac.uk</email>
        <ext-link>https://orcid.org/0000-0001-5090-1040</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Bowers</surname><given-names>David G.</given-names></name>
          
        </contrib>
        <aff id="aff1"><institution>Bangor University, School of Ocean Sciences, Menai Bridge,
Anglesey, LL59 5AB, UK</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">J. A. Mattias Green (m.green@bangor.ac.uk)</corresp></author-notes><pub-date><day>20</day><month>July</month><year>2017</year></pub-date>
      
      <volume>13</volume>
      <issue>4</issue>
      <fpage>599</fpage><lpage>607</lpage>
      <history>
        <date date-type="received"><day>15</day><month>March</month><year>2017</year></date>
           <date date-type="rev-request"><day>7</day><month>April</month><year>2017</year></date>
           <date date-type="rev-recd"><day>3</day><month>June</month><year>2017</year></date>
           <date date-type="accepted"><day>21</day><month>June</month><year>2017</year></date>
      </history>
      <permissions>
<license license-type="open-access">
<license-p>This work is licensed under the Creative Commons Attribution 3.0 Unported License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/3.0/">https://creativecommons.org/licenses/by/3.0/</ext-link></license-p>
</license>
</permissions><self-uri xlink:href="https://os.copernicus.org/articles/.html">This article is available from https://os.copernicus.org/articles/.html</self-uri>
<self-uri xlink:href="https://os.copernicus.org/articles/.pdf">The full text article is available as a PDF file from https://os.copernicus.org/articles/.pdf</self-uri>


      <abstract>
    <p>Doodson proposed a minimum criterion to predict the occurrence of double
high (or double low) waters when a higher-frequency tidal harmonic is added
to the semi-diurnal tide. If the phasing of the harmonic is optimal, the
condition for a double high water can be written <inline-formula><mml:math id="M1" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> where <inline-formula><mml:math id="M2" display="inline"><mml:mi>b</mml:mi></mml:math></inline-formula> is
the amplitude of the higher harmonic, <inline-formula><mml:math id="M3" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> is the amplitude of the
semi-diurnal tide, and <inline-formula><mml:math id="M4" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula> is the ratio of their frequencies. Here we expand
this criterion to allow for (i) a phase difference <inline-formula><mml:math id="M5" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> between the
semi-diurnal tide and the harmonic and (ii) the fact that the double high
water will disappear in the event that <inline-formula><mml:math id="M6" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> becomes large enough for the
higher harmonic to be the dominant component of the tide. This can happen,
for example, at places or times where the semi-diurnal tide is very small.
The revised parameter is <inline-formula><mml:math id="M7" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula>, where <inline-formula><mml:math id="M8" display="inline"><mml:mi>r</mml:mi></mml:math></inline-formula> is a number generally less
than <inline-formula><mml:math id="M9" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, although equal to <inline-formula><mml:math id="M10" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula> when <inline-formula><mml:math id="M11" display="inline"><mml:mrow><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>. The theory predicts that a
double high tide will form when this parameter exceeds 1 and then disappear
when it exceeds a value of order <inline-formula><mml:math id="M12" display="inline"><mml:mrow><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> and the higher harmonic becomes
dominant. We test these predictions against observations at Port Ellen in
the Inner Hebrides of Scotland. For most of the data set, the largest
harmonic of the semi-diurnal tide is the sixth diurnal component, for which
<inline-formula><mml:math id="M13" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula>. The principal lunar and solar semi-diurnal tides are about equal at
Port Ellen and so the semi-diurnal tide becomes very small twice a month at
neap tides (here defined as the smallest fortnightly tidal range). A double
high water forms when <inline-formula><mml:math id="M14" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> first exceeds a minimum value of about 1.5
as neap tides are approached and then disappears as <inline-formula><mml:math id="M15" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> then exceeds a
second limiting value of about 10 at neap tides in agreement with the
revised criterion.</p>
  </abstract>
    </article-meta>
  </front>
<body>
      

      <?xmltex \hack{\newpage}?>
<sec id="Ch1.S1" sec-type="intro">
  <title>Introduction</title>
      <p>Double high waters in the semi-diurnal tide are an intriguing and sometimes
puzzling feature of coastal tidal observations. In a double high water, the
tide rises to a first maximum followed by a short dip in water level before
it rises again to a second maximum and then falls towards the subsequent low
water. This is because the spatial scale of the higher harmonics of the tide
is shorter than that of the dominating semi-diurnal constituent.
Consequently, any location with a double high has a nearby location with a
double low tide. The most famous and best studied example of a double high
water is that at Southampton on the south coast of England, where the
extended period of deep water associated with the double high tide gave the
port a commercial advantage over its rivals. The Southampton double high has
an associated double low in Weymouth, some 50 km away <xref ref-type="bibr" rid="bib1.bibx12" id="paren.1"/>
Other examples of double high waters, and associated double lows, can be found
along the eastern side of the North Sea, e.g. a double in Den Helder in the
Netherlands. On the US east coast, Woods Hole (among other locations)
experiences a double high, whereas a double low water is seen at Providence
in Narragansett Bay (<xref ref-type="bibr" rid="bib1.bibx10 bib1.bibx12" id="altparen.2"/>; see also <xref ref-type="bibr" rid="bib1.bibx5" id="altparen.3"/>, for a
discussion).</p>
      <p>No explanation of double high (or low) waters can be found in the direct
action on the ocean of the tidal forces of the moon and sun. To create a
double high water in the semi-diurnal tide it is necessary to add a higher-frequency oscillation with  an appropriate phase and sufficient amplitude.
Higher harmonics of the semi-diurnal tide (that is oscillations with a
frequency higher than twice per day) are created as the tidal wave enters
shallow water, either through the reduction of water depth, through the
enhanced effects of quadratic bottom friction, or through streamline
curvature around coasts and sandbanks <xref ref-type="bibr" rid="bib1.bibx9" id="paren.4"/>. Transient
seiches created by the tide can also serve to provide the right conditions
for double tides <xref ref-type="bibr" rid="bib1.bibx2" id="paren.5"/>.</p>
      <p>The amplitude and phase of the higher-frequency oscillation must meet certain
conditions to produce a double high water. The simplest of these conditions
was described by <xref ref-type="bibr" rid="bib1.bibx3" id="text.6"/>. If the frequency of the oscillation
is <inline-formula><mml:math id="M16" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula> times that of the semi-diurnal tide and the oscillation has a trough
exactly coinciding with the peak of high water, a double high water will form
if
          <disp-formula id="Ch1.E1" content-type="numbered"><mml:math id="M17" display="block"><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi>b</mml:mi><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow><mml:mi>a</mml:mi></mml:mfrac></mml:mstyle><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M18" display="inline"><mml:mi>b</mml:mi></mml:math></inline-formula> is the amplitude of the oscillation and <inline-formula><mml:math id="M19" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> the amplitude of the
semi-diurnal tide. For example, if the oscillation has the period of a
quarter of a day, then <inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:math></inline-formula> and the amplitude must be at least one-fourth that of
the semi-diurnal tide to produce a double high water. This is a minimum
requirement – in most cases, the quarter diurnal component will not be big
enough on its own to produce a double high water and additional higher
harmonics will be required. As the frequency, <inline-formula><mml:math id="M21" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, of the higher harmonic
increases, Eq. (1) shows that the required amplitude
ratio, <inline-formula><mml:math id="M22" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula>, becomes smaller.</p>
      <p>It is worth noting here that the amplitude, <inline-formula><mml:math id="M23" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula>, of the semi-diurnal tide in
Eq. (1) is that observed on the day – i.e. it is the sum of the
lunar and solar semi-diurnal tidal constituents on that day. The value of it
therefore changes through the fortnightly spring-neap cycle.<fn id="Ch1.Footn1"><p>In the
following we use “neap” to describe the smallest fortnightly tidal range.
Note that this is not necessarily at first or third lunar quadrature.</p></fn> The
semi-diurnal tide on a particular day, changing in amplitude from day to day
in this way is called the D2 tide <xref ref-type="bibr" rid="bib1.bibx8" id="paren.7"/>. Similarly, the amplitude of
the higher-frequency oscillation, <inline-formula><mml:math id="M24" display="inline"><mml:mi>b</mml:mi></mml:math></inline-formula>, is also that observed on the day and
will be different on other days. If the higher frequency oscillation is a
harmonic of the semi-diurnal tide it is referred to as a D4 tide in the case
of a quarter diurnal oscillation, D6 if the period is about 4 h, and so
on. One way to increase the value of the ratio <inline-formula><mml:math id="M25" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> in the Doodson criterion
is to reduce the amplitude, <inline-formula><mml:math id="M26" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula>, of the D2 tide. A good place to look for
double high waters, therefore, is near an amphidrome in the semi-diurnal
tide. In fact, Southampton lies close to the nodal line for the semi-diurnal
tide in the English Channel. Taking another example, in the Irish Sea, there
is a degenerate amphidrome for the semi-diurnal tide close to Courtown on the
Irish coast. Here, the semi-diurnal tide is small and the higher harmonics
can be the dominant components in the tide <xref ref-type="bibr" rid="bib1.bibx8" id="normal.8"/>. Alternatively, in
places where the two principal semi-diurnal constituents – those of the sun
and moon – are about equal, the D2 tide will become small twice a month at
neap tides and any higher-frequency oscillations will assert themselves. In
South Australia, this phenomenon is called the dodge tide
<xref ref-type="bibr" rid="bib1.bibx6" id="paren.9"/>. Similarly, double low waters will occur in locations
near areas with double high waters. This is because we expect double tides
near a semi-diurnal amphidrome (where the semi-diurnal amplitude is small).
Consequently, a short movement across the amphidrome, the phase of the
semi-diurnal tide will change by 180<inline-formula><mml:math id="M27" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>, and if the phase of the overtide
does not change, there will be a switch from a double high to a double low
tide.</p>
      <p>These considerations about small values of <inline-formula><mml:math id="M28" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> lead us to the conclusion that
there must be a second criterion for double high waters. As <inline-formula><mml:math id="M29" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> increases
from the critical value given by Eq. (<xref ref-type="disp-formula" rid="Ch1.E1"/>), a double high water
will first form and then, as <inline-formula><mml:math id="M30" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> continues to increase, the double high
water will disappear as the higher harmonic becomes dominant. The Doodson
criterion in Eq. (<xref ref-type="disp-formula" rid="Ch1.E1"/>) for the production of double high tides
is therefore a necessary requirement but it is not, on its own, a sufficient
requirement, and it explicitly only deals with harmonics smaller than M2. To
develop a more general criterion we need also to impose an upper limit on the
ratio <inline-formula><mml:math id="M31" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> to allow for times and places where the fundamental harmonic (in
Doodson's case, M2) is small compared to its harmonics. We can imagine that a
more general condition for double high waters might take the form
          <disp-formula id="Ch1.E2" content-type="numbered"><mml:math id="M32" display="block"><mml:mrow><mml:mi>x</mml:mi><mml:mo>&gt;</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi>b</mml:mi><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow><mml:mi>a</mml:mi></mml:mfrac></mml:mstyle><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M33" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> is a function to be determined. Our first aim in this paper is thus
to explore the nature of <inline-formula><mml:math id="M34" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula>.</p>
      <p>A further restriction on the Doodson criterion is that it applies only to the
case where the phase of the higher harmonic is optimal – that is, in the case
of double high waters the minimum in the harmonic coincides with the maximum
in the semi-diurnal tide. <xref ref-type="bibr" rid="bib1.bibx2" id="text.10"/> suggested a modification to the
Doodson criterion to allow for a phase difference between the main tide and
its harmonic. A second aim of this paper is to take the opportunity to test
that idea against observations.</p>
</sec>
<sec id="Ch1.S2">
  <title>Theoretical background</title>
      <p>Figure 1 illustrates the formation and subsequent disappearance of a double
high water for the sum of a semi-diurnal D2 tide (amplitude <inline-formula><mml:math id="M35" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula>) and its D6
harmonic (amplitude <inline-formula><mml:math id="M36" display="inline"><mml:mi>b</mml:mi></mml:math></inline-formula>). In Fig. 1a the ratio <inline-formula><mml:math id="M37" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.25</mml:mn></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M38" display="inline"><mml:mrow><mml:mi>a</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> m, and the
phase difference is 0.5 h (i.e. the minimum in D6 occurs half an hour
after the maximum in D2). Adding a D6 curve to D2 creates a double high water
and also a double low water, as the maximum in D6 also coincides with the low
water in D2. In Fig. 1b the amplitude of D2 has been reduced to 0.1 m, keeping
the D6 curve the same. There is now no evidence of a double high water – the
tide is instead best described as sixth diurnal with a weak remnant
semi-diurnal modulation in mean water level.</p>
      <p>In general, the sum of a semi-diurnal tide and a single higher harmonic can be
written as
          <disp-formula id="Ch1.E3" content-type="numbered"><mml:math id="M39" display="block"><mml:mrow><mml:mi>y</mml:mi><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mi>cos⁡</mml:mi><mml:mo>[</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mi>t</mml:mi><mml:mo>]</mml:mo><mml:mo>-</mml:mo><mml:mi>b</mml:mi><mml:mi>cos⁡</mml:mi><mml:mo>[</mml:mo><mml:mi>n</mml:mi><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>)</mml:mo><mml:mo>]</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M40" display="inline"><mml:mi>y</mml:mi></mml:math></inline-formula> is sea level and <inline-formula><mml:math id="M41" display="inline"><mml:mrow><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">2</mml:mn><mml:mi mathvariant="italic">π</mml:mi><mml:mo>/</mml:mo><mml:mn mathvariant="normal">12</mml:mn></mml:mrow></mml:math></inline-formula> lunar hours is the angular
frequency of the semi-diurnal tide. The phase <inline-formula><mml:math id="M42" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> represents the time
difference between the maximum in the semi-diurnal tide and the nearest
minimum in the higher harmonic, so that D2 has a maximum for <inline-formula><mml:math id="M43" display="inline"><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> and the
<inline-formula><mml:math id="M44" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula> harmonic a minimum at <inline-formula><mml:math id="M45" display="inline"><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi></mml:mrow></mml:math></inline-formula>. Time, <inline-formula><mml:math id="M46" display="inline"><mml:mi>t</mml:mi></mml:math></inline-formula>, is measured from zero at the
high water in the semi-diurnal tide. For a small range of times around <inline-formula><mml:math id="M47" display="inline"><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>
and <inline-formula><mml:math id="M48" display="inline"><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi></mml:mrow></mml:math></inline-formula> the cosine curves can be represented by quadratic curves (using
the expansion that <inline-formula><mml:math id="M49" display="inline"><mml:mrow><mml:mi>cos⁡</mml:mi><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn><mml:msup><mml:mi>x</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> for small <inline-formula><mml:math id="M50" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula>):
          <disp-formula id="Ch1.E4" content-type="numbered"><mml:math id="M51" display="block"><mml:mrow><mml:mi>y</mml:mi><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mo>[</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn><mml:msup><mml:mi mathvariant="italic">ω</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:msup><mml:mi>t</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>]</mml:mo><mml:mo>-</mml:mo><mml:mi>b</mml:mi><mml:mo>[</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:msup><mml:mi mathvariant="italic">ω</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>]</mml:mo><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F1"><caption><p><bold>(a)</bold> A double high (and double low) water formed by the addition of a
semi-diurnal curve (dash-dot line, amplitude 1 m) and a sixth-diurnal harmonic
(dashed line, amplitude 0.25 m). There is phase difference of 0.5 h
between the curves, and time is measured from the peak of the semi-diurnal
tide. “D” marks the centre of the dip between the high waters. At D, the
gradient of <inline-formula><mml:math id="M52" display="inline"><mml:mi>y</mml:mi></mml:math></inline-formula> with respect to time is zero. <bold>(b)</bold> As <bold>(a)</bold> but the amplitude of the
semi-diurnal curve is reduced to 0.1 m. The double high and low waters
have disappeared and the tide is effectively sixth diurnal.</p></caption>
        <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://os.copernicus.org/articles/13/599/2017/os-13-599-2017-f01.png"/>

      </fig>

      <p>When there is a double high water, there is a turning point (marked as D in
Fig. 1a) which marks the centre of the dip in sea level between the two high
waters either side. At this turning point, occurring at time <inline-formula><mml:math id="M53" display="inline"><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M54" display="inline"><mml:mrow><mml:mtext>d</mml:mtext><mml:mi>y</mml:mi><mml:mo>/</mml:mo><mml:mtext>d</mml:mtext><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>
and so, for small values of <inline-formula><mml:math id="M55" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula>, from Eq. (<xref ref-type="disp-formula" rid="Ch1.E4"/>),
          <disp-formula id="Ch1.E5" content-type="numbered"><mml:math id="M56" display="block"><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mfenced open="[" close="]"><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfrac></mml:mstyle></mml:mfenced><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>
      <p>Since the term in brackets is always greater than 1, <inline-formula><mml:math id="M57" display="inline"><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo>&gt;</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi></mml:mrow></mml:math></inline-formula>. It is a
feature of double high waters that the turning point lies further away from
the (semidiurnal) high water than does the minimum in the harmonic that
creates it.</p>
      <p>For the dip to be a minimum between the double high waters, a further
condition is that <inline-formula><mml:math id="M58" display="inline"><mml:mrow><mml:msup><mml:mtext>d</mml:mtext><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mi>y</mml:mi><mml:mo>/</mml:mo><mml:mtext>d</mml:mtext><mml:msup><mml:mi>t</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> at the dip. At this point it is necessary
to include the next term (a function of <inline-formula><mml:math id="M59" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">4</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>) in the expansion of the
cosine curve for the higher harmonic. Applying the condition that
<inline-formula><mml:math id="M60" display="inline"><mml:mrow><mml:msup><mml:mtext>d</mml:mtext><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mi>y</mml:mi><mml:mo>/</mml:mo><mml:mtext>d</mml:mtext><mml:msup><mml:mi>t</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:msub><mml:mo>|</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:msub><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> gives the condition for a double high water as
          <disp-formula id="Ch1.E6" content-type="numbered"><mml:math id="M61" display="block"><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi>b</mml:mi><mml:mi>a</mml:mi></mml:mfrac></mml:mstyle><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where
          <disp-formula id="Ch1.E7" content-type="numbered"><mml:math id="M62" display="block"><mml:mrow><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>=</mml:mo><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>[</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mo>(</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>)</mml:mo><mml:msup><mml:mi mathvariant="italic">ω</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>(</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>]</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        which is the result obtained by <xref ref-type="bibr" rid="bib1.bibx2" id="text.11"/>. Note that when the
phase difference <inline-formula><mml:math id="M63" display="inline"><mml:mrow><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> it follows that <inline-formula><mml:math id="M64" display="inline"><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> and therefore that <inline-formula><mml:math id="M65" display="inline"><mml:mrow><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>=</mml:mo><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>; the condition for a double high water then becomes the Doodson
criterion in Eq. (<xref ref-type="disp-formula" rid="Ch1.E1"/>). For all other values of <inline-formula><mml:math id="M66" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula>, <inline-formula><mml:math id="M67" display="inline"><mml:mrow><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>&lt;</mml:mo><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>; then, according to Eq. (<xref ref-type="disp-formula" rid="Ch1.E6"/>), <inline-formula><mml:math id="M68" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> must be larger than the
value required for zero <inline-formula><mml:math id="M69" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> in order to create a double high water. It is
possible for <inline-formula><mml:math id="M70" display="inline"><mml:mrow><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>, and hence <inline-formula><mml:math id="M71" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>, to be negative. This will happen
when either <inline-formula><mml:math id="M72" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> is small or <inline-formula><mml:math id="M73" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> is large. When <inline-formula><mml:math id="M74" display="inline"><mml:mrow><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> is negative, it is
impossible to satisfy Eq. (<xref ref-type="disp-formula" rid="Ch1.E6"/>) whatever the value of <inline-formula><mml:math id="M75" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula>, and no
double high water can be formed (the actual physical interpretation of a
negative <inline-formula><mml:math id="M76" display="inline"><mml:mrow><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> is intriguing, but left for a future paper).</p>
      <p>As an illustration, we show in Fig. 2 the solution to the approximate
parameterization in Eqs. (<xref ref-type="disp-formula" rid="Ch1.E4"/>)–(<xref ref-type="disp-formula" rid="Ch1.E6"/>) for the case of
<inline-formula><mml:math id="M77" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula>, that is the sum of D2 and D6 tides. The solution is symmetrical for
negative and positive values of <inline-formula><mml:math id="M78" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> and so in Fig. 2 we have drawn the
curves as a function of the absolute value of <inline-formula><mml:math id="M79" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula>. It is necessary to
iterate to reach the solution – the critical value of <inline-formula><mml:math id="M80" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> must satisfy both
Eqs. (5) and (6). For D2 and D6 tides, the critical value of <inline-formula><mml:math id="M81" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> for zero
phase is <inline-formula><mml:math id="M82" display="inline"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>/</mml:mo><mml:mn mathvariant="normal">9</mml:mn></mml:mrow></mml:math></inline-formula>, or 0.111. As the phase difference increases from zero, the
critical value of <inline-formula><mml:math id="M83" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> also increases, so that for a phase difference of 1 h, it is necessary for the amplitude of D6 to be at least 0.27 times that
of D2 to produce a double high tide. Also shown in Fig. 2 is the exact
numerical solution to the problem, starting with the cosine curves in Eq. (3).
The approximation of the cosine curves as quadratic curves close to
their maxima and minima means that the analytical solution is not exact, but
it does capture the essential features of the numerical solution, especially
for <inline-formula><mml:math id="M84" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> less than about 1.5 h. In fact, at values of <inline-formula><mml:math id="M85" display="inline"><mml:mrow><mml:mo>|</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>|</mml:mo></mml:mrow></mml:math></inline-formula> greater
than about 1 h the dip becomes so detached from the high water that the
tidal curve is no longer recognizable as a double high water. In practice,
therefore, for a recognizable double high water produced by D2 and D6 tides,
we can limit attention to the region in which <inline-formula><mml:math id="M86" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> is less than about 1 h in Fig.  2.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F2"><caption><p>Theoretical condition for the formation of double high waters
illustrated for the addition of a sixth diurnal (D6) harmonic to a
semi-diurnal D2 tide. The <inline-formula><mml:math id="M87" display="inline"><mml:mi>y</mml:mi></mml:math></inline-formula> axis is the ratio of the amplitudes of the
harmonic to the principal tide. The <inline-formula><mml:math id="M88" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> axis is the time difference between the
minimum in D6 and the maximum in D2 (plus D1 in the case of the observations).
The solid curve shows the critical value of <inline-formula><mml:math id="M89" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> for the formation of double
high waters according to the analytical solution (Eq. <xref ref-type="disp-formula" rid="Ch1.E6"/> in text).
The dotted curve is an exact numerical solution to the problem. The points
show values of <inline-formula><mml:math id="M90" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M91" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> at Port Ellen. Points have been plotted as
solid circles for tides where there is a double high water, grey filled
circles where there is a stand in the tide near high water and open circles
when neither a stand nor a double high water is observed.</p></caption>
        <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://os.copernicus.org/articles/13/599/2017/os-13-599-2017-f02.pdf"/>

      </fig>

      <p>As the amplitude ratio <inline-formula><mml:math id="M92" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> continues to increase above the critical value
for the formation of a double high water, there is a gradual transition
towards a tide dominated by the higher harmonic. As this happens, the level
of water in the dip between the high waters falls towards the level of the
low tide (Fig. 1). There comes a point when the level of water between the
dips is virtually the same as that at the low tide, when we can say that the
transition to the higher harmonic is complete and the double tide has
disappeared. At the time of the minimum in the dip, <inline-formula><mml:math id="M93" display="inline"><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> and sea level is
given by
          <disp-formula id="Ch1.E8" content-type="numbered"><mml:math id="M94" display="block"><mml:mrow><mml:mi>y</mml:mi><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn><mml:msup><mml:mi mathvariant="italic">ω</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:msup><mml:mi>t</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>)</mml:mo><mml:mo>-</mml:mo><mml:mi>b</mml:mi><mml:mo>[</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mo>(</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>)</mml:mo><mml:msup><mml:mi mathvariant="italic">ω</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>]</mml:mo><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
        The first term in the right-hand side of this equation represents the fall in
water level between <inline-formula><mml:math id="M95" display="inline"><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M96" display="inline"><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> due to the semi-diurnal tide; the second
term is the fall in the same time due to the harmonic. The first term will
occur in the absence of the harmonic, so the size of the dip due to the
harmonic is equal to the second term. Referring to Eq. (<xref ref-type="disp-formula" rid="Ch1.E7"/>), this
second term can be written as <inline-formula><mml:math id="M97" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>/</mml:mo><mml:mi>n</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>. The size of the dip relative to the
amplitude of the semi-diurnal tide is therefore <inline-formula><mml:math id="M98" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>/</mml:mo><mml:mi>n</mml:mi><mml:msup><mml:mo>)</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>. The higher
harmonic will become dominant when this ratio exceeds a certain value, of
order 1. That is the upper limit of <inline-formula><mml:math id="M99" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> for a recognizable double high
water is
          <disp-formula id="Ch1.E9" content-type="numbered"><mml:math id="M100" display="block"><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow><mml:mi>f</mml:mi></mml:mfrac></mml:mstyle><mml:mo>&gt;</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi>b</mml:mi><mml:mi>a</mml:mi></mml:mfrac></mml:mstyle><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where the factor <inline-formula><mml:math id="M101" display="inline"><mml:mi>f</mml:mi></mml:math></inline-formula> represents the size of the dip relative to <inline-formula><mml:math id="M102" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula>. For
example, if <inline-formula><mml:math id="M103" display="inline"><mml:mrow><mml:mi>f</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M104" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula>, the upper bound for <inline-formula><mml:math id="M105" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> is <inline-formula><mml:math id="M106" display="inline"><mml:mrow><mml:mn mathvariant="normal">9</mml:mn><mml:mo>/</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>.  In
general, we can write the critical condition for the formation of a
recognizable double high water in the form
          <disp-formula id="Ch1.E10" content-type="numbered"><mml:math id="M107" display="block"><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msup><mml:mi>n</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow><mml:mi>f</mml:mi></mml:mfrac></mml:mstyle><mml:mo>&gt;</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi>b</mml:mi><mml:mi>a</mml:mi></mml:mfrac></mml:mstyle><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        which is the condition we want to test. Note that the theory leads us to
expect that the same parameter, namely <inline-formula><mml:math id="M108" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula>, is important in predicting
the initial onset of the double high water and the disappearance of the high
water as the harmonic becomes dominant.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F3"><caption><p>Location of Port Ellen in the Inner Hebrides, Scotland. The open
circle on the channel marks the M2 amphidrome. The coastline was drawn using
GEBCO2014 <xref ref-type="bibr" rid="bib1.bibx11" id="paren.12"><named-content content-type="pre">see</named-content><named-content content-type="post">and <uri>www.gebco.net</uri></named-content></xref>, available from
<uri>https://www.bodc.ac.uk/data/hosted_data_systems/gebco_gridded_bathymetry_data/</uri>.</p></caption>
        <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://os.copernicus.org/articles/13/599/2017/os-13-599-2017-f03.pdf"/>

      </fig>

</sec>
<sec id="Ch1.S3">
  <title>Observations</title>
      <p>Port Ellen lies on the south coast of Islay, an island which is part of the
Inner Hebrides on the west coast of Scotland (Fig. 3). Tidal data sampled
every 15 min at Port Ellen are available from the UK Tide Gauge Network as
part of The National Tidal and Sea Level Facility (available from
<uri>https://www.bodc.ac.uk/data/hosted_data_systems/sea_level/uk_tide_gauge_network/processed/</uri>),
and Fig. 4 shows a record from Port Ellen for the second half of February,
2010. The 2-week observational period we use began on a new and ended on a
full moon. The semi-diurnal tide is unusual because the two principal
semi-diurnal tidal constituents (M2 and S2) are about equal in amplitude and,
as a consequence, the semi-diurnal tide virtually disappears twice each month
when these two constituents are 180<inline-formula><mml:math id="M109" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula> out of phase. This effect can be
seen for 25 February in Fig. 4. A similar phenomenon is the dodge tide in
South Australia <xref ref-type="bibr" rid="bib1.bibx6" id="paren.13"/>, and the vanishing tide in Courtown,
Ireland <xref ref-type="bibr" rid="bib1.bibx7" id="paren.14"/>.</p>
      <p>The diurnal tide is also important at Port Ellen and produces a significant
diurnal inequality which can be seen at the beginning and end of the record
in Fig. 4. In the figure, the diurnal inequality takes the form of the low
tide in the middle of the day being higher than the low tide at the beginning
and end of the day. The diurnal inequality is important in the selection of
the morning or afternoon tide for the production of double high waters, as we
shall see.</p>
      <p>The double high water is intermittent and occurs most clearly in the morning
high tide in the first part of the record (18 and 19 February) and the
afternoon high tide towards the end of the record (27 and 28 February). There
are also times when there is a stand in water level around the time of high
water – it is sometimes difficult to tell if there is actually a double high
tide present on these occasions or not. These “stand” tides are observed on
all tidal cycles between the morning of 15 February and the morning of 17 February,
and on the mornings of 20, 23, and 27 February.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F4"><caption><p>Water level record for Port Ellen in the second half of February,
2010. Double high waters are seen on the mornings of 18 and 19 February and
the afternoons of 27 and 28 February.</p></caption>
        <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://os.copernicus.org/articles/13/599/2017/os-13-599-2017-f04.png"/>

      </fig>

<sec id="Ch1.S3.SS1">
  <title>Harmonic analysis</title>
      <p>We have analysed short portions of the water level record for the amplitude
and phases of key harmonics using harmonic analysis <xref ref-type="bibr" rid="bib1.bibx4" id="paren.15"><named-content content-type="pre">see e.g.</named-content><named-content content-type="post">for
details</named-content></xref>. The method was first described by <xref ref-type="bibr" rid="bib1.bibx1" id="text.16"/> to
establish the fact that higher harmonics of the semi-diurnal tide were
important at Southampton. Later, <xref ref-type="bibr" rid="bib1.bibx3" id="text.17"/> used the technique to
test their inequality in Eq. (<xref ref-type="disp-formula" rid="Ch1.E1"/>) against observations of the
tide at Southampton.</p>
      <p>Because the diurnal tide is important at Port Ellen, we have analysed the
data for harmonics starting with the diurnal tide as the first harmonic. The
data set of Fig. 4 was broken up into 25 h segments, each starting at
midnight on the chosen day. The selected segments of data therefore run for
1 h into the next day. The harmonic analysis is then applied to each
segment of data to calculate the amplitude and phase of a diurnal constituent
(D1, period 25 h), a semi-diurnal constituent (D2, period 12.5 h), a
quarter-diurnal constituent  (D4, period 6.25 h), and a sixth-diurnal
constituent (D6, period 4.17 h). Note that this analysis produces
amplitudes and phases of the harmonics applicable to that day only. The
harmonics are therefore not the same as the tidal constituents such as M2 and
M4 treated by a full harmonic tidal analysis and which have constant
amplitude and phase on all days. The amplitude and phase of these daily
harmonics changes from day to day. The usefulness of the “D” harmonics has
been discussed by <xref ref-type="bibr" rid="bib1.bibx8" id="text.18"/>. It enables the relative importance of the
principal tide and its harmonics to be established on each day. The
relationship between the relative amplitudes and phases of the harmonics and
the production (or non-production) of double high waters can then be
explored.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F5"><caption><p><bold>(a)</bold> Observed and fitted tidal curves for 19 February. The double
high water on the morning is reproduced by the fitted curve. <bold>(b)</bold> D1 (blue), D2
(red), D4 (green), and D6 (purple) harmonics for 19 February. The high tide
in D1 (black line) occurs in the middle of the day and causes the diurnal
inequality seen in Fig. 4. The D6 harmonic (lightest grey) has a larger
amplitude than D4 and the minimum in D6 occurs closer to the high D2 than
does the minimum in D4. It is the D6 harmonic that is primarily responsible
for the double high water at Port Ellen.</p></caption>
          <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://os.copernicus.org/articles/13/599/2017/os-13-599-2017-f05.png"/>

        </fig>

      <p>Figure <xref ref-type="fig" rid="Ch1.F5"/> shows an example of the curve-fitting for 19 February,
whereas Fig. <xref ref-type="fig" rid="Ch1.F6"/> shows the daily amplitudes of the
harmonics throughout the observational period. Figure <xref ref-type="fig" rid="Ch1.F5"/>a shows
the observed and fitted curve for this day. The double high water on the
morning tide is clearly shown and the fitted curve reproduces this. There is
no double high water in the afternoon, in either the observed or fitted tide.
The fitted curve is not perfect, however, due to for instance meteorological
effects, seiches, and instrument noise. Most importantly, there is a double
low water in the fitted curve, which is not seen in the observations. Figure
5b shows the nature of the harmonics which sum to give the fitted curve on
this day. The semi-diurnal harmonic, D2, has the largest amplitude, followed
by D1, then D6, and finally D4. The low water in the D6 curve is close to the
high water in the D2 curve, a requirement for the formation of double high
waters. In contrast, the amplitude of D4 is smaller than D6 and its phase is
not optimal – the low in D4 occurs about 2 h after the semi-diurnal high
water in Fig. <xref ref-type="fig" rid="Ch1.F5"/>b (see also Fig. <xref ref-type="fig" rid="Ch1.F6"/>).</p>
      <p>The fact that there is a double high water in the morning of 19 February, but
not in the afternoon, can be explained by the effect of the D1 tide. The
diurnal tide is rising in the morning and falling in the afternoon (Figs. 5b
and <xref ref-type="fig" rid="Ch1.F6"/>). This has the effect of pushing the time of high
water forwards in time in the morning, towards noon, and dragging it
backwards in time, also towards noon, in the afternoon. The time interval
between the two high waters produced by D2 and D1 together is therefore less
than 12.5 h. It is thus possible for the low tide in D6 to be close to
the high tide in D1 and D2 in the morning, but not so close in the afternoon.
In fact, the time interval between the D6 low and the (D1 <inline-formula><mml:math id="M110" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula> D2) high on
19 February is zero in the morning and 0.8 h in the afternoon. The close
coincidence of the low water in D6 and the high water in (D1 <inline-formula><mml:math id="M111" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula> D2) produces a
double high water in the morning but not in the afternoon. For this reason,
it is important to consider the time of the high water in D1 <inline-formula><mml:math id="M112" display="inline"><mml:mo>+</mml:mo></mml:math></inline-formula> D2.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F6"><caption><p>The daily amplitudes for fits of D1 (blue), D2 (red), D4
(green), and D6 (purple).</p></caption>
          <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://os.copernicus.org/articles/13/599/2017/os-13-599-2017-f06.pdf"/>

        </fig>

      <p>The harmonic analysis illustrated in Fig. <xref ref-type="fig" rid="Ch1.F5"/> has been applied
to each of the days in the record and the results which we shall use in the
analysis are shown in Table 1. The phase of D6 is expressed as the absolute
value of the time interval between low water in this harmonic and the nearest
high water in D1 and D2 combined.  The semi-diurnal tide virtually disappears
on 24 February, at neap tides, because of the equality in the amplitudes of the
solar and lunar tides.</p>

<?xmltex \floatpos{t}?><table-wrap id="Ch1.T1"><caption><p>The first column shows the day in February 2010 divided into two
tides, one in the morning (a.m.) and a second in the afternoon (p.m.). <inline-formula><mml:math id="M113" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> is
the amplitude of the semi-diurnal tide in metres and is the same in the
morning and in the afternoon. <inline-formula><mml:math id="M114" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> is the ratio of the amplitudes of the
sixth diurnal tide D6 to that of D2. <inline-formula><mml:math id="M115" display="inline"><mml:mrow><mml:mo>|</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>|</mml:mo></mml:mrow></mml:math></inline-formula> is the absolute value of the
time difference in hours between the maximum in D1 and D2 combined and the
minimum in D6; <inline-formula><mml:math id="M116" display="inline"><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> is the time interval calculated with Eq. (<xref ref-type="disp-formula" rid="Ch1.E5"/>),
and <inline-formula><mml:math id="M117" display="inline"><mml:mrow><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> is the parameter in Eq. (<xref ref-type="disp-formula" rid="Ch1.E7"/>). Finally, <inline-formula><mml:math id="M118" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> is the
discriminator for double high waters.</p></caption><oasis:table frame="topbot"><?xmltex \begin{scaleboxenv}{.95}[.95]?><oasis:tgroup cols="7">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="right"/>
     <oasis:colspec colnum="3" colname="col3" align="right"/>
     <oasis:colspec colnum="4" colname="col4" align="right"/>
     <oasis:colspec colnum="5" colname="col5" align="right"/>
     <oasis:colspec colnum="6" colname="col6" align="right"/>
     <oasis:colspec colnum="7" colname="col7" align="right"/>
     <oasis:thead>
       <oasis:row rowsep="1">  
         <oasis:entry colname="col1">Tide</oasis:entry>  
         <oasis:entry colname="col2"><inline-formula><mml:math id="M121" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col3"><inline-formula><mml:math id="M122" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col4"><inline-formula><mml:math id="M123" display="inline"><mml:mrow><mml:mo>|</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>|</mml:mo></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col5"><inline-formula><mml:math id="M124" display="inline"><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M125" display="inline"><mml:mrow><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M126" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>  
         <oasis:entry colname="col1">14 Feb, a.m.</oasis:entry>  
         <oasis:entry colname="col2">0.28</oasis:entry>  
         <oasis:entry colname="col3">0.169</oasis:entry>  
         <oasis:entry colname="col4">0.6</oasis:entry>  
         <oasis:entry colname="col5">1.75</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M127" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">4.57</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M128" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.77</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">14 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.28</oasis:entry>  
         <oasis:entry colname="col3">0.169</oasis:entry>  
         <oasis:entry colname="col4">0.5</oasis:entry>  
         <oasis:entry colname="col5">1.46</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M129" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.43</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M130" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.07</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">15 Feb, a.m.<inline-formula><mml:math id="M131" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.29</oasis:entry>  
         <oasis:entry colname="col3">0.207</oasis:entry>  
         <oasis:entry colname="col4">0.7</oasis:entry>  
         <oasis:entry colname="col5">1.51</oasis:entry>  
         <oasis:entry colname="col6">2.27</oasis:entry>  
         <oasis:entry colname="col7">0.47</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">15 Feb, p.m.<inline-formula><mml:math id="M132" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.29</oasis:entry>  
         <oasis:entry colname="col3">0.207</oasis:entry>  
         <oasis:entry colname="col4">0.4</oasis:entry>  
         <oasis:entry colname="col5">0.86</oasis:entry>  
         <oasis:entry colname="col6">6.8</oasis:entry>  
         <oasis:entry colname="col7">1.41</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">16 Feb, a.m.<inline-formula><mml:math id="M133" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.32</oasis:entry>  
         <oasis:entry colname="col3">0.19</oasis:entry>  
         <oasis:entry colname="col4">0.5</oasis:entry>  
         <oasis:entry colname="col5">1.2</oasis:entry>  
         <oasis:entry colname="col6">3.92</oasis:entry>  
         <oasis:entry colname="col7">0.75</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">16 Feb, p.m.<inline-formula><mml:math id="M134" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.32</oasis:entry>  
         <oasis:entry colname="col3">0.19</oasis:entry>  
         <oasis:entry colname="col4">0.4</oasis:entry>  
         <oasis:entry colname="col5">0.96</oasis:entry>  
         <oasis:entry colname="col6">5.75</oasis:entry>  
         <oasis:entry colname="col7">1.09</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">17 Feb, a.m.<inline-formula><mml:math id="M135" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.31</oasis:entry>  
         <oasis:entry colname="col3">0.194</oasis:entry>  
         <oasis:entry colname="col4">0.5</oasis:entry>  
         <oasis:entry colname="col5">1.17</oasis:entry>  
         <oasis:entry colname="col6">4.4</oasis:entry>  
         <oasis:entry colname="col7">0.85</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">17 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.31</oasis:entry>  
         <oasis:entry colname="col3">0.194</oasis:entry>  
         <oasis:entry colname="col4">0.5</oasis:entry>  
         <oasis:entry colname="col5">1.17</oasis:entry>  
         <oasis:entry colname="col6">4.4</oasis:entry>  
         <oasis:entry colname="col7">0.85</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">18 Feb, a.m.<inline-formula><mml:math id="M136" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">1</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.29</oasis:entry>  
         <oasis:entry colname="col3">0.199</oasis:entry>  
         <oasis:entry colname="col4">0.3</oasis:entry>  
         <oasis:entry colname="col5">0.68</oasis:entry>  
         <oasis:entry colname="col6">7.53</oasis:entry>  
         <oasis:entry colname="col7">1.5</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">18 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.29</oasis:entry>  
         <oasis:entry colname="col3">0.199</oasis:entry>  
         <oasis:entry colname="col4">0.8</oasis:entry>  
         <oasis:entry colname="col5">1.81</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M137" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1.47</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M138" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.29</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">19 Feb, a.m.<inline-formula><mml:math id="M139" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">1</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.25</oasis:entry>  
         <oasis:entry colname="col3">0.199</oasis:entry>  
         <oasis:entry colname="col4">0</oasis:entry>  
         <oasis:entry colname="col5">0</oasis:entry>  
         <oasis:entry colname="col6">9</oasis:entry>  
         <oasis:entry colname="col7">1.79</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">19 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.25</oasis:entry>  
         <oasis:entry colname="col3">0.199</oasis:entry>  
         <oasis:entry colname="col4">0.8</oasis:entry>  
         <oasis:entry colname="col5">1.81</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M140" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1.47</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M141" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.29</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">20 Feb, a.m.</oasis:entry>  
         <oasis:entry colname="col2">0.21</oasis:entry>  
         <oasis:entry colname="col3">0.161</oasis:entry>  
         <oasis:entry colname="col4">0.5</oasis:entry>  
         <oasis:entry colname="col5">1.61</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M142" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3.69</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M143" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.59</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">20 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.21</oasis:entry>  
         <oasis:entry colname="col3">0.161</oasis:entry>  
         <oasis:entry colname="col4">1.2</oasis:entry>  
         <oasis:entry colname="col5">3.87</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M144" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">64.1</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M145" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">10.3</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">21 Feb, a.m.</oasis:entry>  
         <oasis:entry colname="col2">0.19</oasis:entry>  
         <oasis:entry colname="col3">0.167</oasis:entry>  
         <oasis:entry colname="col4">1.5</oasis:entry>  
         <oasis:entry colname="col5">4.48</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M146" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">82</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M147" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">13.7</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">21 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.19</oasis:entry>  
         <oasis:entry colname="col3">0.167</oasis:entry>  
         <oasis:entry colname="col4">2</oasis:entry>  
         <oasis:entry colname="col5">5.98</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M148" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">152</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M149" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">25.5</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">22 Feb, a.m.</oasis:entry>  
         <oasis:entry colname="col2">0.15</oasis:entry>  
         <oasis:entry colname="col3">0.198</oasis:entry>  
         <oasis:entry colname="col4">1.4</oasis:entry>  
         <oasis:entry colname="col5">3.19</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M150" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">23.8</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M151" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">4.71</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">22 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.15</oasis:entry>  
         <oasis:entry colname="col3">0.198</oasis:entry>  
         <oasis:entry colname="col4">0.9</oasis:entry>  
         <oasis:entry colname="col5">2.05</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M152" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">4.56</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M153" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.9</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">23 Feb, a.m.<inline-formula><mml:math id="M154" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.12</oasis:entry>  
         <oasis:entry colname="col3">0.132</oasis:entry>  
         <oasis:entry colname="col4">0.1</oasis:entry>  
         <oasis:entry colname="col5">0.63</oasis:entry>  
         <oasis:entry colname="col6">6.1</oasis:entry>  
         <oasis:entry colname="col7">0.81</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">23 Feb, p.m.<inline-formula><mml:math id="M155" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">1</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.12</oasis:entry>  
         <oasis:entry colname="col3">0.132</oasis:entry>  
         <oasis:entry colname="col4">0.2</oasis:entry>  
         <oasis:entry colname="col5">1.26</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M156" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2.58</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M157" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.34</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">24 Feb, a.m.</oasis:entry>  
         <oasis:entry colname="col2">0.07</oasis:entry>  
         <oasis:entry colname="col3">0.223</oasis:entry>  
         <oasis:entry colname="col4">1.8</oasis:entry>  
         <oasis:entry colname="col5">3.59</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M158" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">23.7</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M159" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">5.28</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">24 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.07</oasis:entry>  
         <oasis:entry colname="col3">0.223</oasis:entry>  
         <oasis:entry colname="col4">1.3</oasis:entry>  
         <oasis:entry colname="col5">2.59</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M160" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">8.06</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M161" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1.8</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">25 Feb, a.m.</oasis:entry>  
         <oasis:entry colname="col2">0.02</oasis:entry>  
         <oasis:entry colname="col3">1.29</oasis:entry>  
         <oasis:entry colname="col4">1.3</oasis:entry>  
         <oasis:entry colname="col5">1.42</oasis:entry>  
         <oasis:entry colname="col6">8.85</oasis:entry>  
         <oasis:entry colname="col7">11.43</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">25 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.02</oasis:entry>  
         <oasis:entry colname="col3">1.29</oasis:entry>  
         <oasis:entry colname="col4">0</oasis:entry>  
         <oasis:entry colname="col5">0</oasis:entry>  
         <oasis:entry colname="col6">9</oasis:entry>  
         <oasis:entry colname="col7">11.63</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">26 Feb, a.m.</oasis:entry>  
         <oasis:entry colname="col2">0.1</oasis:entry>  
         <oasis:entry colname="col3">0.399</oasis:entry>  
         <oasis:entry colname="col4">1.6</oasis:entry>  
         <oasis:entry colname="col5">2.22</oasis:entry>  
         <oasis:entry colname="col6">5.1</oasis:entry>  
         <oasis:entry colname="col7">2.03</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">26 Feb, p.m.</oasis:entry>  
         <oasis:entry colname="col2">0.1</oasis:entry>  
         <oasis:entry colname="col3">0.399</oasis:entry>  
         <oasis:entry colname="col4">0.2</oasis:entry>  
         <oasis:entry colname="col5">0.28</oasis:entry>  
         <oasis:entry colname="col6">8.94</oasis:entry>  
         <oasis:entry colname="col7">3.6</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">27 Feb, a.m.<inline-formula><mml:math id="M162" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.19</oasis:entry>  
         <oasis:entry colname="col3">0.273</oasis:entry>  
         <oasis:entry colname="col4">1.2</oasis:entry>  
         <oasis:entry colname="col5">2.02</oasis:entry>  
         <oasis:entry colname="col6">2.06</oasis:entry>  
         <oasis:entry colname="col7">0.56</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">27 Feb, p.m.<inline-formula><mml:math id="M163" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">1</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.19</oasis:entry>  
         <oasis:entry colname="col3">0.273</oasis:entry>  
         <oasis:entry colname="col4">0.6</oasis:entry>  
         <oasis:entry colname="col5">1.01</oasis:entry>  
         <oasis:entry colname="col6">7.26</oasis:entry>  
         <oasis:entry colname="col7">1.98</oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">28 Feb, a.m.</oasis:entry>  
         <oasis:entry colname="col2">0.36</oasis:entry>  
         <oasis:entry colname="col3">0.224</oasis:entry>  
         <oasis:entry colname="col4">1.1</oasis:entry>  
         <oasis:entry colname="col5">2.18</oasis:entry>  
         <oasis:entry colname="col6"><inline-formula><mml:math id="M164" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col7"><inline-formula><mml:math id="M165" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.67</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>  
         <oasis:entry colname="col1">28 Feb, p.m.<inline-formula><mml:math id="M166" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">1</mml:mn></mml:msup></mml:math></inline-formula></oasis:entry>  
         <oasis:entry colname="col2">0.36</oasis:entry>  
         <oasis:entry colname="col3">0.224</oasis:entry>  
         <oasis:entry colname="col4">0.4</oasis:entry>  
         <oasis:entry colname="col5">0.79</oasis:entry>  
         <oasis:entry colname="col6">7.41</oasis:entry>  
         <oasis:entry colname="col7">1.66</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup><?xmltex \end{scaleboxenv}?></oasis:table><table-wrap-foot><p><inline-formula><mml:math id="M119" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">1</mml:mn></mml:msup></mml:math></inline-formula> Tides in which a double high water is
observed. <inline-formula><mml:math id="M120" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:math></inline-formula> Tides where there is a stand in the
tide, so that a double high water is close to being formed.</p></table-wrap-foot></table-wrap>

</sec>
<sec id="Ch1.S3.SS2">
  <title>Testing the condition for double high waters</title>
      <p>A full analysis of the conditions that produce a double high water would
include all of the relevant harmonics of the semi-diurnal tide. In the case
of the short record at Port Ellen, however, the mean amplitude of D4 is
0.026 m and that of D6 is 0.044 m. The fact that D6 is generally larger,
coupled with the fact that it is of higher frequency and so more potent at
producing double high tides, means that, if we are to consider just a single
harmonic, then that should be D6. The actual advantage for D6 is
<inline-formula><mml:math id="M167" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mn mathvariant="normal">9</mml:mn><mml:mo>/</mml:mo><mml:mn mathvariant="normal">4</mml:mn><mml:mo>)</mml:mo><mml:mo>×</mml:mo><mml:mo>(</mml:mo><mml:mn mathvariant="normal">44</mml:mn><mml:mo>/</mml:mo><mml:mn mathvariant="normal">26</mml:mn><mml:mo>)</mml:mo><mml:mo>∼</mml:mo><mml:mn mathvariant="normal">3.8</mml:mn></mml:mrow></mml:math></inline-formula>, so we limit our analysis in this section to
the production of double high waters by a combination of D6 and D2 tides, and
use this combination to test the theory of Sect. 2. We bear in mind,
however, that the comparison between observations and theory may not be exact
because we are neglecting other harmonics, notably D4 and also D1. It can be
seen in Table 1 that on all days the ratio of amplitudes <inline-formula><mml:math id="M168" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> for D6 and D2
tides is greater than the value of 0.111 required to satisfy the condition (1).
Moreover, the value of <inline-formula><mml:math id="M169" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> is nearly as great on the afternoon of
24 February (0.223, when no double high water is seen) as it is on the afternoon of
the 28 February (0.224, when there is a double high water). The straightforward
Doodson criterion, neglecting the effect of phase difference between the
harmonic and the principal tide, is therefore not the best discriminator for
double high waters. The rows in Table 1 have been marked for tides where a
double high water is observed, or nearly so. Tides in which a double high
water is observed are shaded in the darker grey and those where there is a
stand in the tide, so that a double high water is close to being formed are
marked in lighter grey. On days when there is a clear double high water, that
is on the morning of 18 and 19 February and on the afternoons of 27 and
28 February, the value of <inline-formula><mml:math id="M170" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> is greater than or equal to 1.5. On days when
there is a stand in the tide close to high water, <inline-formula><mml:math id="M171" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> is greater than
about 0.5 and less than 1.5. On all other days, when there is no double high
water (excepting 25 and 26 February) the value of <inline-formula><mml:math id="M172" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> is less than 0.5.
The amended Doodson condition therefore performs well as a predictor of
double tides, although the transitions are not as clear-cut in the
observations as theory would predict. This is likely to be partly because of
difficulty in deciding when a double high water is there by visual inspection
alone and partly because the real Port Ellen tide has harmonics other than
the D2 and D6 components treated in the theory, as well as other sea-level
processes and sampling noise being present.</p>
      <p>On 25 and 26 February, <inline-formula><mml:math id="M173" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> reaches its highest values,
over 11 on 25 February and a more modest level of over 2 on 26 February. There is
no clear double high water on these days, though. Instead the tide is
confused with no clear semi-diurnal pattern (Fig. 4). Closer inspection of
the tides on these days (using plots equivalent to Fig. <xref ref-type="fig" rid="Ch1.F5"/>)
shows that on 25 February the diurnal and quarter diurnal tides are dominant –
there is certainly no double high water in the semi-diurnal tide on this day.
The situation on  26 February is more complicated. The diurnal tide is the
largest component on this day. If the diurnal tide is removed from the
reconstruction and a tidal curve is drawn as the sum of the D2, D4, and D6
harmonics, then a double high water can be seen in both the morning and
afternoon. The predictor is therefore doing a reasonable job, within its
limitations, throughout the fortnight of the observations. When the parameter
<inline-formula><mml:math id="M174" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> has a value less than about 1, no double high water is formed
because the higher harmonic is too weak. When it has a value greater than
about 10, then again no double high water is formed because the higher
harmonic is too strong. Between these limiting values, double high waters, or
at least a tidal stand, are likely to be observed, provided they are not
obscured by the presence of other harmonics.</p>
      <p>We have plotted in Fig. 2 values of <inline-formula><mml:math id="M175" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> for the D2 and D6 tides (from
Table 1) against the phase difference <inline-formula><mml:math id="M176" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> (also from Table 1). The points
have been coded so that those on tides in which a double high tide is
observed are shown as filled circles and those on which no double high tide
is observed are shown as open circles. Tides where there is a point of
inflection, or a stand in the tide, have been marked as grey-filled circles.
The tides of 25 and 26 February, when the semi-diurnal tide is very small,
have been omitted from this diagram. Because Fig. 2 shows a graphical
solution to the inequality <inline-formula><mml:math id="M177" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> the inequality is satisfied and we
expect a double high water to form if we are above the curves (incidentally,
the same is true for double low waters; <inline-formula><mml:math id="M178" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> then becomes the time
difference between the D2 low and the Dn high). This
is very much the case – all four clear double high waters plot above the
theoretical transition curve. Similarly, points representing the tides with
no double high water lie below the critical curve. The grey points,
representing tides with a stand, lie close to the critical lines. It is
illuminating to note that the main spread of the points in Fig. 2 is along
the <inline-formula><mml:math id="M179" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> axis. Because the amplitude of harmonics generated by the
semi-diurnal tide tend to increase with that of D2, the ratio <inline-formula><mml:math id="M180" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> remains
fairly constant, and Fig. <xref ref-type="fig" rid="Ch1.F6"/> confirms that D6 depends on D2 but
D4 does not. The critical factor in deciding whether a double high tide will
form is actually the phase difference. As the time of the dip in D6 moves
close to the high tide in D2, the critical condition for a double high tide
is met. Looking at Fig. 2, most of the spread of the points is along the
<inline-formula><mml:math id="M181" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> axis – the variation of <inline-formula><mml:math id="M182" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> is relatively small (because <inline-formula><mml:math id="M183" display="inline"><mml:mi>b</mml:mi></mml:math></inline-formula> is
proportional to <inline-formula><mml:math id="M184" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> for D6 and D2). Consequently, what brings a point across
the critical line, in practice, is a change in the phase of the constituents.
Although the theory has it that both phase and amplitude ratio are important,
in practice for this data set, the phase difference between the harmonics is
the most important parameter in controlling the formation of a double tide.</p>
</sec>
</sec>
<sec id="Ch1.S4">
  <title>Discussion</title>
      <p>The theoretical considerations presented in this paper, supported by a small
data set, suggest that a single parameter can be used to predict the presence
of a double high water when a higher harmonic is added to a semi-diurnal
tide. As we might expect, this parameter – <inline-formula><mml:math id="M185" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> – depends on the
amplitude and phase of the harmonic (compared to the semi-diurnal tide) and
on the ratio of frequencies of harmonic and main tide. The data in Table 1
can be divided into four categories:
<list list-type="bullet"><list-item>
      <p><inline-formula><mml:math id="M186" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>≤</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>: no double high water is seen.</p></list-item><list-item>
      <p><inline-formula><mml:math id="M187" display="inline"><mml:mrow><mml:mn mathvariant="normal">0.5</mml:mn><mml:mo>≤</mml:mo><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>≤</mml:mo><mml:mo>∼</mml:mo><mml:mn mathvariant="normal">1.5</mml:mn></mml:mrow></mml:math></inline-formula>: a stand in the tide is observed,
but no clear double high water.</p></list-item><list-item>
      <p><inline-formula><mml:math id="M188" display="inline"><mml:mrow><mml:mn mathvariant="normal">1.5</mml:mn><mml:mo>≤</mml:mo><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>≤</mml:mo><mml:mo>∼</mml:mo><mml:mn mathvariant="normal">10</mml:mn></mml:mrow></mml:math></inline-formula>: the regime of double high waters.</p></list-item><list-item>
      <p><inline-formula><mml:math id="M189" display="inline"><mml:mrow><mml:mo>(</mml:mo><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">10</mml:mn></mml:mrow></mml:math></inline-formula>: the harmonic dominates and there is again no
clear double high tide.</p></list-item></list></p>
      <p>In the terms posed in the introduction, it is necessary to place both lower
and upper bounds on the criterion for a double high water to allow for the
higher harmonic being too small and too large.</p>
      <p>A limitation of the theory, as presented here, is that it considers just a
single higher harmonic added to the semi-diurnal tide. There will be places,
and times, when this is appropriate. In the case of the data from Port Ellen
we present here, the theory adequately represents the data for most of the
time. The theory can be extended to include other harmonics in a
straightforward way. The condition for the initial formation of a double high
water by two harmonics added to the semi-diurnal tide, for example, can be
written <xref ref-type="bibr" rid="bib1.bibx3" id="paren.19"><named-content content-type="pre">note that</named-content><named-content content-type="post"> discussed a triple-harmonic case, but for
amplitudes only</named-content></xref>:
          <disp-formula id="Ch1.E11" content-type="numbered"><mml:math id="M190" display="block"><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow><mml:mi>a</mml:mi></mml:mfrac></mml:mstyle><mml:msubsup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow><mml:mi>a</mml:mi></mml:mfrac></mml:mstyle><mml:msubsup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>
      <p>Here, <inline-formula><mml:math id="M191" display="inline"><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mn mathvariant="normal">1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M192" display="inline"><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> are the amplitudes of the two harmonics and <inline-formula><mml:math id="M193" display="inline"><mml:mrow><mml:msubsup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> and
<inline-formula><mml:math id="M194" display="inline"><mml:mrow><mml:msubsup><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> are calculated from Eq. (<xref ref-type="disp-formula" rid="Ch1.E7"/>) with the appropriate phases. The
time at the centre of the dip, <inline-formula><mml:math id="M195" display="inline"><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>, is calculated from a modified form of
Eq. (<xref ref-type="disp-formula" rid="Ch1.E5"/>) which allows for both harmonics acting together. It can
happen that the two harmonics support each other in producing a double high
water or (because the <inline-formula><mml:math id="M196" display="inline"><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> parameter can be negative) that one of the
harmonics suppresses a double high water that would otherwise be formed by a
single harmonic acting on its own.</p>
      <p>Turning now to the observations, it is, in fact, difficult to tell, in
marginal cases, when there is a double high water by visual inspection alone.
This is because, when the double high first forms, the dip between the two
high waters is small – perhaps just a centimetre or so – and this is often at
the noise level of the measurements. In this paper, we have tried to avoid
this problem by referring to a tidal stand – that is, there is clearly a
portion of the tidal record where the water level remains flat for a time,
but it is difficult to say for sure if there is a double high tide.</p>
      <p>For most of the data set, the two most important harmonics at Port Ellen are
the semi-diurnal (D2) and the sixth-diurnal D6. Other harmonics, in
particular D1 and D4 play a role, however. The effect of the diurnal
harmonic, D1, in allowing a double high water in one tide but not in the
second on the same day has already been mentioned. The role of the D4
harmonic at Port Ellen is to suppress the formation of the double low tide,
although it is quite likely that there will be a double low near Port Ellen
<xref ref-type="bibr" rid="bib1.bibx12" id="paren.20"><named-content content-type="pre">cf. Southampton–Weymouth;</named-content></xref>. The sum of D2 and
D6 harmonics tends to produce a double high tide and a double low tide
(Fig. 1a). At Port Ellen, the minimum in D4 occurs at about the time of the
minimum in D2 (see Fig. <xref ref-type="fig" rid="Ch1.F5"/>b). The D4 harmonic, though small,
tends to flatten out the low tide and prevents the sixth diurnal harmonic
(which has a maximum at this time) from producing a double low water. The
opposite is true at high tide, when the minimum in D4 helps the minimum in D6
to produce a double high water.</p>
</sec>
<sec id="Ch1.S5" sec-type="conclusions">
  <title>Conclusions</title>
      <p>The theory presented here has extended the original Doodson criterion and now
includes the phase, <inline-formula><mml:math id="M197" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula>, and it allows for <inline-formula><mml:math id="M198" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> to be large. We have,
however, not covered all situations, and we do not include phases larger than
2 h in Fig. 2. We also show that the phase can be the most important
variable in producing double high waters. That is, <inline-formula><mml:math id="M199" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>/</mml:mo><mml:mi>a</mml:mi></mml:mrow></mml:math></inline-formula> can be large enough
but there is no double high because the phase prevents its formation. This
has been overlooked in previous papers on the subject. Also, we have focused
on double highs in the paper, but the theory is just as applicable for double
lows, where the phase goes from <inline-formula><mml:math id="M200" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula> to <inline-formula><mml:math id="M201" display="inline"><mml:mrow><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">180</mml:mn></mml:mrow></mml:math></inline-formula><inline-formula><mml:math id="M202" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>. There is most
likely a double low in the vicinity of Port Ellen, but locating and
quantifying it is left for a future study.</p>
      <p>Overall, the new theory captures the double highs at Port Ellen very well. It
shows that it is the D6 harmonic which is the dominant one at our location,
rather than the usual D4. The formation of double high and double low waters
in the semi-diurnal tide is a fascinating problem with important practical
implications, e.g. to understand the differences between mean tide level and
mean sea level <xref ref-type="bibr" rid="bib1.bibx12" id="paren.21"><named-content content-type="pre">see the discussion in</named-content></xref>.</p>
</sec>

      
      </body>
    <back><notes notes-type="dataavailability">

      <p>The sea-level data were provided by the UK National Tide
Gauge Network through the National Tidal and Sea Level Facility
(<uri>http://www.ntslf.org/</uri>), sponsored by the UK Environment Agency. The
data for Port Ellen can be downloaded from the BODC archive at
<uri>https://www.bodc.ac.uk/data/hosted_data_systems/sea_level/uk_tide_gauge_network/processed/</uri>.</p>
  </notes><?xmltex \hack{\newpage}?><notes notes-type="competinginterests">

      <p>The authors declare that they have no conflict of
interest.</p>
  </notes><ack><title>Acknowledgements</title><p>Funding was provided by the Natural Environmental Research Council (grant
NE/I030224/1 to JAMG), and from the BurningHam research support foundation.
Constructive comments from Phil Woodworth and two anonymous reviewers greatly
improved the manuscript.<?xmltex \hack{\newline}?><?xmltex \hack{\newline}?> Edited by:
Neil Wells<?xmltex \hack{\newline}?> Reviewed by: Philip Woodworth and two anonymous
referees</p></ack><ref-list>
    <title>References</title>

      <ref id="bib1.bibx1"><label>Airy(1843)</label><mixed-citation>
Airy, G. B.: On the laws of individual tides at Southampton and at Ipswich,
Philos. T. R. Soc. Lond., 133, 45–54, 1843.</mixed-citation></ref>
      <ref id="bib1.bibx2"><label>Bowers et al.(2013)</label><mixed-citation>
Bowers, D., Macdonald, R., McKee, D., Nimmo-Smith, W., and Graham, G.: On the
formation of tide-produced seiches and double high waters in coastal seas,
Estuarine, Coast. Shelf Sci., 134, 108–116, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx3"><label>Doodson and Warburg(1941)</label><mixed-citation>Doodson, A. T. and Warburg, H.: Admiralty manual of tides, Her Majesty's
Stationery Office, London, 1941.  </mixed-citation></ref><?xmltex \hack{\newpage}?>
      <ref id="bib1.bibx4"><label>Emery and Thomson(1996)</label><mixed-citation>
Emery, W. J. and Thomson, R. E.: Data analysis methods in physical
oceanography, 2nd Edn., Pergamon Press, 650 pp., 1996.</mixed-citation></ref>
      <ref id="bib1.bibx5"><label>Godin(1993)</label><mixed-citation>
Godin, G.: An investigation of the phenomenon of double high water or double
low water at some harbours, Deutsche Hydrographische Zeitschrift, 45,
87–106, 1993.</mixed-citation></ref>
      <ref id="bib1.bibx6"><label>Nunes and Lennon(1986)</label><mixed-citation>
Nunes, R. and Lennon, G.: Physical property distributions and seasonal trends
in Spencer Gulf, South Australia, an inverse estuary, Mar. Freshwater Res., 37, 39–53, 1986.</mixed-citation></ref>
      <ref id="bib1.bibx7"><label>Pugh(1982)</label><mixed-citation>
Pugh, D.: A comparison of recent and historical tides and mean sea-levels off
Ireland, Geophys. J. Roy. Astr. S., 71, 809–815, 1982.</mixed-citation></ref>
      <ref id="bib1.bibx8"><label>Pugh(1987)</label><mixed-citation>
Pugh, D. T.: Tides, surges and mean sea level: a handbook for engineers and
scientists, Wiley, Chichester, UK, 1987.</mixed-citation></ref>
      <ref id="bib1.bibx9"><label>Pugh and Woodworth(2014)</label><mixed-citation>
Pugh, D. and Woodworth, P.: Sea-Level Science: Understanding tides, surges,
tsunamis and mean sea-level changes, Cambridge University Press, Cambridge,
UK, 2014.</mixed-citation></ref>
      <ref id="bib1.bibx10"><label>Redfield(1980)</label><mixed-citation>
Redfield, A. C.: The tides of New England and New York, Woods Hole
Oceanographic Institution, USA, 1980.</mixed-citation></ref>
      <ref id="bib1.bibx11"><label>Weatherall et al.(2015)</label><mixed-citation>
Weatherall P., Marks, K. M., Jakobsson, M.,
Schmitt, T., Tani, S., Arndt, J. E., Rovere, M., Chayes, D., Ferrini, V., and Wigley,
R.: A new digital bathymetric model
of the world's oceans, Earth Space
Sci., 2, 331–345, doi:10.1002/2015EA000107, 2015.</mixed-citation></ref>
      <ref id="bib1.bibx12"><label>Woodworth(2017)</label><mixed-citation>
Woodworth, P. L.: Differences between mean tide level and mean sea level,
J. Geodesy, 91, 69–90, 2017.</mixed-citation></ref>

  </ref-list><app-group content-type="float"><app><title/>

    </app></app-group></back>
    <!--<article-title-html>The double high tide at Port Ellen: Doodson's criterion revisited</article-title-html>
<abstract-html><p class="p">Doodson proposed a minimum criterion to predict the occurrence of double
high (or double low) waters when a higher-frequency tidal harmonic is added
to the semi-diurnal tide. If the phasing of the harmonic is optimal, the
condition for a double high water can be written <i>b</i><i>n</i><sup>2</sup>∕<i>a</i> &gt; 1 where <i>b</i> is
the amplitude of the higher harmonic, <i>a</i> is the amplitude of the
semi-diurnal tide, and <i>n</i> is the ratio of their frequencies. Here we expand
this criterion to allow for (i) a phase difference <i>ϕ</i> between the
semi-diurnal tide and the harmonic and (ii) the fact that the double high
water will disappear in the event that <i>b</i>∕<i>a</i> becomes large enough for the
higher harmonic to be the dominant component of the tide. This can happen,
for example, at places or times where the semi-diurnal tide is very small.
The revised parameter is <i>b</i><i>r</i><sup>2</sup>∕<i>a</i>, where <i>r</i> is a number generally less
than <i>n</i>, although equal to <i>n</i> when <i>ϕ</i> = 0. The theory predicts that a
double high tide will form when this parameter exceeds 1 and then disappear
when it exceeds a value of order <i>n</i><sup>2</sup> and the higher harmonic becomes
dominant. We test these predictions against observations at Port Ellen in
the Inner Hebrides of Scotland. For most of the data set, the largest
harmonic of the semi-diurnal tide is the sixth diurnal component, for which
<i>n</i> = 3. The principal lunar and solar semi-diurnal tides are about equal at
Port Ellen and so the semi-diurnal tide becomes very small twice a month at
neap tides (here defined as the smallest fortnightly tidal range). A double
high water forms when <i>b</i><i>r</i><sup>2</sup>∕<i>a</i> first exceeds a minimum value of about 1.5
as neap tides are approached and then disappears as <i>b</i><i>r</i><sup>2</sup>∕<i>a</i> then exceeds a
second limiting value of about 10 at neap tides in agreement with the
revised criterion.</p></abstract-html>
<ref-html id="bib1.bib1"><label>Airy(1843)</label><mixed-citation>
Airy, G. B.: On the laws of individual tides at Southampton and at Ipswich,
Philos. T. R. Soc. Lond., 133, 45–54, 1843.
</mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>Bowers et al.(2013)</label><mixed-citation>
Bowers, D., Macdonald, R., McKee, D., Nimmo-Smith, W., and Graham, G.: On the
formation of tide-produced seiches and double high waters in coastal seas,
Estuarine, Coast. Shelf Sci., 134, 108–116, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>Doodson and Warburg(1941)</label><mixed-citation>
Doodson, A. T. and Warburg, H.: Admiralty manual of tides, Her Majesty's
Stationery Office, London, 1941. 
</mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>Emery and Thomson(1996)</label><mixed-citation>
Emery, W. J. and Thomson, R. E.: Data analysis methods in physical
oceanography, 2nd Edn., Pergamon Press, 650 pp., 1996.
</mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>Godin(1993)</label><mixed-citation>
Godin, G.: An investigation of the phenomenon of double high water or double
low water at some harbours, Deutsche Hydrographische Zeitschrift, 45,
87–106, 1993.
</mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>Nunes and Lennon(1986)</label><mixed-citation>
Nunes, R. and Lennon, G.: Physical property distributions and seasonal trends
in Spencer Gulf, South Australia, an inverse estuary, Mar. Freshwater Res., 37, 39–53, 1986.
</mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>Pugh(1982)</label><mixed-citation>
Pugh, D.: A comparison of recent and historical tides and mean sea-levels off
Ireland, Geophys. J. Roy. Astr. S., 71, 809–815, 1982.
</mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>Pugh(1987)</label><mixed-citation>
Pugh, D. T.: Tides, surges and mean sea level: a handbook for engineers and
scientists, Wiley, Chichester, UK, 1987.
</mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>Pugh and Woodworth(2014)</label><mixed-citation>
Pugh, D. and Woodworth, P.: Sea-Level Science: Understanding tides, surges,
tsunamis and mean sea-level changes, Cambridge University Press, Cambridge,
UK, 2014.
</mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>Redfield(1980)</label><mixed-citation>
Redfield, A. C.: The tides of New England and New York, Woods Hole
Oceanographic Institution, USA, 1980.
</mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>Weatherall et al.(2015)</label><mixed-citation>
Weatherall P., Marks, K. M., Jakobsson, M.,
Schmitt, T., Tani, S., Arndt, J. E., Rovere, M., Chayes, D., Ferrini, V., and Wigley,
R.: A new digital bathymetric model
of the world's oceans, Earth Space
Sci., 2, 331–345, doi:10.1002/2015EA000107, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>Woodworth(2017)</label><mixed-citation>
Woodworth, P. L.: Differences between mean tide level and mean sea level,
J. Geodesy, 91, 69–90, 2017.
</mixed-citation></ref-html>--></article>
