Geology topics

John Townend

Publications and source records attributed to John Townend.

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↗

USGS2018

Frictional properties and 3-D stress analysis of the southern Alpine Fault, New Zealand

New Zealand's Alpine Fault (AF) ruptures quasi-periodically in large-magnitude earthquakes. Paleoseismological evidence suggests that about half of all recognized AF earthquakes terminated at the boundary between the Central and South Westland sections of the fault. There, fault geometry and the polarity of uplift change. The South Westland AF exhibits oblique-normal fault motion on a structure oriented 052°/82°SE that, for at least 35 km along strike, contains saponite-rich principal slip zone gouges. New hydrothermal friction experiments reveal that the saponite fault gouge is frictionally weak, exhibiting friction coefficients between μ = 0.12 and μ = 0.16 for a range of temperatures ( T = 25–210 °C) and effective normal stresses ( σ n ' = 31.2–93.6 MPa). The saponite gouge is rate-strengthening in all velocity steps performed at velocities between 0.01 and 3.0 μm/s, behavior conducive to aseismic creep. A three-dimensional stress analysis shows that the South Westland AF is favorably oriented with respect to the regional stress field for slip within the frictionally weak saponite fault gouge. Geometrically, the fault is severely misoriented for slip in any fault-forming materials with friction coefficients exceeding μ ∼0.5. The combination of weak gouges prone to aseismic creep, strong asperities , and low resolved shear stress may impede earthquake rupture propagation along the South Westland Alpine Fault.

Journal of Structural Geology↗

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↗

USGS2014

Low-frequency earthquakes reveal punctuated slow slip on the deep extent of the Alpine Fault, New Zealand

We present the first evidence of low-frequency earthquakes (LFEs) associated with the deep extension of the transpressional Alpine Fault beneath the central Southern Alps of New Zealand. Our database comprises a temporally continuous 36 month-long catalog of 8760 LFEs within 14 families. To generate this catalog, we first identify 14 primary template LFEs within known periods of seismic tremor and use these templates to detect similar events in an iterative stacking and cross-correlation routine. The hypocentres of 12 of the 14 LFE families lie within 10 km of the inferred location of the Alpine Fault at depths of approximately 20–30 km, in a zone of high P-wave attenuation, low P-wave speeds, and high seismic reflectivity. The LFE catalog consists of persistent, discrete events punctuated by swarm-like bursts of activity associated with previously and newly identified tremor periods. The magnitudes of the LFEs range between M L – 0.8 and M L 1.8, with an average of M L 0.5. We find that the frequency-magnitude distribution of the LFE catalog both as a whole and within individual families is not consistent with a power law, but that individual families' frequency-amplitude distributions approximate an exponential relationship, suggestive of a characteristic length-scale of failure. We interpret this LFE activity to represent quasi-continuous slip on the deep extent of the Alpine Fault, with LFEs highlighting asperities within an otherwise steadily creeping region of the fault.

Alpine Fault↗

USGS2014

Frictional properties of exhumed fault gouges in DFDP-1 cores, Alpine Fault, New Zealand

Principal slip zone gouges recovered during the Deep Fault Drilling Project (DFDP-1), Alpine Fault, New Zealand, were deformed in triaxial friction experiments at temperatures, T , of up to 350°C, effective normal stresses, σ n ′, of up to 156 MPa, and velocities between 0.01 and 3 µm/s. Chlorite/white mica-bearing DFDP-1A blue gouge, 90.62 m sample depth, is frictionally strong (friction coefficient, μ , 0.61–0.76) across all experimental conditions tested ( T = 70–350°C, σ n ′ = 31.2–156 MPa); it undergoes a transition from positive to negative rate dependence as T increases past 210°C. The friction coefficient of smectite-bearing DFDP-1B brown gouge, 128.42 m sample depth, increases from 0.49 to 0.74 with increasing temperature and pressure ( T = 70–210°C, σ n ′ = 31.2–93.6 MPa); the positive to negative rate dependence transition occurs as T increases past 140°C. These measurements indicate that, in the absence of elevated pore fluid pressures, DFDP-1 gouges are frictionally strong under conditions representative of the seismogenic crust.

Alpine Fault↗

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↗