Rupture process of the 19 August 1992 Susamyr, Kyrgyzstan, earthquake

The Susamyr earthquake of August 19, 1992 in Kyrgyzstan is one of the largest events (M_s = 7.4, M_b = 6.8) of this century in this region of Central Asia. We used broadband and long period digital data from IRIS and GEOSCOPE networks to investigate the source parameters, and their space-time distribution by modeling both body and surface waves. The seismic moment (M₀ = 6.8 × 10¹⁹ N m) and the focal mechanism were determined from frequency-time analysis (FTAN) of the fundamental mode of long period surface waves (100–250 s). Then, the second order integral moments of the moment-rate release were estimated from the amplitude spectra of intermediate period surface waves(40–70 s). From these moments we determined a source duration of 11–13 s, major and minor axes of the source of 30 km and 10–22 km, respectively; and an instant centroid velocity of 1.2 km/s. Finally, we performed a waveform inversion of P and SH waves at periods from 5–60 s. We found a source duration of 18–20 s, longer than the integral estimate from surface wave amplitudes. All the other focal parameters inverted from body waves are similar to those obtained by surface waves ( = 87° ± 6°, = 49° ± 6°, = 105° ± 3°, h = 14 ± 2 km, and M₀ = 5.8 ± 0.7 × 10¹⁹ N m). The initial rupture of this shallow earthquake was located at the south-west border of Susamyr depression in the western part of northern Tien Shan. A finite source analysis along the strike suggests a westward propagation of the rupture. The main shock of this event was preceded 2 s earlier by small foreshock. The main event was almost immediately followed by a very strong series of aftershocks. Our surface and body wave inversion results agree with the general seismotectonic features of the region.     

Source parameters of the Pinotepa Nacional, Mexico, earthquake of 27 March, 1996 (Mw = 5.4) estimated from near-field recordings of a single station

We use near-field accelerograms recorded by the very broadband seismographic station of PNIG to locate the Pinotepa Nacional earthquake of 27 March, 1996 (Mw = 5.4) and to determine its source parameters. The data from PNIG on P and S arrival times, the azimuth of the arrival of P wave, and the angle of incidence of the P wave at the free surface permit the determination of the location (16.365° N, 98.303° W, depth = 18 km) and the origin time (12:34;48.35) of the earthquake.The displacement seismograms of the earthquake clearly shows contribution from the near-field terms. We compute a suite of synthetic seismograms for local mechanisms in the vicinity of the mechanism reported by the U.S. Geological Survey (USGS) and compare them with the observed seismograms at PNIG. The point whose synthetics fit the observed records well has the following parameters: seismic moment, M0 = 1.2 × 10²⁴ dyne-cm; source time function: a triangular pulse of 0.9 sec duration; fault plane: strike = 291°, dip = 10°, and rake = 80°. The location and the source parameters obtained from the analysis of PNIG records differ significantly from those reported by the USGS. This demonstrates again, what has been shown by some previous researchers, that high-quality recordings from a single near-field station can considerably improve the estimation of the source parameters of an earthquake.The main event was preceded by a subevent which occurred 0.18 sec before and whose seismic moment was 1% of the main event. It is possible that even this subevent was preceded by a couple of smaller subevents. This earthquake supports the body of evidence showing that an earthquake begins with a sequence of smaller subevents, cascading in the occurrence of the main event.     

Estimates of coda Q attenuation in the João Câmara area (Northeastern Brazil)

S coda wave of seventy-four local earthquakes recorded in a network of ten seismic stations were used to calculate coda Q attenuation (Q_c) in the João Câmara area (northeastern Brazil). The estimates show Q_c as a strong function of frequency in the range from 6.0 to 20.0 Hz. We found out that Q_c in João Câmara has a functional form given by Q_c= Q₀ f, where Q₀= 151 ± 99 and = 0.98 ± 0.05. If the standard deviations are taken into account,we conclude that there are no relevant changes in both Q₀ and values from one station to another. The estimated Q₀ values at the different stations suggest that the Samambaia fault is a boundary between two different seismic attenuation zones. In one side of the fault (left), where stations were installed in Pre-Cambrian terrain and thick sedimentary layer, the seismic attenuation is stronger than in the other side (stations installed in thin sedimentary layer and limestone outcrop).The anomalous Q₀ values in the left side of the Samambaia fault can be explained due to the presence of a shallow conductive layer in the upper crust( 10 km), such as proposed by Padilha et al. (1992). According to our results, if there is a conductive layer in the area, it probably spreads over João Câmara city and surrounding regions.However, more detailed investigation either with seismic methods (seismic attenuation,3D tomography with P and/or S wave velocities) or with other geophysical methods is needed to interpret the observed differences in Q₀ values between the two sides of the Samambaia fault.     

Mining-induced seismicity in faulted geologic structures: An analysis of seismicity-induced slip potential

Relationships between the locations of mining-induced seismic events, local fault structure, and mine geometry were examined in a deep hard-rock mine in northern Idaho. Stopes experiencing rock bursts and other large seismic events were found to fall into two structural regimes: the Silver Vein, and the N48°W Trend, a steeply dipping plane of seismic activity that is subparallel to major local steeply dipping faults which bound blocky structures. The N48°W Trend also intersects a shaft that was seriously damaged when fault gouge was expelled into the opening during a 3-month period of high seismic energy release. Models of stress interaction are used to support the hypothesis that mining-induced deformation was mobilized along a 1.5 km length of the N48°W Trend. Specifically, numerical models are used to simulate rupture of seismic events and estimate induced changes in the quasi-static stress field. A Coulomb failure criterion is used with these results to estimate the spatial variation in potential for slip on planes parallel to local faulting. Increases in the potential for slip on fault planes subparallel to the N48°W Trend are consistent with activation of deformation along its 1.5 km length. For events with constant seismic moment, stress drop is shown to be far more important than source dimension in elevating slip potential along the observed plane of seismic activity     

Analysis of aeromagnetic anomalies over the central part of the Narmada-Son lineament

The Narmada-Son lineament (NSL) forms a major tectonic feature on the Indian subcontinent. The importance of this lineament lies in its evolution as well as its tectonic history. The lineament seems to have been active since Precambrian times. In order to understand the history of its evolution, it is necessary to know what igenous activity has been taking place along this lineament, and how the Deccan trap volcanics, which cover large areas along this lineament, have erupted.For the study of this problem an analysis of the aeromagnetic anomaly map lying between 76°15 to 77°30E and 21°45 to 22°50N has been carried out. Four different profiles (B ₁ B ₁,B ₂ B ₂,B ₃ B ₃ andB ₄ B ₄) have been drawn in N-S direction over this area and interpreted in terms of the intrusive bodies present within or below the surface of Deccan trap exposures. Inversion and forward modelling techniques have been adopted for interpretation purposes. An analysis of frequency spectra along the profiles has also been carried out to estimate the average depth of the different magnetic bodies. These results have been correlated with the available geological information. It has been found that most of the small wavelength anomalies are caused by dyke-like bodies within or below the Deccan trap at a depth of less than 0.5 km.     

Großkreisdifferenzgleichen

Summary Echo soundings of the U.S. Cruiser Milwaukee in the Puerto Rico Trough in 1939 are briefly discussed, and two depths of 30246 feet or 9219 m, found at 19° 36 N, 68° 20.5 W and at 19° 35N, 68° 8.75W, are stated to be the greatest depths which are known so far in the Atlantic Ocean.     

Coseismic slip in the 1964 Prince William Sound earthquake: A new geodetic inversion

The 1964 Prince William Sound earthquake (March 28, 1964;M _w =9.2) caused crustal deformation over an area of approximately 140,000 km² in south central Alaska. In this study geodetic and geologic measurements of this surface deformation were inverted for the slip distribution on the 1964 rupture surface. Previous seismologic, geologic, and geodetic studies of this region were used to constrain the geometry of the fault surface. In the Kodiak Island region, 28 rectangular planes (50 by 50 km each) oriented 218°N, with a dip varying from 8^o nearest the Aleutian trench to 9^o below Kodiak Island, define the rupture surface. In the Prince William Sound region 39 planes with variable dimensions (40 by 50 km near the trench, 64 by 50 km inland) and orientation (218°N in the west and 270°N in the east) were used to approximate the complex faulting. Prior information was introduced to constrain offshore dip-slip values, the strike-slip component, and slip variation between adjacent planes. Our results suggest a variable dip-slip component with local slip maximums occurring near Montague Island (up to 30 m), further to the east near Kayak Island (up to 14 m), and trenchward of the northeast segment of Kodiak Island (up to 17m). A single fault plane dipping 30°NW, corresponding to the Patton Bay fault, with a slip value of 8 m modeled the localized but large uplift on Montague Island. The moment calculated on the basis of our geodetically derived slip model of 5.0×10²⁹ dyne cm is 30% less than the seismic moment of 7.5×10²⁹ dyne cm calculated from long-period surface waves (Kanamori, 1970) but is close to the seismic moment of 5.9×10²⁹ dyne cm obtained byKikuchi andFukao (1987).     

Source characteristics of the 1992 Nicaragua tsunami earthquake inferred from teleseismic body waves

We analyzed the broadband body waves of the 1992 Nicaragua earthquake to determine the nature of rupture. The rupture propagation was represented by the distribution of point sources with moment-rate functions at 9 grid points with uniform spacing of 20 km along the fault strike. The moment-rate functions were then parameterized, and the parameters were determined with the least squares method with some constraints. The centroid times of the individual moment-rate functions indicate slow and smooth rupture propagation at a velocity of 1.5 km/s toward NW and 1.0 km/s toward SE. Including a small initial break which precedes the main rupture by about 10 s, we obtained a total source duration of 110 s. The total seismic moment isM _o =3.4×10²⁰ Nm, which is consistent with the value determined from long-period surface waves,M _o =3.7×10²⁰ Nm. The average rise time of dislocation is determined to be 10 s. The major moment release occurred along a fault length of 160 km. With the assumption of a fault widthW=50 km, we obtained the dislocationD=1.3 m. From andD the dislocation velocity isD=D/0.1 m/s, significantly smaller than the typical value for ordinary earthquakes. The stress drop =1.1 MPa is also less than the typical value for subduction zone earthquakes by a factor of 2–3. On the other hand, the apparent stress defined by 2E _s /M _o , where andE _s are respectively the rigidity and the seismic wave energy, is 0.037 MPa, more than an order of magnitude smaller than . The Nicaragua tsunami earthquake is characterized by the following three properties: 1) slow rupture propagation; 2) smooth rupture; 3) slow dislocation motion.     

Comparison of a complex rupture model with the precursor asperities of the 1975 HawaiiM s-7.2 earthquake

A simplified multiple source model was constructed for the 1975 HawaiiM _s=7.2 earthquake by matching synthetic signals with three component accelerograms at two stations located approximately 45 km from the epicenter. Six major subevents were identified and located approximately. The signals of these are larger by factors of 1.4 to 3.2 than that of theM _L=5.9 foreshock which occurred 70 minutes before the main rupture and also triggered the SAM-1 recorders at the two stations. Dividing the rupture length (40 km) by the duration of strong ground shaking ( 50 sec) an, average rupture velocity of 0.8 km/sec (about 25% of S-velocity) is obtained. Thus it is likely that the rupture stopped between subevents. The approximate epicenters of the 6 major subevents, and of the foreshock, support the hypothesis that they were located in high stress asperities which rupture during the main shock, except for the last events which is interpreted as a stopping phase generated at a barrier. These asperities have been previously defined on the basis of differences in the precursor pattern before the mainshock. Thus, it appears that both the details of the precursors and of the main rupture depended critically on the heterogeneous tress distribution in the source volume. This suggests that main rupture initiation points and locations of high rupture accelerations may be identified before the mainshock occurs, based on precursor anomaly patterns. A satisfactory match of synthetic signals with the observations could be obtained only if the aximuth of the fault plane of subevents was rotated from N60°E to N90°E and back to N30°E. These orientations are approximately parallel to the nearest Kilauea rift segments. Hence the slip directions and greatest principal stresses were oriented perpendicular to the rifts everywhere. From this analysis and other work, it is concluded that this fault surface consisted of three types of segments with different strength: hard asperities (radius 5 km), soft but brittle segments between the asperities (radius 5 km), and a viscous half (10×40 km) which slipped during the mainshock, but where microearthquakes and aftershocks are not common.     

The hydrothermal system of Volcan Puracé, Colombia

This paper presents chemical and isotopic data for thermal waters, gases and S deposits from Volcan Puracé (summit elevation 4600 m) in SW Colombia. Hot gas discharges from fumaroles in and around the summit crater, and thermal waters discharge from three areas on its flanks. The waters from all areas have D values of-75±1, indicating a single recharge area at high elevation on the volcano. Aircorrected values of³He/⁴He in thermal waters range from 3.8 to 6.7 R_A, and approach those for crater fumarole gas (6.1–7.1 R_A), indicating widespread addition of magmatic volatiles. An economic S deposit (El Vinagre) is being mined in the Rio Vinagre fault zone at 3600 m elevation. Sulfur isotopic data are consistent with a magmatic origin for S species in thermal waters and gases, and for the S ore deposit. Isotopic equilibration between S species may have occurred at 220±40°C, which overlaps possible equilibration temperatures (170±40°C) determined by a variety of other geothermometers for neutral thermal waters. Apparent CH₄–CO₂ equilibration temperatures for gases from thermal springs (400±50°C) and crater fumaroles (520±60°C) reflect higher temperatures deeper in the system. Hot magmatic gas ascending through the Rio Vinagre fault zone is though to have precipitated S and generated thermal waters by interaction with descending meteoric waters.     

Fluid infiltration into fault zones: Chemical,isotopic, and mechanical effects

Fluid infiltration into fault zones and their deeper-level counterparts, brittle-ductile shear zones, is examined in diverse tectonic environments. In the 2.7 Ga Abitibi greenstone belt, major tectonic discontinuities, with lateral extents of hundreds of kilometres initiated as listric normal faults accommodating rift extension and acted as sites for komatiite extrusion and locally intense metasomatism. During reverse motion on the structures, accommodating shortening of the belt, these transcrustal faults were utilised as a conduit for the ascent of trondhjemitic magmas from the base of the crust and of alkaline magmas from the asthenosphere and for the discharge of thousands of cubic kilometres of hydrothermal fluids. Such fluids were characterised by ¹⁸O=+6±2, D=–50±20, ¹³C=–4±4, and temperatures of 270 to 450°C, probably derived from devolatilisation of crustal rocks undergoing prograde metamorphism. Hydrothermal fluids were more radiogenic (⁸⁷Sr/⁸⁶Sr=0.7010 to 0.7040) and possessed higher than did contemporaneous mantle, komatiites or tholeiites, and thus carried a contribution from older sialic basement. A provinciality of⁸⁷Sr/⁸⁶Sr and ¹³C is evident, signifying that fault plumbing sampled lower crust which was heterogeneous at the scale of tens of kilometres. Mineralised faults possess enrichments of large ion lithophile (LIL), LIL elements, including K, Rb, Ba, Cs, B, and CO₂, and rare elements, such as Au, Ag, As, Sb, Se, Te, Bi, and W. Fluids were characterised by X_{CO 2}0.1, neutral to slightly acidic pH, low salinity 3 wt-%, K/Na=0.1, they carried minor CH₄, CO, and N₂, and they underwent transient effervescence of CO₂ during decompression. Clastic sediments occupy graben developed at fault flexures. The⁴⁰Ar/³⁹Ar release spectra indicate that fault rocks experienced episodic disturbance on time scales of hundreds of millions of years.At the Grenville front, translation was accommodated along two mylonite zones and an intervening boundary fault. The high-temperature (580°C) and low-temperature (430 to 490°C) mylonite zones, formed in the presence of deep-level crust-equilibrated fluids of metamorphic origin. Late brittle faults contain quartz veins precipitated from fluids with extemely negative ¹⁸O (–14 per mil) at 200 to 300°C. The water may have been derived from downward penetration into fault zones of precipitation of low¹⁸O on a mountain range induced by continental collision, with uplift accommodated at deep levels by the mylonite zones coupled with rebound on the boundary faults.Archean gneisses overlie Proterozoic sediments along thrust surfaces at Lagoa Real, Brazil; the gneisses are transected by brittle-ductile shear zones locally occupied by uranium deposits. Following deformation at 500 to 540°C, in the presence of metamorphic fluids and under conditions of low water-to-rock ratio, shear zones underwent local intense oxidation and desilication. All minerals undergo a shift of –10 per mil, indicating discharge of meteoric-water-recharged formation brines in the underlying Proterozoic sediments up through the Archean gneisses, during overthrusting; 1000 km³ of solutions passed through these structures. The shear zones and Proterozoic sediments are less radiogenic (⁸⁷Sr/⁸⁶Sr=0.720) than contemporaneous Archean gneisses (0.900), corroborating the transport of fluids and solutes through the structure from a large external reservoir.Major crustal detachment faults of Tertiary age in the Picacho Cordilleran metamorphic core complex of Arizona show an upward transition from undeformed granitic basement through mylonitic to brecciated and hydrothermally altered counterparts. The highest tectonic levels are allochthonous, oxidatively altered Miocene volcanics. This transition is accompanied by an increase of 12 per mil in ¹⁸O, from +7 to +19, and a 400°C decrease in temperature. Lower tectonic levels acted as aquifers for the expulsion of large volumes of higher-temperature reduced metamorphic fluids and/or evolved formation brines. The Miocene allochthon was influenced by a lower-temperature reservoir inducing oxidative potassic alteration; mixing occurred between cool downward-penetrating thermal waters and the hot, deeper aqueous reservoir.In general, flow regimes in these fault and shear zones follow a sequence, from conditions of high temperature and pressure with locally derived fluids at low water-to-rock ratios, during initiation of the structures, to high fluxes of reduced formation or metamorphic fluids along conduits as the structures propagate and intersect hydrothermal reservoirs. Later in the tectonic evolution and at shallower crustal levels there was incursion of oxidising fluids from near-surface reservoirs into the faults. In general, magmatism, tectonics, and fluid motion are intimately related.     

The Ms = 6.2, June 15, 1995 Aigion earthquake (Greece): evidence for low angle normal faulting in the Corinth rift

We present the results of a multidisciplinary study of the M_s = 6.2, 1995, June 15, Aigion earthquake (Gulf of Corinth, Greece). In order to constrain the rupture geometry, we used all available data from seismology (local, regional and teleseismic records of the mainshock and of aftershocks), geodesy (GPS and SAR interferometry), and tectonics. Part of these data were obtained during a postseismic field study consisting of the surveying of 24 GPS points, the temporary installation of 20 digital seismometers, and a detailed field investigation for surface fault break. The Aigion fault was the only fault onland which showed detectable breaks (< 4 cm). We relocated the mainshock hypocenter at 10 km in depth, 38 ° 21.7 N, 22 ° 12.0 E, about 15 km NNE to the damaged city of Aigion. The modeling of teleseismic P and SH waves provides a seismic moment M_o = 3.4 10¹⁸ N.m, a well constrained focal mechanism (strike 277 °, dip 33 °, rake – 77°), at a centroidal depth of 7.2 km, consistent with the NEIC and the revised Harvard determinations. It thus involved almost pure normal faulting in agreement with the tectonics of the Gulf. The horizontal GPS displacements corrected for the opening of the gulf (1.5 cm/year) show a well-resolved 7 cm northward motion above the hypocenter, which eliminates the possibility of a steep, south-dipping fault plane. Fitting the S-wave polarization at SERG, 10 km from the epicenter, with a 33° northward dipping plane implies a hypocentral depth greater than 10 km. The north dipping fault plane provides a poor fit to the GPS data at the southern points when a homogeneous elastic half-space is considered: the best fit geodetic model is obtained for a fault shallower by 2 km, assuming the same dip. We show with a two-dimensional model that this depth difference is probably due to the distorting effect of the shallow, low-rigidity sediments of the gulf and of its edges. The best-fit fault model, with dimensions 9 km E–W and 15 km along dip, and a 0.87 m uniform slip, fits InSAR data covering the time of the earthquake. The fault is located about 10 km east-northeast to the Aigion fault, whose surface breaks thus appears as secondary features. The rupture lasted 4 to 5 s, propagating southward and upward on a fault probably outcropping offshore, near the southern edge of the gulf. In the shallowest 4 km, the slip – if any – has not exceeded about 30 cm. This geometry implies a large directivity effect in Aigion, in agreement with the accelerogram aig which shows a short duration (2 s) and a large amplitude (0.5 g) of the direct S acceleration. This unusual low-angle normal faulting may have been favoured by a low-friction, high pore pressure fault zone, or by a rotation of the stress directions due to the possible dip towards the south of the brittle-ductile transition zone. This fault cannot be responsible for the long term topography of the rift, which is controlled by larger normal faults with larger dip angles, implying either a seldom, or a more recently started activity of such low angle faults in the central part of the rift.     

The 1986 Kermadec earthquake and its relation to plate segmentation

To evaluate the tectonic significance of the October 20, 1986 Kermadec earthquake (M _w =7.7), we performed a comprehensive analysis of source parameters using surface waves, body waves, and relocated aftershocks. Amplitude and phase spectra from up to 93 Rayleigh waves were inverted for centroid time, depth, and moment tensor in a two-step algorithm. In some of the inversions, the time function was parameterized to include information from the body-wave time function. The resulting source parameters were stable with respect to variations in the velocity and attenuation models assumed, the parameterization of the time function, and the set of Rayleigh waves included. The surface wave focal mechanism derived (=275°, =61°, =156°) is an oblique-compressional mechanism that is not easy to interpret in terms of subduction tectonics. A seismic moment of 4.5×10²⁰ N-m, a centroid depth of 45±5 km, and a centroid time of 13±3 s were obtained. Directivity was not resolvable from the surface waves. The short source duration is in significant contrast to many large earthquakes.We performed a simultaneous inversion ofP andSH body waves for focal mechanism and time function. The focal mechanism agreed roughly with the surface wave mechanism. Multiple focal mechanisms remain a possibility, but could not be resolved. The body waves indicate a short duration of slip (15 to 20 s), with secondary moment release 60s later. Seismically radiated energy was computed from the body-wave source spectrum. The stress drop computed from the seismic energy is about 30 bars. Sixty aftershocks that occurred within three months of the mainshock were relocated using the method of Joint Hypocentral Determination (JHD). Most of the aftershocks have underthrusting focal mechanisms and appear to represent triggered slip on the main thrust interface. The depth, relatively high stress drop, short duration of slip, and paucity of true aftershocks are consistent with intraplate faulting within the downgoing plate. Although it is not clear on which nodal plane slip occurred, several factors favor the roughly E-W trending plane. The event occurred near a major segmentation in the downgoing plate at depth, near a bend in the trench, and near a right-lateral offset of the volcanic are by 80 km along an E-W direction. Also, all events in the region from 1977 to 1991 with CMT focal mechanisms similar to that of the Mainshock occurred near the mainshock epicenter, rather than forming an elongate zone parallel to the trench as did the aftershock activity. We interpret this event as part of the process of segmentation or tearing of the subducting slab. This segmentation appears to be related to the subduction of the Louisville Ridge, which may act as an obstacle to subduction through its buoyancy.     

Source process of the southeast Beni-Suef, northern Egypt earthquake Using Empirical Green's Function Technique

The rupture process of a moderateearthquake (ML4.9) occurred southeastBeni-Suef, northern Egypt was analyzed. Thecharacterization of the rupture process andsource properties were extracted fromstudying the relative moment rate function(RMRF). RMRFs were retrieved bydeconvolution of small aftershock recordsusing an inversion procedure for onlyP-wave part of the records.Although, this event is a moderateearthquake, its RMRFs exhibitedcomplexity and directivity of the rupturebehavior. The deconvolution pulses reflectthree subevents or more. The detailedanalysis of the distinct subevents revealedgross temporal and spatial characteristicsof the rupture propagation. The azimuthalvariation in the time delay of thesubevents with respect to the initiationrupture indicates that the subevents arelocated at 0.85 ± 0.17 and3.5 ± 0.07 km in directions of320^° ± 45^° and330^° ± 15^° with rupturevelocities 3.4 ± 0.45 and4.0 ± 0.7 km/sec, respectively. Thismeans that the rupture is predominatelypropagated toward the North. Estimation ofthe rupture direction was combined withP-wave focal mechanism to identify thefault plane for the initial rupture ofmainshock.Source parameters were calculated for eachdistinct subevent, including seismicmoments 8.53E14 to 6.80E15 Nm, fault radii713 to 1800 m and stress drops 0.725 to2.932 MPa. The large estimated stress dropfor the main subevent reflects failureasperity.     

Evaluation of gases,condensates, and SO2 emissions from Augustine volcano,Alaska: the degassing of a Cl-rich volcanic system

After the March–April 1986 explosive eruption a comprehensive gas study at Augustine was undertaken in the summers of 1986 and 1987. Airborne COSPEC measurements indicate that passive SO₂ emission rates declined exponentially during this period from 380±45 metric tons/day (T/D) on 7/24/86 to 27±6 T/D on 8/24/87. These data are consistent with the hypothesis that the Augustine magma reservoir has become more degassed as volcanic activity decreased after the spring 1986 eruption. Gas samples collected in 1987 from an 870°C fumarole on the andesitic lava dome show various degrees of disequilibrium due to oxidation of reduced gas species and condensation (and loss) of H₂O in the intake tube of the sampling apparatus. Thermochemical restoration of the data permits removal of these effects to infer an equilibrium composition of the gases. Although not conclusive, this restoration is consistent with the idea that the gases were in equilibrium at 870°C with an oxygen fugacity near the Ni–NiO buffer. These restored gas compositions show that, relative to other convergent plate volcanoes, the Augustine gases are very HCl rich (5.3–6.0 mol% HCl), S rich (7.1 mol% total S), and H₂O poor (83.9–84.8 mol% H₂O). Values of D and ¹⁸O suggest that the H₂O in the dome gases is a mixture of primary magmatic water (PMW) and local seawater. Part of the Cl in the Augustine volcanic gases probably comes from this shallow seawater source. Additional Cl may come from subducted oceanic crust because data by Johnston (1978) show that Cl-rich glass inclusions in olivine crystals contain hornblende, which is evidence for a deep source (>25km) for part of the Cl. Gas samples collected in 1986 from 390°–642°C fumaroles on a ramp surrounding the inner summit crater have been oxidized so severely that restoration to an equilibrium composition is not possible. H and O isotope data suggest that these gases are variable mixtures of seawater, FMW, and meteoric steam. These samples are much more H₂O-rich (92%–97% H₂O) than the dome gases, possibly due to a larger meteoric steam component. The 1986 samples also have higher Cl/S, S/C, and F/Cl ratios, which imply that the magmatic component in these gases is from the more degassed 1976 magma. Thus, the 1987 samples from the lava dome are better indicators than the 1986 samples of degassing within the Augustine magma reservoir, even though they were collected a year later and contain a significant seawater component. Future gas studies at Augustine should emphasize fumaroles on active lava domes. Condensates collected from the same lava-dome fumarole have enrichments ot 10⁷–10² in Cl, Br, F, B, Cd, As, S, Bi, Pb, Sb, Mo, Zn, Cu, K, Li, Na, Si, and Ni. Lower-temperature (200°–650°C) fumaroles around the volcano are generally less enriched in highly volatile elements. However, these lower-termperature fumaroles have higher concentration of rock-forming elements, probably derived from the wall rock.     

Large earthquake doublets and fault plane heterogeneity in the northern Solomon Islands subduction zone

In the Solomon Islands and New Britain subduction zones, the largest earthquakes commonly occur as pairs with small separation in time, space and magnitude. This doublet behavior has been attributed to a pattern of fault plane heterogeneity consisting of closely spaced asperities such that the failure of one asperity triggers slip in adjacent asperities. We analyzed body waves of the January 31, 1974,M _w =7.3, February 1, 1974,M _w =7.4, July 20, 1975 (1437)M _w =7.6 and July 20, 1975 (1945),M _w =7.3 doublet events using an iterative, multiple station inversion technique to determine the spatio-temporal distribution of seismic moment release associated with these events. Although the 1974 doublet has smaller body wave moments than the 1975 events, their source histories are more complicated, lasting over 40 seconds and consisting of several subevents located near the epicentral regions. The second 1975 event is well modeled by a simple point source initiating at a depth of 15 km and rupturing an approximate 20 km region about the epicenter. The source history of the first 1975 event reveals a westerly propagating rupture, extending about 50 km from its hypocenter at a depth of 25 km. The asperities of the 1975 events are of comparable size and do not overlap one another, consistent with the asperity triggering hypothesis. The relatively large source areas and small seismic moments of the 1974 doublet events indicate failure of weaker portions of the fault plane in their epicentral regions. Variations in the roughness of the bathymetry of the subducting plate, accompanying subduction of the Woodlark Rise, may be responsible for changes in the mechanical properties of the plate interface.To understand how variations in fault plane coupling and strength affect the interplate seismicity pattern, we relocated 85 underthrusting earthquakes in the northern Solomon Islands Are since 1964. Relatively few smaller magnitude underthrusting events overlap the Solomon Islands doublet asperity regions, where fault coupling and strength are inferred to be the greatest. However, these asperity regions have been the sites of several previous earthquakes withM _s 7.0. The source regions of the 1974 doublet events, which we infer to be mechanically weak, contain many smaller magnitude events but have not generated any otherM _s 7.0 earthquakes in the historic past. The central portion of the northern Solomon Islands Arc between the two largest doublet events in 1971 (studied in detail bySchwartz et al., 1989a) and 1975 contains the greatest number of smaller magnitude underthrusting earthquakes. The location of this small region sandwiched between two strongly coupled portions of the plate interface suggest that it may be the site of the next large northern Solomon Islands earthquake. However, this region has experienced no known earthquakes withM _s 7.0 and may represent a relatively aseismic portion of the subduction zone.     

Determination of source duration and moment of earthquakes in the Indian Peninsula using love waveforms

Love waves recorded by a long-period seismograph at New Delhi (NDI) from seven earthquakes of magnitude 4.3 to 5.2 in Koyna and Bhatsa on the western coast and one earthquake in Ongole on the eastern coast of the Indian Peninsula have been used to determine the seismic moment for each of the earthquakes by waveform modeling. Transverse component of the synthetic seismogram shows that the maximum amplitude of waveform decreases with an increase of source duration. Thus for an evaluation of the seismic moment by equating the amplitude level of the observed and synthetic waveforms, we must know the source duration. The synthetic seismogram also indicates that a short source duration gives rise to a small but sharp pulse and this pulse is interpreted as anLg wave. Comparison of the observed and synthetic waveforms has been used for a simultaneous evaluation of the source duration and seismic moment. The source durations are found to vary between 2.2 and 4.4 s; for earthquakes with a magnitude range between 4.3 and 5.2 these durations are slightly higher than normal. We obtain moment (M ₀) of Ongole earthquake (M _L =5.1)as 1.7×10²⁴ dyne-cm; moments of Koyna and Bhatsa earthquakes (4.3M _L 5.2) on the western coast lie between 0.7×10²³ and 3.6×10²³ dyne-cm. Moment (M ₀)-magnitude (M _L ) relation logM ₀=1.5M _L +16.0 for the western United States region agrees as well, in general, with the results for the earthquakes in the Indian Peninsula.     

Analysis of State of Stress During The 1997 Earthquake Swarm In Western Bohemia

Focal mechanisms of 70 events of the January 1997 earthquake swarm were analysed using the Gephart and Forsyth method to determine the state of stress in the West Bohemia/Vogtland region and to reveal possible stress variations during the swarm activity. The method was applied to the whole set of 70 focal mechanisms and to several subsets distinguishing magnitudes of events and the spatial and temporal distribution of the swarm. The three representative stress models A ( ₁ ,A: 0/45 (az/dip), ₂ A: 134/35, ₃ A: 243/25), B ( ₁ ,B: 162/38, ₂ B: 263/14, ₃ B: 10/49), and C ( ₁ C: 135/30, ₂ C: 1/50, ₃ C: 239/23) which could characterise the state of stress in the region were found. Model C can be considered to be the most probable stress model because of its consistency with the European stress field. The results of the stress analysis applied to the individual subsets of focal mechanisms indicate that the state of stress was more uniform during the first phase of the swarm. To distinguish between the fault plane and the auxiliary plane of fault plane solutions a statistical approach was used. The results showed that the fault planes striking NW-SE and dipping 45° to NE were activated during the swarm.     

A relationship between phonon conductivity and seismic parameter Φ for silicate minerals

Summary The relationship between the phonon conductivity at room temperature (K _N ) and the seismic parameter () for silicate minerals is suggested. The considerations are based on the Debye model of thermal energy transport phenomena in solids and on the seismic equation of state for silicates and oxides given byAnderson (1967). The semiempirical relationship is the formK _N = 0.43^0.82 where is in km²/s² andK _N in mcal/cm s K, and the empirical relationship isK _N =(0.528±0.006) –(8.18±2.11). The laboratory data on thermal and elastic properties for several silicates were taken fromHorai andSimmons (1970).     

Source Parameters of the ML 3.8 Earthquake on January 20, 2000 near Meckenheim, Germany

In the early morning hours on January 20,2000 at 03:03:17 (UTC) an M_L = 3.8earthquake occurred east of the city ofMeckenheim. It awakened many people in theMeckenheim-Bonn-Siegburg area. At least 200people called the police or fire brigade inBonn. The earthquake was felt as far asLimburg (Westerwald) and east Belgium. Themaximum intensity (EMS) in theMeckenheim-Adendorf area was 5. Theinstrumental epicenter was located at50.60^° N latitude and 7.08^° Elongitude. Source depth was 9.5 km. Localmagnitude M_L was 3.8 ± 0.3 andseismic moment was 1.86*10¹⁴ Nm,corresponding to a momentmagnitude M_W of 3.4± 0.2.Following the Brune model, the sourceradius was 0.45 km with an averagedislocation of 1.0 cm and a 0.8 MPa stressdrop. The source mechanism from a firstmotion fault plane solution shows a normalfault mechanism with a significant strikeslip contribution. Trend and plunge of themajor and minor principal axes P/T are168^°/54^° and275^°/12^°, respectively. Aninversion of the moment tensor fromamplitudes of direct waves reveals a doublecouple mechanism with a smaller strike slipcontribution than the first motion faultplane solution, but also shows asignificant 16% non-double couplecomponent of the total moment. Theintensity epicenter was determined by agrid search algorithm and was found to be14 km east from the instrumental epicenter.The macroseismic magnitude was 3.6.