[ 505 ]DIRECTIONAL RECORDING OF SWELL FROM DISTANT STORMSBY W. H. MUNK, G. R. MILLER, F. E. SNODGRASS AND N. F. BARBERInstitute <strong>of</strong> Geophysics and Planetary Physics, University <strong>of</strong> California, La Jolla(Communicated by G. E. R. Deacon, FER.S.-Received 19 April 1962)PAGE1. INTRODUCTION 5062. INSTALLATION 507(a) Site 507(b) Method <strong>of</strong> installation 509(c) The Vibrotron pressure transducer509(d) <strong>Recording</strong> 5103. SAMPLE RECORDS 5124. DIGITAL ANALYSIS 512(a) The detection <strong>of</strong> errors 512(b) 'Aliasing' and hydrodynamicfiltering 512(c) Suppression <strong>of</strong> tides 513(d) The spectral matrix 513(e) The reliability <strong>of</strong> the computedspectra 5145. FIRST INFERENCE 515(a) Sample spectrum 515(b) Power spectra 515(c) Direction finding 515(d) Beam width 517(e) Dispersion 518(f) Effect <strong>of</strong> finite depth on ridgelines 5196. GENERAL THEORY OF ARRAYS 520(a) Elementary wave trains 520(b) The spectrum 5217. INTERPRETATION OF THE CROSS-SPECTRALMATRIX 523(a) Conventional treatment 523(b) A best-fitting single wave train 524(c) A best-fitting pair <strong>of</strong> wave trains 526(d) Experimental test <strong>of</strong> the dispersionlaw 527[Plate 7]CONTENTSPAGE(e) Analysis by iteration 528(f) Fourier-Bessel analysis 531(g) Other methods 5338. RESULTS 5409. REFRACTION 54310. REFLEXION 54911. IDENTIFICATION OF INDIVIDUAL EVENTS 551(a) The hurricane <strong>of</strong> 23 October 556(b) Antipodal swell 557(c) The quickly moving storm <strong>of</strong> 12-4May 558(d) The end <strong>of</strong> the storm season 561(e) Narrow beams and sharp ridges 561(f) A discrepancy in wave direction? 56212. AN EXPERIMENT IN HETERODYNING 562(a) Removal <strong>of</strong> frequency shift 562(b) Heterodyne theory 563(c) The results 56413. GENERATION AND DECAY 565(a) Ridge cuts 565(b) Geometric spreading 566(c) Attenuation 56714. MEAN MONTHLY SPECTRA 56715. SCATTER AND ABSORPTION 571(a) The evidence 571(b) Scattering by waves 572(c) Scattering by turbulence 575(d) Island scattering 576(e) Absorption 58016. THE SOUTHERN SWELL 581REFERENCES 583VOL. 255 A. io62. (Price $x 5s.) 62 [Published i8 April I963
506 W. H. MUNK, G. R. MILLER, F. E. SNODGRASS AND N. F. BARBERWe have measured the distribution <strong>of</strong> wave energy with frequency and direction for several months,and attempted to interpret the resulting field E(f, 0, t) in terms <strong>of</strong> pertinent geophysical processes. Thefluctuating pressure on the sea bottom was measured with a triangular array <strong>of</strong> sensitive transducerslocated 2 miles <strong>of</strong>f shore <strong>from</strong> San Clemente Island, California, at a depth <strong>of</strong> 100 m. The transduceroutput was brought ashore by undersea cable and telemetered to the laboratory at La Jolla,California, where it was digitally recorded. The analysis by means <strong>of</strong> high-speed computers consists<strong>of</strong> two distinct steps: (i) computing the spectral matrix Sijk = Cijk+iQijk between the records xi(t)and xj(t) (i,j = 1, 2, 3) at frequency k (k- 1, 2, ..., 1]00); Cijk, Qijk are the co- and quadraturespectra, and -Siik CLiik the ordinary power spectra <strong>of</strong>' the three records; (ii) interpreting Sijk interms <strong>of</strong> the energy spectruim in frequency and direction (or equivalently in horizontal wave-numberspace). 'T'he problem is discussed in some detail. The method adapted is to find the energy anddirection <strong>of</strong> a point source which can best account (in the least-square sense) for Sijk; to computethe matrix S'k associated with this optimum point source; then to find the optimum point source <strong>of</strong>the 'residual matrix' Sijk.- YS'J (we chose y = 0-1), etc, for 100 iterations, or until the residualmatrix vanishes. The calculation is made independently for each frequency. The iterativepoint-source method is satisfactory if the incoming radiation is associated with one or two strongsources.From daily measurements over several months we obtain E( f, 0, t). The non-directional spectra,E( f, t) =fE(f, 0, t) dO are contoured, and show pronounced slanting ridges associated with thedispersive arrivals <strong>from</strong> individual storms. Their slope and intercept give the distance and time <strong>of</strong> thesource, respectively; the energy density along the ridge is associated with the wind speed, the width<strong>of</strong> the ridge with the storm duration, and the variation <strong>of</strong> width with the rate at which the stormapproaches the station. The distribution <strong>of</strong> energy with direction along these ridge lines give thedirection <strong>of</strong> the source (after correction for local refraction), some upper limit to source apertureand an indication <strong>of</strong> its motion normal to the line-<strong>of</strong>-site. With distance, direction, time and fiveother source characteristics thus derived <strong>from</strong> the wave records, the sources are fairly well specified.Comparison with weather maps leads to fair agreement in most instances. FromJune to Septembermost sources lie in the New Zealand-Australia-Antarctic section or in the Ross Sea. In threeinstances the swell was generated in the Indian Ocean near the antipole and entered the Pacificalong the great circle route between New Zealand and Antarctica. By October the southern winterhas come to an end and the northern hemisphere takes over as the chief region <strong>of</strong> swell generation.Antarctic pack ice may be a factor, by narrowing the 'window' between Antarctica and NewZealand or by impairing storm fetches in the Ross Sea. The northward travelling swell is furtherimpaired by densc island groups in the South Pacific, and an attempt is made to estimate theresulting scattering. The recorded spectrum at frequencies below 0 05 c/s is in rough accord withwhat would be expected <strong>from</strong> typical storms after due allowance for geometric spreading; at higherfrequencies the spectrum fails to rise in the expected manner. A resonant interaction between theswell and the trade wind sea may be responsible for back-scattering frequencies above 0 05 c/s. Theobserved widths <strong>of</strong> the beam (< 0-1 radian) ancd <strong>of</strong> the spectral peak (tAf/f < 0.1) put an upper limiton forward scattering over the transmission path <strong>of</strong> 104 wavelengths, and thus place some restrictionon the energy <strong>of</strong> turbulent eddies with dimensions <strong>of</strong> the order <strong>of</strong> 1 km in the shallow layers <strong>of</strong> thesea. Finally, a comparison <strong>of</strong> mean monthly wave spectra with those <strong>from</strong> individual storm permitsus to estimate the absorption time <strong>of</strong> the North Pacific basin. The resultant value <strong>of</strong> j week is notinconsistent with the assumption that most <strong>of</strong> the swell energy is absorbed along the boundaries.The inferred coastal absorption is roughly in accord with what we found at San Clemente Island(beach slope 1:30): the coast line changes <strong>from</strong> a predominantly reflecting to a predominantlyabsorbing boundary as the frequency increases <strong>from</strong> 0 03 to 0 05 c/s.1. INTRODUCTIONAlmost fifteen years ago in the pages <strong>of</strong> this Journal, one <strong>of</strong> us presented power spectra <strong>of</strong>ocean waves and swell <strong>of</strong>f Pendeen and Perranporth in north Cornwall (Barber & Ursell1948). The outstanding feature <strong>of</strong> these spectra is the successive shift <strong>of</strong> peaks towardhigher frequencies. This is the expected behaviour <strong>of</strong> dispersive wave trains <strong>from</strong> rather
- Page 1: Directional Recording of Swell from
- Page 5 and 6: 508 W. H. MUNK, G. R. MILLER, F. E.
- Page 7 and 8: 510 W. H. MUNK, G. R. MILLER, F. E.
- Page 9 and 10: 512 W. H. MUNK, G. R. MILLER, F. E.
- Page 11 and 12: 514 W. H. MUNK, G. R. MILLER, F. E.
- Page 13 and 14: 516 W. H. MUNK, G. R. MILLER, F. E.
- Page 15 and 16: 518 W. H. MUNK, G. R. MILLER, F. E.
- Page 17 and 18: 520 W. H. MUNK, G. R. MILLER, F. E.
- Page 19 and 20: 522 W. H. MUNK, G. R. MILLER, F. E.
- Page 21 and 22: 524 W. H. MUNK, G. R. MILLER, F. E.
- Page 23 and 24: 526 W. H. MUNK, G. R. MILLER, F. E.
- Page 25 and 26: 528 W. H-. MUNK, G. R. MILLER, F. E
- Page 27 and 28: 530 W. H. MUNK, G. R. MILLER, F. E.
- Page 29 and 30: 532 W. H. MUNK, G. R. MILLER, F. E.
- Page 31 and 32: 534 W. H. MUNK, G. R. MILLER, F. E.
- Page 33 and 34: 00207'-0.C-0-01~~~~~~~~~~~~~~~~~~~~
- Page 35 and 36: 538714Cld0,9~~~~~~~U0CldCIO.0a.)a.)
- Page 37 and 38: 540 W. H. MUNK, G. R. MILLER, F. E.
- Page 39 and 40: 542 W. H. MUNK, G. R. MILLER, F. E.
- Page 41 and 42: 544 W. H. MUNK, G. R. MILLER, F. E.
- Page 43 and 44: 546 W. H. MUNK, G. R. MILLER, F. E.
- Page 45 and 46: 548 W. H. MUNK, G. R. MILLER, F. E.
- Page 47 and 48: 550 W. H. MUNK, G. R. MILLER, F. E.
- Page 49 and 50: 140? 1300 1200 1100 1000 900 80016t
- Page 51 and 52: ;554 W. H. MUNK, G. R. MILLER, F. E
- Page 53 and 54:
556 W. H. MUNK, G. R. MILLER, F. E.
- Page 55 and 56:
558 W. H. MUNK, G. R. MILLER, F. E.
- Page 57 and 58:
560 W. H. MUNK, G. R. MILLER, F. E.
- Page 59 and 60:
562 W. H. MUNK, G. R. MILLER, F. E.
- Page 61 and 62:
564 W. H. MUNK, G. R. MILLER, F. E.
- Page 63 and 64:
566 W. H. MUNK, G. R. MILLER, F. E.
- Page 65 and 66:
568 W. H. MUNK, G. R. MILLER, F. E.
- Page 67 and 68:
570 W. H. MUNK, G. R. MILLER, F. E.
- Page 69 and 70:
572 W. H. MUNK, G. R. MILLER, F. E.
- Page 71 and 72:
574 W. H. MUNK, G. R. MILLER, F. E.
- Page 73 and 74:
576 W. H. MUNK, G. R. MILLER, F. E.
- Page 75:
578 W. H. MUNK, G. R. MILLER, F. E.
- Page 78 and 79:
Munk et al. Phil. Trans A, volume 2
- Page 80 and 81:
582 W. H. MUNK, G. R. MILLER, F. E.
- Page 82:
584 W. H. MUNK, G. R. MILLER, F. E.