Joan L. Latchman

Soufrière Hills Eruption, Montserrat, 1995–1997: Volcanic earthquake locations and fault plane solutions

A total of 9242 seismic events, recorded since the start of the eruption on Montserrat in July 1995, have been uniformly relocated with station travel-time corrections. Early seismicity was generally diffuse under southern Montserrat, and mostly restricted to depths less than 7 km. However, a NE-SW alignment of epicentres beneath the NE flank of the volcano emerged in one swarm of volcano-tectonic earthquakes (VTs) and later nests of VT hypocentres developed beneath the volcano and at a separated location, under St. Georges Hill. The overall spatial distribution of hypocentres suggests a minimum depth of about 5 km for any substantial magma body. Activity associated with the opening of a conduit to the surface became increasingly shallow, with foci concentrated below the crater and, after dome building started in Fall 1995, VTs diminished and repetitive swarms of ‘hybrid’ seismic events became predominant. By late-1996, as magma effusion rates escalated, most seismic events were originating within a volume about 2 km diameter which extended up to the surface from only about 3 km depth - the diminution of shear failure earthquakes suggests the pathway for magma discharge had become effectively unconstricted. Individual and composite fault plane solutions have been determined for a few larger earthquakes. We postulate that localised extensional stress conditions near the linear VT activity, due to interaction with stresses in the overriding lithospheric plate, may encourage normal fault growth and promote sector weaknesses in the volcano.

Seismicity associated with dome growth and collapse at the Soufriere Hills Volcano, Montserrat

Varied seismicity has accompanied growth and collapse of the lava dome of the Soufriere Hills Volcano, Montserrat. Earthquakes have been classified as either volcano-tectonic, long-period or hybrid, and daily variations in the numbers of events have mapped changes in the style of eruption. Repetitive hybrid earthquakes were common during the first months of dome growth. In July 1996 the style of seismicity changed and regular, short-lived hybrid earthquake swarms became common. This change was probably caused by an increase in the magma flux. Earthquake swarms have preceded almost all major dome collapses, and have accompanied cyclical deformation, thought to be due to a build-up of pressure in the upper conduit which is later released by magma moving into the dome.

Spatial variations in the frequency‐magnitude distribution of earthquakes at Soufriere Hills Volcano, Montserrat, West Indies

The frequency-magnitude distribution of earthquakes measured by the b-value is determined as a function of space beneath Soufriere Hills Volcano, Montserrat, from data recorded between August 1, 1995 and March 31, 1996. A volume of anomalously high b-values (b > 3.0) with a 1.5 km radius is imaged at depths of 0 and 1.5 km beneath Englishs Crater and Chances Peak. This high b-value anomaly extends southwest to Gages Soufriere. At depths greater than 2.5 km volumes of comparatively low b-values (b∼1) are found beneath St. Georges Hill, Windy Hill, and below 2.5 km depth and to the south of Englishs Crater. We speculate the depth of high b-value anomalies under volcanoes may be a function of silica content, modified by some additional factors, with the most siliceous having these volumes that are highly fractured or contain high pore pressure at the shallowest depths.

Probabilistic Seismic Hazard Assessment at the Eastern Caribbean Islands

A probabilistic seismic hazard analysis has been performed to compute probabilistic seismic hazard maps for the eastern Caribbean region (10° N-19° N, 59° W-64° W), which includes in the north the Leeward Islands (from Anguilla to Dominica) and in the south the Windward Islands (from Martinique to Grenada), Barbados, Trinidad, and Tobago. The analysis has been conducted using a standard logic-tree approach that allowed systematically taking into account the model-based (i.e., epistemic) uncertainty and its influence on the computed ground motion param- eters. Hazard computations have been performed using a grid of sites with a space resolution of 0.025 degrees covering the territory of the considered islands. Two different computation methodologies have been adopted: the standard Cornell-McGuire approach (Cornell, 1968; McGuire, 1976) based on the definition of appropriate seismogenic zones (SZ), and the zone-free approach developed by Woo (1996), which overcomes the ambiguities related with the definition of seis- mic sources. The interplay and complexities between shallow crustal, intraplate, and interface subduction seismicity of the Caribbean region have been thoroughly investigated. By merging all available databases, a comprehensive and updated earthquake catalog for the region has been compiled. Also, a thorough investigation has been undertaken to identify the most suitable ground motion prediction equa- tions to be used in the analyses. Uniform hazard spectra have been calculated for the horizontal component of ground motion (rock and level site conditions), 4 return periods (RP) (95-, 475-, 975-, and 2475-yr), and 22 spectral accelerations (SA) with structural periods ranging from 0 to 3 s. SAs at 0.2 and 1.0 s for 2475-yr RP have been calculated to allow the definition of seismic hazard in the region of study according to the International Building Code (IBC, International Code Council (ICC), 2009).

Attenuation of seismic waves in the Trinidad and Tobago area

Abstract The attenuation of seismic waves from earthquakes located within the area bounded by 9–12°N and 60–63°W was estimated from short-period seismograms. Coda- Q , Q c , determinations were made for each of the six seismograph stations within the area, while spectral Q values from P-phases, Q α , were estimated for station TRN. The S-S single-scattering model was assumed for coda generation, and the ω −2 source model was assumed for the spectral Q determinations. The Q c values show a strong frequency dependence in the frequency range 1.5–12 Hz. The value of Q at 1 Hz, Q o , was found to lie within the range 107–132, while the rate of frequency dependence, n , extends from 0.80 to 1.06 for shallow events. For intermediate-depth events, Q o varies from 101 to 173 and n from 0.80 to 1.02. The Q α values obtained show a spatial variation within the region, the highest attenuation was obtained on land Trinidad.

Tectonic inversion in the Caribbean‐South American plate boundary: GPS geodesy, seismology, and tectonics of the Mw 6.7 22 April 1997 Tobago earthquake

On 22 April 1997 the largest earthquake recorded in the Trinidad-Tobago segment of the Caribbean-South American plate boundary zone (Mw 6.7) ruptured a shallow (~9 km), ENE striking (~250° azimuth), shallowly dipping (~28°) dextral-normal fault ~10 km south of Tobago. In this study, we describe this earthquake and related foreshock and aftershock seismicity, derive coseismic offsets using GPS data, and model the fault plane and magnitude of slip for this earthquake. Coseismic slip estimated at our episodic GPS sites indicates movement of Tobago 135 ± 6 to 68 ± 6 mm NNE and subsidence of 7 ± 9 to 0 mm. This earthquake was anomalous and is of interest because (1) its large component of normal slip and ENE strike are unexpected given the active E-W dextral shearing across the Caribbean-South American plate boundary zone, (2) it ruptured a normal fault plane with a low (~28°) dip angle, and (3) it reactivated and inverted the preexisting Tobago terrrane-South America ocean-continent (thrust) boundary that formed during early Tertiary oblique plate convergence.

A microprocessor-based system for digitizing seismic events from magnetic-tape recordings

Abstract The task of converting analogue recordings of seismic events into digital form for computer input can be rendered faster and more reliable than hand-tracing on a digitizing table by the use of a simple microprocessor-based automatic system. We describe the circuitry and software of a system which takes 4 sec to digitize an event and 40 min to punch a card-deck containing 4096 data points; hand-tracing the same record on a digitizing table takes 8 hr. In addition, the automatic system removes the needs of correcting for duplicate points and axis alignment, and of interpolating to obtain equispaced data points. The system, as described, has been configured at minimal cost for input to an outdated mainframe facility; it demonstrates that microprocessor technology can be used to produce time- and labour-saving gains.