. . M; £.;’.;:., 31.3. . ->.Il «‘t i . ‘v‘\ H I “,1 A _ , n. - b . . ., . ' L 'f' IA: 3‘. .1- : . . And”, I- ..v I , , .' x':.A‘ .A . ,. I A ' " h“ . _ _ . . v A ., _A..A'"A . A _ 1 4. a , . ‘ I . .A . f . . A I ‘ ‘l' "ZN‘" . u A i - .. I I .. . +1.“ TALE??? ’ A I ‘JR'J E “.531. IuAAr-c : A; A‘ I ‘ . C ~ ILQ‘ 3...“..AA I ‘« ' ‘44:!“1; .. I. a .: -..... 13;. . . ‘l I. A .. -. A. ..-‘ .I ' ' A ' ‘ ‘ ”LEVEL; 1 P .. L . U. _ v- 5* ..A. ‘ r: .z-“Jt' ” ' . I.» _ > . . Ab _, . . ‘ gm.‘ I' A” A , . A 3.- . ‘ 1T. ‘ '1‘ II.AAA;~IW“H'...‘I .. 1 I. ~. . .4. “s A . L A... , . ‘ _ k . . A. .. - . .. I z * ‘ ~‘A"-" A A ; v4. ~ - ' ‘ A ,\ 49.1“ w ‘l M I .A.., r? ,, AA.;'A~..,. ' “'vr '1\I:" A. .m .. ..A m ... “I... “I .u. "ww- ~~L~Jn. .. u Tm L.',. w \s “ J... 45:", .4. H ‘ Ii; n.“ 311$: A. :1 N. ‘ .. Q _ 4' .241...’_‘.‘ .ZZT‘.I_.A‘.I.F.‘.‘. A ' s a» . A . w . ~ ' ‘2‘ .‘"‘ Mu.)- 53:3,:2‘. ~.. a v“. I. \.. ' .1. L3,... -,A. “a; «~_ .,.. .I-u. mg.- .5: .... as. v Q I. .3.'.A~. 4;- ". . . IA , n.>>~‘lv.‘r¢ n.- n'. «11> J' v 9.. 1;..- ... A A. .I m-«a 7...-..5'3. [3:55-- ‘.-A~\ m “V 1;; . q: ~ . a I. I .A n. .5" gum r5 (- 5|Lu323Iu13'. .w A I ' .-:=.--;.-’~.é . 2',“ in; fid‘ . \. . .. IA ., “Ti-"Au-“ II. - - ”ILA“; u“ .. 'u A . . seem.“ \ Nd A. . l ‘ A» 34:..- .u~ll -‘~~- l ¢A~.1A_(“ A .. .AI-v foam“. D\ I. .2- ":' .‘I. . V I... , '5 4-. - s r In“: 7C4.“ .. I», I 'lfihfi IVERSITY LIBRARIES TWWWWWM 3 1293 00891 4230 l This is to certify that the thesis entitled The Eurasian-North American Plate Boundary Through the Area of the Laptev Sea presented by Daniel Richard Olson has been accepted towards fulfillment of the requirements for Masters degree in Geology Major p fessor Date 6/l‘f/90 0-7639 MS U is an Affirmative Action/Equal Opportunity Institution h f mum Mkmzan State Unlve l to my J *— fl PLACE IN RETURN BOX to remove this checkout from your record. TO AVOID FINES return on or before date due. DATE DUE DATE DUE DATE DUE II ‘.——|—__ II II T r—~ 7 MSU Is An Afllrmative Action/Equal Opportunity Institution l __ __ __ Cimm, pita-9.1 THE EURASIAN - NORTH AMERICAN PLATE BOUNDARY THROUGH THE AREA OF THE LAPTEV SEA By Daniel Richard Olson A THESIS Submitted to Michi an State Universi in partial fu llment of the requirements for the degree of MASTER OF SCIENCE Department of Geological Sciences 1990 ABSTRACT The Eurasian -North American Plate Boundary Through the Area of the Laptev Sea By Daniel Richard Olson 24 teleseismic events occurring between 1962 - 1988 in the Laptev Sea region were examined to determine focal depths and focal mechanisms using short-period and long-period synthetic waveform modelling and first-motion analysis. Focal depths are calculated to be 11 _+_ 1 km for earthquakes on a trend along the 130°E meridian, between 77° - 710 North latitude. Focal depths increase by 7 to 14 km to the east and west of 130°E, respectively. This distribution is consistent with a proposed rift margin extending north-south through this region, with focal depths limited by high heat flow along the rift axis. Focal mechanisms are predominately extensional to the north and have greater strike-slip components to the south, with compression being dominant at 71°N. This is consistent with placement of the Eurasian - North American pole of rotation near the Lena River Delta by Cook et a1. (1986). The data are insufficient to further resolve a definitive plate margin. ACKNOWLEDGEMENTS My thanks go to Kazuya Fujita, Michael Velbel, Hugh Bennett and F. William Cambray, and all others within Geological Sciences for making this work possible. Special thanks go to Glenn Kroeger for use of his synthetic seismogram package. This research was supported in part by NSF grant DPP 87-21119, and the Department of Geological Sciences and is gratefully acknowledged. TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES INTRODUCTION THIS STUDY METHODOLOGY INDIVIDUAL SEISMIC EVENTS NOTATION ON FIGURES April 19, 1962 ' November 22, 1984 January 21, 1976 February 22, 1987 August 25, 1964 July 20, 1976 March 21, 1988 April 7, 1969 September 22, 1987 August 5, 1986 April 23, 1977 June 10, 1983 January 1, 1988 ' May 20, 1963 July 21, 1964 iv 12 12 13 16 20 23 23 29 32 35 39 39 42 46 46' . 50. 50 February 1, 1980 DISCUSSION GEOLOGIC INTERPRETATIONS TECTONIC IMPLICATIONS CONCLUSIONS REFERENCES CITED 53 59 . 60 61 62 63 LIST OF TABLES Table 1. Compiled focal information for the present study area. vi Figure 1. Figure 2. Figure 3. Figure 4. Figure 5. Figure 6. Figure 7. Figure 8. Figure 9. Figure 10. Figure 11. Figure 12. Figure 13. Figure 14. Figure 15. Figure 16. Figure 17. Figure 18. Figure 19. Figure 20. LIST OF FIGURES Focal mechanisms , focal depths, and inferred grabens (Fujita et al., 1990) in the Laptev Sea area. Epicentral locations and dates for events studied. Poles of rotation for the Eurasian-North American plates within northeastern Asia Short Period waveform fits for July 21, 1964 event (Fujita et al., 1990a) Solutions for the April 19, 1962 event. Synthetics for the April 19, 1962 event. Solutions for'the November 22, 1984 event. Synthetics for the November 22, 1984 event. Synthetics for the Dziewonski et a1. (1985) mechanism for the November 22, 1984 event. Solutions for the January 21, 1976 event. Synthetics for the January 21, 1976 event. Solutions for the February 22, 1987 event. Synthetics for the February 22, 1987 event. Solutions for the August 25, 1964 event. Synthetics for the August 25, 1964 event. Synthetics for the August 25, 1964 event. Solution for the July 20, 1976 event. Synthetics for the July 20, 1976 event. -. Solution. for the March 21, 1988 event. ' Synthetics for the March 21, 1988 event. vii 11 14 15 17 18 19 21 22 24 25 26 27 28 30 3 1 '33 34 Figure 21. Figure 22. Figure 23. Figure 24. Figure 25. Figure 26. Figure 27. Figure 28. Figure 29. Figure 30. Figure 31. Figure 32. Figure 33. Figure 34. Figure 35. Figure 36. Figure 37. . Figure 38. Solutions for the April 7, 1969 event. Synthetics for the April 7, 1969 event. Synthetics for the April 7, 1969 event. Solution for the September 22, 1987 event. Synthetics for the September 22, 1987 event. Observed records for the Au t5, 1986 event with suspected small P and p phases visible. Solutions for the April 23, 1977 event. Synthetics for the April 23, 1977 event. Solutions for the June 10, 1983 event. Synthetics for the June 10, 1983 event. Synthetics for the June 10, 1983 event. Solutions for the May 20, 1963 event. Synthetics for the May 20, 1963 event. Solutions for the July 21, 1964 event. Synthetics for the July 21, 1964 event. Synthetics for the July 21, 1964 event. Solutions for the February 1, 1980 event. Synthetics for the February 1, 1980 event. viii 36 37 38 40 41 43 45 47 48 49 51 52 54 55 56 57 58 INTRODUCTION The Arctic Mid-Ocean ridge defines one segment of the boundary between the North American and Eurasian plates. This segment terminates at the slope of the Eurasian continent north of the Laptev Sea, northeastern Siberia. South of the ridge terminus, seismicity continues but occurs over a more diffuse zone which does not define a sharp continuation of the plate boundary. The Laptev Sea covers a section of the extensive continental shelf of Eurasia, north of northeastern Siberia (Figure 1). It is bounded to the east by the New Siberian Islands, and to the west by Taimyr Peninsula. Prominent features on the south include the Lena River delta and Bour Khaya Bay, and to the north is the Arctic abyssal plain. Tectonically, this area is of interest because an oceanic ridge enters a continental margin at almost a right angle, rather than at a low angle, as in the Gulf of California. Also a change in stress field on the Eurasian-North American plate boundary occurs in the area, extension to the north, compression to the south (Cook et al., 1986; Cook, 1988; Parfenov et al., 1988), based on focal mechanism solutions for events along the boundary. Several grabens (Figure 1) are inferred to be below the Laptev Sea based on magnetic and gravity surveys (Gaponenko et al., 1968; Genin et al., 1977; Kim, 1986) and refraction lines (Kogan, 1974; Vinogradov, 1984). To the southeast is the Moma rift system of northwest-southeast trending grabens. The grabens are presumed 'tO be of Cenozoic age, are‘ sediment filled, ‘with infered time of opening being earlier (Eocene earliest) for those grabens lying 1 .88 8w >653 2: s 82,: .._m 5 arse 289a 3:25 new 23% 302 .mEmEmsooE .moou .p 239”. woo”. zoOn. WON. ncurugdi Z. L Ions: 332.22 :12: 20°F Iii i Ii ii i/ co..e¢e:4n 'l’l L 200» (um 25...; I O z 4.5 occurring between 1962 and 1988 were choosen for study (Figure 2; Table 1) which occurred within 65° and 79°N. latitude and 118° and 147°E. longitude. METHODOLOGY Waveform modelling of station records is used to determine various focal parameters associated with an earthquake. There are six such parameters that .u_m< anmoozto: £53. «221 canto—5‘ ctoz-ca_ma.5m 9: .2 20222 CD «20¢ .m 230E u o: O 2.2.86.8 u on. u on. m o: D co. m on 14 J / I 1’ 050. 62.3.0» see 5:...ch w on I . . 32.... 8 .838 z o . .8. .Mfinfififi 0 32. 2.4.5 z<_m<§u 4 82 9.2. ssoxom 522. see 8.2:: 3.2 ...o t Eng-cecoN . . 1b. :3 .o .22.: O O 95. .2933... 93925... . u 3.1 . \ com. :2. .38 2! it .323. «um it: Q r “#4.: 9% £2 8 zem i .M one .m mama ..ae no xemEen i .Eeb mama .xooo a gooo i meonenemom numeo aeoom i n wonum mane oa .omm .wm .mvm o.m o.o mv.m~a oo.hh Sam mmm am He: mama mm .NQ .amm .av .w monum man» ma .onm .mv .mma m.m m.m mm.4ma N4.mh 2mo mmm mm mew hmma 9mm .NQ .mmm .am .mma honum many ma .ovm .om .om o.m N.m o.oma mm.mh ZNN mao mm new hmma mm .ND .mv .mh .nm #000 .m .om .mm monum many ma .o .om .mm m.m m.m mm.ova mm.mm Smm mma mm >oz 4mma mm .NQ .mvm .mm .Nva .Eeh .mmm .Nn .44a xoou .vmm .oh .oma Naonum many mm .mmm .mn .mma v.m m.m m>.mma mm.mh Ema mmo oa non mmma emm .NQ .mmm .mm .mam xoou .mmm .om .mwa monum man» mm .mam .om .moa m.m 4.m mm.mma mo.mn £om mha as new omma maam mao .Hum .hmm a ZmHdeumz S E m 020A 2 Ba MZHB mafia ao.unoov .a wanes 9 are of primary interest. These are 1) the strike (azimuth) of the fault plane, 2) the dip of the fault plane, 3) the angle of motion (slip) of the hanging wall block with respect to the foot wall block, 4) the focal depth of the event, 5) the source-time function (STF, rise time, flat time, fall time), which is basically a record of the faulting process (which for simplicity is approximated as a trapizoidal or pyramidal function). STF is also a measurement of event duration and 6) crustal structure in the near source region. Synthetic seismograms are computed for stations in a distance range of 30° - 90° (Helmberger and Burdick, 1979). Stations at less that 30° are generally not used due to the possibility that un-modellable triplication phases may appear in the records, stations at greater than 90° are not used because of effects of the core. First-motions can constrain the focal mechanism of the event, and waveform modelling can be used to estimate all six parameters. Some crustal structures for areas under study have been calculated based on two refraction lines, one shot across Bour Khaya Bay (Kogan, 1974) and another shot to the west of the Lena River delta (Vinogradov, 1984). The interpreted results of gravity and magnetic surveys are presented in Gaponenko et a1. (1968). From this information sediment-filled grabens below the Laptev Sea, are inferred (Vinogradov, 1984). When no estimates of the crustal structure are known, a reasonable simplified crustal model is assumed. For some of the largest (mb 6.0-6.2) well-studied events, a reliable source-time function exists. Jemsek et al., (1986) is notable for producing source-time functions for the August 25, 1964, April 7, 1969, and June 10, 1983, Laptev Sea events. These source-time functions are used as initial estimates for STF for these events. 1 0 To create a model waveform, an initial estimate for all focal parameters for the event is needed. Very often the results of a first motion analysis and the geology/structure of the area (if known) will be the basis for the first estimate in waveform modelling. Waveforms are created through forward modelling using the algorithm of Kroeger (1987). The algorithm uses a spike with an amplitude determined by the radiation pattern of the assumed focal mechanism. This spike is then reflected and transmitted within the near-source crustal layers. The resulting spikes, with various amplitudes and time delays, are convolved with the source-time function for the event, the Futterrnan Q operator for attenuation and the instrument response, to create the synthetic waveform. The remainder of the waveform modelling technique involves changing the six parameters until an acceptable fit for all observed station records for that event is obtained. A good fit requires that all phases in the observed records (discounting high frequency and incoherent noise) are matched with phases of similar amplitude, polarity and time delay in the synthetic waveform. An example from Fujita et al. (1990a) (Figure 4) shows four waveform fits for the July 21, 1964 event. The upper three would be considered "good" fits, the fourth is not good. This process requires that synthetic phase amplitudes approach those of the observed records, and that polarities and time delays are correct. This is a trail and error procedure. Both long-period and short-period station records are generally available for study. Short-period records are useful for small events, long-period records are most useful for studying larger events. For the largest events (111}, _2_ 6.0), short-period records tend to record the motion of the initial break of the event, long-period records record the overall motion. This can allow - two different solutions to be created from data for one event. When long and short 11 Anomaa .am am muaasnv .Eo>o Ema .& 22. .2 9: 8.2962, oozed :96 .e 939“. ma a onm m a no ea a ona a4ao_o sumo: uuuuu emanaeo: "emamo _ mozooumnx. . m.o n.o a.o “new 1;; Acmma. 2H2.NOX Ammaa. x000 am 09 Bah 9mmm Aomma. ZOmAO Basso . 515:} gig EE 5:? éT am mo coma >mv. azw Nom Dom OOOCD 3 1 2 period solutions exist and are in conflict, the long-period record solution is accepted as a better record of overall tectonic stress release. Errors can arise in determining all of the above parameters. Error sources include: noise in the observed record, over-simplified near-source crustal model, or an invalid crustal model, the source-time function was incorrect, or the focal mechanism or focal depth were wrong. These parameters are discussed in Stein and Kroeger (1980). Attenuation and magnification factors have relatively little effect on waveform appearance, and are generally held constant during modelling. Confidence is highest for a solution which produces good waveform fits (similar phase amplitudes, polarities and time delays in the observed and synthetic waveforms) for several stations with wide azimuthal coverage. The solution may not be unique. Always a survey of possible values is made on the chance of detecting another solution. It is possible to quantitatively evaluate the error in waveform matching as is done for long period records, it is not practical at present for short period synthetics. Error margins are discussed event-by-event. INDIVIDUAL SEISMIC EVENTS NOTATION ON FIGURES Stations plotted on a focal mechanism may appear as one of three symbols: open circle, denoting a dilatational first motion; filled circle, denoting a compression; open triangle, denoting no first motion pick was 'made. , Larger symbols with station abbreviations are first mofion picks of this study. Smaller symbols with no station abbreviations indicate first motion picks 1 3 compiled from the 188 or ISC bulletins. Source-time functions in this study are comprised of the rise time, flat time, and fall time of the event. The units of the source-time function are seconds. Crustal structures are represented by four columns of numbers. Each column is labelled below with: column 1: Alpha, indicating P wave velocity (km/sec), column 2: Beta, indicating S wave velocity (km/sec), column 3: Rho, indicating density (g/cm3), column 4: T, the thickness of that layer in kilometers. Focal mechanisms are given as angles of azimuth, dip and slip, respectively. Focal depths are in kilometers. In presenting comparisions of observed and synthetic station records, these records are placed .in pairs, with the observed station record above and the synthetic record below. sz and lpz refer to short period vertical and long period vertical records respectively. April 19, 1962 The April 19, 1962 event (mb 6.2) is located near the northern terminus of the Cherskii mountains, in the south-east section of the study area (Figure 1). First-motions for this event are all compressions. A previous solution by Cook, (1988) (Figure 5) based on 188 reported first-motions finds a poorly constrained mechanism which is purely thrust on a fault striking nearly north-south (strike 170 dip 40 slip 90). First-motions were read and synthetic waveforms calculated for six long-period records - (FigUre 6) and good matches were obtained for the mechanism (strike 80 dip 39 slip 58). This mechanism is rotated by roughly 90 14 .226 $2 .2 .22 6.: .2 26.5.8 .m 2:9“. cm n=._w 0* n=0 ON... mx_m._.m mmm— .xooo ma .=mm< Nmma an aim... an $0 on 5:15 >onkm m_I._. 9. Din? Nmma 15 .296 $2 .2 .294 65 .2 8.55:6 .o 2:2". on mam an as on meB 535 mg a 8:88 8 _ m. DEA? «mm. E a. .22 N 3 Don eon we 3 on ||\l/>\»< ”528 .2 .o .2 EMS: momaom i .5. 2 Irma Na_ mn< .wm .mn .Om Em_Zo 33 .NN .3826: 2: .2 23.5.90; K 0.59“. afim om n=o mm wx_m._.m mmma .XOOU l Duh—.00 n=i_m wk. n=o hm mx.m._.m ._o a0 _xm203u_No l owImoz 44mm— m n__..m om 15 mm MXEHM >095 m_I._. . 42: .Nm Emzm>oz 18 .— .o a .E mic. womDOm .296 39 «a 32532 2: .2 mozofigm .m 952”. m ajm om a5 mm meB i( (Ll.\c< .535 Exp *mma NN mmm2m>oz _ mozoowm on Nom_oz Z ._. Q Q o n.n msv 0d oda 5N m6 0.0 Ca .v.N Tn m6 ”HmDmo Ex 2 Irma New 0x» .m .Om .mm 2m_z<10m2 4?: ll}: g 42:: New 00x ll? 4:33.): Nam .50 19 4:26 39 .NN .38262 2: .3. 62:23:. 68: .a .6 326328 65 .2 3.5554,. .m 2:2“. _ 1:12:44 3.4. New 00x New 0x.» mozoomm on Nom_oz m¢ n._._m 0% m5 hm mx_m._.m mmm— .._D we _¥mzo>>m_No . 20:5qu mOmZMF .FZMEOE QOEZMU mom mo_._..mI._.Z>m OOEMQ .EOIm $2 .8 mmm2u>oz Nam Odd Nam .50 20 these fairly low quality of waveform matches, or by noting that the centroid- moment tensor technique used by Dziewonski et al. (1985) becomes more unreliable for events at shallower depth, longer wavelengths being less reliable (Michael and Geller, 1984). The depth for the event is calculated at 11 km. It is possible that both of the above possiblities contribute to the difference observed. Synthetics were also created for the Dziewonski mechanism (Figure 9), showing relatively poor short-period waveform fits. January 21, 1976 The January 21, 1976 event occured in the northern Cherskii mountains (Figure 2). Previous solutions were created by Savostin and Karasik (1981), Koz’min (1984) and Cook (1988), based on first-motions from WWSSN station records (Figure 10). The mechanism of Koz’min (1984), strike 88 dip 58 slip 354, has a dilatational component and a large strike-slip component. That of Savostin and Karasik (1981), strike 106 dip 34 slip 32, and Cook (1988), strike 90 dip 45 slip 29, both have a slight thrust component and a large strike-slip component. In this study first motions from WWSSN stations were picked and three short-period station records were modelled (Figure 11). Waveform fits of good quality were produced using the mechanism (strike 56 dip 70 slip 11), this mechanism being very similar to the thrust mechamisms of Savostin and Karasik (1981) and Cook (1988). The crustal model used was a half-space of typical continent material, overlain by a layer of lighter material with an arbitrary thickness of one kilometer. The calculated fault plane can be reconciled with faults in the area reported by Zonenshain et al. (1978), Savostin and Karasik (1981), and with fault traces observed on LANDSAT imagery by CDDk (1988). Poor constraint of this mechanism exists due to very few stations used for 21 pump . .826 2.9 .5 33:2. 2: .8 23.5.8 .3 952". mu mfim me 1.0 cm 9315 coma .xooo l 030m #mn n=n_m mm $0 on w¥.m._.m vmmp .z_2.NO¥ l ouImm095 m_I._. mhma .aN >m<32m<3mmw. nmma 8mm : ._xmzo;u.~a onu mam E .5 2... 2.25 20:33 .5sz Ego: 90228 N 25 .296 B2 .2 29.26... 65 .2 8222.6 .2 9:9. mm >m<3mmmu Emm— _ mozooum on _ 2W dune O 8.» Now use . . u . 2 mm 22 52%;? é a Q n a . 5N 0.0 0.0 0._. TN Nd n...“ 0.n 0; 0.0 m; _ ”hm—4&0 . .a .0 a .E mic. momDOm . 20. 2 Same .24.. .2 .62 292.45% Nam |_Im Nam 0v; 26 .296 462 .3 8:34 65 .2 26.2.8 .3 9:9. onN mflm 50 m5 0s... $.st 000— .xooo l Im095 91.. +009 .mm .5303. 27 .805 39 .8 5:93 2: .2 6.2.6556 .2. 2:9“. :8 2.2 on 6.: :3 2.25 _ 2286» on a 5:5 m...» 4om_.mu hmnoa< 28 .296 462 .8 2:63 9: .2 8.622% .2 9:9. one .3. on .5 Sn .2... .62.... 2: _l 62.88 on _ *00~.m~ km303< 2% «ans: .2 q u a ma fin 4.. o.~ ca 4.. ma 2.. m 5.2. .03 en .03 292.418: N9 Dem m... o; 0.: ma . ”5:5 \/\.\/\\, .. Q: o.~ as t .2: momaom . 29 three fault strikes of between 338 and 356 degrees. The dip of these fault planes varies by 30 degrees, while the amount of strike-slip component, that is the slip angle, varies from 270 degrees, (no strike slip motion, Jemsek et al. (1986)), to 247 degrees (slight strike slip motion, Savostin and Karasik (1981)). These mechanisms overall are fairly consistent with each other. In this study first motions were picked and long-period waveform modelling was performed for ten WWSSN stations (Figures 15 and 16). The mechanism produced: strike 310 dip 36 slip 250. provides good waveform fits to the long-period station records. This mechanism resembles most closely that of J emsek et a1. (1986), but with a greater strike-slip component. The crustal model and the source-time function used in studying this event were calculated by Jemsek et al. (1986). This crustal model was somewhat simplified for this study due to software limitations. This mechanism, like the previous mechanisms, fulfills the general expectation that a ridge earthquake have a large extensional component and act on faults aligned sub-parallel with the ridge. July 20, 1976 The July 20, 1976 (mb 4.8) event took place on the Arctic Mid-Ocean ridge very close to where the August 1964 event occured (Figure 1). This event has not been studied previously. First motions were picked and waveform modelling was done for two short period station records (Figures 17 and 18). The resulting fits were fair. From this information, a mechanism can be roughly constrained: strike 320 dip 50 slip 210. This is an extensional mechanism with a considerable strike-slip cOmponent. The error margins 'are wide, with so few stations modelled, that large changes, for example ,1 30 degrees in azimuth 30 .526 2.9 6“ 23.. 2: .2 cozaom .5 059“. >035 mi. ohm? .ON >42. fwd 40m: 31 .225 29 .8 22. 2: .2 822.2% .2 2:9“. 52.5 m=¢ who _. ON >43... mozooum on _ fwd 0m: 4 06 I 32 _ a. q n a ‘ t (:Efi ( z n.n 3. od . _ 0.2 3 Qn 3 o; in fin on ma 3 0.0 m... E; “525 _ to 0.0 to t “5: momaom 5.. 8 Emma Nam .Em Nam om: .OPN .on .ONn 2m_Z035 mih mm? .E 10%.: .226 82. a 5.22 a... .2 8:»...5w .8 2:9“. 52m 9.: _ mam. .R 10%.: 828% on _\<<>s>> 34 oduoE hand 3 m.» Om 0N flu n2 2 g 3 82 0.0 No.2 .5st ca o.m 0N tuzp momaom c... 2 Iran . .03 .8 .22” 2925.82 N3 132 NE O._.n_ 35 April 7, 1969 The event on April 7, 1969 (mb 5.4) occurred on the Eurasian continental shelf, near the northern end of the Omoloi graben, below the northern Laptev Sea (Figure 1). This event has been well-studied by several previous workers (e.g., Chapman and Solomon, 1976; Savostin and Karasik, 1981; Jemsek et al., 1986; Koz’min, 1984; Cook, 1988), and calculated mechanisms agree only that some extensional component is present (Figure 21). It appears that centroid— moment tensor studies (i.e., Chapman and Solomon, 1976) favor a dilatational mechanism with a northwest - southeast striking plane, while first-motion studies favor mechanisms with more strike-slip component, and a more north-south striking fault plane. Faults of the western side of the Omoloi graben may be reconciled with either of these fault plane trends. This study picks first-motion information and matches seven long-period and six short-period station records in order to constrain a focal mechanism (strike 314 dip 62 slip 301; Figures 22 and 23). Waveform fits were very good for both long-period and short-period station records. The source time function used (1.0 0.0 2.0) was adopted as closely as possible from Jemesk et al. (1986), who also modelled long-period records and produced good fits for this event. The crustal model used was a half-space of continental material overlain by a one kilometer thick layer. The good matches for so many station records suggests the error margin for this mechanism is low. The reason for the inconsistency among the previous solutions for this event is not clear. The mechanism of this study generally resembles other solutions derived from first-motion and waveform-matching studies when comparing dip and slip angle. The strike of the present mechanism agrees. with the Jemsek et al. (1986) mechanism (314° azimuth), but is rotated by as much as 66° from that of other 36 . . .ao . 0:6 «mm .5m 2. .3 ca 9.25 .226 89 5:23 on. .8 ace.“ _ m E E moms .xooo I Duh—.00 non mam 0+ n=o up 5:35 pour .v=m(¢(x 02( 3.596 I 813 .vnN mfim to ma con “xi—h chm. .ZOZOJOm 02( 253.210 I 030m h 4.1% mmmp nap cam on m5 2. 9:35 #00" .zSNox I auxmoam m=.:. mmmp .m. ..=mm< z... 37 an aim S “.5 in 9..um .625 mi» mom? K nan? ._. Q n 6 EN m.n n6 0.? TN in Qm "kmamo O.N 0.0 o; E mic. momaom Ev. : Ikmmo .Pon .Nm sin 2923.552 .225 82 K __.._< 2: .2 8:025»... .«u 2:9". _ avenue» om _ .226 23 K 22 2: .2 62622.6 .8 2:2". 5n ajm um a5 In meB >035 m=.:. mwmp K .=mn_< . 828% on _ 39 mechanisms. The consistency of the fault plane between this study and that of Jemsek et al. (1986) suggests this fault plane may be correct. September 22, 1987 The September 22, 1987 event (mb 5.6) occurred on the continental shelf of the Eurasian continent, below the western side of the Bel’kov graben, in the northern Laptev Sea (Figure 1). Dziewonski (1988) publishes a centroid-moment tensor solution for this event: strike 6 dip 41 slip 291 (Figure 24). This study picks first-motion information and matches two long-period and four short-period station records (Figure 25). First-motions were dilatations. Good waveform fits were found for both long-period and short-period records using the mechanism (strike 189 dip 45 slip 270). Either this or the CMT mechanism can be reconciled with normal faults of the Bel’kov graben (Figure 1). Short period modelling was possible using the source time function (1.0 0.2 1.0). August 5, 1986 A small event (mb 4.7) occurred on the continental shelf of Eurasia, north of the New Siberian Islands on August 5, 1986 (Figure 1). This event is the most eastern event in the present study. Little information from this event is found on WWSSN station records. Few first-motion picks were made and no waveform modelling was successful. The focal mechanism is unconstrained. It was observed on some station records that two phases appear consistently and with a fairly consistent difference in time between them. These phases appear ' On the statiOn- records at about the right time to be part of the P wavetrain from this event. These otherwise poor records were digitized on the hope that 40 .226 .82 .«~ .3528 2: .2 5:28 SN. mam .3 ma .m “22.5 38: 29026.3 .. 20:38 :8 ‘w .3. 2:9". >035 m_Ih mm? .mm mmmfimEmm .u:w>m hmmfi .mm Hmnamummm on» How moflumnucmm .mm wusmfim >0 ka m_I._. Ba. .3 mmmzwium I\<, _ 2...... _ .5 ., 3. use 41 z Nam 2.4. .P Q n 6 . man we: 0d 0; TN in m.m ”.528 o; g o; E m:: womaom 22:22.22: 3% 2.. m: Irma Nam .03 .2. .2: 2924.18: .3... mac 2m 42 these phases could later be named as the phases P and pP for this event, and thereby be useful to determine focal depth for this event (Figure 26). If one does assume these phases are P and pP, then assuming a reasonable crustal model (mantle half-space, continental layer), a focal depth of 18 km is calculated. More detailed crustal structure information is available for the New Siberian Islands region (Avetisov, 1983), but is not implemented here because the focal mechanism is unconstrained regardless. No other useful information was obtained for this event. April 23, 1977 An event occurred on April 23, 1977 (mb 5.0) in the eastern Laptev Sea, below the eastern margin . of the Kykov graben (Figure 1). A centroid-moment tensor solution from Dziewonski, (1987) is strike 0 dip 45 slip 270 for this event (Figure 27). In this study, first-motion picks are made and waveform matching of five short-period station records was done (Figure 28). All first-motions are dilatations. The waveform modelling produced fair fits for the mechanism (strike 30 dip 45 slip 270). The observed records were generally noisy. Either this or the CMT mechanisms could be easily reconciled with the normal faults of the north-northeast trending Kykov graben. The crustal model used for study of this event is a half space of continental material created by Jemsek et a1. (1986). The crustal model for the New Siberian Islands region from Avetisov, (1983) is not used. The April, 1977, event is small, so that regardless of crustal structure refinement adequate fits should be obtained. The Jemsek et al. (1986) half-space crustal model was used. 43 .6333 momma: an new a =uEm uofionmam 53. 2.6.6 83 .m 5:92 65 .2 3.80.. 603030 .3 232m _ mozoowm on _ 1 N66 3.. $3. 6; Nam 0.40 , Nam .5: Ex 9 ”Emma .200... 82280 mums... an 024 2 853mg 1:; momooum omzzEmzooz: 2m_z35 9.: hum? mm .554. 2 mmzooum on 2 22.: $2222 Nam 41m Nam 030 mum 2mm. ow i: ( C #mDmo o; 0.0 0.: E 24: momaom 2.. 2 Emma é 2:; . . Nam cum 9» .on 2m_z095 m_I._. mmmr .oe mZDw 48 .225 22 .2 2.2. 2: .2 «25556 .8 2:2“. mnu 2.5 mm ma mm. 9..um 52.5 mi» mm? 0, uzan z p q n a n.n 0... 9m 23.. 0.». o... o... 0:... 3 man o.m ow ON 3 3 od 3 3 Wm upmamo .N .o 2 big momaom 5.. mm :58 .mnN .mm. .mmF 2m_zoDhm mi» _ nmmp .OF mZDfi mozooum on ‘ 50 recorded on at least some stations and in fact, other events of smaller magnitude have been better recorded. Reasons for this difference are not clear. The remaining events under study occurred in the vicinity of the Lena River delta. May 20, 1963 The event of May 20, 1963 (mb 5.0) occurred below the Lena River delta (Figure 1). It was studied previously by Cook (1988) using first-motion information (Figure 32). He offers two possible mechanism solutions (strike 174 dip 60 slip 349, and strike 12 dip 80 slip 180), both mechanisms having large strike-slip components. This study considers first-motion picks and matches four short-period station .records for this event (Figure 33). The resulting mechanism (strike 270 dip 75 slip 339) has a large strike-slip component and is very similar to one of the Cook (1988) solutions. Waveform fits are good. The crustal model assumed for the epicentral region, simplified from Vinogradov (1984) consists of a kilometer-thick layer over a half space. The east-west striking fault planes of these mechanism are not aligned with a prominent graben trend. July 21, 1964 A fairly large event (mb 5.4) occured below Bour Khaya Bay (Bour Khaya graben; Figure 1), to the east of the Lena River delta on July 21, 1964. This event has been studied by Cook (1988) using first-motion information, surface-wave radiation studies and short-period aneform modelling of various station records. The mechanism calculated is (strike 170 dip 45 slip 286), which 51 .596 near .0“ >22 05 .2 325.2% .3 052“. can djm om n:n_ mmw “.3:um 2m mam om a5 03 9..um mnn mam nu as OR 9.5.5 2.2 .xooo >m mzoEjom 53.5 m: 0m .22 mm? BE ON .25. 52 .226 82 .8 .6: o... .2 8.652“... .3 2:2“. mun him 2 .5 2a 9.55 .696 95 fl nmme ON >32 mozoowm on z . F q n a - ....N m.n 0m oé ON ON 0% upmamo Do 0.0 00 t 2...: momaom . c... 2 IE8 3-2,. .24.. .2 .08 295182 Nam 0.3. 53 has a large dilatational component, with a slight strike-slip motion (Figure 34). He utilizes crustal model information from the area of Bour Khaya Bay obtained from a refraction line shot across that area by Kogan (1974). In this study first motion picks are considered and waveform modelling of one long period and six short period station records was done (Figures 35 and 36). The waveform fits are generally good and the mechanism obtained (strike 177 dip 44 slip 255) is in agreement with the Cook (1988) mechanism. The north-south striking fault plane aligns well with the bounding faults of the Bour Khaya graben. February 1, 1980 The final event under study occured on February 1, 1980 (mb 5.4) below the western edge of the Lena River delta (Figure 1), below the Olenek aulacogen. This event was studied previously by Cook (1988) (Figure 37), using first-motion information and surface-wave radiation patterns, producing two mechanisms: strike 274 dip 71 slip 222, and strike 274 dip 72 slip 344. The strike of the fault planes for these mechanisms align well with the east-west striking faults associated with the aulacogen. A centroid-moment tensor solution, strike 315 dip 55 slip 282 by Dziewonski et al. (1988a) (Figure 37) could also be reconciled with probable faults in the area. This study considers first motion picks and models six short-period station records for this event (Figure 38). The resulting mechanism (strike 288 dip 60 slip 228) is in general agreement with the previously offered mechanisms. Waveform fits are good. The crustal model used includes a sediment layer to approximate ' the effect of the nearby Lena River delta. The source-time function used (1.0 0.0 2.0) was adopted from Jemsek et al. (1986). 54 .805 39 .Fu :2. 2: .2 «cozaom .3 2:9“. mmN mjm 3V n:0 E; mv._m:.m >0D._.m m_I._. #mmw .PN >42. mmm n__.._m m¢ n:0 0m; meHm mmm— .xooo HZO_._._.jOm mDO_>mmn_ 55 .526 82 .5 22. 2: .2 8.85:6 .8 2:9". non mam>M¢Ewnth $.55 _ mozoomm on _ {. temp .3 >40... 3 f Nam 42< _ nucooou om _ n.n mé od . oén ad m6 m6 , od 3 ma 3 . o.n Q. 5.. ma (/\/\/)\/ “5an . n.o 0.0 nd tug... momaom 29.22 Q g. N. a... N m >9. N9 205 .mmN .1. NC. 2m_Z0D._.m m_I._. $2 .a .52. _ mozoowm on Nam o>m_N0 .520 I 0.40m Kym. 13m NB n:0 .vhN mxfifim NNN n:u_m K n:0 .vhN mxfifim mmmr .xooo HmZOEjOm mDO.>umn_ F >m02.m m_I._. P >mm<3mmmm owmp mozooumon : : Nam 1.00 Nam Jim I? It? 5:255 Nam 030 Nam Sm 2...: ad fin to 3 ca NWDWM 3 od 0.. E m2: .258 E _ 5:? ;( 5.. «a Emma 63 .8 .mnw 2925.8: Nam 00.. Nam ozm 59 DISCUSSION Calculated and plotted focal depths (Figure 2), show a pattern of shallow focal depth across the Laptev Sea region. Events of the north-central and south-central Laptev Sea have focal depths of 11 km (April, 1969), 12 km (July, 1964), and 10 km (May, 1963). The events of the northwestern Laptev Sea, excluding those north of the continental slope, and of the southwestern Laptev Sea have focal depths of 25 km (June, 1983) and 22 km (February, 1980). Events near the New Siberian Islands (eastern Laptev Sea) have focal depths, given in order from east to west, of 15 km (September, 1987), 17 km (April, 1977) and 18 km (August, 1986). There is a pattern of increasing focal depths both toward the east and. west, away from the 130° meridian. The grabens in the 130° East area and aligned north-south are the Omoloi graben to the north and the Bour Khaya graben to the south (Figure 38). It is interpreted that these represent the current locus of spreading in the Laptev Sea area as suggested by Fujita et al. (1990b) and Kim (1986). The data are sufficient to resolve the current locus of spreading, but do not resolve which style of rifting (pure, simple, or a combination of these) is presently active in the area. This is due to the similarity of heat flow regimes (and therefore shallow earthquake trends) arising from both types of rifting, as modelled by Buck et al. (1988). All events studied occurring beneath the Laptev Sea or to the north include an extensional component, with some strike-slip motion (Figure 1). The events occurring in the Cherskii mountains have considerable strike-slip and compressional components (January, 1976: strike 56 dip 70 slip 11; April, 1962: strike 80 dip 39 slip 58; November, 1984: strike 58 dip 90 slip 6)., This change ' in stress field on going north to south within the study area indicates that 60 some events occurred north of the pole of rotation in this area, and other events occurred south of the pole. However, grabens or fault-bounded depressions extend through the Laptev Sea area and through the Cherskii Range (extending out of the study area). These features, as well as high heat flow, have been cited by Savostin and Karasik (1981) as the basis for concluding that this entire region is presently undergoing extension. This is not the case, based on seismic information. Compression is probably a fairly recent development on a previously extensional margin. The style of rifting in the region of the Laptev Sea is discussed by Fujita et al. (1990b). They describe the characteristics of both the pure and simple shear rifting models. The pure shear model involves a rift which is separating symmetrically about the spreading axis, with faults on either side extending into the ductile zone of the crust (McKenzie, 1978). The simple shear model rift (Wernicke, 1985; Lister et al., 1986) is distinctly asymmetric, separating on a single detachment fault which extends through the crust. This detachment produces an upper and lower plate in the rifting margin. The polarity of this detachment plane may reverse across transfer faults in the region. Fujita et al. (1990b) determine that the simple shear model of rifting is most appropriate for describing the rifting of the Laptev Sea region. GEOLOGIC INTERPRETATIONS In the Laptev Sea region, the Eurasian and North American plate form a two-plate system, with the North American plate rotating clockwise with respect to the Eurasian plate about a pole lying near Bour Khaya Bay, southern Laptev Sea (Cook et al., 1986). Rifting is in progress north of the pole. The northern I Moma Rift system, to the south of the Laptev Sea, is currently undergoing 61 compression. The presence of the Moma Rift system, with high heat flow and Cenozoic bimodal volcanism suggest the northern Cherskii Range area (Moma rift) was recently under extension and possibly connected with the Laptev Sea rift system through the Ust’ Yana graben (Fujita et al., 1990b). The high heat flow of the Moma rift system should be dissipating after extension ended in the area. TECTONIC IMPLICATIONS From this study, all Laptev Sea events occur on high angle faults. Because no low angle normal faulting (detachment plane) events are found, it seems unlikely that the simple shear model is presently active in the area. However arguments from Fujita et al. ( 1990b), based on topography, etc., suggest that simple shear rifting was active in the area at one time. It is possible that rifting in the area began as simple shear motion and has since either changed to pure shear motion or has developed beyond the point where any single mode of rifting is recognizable. The deep focal depths occurring in the western Laptev Sea (22 km, 25 km) are difficult to explain with either pure or simple shear rifting and may support the idea that continental material in the area is being broken up by the extensional stress, and that the Laptev Sea rift system, which is a system in transition from oceanic to continental affinity, may be complex and un-modellable with a single rifting model. The presence of more than one north-south trend in seismicity in the area may be delimiting the edges of various continental blocks. The trend of seismicity farthest to the east, just west of the New Siberian Islands (Bel’kov, Kykov grabens), include nearly pure normal faulting mechanisms and may be occurring on a major . bounding normal fault of the "solid" continent to the east. Other events in the 62 Laptev Sea have greater strike-slip components and may be occurring between various smaller continental blocks in the rift zone. CONCLUSIONS A trend exists in the calculated focal depths for seismic events occurring in the region of the Laptev Sea, with shallower events occuring along the 130°E meridian. This corresponds with and aligns with the entry point of the Mid-Arctic ridge into the Eurasian continent and is thought to be caused by the propagation of the ridge extensional regime into the continent. It is interpreted that the Omoloi graben (north Laptev Sea) and the Bour Khaya graben (south Laptev Sea) comprise the current locus of spreading in this area, and that the Eurasian-North American plate boundary lies along the 130° meridian between 77° and 72° north latitude. Focal mechanisms are extensional in the north of the study area and gradually have more strike-slip motion to the south, with a compressive component existing in the southern-most events. This is interpreted as confirmation that the pole of rotation for Eurasia-North America lies in the study area. 53 REFERENCES CITED Avetisov, G. P., 1983, Seismologic data on the deep structure of the New Siberian Islands and adjacent sea areas, International Geology Review, v. 25, p. 651-660. Buck, W. R., Martinez, F., Steckler, M. S., and Cochran, J. R., 1988. Thermal cogsqugzences of lithospheric extension: Pure and Simple, Tectonics, v. 7, p. 1 - 4. Chapman, M. E., and Solomon, S. C., 1976. North American-Eurasian plate ggpréggry in Northeast Asia, Journal of Geophysical Research, v. 81, p. Chase, C. G., 1978. Plate kinematics: The Americas, East Africa, and the rest of the world, Earth and Planetary Science Letters, v. 37, p. 355-368. Chen, W.-P., and Molnar, P., 1983. Focal depths of intracontinent and intraplate earthquakes and their implications for the thermal and mechanical roperties of the lithosphere, Journal of Geophysical Research, v. 88, p.418 -4214. Cook, D. B., Fujita, K, and McMullen, c. A., 1986. Present-Dagkgflate interactions in Northeast Asia: North American, Eurasian, and otsk plates, Journal of Geodynamics, v. 6, p. 33-51. Cook, D. B., 1988. Seismology and tectonics of the North American late in the Arctic: Northeast Siberia and Alaska, Ph.D. Dissertation, ichigan State University, 250 p. DeMets, C., Gordon, R., Argus, D., and Stein, S., 1990, Current plate motions, Geophysical Journal International, in press. Dziewonski, A. M., Franzen, J. E., and Woodhouse, J. H., 1983. Centroid-moment tensor solutions for April - June, 1983, Physics of the Earth and Planetary Interiors, v. 33, p. 243-249. Dziewonski, A. M., Franzen, J. E., and Woodhouse, J. H., 1985. Centroid-moment tensor solutions for October - December, 1984, Physics of the Earth and Planetary Interiors, v. 39, p. 147-156. Dziewonski, A. M., Ekstrom, G., Franzen, J. E., and Woodhouse, J. H., 1987. Global seismicity of 1977: centroid-moment tensor solutions for 471 earthquakes, Physics of the Earth and Planetary Interiors, v. 45, p. 11-36. Dziewonski, A. M., Ekstrom, G., Franzen, J. E., and Woodhouse, J. H., 1988a. Global seismicity of 1980: centroid-moment tensor solutions for 515 earthquakes, Physics .~ of the Earth and Planetary Interiors, v.‘ 50, p. 127-1 4. ' ‘ 64 Dziewonski, A. M., Ekstrom, G., Woodhouse, J. H., and Zwart, G., 1988b. Centroid-moment tensor solutions for J anua - March 1987, Physics of the Earth and Planetary Interiors, v. 50, p. 116-12 . Dziewonski, A. M., 1988. Centroid-moment tensor solutions for July-September 1987, Physics of the Earth and Planetary Interiors, v. 53, p. 1-11. Fujita, K., Olson, D. R., and Kroeger, G. C., 1990a. The transition from extension to compression in the southern Laptev Sea, Northeast Siberia, Eos (Transactions of the American Geophysical Union), v. 71, p. 641. Fujita, K., Cambray, F. W., and Velbel, M. A., 1990b. Tectonics of the La tev Sea and Moma rift systems, Northeastern USSR, Marine Geology, v. 9 , in press. Gaponenko, G. I., Litinskii, V. A., Levin, D. V., Orlov, A. N., and Pol’kin, Y. I., 1968. Geolo 'c-tectonic structure of the bottom of the Laptev and western part of the ast Siberian seas according to eophysical data (in Russian), Geofizicheskie Metody Razvedki v Arktiki, v. , p. 20-32. Genin, B. L., Lipkov, L. Z., and Piskarev, A. L, 1977. On the structure of the basement of the East Siberian shelf as exemplified b the New Siberian Islands archipelago (in Russian). In: Y. E. Pogrebitskii and M. K. Kos’ko (Editors), Arctic Tectonics, The Folded Basement of Sedimentary Basins. NIIGA, Leningrad, p. 85-97. Helrnberger, D. V., and Burdick, L. J., 1979. Synthetic seismograms, Annual Rev1ew of Earth and Planetary Science, v. 7, p. 417-442. Jemsek, J. P., Bergman, E. A., Nabelek, J. L., and Solomon, S. C., 1986. Focal depths and mechanisms of large earthquakes on the Arctic Mid-Ocean Ridge system, Journal of Geophysical Research, v. 91, p. 13993-14005. Kim, B. 1., 1986. Structural continuation of the rift valley of Gakkel ridge on the Laptev shelf (in Russian). In: B. K. Egiazarov and Y. B. Kazmin (Editors), Structure and Histo of Development of the Arctic Ocean. Sevmorgeologiya, Leningrad, p. 1 3-139. Kogan, A. L., 1974. Performing seismic work of the CRMW-DSS method on oceanic ice on the shelves of Arctic Seas (example of work in the Laptev Sea) (in Russian), Geofizicheskie Metody Razvedki v Arktiki, v. 9, p. 33-38. Koz’min, B. M., 1984. Seismic zones of Yakutia and the focal mechanism of their earthquakes (in Russian). N auka, Moscow, 126 p. Kroeger, G. C., 1987. Synthesis and analysis of teleseismic body wave seismograms, Ph.D. Dissertation, Stanford niversity, 135 p. Lander, A. V., 1990. Dispersion characteristics of Rayleigh waves on the Asiatic Arctic shelf and features of the structure. of the Laptev Sea, Vychilitel’naya Seismologiya, v. 21, p. 152- 157. ' 65 Lister, G. S., Etheridge, M. A., and Symonds, P. A., 1986. Detachment faulting and the evolution of passive continental margins, Geology, v. 14, p. 246-250. McKenzie, D. P., 1978. Some remarks on the develo ment of sedimentary basins, Earth and Planetary Science Letters, v. 40, p. -32. Michael, A. J., and Geller, R. J., 1984. Linear moment tensor inversion for shallow thrust earthquakes combining first-motion and surface wave data, Journal of Geophysical Research, v. 89, p. 1889-1897. Minster, J. B., and Jordan, T. H., 1978. Present-day plate motions, Journal of Geophysical Research, v. 83, p. 5331-5354. Parfenov, L. M., Koz’min, B. M., Grinenko, O. V., Imaev, V. S., and Imaeva, L P., 1988. Geodynamics of the Chersky seismic belt. Journal of Geodynamics, v. 9, p. 15-37. Savostin, L A., and Karasik, A. M., 1981. Recent plate tectonics of the Arctic basin and of Northeastern Asia, v. 74, p. 111-145. Stein, S., and Kroeger, G., 1980. Estimating earthquake source parameters from seismolo 'cal data. In: S. Nemat-Nasser (Editor), Solid Earth Geophysics an Geotechnology, AMD, v. 42, p. 61-71. Stein, S., and Wiens, D. A., 1986. Depth determination for shallow teleseismic gb‘arthé sakes: methods and results, Review of Geophysics, v. 24, p. Sykes, L. R., 1967. The Seismici of the Arctic, Bulletin of the Seismological Society of America, v. 55, 501- 18. Vino adov, V. A., 1984. Laptev Sea (in Russian). In: I. S. Gramber and Y. E. o rebitskii (Editors), Geologic Structure of the USS and its Re ationshi to the Distribution of Mineral Resources, v. 9, Seas of the Soviet Al’Cth. Nedra, Leningrad, p. 51-60. Wernicke, B., 1985. Low-angle normal faults in the basin and ran e province: nappe tectonics in an extending orogen, Nature, v. 291, p. 645-64 . Wiens, D. A., 1987. Effects of near source bathymetry on teleseismic P waveforms, Geophysical Research Letters, v. 14, p. 761-764. Zonenshain, L P., Natapov, L M., Savostin, L A., and Stavsky, A. P., 1978. Recent plate tectonics of Northeastern Asia in connection with the opening of the Northern Atlantic and Arctic Basins, Oceanology, v. 18, p. 846-85 . NSTRTE UNI V. IIIIIIIIIIIIIII I IIII I II IIIIIIIIIIIIIIIII