Geology topics

Martha K. Savage

Publications and source records attributed to Martha K. Savage.

6 recordsLinked to original sources

USGS2023

Ambient seismic noise tomography of the Kingdom of Saudi Arabia

Harrat Rahat is a Cenozoic volcanic field in the west-central part of the Kingdom of Saudi Arabia, 150 kilometers east of the Red Sea, and is the site of the most recent eruption in the country (1256 C.E.; 654 in the year of the Hijra). The city of Al Madīnah lies at the north end of Harrat Rahat, and its volcanic and seismic risks are frequently reassessed. In 2009 C.E. an earthquake swarm at Harrat Lunayyir, 200 km west-northwest of Al Madīnah, produced significant shaking and some building damage in nearby communities, which prompted a revision of seismic hazard models for the region. A network of seismic stations on this volcanically active western side of the Arabia Plate was installed, and stations were also added in the tectonically active northern part of the country. Although regional earthquakes may be used to determine the crustal structure of the western Arabia Plate, such crustal models are hindered by insufficient numbers of earthquakes in the stable plate interior. Tomography studies can be used to infer material properties of the subsurface, such as presence of partial melt, and are beneficial for volcanic hazard assessments. We use ambient seismic noise to compute Rayleigh and Love surface-wave dispersion maps between 5 and 12 second periods for a subset of seismic stations within and near northern Harrat Rahat. The surface-wave maps are inverted to produce shear-wave velocities using a neighborhood algorithm and interpolated into a pseudo three-dimensional model. The distributions of surface-wave and shear-wave velocities are heterogenous, varying from ±3–8 percent. However, low velocities are not restricted to Harrat Rahat. We observed a difference between Rayleigh- and Love-wave velocities that extends north of the site of the 1256 C.E. eruption and coincides with a low gravity anomaly. We obtain a shear-wave velocity increase of 10–15 percent between 15- and 25-kilometers depth, which is consistent with the presence of a transition between the felsic upper crust and the mafic lower crust of the Arabian Shield. The average shear-wave velocities of the upper and lower crust are estimated to be 3.64 and 3.95 kilometers per second using Rayleigh waves and 3.53 and 4.16 kilometers per second using Love waves, which are in good agreement with the results of other geophysical surveys in this area. The modest magnitude of the low-velocity anomalies within the crust and their locations extending well beyond the limits of Harrat Rahat indicate that they are not caused by a crustal magma chamber. If magma chambers exist, they are smaller than can be imaged with our seismological method (resolution corresponds to a 15-kilometer wavelength), deeper than 30 kilometers, or shallower than 5 kilometers with a small velocity contrast. The last possibility is discounted by weak-to-absent surface geothermal activity and by a lack of shallow seismicity that characterizes other areas with known shallow magmas, such as Hawaiʻi. We then expanded our work to the entire Arabian Shield with principally the same methodology to look at large-scale regional patterns. We found modest shear-wave velocity deviations on the order of only ±3 percent, which are within expected ranges for lithological variation.

Professional Paper↗

USGS2020

Global quieting of high-frequency seismic noise due to COVID-19 pandemic lockdown measures

Human activity causes vibrations that propagate into the ground as high-frequency seismic waves. Measures to mitigate the COVID-19 pandemic caused widespread changes in human activity, leading to a months-long reduction in seismic noise of up to 50%. The 2020 seismic noise quiet period is the longest and most prominent global anthropogenic seismic noise reduction on record. While the reduction is strongest at surface seismometers in populated areas, this seismic quiescence extends for many kilometers radially and hundreds of meters in depth. This provides an opportunity to detect subtle signals from subsurface seismic sources that would have been concealed in noisier times and to benchmark sources of anthropogenic noise. A strong correlation between seismic noise and independent measurements of human mobility suggests that seismology provides an absolute, real-time estimate of population dynamics.

Science↗

USGS2017

Bedrock geology of DFDP-2B, central Alpine Fault, New Zealand

During the second phase of the Alpine Fault, Deep Fault Drilling Project (DFDP) in the Whataroa River, South Westland, New Zealand, bedrock was encountered in the DFDP-2B borehole from 238.5–893.2 m Measured Depth (MD). Continuous sampling and meso- to microscale characterisation of whole rock cuttings established that, in sequence, the borehole sampled amphibolite facies, Torlesse Composite Terrane-derived schists, protomylonites and mylonites, terminating 200–400 m above an Alpine Fault Principal Slip Zone (PSZ) with a maximum dip of 62°. The most diagnostic structural features of increasing PSZ proximity were the occurrence of shear bands and reduction in mean quartz grain sizes. A change in composition to greater mica:quartz + feldspar, most markedly below c. 700 m MD, is inferred to result from either heterogeneous sampling or a change in lithology related to alteration. Major oxide variations suggest the fault-proximal Alpine Fault alteration zone, as previously defined in DFDP-1 core, was not sampled.

New Zealand Journal of Geology and Geophysics↗

USGS2010

Anisotropy, repeating earthquakes, and seismicity associated with the 2008 eruption of Okmok Volcano, Alaska

We use shear wave splitting (SWS) analysis and double-difference relocation to examine temporal variations in seismic properties prior to and accompanying magmatic activity associated with the 2008 eruption of Okmok volcano, Alaska. Using bispectrum cross-correlation, a multiplet of 25 earthquakes is identified spanning five years leading up to the eruption, each event having first motions compatible with a normal fault striking NE–SW. Cross-correlation differential times are used to relocate earthquakes occurring between January 2003 and February 2009. The bulk of the seismicity prior to the onset of the eruption on 12 July 2008 occurred southwest of the caldera beneath a geothermal field. Earthquakes associated with the onset of the eruption occurred beneath the northern portion of the caldera and started as deep as 13 km. Subsequent earthquakes occurred predominantly at 3 km depth, coinciding with the depth at which the magma body has been modeled using geodetic data. Automated SWS analysis of the Okmok catalog reveals radial polarization outside the caldera and a northwest-southeast polarization within. We interpret these polarizations in terms of a magma reservoir near the center of the caldera, which we model with a Mogi point source. SWS analysis using the same input processing parameters for each event in the multiplet reveals no temporal changes in anisotropy over the duration of the multiplet, suggesting either a short-term or small increase in stress just before the eruption that was not detected by GPS, or eruption triggering by a mechanism other than a change of stress in the system.

Alaska↗

USGS2010

Earthquake refraction profiles of the root of the Sierra Nevada

We examine the seismic structure of the Sierra Nevada using records of nine earthquakes and one explosion in and near the Sierra, recorded on stations in the Sierra. We first interpret travel times from these paths, which are confined to a single tectonic block, in terms of one‐dimensional structures. The most nearly reversed pair of earthquakes, the 1966 Truckee and 1983 Durrwood Meadows earthquakes, share refracted ( P n ) arrival times (corrected to surface focus) along a line t =8.75±0.25+Δ/8.0, suggesting that a nearly flat layer of 8.0 km/s mantle material lies at depths of 46–48 km. First arrivals from these events do not constrain velocities from ≈30 to 45 km depth. Secondary arrivals and some first arrivals from other earthquakes suggest that velocities in part of this region range between 6.9 and 7.8 km/s. The presence of this “7.x‐km/s” layer can help to explain previous contradictory observations. The 7.x‐km/s layer could be interpreted as either the mafic bottom of a silicic, Mesozoic magmatic arc or as accreted mafic underplating or rejuvenated mantle related to Cenozoic arc volcanism or Basin and Range spreading. Arrivals at stations in the foothills and the crest of the Sierra cannot be fit with a single longitudinal structure, indicating a lateral variation of velocity structure. These variations support previously inferred variations of lithospheric structure, with higher‐velocity, thinner crust to the west beneath the Sierran foothills and slower‐velocity crust (or possibly upper mantle) beneath the high mountains in the eastern Sierra. Rapid changes in arrival times between stations separated by short distances in the eastern Sierra suggest that a sharp boundary exists between the Sierra and the Basin and Range at Moho depths. We also present fresh evidence of the asymmetry of the root of the Sierra, wherein arrivals from earthquakes on the west of the Sierra are delayed within the Sierra and return to original values in the Basin and Range, while arrivals from earthquakes and explosions from the Sierra into the Great Valley. We suggest that if the 7.x‐km/s material occurs in a wedge above the Moho, then the asymmetry can be explained by arrivals from the west being delayed by the dipping 8.0‐km/s Moho, while those from the east may be entering the root along a 7.x‐km/s layer that is near the depth of the Basin and Range Moho.

New Hampshire↗

USGS1994

Earthquakes in California and Nevada

Each year seismic networks record more than 30,000 earthquakes with magnitudes (M) as small as 1.0 in California and Nevada. Historic accounts indicate more than 360 earthquakes of M ≥5.5 have occurred in this region between 1836 and 1994. This map depicts the recent and historic seismicity in three ways. The small red dots represent epicenters of earthquakes between 1 January 1980 and 31 January 1994. The black circles depict all known epicenters of M≥5.5 earthquakes since 1836. A number adjacent to an epicenter location refers to an entry in the table "Known M≥6.5 Earthquakes since 1836." The hachured regions indicate the surface projection of fault rupture for significant earthquakes as inferred from geologic observations of surface deformation and estimates of subsurface faulting from geophysical methods. Tightly spaced hachures represent vertical, strike-slip faults, and wider spaced hachures designate slipping faults.

California, Nevada↗