Search

An eruption of the Soufrière Hills Volcano (SHV) on the eastern Caribbean island of Montserrat began on 18 July 1995 and continued until February 2010. Within nine days of the eruption onset, an existing four-station analog seismic network (ASN) was expanded to 10 sites. Telemetered data from this network were recorded, processed, and archived locally using a system developed by scientists from the U.S. Geological Survey (USGS) Volcano Disaster Assistance Program (VDAP). In October 1996, a digital seismic network (DSN) was deployed with the ability to capture larger amplitude signals across a broader frequency range. These two networks operated in parallel until December 2004, with separate telemetry and acquisition systems (analysis systems were merged in March 2001). Although the DSN provided better quality data for research, the ASN featured superior real-time monitoring tools and captured valuable data including the only seismic data from the first 15 months of the eruption. These successes of the ASN have been rather overlooked. This article documents the evolution of the ASN, the VDAP system, the original data captured, and the recovery and conversion of more than 230,000 seismic events from legacy SUDS, Hypo71, and Seislog formats into Seisan database with waveform data in miniSEED format. No digital catalog existed for these events, but students at the University of South Florida have classified two-thirds of the 40,000 events that were captured between July 1995 and October 1996. Locations and magnitudes were recovered for ∼10,000 of these events. Real-time seismic amplitude measurement, seismic spectral amplitude measurement, and tiltmeter data were also captured. The result is that the ASN seismic dataset is now more discoverable, accessible, and reusable, in accordance with FAIR data principles. These efforts could catalyze new research on the 1995–2010 SHV eruption. Furthermore, many observatories have data in these same legacy data formats and might benefit from procedures and codes documented here.

SUMMARYEstimating subsurface density is important for imaging various geologic structures such as volcanic edifices, reservoirs and aquifers. Muon tomography has recently been used to complement traditional gravity measurements as a powerful method for probing shallow subsurface density structure beneath volcanoes. Gravity and muon data have markedly different spatial sensitivities and, as a result, the combination is useful for imaging structures on spatial scales that are larger than the area encompassed by crossing muon trajectories. Here we explore and test a joint inversion of gravity and muon data in a study area where there is an independently characterized target anomaly: a regionally extensive, high-density layer beneath Los Alamos, New Mexico, USA. We resolve the nearly flat-lying structure using a unique experimental set-up wherein surface and subsurface gravity and muon measurements are obtained above and below the target volume. Our results show that with minimal geologic (prior) constraints, the joint inversion correctly recovers salient features of the expected density structure. The results of our study illustrate the potential of combining surface and subsurface (e.g. borehole) gravity and muon measurements to invert for shallow geologic structures.

Imaging shallow subsurface density structure is an important goal in a variety of applications, from hydrogeology to seismic and volcanic hazard assessment. We assess the effectiveness of surface and subsurface gravity measurements in estimating the density structure of a well-characterized rock volume: the mesa (a small, flat-topped plateau) upon which the town of Los Alamos, New Mexico, USA is located. Our gravity measurements were made on the mesa surface above a horizontal tunnel and underground, within the tunnel. We demonstrate that, in the absence of other geophysical data such as seismic data or muon attenuation, subsurface (tunnel) gravity measurements are critical to accurately recovering geologic structure. Without the tunnel data, our resolution is limited to roughly the surface gravity station spacing, but by including the tunnel data we can resolve structure to a depth of ~ 10 times the surface gravity station spacing. Densities were obtained using both forward modeling and a Bayesian inverse modeling approach, incorporating relevant constraints from geologic observations. We find that Bayesian inversion, with geologically relevant prior, is a superior approach to the forward models in terms of both robustness and efficiency and correctly predicts the orientation and elevation of important geologic features.

Muons can be used to image the density of materials through which they pass, including geological structures. Subsurface applications of the technology include tracking fluid migration during injection or production, with increasing concern regarding such timely issues as induced seismicity or chemical leakage into aquifers. Current density monitoring options include gravimetric data collection and active or passive seismic surveys. One alternative, or complement, to these methods is the development of a muon detector that is sufficiently compact and robust for deployment in a borehole. Such a muon detector can enable imaging of density structure to monitor small changes in density – a proxy for fluid migration – at depths up to 1500 m. Such a detector has been developed, and Monte Carlo modeling methods applied to simulate the anticipated detector response. Testing and measurements using a prototype detector in the laboratory and shallow underground laboratory demonstrated robust response. A satisfactory comparison with a large drift tube-based muon detector is also presented.

We present an underground cosmic ray muon tomographic experiment imaging 3D density of overburden, part of a joint study with differential gravity. Muon data were acquired at four locations within a tunnel beneath Los Alamos, New Mexico, and used in a 3D tomographic inversion to recover the spatial variation in the overlying rock-air interface, and compared with a priori knowledge of the topography. Densities obtained exhibit good agreement with preliminary results of the gravity modeling, which will be presented elsewhere, and are compatible with values reported in the literature. The modeled rock-air interface matches that obtained from LIDAR within 4 m, our resolution, over much of the model volume. This experiment demonstrates the power of cosmic ray muons to image shallow geological targets using underground detectors, whose development as borehole devices will be an important new direction of passive geophysical imaging., 17 pages, 8 figures

Between May and June 2006, an earthquake swarm occurred near Volcan Tancitaro in Mexico, which was recorded by a temporary seismic deployment known as the MARS network. We located ∼1000 events from this seismic swarm. Previous earthquake swarms in the area were reported in the years 1997, 1999 and 2000. We relocate and analyze the evolution and properties of the 2006 earthquake swarm, employing a waveform cross-correlation-based phase repicking technique. Hypocenters from 911 events were located and divided into eighteen families having a correlation coefficient at or above 0.75. 90% of the earthquakes provide at least sixteen phase picks. We used the single-event location code Hypo71 and the P-wave velocity model used by the Jalisco Seismic and Accelerometer Network to improve hypocenters based on the correlation-adjusted phase arrival times. We relocated 121 earthquakes, which show clearly two clusters, between 9–10 km and 3–4 km depth respectively. The average location error estimates are

The internal structure of the Moon has been investigated over many years using a variety of seismic methods, such as travel time analysis, receiver functions, and tomography. Here we propose to apply body-wave seismic interferometry to deep moonquakes in order to retrieve zero-offset reflection responses (and thus images) beneath the Apollo stations on the nearside of the Moon from virtual sources colocated with the stations. This method is called deep-moonquake seismic interferometry (DMSI). Our results show a laterally coherent acoustic boundary around 50 km depth beneath all four Apollo stations. We interpret this boundary as the lunar seismic Moho. This depth agrees with Japan Aerospace Exploration Agency's (JAXA) SELenological and Engineering Explorer (SELENE) result and previous travel time analysis at the Apollo 12/14 sites. The deeper part of the image we obtain from DMSI shows laterally incoherent structures. Such lateral inhomogeneity we interpret as representing a zone characterized by strong scattering and constant apparent seismic velocity at our resolution scale (0.2–2.0 Hz).

Imaging subsurface rock formations or geological objects like oil and gas reservoirs, mineral deposits, cavities or even magmatic plumbing systems under active volcanoes has been for many years a major quest of geoscientists. Since these subsurface objects cannot be observed directly, different indirect methods have been developed. These methods are all based on variations of certain physical properties of the subsurface materials that can be detected from the ground surface or from boreholes. To determine the density distribution, a new imaging technique using cosmic-ray muon detectors deployed in a borehole has been developed and a first prototype of a borehole muon detector successfully tested. In addition to providing a static image of the subsurface density in three dimensions (or three-dimensional tomography), borehole muography can also inform on the variations of density with time, which recently became of major importance with the injection of large volumes of fluids, mainly water and CO 2 , in porous subsurface reservoirs (e.g. aquifer storage and recovery, wastewater disposal, enhanced oil recovery and carbon sequestration). This raises several concerns about the risk of leakage and the mechanical integrity of the reservoirs. Determining the field scale induced displacement of fluids by geophysical methods like muography is thus a priority. This article is part of the Theo Murphy meeting issue ‘Cosmic-ray muography'.

No conclusive evidence has been presented to date for tectonic tremor (TT) in the vicinity of central Chile, where the Nazca Plate is subducting beneath the South American Plate. Subduction in our experimental location (roughly 35.5° S, 70.5° W) is steep and fairly unobstructed compared to the flattened and more seismogenic behavior to the north. We seek to identify TT in our experimental area, whose geodynamics are comparable to tremor-rich subduction zones such as Cascadia and the Nankai Trough. Our method combines time-series visual inspection, frequency-spectrum analysis, waveform cross-correlation, and 3-component (3C) signal covariance to explore the presence of TT in this region. We have identified TT using stations in central Chile and the Malargue region, Argentina. The TT exhibits similar features to other TT observations worldwide. Waveform characteristics for the TT in our study, particularly dimension of the 3C signal covariance, vary as a function of apparent source location. The duration of one episode of identified TT was about 10 h, which may indicate that the plate interface where tremor generates is strongly coupled. We conclude that our observations reflect features of the local propagation, rather than the tremor source itself.

[1] The 15 June 1991 climactic eruption of Mount Pinatubo (Philippines) was followed by intense seismicity that remained at a high level for several months. We located 10,839 events recorded between 1 July and mid-December 1991. In contrast to the preeruptive seismicity which was focused in two groups below the summit area, posteruptive events were widely distributed below and around the volcano. The classification of the events indicates the presence of several large multiplets, and the application of relative relocation techniques to the similar events by calculating high-precision delays between traces outlines a number of clear seismogenic structures. We used different methods to confirm the validity of our results; these tests indicate that reliable features can be detected with a small monitoring network. While the main cluster of activity can be attributed to an intrusive process starting from below the 15 June crater, the volcanic origin of the seismic activity in the other areas is more difficult to establish. Away from the summit, relocations define streaks or planes which are oriented predominantly southwest-northeast, with in several cases the presence of northwest-southeast conjugate structures. Most of the composite focal mechanisms that we could determine indicate predominantly strike-slip, right-lateral faulting. Our results indicate that most of the seismicity that occurred after the 15 June eruption is related to the east-west regional compressional stress field related to the subduction. We suggest that the regional stress field induces seismicity along new or preexisting faults in the medium surrounding the volcano where the stress field was locally disturbed by the volcanic eruption.