Surabaya, Jawa Timur
Mixed, mainly semidiurnal- Period
- 01 Jan 2026 — 01 Agu 2026212 days
- Datum
- Nol sensor stasiun (bukan MSL)
- gaps
- 7891 jam
- Source
- IOC Sea Level Station Monitoring Facility (UNESCO/IOC & VLIZ)
Source, licence and datum notes
IOC menyajikan data mentah terhadap nol sensor masing-masing stasiun. Tinggi di sini tidak merujuk MSL, LAT, maupun chart datum.
Licence: Akses terbuka dengan kewajiban sitasi (VLIZ/IOC, DOI 10.14284/482)
Flanders Marine Institute (VLIZ); Intergovernmental Oceanographic Commission (IOC) (2026): Sea level station monitoring facility. https://www.ioc-sealevelmonitoring.org — DOI 10.14284/482
Not for navigationThe official Indonesian tide tables are published by Pushidrosal.
The chart
The record
- Observed the water level the gauge actually recorded
- Predicted what the fitted constituents say it should be
- Residual the difference — weather, wind and instrument trouble
The rest is held out and predicted. Shorten the window and constituents start dropping out — the same Rayleigh criterion, moved by the window rather than by the length of the record.
Result
Request refused
S2 and K2 cannot be separated on a record of 141.3 days.
- S2 and K2 differ by only 0.0821°/h: separating them needs a record of 182.6 days, and 141.3 are available.
- K1 and P1 differ by only 0.0821°/h: separating them needs a record of 182.6 days, and 141.3 are available.
Record length: 141.3 days · Length required: 182.6 days
What follows is the largest set this window can still support: M2, S2, N2, K1, O1, Q1, M4, MS4.
- Condition number κ
- 1.44good
- Residual RMS
- 0.1498 m
- Mean level Z₀
- 1.7488 m
- Record length
- 141.3 days
Unexplained, over the fitted part
0.1498m
Mean distance between the black line and the blue one, over the part used to fit.
Unexplained, over the held-out part
0.1579m
Over the part deliberately not seen while fitting. This is the real test.
Constituents
Harmonic constituents
| Constituent | Amplitude H (m) | Phase g (°) | Lag (h) | Speed (°/h) | Period (h) | Nodal factor f | Nodal correction u (°) |
|---|---|---|---|---|---|---|---|
| K1165.555 | 0.4554±0.0033 | 25.1±0.4 | 1.67 | 15.041069 | 23.934 | 1.1068 | 2.86 |
| M2255.555 | 0.3786±0.0038 | 118.5±0.6 | 4.09 | 28.984104 | 12.421 | 0.9659 | 0.79 |
| O1145.555 | 0.2750±0.0031 | 346.1±0.7 | 24.82 | 13.943036 | 25.819 | 1.1728 | -3.23 |
| S2273.555 | 0.2538±0.0037 | 132.2±0.8 | 4.41 | 30.000000 | 12.000 | 1.0000 | 0.00 |
| N2245.655 | 0.0659±0.0038 | 107.3±3.3 | 3.77 | 28.439730 | 12.658 | 0.9659 | 0.79 |
| Q1135.655 | 0.0509±0.0031 | 335.8±3.5 | 25.06 | 13.398661 | 26.868 | 1.1728 | -3.23 |
| MS4473.555 | 0.0223±0.0038 | 196.0±9.8 | 3.32 | 58.984104 | 6.103 | 0.9659 | 0.79 |
| M4455.555 | 0.0211±0.0039 | 199.0±10.7 | 3.43 | 57.968208 | 6.210 | 0.9329 | 1.57 |
H is half the height of that constituent’s wave; g is how late it arrives behind the Moon or Sun, and the column beside it says the same thing in hours — g divided by the constituent’s speed. f and u are the 18.6-year corrections, applied and still shown.
The ± figure is one standard error, taken from the diagonal of the solve’s covariance matrix. It assumes the record’s noise is independent from hour to hour, and tidal residuals are not: weather lasts for days, so neighbouring hours are wrong in the same direction. The true uncertainty is larger than what is printed here — by how much, this site does not compute.
A property of the place, or of these months?
The same numbers, four different stretches of the record
This whole site rests on one claim: that amplitude and phase are properties of the place — what a coastline does to the same forcing. But every number above came from one particular stretch of 2026. Here the record is cut into 4 equal, non-overlapping pieces and each is fitted on its own. Where an amplitude barely moves between them it is behaving like a property of the harbour. Where it swings, the record is doing the work rather than the place.
- Stretch 1 01 Jan 2026 — 23 Feb 2026
- Stretch 2 23 Feb 2026 — 17 Apr 2026
- Stretch 3 17 Apr 2026 — 09 Jun 2026
- Stretch 4 09 Jun 2026 — 01 Agu 2026
| Constituent | Mean H (m) | Lowest — highest | Swing | Phase swing |
|---|---|---|---|---|
| K1 | 0.4903 | 0.3796 — 0.5491 | 35% | ±13.8° |
| M2 | 0.3781 | 0.3712 — 0.3854 | 4% | ±2.0° |
| O1 | 0.2733 | 0.2671 — 0.2782 | 4% | ±4.0° |
| S2 | 0.2258 | 0.1493 — 0.3024 | 68% | ±23.1° |
| N2 | 0.0741 | 0.0661 — 0.0813 | 21% | ±15.0° |
| Q1 | 0.0508 | 0.0427 — 0.0569 | 28% | ±3.1° |
| MS4 | 0.0227 | 0.0072 — 0.0370 | 131% | ±27.8° |
| M4 | 0.0202 | 0.0182 — 0.0219 | 18% | ±9.2° |
One constituent, from the start
How K1 came out of this record
The table above is the answer. This is the route to it, in this record’s own numbers — nothing simplified or rounded for the example.
- Step 1 — Where the Moon and Sun are. Every constituent has six Doodson numbers saying how many times each astronomical element enters its angle. For K1 they are 1 1 0 0 0 0, so its equilibrium argument is V(t) = τ + s − 90°. The element polynomials are from Meeus, Astronomical Algorithms (2nd ed.), chapters 22 and 47.
τ = mean lunar time · s = the Moon’s mean longitude · h = the Sun’s mean longitude · p = lunar perigee · N = the Moon’s ascending node · p′ = solar perigee.
Step 2 — Two columns, hour by hour. cos(V + u) and sin(V + u) are evaluated at every sample in the record. These are K1’s two columns of the design matrix; four consecutive hours of this record:
How K1 came out of this record Time (UTC) V(t)° cos(V+u) sin(V+u) K1’s share (m) Recorded (m) 12 Mar 2026 12:00 UTC 80.15 0.1217 0.9926 +0.2678 2.0620 12 Mar 2026 13:00 UTC 95.20 -0.1401 0.9901 +0.1479 1.9670 12 Mar 2026 14:00 UTC 110.24 -0.3922 0.9199 +0.0177 1.7650 12 Mar 2026 15:00 UTC 125.28 -0.6175 0.7866 -0.1136 1.5580 - Step 3 — Least squares picks one pair of numbers. The solve looks for the a and b that bring Σ (a·cos + b·sin) as close to the record as it can get, for every constituent at once. For K1 it returned a = 0.4564 and b = 0.2139. The last column above is a·cos + b·sin: what K1 alone contributes at that hour.
- Step 4 — Two lines of trigonometry. The pair (a, b) is amplitude and phase in another form: H = √(a² + b²) ÷ f = 0.4554 m, and g = atan2(b, a) = 25.1°. Exactly the numbers on K1’s row above. The division by f takes out the 18.6-year nodal correction, so the constant does not depend on which year you happened to observe.