<<

View metadata, citation and similar papers at core.ac.uk brought to you by CORE

provided by Kyoto University Research Information Repository

Characteristics of Seismicity Distribution along the : Title Some New Observations

Author(s) GHOSE, Ranajit; OIKE, Kazuo

Bulletin of the Disaster Prevention Research Institute (1988), Citation 38(2): 29-48

Issue Date 1988-06

URL http://hdl.handle.net/2433/124954

Right

Type Departmental Bulletin Paper

Textversion publisher

Kyoto University Bull. Disas. Prey. Res. Inst., Kyoto Univ., Vol. 38, Part 2, No. 332, June, 1988 29

Characteristics of Seismicity Distribution along the Sunda Arc: Some New Observations

By Ranajit GHOSEand Kazuo OIKE

(Manuscript received March 7, 1988)

Abstract

Spatio-temporal variations of activity along the Sunda arc were investigated. We prepared a strain release map for this century. Adjacent to the zones of high strain release, presence of seismically quiet zones was noted. A careful inspection of the depth distribution of the revealed that in the eastern Sunda arc, possibly there exists a zone of scarce seismicity at an interme- diate depth. We discussed the probable implications. We also analysed the patterns of temporal distributions of earthquakes at the three different seismotectonic provinces of the Sunda arc—, , and the Lesser Sunda . We could clearly see that, although the causative geodynamic situations for seismicity vary significantly in space along the length of the arc, the period of increase or decrease in seismicity is largely space invariant. The locally differing levels of seismicity are superposed on the common background of long period seismicity fluctuation. Finally, clustering of seismicity at some patches along the Sunda arc was studied with respect to the altimetric gravity anomaly data. We noted some apparent conformities. For some of our present observations, there exists much scope for detailed future investigations, which we discussed.

Introduction

We analyzed the seismicity data along the entire stretch of the Sunda arc which is unfurled in the southeast Asean region between Andaman-Nicobar Islands in the northwest and in the east. The total length of the arc is over 5600 km and it represents all the features of a typical arc where an oceanic plate (here the -Australian plate) subducts below a continental plate (here the south- east Asean flank of the Eurasian plate). Either in a local scale or regionally, the nature of distribution of earthquakes, volcanoes, gravity and magnetic anomalies, petrology, geochronology and palaeomagnetism of rock samples, crustal seismic velocity structures determined using artificial and/or natural seismic sources, as well as the focal mechanisms of earthquakes and the derived stress distribution in the descending oceanic slab and adjacent regions have been investigated and reported by many researchers, references of any or all of which we refrain from citing here for the sake of brevity. Recently, Newcomb and McCann (1987)1' has reported of the seismic history of the Sunda arc region. Fig. 1 shows the plate tectonic setting and bathymetry of the Sunda arc region. Many of the morphologic and geodynamic features of the Sunda arc change sig- nificantly along its strike. Minster and Jordan's (1978)2' RM2 model proposed a plate convergence rate of 6 cm/yr near the northern tip of Sumatra to 7.8 cm/yr near in the east. The age of the youngest subducting at Sumatra is 46 30 R. GHOSE and K. OIKE

/00t.f:i A ' t:. 3.. 1•. E.._ s..__AN‘.:1-ITT__, .' f ---- .13gal"v...,Arr-F.ioh-FLAT '41 i J:: -• '''`k 5cV_S.F.S. .. ,, ' .•: •. A .. ,....41._E2C . NcL3...Ai'. ,,,; . SIMEULUE,"''•<'..) .:.4(,c,--- -BIM ,. 'v.. . ,0,`....(;) cja_NIASBARI0 PINT• 1..4. ,„ ••T CEANA... NT) ------\--°.'"r----"";:-."

.. .

SIKH'-k,,1...C).R\PLATET...,.. ... ENGOAN° -- --\\\ •`-, i 0,-, -'--? .,.?.' N'Baia_a a.

59,------•-VIi-"?)•„„,6r{ Ø indian Ocean n vit.,,,cr__,,A,,,s3,--,%q r _, ..k .:,,•4111166'.4,-,!t_-6/=, -V — li a1,ft VcLN 1. ----P,r---i?':"44S1%11Wri- C...t,6-----.\... .,-,Lfia,-5 .11 if 11 0:--, _ ,__2!...,,-.rf---"- 5' 00.-0,e'' 544._zs_\j,,)• lt2J 53 1 1 5(ib '') fi t.)31 [ [ • [ .R. 1 „CAL--I I/ 90° 100° 110° 120° Fig. 1, Plate tectonic setting and bathymetry in the Sunda arc region (modified after New- comb and McCann, 19871)). Black dots denote the sesimic refraction profile of Curray et al. (1977)0. The fine dotted line shows one seismic reflection profile from New- comb and Mccann (1987). S.F.S. stands for Sumatra Fault System. C. IS., R.R., and J.R. are Christmas Island, Roo Rise and Joey Rise respectively. The bathy- metric contours are in km.

Ma while the oldest oceanic crust near Java and Lesser is 152 Ma old. The major morphostructural units can be recognised all along the margin of Sumatra and Java. Fore arc geometry systematically varies from west to east (the depth of the fore arc basins, the trench slope break and the trench increase towards Java) as the sedi- ment thickness on the subducting oceanic plate decreases. A relatively simple style of occurs where the oceanic crust is subducted beneath the continental platform of Sumatra and western Java as well as the of eastern Java and the western (Newcomb and McCann, 1987)". More complex tectonic systems exist in the Andaman Islands, where convergence is essentially subparallel to the arc (Curray et al., 1978)", and east of Java, where collision of the island arc of the Lesser Sunda Islands and the Australian occurs (Silver et al., 1983)4'. , which separates Sumatra from Java serves as a major boundary between two different geodynamic regimes. The trench itself is curved landward west of the Sunda Strait, so that the Strait appears to be a transition zone between two different regions (Budiarto, 1976" ; Ranneft, 19796' ; Zen, 1983"). An abrupt change in the depth of penetration of the Benioff zone between east and west of the Strait is

_ o Characteristicsof Seismicity Distribution along the Sunda Arc: Some New Observations 31 very significant. West of Sunda Strait the seismic activity does not exceed below 200 km, while in the east there is a sudden increase in the maximum depth of seismicity to 650 km with a gap between 300 and 500 km. Le Pichon and Hiertzler (1968)9', Fitch (1970)10),and Fitch and Molner (1970)") discussed the possible explanations for this; however for lack of further substancial evidences equivocality still remains. Interplate motion, normal to the island arc near Java, becomes oblique at Sumatra, where motion parallel to the arc is accomodated by dextral strike-slip displacement along the 1600 km long Sumatra (also known as Semangko or Bati) Fault System, shown in Fig. I_ (for details, refer to Fitch, 197211'; Katili, 197419); Tija, 1978'4' ; Curray et al., 19783); Huchon and Le Pichon, 1984"), among others). With this background knowledge and using updated data base we assessed the seismicity in the whole of the Sunda arc region. Spatio-temporal variation in - quake activity was investigated. In order to obtain an insight into the process of how the heterogeneity in the horizontal compressive stress (coupling) between plates across the subduction boundary is related to the local or regional morphologic or tectonic features, we prepared a regional strain release map, using earthquake data of this century. A careful reinspection of the depth distribution of earthquakes revealed some hitherto unreported peculiarities. We also noticed a remarkable resemblance in the patterns of temporal variations of seismicity in the three different seismotectonic provinces of the Sunda arc Sumatra, Java and the Lesser Sunda Islands. It clearly shows that although the causative geodynamic situations for seismicity vary significantly in space, for the entire Sunda convergence zone, the period of increase or decrease in seismic activity is lagely space invariant. As if, the varying local levels of seismicity are superposed on a common background of long period seismicity fluctuation. Further to these observations, clustering of seismicity at some patches along the Sunda arc has been studied with respect to the gravity anomaly data extracted from the satellite altimetry. Some interesting conformities could be noticed. In this paper, we report of all these observations, and discussed some probable implications. It is expected that enhanced understanding in spatial and temporal changes in seismic activity along the Sunda arc will also provide a better background to identify the regional distinctive features of .

Data

The aim of our present study was to assess the spatial and temporal distribution of seismic activity along the Sunda arc. For this purpose we used two data bases, viz. the NOAA Hypocentre Data File (from Jan., 1900 to May, 1981) and the 1SC Data File (from Jan., 1971 to Dec., 1983). The NOAA Hypocentre Data File includes all the earthquakes recorded by at least five stations. The coverage of the earthquake activity in this data file is somewhat uneven. Only the large events have been reported prior to 1922 because of poor seismographic coverage, while since 1961 most earthquakes with magnitude > 4.5 are included. The completeness of the 32 R. CHOSE and K. DIKE

NOAA hypocentre Data File was investigated by Perez and Scholz (1984)16'. In order to ensure the possibility that data coverage is sufficient to identify temporal changes in seismicity, we used only large events with magnitude > 6, for which the NOAA Data File, excepting possibly a few earlier years of this century, is more or less com- plete. For the large events we also consulted the recent catalogue of Newcomb and McCann (1987)1', with some of the epicentres relocated. For examining the spatial distribution of earthquakes, in relation to volcanicity, bathymetry, and other tectonic and structural features, we used the ISC Data File. We used only those events located by 15 or more station reportings. Out of all hypocentral parameters, the depth estimates are perennially most unreliable. We accepted all the reported pP-P depth estimates to minimize any major depth bias. As the possibility of spurious identification of pwP phases as pP (resulting in deeper depths) still remains, we were careful to keep enough allowance for errors in depth estimates before reaching any conclusion based on them.

Seismicity and Regional Strain Release

In Fig. 2a and 2b the earthquakes from respectively the NOAA and the ISC Data Files were plotted. The ISC events of Fig. 2b are those events which have been located by 25 or more station readings, and hence the epicentral locations are reliable. We can easily recognize some seismically more active patches. We reanalysed the regional seismicity distribution with the background of the hitherto known facts. In Sumatra, the subduction process has a much younger history than Java and Lesser Sunda Islands, where the subducting oceanic crust is very old, and the slab has penetrated much deeper down, to generate some of the deepest earthquake events with downdip compressional mechanism. The shallowest earth- quakes (0-20 km) immediately adjacent to the are commonly associated with strain release occurring within the subducted lithosphere. At Sumatra, the interplate seismic coupling is very high, and thus causes the large thrust events. The majority of strong earthquakes in both the historic and instrumental catalogues of the Sunda arc are located in the fore arc of Sumatra. While Nias and Enggano (shown in Fig. 1) were documented as the locus of activity in the historic record, moderate and large instrumentally recorded earthquakes group principally in central Sumatra, on the Pini Arch, and in southern Sumatra, at Enganno (Newcomb and McCann, 1987", and also visible in Fig. 2c here). The region near Enggano seems to have experienced a large earthquake every 27 years, on an average, since 1756. Besides the established evidence of seismicity related to the back arc thrusting north of fibres, in recent years support continues to grow for the importance of back arc thrusting north of the trough, along which the Australian is being subducted northward beneath the eastern Sunda arc (Silver et al., 19834'; McCaffrey et al., 1985'7' ; McCaffrey and Nabelek, 1984"), 1986m). Another recent explanation of a local seismicity has been presented by Spence (1986)2°' who argues Characteristicsof SeismicityDistribution along the Sunda Arc: Some New Observations 33

'* 4NOAA -., ..•... ••--"•..,.• t,i. *'9 .. t•,, .'''`.--,.7'i.•** ' "t ,„, i.4.*:.0--,,i'----11—',,.• ..••\\\\C"."\\,.\\\ .,.,:i 4.-L.- • ,.) .'''' Fa-3'4'."-,fit; 0,40;„,,„4.,,,d, .f4, -iirI,,...2v,....,,. ._,. _;, I.-...... '...... _4,..4. ..‹, , ,_..,C,„(°*—*)g po\----.1_,,j n.•.4„. „..,.. • .. o r• , # ;14 • Q, •.6, , tif' d3It . ''^.. ,....; L.- • • 4i •..r. 4., a.7.4;r . . • fr _ ,r..,t # • ..4.*.p..••,. 1'' +-*•L...A•7„.:,),0...,,,AG,..,7.. 10 ',., N,-- ,'3,••if-f,... •.Iry, ..-e\--..., =.--..-4,,,, ,,;c---4.'-3, .- .,".." .., ;to* a .• e ' 100 120

.• 46 it,. , =le o o ) 7 ,,, pt ° : •. -• . „..., 401,, ,.?i.„ --V, ,,,,, 1 ..9 .1:''.,-..li..,!Aoaloe 0-9''''Ala 411 ill -I ,;;C° .•4:siti.•' ...',%4.4 4.... • ' .. . , 01 ..... * .°. o .1 Cji 1 CZ: Qr't* .° a ';”:-. . ' C‘t . • 41 .p . . .„ • .. . t° !..1' . , It.....° • All, . • • " , 4, ' • ,.. 4. 9a0,. . • 6 i...• --.:.,4..r*"i;;C--. .:.,,,,„4:•..:i....'= 1.•°4r. I.lt.la,iiiii..',...4,...... ,1. ..--—:::** 10.-...-4. --^ -.--- „ b „r.,44.,.... ;;.wi.- : ejoist.. . • 084 100 120 Fig. 2a. Distributionof earthquakesin the Sunda arc regionduring 1900—May,1981, from NOAAHypocentre Data File. Fig. 2b. Distributionof earthquakes in the Sunda arc regionduring 1971-1983,located by 25 or more stationreadings, from ISC Data File. that the great (M 7.9) Sumba earthquake of August 9, 1977, which is the largest normal faulting event to occur since the Sanriku event of 1933, supports slab-pull as the dominant process for subduction of the oceanic crust beneath the easternmost section of the Sunda trench. In Figs. 2a and 2b we can identify the localities with maximum earthquake frequencies. They are located at the northern tip of Sumatra, north of Simeulue island, southeast of Nias (near Pini), southern tip of Siberut island, near Enggano, southern end of Sumatra near Sunda Strait; and at Java and Lesser Sunda Islands,

• 34 R. GHOSE and K. DIKE

. \ l 1 K CI NOAAIi''''atr, iziV-A.N,y4-.*::..;--.4 *CI%I's1s-\,/v\-1_, M>/6vice•-,8 _,a-br:,t.',7'',4 0 '4PI.,.1.*\I-(' ,-/'Se.1-ft,(-4Th ',, .f:46,:T•111 1i..,-,,..lirwi ,*Da .__,, ,-,,,,,N. '•7APO:M.Eigir-P'1•N7-80 ...-ii-mac+ -.3 J,,,, \illfflairm,:

i....IPIke.,-..1vetoc, -o';do' •^^(-)1,. ° .n Iav-alrimmitli.,'A't6 l'4/05.161r,:, msi,!-,:airExtv.4f----mlislo:..,;i-'.•-;At —„.••.161i 0 -7 KM 0 4',',*t'''---:',^..L,.;-:,\,..,'--),.:,C)_,Ei0, \t*elkdill&4-:-eilWs. ._7--„_,.--5.),,'1k'7 .j41,u.-girt V/ 7T'ttAV/41% — 140,ni 0-T T. r)itx...,,,,....INI - A7 .•itAv,7:4280r-' ( C -7,;=7,wekli',i& 11;: ',;.' <> .-/...‘ - . ,e,14* \ AI OnV01-4N1*":ik.t;al 420k, ()•Z-WiE1106' 10. _ _,./._, r-d-okv.aiaiett...:,..itto#4000,.-..-',aiir560-1L'i c ,Di ,0,z,,,L7 .._L.,•4,,,a, ioo-I

lot] 120 U ', --NOAAA^,-.<,•4 ', r4%-;-7-Jr r\---') )'-,8 h150KM,_/,--._-.--..i.,.o.w,,•.14...,„M1-\ ) N.. , _-,47n ,-,(it\,, \-,v.r,,-„..,+ c,,---),,s,1i- ii..L,,ig,I7,. , 7-8C) ,--4CI, ),-- \4 -1\ Cl.,-,,, . - i'---",,/,?\IAnit4r,.,--4, .••71Ms 11 in,^11.0 6-7 0-.°--,1Av.,'.may; ,,e I ,, _.., ,„,.....:,,,17,-„?..r,; 5-6 0 sq,_)----,\4-5 ,)w...6p L1.-',I:)- , c'-'9'' .,1-- ,.•-''' .---, C''0-4 ,'----•-•_4(.L,r,u r)717 ' .fr,...111^454 _, KM150 co .`r;j,..,6./;,.—)I,. ___,))- a ;.',44•;,;,'YA: / *?f, a ---,ri.11 --11 f - ...5r vV"..,'"I,260-_i ....1.,i0.4-0.....-. i!....*,v•404.*,41127,4ii.P.,/ '- > ID .--.. „*„.1---,-,..•370--480j -,..._., ., d590- 100 120 HO: V

Fig. 2c. Distribution of large earthquakes with magnitude 6 in the Sunda arc region during 1900-May, 1981, from NOAA Data File. Fig. 2d. Distribution of earthquakes deeper than 150 km in the Sunda arc region during 1900-May, 1981, from NOAA Data F le.

near the western tip of Java, south off Java around 107°E and 110°.5E, and southwest of Sumba (the region of locally enhanced slab-pull force). It is striking that at Sumatra, although the plates are strongly coupled, the shallow seismicity does not extend much into the continental platform, whereas whole of Java experiences many such shallow events. The difference in slab penetration between Sumatra and Java cannot completely account for this difference in shallow seismicity inlands of Sumatra and Java. On the basis of our observation, we suggest two possible factors. Newcomb and McCann (1987)" proposed seismic and aseismic modes of subduction for Sumatra and Java respectively. We speculate that as most of the interplate stress build-up is released right at the outer arc thrust zone for Sumatra, and as much of the seismicity observed inland Sumatra is caused by the movement of the Sumatra Fault System,

. I : , Characteristicsof Seismicity Distribution along the Sunda Arc : SomeNew Observations 35 there is possibly not much residual interplate stress related to subduction, left to be transmitted for causing the shallow events in the overriding plate, north of the volcanic front. In contrary, in case of a decoupled arc segment like Java, although the inter- plate stress accumulation is less, the stress release in the thrust zone is also relatively much less. As a result, the net residual stress transmitted to the continental side might be greater in the decoupled portion (Java) than in the coupled one (Sumatra). Secondly, the distance from the trench to inland Sumatra exceeds that to Java. As a matter of fact, the line of active volcanoes (which usually corresponds to the 100 km depth contour for the inclined seismic zone, notable in Fig. 4) passes through the centre (E—W)of Java, whereas for Sumatra it runs close to the Indian ocean coast. This also contributes to the greater seismic activity in the shallow continental part north of Java, compared to Sumatra. Another feature which draws attention in Figs. la and lb is that, out of many scattered localities of enhanced seismic activity and of relative seismic quiescence, just west of 110°E longitude, there is a remarkably silent zone, south off Java. This seismically quiet zone is about 75 km in width and trends almost N—Sfrom the trench to the inland Java (Figs. la and lb). Although minor seismic activity (not visible in our data base) might be present here, absence of earthquake with magnitude > 4, even in the inner thrust zone is certainly noteworthy. Relevantly, a buoyant oceanic high—the Roo Rise is subducting just south of this locality, and Curray et al. (1977)83 also noticed, in their crustal refraction section, some irregularities here. In Fig. 2c all events with magnitude  6 from the NOAA Data File have been plotted. Fig. 2d shows all the events with focal depth  150 km. We can clearly see that the number of large interplate events in the thrust zone off Sumatra is much greater than that for Java. This certainly reflects greater coupling at Sumatra. However, two earthquakes with magnitude _ 8 occurred south of Java, indicating local zones of high horizontal compressive stress. For a better quantitative estimate of the seismic activity along the Sunda arc we prepared a strain release map for the entire region (numerical estimates are presented on maps in the Appendix, in four parts—A, B, C and D). We followed the method of Allen et al. (1965) "'. Each 2° X 2° square was divided into 25 small squares, each of 24' x 24' dimension. The map area (denoted by the boundary of Fig. 3) was divided into 2850 such twentyfour minute squares, each of which was treated as a unit. The number of earthquakes occurring within a unit block was summed up for the entire period. The strain energy released by all the earthquakes in each 24' x 24' square was calculated and converted into an equivalent number of earthquakes (N4) of magnitude M-4. The equation used for this purpose 1S: N4= 100•""4-4) where M is the magnitude of the corresponding earthquake. These numbers were plotted at the centre of each small square. The strain energy release map thus pre- pared is given in the Appendix (in four parts). In stead of contouring the values, we 36 R. GHOSE and K. OIKE

10 STRAIN RELEASE MAP

91 5— Symbols: ^ VERYLOW STRAIN RELEASE ^ ^ 175 < N4 < 1000 0— : ®N4 1000 41 0 1.111111111111,

5- : • 95 0 0 MO0 01 ^ ^ 00,, 10 • ":. 0.• 0 • 0 0• 0 (E1 100 105^ ^mg

14 110 120 r 130 Fig. 3. Strainrelease map of the Sundaarc region. The dottedline demarcatethe boundary of magnitude4 earthquakes. Detailedstrain release figuresare presented in the Appendix. plotted in Fig. 3 the extreme values in symbols. For computing the strain release figures, we used the NOAA (1900 to May, 1981) and the ISC (1971 to 1983, for events with at least 15 station reportings) Data Files. The accuracy of NOAA epicentral locations is not always very good. Also for the earlier years of this century, only the very large events have been reported. Further, the magnitude scale has also changed with time. For these difficulties, we carefully examined each and every event within a unit block. We checked with the recently developed catalogue (with some of the large events relocated and magnitudes homo- geneized) of Newcomb and McCann (1987)1)and minimized other data base incoher- encies. For the smaller events the strain release is much less. For example, an earthquake of magnitude 8 releases strain energy equivalent to almost 1000 magnitude 4 earthquakes, for an event of magnitude 7 the equivalent number of magnitude 4 event is about 178, for magnitude 6 it is about 32, and so on. So missing of a few smaller events will change the strain release picture very little. For the large events our data base was more or less complete. The resulting strain release map is thus quite realistic. In Fig. 3 the hollow squares represent locations with very low strain release (say, less than 10 to 15 N4). Areas within the dotted boundaries, but with no symbol are those of intermediate strain release. A very interesting pattern which emerges in Fig. 3 (see also the Appendix) is that, adjacent to a very large strain release, there is almost always a zone of low strain release. This is more evident in Sumatra where more energy is being released seismically. Whenever a large earthquake occurs, a Characteristicsof Seismicity Distribution along the Sunda Arc: SomeNew Observations 37 significant rupture accompanies. Tocher (1958)22', Iida (1965)23', Ambraseys and Tchalenko (1968) 24), Bonilla and Buchanan (1970)2", and Slemmons (1977)28),among many others, have studied the relationships between fault length and/or surface area and, magnitudes. Acharya (1979)") investigated the regional variations in the rupture- length to magnitude relationships. For Sunda arc no such empirical relationship is yet established. However, Newcomb and McCann (1987)1) delineated from the intensity maps the rupture zones for some large events along the Sunda arc. It appears that many of the very low strain release localities demarcated in Fig. 3 may be located within such rupture zones of very large events, and hence now being in a stage of strain accumulation, are quiet. As a few of these quiet zones may, however, be potential seismic gaps, more detailed and accurate spatio-temporal investigations are essential for understanding their true natures. Habermann et al. (1986)2" explained the low seismicity at two regions south of Sumatra (note in Fig. 3 at around 2°N, 97°.5E and 3°S, 100°.5E) to be possibly due to the interaction with and subduction below Sumatra of a migrating oceanic ridge (Investigator Ridge). The aseismicity of Java compared to Sumatra, excepting for two restricted spots (where very large strain energy was released in the past), is more prominent in our strain release map. In Lesser Sunda Islands, the strain release is again quite high, possibly in part due to the gradually increasing influence of the arc-continent collision. The very high strain release at the extreme right corner of Fig. 3 is associated with the active subduction in the Banda arc.

Depth Distribution of Earthquakes

From ISC Data File, earthquake events between 1971 and 1983 with 15 or more station reportings were plotted in vertical sections of Fig. 4 using the computer pro- gram of Ishikawa et al. (1985)29'. The profile lines are shown in the upper part. All the previously reported features of the Benioff zone in this region are present in these updated depth sections. We also constructed the vertical sections perpendicular to the profile lines and looked from the oceanic side to the island. For the last three profiles —MM', NN' and 00', a notable feature could be traced, which is shown in Fig. 5. In this region the lack of continuity in seismicity between 300 and 500 km is well known since long time. However, as can be seen in Fig. 5, at around a depth range of 120 to 150 km, a zone of scarce seismicity exists. The feature is a consistent one. The density of earthquakes above and below this zone is somewhat greater. The feature does not appear to be just an artifact of inaccurate depth estimates. We suggest the possibility of a zone of stress transition at this depth. In our separate work (Ghose and Oike, in preparation), we have noticed that some of the earthquake focal mechanisms around this depth range show almost opposite orientations of the axes of maximum compressive stress, particularly in the eastern side. Isacks and Molner (1969)3", at first, talked about the possible existence of a stress minima at an intermediate depth 38 R. GHOSE and K. OIKE

(at a different depth range and situation than the present feature). The present feature deserves attention. However, we must admit that our observation, at present,

-..\ \ 1 1 ri ..Th51/_ * 441t, c. ' ,1 flv,1:4-01...... , 0 I 1 r ' iR17di 4004 # ,4146.41iNgisti'y l'1 , . ' C:i'e 10 —wilimi. •dr/19'IIot>"yr _ 1 I ' I _ 100 120 0Ay vA' E y v_ E' Iy vI' Myv M' 140.g, .,44,11'*Ifer" T - N, 280-n, o.. I...114?.41*74T— , 420 - 560 - MO- i 4., B y vB' F y v F' 4y 4`"gt v 41iv_ v11' !::,. o 0:l . Fr- Iv1 1.f,Ar' i C Y C' Gy v G'Ky v K' 0 y v 01 ado '12 -1 t,r,:''''''t4/7'',°0,'•`IV: . c •- 'cil

J

Dy D' Hy v HiLy v V SCALE

- HORIZ.:VERT. -1:1 ,,

. , V —VOLCANICFRONT - Y —TRENCH LOCATION

. ' Fig. 4. Benioff zone structure along the Sunda arc. Earthquakes reported by 15 or more stations from the ISC Data File (1971-1983) are plotted against depth. All reported pP—P depths are used.

f 4 , G 1 1 0 ' Characteristicsof Seismicity Distribution along the SundaArc: SomeNew Observations 39

KM M N 0

o, • r 0— ••• - . • 016o •f.vie•a4.•.11.. - . , 0)00 ;'rfCC4tSi `0. 0ce.0"*e 00 0 6;o 060) • a 0 ° 140-0ef420 0 °o 0 0 00 0 0.0— ••• .•*0°CP0vo0

280-40°C0 ,8D 0oco% ° 0 0

0 0 420- 0 0 0 560- 0 0 p*: 0 op 0 coo;. cop0‘°0 0 60 0 700-

Fig. 5. Seismicitydepth sectionsfor the three easternmostblocks (with profiles MM', NN' and 00') of Fig. 4. The depth distributionsare whilelooking from the trench to- wardsthe islands,along the respectiveprofile lines of Fig. 4. Positionsof thesouthern ends of the profileslines—M, N and 0 are shown. The scarceseismicity at an inter- mediaiedepth rangeis marked by the tiny arrows. is a preliminary one which awaits further careful research.

Temporal Distribution of Seismicity

It is well known that the subduction dynamics, interplate coupling, as well as the resulting seismicity exhibit great variation along the length of the Sunda arc. Local heterogeneities, bathymetric undulations on the subducting oceanic plate and its age also vary significantly. Plate tectonic set ups also differ; at Sumatra, oblique subduction between two strongly coupled plates is believed to be occurring, causing transcurrent faulting and much of the slip seismically; at Java, it is normal subduction between two very weakly coupled plate segments, causing more aseismic slip; whereas at Lesser Sunda Islands, the effect of arc-continent collision is superposed on that of the subduction. In this situation, we were interested to see the patterns of temporal variation of seismicity in these three regions, in order to understand whether the locally varying subduction dynamics and tectonic factors dominate the pattern, or the general regional subductional stress. Fig. 6 shows the temporal distribution of earthquakes with magnitude 6 or over of this century (1900—May, 1981). We can clearly see that the overall patterns of distribution of number (N) of events with magnitude > 6 against time have a remarkable similarity for the three different tectonic provinces. Till 1935-40, the trend was rising. Then from 1935-40 to 40 R. CHOSE and K. OIKE

97 g LTMA.T'' It A.a A."V..z, LESSER 8'8;81 T31,4r) A MM 1,IEMI, 1II 17'1 7' '111'1)1111I11,i11 , 1 1 1 , 1 e II1III,1111111;1 IIII6' 6'1 1 111 L'i 1111 1900 1920 1940 1960 1980 1900 1920 1940 1960 19 OV001920 1940 1960 19:0 80TSUMATRANA'1'FtV A71/40FJAVAV'A.5U7. LESSER I ,-NI1---fIMInp.- N--r401S N I 50 6030i.-' ,,30_ , 40 ,I 20 _,-, , ,--, x 30- /-,1 In 201,-Fr. 20 i 10, _I ,i----- 1 rri',-,- 10I101 1 _CI' 1900 1920 1940' 1960 198019001920 1940 1960 198019001920 1940 1960 1980 Year Year Year Fig. 6. Temporaldistribution of seismicity in the Sunda arc region,using NOAA Data File (1900—May,1981). Only the largeevents with magnitude >6 have beenconsidered. In the upper part, magnitudeis plotted againsttime, while in the lower figures, cummulativenumber of the eventsis plottedagainst time. around 1960, for approximately 20 years, there were few large earthquakes in all of the three segments of Sunda arc. Subsequent to this till May, 1981, the seismicity was again high. (For Lesser Sunda Islands, the sharp rise around 1977 is obviously because of the large Sumba earthquake series.) From this marked similarity in the temporal pattern of seismicity, it is evident that the variation of the common tectonic force (viz., that due to the subduction of the Indian Ocean- Australian plate beneath the southeast Asia) with time, for all of the three segments of the Sunda arc was same, and that dictated the overall long period temporal pattern. Over this common long period fluctuation, the locally varying effects of subbduction dynamics, local plate tectonics and interplate coupling at the three individual segments are superposed, to give the ultimate picture (Fig. 6). It is possible to conclude that, although the causative geodynamic situations vary significantly in space, for the entire Sunda arc, the periods of increase or decrease in seismicity is largely space invariant. As if the varying local levels of seismicity mingles with a common long period seismicity fluctuation.

Seismicity and Altimetric Gravity Anomaly In order to examine the seismicity of the Sunda arc in relation to the gravity anomaly, we made use of the satellite altimetric free air gravity anomaly data of Segawa and Matsumoto (1987)3". We traced the free air gravity anomaly map derived

_ - Characteristicsof SeismicityDistribution along the Sunda Arc : Some New Observations 41

10ity . fC1 ,4) ..;,4,(,_,- ,-,`;Illb Free Air Gravity .) i,,7:. AnomalyDerived From GEOS-3 41+k,,.;;Orlr,,-',/7. (),4--i;;,.;,,,L,...cv.----,,,'4_,.__ andSEASAT-1 Altimetry -' ,"4•..-' 'S.' _. .., :',_-., w ,,

Q ,j,i...•'- \t'P'qtyf' ,,I,J,:,: ..7 n,,.,''''‘..': ,...- l'._... 4 "f..--- 114 ,:',;0\^;:s.'- ess,; .-•,....,,,,„, ,.„,...,_.,,,,,,.,::,,,,,,,•,,,:,,,, 4,'' ,,' 0\d--',,,()sv.,13/4,..,---),,„,'r.-- ,,',,"1_,,,,,,,„.,,_,.,_ ,,.„4:,'",!' ,,; ----,,,J,,, . •,..- i,‘,---0).:, .-:,_-,-,--,,,,,,,,-•,-;-)i--,,,,,,,-__-_,!:,1-,, i?.,,,,,-, J , c'sc:' "c' 11'4•:1=.;,-7,-\t'<-'`'''',t?,-',':,. ter.!,''^';,-.:;;., t:,.,A:,_ 0. ,.,-1-..v.,:-(41-`: __;;'4'_,,,...."7

120 i' ', i, ' ' ,..),,'', , ',',;', '144144",..._4 ,...',,„1''''' Ir.,^,...-- —*-..----'---:-'-`--._:---,•_'-'",'"1,----:-al'''5,t,----,..-- , _ _.,„ , qi. ''w ;„--1, -..----,-›;:',ts—„, ---:=... _.--4.1 iC20 C.; i::;:,' IT ' ^)0 - ,-,.,,.. ,','\ ; ^.,--‘,Ile-''''‘'‘,„7-_.41.."-7.—,,,•.=--_---;----,riiit ---- _ci.,-rati.i..-_:.:=.--. s---. ----_-_,:niki1Nir -- CL-1.-..". -,-Pai,, ( ,, ,,, t ,',.:'-‘‘-_, .`"-:-,•-----:--"t77.---,--Z"-"mtissiS .^, h,.. it',' 0pi:,)''',-r,',Y1„,<',..-i,,,,er;),:-----:‘,;It),i''.11.e• , it..- ii...,,,„...,.,-),-,_„).,.2.,,,---..,---_,...,, . ::-.". I s 110 120 10 °°•.i. *-.I .Aa (.73\3 -V \-- .,.,.,.... b-\Seismicity e . = \\ (SC: 1971-1983 rr-Vb

. .0*/. . 0 . III A;kitia .t.4% %. '' ...... ',i.'••:.'s• *7 °•••.''L .o.„ .7t*-''1.' 2)6) ...... ,...i...... 0 •46.. •••: '....:':'• • .g7 ° ,eZ1( --. .4%,* •#.1 t•• • . ; . . •k :Is, ...1 t. h •,:i.,.. *.J.: •---,. T . ..t.„..;7 7440°4.' ••% -it . .,Q,-.. .:L4;..i.v., -'. .* diiiii„.iPIC 10 . ..•: %.• . - ... : t • • . :'.;-.. 'Aphi, • . 47 . • 40 :I ;ti, It i,.,011; .. ..,:.

100 110 120 Fig. 7a. Free air gravity anomaly map of the Sunda arc region, as derived from satellite altimeter data (after Segawa and Matsumoto, 198731)).Solid and dashed contours represent negative and positive anomalies respectively. Arrows point to the locations (Fig. 7b) corresponding to high earthquake activity. Fig. 7b. Seismicity distribution map of the Sunda arc (ISC Data File, 1971-1983,with events located by 15 or more stations), for comparing with gravity anomaly. Arrows show the patches of high seismicity.

i 4 42 R. GHOSE and K. OIKE from GEOS-3 and SEASAT-1 altimeter data (Fig. 7a) . Segawa and Matsumoto (1987)3" converted analytically the satellite altimeter data to free air gravity anomaly for the oceanic regions. Although differences exist between this alimetric gravity anomaly and true gravity anomaly particularly on seamounts and trenches, for preliminary investigation the accuracy is quite acceptable. In Fig. 7, the contour interval is of 20 mgal. Solid and dashed contours are respectively for negative and positive anomalies. (Although we did not draw the 0 mgal contour for the sake of simplicity of the map, the smallest radius solid line contour represents the negative minima). For comparision, in Fig. 7b we presented a plot of the ISC earthquake epicentres (1971-1983) located by 15 or more stations. The arrows show the zones of high earthquake frequency. The corresponding localities are also indicated by arrows in the free air gravity anomaly map. One can see a good conformity between the seismicity and the altimetric free air gravity anomaly in that, most of the earth- quake clusters are located where the negative free air gravity anomaly contours are rather dense (i.e., deeper parts). Excepting a few (viz., d, g and h), for most of the seismicity clusters the matching is fairly good. Even if the possibilities of some earthquake mislocations and the deviation of the altimetric gravity anomaly from the true anomaly are taken into account, the picture does not appear to become very much different. Morphologically, the earthquakes tend to concentrate in the deeper regions of the inter-arc basin or its slope regions. We noticed, in vertical seismicity profiles taken parallel to the trench that, below these patches of earthquake clustering, the maximum depths of seismicity are greater compared to the surrounding. At the present moment, our evidences are not conclusive. However, the observation invites attention for future investigation.

Conclusions Based on our present investigation of the spatio-temporal variations of seismicity the following conclusions could be reached : (a) From the epicentral distribution maps, as well as the strain release map, we could notice that adjacent to the zones of high strain release, seismically quiet zones are often present. Some of these low strain release localities are possibly located within the rupture zones of the adjacent large events, and hence being in a stage of strain accumulation, are presently quiet. However, as the possibility of a few of these quiet zones to be potential seismic gaps certainly exists, more detailed and accurate spatio-temporal analysis is very necessary to ascertain their true natures. (b) A careful inspection of the depth distribution of seismicity revealed that in a part of the eastern Sunda arc, possibly a zone of scarce seismicity at an intermediate depth exists. We discussed the implications, and recommended further investigations. (c) The patterns of temporal variation of seismicity at the three different seismotectonic provinces of the Sunda arc—Sumatra, Java, and the Lesser Sunda Islands, were analysed. It could be observed that, although the causative geodynamic Characteristicsof Seismicity Distribution along the Sunda Arc: SomeNew Observations 43 situations for seismicity vary significantly in space along the length of the arc, the period of increase or decrease in seismicity is largely space invariant. The locally differing levels of seismicity are superposed on the common background of long period seismicity fluctuation. (d) Finally, a comparision of the locations of seismicity clustering at the Sunda arc with the contours of the altimetric free air gravity anomaly revealed some positive comformities. We recognised this observation as a preliminary one, for the purpose of inviting more devoted research attention in future.

Acknowledgements We are grateful to Dr. Y. Kishimoto, Mr. K. Watanabe, Mr. F. Takeuchi and Dr. K. Matsumura of the D.P.R.I., Kyoto University for their valuable comments and cooperations. Dr. J. Segawa of the Ocean Research Institute, University of Tokyo, kindly provided the altimetric gravity anomaly maps. R.G. acknowledges thankfully the scholarship awarded to him by the Ministry of Science, Education and Culture, Govt. of Japan.

References

1) Newcomb, K.R., and W.R. McGann: Seismic history and seismotectonics of the Sunda arc, J. Geophys. Res., Vol. 92, 1987, pp. 421-439. 2) Minster, J.B., and T. Jordan: Present-day plate motions, J. Geophys. Res., Vol. 83, 1978, pp. 5331-5354. 3) Curray, J.R., D.G. Moore, L.A. Lawyer, F.J. Emmel, R.W. Raitt, M. Henry, and R. Kieckhefer: Tectonics of the Andaman Sea and Burma, Am. Assoc. Pet. Geol. Mem., Vol. 29, 1978, pp. 189- 198. 4) Silver, E.A., D. Reed, R. McCaffrey, and Y. Joyodiwiryo: Back arc thrusting in the eastern Sunda Arc, : a consequence of arc-continent collision, J. Geophys. Res., Vol. 88, 1983, pp. 7429-7448. 5) Budiarto, R.: Sunda Strait: a dividing line between structural patterns in Sumatra and Java islands, Geol. Indonesia, Vol. 3, 1976, pp. 11-20. 6) Ranneft, T.S.M.: Segmentation of island arcs and application to Petroleum Geology, J. Petro- leum Geology, Vol. 1, 1979, pp. 35-53. 7) Zen, M.T.: Krakatau and the tectonic importance of Sunda Strait, Buletin Jurusan geologi, Vol. 12, 1983, pp. 9-22. 8) Curray, J.R., G.G. Shor, Jr., R.W. Raitt, and M. Henry: Seismic refraction and reflection studies of crustal structure of the eastern Sunda and western Banda arcs, J. Geophys. Res., Vol. 82, 1977, pp. 1479-1489. 9) Le Pichon, X., and J.R. Heirtzler: Magnetic anomalies in the Indian Ocean and sea-floor spreading , J. Geophys. Res., Vol. 73, 1968, pp. 2101-2118. 10) Fitch, T.J.: Earthquake mechanisms and island arc tectonics in the Indonesian-Philippine region, Bull. Seis. Soc. Am., VoI. 60, 1970, pp. 565-591. 11) Fitch, T.J., and P. Molner: Focal mechanisms along inclined earthquake zones in the Indonesia- Philippine region, J. Geophys. Res., Vol. 75, 1970, pp. 1431-1444. 12) Fitch, T.J.: Plate convergence, transcurrent faults, and internal deformation adjacent to south- east Asia and the western Pacific, J. Geophys. Res., Vol. 77, 1972, pp. 4432-4460. 13) Katili, J.A.: Sumatra, in - Orogenic Belts, Data for Orogenic Studies, Spec. 44 R. GHOSEand K. DIKE

Publ. 4, ed. A.M. Spencer,Geological Society, London, 1974, pp. 317-331. 14) Tija, H.D.: Activefaults in Indonesia,Geol. Soc. Malaysia, Vol. 10, 1978,pp. 73-92. 15) Huchon, P., and X. Le Pichon: Sunda Strait and Central Sumatra Fault, Geology,Vol. 12, 1984,pp. 668-672. 16) Perez, 0., and C. Scholz: Heterogeneitiesof the instrumentalseismicity catalogue (1904-1980) for strongshallow earthquakes, Bull. Seis. Soc. Am., Vol. 74, 1984,pp. 669-686. 17) McCaffrey,R., P. Molner, and S.W. Roecker: Microearthquakeseismicity and fault plane solutionsrelated to arc-continentcollision in the easternSunda arc, Indonesia,J. Geophys. Res., VoI.90, 1985,pp. 4511-4528. 18) McCaffrey,R., and J. Nabel6k: The geometryof back arc thrustingalong the easternSunda arc, Indonesia:constraints from earthquakeand gravitydata, J. Geophys.Res., Vol. 89, 1984, 6171-6179. 19) McCaffrey,R., and J. Nabelék: Seismologicalevidence for shallowthrusting north of the Timor Trough,Geophys. J.R. astr. Soc., Vol. 85, 1986,pp. 365-381. 20) Spence,W.: The 1977Sumba earthquakeseries: evidence for slab pull forceacting at a sub- ductionzone, J. Geophys.Res., Vol. 91, 1986,pp. 7225-7239. 21) Allen,C.R., P. St. Amand,C.F. Richter,and J.M. Nordquist: Relationshipbetween seismicity and geologicstructure in the southernCalofornia region, Bull. Seism. Soc. Am., Vol. 55, 1965, pp. 753-797. 22) Tocher, D.: Earthquakeenergy and ground breakage,Bull. Seis. Soc. Am., Vol. 48, 1958, pp. 147-153. 23) Iida,K.: Earthquakemagnitude, earthquake fault and sourcedimensions, J. Earth Sci.,Nagoya Univ., Vol. 13, 1965,pp. 115-132. 24) Ambraseys,N., and J. Tchalenko: Documentationof faulting associatedwith earthquakes, Part I. Unpublishedmanutcript, Dept. of Civil Engineering,Empirical College of Science, London,1968. 25) Bonilla,M.G., and J.M. Buchanan: Interimreport on worldwidehistoric surface faulting, U.S. Geol.Surv., Open File Report, 1970. 26) Slemmons,D.B.: Faults and earthquakemagnitude, Rept. 6, State of the Art for Assessing EarthquakeHazards in the United States, MiscellaneousPaper S-73-I, U.S. Army Engineer WaterwaysExperiment Station, Vicksburg, Miss., 1977, 129 pp. 27) Acharya, H.K.: Regionalvariations in the rupture-lengthmagnitude relationships and their dynamcalsignificance, Bull. Seism. Soc. Am., Vol. 69. 1979,pp. 2063-2084. 28) Habermann,R.E., W.R. McCann,and B.Perin: Spatialseismicity variations along convergent plate boundaries,Geophys. J.R. astr. Soc.,Vol. 5, 1986,pp. 43-68. 29) Ishikawa,Y., K. Matsumura,H. Yokoyama,and H. Matsumoto: SEIS-PC—itsoutline, Geol. Data Proc.,Vol. 10, 1985,pp. 19-34. 30) Isacks,B., and P. Molner: Mantleearthquake mechanisms and the sinkingof the lithosphere, Nature,Vol. 223, 1969,pp. 1121-1124. 31) Segawa,J., and T. Matsumoto: Free air gravityanomaly of the world oceanas derivedfrom satellitealtimeter data, Bull. OceanRes. Inst., Univ. of Tokyo,No. 25, 1987.

Appendix Strain Release Map of the Sunda Arc Region (presented in four parts A, B, C and D) Characteristicsof SeismicityDistribution along the Sunda Arc: SomeNew Observations 45

11) .- 0 ,-- 1.1.1 a -..... a

Cr) -0 .r-

\

isk 11. PAP'''-r\b*:)1/4C.414' Lg "r

o .- en en .- —en Nco coN ea,- cola -0) 0) V- 'll- l'-• N , V1 t... 10, ,-Net ,,,,46 - ,-01,-C-J.^)"r-r0l. .--1c„

i,... ".--P-,-..,-. °'Eli.." Z; 2 4inZ.;o..CI ,9. ;-C7~17_

8 NNN•-•N NN4CDsr c.0 • CI ..- .-- 9. / iv -,-io

co CO Pc 0N N KI el0 l•-•el CIb w sr co t•D Lc) r- ..----.7/ sren en .-

N o N NCI '''.N;1.3- N;7,(9, C. CO•-• in 0) V t•-•CO.3COLo.0 o3 N cn r-c. - NN. co chm N N NN a.- 0- N - - of co a -,- v - t:.ti ;Z• T a> =(i:CP,_.to N. to F. co (c. 4 CO . atr, N (.7 N 1 I 1,....0) izt.. ,..) r-

4 sr co sr 1.- 0 ...- 2 a ' co coo

( 46 R. CHOSE and K. OIKE

2 B 17 2421 9 10113 13 17 18 25 l,j (N) 9 2ts3219119 61)2.3 45\9\.) 9734 5 7 93 C' 9 8 25 3 11614 18 28 0 B 10 1198 16 126 7 7 — 6 8 8 22 32 133183 38 if 3 8 11 26 7 21 17 E 17 17 ..19 5111061ro. (S) 8 2 9 \32252 13\2 22'$) 2 33 177 15 150 5 106 2 1 93 27 8 21 6 28 23 \7 3 28 16 4 218 5 22 1341 10 14 38 32 7 40 37 35 95 31 3 4 21 10 20025 7 6 9 7 29 32\3 15 114 78 10717305514 11 7 15 25 6 61 214 48194\-,78 6 1 .. .,.. 27 245 99 156 5216,10 \ 9 5 5 152 6 47 137 56 13 38 .C\ 251 6 16 6 2 11 4 60 77 90 50 91,-8188 31 949498 1 43 2 15 31 13 99 12276 10792 27 6 38 62 69 7 98406-8 6 15 7 2 13 73 26 6 16 5 12 3 3 3 3 6 17 2

10 1 1 102 (°E) 106

_ Characteristicsof SeismicityDistribution along the Sunda Arc: SomeNew Observations 47

2 C C3 BORNEO3 C 32 rS) 431 63 $8122 32 166 10 8 3 7 32 62 92 4 149 91 172 63 75 6 11 7 19 15 96 113 13 251 4 9 35 50 32 6 16 102 3 27 10 9737 47 3 19 15 15 7 8 (-) P 6 11 9 7232

3 61-26_73"-- 2 30 a 13 2 JAVA 11 3 9 182 1897 --..,..---- 5 76 41 11 444683-510 43 6 68 3,5417 610''''-----'13 21 7 _3-22 5 17 7 1• 6 3 1370 35 1„ 49 39 30 13 8 6 25 9 8 67 68 1541 5 5 53 7 73i11923515-'.."3-5--\. 8 15 9 2 8 25 32 221 12 33 168438 7 11 cr....)3175 48 ..------352526 7 12 30 36 118912 28 5 35 21 15 12 9 23 11 2 31 141 103 167 34 11 3 7 9 22 11 36 4 6 69 16 5 2 19 35 8 5273 147 1471724 7 8 63 11 40 3 11 4 2 49 6 53 53 7 5 9 6 4 13 90 38 10 1 2 8 14 4 13 35 10732 47 65 4 10 6 7 27 8 9 7 INDIAN OCEAN 3

14 I 1 1I I 106 110114 (°E) 118 4 8 R. CHOSE and K. DIKE

a ....., CI a,et.1..1 ...... \ ....cl la v a. ..r10I-el0 —01 C.CY In ID 10 1- CII 0 CD 4- s- s- -N-_.4.IDCP,„0NNelID1.".0

CO g". 'P; ; ?I

SO Of0. ,-.0 0 100400 - O.1.".8 .r, N NTr.,_ IN 0 0 " el C. Cl c'. 10 •-• el ILI N 0 0...... ,.-N03...-•-'- .- Cl.C4 r- et r)ci0lo. o-et CD0 f,C.CD I cO. GI N10...*-1.... cvInes slo O" 0.er •Ng1N 00 .a.,1-, ,,-10 ,_ VI41,72I0 I-es r) 0 03 ...-^ C. ..- c^ICII a' CO CO

CO .1. 44.. , 0 N l•-• ,NDI. srcmro)N 0 .-,- INs-•—...„,,, w i„,in.., ._„.20—mrricv0.Coo,s,4-s...._z_r....._(Cl et roetro ..- cm<0o-r-.I.Oi0/.-z _r\ 0 .,.a20'. s „.._ ,.., c,co Cl)IAsN

C,,,-0, ,-GO,••, to tN —01 IN0etrto-.1.- m ,.-et N0 N117V0C. 1043 03 WI coel _ on 0, N Lad C., C., CO 10 N ,C.CO,.CO.4:.---.)et CC Ner....1-CO1-•^•... 4,, Q. CONCO - ea.NN ,-.. V3I, 0 t 4t, /..- VI 040s- ,- .- 0soon ej,,c-,or - 0 eg r- .3,-.Ns-i•••0s-.1. N^•••.4. = II .CC IN CI V ‘.. or-aonga .I. Cl co . \--) r-es.v. ae C.,N.;F.-. ,9. 2 ., .cc 0 ozyN COt,- ...,,.,, ,-r- ,...so ,.., CIsO '-'`,,,, bet9 so N..,CONS V10CONs-10 . 00 I- sO IN CI = ..... /r/ 03 I -- I II.- .- (0CA 0'I a 1, 1"...

. -