SEARCH · USGS Science

Results for “Pure & Applied Geophysics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7Linked to original sources

USGS2010

Combined effects of tectonic and landslide-generated Tsunami Runup at Seward, Alaska during the Mw 9.2 1964 earthquake

We apply a recently developed and validated numerical model of tsunami propagation and runup to study the inundation of Resurrection Bay and the town of Seward by the 1964 Alaska tsunami. Seward was hit by both tectonic and landslide-generated tsunami waves during the M W "> M W MW 9.2 1964 megathrust earthquake. The earthquake triggered a series of submarine mass failures around the fjord, which resulted in landsliding of part of the coastline into the water, along with the loss of the port facilities. These submarine mass failures generated local waves in the bay within 5 min of the beginning of strong ground motion. Recent studies estimate the total volume of underwater slide material that moved in Resurrection Bay to be about 211 million m 3 (Haeussler et al. in Submarine mass movements and their consequences, pp 269–278, 2007 ). The first tectonic tsunami wave arrived in Resurrection Bay about 30 min after the main shock and was about the same height as the local landslide-generated waves. Our previous numerical study, which focused only on the local landslide-generated waves in Resurrection Bay, demonstrated that they were produced by a number of different slope failures, and estimated relative contributions of different submarine slide complexes into tsunami amplitudes (Suleimani et al. in Pure Appl Geophys 166:131–152, 2009 ). This work extends the previous study by calculating tsunami inundation in Resurrection Bay caused by the combined impact of landslide-generated waves and the tectonic tsunami, and comparing the composite inundation area with observations. To simulate landslide tsunami runup in Seward, we use a viscous slide model of Jiang and LeBlond (J Phys Oceanogr 24(3):559–572, 1994 ) coupled with nonlinear shallow water equations. The input data set includes a high resolution multibeam bathymetry and LIDAR topography grid of Resurrection Bay, and an initial thickness of slide material based on pre- and post-earthquake bathymetry difference maps. For simulation of tectonic tsunami runup, we derive the 1964 coseismic deformations from detailed slip distribution in the rupture area, and use them as an initial condition for propagation of the tectonic tsunami. The numerical model employs nonlinear shallow water equations formulated for depth-averaged water fluxes, and calculates a temporal position of the shoreline using a free-surface moving boundary algorithm. We find that the calculated tsunami runup in Seward caused first by local submarine landslide-generated waves, and later by a tectonic tsunami, is in good agreement with observations of the inundation zone. The analysis of inundation caused by two different tsunami sources improves our understanding of their relative contributions, and supports tsunami risk mitigation in south-central Alaska. The record of the 1964 earthquake, tsunami, and submarine landslides, combined with the high-resolution topography and bathymetry of Resurrection Bay make it an ideal location for studying tectonic tsunamis in coastal regions susceptible to underwater landslides.

Alaska↗

USGS2011

Book review: Extreme ocean waves

‘‘Extreme Ocean Waves’’ is a collection of ten papers edited by Efim Pelinovsky and Christian Kharif that followed the April 2007 meeting of the General Assembly of the European Geosciences Union. A note on terminology: extreme waves in this volume broadly encompass different types of waves, includ- ing deep-water and shallow-water rogue waves (alternatively termed freak waves), storm surges from cyclones, and internal waves. Other types of waves such as tsunamis or rissaga (meteotsunamis) are not discussed in this volume. It is generally implied that ‘‘extreme’’ has a statistical connotation relative to the average or significant wave height specific to each type of wave. Throughout the book, in fact, the reader will find a combination of theoretical and statistical/ empirical treatment necessary for the complete examination of this subject. In the introduction, the editors underscore the importance of studying extreme waves, documenting several dramatic instances of damaging extreme waves that occurred in 2007. Review info: Extreme Ocean Waves . By E. Pelinovsky and C. Kharif (eds), 2008. ISBN: 978-1402083136, xiii, 196 pp.

Pure and Applied Geophysics↗

USGS2011

Book review: Nonlinear ocean waves and the inverse scattering transform

Nonlinear Ocean Waves and the Inverse Scattering Transform is a comprehensive examination of ocean waves built upon the theory of nonlinear Fourier analysis. The renowned author, Alfred R. Osborne, is perhaps best known for the discovery of internal solitons in the Andaman Sea during the 1970s. In this book, he provides an extensive treatment of nonlinear water waves based on a nonlinear spectral theory known as the inverse scattering transform. The writing is exceptional throughout the book, which is particularly useful in explaining some of the more difficult mathematical concepts. Review info: N onlinear Ocean Waves and the Inverse Scattering Transform. B y Alfred R. Osborne, 2010. ISBN: 978-125286299, 917 pp.

Pure and Applied Geophysics↗

USGS2011

Nearshore Tsunami Inundation Model Validation: Toward Sediment Transport Applications

Model predictions from a numerical model, Delft3D, based on the nonlinear shallow water equations are compared with analytical results and laboratory observations from seven tsunami-like benchmark experiments, and with field observations from the 26 December 2004 Indian Ocean tsunami. The model accurately predicts the magnitude and timing of the measured water levels and flow velocities, as well as the magnitude of the maximum inundation distance and run-up, for both breaking and non-breaking waves. The shock-capturing numerical scheme employed describes well the total decrease in wave height due to breaking, but does not reproduce the observed shoaling near the break point. The maximum water levels observed onshore near Kuala Meurisi, Sumatra, following the 26 December 2004 tsunami are well predicted given the uncertainty in the model setup. The good agreement between the model predictions and the analytical results and observations demonstrates that the numerical solution and wetting and drying methods employed are appropriate for modeling tsunami inundation for breaking and non-breaking long waves. Extension of the model to include sediment transport may be appropriate for long, non-breaking tsunami waves. Using available sediment transport formulations, the sediment deposit thickness at Kuala Meurisi is predicted generally within a factor of 2.

Pure and Applied Geophysics↗

USGS2011

The role of adsorbed water on the friction of a layer of submicron particles

Anomalously low values of friction observed in layers of submicron particles deformed in simple shear at high slip velocities are explained as the consequence of a one nanometer thick layer of water adsorbed on the particles. The observed transition from normal friction with an apparent coefficient near μ = 0.6 at low slip speeds to a coefficient near μ = 0.3 at higher slip speeds is attributed to competition between the time required to extrude the water layer from between neighboring particles in a force chain and the average lifetime of the chain. At low slip speeds the time required for extrusion is less than the average lifetime of a chain so the particles make contact and lock. As slip speed increases, the average lifetime of a chain decreases until it is less than the extrusion time and the particles in a force chain never come into direct contact. If the adsorbed water layer enables the otherwise rough particles to rotate, the coefficient of friction will drop to μ = 0.3, appropriate for rotating spheres. At the highest slip speeds particle temperatures rise above 100°C, the water layer vaporizes, the particles contact and lock, and the coefficient of friction rises to μ = 0.6. The observed onset of weakening at slip speeds near 0.001 m/s is consistent with the measured viscosity of a 1 nm thick layer of adsorbed water, with a minimum particle radius of approximately 20 nm, and with reasonable assumptions about the distribution of force chains guided by experimental observation. The reduction of friction and the range of velocities over which it occurs decrease with increasing normal stress, as predicted by the model. Moreover, the analysis predicts that this high-speed weakening mechanism should operate only for particles with radii smaller than approximately 1 μm. For larger particles the slip speed required for weakening is so large that frictional heating will evaporate the adsorbed water and weakening will not occur.

Pure and Applied Geophysics↗

USGS2011

Earthquake rupture at focal depth, part II: mechanics of the 2004 M2.2 earthquake along the Pretorius Fault, TauTona Mine, South Africa

We analyze here the rupture mechanics of the 2004, M2.2 earthquake based on our observations and measurements at focal depth (Part I). This event ruptured the Archean Pretorius fault that has been inactive for at least 2 Ga, and was reactivated due to mining operations down to a depth of 3.6 km depth. Thus, it was expected that the Pretorius fault zone will fail similarly to an intact rock body independently of its ancient healed structure. Our analysis reveals a few puzzling features of the M2.2 rupture-zone: (1) the earthquake ruptured four, non-parallel, cataclasite bearing segments of the ancient Pretorius fault-zone; (2) slip occurred almost exclusively along the cataclasite-host rock contacts of the slipping segments; (3) the local in-situ stress field is not favorable to slip along any of these four segments; and (4) the Archean cataclasite is pervasively sintered and cemented to become brittle and strong. To resolve these observations, we conducted rock mechanics experiments on the fault-rocks and host-rocks and found a strong mechanical contrast between the quartzitic cataclasite zones, with elastic-brittle rheology, and the host quartzites, with damage, elastic–plastic rheology. The finite-element modeling of a heterogeneous fault-zone with the measured mechanical contrast indicates that the slip is likely to reactivate the ancient cataclasite-bearing segments, as observed, due to the strong mechanical contrast between the cataclasite and the host quartzitic rock.

Tautona Mine↗

USGS2012

Near-field tsunami edge waves and complex earthquake rupture

The effect of distributed coseismic slip on progressive, near-field edge waves is examined for continental shelf tsunamis. Detailed observations of edge waves are difficult to separate from the other tsunami phases that are observed on tide gauge records. In this study, analytic methods are used to compute tsunami edge waves distributed over a finite number of modes and for uniformly sloping bathymetry. Coseismic displacements from static elastic theory are introduced as initial conditions in calculating the evolution of progressive edge-waves. Both simple crack representations (constant stress drop) and stochastic slip models (heterogeneous stress drop) are tested on a fault with geometry similar to that of the M w = 8.8 2010 Chile earthquake. Crack-like ruptures that are beneath or that span the shoreline result in similar longshore patterns of maximum edge-wave amplitude. Ruptures located farther offshore result in reduced edge-wave excitation, consistent with previous studies. Introduction of stress-drop heterogeneity by way of stochastic slip models results in significantly more variability in longshore edge-wave patterns compared to crack-like ruptures for the same offshore source position. In some cases, regions of high slip that are spatially distinct will yield sub-events, in terms of tsunami generation. Constructive interference of both non-trapped and trapped waves can yield significantly larger tsunamis than those that produced by simple earthquake characterizations.

Pure and Applied Geophysics↗

USGS2013

What do data used to develop ground-motion prediction equations tell us about motions near faults?

A large database of ground motions from shallow earthquakes occurring in active tectonic regions around the world, recently developed in the Pacific Earthquake Engineering Center’s NGA-West2 project, has been used to investigate what such a database can say about the properties and processes of crustal fault zones. There are a relatively small number of near-rupture records, implying that few recordings in the database are within crustal fault zones, but the records that do exist emphasize the complexity of ground-motion amplitudes and polarization close to individual faults. On average over the whole data set, however, the scaling of ground motions with magnitude at a fixed distance, and the distance dependence of the ground motions, seem to be largely consistent with simple seismological models of source scaling, path propagation effects, and local site amplification. The data show that ground motions close to large faults, as measured by elastic response spectra, tend to saturate and become essentially constant for short periods. This saturation seems to be primarily a geometrical effect, due to the increasing size of the rupture surface with magnitude, and not due to a breakdown in self similarity.

Pure and Applied Geophysics↗

USGS2014

Laboratory generated M -6 earthquakes

We consider whether mm-scale earthquake-like seismic events generated in laboratory experiments are consistent with our understanding of the physics of larger earthquakes. This work focuses on a population of 48 very small shocks that are foreshocks and aftershocks of stick–slip events occurring on a 2.0 m by 0.4 m simulated strike-slip fault cut through a large granite sample. Unlike the larger stick–slip events that rupture the entirety of the simulated fault, the small foreshocks and aftershocks are contained events whose properties are controlled by the rigidity of the surrounding granite blocks rather than characteristics of the experimental apparatus. The large size of the experimental apparatus, high fidelity sensors, rigorous treatment of wave propagation effects, and in situ system calibration separates this study from traditional acoustic emission analyses and allows these sources to be studied with as much rigor as larger natural earthquakes. The tiny events have short (3–6 μs) rise times and are well modeled by simple double couple focal mechanisms that are consistent with left-lateral slip occurring on a mm-scale patch of the precut fault surface. The repeatability of the experiments indicates that they are the result of frictional processes on the simulated fault surface rather than grain crushing or fracture of fresh rock. Our waveform analysis shows no significant differences (other than size) between the M -7 to M -5.5 earthquakes reported here and larger natural earthquakes. Their source characteristics such as stress drop (1–10 MPa) appear to be entirely consistent with earthquake scaling laws derived for larger earthquakes.

Pure and Applied Geophysics↗

USGS2014

Tsunami forecast by joint inversion of real-time tsunami waveforms and seismic of GPS data: application to the Tohoku 2011 tsunami

Correctly characterizing tsunami source generation is the most critical component of modern tsunami forecasting. Although difficult to quantify directly, a tsunami source can be modeled via different methods using a variety of measurements from deep-ocean tsunameters, seismometers, GPS, and other advanced instruments, some of which in or near real time. Here we assess the performance of different source models for the destructive 11 March 2011 Japan tsunami using model–data comparison for the generation, propagation, and inundation in the near field of Japan. This comparative study of tsunami source models addresses the advantages and limitations of different real-time measurements with potential use in early tsunami warning in the near and far field. The study highlights the critical role of deep-ocean tsunami measurements and rapid validation of the approximate tsunami source for high-quality forecasting. We show that these tsunami measurements are compatible with other real-time geodetic data, and may provide more insightful understanding of tsunami generation from earthquakes, as well as from nonseismic processes such as submarine landslide failures.

Honshu↗

USGS2014

Evolution of wear and friction along experimental faults

We investigate the evolution of wear and friction along experimental faults composed of solid rock blocks. This evolution is analyzed through shear experiments along five rock types, and the experiments were conducted in a rotary apparatus at slip velocities of 0.002–0.97 m/s, slip distances from a few millimeters to tens of meters, and normal stress of 0.25–6.9 MPa. The wear and friction measurements and fault surface observations revealed three evolution phases: A) An initial stage (slip distances <50 mm) of wear by failure of isolated asperities associated with roughening of the fault surface; B) a running-in stage of slip distances of 1–3 m with intense wear-rate, failure of many asperities, and simultaneous reduction of the friction coefficient and wear-rate; and C) a steady-state stage that initiates when the fault surface is covered by a gouge layer, and during which both wear-rate and friction coefficient maintain quasi-constant, low levels. While these evolution stages are clearly recognizable for experimental faults made from bare rock blocks, our analysis suggests that natural faults “bypass” the first two stages and slip at gouge-controlled steady-state conditions.

Pure and Applied Geophysics↗

USGS2014

Book review: Three great tsunamis: Lisbon (1755), Sumatra-Andaman (2004), and Japan (2011)

“Three Great Tsunamis: Lisbon (1755), Sumatra–Andaman (2004), and Japan (2011)” is published in Springer’s new series SpringerBriefs. According to Springer’s website, the SpringBriefs volumes are intended to provide “concise summaries of cutting-edge research and practical applications across a wide spectrum of fields”. Among the several categories considered for SpringerBriefs are in-depth case studies, for which this volume is most closely aligned. This 89-page volume comprises six chapters. Chapter 1 describes the fundamentals of tsunamis at a basic level, including a selective list of major historical tsunamis. Chapter 2 describes the 1755 Lisbon earthquake and tsunami, although no tectonic background is given for this event, as it is for the other two events considered in this volume. Included in this chapter, however, is a brief but fascinating account of how this earthquake affected scientific and philosophical thought during the Age of Enlightenment. The authors’ description of the 2004 Sumatra–Andaman earthquake and tsunami in Chapter 3 includes instrumental and post-tsunami measurements of the tsunami, ground shaking and subsidence, and other effects such as mud volcanoes. Primarily, tsunami observations are described for the 2011 Japan event, including a brief summary of the nuclear disaster in the Fukushima region. The fifth chapter is on tsunami warning systems, including new and emerging systems around the world in the twenty-first century. The final chapter provides perspectives on why these events were so destructive and how future hazards may be mitigated. Review info: Three Great Tsunamis: Lisbon (1755), Sumatra–Andaman (2004), and Japan (2011). By Harsh K. Gupta and Vineet K. Gahalaut, 2013. ISBN: 978-9400765757

Pure and Applied Geophysics↗

USGS2014

Joint inversion of seismic and magnetotelluric data in the Parkfield Region of California using the normalized cross-gradient constraint

We present jointly inverted models of P-wave velocity (Vp) and electrical resistivity for a two-dimensional profile centered on the San Andreas Fault Observatory at Depth (SAFOD). Significant structural similarity between main features of the separately inverted Vp and resistivity models is exploited by carrying out a joint inversion of the two datasets using the normalized cross-gradient constraint. This constraint favors structurally similar Vp and resistivity images that adequately fit the seismic and magnetotelluric (MT) datasets. The new inversion code, tomoDDMT, merges the seismic inversion code tomoDD and the forward modeling and sensitivity kernel subroutines of the MT inversion code OCCAM2DMT. TomoDDMT is tested on a synthetic dataset and demonstrates the code’s ability to more accurately resolve features of the input synthetic structure relative to the separately inverted resistivity and velocity models. Using tomoDDMT, we are able to resolve a number of key issues raised during drilling at SAFOD. We are able to infer the distribution of several geologic units including the Salinian granitoids, the Great Valley sequence, and the Franciscan Formation. The distribution and transport of fluids at both shallow and great depths is also examined. Low values of velocity/resistivity attributed to a feature known as the Eastern Conductor (EC) can be explained in two ways: the EC is a brine-filled, high porosity region, or this region is composed largely of clay-rich shales of the Franciscan. The Eastern Wall, which lies immediately adjacent to the EC, is unlikely to be a fluid pathway into the San Andreas Fault’s seismogenic zone due to its observed higher resistivity and velocity values.

California↗

USGS2016

Reconstruction of far-field tsunami amplitude distributions from earthquake sources

The probability distribution of far-field tsunami amplitudes is explained in relation to the distribution of seismic moment at subduction zones. Tsunami amplitude distributions at tide gauge stations follow a similar functional form, well described by a tapered Pareto distribution that is parameterized by a power-law exponent and a corner amplitude. Distribution parameters are first established for eight tide gauge stations in the Pacific, using maximum likelihood estimation. A procedure is then developed to reconstruct the tsunami amplitude distribution that consists of four steps: (1) define the distribution of seismic moment at subduction zones; (2) establish a source-station scaling relation from regression analysis; (3) transform the seismic moment distribution to a tsunami amplitude distribution for each subduction zone; and (4) mix the transformed distribution for all subduction zones to an aggregate tsunami amplitude distribution specific to the tide gauge station. The tsunami amplitude distribution is adequately reconstructed for four tide gauge stations using globally constant seismic moment distribution parameters established in previous studies. In comparisons to empirical tsunami amplitude distributions from maximum likelihood estimation, the reconstructed distributions consistently exhibit higher corner amplitude values, implying that in most cases, the empirical catalogs are too short to include the largest amplitudes. Because the reconstructed distribution is based on a catalog of earthquakes that is much larger than the tsunami catalog, it is less susceptible to the effects of record-breaking events and more indicative of the actual distribution of tsunami amplitudes.

Hawaii↗

USGS2016

A study of the 2015 Mw 8.3 Illapel earthquake and tsunami: Numerical and analytical approaches

The September 16, 2015 Illapel, Chile earthquake triggered a large tsunami, causing both economic losses and fatalities. To study the coastal effects of this earthquake, and to understand how such hazards might be accurately modeled in the future, different finite fault models of the Illapel rupture are used to define the initial condition for tsunami simulation. The numerical code Non-hydrostatic Evolution of Ocean WAVEs (NEOWAVE) is employed to model the tsunami evolution through the Pacific Ocean. Because only a short time is available for emergency response, and since the earthquake and tsunami sources are close to the coast, gaining a rapid understanding of the near-field run-up behavior is highly relevant to Chile. Therefore, an analytical solution of the 2 ? 1 D shallow water wave equations is considered. With this solution, we show that we can quickly estimate the run-up distribution along the coastline, to first order. After the earthquake and tsunami, field observations were measured in the surrounded coastal region, where the tsunami resulted in significant run-up. First, we compare the analytical and numerical solutions to test the accuracy of the analytical approach and the field observations, implying the analytic approach can accurately model tsunami run-up after an earthquake, without sacrificing the time necessary for a full numerical inversion. Then, we compare both with field run-up measurements. We observe the consistency between the two approaches. To complete the analysis, a tsunami source inversion is performed using run-up field measurements only. These inversion results are compared with seismic models, and are shown to capture the broad-scale details of those models, without the necessity of the detailed data sets they invert.

Illapel↗

USGS2016

Source characterization and tsunami modeling of submarine landslides along the Yucatán Shelf/Campeche Escarpment, southern Gulf of Mexico

Submarine landslides occurring along the margins of the Gulf of Mexico (GOM) represent a low-likelihood, but potentially damaging source of tsunamis. New multibeam bathymetry coverage reveals that mass wasting is pervasive along the Yucatán Shelf edge with several large composite landslides possibly removing as much as 70 km 3 of the Cenozoic sedimentary section in a single event. Using GIS-based analysis, the dimensions of six landslides from the central and northern sections of the Yucatán Shelf/Campeche Escarpment were determined and used as input for preliminary tsunami generation and propagation models. Tsunami modeling is performed to compare the propagation characteristics and distribution of maximum amplitudes throughout the GOM among the different landslide scenarios. Various factors such as landslide geometry, location along the Yucatán Shelf/Campeche Escarpment, and refraction during propagation result in significant variations in the affected part of the Mexican and US Gulf Coasts. In all cases, however, tsunami amplitudes are greatest along the northern Yucatán Peninsula.

Pure and Applied Geophysics↗

USGS2016

Introduction to “Global tsunami science: Past and future, Volume I”

Twenty-five papers on the study of tsunamis are included in Volume I of the PAGEOPH topical issue “Global Tsunami Science: Past and Future”. Six papers examine various aspects of tsunami probability and uncertainty analysis related to hazard assessment. Three papers relate to deterministic hazard and risk assessment. Five more papers present new methods for tsunami warning and detection. Six papers describe new methods for modeling tsunami hydrodynamics. Two papers investigate tsunamis generated by non-seismic sources: landslides and meteorological disturbances. The final three papers describe important case studies of recent and historical events. Collectively, this volume highlights contemporary trends in global tsunami research, both fundamental and applied toward hazard assessment and mitigation.

Pure and Applied Geophysics PAGEOPH↗

USGS2016

Solving for source parameters using nested array data: A case study from the Canterbury, New Zealand earthquake sequence

The seismic spectrum can be constructed by assuming a Brune spectral model and estimating the parameters of seismic moment ( M 0 ), corner frequency ( f c ), and high-frequency site attenuation ( κ ). Using seismic data collected during the 2010–2011 Canterbury, New Zealand, earthquake sequence, we apply the non-linear least-squares Gauss–Newton method, a deterministic downhill optimization technique, to simultaneously determine the M 0 , f c , and κ for each event-station pair. We fit the Brune spectral acceleration model to Fourier-transformed S-wave records following application of path and site corrections to the data. For each event, we solve for a single M 0 and f c , while any remaining residual kappa, κ r "> κ r κr , is allowed to differ per station record to reflect varying high-frequency falloff due to path and site attenuation. We use a parametric forward modeling method, calculating initial M 0 and f c values from the local GNS New Zealand catalog M w, GNS magnitudes and measuring an initial κ r "> κ r κr using an automated high-frequency linear regression method. Final solutions for M 0 , f c , and κ r "> κ r κr are iteratively computed through minimization of the residual function, and the Brune model stress drop is then calculated from the final, best-fit f c . We perform the spectral fitting routine on nested array seismic data that include the permanent GeoNet accelerometer network as well as a dense network of nearly 200 Quake Catcher Network (QCN) MEMs accelerometers, analyzing over 180 aftershocks M w,GNS ≥ 3.5 that occurred from 9 September 2010 to 31 July 2011. QCN stations were hosted by public volunteers and served to fill spatial gaps between existing GeoNet stations. Moment magnitudes determined using the spectral fitting procedure ( M w,SF ) range from 3.5 to 5.7 and agree well with M w,GNS , with a median difference of 0.09 and 0.17 for GeoNet and QCN records, respectively, and 0.11 when data from both networks are combined. The majority of events are calculated to have stress drops between 1.7 and 13 MPa (20th and 80th percentile, correspondingly) for the combined networks. The overall median stress drop for the combined networks is 3.2 MPa, which is similar to median stress drops previously reported for the Canterbury sequence. We do not observe a correlation between stress drop and depth for this region, nor a relationship between stress drop and magnitude over the catalog considered. Lateral spatial patterns in stress drop, such as a cluster of aftershocks near the eastern extent of the Greendale fault with higher stress drops and lower stress drops for aftershocks of the 2011 M w,GNS 6.2 Christchurch mainshock, are found to be in agreement with previous reports. As stress drop is arguably a method-dependent calculation and subject to high spatial variability, our results using the parametric Gauss–Newton algorithm strengthen conclusions that the Canterbury sequence has stress drops that are more similar to those found in intraplate regions, with overall higher stress drops that are typically observed in tectonically active areas.

Christchurch↗