Using surface waves recorded by a large mesh of three-element arrays to detect and locate disparate seismic sources

Author Posting. © The Authors, 2018. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 215 (2018): 942–958, doi:10.1093/gji/ggy316.Surface waves recorded by global arrays have proven useful for locating tectonic earthquakes and in detecting slip events depleted in high frequency, such as glacial quakes. We develop a novel method using an aggregation of small- to continental-scale arrays to detect and locate seismic sources with Rayleigh waves at 20–50 s period. The proposed method is a hybrid approach including first dividing a large aperture aggregate array into Delaunay triangular subarrays for beamforming, and then using the resolved surface wave propagation directions and arrival times from the subarrays as data to formulate an inverse problem to locate the seismic sources and their origin times. The approach harnesses surface wave coherence and maximizes resolution of detections by combining measurements from stations spanning the whole U.S. continent. We tested the method with earthquakes, glacial quakes and landslides. The results show that the method can effectively resolve earthquakes as small as ∼M3 and exotic slip events in Greenland. We find that the resolution of the locations is non-uniform with respect to azimuth, and decays with increasing distance between the source and the array when no calibration events are available. The approach has a few advantages: the method is insensitive to seismic event type, it does not require a velocity model to locate seismic sources, and it is computationally efficient. The method can be adapted to real-time applications and can help in identifying new classes of seismic sources.WF is currently supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship. This work was supported by National Science Foundation grant EAR-1358520 at Scripps Institution of Oceanography, UC San Diego

Relationships between gravity waves observed at Earth's surface and in the stratosphere over the central and eastern United States

Observations of tropospheric gravity waves (GWs) made by the new and extensive USArray Transportable Array (TA) barometric network located east of the Rockies, in the central and eastern United States and of stratospheric (30–40 km above sea level) GWs made by the Atmospheric Infrared Sounder (AIRS) are compared over a 5 year time span from 2010 through 2014. GW detections in the period band from 2 to 6 h made at the Earth's surface during the thunderstorm season from May through August each year exhibit the same broad spatial and temporal patterns as observed at stratospheric altitudes. At both levels, the occurrence frequency of GWs is higher at night than during the day and is highest to the west of the Great Lakes. Statistically significant correlations between the variance of the pressure at the TA, which is a proxy for GWs at ground level, with 8.1 μm brightness temperature measurements from AIRS and rain radar precipitation data, which are both proxies for convective activity, indicate that GWs observed at the TA are related to convective sources. There is little, if any, time lag between the two. Correlations between GWs in the stratosphere and at ground level are weaker, possibly due to the AIRS observational filter effect, but are still statistically significant at nighttime. We conclude that convective activity to the west of the Great Lakes is the dominant source of GWs both at ground level and within the stratosphere

EON-ROSE and the Canadian Cordillera Array – Building Bridges to Span Earth System Science in Canada

EON-ROSE (Earth-System Observing Network - Réseau d’Observation du Système terrestrE) is a new initiative for a pan-Canadian research collaboration to holistically examine Earth systems from the ionosphere into the core. The Canadian Cordillera Array (CC Array) is the pilot phase, and will extend across the Cordillera from the Beaufort Sea to the U.S. border. The vision for EON-ROSE is to install a network of telemetered observatories to monitor solid Earth, environmental and atmospheric processes. EON-ROSE is an inclusive, combined effort of Canadian universities, federal, provincial and territorial government agencies, industry, and international collaborators. Brainstorming sessions and several workshops have been held since May 2016. The first station will be installed at Kluane Lake Research Station in southwestern Yukon during the summer of 2018. The purpose of this report is to provide a framework for continued discussion and development.RÉSUMÉEON-ROSE (Earth-System Observing Network - Réseau d’Observation du Système terrestrE) est une nouvelle initiative de collaboration de recherche pancanadienne visant à étudier de manière holistique les systèmes terrestres, depuis l’ionosphère jusqu’au noyau. Le Réseau canadien de la cordillère (CC Array) en est la phase pilote, laquelle couvrira toute la Cordillère, de la mer de Beaufort jusqu’à la frontière étasunienne. L’objectif d’EON-ROSE est d’installer un réseau d’observatoires télémétriques pour suivre en continu les processusterrestres, environnementaux et atmosphériques. EON-ROSE est un effort combiné et inclusif des universités canadiennes, des organismes gouvernementaux fédéraux, provinciaux et territoriaux, de l’industrie et de collaborateurs internationaux. Des séances de remue-méninges et plusieurs ateliers ont été tenus depuis mai 2016. La première station sera installée à la station de recherche du lac Kluane, dans le sud-ouest du Yukon, au cours de l’été 2018. Le but du présent rapport est de fournir un cadre de discussion et de développement continu

Seafloor seismic monitoring of an active submarine volcano : local seismicity at Vailulu'u Seamount, Samoa

Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 5 (2004): Q06007, doi:10.1029/2004GC000702.We deployed five ocean bottom hydrophones (OBHs) for a 1-year seismic monitoring study of Vailulu'u Seamount, the youngest and easternmost volcano in the Samoan Archipelago. Four instruments were placed on the summit crater rim at 600–700 m water depth, and one was placed inside the crater at 1000 m water depth. An analysis of the first 45 days of records shows a very large number of seismic events, 211 of them local. These events define a steady background activity of about four seismic events per day, increasing to about 10 events per day during a week of heightened seismic activity, which peaked at 40 events during 1 day. We identified 107 earthquakes, whose arrivals could be picked on all five stations and that are likely located within the seamount, based on their similar waveforms. Two linear trends are defined by 21 of these events. These are extremely well correlated and located, first downward then upward on a steeply inclined plane that is close to the axial plane of the southeast rift as it emerges from the main summit of Vailulu'u. These events resemble volcanotectonic earthquakes from subaerial volcanoes in displaying very coherent seismic waveforms and by showing systematic, narrowly defined progressions in hypocenter locations. We propose that these events reflect brittle rock failure due to magma redistribution in or near a central magma reservoir.The bulk of this work was funded by NSF-OCE, in grants to HS and SRH and the OBSIP facility at Scripps

Optimization and geophysical inverse problems

A fundamental part of geophysics is to make inferences about the interior of the earth on the basis of data collected at or near the surface of the earth. In almost all cases these measured data are only indirectly related to the properties of the earth that are of interest, so an inverse problem must be solved in order to obtain estimates of the physical properties within the earth. In February of 1999 the U.S. Department of Energy sponsored a workshop that was intended to examine the methods currently being used to solve geophysical inverse problems and to consider what new approaches should be explored in the future. The interdisciplinary area between inverse problems in geophysics and optimization methods in mathematics was specifically targeted as one where an interchange of ideas was likely to be fruitful. Thus about half of the participants were actively involved in solving geophysical inverse problems and about half were actively involved in research on general optimization methods. This report presents some of the topics that were explored at the workshop and the conclusions that were reached. In general, the objective of a geophysical inverse problem is to find an earth model, described by a set of physical parameters, that is consistent with the observational data. It is usually assumed that the forward problem, that of calculating simulated data for an earth model, is well enough understood so that reasonably accurate synthetic data can be generated for an arbitrary model. The inverse problem is then posed as an optimization problem, where the function to be optimized is variously called the objective function, misfit function, or fitness function. The objective function is typically some measure of the difference between observational data and synthetic data calculated for a trial model. However, because of incomplete and inaccurate data, the objective function often incorporates some additional form of regularization, such as a measure of smoothness or distance from a prior model. Various other constraints may also be imposed upon the process. Inverse problems are not restricted to geophysics, but can be found in a wide variety of disciplines where inferences must be made on the basis of indirect measurements. For instance, most imaging problems, whether in the field of medicine or non-destructive evaluation, require the solution of an inverse problem. In this report, however, the examples used for illustration are taken exclusively from the field of geophysics. The generalization of these examples to other disciplines should be straightforward, as all are based on standard second-order partial differential equations of physics. In fact, sometimes the non-geophysical inverse problems are significantly easier to treat (as in medical imaging) because the limitations on data collection, and in particular on multiple views, are not so severe as they generally are in geophysics. This report begins with an introduction to geophysical inverse problems by briefly describing four canonical problems that are typical of those commonly encountered in geophysics. Next the connection with optimization methods is made by presenting a general formulation of geophysical inverse problems. This leads into the main subject of this report, a discussion of methods for solving such problems with an emphasis upon newer approaches that have not yet become prominent in geophysics. A separate section is devoted to a subject that is not encountered in all optimization problems but is particularly important in geophysics, the need for a careful appraisal of the results in terms of their resolution and uncertainty. The impact on geophysical inverse problems of continuously improving computational resources is then discussed. The main results are then brought together in a final summary and conclusions section

Plate-boundary deformation associated with the great Sumatra–Andaman earthquake

The Sumatra–Andaman earthquake of 26 December 2004 is the first giant earthquake (moment magnitude M_w > 9.0) to have occurred since the advent of modern space-based geodesy and broadband seismology. It therefore provides an unprecedented opportunity to investigate the characteristics of one of these enormous and rare events. Here we report estimates of the ground displacement associated with this event, using near-field Global Positioning System (GPS) surveys in northwestern Sumatra combined with in situ and remote observations of the vertical motion of coral reefs. These data show that the earthquake was generated by rupture of the Sunda subduction megathrust over a distance of >1,500 kilometres and a width of <150 kilometres. Megathrust slip exceeded 20 metres offshore northern Sumatra, mostly at depths shallower than 30 kilometres. Comparison of the geodetically and seismically inferred slip distribution indicates that ~30 per cent additional fault slip accrued in the 1.5 months following the 500-second-long seismic rupture. Both seismic and aseismic slip before our re-occupation of GPS sites occurred on the shallow portion of the megathrust, where the large Aceh tsunami originated. Slip tapers off abruptly along strike beneath Simeulue Island at the southeastern edge of the rupture, where the earthquake nucleated and where an M_w = 7.2 earthquake occurred in late 2002. This edge also abuts the northern limit of slip in the 28 March 2005 M_w = 8.7 Nias–Simeulue earthquake