Columbia ainlbersstti?
PALISADES, NEW YORK
Lamont Geological Observatory
Technical Report on Seismology
No. 19
A Study and Evaluation of The Tripartite Seismic Method of Locating Hurricanes
Contract N6 onr 27133
LAMONT GEOLOGICAL OBSERVATORY
(Columbia University) Palisades, New York
A Study and Evaluation of
The Tripartite Seismic Method of Locating Hurricanes
Technical Report No. 19 by
William L. Donn and Maurice Blaik
The research reported in this document was supported by Contract N6-onr-27133 with the Office of Naval Research of the United States Navy Department.
February 1952
ABSTRACT
Tripartite records from the U. 5. Navy stations at Bermuda, Cherry Point and Miami were studied in detail for the 1950 hurricane season. Following a discussion of the theoretical error expected in this study, the results of single and average azimuth computation are given. Since the error between computed and observed storm azimuths exceeds the theoretical by a considerable amount, a study of the causes of the errors was undertaken. The difficulties are considered to result from instrumentations and procedure, lack of wave coherence at the three elements of the tripartite net, and refraction or multiple wave paths.
A large range in velocities was observed with indications that the lower values are the more reliable. Selection on a veloc¬ ity basis gives somewhat better success than averaging all readings over the chosen interval of time. Suggestions for improvement of the instrumentation program are given which may help reduce exist¬ ing errors and give storm location with greater certainty. Opera¬ tional success may then be obtained more frequently.
INTRODUCTION
The use of more than one seismic station in determining storm position or direction of wave approach has been described by a number of investigators, chiefly Hecker (1), Shaw (2), Krug (3), Trommsdorf (4), Ramirez (5), Gilmore (6, 7, £), lynch (9) and Kammer and Dinger (10). Owing to the diversity of results obtained in the application of this procedure to the location and tracking of hurricanes, an intensive study and evaluation of the tripartite method was undertaken. Seismograms from U.S. Navy Hurricane Tracking stations were studied for most of the stations for the five hurricanes in the Western North Atlantic Ocean for 1950. The data given here are for the stations at Bermuda, Cherry Point and Miami ("B", "CP" and "M" respectively in Figures 1 and 2). The instruments are reported as being "always in phase", and having a galvanometer-seismometer system period of 7 seconds (with a possi¬ ble error of 0.1 to 0.2 seconds). The magnification is reported as being 5000 for all instruments. All instruments were oriented N-S, with the free end of the pendulum to the north. These stations were selected owing to their disposition in a great triangle, and the fact that one is an island station. The azimuths computed from seismic data are compared with storm azimuths determined from marine weather charts. The tracks of the hurricanes used in this study are shown in Figures 1 and 2.
- 1
STORM AZIMUTH DETERMINED BY AVERAGE SEISMIC AZIMUTH COMPUTATION
Average Azimuth Computations
The average azimuth of wave approach at any time was determined by selecting the sixteen most regular waves during a six minute interval and finding the arithmetic mean for the six¬ teen individual directions computed for these waves. The aver¬ age deviation was then computed for each set of sixteen waves. Individual directions were obtained by use of the formula develop¬ ed by Gilmore (6) after Krug (3) and Ramirez (5). Table I compares computed average azimuth of wave approach with hurricane azimuth for the stations referred to above. The observed storm azimuth given is for a line to the storm center, but the tables also give the angle subtended at each of the stations by the effective wind area of the storm. This was usually symmetrical about the line to the center. Two observed azimuths are given where obvious ambigu¬ ity existed owing to the presence of simultaneous hurricanes. Discussion of Procedure
The traces from the three elements of the tripartite station were recorded on a single drum when azimuths were to be determined, and drum speeds we re increased from the normal l/2 mm per sec to 5 mm per sec at Bermuda and Miami, and to 2 mm per sec at Cherry Point. Then the interval from an arbitrary time to the crest of the same wave on each trace was measured. From these measurements are determined arrival order (which is often obvious from visual examination) and the arrival time differences between the waves of first and last arrival(A t]_), and the wave of first
|
Date 1950 |
G.G.T . |
TABLE I. COMPARISON OF COMPUTED AVERAGE AZIMUTHS WITH OBSERVED STORM AZIMUTHS BERMUDA Angle sub¬ tended at Observed Station by Computed Storm Storm Date Computed Az im- A. D. Azimuth (degrees) 1950 G.C.T. Azim-A.D. |
Observed Storm Azimuth |
2 Angle sub¬ tended at Station by Storm (degrees) |
|||||||
|
Aug. 15 |
0700 |
030 |
3 |
021/194 |
N* |
Aug. 30 |
0600 |
060 |
4 |
095 |
45 |
|
1200 |
020 |
6 |
021/196 |
25/38 |
1200 |
020 |
8 |
092 |
40 |
||
|
2100 |
040 |
4 |
016/202 |
20/32 |
2200 |
050 |
2 |
087 |
54 |
||
|
Aug. 16 |
0300 |
050 |
5 |
018/205 |
14/34 |
Sept.l |
1500 |
060 |
6 |
089 |
80 |
|
0600 |
060 |
5 |
018/207 |
10/38 |
2400 |
080 |
6 |
082 |
70 |
||
|
1200 |
050 |
3 |
019/209 |
24/38 |
Sept . 2 |
0600 |
060 |
9 |
084 |
80 |
|
|
1S00 |
210 |
17 |
213 |
38 |
1200 |
050 |
5 |
O63 |
75 |
||
|
2400 |
220 |
11 |
217 |
40 |
1800 |
050 |
8 |
055 |
75 |
||
|
Aug. 17 |
0600 |
240 |
14 |
220 |
45 |
2400 |
120 |
11 |
052 |
75 |
|
|
1200 |
210 |
2 |
226 |
36 |
Sept .3 |
1200 |
040 |
5 |
049 |
80 |
|
|
1S00 |
200 |
6 |
230 |
24 |
1800 |
050 |
5 |
046 |
60 |
||
|
Aug. IS |
0300 |
210 |
4 |
237 |
28 |
2400 |
050 |
6 |
040 |
50 |
|
|
0900 |
230 |
11 ' |
242 |
31 |
Sept. 4 |
0600 |
020 |
11 |
034 |
40 |
|
|
1500 |
200 |
6 |
247 |
32 |
1200 |
030 |
5 |
038 |
70 |
||
|
1S00 |
190 |
4 |
250 |
33 |
1800 04(/L80 10/11 039/191 |
40/35 |
|||||
|
2400 |
220 |
S |
253 |
37 |
2400 190/010 |
7/16 |
194/041 |
50/30 |
|||
|
Aug. 19 |
0900 |
230 |
6 |
259 |
30 |
Sept . 5 |
1200 |
330 |
11 |
045/199 |
30 |
|
1500 |
230 |
4 |
265 |
32 |
2400 |
010 |
6 |
054/208 |
13 |
||
|
1S00 |
220 |
6 |
267 |
37 |
Sept . 6 |
1200 |
340 |
10 |
213 |
60 |
|
|
2400 |
220 |
6 |
275 |
40 |
2400 |
250 |
10 |
220 |
60 |
||
|
Aug. 20 |
0600 |
200 |
5 |
285 |
45 |
Sept. 7 |
1200 |
340 |
5 |
234 |
90 |
|
1200 |
240 |
5 |
300 |
45 |
2400 |
340 |
3 |
234 |
120 |
||
|
1S00 |
240 |
5 |
318 |
45 |
Sept. 8 |
0700 |
330 |
8 |
236 |
115 |
|
|
2400 |
220 |
S |
333 |
65 |
1200 |
290 |
5 |
239 |
120 |
||
|
Aug. 21 |
0600 |
230 |
6 |
347 |
45 |
1800 |
280 |
8 |
239 |
135 |
|
|
1200 |
240 |
7 |
167/003 |
12/47 |
Sept. 9 |
1200 |
290 |
6 |
245 |
130 |
|
|
1S00 |
240 |
9 |
168/008 |
13/32 |
2400 |
290 |
8 |
250 |
100 |
||
|
2400 |
200 |
5 |
I69/OII |
12/20 |
Oct .13 |
2300 |
320 |
12 |
134 |
30 |
|
|
Aug. 22 |
0600 |
190 |
2 |
170/016 |
13/20 |
Oct .14 |
2300 |
300 |
12 |
108 |
60 |
|
1200 |
200 |
2 |
171/024 |
N |
Oct. 15 |
1200 |
240 |
6 |
081 |
70 |
|
|
2200 |
320 |
17 |
175/029 |
N |
2300 |
230 |
18 |
068 |
50 |
||
|
Aug. 29 |
2400 |
~0§0~ |
S |
100 |
45 |
||||||
|
CHERRY POINT |
|||||||||||
|
Aug. 16 |
1200 |
090 |
2 |
146 |
26 |
Aug. 19 |
0300 |
100 |
4 |
165 |
50 |
|
1S00 |
100 |
2 |
149 |
28 |
0600 |
100 |
4 |
165 |
50 |
||
|
2400 |
100 |
1 |
152 |
34 |
0900 |
100 |
2 |
163 |
58 |
||
|
Aug. 17 |
0600 |
090 |
1 |
154 |
35 |
1500 |
100 |
3 |
158 |
75 |
|
|
1200 |
100 |
2 |
157 |
28 |
1800 |
110 |
3 |
156 |
85 |
||
|
1S00 |
090 |
3 |
159 |
28 |
2400 |
100 |
5 |
143 |
140 |
||
|
2400 |
090 |
2 |
163 |
32 |
Aug. 20 |
0600 |
090 |
3 |
108 |
120 |
|
|
Aug. IS |
0600 |
100 |
3 |
166 |
35 |
1200 |
090 |
5 |
072 |
95 |
|
|
0900 |
100 |
3 |
167 |
40 |
1800 |
090 |
3 |
060 |
80 |
||
|
1500 |
100 |
2 |
168 |
57 |
2400 |
090 |
2 |
055 |
54 |
2100 090 3 167 60
( cont 1 d )
-3
TABLE I.-(cont'd)
Cherry Point
Angle sub- Angle sub¬ tended at tended at
Observed Station by Observed Station by
Date Computed Storm Storm Date Computed Storm Storm
|
1950 |
G.C.T. |
Azim- |
-A.D |
.Azimuth |
(degrees ) |
1950 |
G.C.T. |
Azim- |
-A.D. |
Azimuth |
(degrees) |
|
Aug. 21 |
0600 |
080 |
2 |
052 |
41 |
Sept .9 |
0600 |
100 |
3 |
116 |
35 |
|
1200 |
090 |
3 |
050/142 |
37/9 |
1200 |
100 |
5 |
118 |
50 |
||
|
1800 |
080 |
2 |
048/142 |
32/9 |
Sept .10 |
2400 |
100 |
2 |
108 |
80 |
|
|
2400 |
080 |
3 |
045/143 |
20/9 |
Sept. 11 |
0600 |
100 |
4 |
090 |
90 |
|
|
Aug. 31 |
1200 |
110 |
7 |
“096 |
23 |
0900 |
080 |
2 |
079 |
90 |
|
|
2400 |
310 |
7 |
097 |
N |
1200 |
080 |
3 |
071 |
90 |
||
|
Sept.l |
1200 |
100 |
3 |
100 |
30 |
1800 |
080 |
2 |
059 |
75 |
|
|
Sept. 4 |
0020 |
~Wo~ |
7 |
086/140 |
30/15 |
2400 |
100 |
3 |
054 |
55 |
|
|
1200 |
130 |
3 |
076/141 |
35/25 |
Sept. 12 |
0600 |
050 |
1 |
053 |
50 |
|
|
1800 |
no |
4 |
074/142 |
18/23 |
1200 |
100 |
4 |
052 |
45 |
||
|
2400 |
090 |
3 |
072/143 |
25/30 |
1800 |
120 |
2 |
054 |
40 |
||
|
Sept . 5 |
0600 |
090 |
5 |
143 |
30 |
2400 |
070 |
3 |
056 |
35 |
|
|
1200 |
no |
6 |
143 |
26 |
Sept. 29 |
1800 |
090 |
3 |
ll6 |
N |
|
|
1800 |
080 |
2 |
143 |
30 |
Oct .1 |
1200 |
100 |
4 |
116 |
25 |
|
|
2400 |
100 |
3 |
140 |
30 |
1800 |
no |
5 |
118 |
30 |
||
|
Sept • 6 |
0600 |
100 |
2 |
138 |
30 |
2400 |
100 |
3 |
117 |
35 |
|
|
1200 |
100 |
3 |
135 |
40 |
Oct. 2 |
0600 |
098 |
2 |
114 |
40 |
|
|
1800 |
110 |
4 |
132 |
37 |
1200 |
no |
6 |
112 |
35 |
||
|
2400 |
090 |
5 |
128 |
30 |
1800 |
100 |
2 |
113 |
40 |
||
|
Sept. 7 |
0600 |
100 |
5 |
124 |
30 |
2400 |
100 |
4 |
113 |
40 |
|
|
1200 |
090 |
4 |
122 |
40 |
Oct. 3 |
1200 |
100 |
3 |
112 |
45 |
|
|
1800 |
080 |
2 |
121 |
35 |
1800 |
090 |
3 |
105 |
45 |
||
|
2400 |
090 |
5 |
120 |
50 |
2400 |
100 |
3 |
098 |
50 |
||
|
Sept. 8 |
0600 |
090 |
5 |
118 |
35 |
Oct. 4 |
0600 |
100 |
3 |
087 |
55 |
|
1200 |
080 |
3 |
117 |
40 |
1200 |
090 |
3 |
075 |
35 |
||
|
1800 |
080 |
2 |
116 |
35 |
1800 |
100 |
4 |
067 |
35 |
||
|
2400 |
100 |
5 |
115 |
40 |
2400 |
100 |
2 |
063 |
30 |
||
|
MIAMI |
|||||||||||
|
Aug. 18 |
1200 |
090 |
4 |
064 |
60 |
Aug. 20 |
0500 |
070 |
3 |
032 |
50 |
|
1800 |
100 |
2 |
055 |
64 |
1200 |
070 |
4 |
032 |
36 |
||
|
Aug. 19 |
0300 |
080 |
4 |
048 |
54 |
1800 |
060 |
8 |
034 |
26 |
|
|
1400 |
080 |
2 |
039 |
49 |
Aug. 21 |
0500 |
070 |
8 |
036 |
21 |
|
|
2400 |
070 |
2 |
033 |
34 |
* - N = insufficient data
*#**##****■» f
-4
and second arrival (At2). The maximum time difference possible occurs for a wave traveling parallel to a leg of the tripartite triangle. If a wave should arrive from a direction 180 degrees away, this maximum time difference would be the same, but the arrival order would be reversed.
Time measurements were made to one-tenth mm (0.02 sec for Bermuda and Miami, and 0.05 sec for Cherry Point) and were only converted to seconds for purposes of velocity computations.
The physical nature of the records did not warrant the use of more refined measurements. This was verified by trial measurements.
Azimuths were computed from the formula:
sin B
tan A a -
Rt Rs - Cos B
where: A = direction angle between the wave front and the leg of the tripartite station connecting the elements of first and second arrivals, and thus refers to different legs for different arrival orders;
B = the vertex angle of the triangle at the element of first arrival; . .
Rt « ^ or the ratio of the time difference between
A
the first and last arrivals to that between first and second arrivals.
Rs = the ratio of the length of the leg connecting the elements of first and last arrivals to the length of the leg connecting the elements of first and second arrival.
It is important to realize the theoretical accuracy that can be obtained for the tripartite stations used. In making these determinations it is assumed that.no significant velocity differences exist for the microseism periods observed. The Bermuda net will be taken as an example. This triangle is essentially equilateral with legs of 1800 feet. The angle A as defined, can only have values from zero to sixty degrees for any arrival order. With the further assump¬ tion that a surface wave velocity of 10,000 feet per sec exists at
-5
Bermuda (justification for this will be given later), then the largest time differences (Atq) for a given arrival order will be 0.9 mm (0.18 sec) in view of the precision obtainable here. Again in view of the precision, A t2 can vary from 0.0 mm (0.00 sec) to 0.9 mm (0.18 sec) in steps of 0.1 ram (0.02 sec). Hence the angle A for any arrival order is determined by one of the ten possible ratios ( A t^/At2), and gives one of ten possible azimuth sectors whose size is 60/10 degrees. However, considera¬ tion of the significance to be attached to the ratios based on trial measurements shows that two consecutive ratios may not be truly distinct. Consequently the best theoretical accuracy would be double the above sector (12 degrees) with half the number of possible sectors (30).
Several improvements in instrumentation are immediately suggested in view of the above discussion: a) an increase in the size of the network to a limit imposed by the need for recognition of similar waves, and by practical considerations; b) an increase in drum speed always accompanied by an increase in magnification to maintain sharp wave crests; c) an improvement in the quality of the records to permit greater precision of measurement; d) an increase of the number of instruments used in the net to better define the wave motion and give additional data for computations. The combined advantage of the first three suggestions would be to permit greater precision in measurements. This would reduce the size of the theoretical sector of error for each station, and increase the possible number of azimuths obtainable.
-6
Discussion of Results
Table I shows that agreement between the computed and observed azimuths occurred only 13 times for the 148 sets of computations. Agreement is here considered to be cases where the computed azimuth sector (determined by the A.D.) includes the azimuth of the storm center. Poor agreement still occurs if cases are considered where any part of the computed sector overlaps the sector to the effective wind area. Table II summarizes the re¬ sults for both cases.
TABLE II.
Bermuda Cherry Point Miami
Total Computations 53 74 9
Success using azimuth of center 8 5 0
Success using sector of effective 24 32 1
wind area
Striking negative correlation between computed and ob¬ served azimuths occurred for the Bermuda station for the hurricane of October 13-17. Table I indicates computed azimuths to be approx¬ imately 180 degrees in error even when the storm made its closest approach, with coincident maximum amplitudes occurring. There was clearly no meteorological ambiguity at the time. Similar results have been reported at other stations during other hurricane seasons. Careful examination of the responses and the physical characteristics of each trace gives no obvious indication that the cause is instru¬ mental.
The above data and discussion is based on average azimuths, following the standard procedure for such determinations. However, if individual wave azimuths are considered, the data can be given as
-7
in Figure 3, showing angle of error against frequency of occurr¬ ence. (The data for the clearly anomalous case for Bermuda - October 13 to 17 - were not used.) Smoothed frequency data are shown since each point includes values spread on both sides of the designated ordinates. The curve for Bermuda is based on 669 computed azimuths, for Cherry Point 749 and Miami 116. The dis¬ tribution of the angle error is far from random, showing definite modes. The deviations of the Cherry Point and Miami curves in¬ dicate a systematic error very possibly a result of refraction. Despite the peak near zero for Bermuda, the shape of the curve indicates poor accuracy was obtained. It appears that the method does give azimuths although not yet accurate enough for operation¬ al purposes.
This leads to a further consideration of the causes of error and possible remedies. The sources of error resulting from procedure account for only part of the total difficulty. Other, and possibly more significant causes are indicated from the study, especially when individual waves are considered. Individual wave azimuths were based on the ratio of the time differences Atq and A tp. The observed time differences for Atq showed consider¬ able variation, which for Bermuda ranged from 0.1 to 1.1 mm (0.02 to 0.22 sec). From the theoretical conditions given in the pre¬ ceding section this quantity should have a very small range. Even allowing reasonable velocity variations this tenfold range seems far too great. Careful examination and measurement of similar
-8
waves on each trace showed that they are rarely coherent parti¬ cularly with respect to period. These period differences are often adequate to account for the anomalous values of At]_.
STUDY OF WAVE COHERENCE
The computation of wave direction is based on the assump¬ tion that the particular wave is coherent at all three elements. Measurements of the periods of supposedly identical waves almost always indicated differences in period of the order of magnitude of the arrival time differences used in the computations. These period differences are sufficient to account for the large sector over which individual wave azimuths varied during a particular set of observations (six minutes). This sector was often as large as 90 degrees even in cases of good agreement between computed aver¬ age azimuth and observed azimuth. This effect, and the usual occurrence of microseisms in beat patterns suggest that the in¬ coherence of the individual waves over the small distance separat¬ ing the elements is a result mainly of waves arriving at the same time from different directions. Similar observations and conclus¬ ions were noted by Kamraer and Dingep (10), Leet (11) and have been developed theoretically by Bungers (12). Velocity considerations, to be given below, further support this conclusion.
The approach of microseisms from different directions at the same time may be due to a combination of factors, namely, simultaneous sources at different parts of the storm area, two or more storms at different azimuths from the station, and refraction
-9
at the continental borders, the latter being admitted for earth¬ quake Rayleigh waves. These results indicate two additional and possibly correctible sources of error. First, the instrumental frequency response at present is so broad that no single source of microseisms can be studied. It has already been shown by Donn (13) that microseism period is apparently a function of water depth conditions in the generating area. Secondly, the separation of the seismographs is not suitable to the order of magnitude of the time measurements necessary. Additional elements in the net would further correct this.
Since lack of coherence appears to be a major source of error in azimuth determinations, azimuths were computed using individual waves and the results for each wave were compared with its coherence at the three elements. The parameters of amplitude and period were used as the measure of coherence, and only regular waves that showed no obvious incoherence were selected. Eight observation times covering six minutes were used and were taken from the previous data. Only 5 of 131 waves had the same period at each element for the precision used, and these show no correla¬ tion with success. Differences in period for the others varied from 0.1 mm to 1.4 mm. None of the waves showed constant -ampli¬ tude. The lack of success for the five coherent waves, assuming this small number to be significant, may be explained by (a) the presence of composite waves formed by waves of the same period but traveling along different paths and (b) refraction. By analogy with earthquake seismology the latter may be assumed to
-10
be often of large magnitude, especially for the short-period waves studied.
Considerable variations in amplitude existed at the three elements, with a systematic but not constant difference among them. For example, Element A at Cherry Point always show¬ ed much higher values than the other two elements, etc. Since the instruments have been described as being nearly identical in response and orientation, the systematic amplitude diversity is possibly the result of differences in anchoring or in local sur- ficial geology. This effect is significant if also accompanied by phase differences which would effect azimuth computations.
I
STUDY OF VELOCITIES
Velocity Data
The method which was used by Kammer and Dinger (10) with some indications of success was applied in this study. The procedure consisted of computing azimuths and velocities for a series of waves, and attempting to reduce the angular spread of results by considering only azimuths determined from waves show¬ ing velocities below 11,000 feet per second. Encouraging but unsatisfactory results were obtained by the use of this proced¬ ure. However, an analysis of the velocity data is given below since it reveals significant information bearing on the problem.
The curves in Figure U show the frequency distribution of velocities on a logarithmic scale for Bermuda, Cherry Point and Miami for most of the hurricanes studied. A total of 719 individual wave velocity determinations are used for Bermuda,
-11
731 for Cherry Point and 116 for Miami. The points shown are
t
plotted at the median values of the populations in non-overlapp¬ ing velocity sequences. Each curve reveals a considerable spread in velocities. Table III summarizes the minimum, modal, mean and maximum velocity values for each station. The maximum values have much less reliability than the others since they depend on the smallest time differences that can be measured (0.1 mm on the records). Hence variations are difficult to distinguish with existing drum speeds. By the same reasoning the reliability is greatest for the minimum values.
TABLE III.
|
Bermuda |
Cherry Point |
Miami |
|
|
Minimum Velocity |
7,500 |
2,100 |
5,100 |
|
Modal Velocity |
17,000 |
3,800 |
11,000 |
|
Mean Velocity |
17,000 |
6,200 |
13,000 |
|
Maximum Velocity |
80,000 |
20,000 |
60,000 |
It is obvious that velocities in general are intermed¬ iate for Miami, and distinctly lowest for Cherry Point. It is further obvious, and considered of significance that each curve shows a decided concentration of velocities even though not at the same values.
To study the significance of the velocity data given, an analysis was made of the success obtained for waves of differ¬ ent velocities. Cumulative frequency curves are given for Bermuda and Cherry Point in Figures 5 and 6, respectively. These show angle of error (to storm center) against cumulative frequency of error for the velocities given by the curves. All velocities lie
-12
in a rather narrow band. However, to consider Figure 5, for Bermuda, first, it is seen that the curves of lower velocity (8,500, 10,000, and 11,000 feet per second) lie above the others, almost overlapping, and indicate greater relative success. It is of note that these values are represented by ’relatively few observations and are at the low velocity end of the curve. It should be noted that "success" here is quite rela¬ tive since forty to sixty percent of the best observations still show errors of 20 to 30 degrees.
The cumulative frequency curves for Cherry Point (Figure 6) are in general similar to and no better than Bermuda. However, in this case the higher velocities show somewhat better success, with best success given by the 7,000 ft ./sec curve. The low and high portions of all of the curves for both Bermuda and Cherry Point are of lesser reliability than the central portions owing to a very irregular and very sparse distribution of velocities at low and high angles of error, respectively.
The Miami data were too few for analysis in this manner.
It is apparent that although the unique frequency curves show concentrations of velocities, these most frequent velocities do not give the best relative success. This suggests that determina¬ tions of azimuths by averaging data from the entire record should give generally poorer results than determinations based on selection. The broader distribution of Bermuda velocities may reflect the great¬ er potential sources referred to earlier in connection with the A.D.
differences at the stations
-13
A further analysis has been made in order to note any possible trends in average velocity with time. The data from Bermuda were used for this purpose. It was discovered that very short-time variations occurred which could not be related to any obvious causes. However a definite trend was apparent from September 3 through September 8. The data for this trend, together with simultaneous average micro seism amplitude and period are shown in Figure 7. Velocities, which showed about one hundred percent variation, appear to have been independent of period. The ampli¬ tude curve shows a distinct maximum at the time of maximum velocity. An amplitude minimum occurs about 1200 on September 6. During the time of maximum velocity and amplitude, the records show that most of the waves appeared to arrive simultaneously at all three elements. In computing velocities only the waves (the minority) showing measurable time differences were used. This condition existed for a day or more and correlates with the presence of two hurricanes about 180 degrees apart, and approximately equidistant from the station. The tracks of these storms are shown in Figure 2. The southern storm was approaching as the northern storm receded. Hence the amplitude high on September 5 is interpreted as marking the time when the combined effect of both storms was at a maximum. During this intens¬ ity increase, time differences between "unique" waves at the Bermuda tripartite elements diminished, until a standing wave effect occurr¬ ed, the time of which corresponds to the time of amplitude maximum on September 5. Velocities, which vary inversely with time differ¬ ences increased to a maximum at this time. Table I shows that at
-14
times during this interval selection on the records was possible giving directions roughly toward both storms. Velocities decreas¬ ed as time differences increased from the decreasing effect of the northern storm. A second and larger amplitude maximum occurred with the close approach of the southern storm, however velocities continued to fall to a low level. The above case suggests again that velocity determinations, using tripartite stations, are often velocities of composite waves, and in such cases storm azimuth computations must be erroneous. Further, in this case, the cause of the composite waves can be ascribed to the presence of wave paths from two distinct source areas.
Discussion
The large variations of velocity for waves of the same character, plus the occurrence of beats, suggest that the record¬ ed microseisms are frequently caused by the superposition of two or more pure waves approaching from different directions. With this assumption, analysis of these interference beats revealed that frequently, only two wave trains differing by approximately ten percent in period, or a continuous disturbance over this range, could have caused the microseism patterns for intervals of at least thirty seconds. This together with observed velocity variations, leads to an interpretation of the apparently anomalous velocities and directions computed from a "unique" wave.
The assumption is made that the velocities of the compon¬ ent pure waves are practically independent of wave period for the
periods and tripartite distances involved. Then for the case of two such wave trains arriving at a station with the phase velocity v
-15
and the angle C between their paths, the direction of the compos¬ ite wave will be that of the bisector of C, and its phase velocity will be:
V ~ Cos c/2 ^
If significant phase differences exist appropriate modifications can be made to this simple formula.
This interpretation shows that the apparent velocity measured by a tripartite system will vary from a minimum value (the true phase velocity of the pure wave) to an infinite velocity (giving a standing wave) as C varies from 0 to 180 degrees. Thus the modal velocities shown in Figure 4 would be a function of the most frequent angular separation between the paths of waves arriving simultaneously at a tripartite station. According to the formula (1) this separation must be of the order of 90 degrees or more to account for these modal values. The true surface wave velocities should thus be close to the lowest velocities computed, or 2100,
5100 and 7500 feet per sec for Cherry Point, Miami and Bermuda respectively. These are infrequent values which is expected from a consideration of the factors of origin and propagation given.
No satisfactory results in azimuth determinations were obtained on the basis of the lowest velocity. However these values were very infrequent and were probably associated with serious refraction.
-16
SUMMARY AMD CONCLUSIONS
1. Theoretically a given sequence of arrivals at the elements of a tripartite network indicates qualitatively a sixty degree sector of possible azimuths. This may be narrowed by con¬ sidering the ratio A t]/ A tp. These time differences depend upon wave velocity, direction of wave approach, and separation between seismographs, with the factor of drum speed further affect¬ ing the precision of the results. In most cases only one signifi¬ cant figure for time differences could be carried from the measure¬ ments, which permitted a theoretical reduction of the qualitatively determined sixty-degree sector to no better than about ten degrees.
2. Empirical studies do not support the theoretical con¬ clusions as to accuracy since a much greater error occurred for azimuths computed on both an average and a selective basis.
Further, observed maximum arrival time differences ( A t^) show a far greater range than is expected for the tripartite station and wave velocities used.
3* Directions to storm centers based on average azimuth computations gave accuracy too poor for operational purposes. How¬ ever rough azimuths were obtainable with angles of error from 20 to 40 degrees.
4. Several factors appear to contribute to the lack of success in locating storm areas. These are summarized as (a) errors resulting from procedure of measuring and computing, (b) errors resulting from the incoherence of waves recorded at the elements of
-17
the tripartite station, and (c) errors resulting from refraction and multiple wave paths.
5. To explain the discrepancies noted, attention was directed to the study of individual waves recorded at each element of the tripartite nets. Definite, often pronounced differences in period and amplitude of "unique" waves at each of the seismograph elements indicate that the waves and groups measured are incoherent, with the period differences being of the same order as arrival time differences. This precludes any accuracy in computational results depending on such time differences. In some cases the presence of apparently identical waves arriving simultaneously suggests the existence of standing waves from opposite sources.
6. The above led to a study of individual and average wave velocities for the three stations used. In general wave velocities were lowest for Cherry Point and highest for Bermuda with intermediate values for Miami. These might be a result of local geology, and are in agreement with the known geologic relations among the stations. Anomalously low and high velocities are observ¬ ed at all stations, although a definite concentration is noted for each. This study suggests that present tripartite records show the progress of a composite wave form across the net rather than a unique microseism wave traveling a unique path. The study of average velocities for a particular case furthers this view and permits the distinction between two source areas, for the case given.
7. The cause of the ambiguity in recognizing pure waves is considered the presence of multiple paths. These in turn prob¬ ably originate from a combination of refraction at coastal zones.
-IS
two or more source areas, and broad source areas. Based on the assumptions made, both empirical and theoretical results suggest that the lowest velocity values observed at a station most nearly approach the velocities of the component waves forming the record¬ ed microseisms, which is similar but more specific than the find¬ ings of Kammer and Dinger (10). Reliable azimuths should then be computed from the pure waves and not the composite microseisms.
Such waves are difficult to distinguish with present instrumenta¬ tion. Even if unique wave paths can be recognized, they may still differ from true storm azimuths owing to refraction.
8. This study has suggested several improvements in instrumentation, listed below, which may increase the operational value of seismic storm location and reveal further significant data on the basic nature of microseisms.
(a) Use more than three instruments in a net for operational purposes and as many as possible in a research net for further study.
(b) All of the instruments should be vertical components, with at least one instrument having two matched horizontals associated.
(c) The instruments should be sharply tuned to minimize inter¬ ference of waves of different period.
(d) Instruments should be spaced further apart than at present for operational purposes and at variable distances for further re¬ search. An array of numerous instruments along intersecting lines
at right angles would give more information on wave propagation,
«
and would also provide several networks of different spacing for study of azimuths.
-19
(e) Amplitudes should be increased in proportion to increas¬ ing drum speeds.
(f ) An improvement in record quality should be made to permit greater measurement precision, for example frequent simultaneous brief interruptions of all light beams and finer line reproduction.
9. It is considered at present that the operational value of tripartite stations in locating and tracking storms is small, and that almost as much can be determined from a qualitative appraisal of the records as from time-consuming measurements and computations. It is further believed that attention should be concentrated in research, both experimental and theoretical, on the origin and propagation of microseisms before attempting operational application.
-20
ACKNOWLEDGMENTS
This research was supported by Contract N6-onr-27133 between Columbia University and the Office of Naval Research.
The seismograms used in the study were loaned by M. H. Gilmore from the records of the Navy Microseismic Project in the Hurricane Weather Central, Miami, Florida. The weather data was obtained from marine weather charts supplied by the United States Weather Bureau Office at La Guardia Field, Long Island.
The writers gratefully acknowledge the help and criticism given by Doctors Maurice Swing and Frank Press of the Lamont Geological Observatory in the preparation of the manuscript.
>L Ji. X. -M-
J\ /> /\ /> /\
BIBLIOGRAPHY
(1) Hecker, 0.,nVersuche zur Bestiramung der Fortpflanzungs- geschwindigkeit der Bodenbewegung bei der mikr. Unruhe", Mitt.
Zent. Int. Seism. Ass., Gerlands Beitr, Geophysik V14, 28-33,
1915.
(2) Shaw, J. J. , "Communication de M.J.J. Shaw sur les mouve- ments microseismiques", C.R. Union Geod. et Geoph. Int., Rome,
52-53, 1922.
(3) Krug, H.D., "Ausbreitung der naturlichen Bodenruhe nach Aufzeichnungen mit transportablen Horizontal^eismographen", Z.
F. Geophysik V.13, 328-348, 1937.
(4) Trommsdorf, F., "Untersuchungen uber die naturliche Boden¬ ruhe (Mikroseismik) mit transportablen Dreikoraponentenstationen", Z.F. Geophysik,. V.15, 304-320, 1939.
(5) Ramirez, J.E., "An experimental investigation of the nature and origin of microseisms at St. Louis, Mo., " Bull Seism. Soc. Amer., V.30, 35-84, 139-178, 1940.
(6) Gilmore, M. H., "Microseisms and ocean storms". Bull. Seism. Soc. Amer., V.36, 73-85, 1946.
(7) - "Tracking ocean storms with the seismograph,"
Bull. Am. Met. Soc., V.28, 73-85, 1947.
(8) Gilmore, M.H., and Hubert, W.E., "Microseisms and Pacific Typhoons", Bull. Seism. Soc. Amer., V.38, 195-122, 1948.
(9) Lynch, J. J., "An investigation of two-second microseisms associated with cold fronts and a new method for tracking the cold front center", Trans. Amer. Geoph. Un., V.31, 525-528, 1950.
(10) Kammer, E. W., and Dinger, J.E., "Hurricane swell as a gener¬ ator of microseisms", J. Meteor., V.8, 347-353, 1951.
(11) Leet, L.D., "Discussion of tripartite microseism measurements", Bull. Seism. Soc. Amer., V.39, 249-255, 1949.
(12) Bungers, R., "Die Uberlagerung zweier Wellen verschiedener Herkunftsrichtung", Z. F. Geophysik, V.15, 321-332, 1939.
(13) Donn, W. L., "Frontal microseisms generated in the western North Atlantic Ocean", J. Meteor., V.8, 406-415, 1951.
-H-
-*
_v/
a"
Figure 1. Paths of hurricanes of August 13-24; September 30-0ctober 6; and October 12-18,1950*
Figure 2. Paths of hurricanes of August 27- ^ept ember 6 and August 31-September 14, 1950.
FREQUENCY (PERCENT)
Figure 3.
FREQUENCY (No. of Observations)
25
Figure 4.
CUMULATIVE FREQUENCY (PERCENT)
Figure 5. Graph of cumulative frequency of angles of error associated with velocities at Bermuda. Curves show velocity in feet per second.
CUMULATIVE FREQUENCY (PERCENT)
Figure 6, Graph of cumulative frequency of angles of error associated with velocities at Cherry Point.
1
SEPTEMBER 1950
Figure 7. Velocity, period and amplitude data for Bemruda for September 3-9, 1950.
COLUMBIA LIBRARIES OFFSITE
L -
90645774
'