USGS Science⌕ Search

SEARCH · USGS Science

Results for “Pure and Applied Geophysics”

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 19 recordsLinked to original sources
USGS2017

Performance of Irikura recipe rupture model generator in earthquake ground motion simulations with Graves and Pitarka hybrid approach

We analyzed the performance of the Irikura and Miyake (Pure and Applied Geophysics 168(2011):85–104, 2011) (IM2011) asperity-based kinematic rupture model generator, as implemented in the hybrid broadband ground motion simulation methodology of Graves and Pitarka (Bulletin of the Seismological Society of America 100(5A):2095–2123, 2010), for simulating ground motion from crustal earthquakes of intermediate size. The primary objective of our study is to investigate the transportability of IM2011 into the framework used by the Southern California Earthquake Center broadband simulation platform. In our analysis, we performed broadband (0–20 Hz) ground motion simulations for a suite of M6.7 crustal scenario earthquakes in a hard rock seismic velocity structure using rupture models produced with both IM2011 and the rupture generation method of Graves and Pitarka (Bulletin of the Seismological Society of America, 2016) (GP2016). The level of simulated ground motions for the two approaches compare favorably with median estimates obtained from the 2014 Next Generation Attenuation-West2 Project (NGA-West2) ground motion prediction equations (GMPEs) over the frequency band 0.1–10 Hz and for distances out to 22 km from the fault. We also found that, compared to GP2016, IM2011 generates ground motion with larger variability, particularly at near-fault distances (<12 km) and at long periods (>1 s). For this specific scenario, the largest systematic difference in ground motion level for the two approaches occurs in the period band 1–3 s where the IM2011 motions are about 20–30% lower than those for GP2016. We found that increasing the rupture speed by 20% on the asperities in IM2011 produced ground motions in the 1–3 s bandwidth that are in much closer agreement with the GMPE medians and similar to those obtained with GP2016. The potential implications of this modification for other rupture mechanisms and magnitudes are not yet fully understood, and this topic is the subject of ongoing study. We concluded that IM2011 rupture generator performs well in ground motion simulations using Graves and Pitarka hybrid method. Therefore, we recommend it to be considered for inclusion into the framework used by the Southern California Earthquake Center broadband simulation platform.

Pure and Applied Geophysics↗
USGS2019

High-frequency ground motion and source characteristics of the 2008 Wenchuan and 2013 Lushan, China, earthquakes

The 2008 M W 7.9 Wenchuan and the 2013 M W 6.6 Lushan earthquakes, which both occurred on the Longmen Shan thrust belt, show some interesting similarities and differences. Whereas the Wenchuan earthquake entailed a rupture zone that extended about 300 km northeastward, with fault slip extending to the surface, the Lushan earthquake was the result of a buried and much more compact zone of rupture. The high-frequency ground motions, however, for these two earthquakes, as measured by the peak ground acceleration, were evidently influenced by neither the extent of rupture nor the presence or absence of surface rupture. The source parameters for these two earthquakes tend to confirm the idea that high-frequency ground motion is controlled by stress changes in the rupture zone that give rise to the radiated ground acceleration. The apparent stresses for the Wenchuan and Lushan earthquakes are about 0.5 and 0.75 MPa, respectively, and the stress drops, in the same order, are about 2.5 and 3.5 MPa. The ratios of average stress drop to apparent stress are in the range 4.5–5 for both events, consistent with expectations based on the Brune (J Geophys Res 75(26):4997–5009, 1970) source model.

Lushan, Wenchuan↗
USGS2005

The inverse problem of refraction travel times, part I: Types of Geophysical Nonuniqueness through Minimization

In a set of two papers we study the inverse problem of refraction travel times. The purpose of this work is to use the study as a basis for development of more sophisticated methods for finding more reliable solutions to the inverse problem of refraction travel times, which is known to be nonunique. The first paper, "Types of Geophysical Nonuniqueness through Minimization," emphasizes the existence of different forms of nonuniqueness in the realm of inverse geophysical problems. Each type of nonuniqueness requires a different type and amount of a priori information to acquire a reliable solution. Based on such coupling, a nonuniqueness classification is designed. Therefore, since most inverse geophysical problems are nonunique, each inverse problem must be studied to define what type of nonuniqueness it belongs to and thus determine what type of a priori information is necessary to find a realistic solution. The second paper, "Quantifying Refraction Nonuniqueness Using a Three-layer Model," serves as an example of such an approach. However, its main purpose is to provide a better understanding of the inverse refraction problem by studying the type of nonuniqueness it possesses. An approach for obtaining a realistic solution to the inverse refraction problem is planned to be offered in a third paper that is in preparation. The main goal of this paper is to redefine the existing generalized notion of nonuniqueness and a priori information by offering a classified, discriminate structure. Nonuniqueness is often encountered when trying to solve inverse problems. However, possible nonuniqueness diversity is typically neglected and nonuniqueness is regarded as a whole, as an unpleasant "black box" and is approached in the same manner by applying smoothing constraints, damping constraints with respect to the solution increment and, rarely, damping constraints with respect to some sparse reference information about the true parameters. In practice, when solving geophysical problems different types of nonuniqueness exist, and thus there are different ways to solve the problems. Nonuniqueness is usually regarded as due to data error, assuming the true geology is acceptably approximated by simple mathematical models. Compounding the nonlinear problems, geophysical applications routinely exhibit exact-data nonuniqueness even for models with very few parameters adding to the nonuniqueness due to data error. While nonuniqueness variations have been defined earlier, they have not been linked to specific use of a priori information necessary to resolve each case. Four types of nonuniqueness, typical for minimization problems are defined with the corresponding methods for inclusion of a priori information to find a realistic solution without resorting to a non-discriminative approach. The above-developed stand-alone classification is expected to be helpful when solving any geophysical inverse problems. ?? Birkha??user Verlag, Basel, 2005.

Pure and Applied Geophysics↗
USGS1985

Prediction of central California earthquakes from soil-gas helium fluctuations

The observations of short-term decreases in helium soil-gas concentrations along the San Andreas Fault in central California have been correlated with subsequent earthquake activity. The area of study is elliptical in shape with radii approximately 160??80 km, centered near San Benito, and with the major axis parallel to the Fault. For 83 percent of the M>4 earthquakes in this area a helium decrease preceded seismic activity by 1.5 to 6.5 weeks. There were several earthquakes without a decrease and several decreases without a corresponding earthquake. Owing to complex and unresolved interaction of many geophysical and geochemical parameters, no suitable model is yet developed to explain the observations. ?? 1985 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗
USGS2011

Book review: Rogue waves in the ocean

Rogue Waves in the Ocean (2009) is a follow-on text to Extreme Ocean Waves (2008) edited by Pelinovsky and Kharif, both published by Springer. Unlike the earlier text, which is a compilation of papers on a variety of extreme waves that was the subject of a scientific conference in 2007, Rogues Waves in the Ocean is written, rather than edited, by Kharif, Pelinovsky, and Slunyaev and is focused on rogue waves in particular. The book consists of six chapters covering 216 pages. As the subject matter of each chapter is distinct, references appear at the end of each chapter rather than at the end of the book. The preface shows how each of the chapters relates to the larger study of rogue waves. The result is a book with a nice mix of eyewitness observations, physical theory, and statistics. Review info: Rogue Waves in the Ocean. Advances in Geophysical and Environmental Mechanics and Mathematics. B y Christian Kharif, Efim Pelinovsky and Alexey Slunyaev, 2009. ISBN: 978-3540884187, xiii, 216 pp.

Pure and Applied Geophysics↗
USGS2008

Interactive visualization to advance earthquake simulation

The geological sciences are challenged to manage and interpret increasing volumes of data as observations and simulations increase in size and complexity. For example, simulations of earthquake-related processes typically generate complex, time-varying data sets in two or more dimensions. To facilitate interpretation and analysis of these data sets, evaluate the underlying models, and to drive future calculations, we have developed methods of interactive visualization with a special focus on using immersive virtual reality (VR) environments to interact with models of Earth's surface and interior. Virtual mapping tools allow virtual "field studies" in inaccessible regions. Interactive tools allow us to manipulate shapes in order to construct models of geological features for geodynamic models, while feature extraction tools support quantitative measurement of structures that emerge from numerical simulation or field observations, thereby enabling us to improve our interpretation of the dynamical processes that drive earthquakes. VR has traditionally been used primarily as a presentation tool, albeit with active navigation through data. Reaping the full intellectual benefits of immersive VR as a tool for scientific analysis requires building on the method's strengths, that is, using both 3D perception and interaction with observed or simulated data. This approach also takes advantage of the specialized skills of geological scientists who are trained to interpret, the often limited, geological and geophysical data available from field observations. ?? Birkhaueser 2008.

Pure and Applied Geophysics↗
USGS1976

Relative influence upon microwave emissivity of fine-scale stratigraphy, internal scattering, and dielectric properties

The microwave emissivity of relatively low-loss media such as snow, ice, frozen ground, and lunar soil is strongly influenced by fine-scale layering and by internal scattering. Radiometric data, however, are commonly interpreted using a model of emission from a homogeneous, dielectric halfspace whose emissivity derives exclusively from dielectric properties. Conclusions based upon these simple interpretations can be erroneous. Examples are presented showing that the emission from fresh or hardpacked snow over either frozen or moist soil is governed dominantly by the size distribution of ice grains in the snowpack. Similarly, the thickness of seasonally frozen soil and the concentration of rock clasts in lunar soil noticeably affect, respectively, the emissivities of northern latitude soils in winter and of the lunar regolith. Petrophysical data accumulated in support of the geophysical interpretation of microwave data must include measurements of not only dielectric properties, but also of geometric factors such as finescale layering and size distributions of grains, inclusions, and voids.

Pure and Applied Geophysics PAGEOPH↗
USGS1985

A study of the depth of weathering and its relationship to the mechanical properties of near-surface rocks in the Mojave Desert

Weathered granite extends 70 m deep at Hi Vista in the arid central Mojave Desert of southern California. The low strength of this granite is due to the alteration of biotite and chlorite montmorillonite. Deep weathering probably occurs in most granites, although we cannot rule out some anomalous mechanisms at Hi Vista. Geophysical instruments set in these slightly altered rocks are limited by the unstable behavior of the rocks. Thus, tectonic signals from instruments placed in shallow boreholes give vague results. Geophysical measurements of these weathered rocks resemble measurements of granitic rocks near major faults. The rheology of the rocks in which instruments are placed limits the useful sensitivity of the instruments. ?? 1985 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗
USGS2008

Data-resolution matrix and model-resolution matrix for Rayleigh-wave inversion using a damped least-squares method

Inversion of multimode surface-wave data is of increasing interest in the near-surface geophysics community. For a given near-surface geophysical problem, it is essential to understand how well the data, calculated according to a layered-earth model, might match the observed data. A data-resolution matrix is a function of the data kernel (determined by a geophysical model and a priori information applied to the problem), not the data. A data-resolution matrix of high-frequency (>2 Hz) Rayleigh-wave phase velocities, therefore, offers a quantitative tool for designing field surveys and predicting the match between calculated and observed data. We employed a data-resolution matrix to select data that would be well predicted and we find that there are advantages of incorporating higher modes in inversion. The resulting discussion using the data-resolution matrix provides insight into the process of inverting Rayleigh-wave phase velocities with higher-mode data to estimate S-wave velocity structure. Discussion also suggested that each near-surface geophysical target can only be resolved using Rayleigh-wave phase velocities within specific frequency ranges, and higher-mode data are normally more accurately predicted than fundamental-mode data because of restrictions on the data kernel for the inversion system. We used synthetic and real-world examples to demonstrate that selected data with the data-resolution matrix can provide better inversion results and to explain with the data-resolution matrix why incorporating higher-mode data in inversion can provide better results. We also calculated model-resolution matrices in these examples to show the potential of increasing model resolution with selected surface-wave data. ?? Birkhaueser 2008.

Pure and Applied Geophysics↗
USGS2016

A possible source mechanism of the 1946 Unimak Alaska far-field tsunami, uplift of the mid-slope terrace above a splay fault zone

In 1946, megathrust seismicity along the Unimak segment of the Alaska subduction zone generated the largest ever recorded Alaska/Aleutian tsunami. The tsunami severely damaged Pacific islands and coastal areas from Alaska to Antarctica. It is the charter member of “tsunami” earthquakes that produce outsized far-field tsunamis for the recorded magnitude. Its source mechanisms were unconstrained by observations because geophysical data for the Unimak segment were sparse and of low resolution. Reprocessing of legacy geophysical data reveals a deep water, high-angle reverse or splay thrust fault zone that leads megathrust slip upward to the mid-slope terrace seafloor rather than along the plate boundary toward the trench axis. Splay fault uplift elevates the outer mid-slope terrace and its inner area subsides. Multibeam bathymetry along the splay fault zone shows recent but undated seafloor disruption. The structural configuration of the nearby Semidi segment is similar to that of the Unimak segment, portending generation of a future large tsunami directed toward the US West coast.

Alaska↗
USGS2010

A trade-off solution between model resolution and covariance in surface-wave inversion

Regularization is necessary for inversion of ill-posed geophysical problems. Appraisal of inverse models is essential for meaningful interpretation of these models. Because uncertainties are associated with regularization parameters, extra conditions are usually required to determine proper parameters for assessing inverse models. Commonly used techniques for assessment of a geophysical inverse model derived (generally iteratively) from a linear system are based on calculating the model resolution and the model covariance matrices. Because the model resolution and the model covariance matrices of the regularized solutions are controlled by the regularization parameter, direct assessment of inverse models using only the covariance matrix may provide incorrect results. To assess an inverted model, we use the concept of a trade-off between model resolution and covariance to find a proper regularization parameter with singular values calculated in the last iteration. We plot the singular values from large to small to form a singular value plot. A proper regularization parameter is normally the first singular value that approaches zero in the plot. With this regularization parameter, we obtain a trade-off solution between model resolution and model covariance in the vicinity of a regularized solution. The unit covariance matrix can then be used to calculate error bars of the inverse model at a resolution level determined by the regularization parameter. We demonstrate this approach with both synthetic and real surface-wave data. ?? 2010 Birkh??user / Springer Basel AG.

Pure and Applied Geophysics↗
USGS2005

The inverse problem of refraction travel times, part II: Quantifying refraction nonuniqueness using a three-layer model

This paper is the second of a set of two papers in which we study the inverse refraction problem. The first paper, "Types of Geophysical Nonuniqueness through Minimization," studies and classifies the types of nonuniqueness that exist when solving inverse problems depending on the participation of a priori information required to obtain reliable solutions of inverse geophysical problems. In view of the classification developed, in this paper we study the type of nonuniqueness associated with the inverse refraction problem. An approach for obtaining a realistic solution to the inverse refraction problem is offered in a third paper that is in preparation. The nonuniqueness of the inverse refraction problem is examined by using a simple three-layer model. Like many other inverse geophysical problems, the inverse refraction problem does not have a unique solution. Conventionally, nonuniqueness is considered to be a result of insufficient data and/or error in the data, for any fixed number of model parameters. This study illustrates that even for overdetermined and error free data, nonlinear inverse refraction problems exhibit exact-data nonuniqueness, which further complicates the problem of nonuniqueness. By evaluating the nonuniqueness of the inverse refraction problem, this paper targets the improvement of refraction inversion algorithms, and as a result, the achievement of more realistic solutions. The nonuniqueness of the inverse refraction problem is examined initially by using a simple three-layer model. The observations and conclusions of the three-layer model nonuniqueness study are used to evaluate the nonuniqueness of more complicated n-layer models and multi-parameter cell models such as in refraction tomography. For any fixed number of model parameters, the inverse refraction problem exhibits continuous ranges of exact-data nonuniqueness. Such an unfavorable type of nonuniqueness can be uniquely solved only by providing abundant a priori information. Insufficient a priori information during the inversion is the reason why refraction methods often may not produce desired results or even fail. This work also demonstrates that the application of the smoothing constraints, typical when solving ill-posed inverse problems, has a dual and contradictory role when applied to the ill-posed inverse problem of refraction travel times. This observation indicates that smoothing constraints may play such a two-fold role when applied to other inverse problems. Other factors that contribute to inverse-refraction-problem nonuniqueness are also considered, including indeterminacy, statistical data-error distribution, numerical error and instability, finite data, and model parameters. ?? Birkha??user Verlag, Basel, 2005.

Pure and Applied Geophysics↗
USGS1977

Earthquake stress drops, ambient tectonic stresses and stresses that drive plate motions

A variety of geophysical observations suggests that the upper portion of the lithosphere, herein referred to as the elastic plate, has long-term material properties and frictional strength significantly greater than the lower lithosphere. If the average frictional stress along the non-ridge margin of the elastic plate is of the order of a kilobar, as suggested by the many observations of the frictional strength of rocks at mid-crustal conditions of pressure and temperature, the only viable mechanism for driving the motion of the elastic plate is a basal shear stress of several tens of bars. Kilobars of tectonic stress are then an ambient, steady condition of the earth's crust and uppermost mantle. The approximate equality of the basal shear stress and the average crustal earthquake stress drop, the localization of strain release for major plate margin earthquakes, and the rough equivalence of plate margin slip rates and gross plate motion rates suggest that the stress drops of major plate margin earthquakes are controlled by the elastic release of the basal shear stress in the vicinity of the plate margin, despite the existence of kilobars of tectonic stress existing across vertical planes parallel to the plate margin. If the stress differences available to be released at the time of faulting are distributed in a random, white fasbion with a mean-square value determined by the average earthquake stress drop, the frequency of occurrence of constant stress drop earthquakes will be proportional to reciprocal faulting area, in accordance with empirically known frequency of occurrence statistics. ?? 1977 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗
USGS2009

Dipping-interface mapping using mode-separated Rayleigh waves

Multichannel analysis of surface waves (MASW) method is a non-invasive geophysical technique that uses the dispersive characteristic of Rayleigh waves to estimate a vertical shear (S)-wave velocity profile. A pseudo-2D S-wave velocity section is constructed by aligning 1D S-wave velocity profiles at the midpoint of each receiver spread that are contoured using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. Based on the assumption that a dipping-layer model can be regarded as stepped flat layers, high-resolution linear Radon transform (LRT) has been proposed to image Rayleigh-wave dispersive energy and separate modes of Rayleigh waves from a multichannel record. With the mode-separation technique, therefore, a dispersion curve that possesses satisfactory accuracy can be calculated using a pair of consecutive traces within a mode-separated shot gather. In this study, using synthetic models containing a dipping layer with a slope of 5, 10, 15, 20, or 30 degrees and a real-world example, we assess the ability of using high-resolution LRT to image and separate fundamental-mode Rayleigh waves from raw surface-wave data and accuracy of dispersion curves generated by a pair of consecutive traces within a mode-separated shot gather. Results of synthetic and real-world examples demonstrate that a dipping interface with a slope smaller than 15 degrees can be successfully mapped by separated fundamental waves using high-resolution LRT. ?? Birkh??user Verlag, Basel 2009.

Pure and Applied Geophysics↗
USGS1994

Ductile creep and compaction: A mechanism for transiently increasing fluid pressure in mostly sealed fault zones

A simple cyclic process is proposed to explain why major strike-slip fault zones, including the San Andreas, are weak. Field and laboratory studies suggest that the fluid within fault zones is often mostly sealed from that in the surrounding country rock. Ductile creep driven by the difference between fluid pressure and lithostatic pressure within a fault zone leads to compaction that increases fluid pressure. The increased fluid pressure allows frictional failure in earthquakes at shear tractions far below those required when fluid pressure is hydrostatic. The frictional slip associated with earthquakes creates porosity in the fault zone. The cycle adjusts so that no net porosity is created (if the fault zone remains constant width). The fluid pressure within the fault zone reaches long-term dynamic equilibrium with the (hydrostatic) pressure in the country rock. One-dimensional models of this process lead to repeatable and predictable earthquake cycles. However, even modest complexity, such as two parallel fault splays with different pressure histories, will lead to complicated earthquake cycles. Two-dimensional calculations allowed computation of stress and fluid pressure as a function of depth but had complicated behavior with the unacceptable feature that numerical nodes failed one at a time rather than in large earthquakes. A possible way to remove this unphysical feature from the models would be to include a failure law in which the coefficient of friction increases at first with frictional slip, stabilizing the fault, and then decreases with further slip, destabilizing it. ?? 1994 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗
USGS1994

Direct observation of frictional contacts: New insights for state-dependent properties

Rocks and many other materials display a rather complicated, but characteristic, dependence of friction on sliding history. These effects are well-described by empirical rate- and state-dependent constitutive formulations which have been utilized for analysis of fault slip and earthquake processes. We present a procedure for direct quantitative microscopic observation of frictional contacts during slip. The observations reveal that frictional state dependence represents an increase of contact area with contact age. Transient changes of sliding resistance correlate with changes in contact area and arise from shifts of contact population age. Displacement-dependent replacement of contact populations is shown to cause the diagnostic evolution of friction over a characteristic sliding distance that occurs whenever slip begins or sliding conditions change. ?? 1994 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗
USGS1994

Earthquake mechanism and predictability shown by a laboratory fault

Slip events generated in a laboratory fault model consisting of a circulinear chain of eight spring-connected blocks of approximately equal weight elastically driven to slide on a frictional surface are studied. It is found that most of the input strain energy is released by a relatively few large events, which are approximately time predictable. A large event tends to roughen stress distribution along the fault, whereas the subsequent smaller events tend to smooth the stress distribution and prepare a condition of simultaneous criticality for the occurrence of the next large event. The frequency-size distribution resembles the Gutenberg-Richter relation for earthquakes, except for a falloff for the largest events due to the finite energy-storage capacity of the fault system. Slip distributions, in different events are commonly dissimilar. Stress drop, slip velocity, and rupture velocity all tend to increase with event size. Rupture-initiation locations are usually not close to the maximum-slip locations. ?? 1994 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗