scholarly journals Microseismicity, tectonic stress, and exhumation rates near the central Alpine Fault, New Zealand

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):  
◽  
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):  
K Michailos ◽  
Rupert Sutherland ◽  
John Townend ◽  
Martha Savage

© 2020. American Geophysical Union. All Rights Reserved. We investigate orogenic uplift rates and the thermal structure of the crust in the hanging wall of the Alpine Fault, New Zealand, using the hypocenters of 7,719 earthquakes that occurred in the central Southern Alps between late 2008 and early 2017, and previously published thermochronological data. We assume that the base of the seismogenic zone corresponds to a brittle-ductile transition at some fixed temperature, which we estimate by fitting the combined thermochronological data and distribution of seismicity using a multi-1-D approach. We find that exhumation rates vary from 1 to 8 mm/yr, with maximum values observed in the area of highest topography near Aoraki/Mount Cook, a finding consistent with previous geologic and geodetic analyses. We estimate the temperature of the brittle-ductile transition beneath the Southern Alps to be 410–430°C, which is higher than expected for Alpine Fault rocks whose bulk lithology is likely dominated by quartz. The high estimated temperatures at the base of the seismogenic zone likely reflect the unmodeled effects of high fluid pressures or strain rates.


2021 ◽  
Author(s):  
K Michailos ◽  
Rupert Sutherland ◽  
John Townend ◽  
Martha Savage

© 2020. American Geophysical Union. All Rights Reserved. We investigate orogenic uplift rates and the thermal structure of the crust in the hanging wall of the Alpine Fault, New Zealand, using the hypocenters of 7,719 earthquakes that occurred in the central Southern Alps between late 2008 and early 2017, and previously published thermochronological data. We assume that the base of the seismogenic zone corresponds to a brittle-ductile transition at some fixed temperature, which we estimate by fitting the combined thermochronological data and distribution of seismicity using a multi-1-D approach. We find that exhumation rates vary from 1 to 8 mm/yr, with maximum values observed in the area of highest topography near Aoraki/Mount Cook, a finding consistent with previous geologic and geodetic analyses. We estimate the temperature of the brittle-ductile transition beneath the Southern Alps to be 410–430°C, which is higher than expected for Alpine Fault rocks whose bulk lithology is likely dominated by quartz. The high estimated temperatures at the base of the seismogenic zone likely reflect the unmodeled effects of high fluid pressures or strain rates.


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>


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;


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>


2021 ◽  
Author(s):  
Duna Roda-Boluda ◽  
Taylor Schildgen ◽  
Hella Wittmann-Oelze ◽  
Stefanie Tofelde ◽  
Aaron Bufe ◽  
...  

&lt;p&gt;The Southern Alps of New Zealand are the expression of the oblique convergence between the Pacific and Australian plates, which move at a relative velocity of nearly 40 mm/yr. This convergence is accommodated by the range-bounding Alpine Fault, with a strike-slip component of ~30-40 mm/yr, and a shortening component normal to the fault of ~8-10 mm/yr. While strike-slip rates seem to be fairly constant along the Alpine Fault, throw rates appear to vary considerably, and whether the locus of maximum exhumation is located near the fault, at the main drainage divide, or part-way between, is still debated. These uncertainties stem from very limited data characterizing vertical deformation rates along and across the Southern Alps. Thermochronology has constrained the Southern Alps exhumation history since the Miocene, but Quaternary exhumation is hard to resolve precisely due to the very high exhumation rates. Likewise, GPS surveys estimate a vertical uplift of ~5 mm/yr, but integrate only over ~10 yr timescales and are restricted to one transect across the range.&lt;/p&gt;&lt;p&gt;To obtain insights into the Quaternary distribution and rates of exhumation of the western Southern Alps, we use new &lt;sup&gt;10&lt;/sup&gt;Be catchment-averaged erosion rates from 20 catchments along the western side of the range. Catchment-averaged erosion rates span an order of magnitude, between ~0.8 and &gt;10 mm/yr, but we find that erosion rates of &gt;10 mm/yr, a value often quoted in the literature as representative for the entire range, are very localized. Moreover, erosion rates decrease sharply north of the intersection with the Marlborough Fault System, suggesting substantial slip partitioning. These &lt;sup&gt;10&lt;/sup&gt;Be catchment-averaged erosion rates integrate, on average, over the last ~300 yrs. Considering that the last earthquake on the Alpine Fault was in 1717, these rates are representative of inter-seismic erosion. Lake sedimentation rates and coseismic landslide modelling suggest that long-term (~10&lt;sup&gt;3&lt;/sup&gt; yrs) erosion rates over a full seismic cycle could be ~40% greater than our inter-seismic erosion rates. If we assume steady state topography, such a scaling of our &lt;sup&gt;10&lt;/sup&gt;Be erosion rate estimates can be used to estimate rock uplift rates in the Southern Alps. Finally, we find that erosion, and hence potentially exhumation, does not seem to be localized at a particular distance from the fault, as some tectonic and provenance studies have suggested. Instead, we find that superimposed on the primary tectonic control, there is an elevation/temperature control on erosion rates, which is probably transient and related to frost-cracking and glacial retreat.&lt;/p&gt;&lt;p&gt;Our results highlight the potential for &lt;sup&gt;10&lt;/sup&gt;Be catchment-averaged erosion rates to provide insights into the magnitude and distribution of tectonic deformation rates, and the limitations that arise from transient erosion controls related to the seismic cycle and climate-modulated surface processes.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;


Geomorphology ◽  
2016 ◽  
Vol 263 ◽  
pp. 71-86 ◽  
Author(s):  
T.R. Robinson ◽  
T.R.H. Davies ◽  
T.M. Wilson ◽  
C. Orchiston

Sign in / Sign up

Export Citation Format

Share Document