USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geophysics Journal International”

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

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

At least 181 records · Page 10Linked to original sources
USGS2007

Spectral element modelling of fault-plane reflections arising from fluid pressure distributions

The presence of fault-plane reflections in seismic images, besides indicating the locations of faults, offers a possible source of information on the properties of these poorly understood zones. To better understand the physical mechanism giving rise to fault-plane reflections in compacting sedimentary basins, we numerically model the full elastic wavefield via the spectral element method (SEM) for several different fault models. Using well log data from the South Eugene Island field, offshore Louisiana, we derive empirical relationships between the elastic parameters (e.g. P-wave velocity and density) and the effective-stress along both normal compaction and unloading paths. These empirical relationships guide the numerical modelling and allow the investigation of how differences in fluid pressure modify the elastic wavefield. We choose to simulate the elastic wave equation via SEM since irregular model geometries can be accommodated and slip boundary conditions at an interface, such as a fault or fracture, are implemented naturally. The method we employ for including a slip interface retains the desirable qualities of SEM in that it is explicit in time and, therefore, does not require the inversion of a large matrix. We performa complete numerical study by forward modelling seismic shot gathers over a faulted earth model using SEM followed by seismic processing of the simulated data. With this procedure, we construct post-stack time-migrated images of the kind that are routinely interpreted in the seismic exploration industry. We dip filter the seismic images to highlight the fault-plane reflections prior to making amplitude maps along the fault plane. With these amplitude maps, we compare the reflectivity from the different fault models to diagnose which physical mechanism contributes most to observed fault reflectivity. To lend physical meaning to the properties of a locally weak fault zone characterized as a slip interface, we propose an equivalent-layer model under the assumption of weak scattering. This allows us to use the empirical relationships between density, velocity and effective stress from the South Eugene Island field to relate a slip interface to an amount of excess pore-pressure in a fault zone. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗
USGS2007

Lithology-derived structure classification from the joint interpretation of magnetotelluric and seismic models

Magnetotelluric and seismic methods provide complementary information about the resistivity and velocity structure of the subsurface on similar scales and resolutions. No global relation, however, exists between these parameters, and correlations are often valid for only a limited target area. Independently derived inverse models from these methods can be combined using a classification approach to map geologic structure. The method employed is based solely on the statistical correlation of physical properties in a joint parameter space and is independent of theoretical or empirical relations linking electrical and seismic parameters. Regions of high correlation (classes) between resistivity and velocity can in turn be mapped back and re-examined in depth section. The spatial distribution of these classes, and the boundaries between them, provide structural information not evident in the individual models. This method is applied to a 10 km long profile crossing the Dead Sea Transform in Jordan. Several prominent classes are identified with specific lithologies in accordance with local geology. An abrupt change in lithology across the fault, together with vertical uplift of the basement suggest the fault is sub-vertical within the upper crust. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗
USGS2007

Upper-crustal structure beneath the strait of Georgia, Southwest British Columbia

We present a new three-dimensional (3-D) P -wave velocity model for the upper-crustal structure beneath the Strait of Georgia, southwestern British Columbia based on non-linear tomographic inversion of wide-angle seismic refraction data. Our study, part of the Georgia Basin Geohazards Initiative (GBGI) is primarily aimed at mapping the depth of the Cenozoic sedimentary basin and delineating the near-surface crustal faults associated with recent seismic activities (e.g. M = 4.6 in 1997 and M = 5.0 in 1975) in the region. Joint inversion of first-arrival traveltimes from the 1998 Seismic Hazards Investigation in Puget Sound (SHIPS) and the 2002 Georgia Basin experiment provides a high-resolution velocity model of the subsurface to a depth of ∼7 km. In the southcentral Georgia Basin, sedimentary rocks of the Cretaceous Nanaimo Group and early Tertiary rocks have seismic velocities between 3.0 and 5.5 km s −1 . The basin thickness increases from north to south with a maximum thickness of 7 (±1) km (depth to velocities of 5.5 km s −1 ) at the southeast end of the strait. The underlying basement rocks, probably representing the Wrangellia terrane, have velocities of 5.5–6.5 km s −1 with considerable lateral variation. Our tomographic model reveals that the Strait of Georgia is underlain by a fault-bounded block within the central Georgia Basin. It also shows a correlation between microearthquakes and areas of rapid change in basin thickness. The 1997/1975 earthquakes are located near a northeast-trending hinge line where the thicknesses of sedimentary rocks increase rapidly to the southeast. Given its association with instrumentally recorded, moderate sized earthquakes, we infer that the hinge region is cored by an active fault that we informally name the Gabriola Island fault. A northwest-trending, southwest dipping velocity discontinuity along the eastern side of Vancouver Island correlates spatially with the surface expression of the Outer Island fault. The Outer Island fault as mapped in our seismic tomography model is a thrust fault that projects directly into the Lummi Island fault, suggesting that they are related structures forming a fault system that is continuous for nearly 90 km. Together, these inferred thrust faults may account for at least a portion of the basement uplift at the San Juan Islands.

British Columbia↗
USGS2007

Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity

The GRACE satellite mission has been measuring the Earth's gravity field and its temporal variations since 2002 April. Although these variations are mainly due to mass transfer within the geofluid envelops, they also result from mass displacements associated with phenomena including glacial isostatic adjustment and earthquakes. However, these last contributions are difficult to isolate because of the presence of noise and of geofluid signals, and because of GRACE's coarse spatial resolution (>400 km half-wavelength). In this paper, we show that a wavelet analysis on the sphere helps to retrieve earthquake signatures from GRACE geoid products. Using a wavelet analysis of GRACE geoids products, we show that the geoid variations caused by the 2004 December ( M w = 9.2) and 2005 March ( M w = 8.7) Sumatra earthquakes can be detected. At GRACE resolution, the 2004 December earthquake produced a strong coseismic decrease of the gravity field in the Andaman Sea, followed by relaxation in the area affected by both the Andaman 2004 and the Nias 2005 earthquakes. We find two characteristic timescales for the relaxation, with a fast variation occurring in the vicinity of the Central Andaman ridge. We discuss our coseismic observations in terms of density changes of crustal and upper-mantle rocks, and of the vertical displacements in the Andaman Sea. We interpret the post-seismic signal in terms of the viscoelastic response of the Earth's mantle. The transient component of the relaxation may indicate the presence of hot, viscous material beneath the active Central Andaman Basin.

Geophysical Journal International↗
USGS2007

Finite-frequency traveltime tomography of San Francisco Bay region crustal velocity structure

Seismic velocity structure of the San Francisco Bay region crust is derived using measurements of finite-frequency traveltimes. A total of 57 801 relative traveltimes are measured by cross-correlation over the frequency range 0.5–1.5 Hz. From these are derived 4862 ‘summary’ traveltimes, which are used to derive 3-D P -wave velocity structure over a 341 × 140 km 2 area from the surface to 25 km depth. The seismic tomography is based on sensitivity kernels calculated on a spherically symmetric reference model. Robust elements of the derived P -wave velocity structure are: a pronounced velocity contrast across the San Andreas fault in the south Bay region (west side faster); a moderate velocity contrast across the Hayward fault (west side faster); moderately low velocity crust around the Quien Sabe volcanic field and the Sacramento River delta; very low velocity crust around Lake Berryessa. These features are generally explicable with surface rock types being extrapolated to depth ∼10 km in the upper crust. Generally high mid-lower crust velocity and high inferred Poisson's ratio suggest a mafic lower crust.

California↗
USGS2007

A crustal seismic velocity model for the UK, Ireland and surrounding seas

A regional model of the 3-D variation in seismic P -wave velocity structure in the crust of NW Europe has been compiled from wide-angle reflection/refraction profiles. Along each 2-D profile a velocity–depth function has been digitised at 5 km intervals. These 1-D velocity functions were mapped into three dimensions using ordinary kriging with weights determined to minimise the difference between digitised and interpolated values. An analysis of variograms of the digitised data suggested a radial isotropic weighting scheme was most appropriate. Horizontal dimensions of the model cells are optimised at 40 × 40 km and the vertical dimension at 1 km. The resulting model provides a higher resolution image of the 3-D variation in seismic velocity structure of the UK, Ireland and surrounding areas than existing models. The construction of the model through kriging allows the uncertainty in the velocity structure to be assessed. This uncertainty indicates the high density of data required to confidently interpolate the crustal velocity structure, and shows that for this region the velocity is poorly constrained for large areas away from the input data.

Geophysical Journal International↗
USGS2008

Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M = 9.2 Sumatra earthquake

The 2004 M=9.2 Sumatra-Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ???1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress-strain relation; an aspherical perturbation in viscoelastic structure; a 'static'mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra-Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation. ?? Journal compilation ?? 2008 RAS.

Geophysical Journal International↗
USGS2008

The role of shear and tensile failure in dynamically triggered landslides

Dynamic stresses generated by earthquakes can trigger landslides. Current methods of landslide analysis such as pseudo-static analysis and Newmark's method focus on the effects of earthquake accelerations on the landslide mass to characterize dynamic landslide behaviour. One limitation of these methods is their use Mohr-Coulomb failure criteria, which only accounts for shear failure, but the role of tensile failure is not accounted for. We develop a limit-equilibrium model to investigate the dynamic stresses generated by a given ground motion due to a plane wave and use this model to assess the role of shear and tensile failure in the initiation of slope instability. We do so by incorporating a modified Griffith failure envelope, which combines shear and tensile failure into a single criterion. Tests of dynamic stresses in both homogeneous and layered slopes demonstrate that two modes of failure exist, tensile failure in the uppermost meters of a slope and shear failure at greater depth. Further, we derive equations that express the dynamic stress in the near-surface in the acceleration measured at the surface. These equations are used to approximately define the depth range for each mechanism of failure. The depths at which these failure mechanisms occur suggest that shear and tensile failure might collaborate in generating slope failure. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗
USGS2008

Effects of topography and crustal heterogeneities on the source estimation of LP event at Kilauea volcano

The main goal of this study is to improve the modelling of the source mechanism associated with the generation of long period (LP) signals in volcanic areas. Our intent is to evaluate the effects that detailed structural features of the volcanic models play in the generation of LP signal and the consequent retrieval of LP source characteristics. In particular, effects associated with the presence of topography and crustal heterogeneities are here studied in detail. We focus our study on a LP event observed at Kilauea volcano, Hawaii, in 2001 May. A detailed analysis of this event and its source modelling is accompanied by a set of synthetic tests, which aim to evaluate the effects of topography and the presence of low velocity shallow layers in the source region. The forward problem of Green's function generation is solved numerically following a pseudo-spectral approach, assuming different 3-D models. The inversion is done in the frequency domain and the resulting source mechanism is represented by the sum of two time-dependent terms: a full moment tensor and a single force. Synthetic tests show how characteristic velocity structures, associated with shallow sources, may be partially responsible for the generation of the observed long-lasting ringing waveforms. When applying the inversion technique to Kilauea LP data set, inversions carried out for different crustal models led to very similar source geometries, indicating a subhorizontal cracks. On the other hand, the source time function and its duration are significantly different for different models. These results support the indication of a strong influence of crustal layering on the generation of the LP signal, while the assumption of homogeneous velocity model may bring to misleading results.

Geophysical Journal International↗
USGS2008

A deep crustal fluid channel into the San Andreas Fault system near Parkfield, California

Magnetotelluric (MT) data from 66 sites along a 45-km-long profile across the San Andreas Fault (SAF) were inverted to obtain the 2-D electrical resistivity structure of the crust near the San Andreas Fault Observatory at Depth (SAFOD). The most intriguing feature of the resistivity model is a steeply dipping upper crustal high-conductivity zone flanking the seismically defined SAF to the NE, that widens into the lower crust and appears to be connected to a broad conductivity anomaly in the upper mantle. Hypothesis tests of the inversion model suggest that upper and lower crustal and upper-mantle anomalies may be interconnected. We speculate that the high conductivities are caused by fluids and may represent a deep-rooted channel for crustal and/or mantle fluid ascent. Based on the chemical analysis of well waters, it was previously suggested that fluids can enter the brittle regime of the SAF system from the lower crust and mantle. At high pressures, these fluids can contribute to fault-weakening at seismogenic depths. These geochemical studies predicted the existence of a deep fluid source and a permeable pathway through the crust. Our resistivity model images a conductive pathway, which penetrates the entire crust, in agreement with the geochemical interpretation. However, the resistivity model also shows that the upper crustal branch of the high-conductivity zone is located NE of the seismically defined SAF, suggesting that the SAF does not itself act as a major fluid pathway. This interpretation is supported by both, the location of the upper crustal high-conductivity zone and recent studies within the SAFOD main hole, which indicate that pore pressures within the core of the SAF zone are not anomalously high, that mantle-derived fluids are minor constituents to the fault-zone fluid composition and that both the volume of mantle fluids and the fluid pressure increase to the NE of the SAF. We further infer from the MT model that the resistive Salinian block basement to the SW of the SAFOD represents an isolated body, being 5-8km wide and reaching to depths >7km, in agreement with aeromagnetic data. This body is separated from a massive block of Salinian crust farther to the SW. The NE terminus of resistive Salinian crust has a spatial relationship with a near-vertical zone of increased seismic reflectivity ???15km SW of the SAF and likely represents a deep-reaching fault zone. ?? 2008 The Authors Journal compilation ?? 2008 RAS.

Geophysical Journal International↗
USGS2009

The postseismic response to the 2002 M 7.9 Denali Fault earthquake: Constraints from InSAR 2003-2005

InSAR is particularly sensitive to vertical displacements, which can be important in distinguishing between mechanisms responsible for the postseismic response to large earthquakes (afterslip, viscoelastic relaxation). We produce maps of the surface displacements resulting from the postseismic response to the 2002 Denali Fault earthquake, using data from the Canadian Radarsat-1 satellite from the periods summer 2003, summer 2004 and summer 2005. A peak-to-trough signal of amplitude 4 cm in the satellite line of sight was observed between summer 2003 and summer 2004. By the period between summer 2004 and summer 2005, the displacement rate had dropped below the threshold required for observation with InSAR over a single year. The InSAR observations show that the principal postseismic relaxation process acted at a depth of ∼50 km, equivalent to the top of the mantle. However, the observations are still incapable of distinguishing between distributed (viscoelastic relaxation) and localized (afterslip) deformation. The imposed coseismic stresses are highest in the lower crust and, assuming a Maxwell rheology, a viscosity ratio of at least 5 between lower crust and upper mantle is required to explain the contrast in behaviour. The lowest misfits are produced by mixed models of viscoelastic relaxation in the mantle and shallow afterslip in the upper crust. Profiles perpendicular to the fault show significant asymmetry, which is consistent with differences in rheological structure across the fault.

Geophysical Journal International↗
USGS2009

Developing framework to constrain the geometry of the seismic rupture plane on subduction interfaces a priori - A probabilistic approach

A key step in many earthquake source inversions requires knowledge of the geometry of the fault surface on which the earthquake occurred. Our knowledge of this surface is often uncertain, however, and as a result fault geometry misinterpretation can map into significant error in the final temporal and spatial slip patterns of these inversions. Relying solely on an initial hypocentre and CMT mechanism can be problematic when establishing rupture characteristics needed for rapid tsunami and ground shaking estimates. Here, we attempt to improve the quality of fast finite-fault inversion results by combining several independent and complementary data sets to more accurately constrain the geometry of the seismic rupture plane of subducting slabs. Unlike previous analyses aimed at defining the general form of the plate interface, we require mechanisms and locations of the seismicity considered in our inversions to be consistent with their occurrence on the plate interface, by limiting events to those with well-constrained depths and with CMT solutions indicative of shallow-dip thrust faulting. We construct probability density functions about each location based on formal assumptions of their depth uncertainty and use these constraints to solve for the ‘most-likely’ fault plane. Examples are shown for the trench in the source region of the M w 8.6 Southern Sumatra earthquake of March 2005, and for the Northern Chile Trench in the source region of the November 2007 Antofagasta earthquake. We also show examples using only the historic catalogues in regions without recent great earthquakes, such as the Japan and Kamchatka Trenches. In most cases, this method produces a fault plane that is more consistent with all of the data available than is the plane implied by the initial hypocentre and CMT mechanism. Using the aggregated data sets, we have developed an algorithm to rapidly determine more accurate initial fault plane geometries for source inversions of future earthquakes.

Geophysical Journal International↗
USGS2009

Validation of the rupture properties of the 2001 Kunlun, China (Ms = 8.1), earthquake from seismological and geological observations

We determine the finite-fault slip distribution of the 2001 Kunlun earthquake ( M s = 8.1) by inverting teleseismic waveforms, as constrained by geological and remote sensing field observations. The spatial slip distribution along the 400-km-long fault was divided into five segments in accordance with geological observations. Forward modelling of regional surface waves was performed to estimate the variation of the speed of rupture propagation during faulting. For our modelling, the regional 1-D velocity structure was carefully constructed for each of six regional seismic stations using three events with magnitudes of 5.1–5.4 distributed along the ruptured portion of the Kunlun fault. Our result shows that the average rupture velocity is about 3.6 km s −1 , consistent with teleseismic long period wave modelling. The initial rupture was almost purely strike-slip with a rupture velocity of 1.9 km s −1 , increasing to 3.5 km s −1 in the second fault segment, and reaching a rupture velocity of about 6 km s −1 in the third segment and the fourth segment, where the maximum surface offset, with a broad fault zone, was observed. The rupture velocity decelerated to a value of 3.3 km s −1 in the fifth and final segment. Coseismic slip on the fault was concentrated between the surface and a depth of about 10 km. We infer that significant variations in rupture velocity and the observed fault segmentation are indicative of variations in strength along the interface of the Kunlun fault, as well as variations in fault geometry.

Qinghai Province↗
USGS2009

Stress evolution following the 1999 Chi-Chi, Taiwan, earthquake: Consequences for afterslip, relaxation, aftershocks and departures from Omori decay

We explore how Coulomb stress transfer and viscoelastic relaxation control afterslip and aftershocks in a continental thrust fault system. The 1999 September 21 M w = 7.6 Chi-Chi shock is typical of continental ramp-d??collement systems throughout the world, and so inferences drawn from this uniquely well-recorded event may be widely applicable. First, we find that the spatial and depth distribution of aftershocks and their focal mechanisms are consistent with the calculated Coulomb stress changes imparted by the coseismic rupture. Some 61 per cent of the M ??? 2 aftershocks and 83 per cent of the M ??? 4 aftershocks lie in regions for which the Coulomb stress increased by ???0.1 bars, and there is a 11-12 per cent gain in the percentage of aftershocks nodal planes on which the shear stress increased over the pre-Chi Chi control period. Second, we find that afterslip occurred where the calculated coseismic stress increased on the fault ramp and d??collement, subject to the condition that friction is high on the ramp and low on the d??collement. Third, viscoelastic relaxation is evident from the fit of the post-seismic GPS data on the footwall. Fourth, we find that the rate of seismicity began to increase during the post-seismic period in an annulus extending east of the main rupture. The spatial extent of the seismicity annulus resembles the calculated ???0.05-bar Coulomb stress increase caused by viscoelastic relaxation and afterslip, and we find a 9-12 per cent gain in the percentage of focal mechanisms with >0.01-bar shear stress increases imparted by the post-seismic afterslip and relaxation in comparison to the control period. Thus, we argue that post-seismic stress changes can for the first time be shown to alter the production of aftershocks, as judged by their rate, spatial distribution, and focal mechanisms. ?? Journal compilation ?? 2009 RAS.

Geophysical Journal International↗
USGS2009

Offshore double-planed shallow seismic zone in the NE Japan forearc region revealed by sP depth phases recorded by regional networks

We detected the sP depth phase at small epicentral distances of about 150 km or more in the seismograms of shallow earthquakes in the NE Japan forearc region. The focal depths of 1078 M > 3 earthquakes that occurred from 2000 to 2006 were precisely determined using the time delay of the sP phase from the initial P-wave arrival. The distribution of relocated hypocentres clearly shows the configuration of a double-planed shallow seismic zone beneath the Pacific Ocean. The upper plane has a low dip angle near the Japan Trench, increasing gradually to ???30?? at approximately 100 km landward of the Japan Trench. The lower plane is approximately parallel to the upper plane, and appears to be the near-trench counterpart of the lower plane of the double-planed deep seismic zone beneath the land area. The distance between the upper and lower planes is 28-32 km, which is approximately the same as or slightly smaller than that of the double-planed deep seismic zone beneath the land area. Focal mechanism solutions of the relocated earthquakes are determined from P-wave initial motion data. Although P-wave initial motion data for these offshore events are not ideally distributed on the focal sphere, we found that the upper-plane events that occur near the Japan Trench are characterized by normal faulting, whereas lower-plane events are characterized by thrust faulting. This focal mechanism distribution is the opposite to that of the double-planed deep seismic zone beneath the land area. The characteristics of these focal mechanisms for the shallow and deep doubled-planed seismic zones can be explained by a bending-unbending model of the subducting Pacific plate. Some of relocated earthquakes took place in the source area of the 1933 Mw8.4 Sanriku earthquake at depths of 10-23 km. The available focal mechanisms for these events are characterized by normal faulting. Given that the 1933 event was a large normal-fault event that occurred along a fault plane dipping landward, the earthquakes that currently occur just beneath or oceanwards of the Japan Trench are probably its aftershocks, suggesting that aftershock activity continues to the present day, 70 years after the main shock. ?? 2009 The Authors, Journal compilation ?? 2009 RAS.

Geophysical Journal International↗
USGS2009

Shallow seismic structure of Kunlun fault zone in northern Tibetan Plateau, China: Implications for the 2001 M s8.1 Kunlun earthquake

The shallow seismic velocity structure of the Kunlun fault zone (KLFZ) was jointly deduced from seismic refraction profiling and the records of trapped waves that were excited by five explosions. The data were collected after the 2001 Kunlun M s8.1 earthquake in the northern Tibetan Plateau. Seismic phases for the in-line record sections (26 records up to a distance of 15 km) along the fault zone were analysed, and 1-D P - and S -wave velocity models of shallow crust within the fault zone were determined by using the seismic refraction method. Sixteen seismic stations were deployed along the off-line profile perpendicular to the fault zone. Fault-zone trapped waves appear clearly on the record sections, which were simulated with a 3-D finite difference algorithm. Quantitative analysis of the correlation coefficients of the synthetic and observed trapped waveforms indicates that the Kunlun fault-zone width is 300 m, and S -wave quality factor Q within the fault zone is 15. Significantly, S -wave velocities within the fault zone are reduced by 30–45 per cent from surrounding rocks to a depth of at least 1–2 km, while P -wave velocities are reduced by 7–20 per cent. A fault-zone with such P - and S -low velocities is an indication of high fluid pressure because Vs is affected more than Vp . The low-velocity and low- Q zone in the KLFZ model is the effect of multiple ruptures along the fault trace of the 2001 M s8.1 Kunlun earthquake.

Tibetan Plateau↗
USGS2009

A grid-doubling finite-element technique for calculating dynamic three-dimensional spontaneous rupture on an earthquake fault

We present a new finite-element technique for calculating dynamic 3-D spontaneous rupture on an earthquake fault, which can reduce the required computational resources by a factor of six or more, without loss of accuracy. The grid-doubling technique employs small cells in a thin layer surrounding the fault. The remainder of the modelling volume is filled with larger cells, typically two or four times as large as the small cells. In the resulting non-conforming mesh, an interpolation method is used to join the thin layer of smaller cells to the volume of larger cells. Grid-doubling is effective because spontaneous rupture calculations typically require higher spatial resolution on and near the fault than elsewhere in the model volume. The technique can be applied to non-planar faults by morphing, or smoothly distorting, the entire mesh to produce the desired 3-D fault geometry. Using our FaultMod finite-element software, we have tested grid-doubling with both slip-weakening and rate-and-state friction laws, by running the SCEC/USGS 3-D dynamic rupture benchmark problems. We have also applied it to a model of the Hayward fault, Northern California, which uses realistic fault geometry and rock properties. FaultMod implements fault slip using common nodes, which represent motion common to both sides of the fault, and differential nodes, which represent motion of one side of the fault relative to the other side. We describe how to modify the traction-at-split-nodes method to work with common and differential nodes, using an implicit time stepping algorithm.

Geophysical Journal International↗
USGS2009

Rayleigh-wave mode separation by high-resolution linear radon transform

Multichannel analysis of surface waves (MASW) method is an effective tool for obtaining vertical shear wave profiles from a single non-invasive measurement. One key step of the MASW method is generation of a dispersion image and extraction of a reliable dispersion curve from raw multichannel shot records. Because different Rayleigh-wave modes normally interfere with each other in the time and space domain, it is necessary to perform mode separation and reconstruction to increase the accuracy of phase velocities determined from a dispersion image. In this paper, we demonstrate the effectiveness of high-resolution linear Radon transform (LRT) as a means of separating and reconstructing multimode, dispersive Rayleigh-wave energy. We first introduce high-resolution LRT methods and Rayleigh-wave mode separation using high-resolution LRT. Next, we use synthetic data and a real-world example to demonstrate the effectiveness of Rayleigh-wave mode separation using high-resolution LRT. Our synthetic and real-world results demonstrate that (1) high-resolution LRT successfully separates and reconstructs multimode dispersive Rayleigh-wave energy with high resolution allowing the multimode energy to be more accurately determined. The horizontal resolution of the Rayleigh-wave method can be increased by extraction of dispersion curves from a pair of traces in the mode-separated shot gather and (2) multimode separation and reconstruction expand the usable frequency range of higher mode dispersive energy, which increases the depth of investigation and provides a means for accurately determining cut-off frequencies. ?? 2009 The Authors Journal compilation ?? 2009 RAS.

Geophysical Journal International↗