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

2021 ◽  
Author(s):  
Calum Chamberlain ◽  
D Shelly ◽  
John Townend ◽  
Timothy Stern

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 ML - 0.8 and ML 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. © 2014. American Geophysical Union. All Rights Reserved.

2021 ◽  
Author(s):  
Calum Chamberlain ◽  
D Shelly ◽  
John Townend ◽  
Timothy Stern

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 ML - 0.8 and ML 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. © 2014. American Geophysical Union. All Rights Reserved.


2021 ◽  
Author(s):  
V Lay ◽  
S Buske ◽  
SB Bodenburg ◽  
John Townend ◽  
R Kellett ◽  
...  

No description supplied


2021 ◽  
Author(s):  
◽  
Konstantinos Michailos

<p>This thesis documents a detailed examination of the seismic activity and characteristics of crustal deformation along the central Alpine Fault, a major obliquely convergent plate-boundary fault. Paleoseismic evidence has established that the Alpine Fault produces large to great (M7−8) earthquakes every 250−300 years, in a quasi-periodic manner, with the last surface-rupturing earthquake occurring in 1717. This renders the fault late in its typical earthquake cycle, posing substantial seismic risk to southern and central New Zealand. Understanding the seismic and tectonic character of this fault may yield information of both societal and scientific significance regarding seismic hazard and late-interseismic processes leading up to a large earthquake. However, the central Alpine Fault is currently seismically quiescent when compared to adjacent regions, and therefore requires detailed, long-duration observations to study seismotectonic processes. The work in this thesis addresses the need for a greater understanding of along-strike variations in seismic character of the Alpine Fault ahead of an anticipated large earthquake.  To achieve observations with high spatial and temporal resolution across the length of the central Alpine Fault, I use 8.5 years of continuous seismic data from the Southern Alps Microearthquake Borehole Array (SAMBA), and data from four other temporary seismic networks and five local GeoNet permanent sites. Incorporating all of these temporary and permanent seismic sites provides us with a dense composite network of seismometers. Without such a dense network, homogeneous examination of the characteristics of low-magnitude seismicity near the Alpine Fault would be impossible.  Using this dataset, I have constructed the most extensive microearthquake catalog for the central Alpine Fault region to date, containing 9,111 earthquakes and covering the time between late 2008 and early 2017. To construct this catalog I created an objective workflow to ensure catalog uniformity. Overall, 7,719 earthquakes were successfully relocated with location uncertainties generally ≤ 0.5 km in both the horizontal and vertical directions. The majority of the earthquakes were found to occur southeast of the Alpine Fault (i.e. in the hanging-wall). I observed a lack of seismicity beneath Aoraki/Mount Cook that has previously been shown to be associated with locally high uplift rates (6–10 mm/yr) and high geothermal gradients (∼60◦C/km). Seismogenic cut-off depths were observed to significantly vary along the strike of the Alpine Fault, ranging from 8 km beneath the highest topography to 20 km in the adjacent areas.  To quantify the scale of the seismic deformation, a new local magnitude scale was also derived, corrected for geometric spreading, attenuation and site terms based on individually calculated GeoNet moment magnitude (Mw) values. Earthquake local magnitudes range between ML –1.2 and 4.6 and the catalog is complete above ML 1.1.  To examine the stress regime near the central Alpine Fault, I built a new data set of 845 focal mechanisms from earthquakes in our catalog. This was achieved by manually determining P wave arrival polarity picks from all earthquakes larger than ML 1.5. In order to determine the orientations and characteristics of the stress parameters, I grouped these focal mechanisms and performed stress inversion calculations that provided an average maximum horizontal compressive stress orientation, SHmax, of 121±11◦ , which is uniform within uncertainty along the length of the central Southern Alps. I observed an average angle of 65◦ between the SHmax and the strike of the Alpine Fault, which is consistent with results from similar previous studies in the northern and southern sections of the Alpine Fault. This implies that the Alpine Fault is misoriented for reactivation, in the prevailing stress field.  Using a 1-D steady-state thermal structure model constrained by seismicity and thermochronology data, I investigated the crustal thermal structure and vertical kinematics of the central Southern Alps orogen. The short-term seismicity data and longer-term thermochronology data impose complementary constraints on the model. I observed a large variation in exhumation rate estimates (1–8 mm/yr) along the length of the Alpine Fault, with maximum calculated values observed near Aoraki/Mount Cook. I calculated the temperature at the brittle-ductile transition zone, which ranges from 440 to 457◦C in the different models considered. This temperature is slightly hotter than expected for crust composed by quartz-rich rocks, but consistent with the presence of feldspar-rich mafic rocks in parts of the crust.</p>


2021 ◽  
Author(s):  
◽  
Laura-May Baratin Wachten

<p>This thesis involves the study of low-frequency earthquakes (LFEs) in the central Southern Alps. The Alpine Fault is the principal locus of deformation within the Australia–Pacific plate boundary in the South Island of New Zealand and it is late in its typical ∼300-year seismic cycle. Surveying the seismicity associated with slow deformation in the vicinity of the Alpine Fault may provide constraints on the stresses acting on a major transpressive margin prior to an anticipated great (≥M8) earthquake. Here, we use 8 years of data from the Southern Alps Microearthquake Borehole Array (SAMBA) (amongst those, 3 years of data were collected as part of this project) in order to: (1) generate an updated LFE catalogue using an improved matched-filter technique that incorporates phase-weighted stacking; (2) compute LFE focal mechanisms and invert them to infer the crustal stress field on the deep extent of the Alpine Fault; (3) expand the LFE catalogue to cover a wider range of spatial/temporal behaviours; (4) study LFE families’ characteristics to identify periods where slow slip might happen.  We first use fourteen primary LFE templates in an iterative matched-filter and stacking routine, which allows the detection of similar signals and produces LFE families sharing common locations. We generate an 8-yr catalogue containing 10,000 LFEs that are combined for each of the 14 LFE families using phase-weighted stacking to produce signals with the highest possible signal-to-noise ratios. We find LFEs to occur almost continuously during the 8-yr study period and we highlight two types of LFE distributions: (1) discrete behaviour with an inter-event time exceeding 2 minutes; (2) burst-like behaviour with an inter-event time below 2 minutes. The discrete events are interpreted as small-scale frequent deformation on the deep extent of the Alpine Fault and the LFE bursts (corresponding in most cases to known episodes of tremor or large regional earthquakes) are interpreted as brief periods of increased slip activity indicative of slow slip. We compute improved non-linear earthquake locations using a 3D velocity model and find LFEs to occur below the seismogenic zone at depths of 17–42 km, on or near the hypothesised deep extent of the Alpine Fault. We then compute the first estimates of LFE focal mechanisms associated with continental faulting. Focal mechanisms, in conjunction with recurrence intervals, are consistent with quasi-continuous shear faulting on the deep extent of the Alpine Fault.  We then generate a new catalogue that regroups hundreds of LFE families. This time 638 synthetic LFE waveforms are generated using a 3D grid and used as primary templates in a matched-filter routine. Of those, 529 templates yield enough detections during the first iteration of the matched-filter routine (≥ 500 detections over the 8-yr study period) and are kept for further analysis. We then use the best 25% of correlated events for each LFE family to generate linear stacks which create new LFE templates. From there, we run a second and final iteration of the matched-filter routine with the new LFE templates to obtain our final LFE catalogue. The remaining 529 templates detect between 150 and 1,671 events each totalling 300,996 detections over the 8-yr study period. Of those 529 LFEs, we manage to locate 378 families. Their depths range between 11 and 60 km and LFEs locate mainly in the southern part of the SAMBA network. We finally examine individual LFE family rates and occurrence patterns. They indicate that LFE sources seem to evolve from an episodic or ‘stepped’ to a continuous behaviour with depth. This transition may correspond to an evolution from a stick-slip to a stable-sliding slip regime. Hence, we propose that the distinctive features of LFE occurrence patterns reflect variations in the in-situ stress and frictional conditions at the individual LFE source locations on the Alpine Fault.  Finally, we use this new extensive catalogue as a tool for in-depth analyses of the deep central Alpine Fault structure and its slip behaviour. We identify eight episodes of increased LFE activity between 2009 and 2017 and provide time windows for further investigations of tremor and slow slip. We also study the spatial and temporal behaviours of LFEs and find that LFEs with synchronous occurrence patterns tend to be clustered in space. We thus suggest that individual LFE sources form spatially coherent clusters that may represent localised asperities or elastic patches on the deep Alpine Fault interface. We infer that those clusters may have a similar rheological response to tectonic forcing or to potential slow slip events. Eventually, we discover slow (10km/day) and rapid (∼20-25km/h) migrations of LFEs along the Alpine Fault. The slow migration might be controlled by slow slip events themselves while the rapid velocities could be explained by the LFE sources’ intrinsic properties.</p>


1980 ◽  
Vol 70 (5) ◽  
pp. 1607-1633
Author(s):  
Christopher D. Stephens ◽  
John C. Lahr ◽  
Kent A. Fogleman ◽  
Robert B. Horner

abstract The St. Elias, Alaska, earthquake (Ms 7.1) of February 28, 1979 occurred beneath the Chugach and St. Elias Mountains of southeastern Alaska and southwestern Yukon Territory. The main shock and aftershocks were recorded at regional high-gain, high-frequency seismographs operated by the U.S. and Canada. Hypocenters and magnitudes are presented for 308 aftershocks that occurred between February 28 and March 31, 1979. These data contain a nearly complete record of events of magnitude 3.5 and larger starting about 20 min after the main shock. The largest aftershock has a poorly determined magnitude slightly above 5, and the frequency-magnitude distribution has a b value of 1.36. A t−p inverse power law with an unusually low value of 0.93 for p adequately describes the decay with time in the frequency of occurrence of large aftershocks. The aftershocks occurred in a broad zone that extends about 115 km southeast from the epicenter of the main shock. Events tend to form clusters within this zone. One of the most remarkable features in the distribution of epicenters is that relatively few aftershocks were located near the epicenter of the main shock, and the highest rate of activity was centered about 50 km southeast of the epicenter of the main shock. Within the accuracy of the data, the depths of the aftershocks are all less than about 20 km. In the few areas where good depth control is available, the seismicity appears not to extend to the Earth's surface. Additional data from temporary stations operated in the aftershock zone during July and August 1979 indicate that the seismicity in some areas may be confined to a zone less than 6 km in vertical thickness. Focal mechanisms determined from p-wave first motions for some of the larger aftershocks all indicate northward-directed compression, which is consistent with the focal mechanism of the main shock. A review of the regional seismicity during the 6-month period preceding the St. Elias earthquake indicates that, relative to a comparable 6-month period 1 yr earlier, there was a 45 per cent increase in the rate of activity for events of magnitude 1.8 and larger, and possibly a decrease in the b value during the same period. Also, a prominent cluster of events with magnitudes less than 3.3 occurred at the southeast corner of the aftershock zone during the 6 months prior to the earthquake. The seismic record from the USGS network since 1974 is not yet complete in time, so it is not possible to determine how unusual the seismic activity preceding this earthquake has been.


2020 ◽  
Author(s):  
Vera Lay ◽  
Stefan Buske ◽  
Sascha Barbara Bodenburg ◽  
Franz Kleine ◽  
John Townend ◽  
...  

&lt;p&gt;The Alpine Fault along the West Coast of the South Island (New Zealand) is a major plate boundary that is expected to rupture in the next 50 years, likely as a magnitude 8 earthquake. The Deep Fault Drilling Project (DFDP) aims to deliver insight into the geological structure of this fault zone and its evolution by drilling and sampling the Alpine Fault at depth. &amp;#160;&lt;/p&gt;&lt;p&gt;Here we present results from a 3D seismic survey around the DFDP-2 drill site in the Whataroa Valley where the drillhole penetrated almost down to the fault surface. Within the glacial valley, we collected 3D seismic data to constrain valley structures that were obscured in previous 2D seismic data. The new data consist of a 3D extended vertical seismic profiling (VSP) survey using three-component receivers and a fibre optic cable in the DFDP-2B borehole as well as a variety of receivers at the surface.&lt;/p&gt;&lt;p&gt;The data set enables us to derive a reliable 3D P-wave velocity model by first-arrival travel time tomography. We identify a 100-460 m thick sediment layer (average velocity 2200&amp;#177;400 m/s) above the basement (average velocity 4200&amp;#177;500 m/s). Particularly on the western valley side, a region of high velocities steeply rises to the surface and mimics the topography. We interpret this to be the infilled flank of the glacial valley that has been eroded into the basement. In general, the 3D structures implied by the velocity model on the upthrown (Pacific Plate) side of the Alpine Fault correlate well with the surface topography and borehole findings.&lt;/p&gt;&lt;p&gt;A reliable velocity model is not only valuable by itself but it is also required as input for prestack depth migration (PSDM). We performed PSDM with a part of the 3D data set to derive a structural image of the subsurface within the Whataroa Valley. The top of the basement identified in the P-wave velocity model coincides well with reflectors in the migrated images so that we can analyse the geometry of the basement in detail.&lt;/p&gt;


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. D73-D79 ◽  
Author(s):  
Qiaomu Qi ◽  
Arthur C. H. Cheng ◽  
Yunyue Elita Li

ABSTRACT Formation S-wave attenuation, when combined with compressional attenuation, serves as a potential hydrocarbon indicator for seismic reservoir characterization. Sonic flexural wave measurements provide a direct means for obtaining the in situ S-wave attenuation at log scale. The key characteristic of the flexural wave is that it propagates at the formation shear slowness and experiences shear attenuation at low frequency. However, in a fast formation, the dipole log consists of refracted P- and S-waves in addition to the flexural wave. The refracted P-wave arrives early and can be removed from the dipole waveforms through time windowing. However, the refracted S-wave, which is often embedded in the flexural wave packet, is difficult to separate from the dipole waveforms. The additional energy loss associated with the refracted S-wave results in the estimated dipole attenuation being higher than the shear attenuation at low frequency. To address this issue, we have developed a new method for accurately determining the formation shear attenuation from the dipole sonic log data. The method uses a multifrequency inversion of the frequency-dependent flexural wave attenuation based on energy partitioning. We first developed our method using synthetic data. Application to field data results in a shear attenuation log that is consistent with lithologic interpretation of other available logs.


2016 ◽  
Vol 59 (4) ◽  
pp. 483-495 ◽  
Author(s):  
J Feenstra ◽  
C Thurber ◽  
J Townend ◽  
S Roecker ◽  
S Bannister ◽  
...  

Geophysics ◽  
1993 ◽  
Vol 58 (1) ◽  
pp. 20-29 ◽  
Author(s):  
Nabil Akbar ◽  
Jack Dvorkin ◽  
Amos Nur

To relate P‐wave attenuation to permeability, we examine a three‐dimensional (3-D) theoretical model of a cylindrical pore filled with viscous fluid and embedded in an infinite isotropic elastic medium. We calculate both attenuation and permeability as functions of the direction of wave propagation. Attenuation estimates are based on the squirt flow mechanism; permeability is calculated using the Kozeny‐Carman relation. We find that in the case when a plane P‐wave propagates perpendicular to the pore orientation [Formula: see text], attenuation is always higher than when a wave propagates parallel to this orientation [Formula: see text]. The ratio of these two attenuation values [Formula: see text] increases with an increasing pore radius and decreasing frequency and saturation. By changing permeability, varying the radius of the pore, we find that the permeability‐attenuation relation is characterized by a peak that shifts toward lower permeabilities as frequency decreases. Therefore, the attenuation of a low‐frequency wave decreases with increasing permeability. We observe a similar trend on relations between attenuation and permeability experimentally obtained on sandstone samples.


Sign in / Sign up

Export Citation Format

Share Document