Search USGS⌕ Search

SEARCH · Search USGS

Results for “Geophysical Journal”

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
USGS2011

Geophysical Research Letters: New policies improve top-cited geosciences journal

Geophysical Research Letters (GRL) is the American Geophysical Union's premier journal of fast, groundbreaking communication. It rapidly publishes high- impact,letter-length articles, and it is the top-cited multidisciplinary geosciences journal over the past 10 years, with an impact factor that increased again in 2009, to 3.204. For manuscripts submitted to GRL, the median time to first and final decision is 23 and 27 days, respectively—a 35% improvement since 2007—and the median time from submission to publication is 13 weeks for 90% of GRL papers—a 25% improvement since 2007. Among high-impact publications in the geosciences, GRL has the fastest turnaround.

Eos, Transactions, American Geophysical Union↗
USGS2011

Lithosphere-asthenosphere interaction beneath the western United States from the joint inversion of body-wave traveltimes and surface-wave phase velocities

The relation between the complex geological history of the western margin of the North American plate and the processes in the mantle is still not fully documented and understood. Several pre-USArray local seismic studies showed how the characteristics of key geological features such as the Colorado Plateau and the Yellowstone Snake River Plains are linked to their deep mantle structure. Recent body-wave models based on the deployment of the high density, large aperture USArray have provided far more details on the mantle structure while surface-wave tomography (ballistic waves and noise correlations) informs us on the shallow structure. Here we combine constraints from these two data sets to image and study the link between the geology of the western United States, the shallow structure of the Earth and the convective processes in mantle. Our multiphase DNA10-S model provides new constraints on the extent of the Archean lithosphere imaged as a large, deeply rooted fast body that encompasses the stable Great Plains and a large portion of the Northern and Central Rocky Mountains. Widespread slow anomalies are found in the lower crust and upper mantle, suggesting that low-density rocks isostatically sustain part of the high topography of the western United States. The Yellowstone anomaly is imaged as a large slow body rising from the lower mantle, intruding the overlying lithosphere and controlling locally the seismicity and the topography. The large E-W extent of the USArray used in this study allows imaging the 'slab graveyard', a sequence of Farallon fragments aligned with the currently subducting Juan de Fuca Slab, north of the Mendocino Triple Junction. The lithospheric root of the Colorado Plateau has apparently been weakened and partly removed through dripping. The distribution of the slower regions around the Colorado Plateau and other rigid blocks follows closely the trend of Cenozoic volcanic fields and ancient lithospheric sutures, suggesting that the later exert a control on the locus of magmato-tectonic activity today. The DNA velocity models are available for download and slicing at http://dna.berkeley.edu. ?? 2011 The Authors Geophysical Journal International ?? 2011 RAS.

Geophysical Journal International↗
USGS2010

Locating non-volcanic tremor along the San Andreas Fault using a multiple array source imaging technique

Non-volcanic tremor (NVT) has been observed at several subduction zones and at the San Andreas Fault (SAF). Tremor locations are commonly derived by cross-correlating envelope-transformed seismic traces in combination with source-scanning techniques. Recently, they have also been located by using relative relocations with master events, that is low-frequency earthquakes that are part of the tremor; locations are derived by conventional traveltime-based methods. Here we present a method to locate the sources of NVT using an imaging approach for multiple array data. The performance of the method is checked with synthetic tests and the relocation of earthquakes. We also applied the method to tremor occurring near Cholame, California. A set of small-aperture arrays (i.e. an array consisting of arrays) installed around Cholame provided the data set for this study. We observed several tremor episodes and located tremor sources in the vicinity of SAF. During individual tremor episodes, we observed a systematic change of source location, indicating rapid migration of the tremor source along SAF. ?? 2010 The Authors Geophysical Journal International ?? 2010 RAS.

Geophysical Journal International↗
USGS2011

Critical steps for the continuing advancement of hydrogeophysics

Special hydrogeophysics issues published by hydrology and geophysics journals, special sessions and workshops at conferences, and an increasing number of short courses demonstrate the growing interest in the use of geophysics for hydrologic investigations. The formation of the hydrogeophysics technical subcommittee of AGU's Hydrology section adds further evidence of the recognized significance of this growing interdisciplinary field. Given the clear value of nondestructive and nonintrusive imaging for subsurface investigations, we believe the advances in the adoption of existing geophysical methods, the development of novel methods, and the merging of geophysical and other data made in hydrogeophysics could be applied to a wide range of geological, environmental, and engineering applications.

Eos Science News↗
USGS2011

Coulomb stress change sensitivity due to variability in mainshock source models and receiving fault parameters: A case study of the 2010-2011 Christchurch, New Zealand, earthquakes

Strong aftershocks following major earthquakes present significant challenges for infrastructure recovery as well as for emergency rescue efforts. A tragic instance of this is the 22 February 2011 M w 6.3 Christchurch aftershock in New Zealand, which caused more than 100 deaths while the 2010 M w 7.1 Canterbury mainshock did not cause a single fatality (Figure 1). Therefore, substantial efforts have been directed toward understanding the generation mechanisms of aftershocks as well as mitigating hazards due to aftershocks. Among these efforts are the prediction of strong aftershocks, earthquake early warning, and aftershock probability assessment. Zhang et al. (1999) reported a successful case of strong aftershock prediction with precursory data such as changes in seismicity pattern, variation of b-value, and geomagnetic anomalies. However, official reports of such successful predictions in geophysical journals are extremely rare, implying that deterministic prediction of potentially damaging aftershocks is not necessarily more scientifically feasible than prediction of mainshocks.

Christchurch↗
USGS1990

Sediment movement along the U.S. east coast continental shelf-I. Estimates of bottom stress using the Grant-Madsen model and near-bottom wave and current measurements

Bottom stress is calculated for several long-term time-series observations, made on the U.S. east coast continental shelf during winter, using the wave-current interaction and moveable bed models of Grant and Madsen (1979, Journal of Geophysical Research, 84, 1797-1808; 1982, Journal of Geophysical Research, 87, 469-482). The wave and current measurements were obtained by means of a bottom tripod system which measured current using a Savonius rotor and vane and waves by means of a pressure sensor. The variables were burst sampled about 10% of the time. Wave energy was reasonably resolved, although aliased by wave groupiness, and wave period was accurate to 1-2 s during large storms. Errors in current speed and direction depend on the speed of the mean current relative to the wave current. In general, errors in bottom stress caused by uncertainties in measured current speed and wave characteristics were 10-20%. During storms, the bottom stress calculated using the Grant-Madsen models exceeded stress computed from conventional drag laws by a factor of about 1.5 on average and 3 or more during storm peaks. Thus, even in water as deep as 80 m, oscillatory near-bottom currents associated with surface gravity waves of period 12 s or longer will contribute substantially to bottom stress. Given that the Grant-Madsen model is correct, parameterizations of bottom stress that do not incorporate wave effects will substantially underestimate stress and sediment transport in this region of the continental shelf.

Georges Bank, Mid-Atlantic Bight↗
USGS2011

A trans-dimensional Bayesian Markov chain Monte Carlo algorithm for model assessment using frequency-domain electromagnetic data

A meaningful interpretation of geophysical measurements requires an assessment of the space of models that are consistent with the data, rather than just a single, ‘best’ model which does not convey information about parameter uncertainty. For this purpose, a trans-dimensional Bayesian Markov chain Monte Carlo (MCMC) algorithm is developed for assessing frequencydomain electromagnetic (FDEM) data acquired from airborne or ground-based systems. By sampling the distribution of models that are consistent with measured data and any prior knowledge, valuable inferences can be made about parameter values such as the likely depth to an interface, the distribution of possible resistivity values as a function of depth and non-unique relationships between parameters. The trans-dimensional aspect of the algorithm allows the number of layers to be a free parameter that is controlled by the data, where models with fewer layers are inherently favoured, which provides a natural measure of parsimony and a significant degree of flexibility in parametrization. The MCMC algorithm is used with synthetic examples to illustrate how the distribution of acceptable models is affected by the choice of prior information, the system geometry and configuration and the uncertainty in the measured system elevation. An airborne FDEM data set that was acquired for the purpose of hydrogeological characterization is also studied. The results compare favorably with traditional least-squares analysis, borehole resistivity and lithology logs from the site, and also provide new information about parameter uncertainty necessary for model assessment.

Geophysical Journal International↗
USGS2011

High-frequency Born synthetic seismograms based on coupled normal modes

High-frequency and full waveform synthetic seismograms on a 3-D laterally heterogeneous earth model are simulated using the theory of coupled normal modes. The set of coupled integral equations that describe the 3-D response are simplified into a set of uncoupled integral equations by using the Born approximation to calculate scattered wavefields and the pure-path approximation to modulate the phase of incident and scattered wavefields. This depends upon a decomposition of the aspherical structure into smooth and rough components. The uncoupled integral equations are discretized and solved in the frequency domain, and time domain results are obtained by inverse Fourier transform. Examples show the utility of the normal mode approach to synthesize the seismic wavefields resulting from interaction with a combination of rough and smooth structural heterogeneities. This approach is applied to an ∼4 Hz shallow crustal wave propagation around the site of the San Andreas Fault Observatory at Depth (SAFOD).

Geophysical Journal International↗
USGS2006

Constraints on the mechanism of long-term, steady subsidence at Medicine Lake volcano, northern California, from GPS, leveling, and InSAR

Leveling surveys across Medicine Lake volcano (MLV) have documented subsidence that is centered on the summit caldera and decays symmetrically on the flanks of the edifice. Possible mechanisms for this deformation include fluid withdrawal from a subsurface reservoir, cooling/crystallization of subsurface magma, loading by the volcano and dense intrusions, and crustal thinning due to tectonic extension (Dzurisin et al., 1991 [Dzurisin, D., Donnelly-Nolan, J.M., Evans, J.R., Walter, S.R., 1991. Crustal subsidence, seismicity, and structure near Medicine Lake Volcano, California. Journal of Geophysical Research 96, 16, 319-16, 333.]; Dzurisin et al., 2002 [Dzurisin, D., Poland, M.P., Bürgmann, R., 2002. Steady subsidence of Medicine Lake Volcano, Northern California, revealed by repeated leveling surveys. Journal of Geophysical Research 107, 2372, doi:10.1029/2001JB000893.]). InSAR data that approximate vertical displacements are similar to the leveling results; however, vertical deformation data alone are not sufficient to distinguish between source mechanisms. Horizontal displacements from GPS were collected in the Mt. Shasta/MLV region in 1996, 1999, 2000, 2003, and 2004. These results suggest that the region is part of the western Oregon block that is rotating about an Euler pole in eastern Oregon. With this rotation removed, most sites in the network have negligible velocities except for those near MLV caldera. There, measured horizontal velocities are less than predicted from ∼10 km deep point and dislocation sources of volume loss based on the leveling data; therefore volumetric losses simulated by these sources are probably not causing the observed subsidence at MLV. This result demonstrates that elastic models of subsurface volume change can provide misleading results where additional geophysical and geological constraints are unavailable, or if only vertical deformation is known. The deformation source must be capable of causing broad vertical deformation with comparatively smaller horizontal displacements. Thermoelastic contraction of a column of hot rock beneath the volcano cannot reproduce the observed ratio of vertical to horizontal surface displacements. Models that determine deformation due to loading by the volcano and dense intrusions can be made to fit the pattern of vertical displacements by assuming a weak upper crust beneath MLV, though the subsidence rates due to surface loading must be lower than the observed displacements. Tectonic extension is almost certainly occurring based on fault orientations and focal mechanisms, but does not appear to be a major contributor to the observed deformation. We favor a model that includes a combination of sources, including extension and loading of a hot weak crust with thermal contraction of a cooling mass of rock beneath MLV, which are processes that are probably occurring at MLV. Future microgravity surveys and the planned deployment of an array of continuous GPS stations as part of a Plate Boundary Observatory volcano cluster will help to refine this model.

California↗
USGS2006

Body-wave traveltime and amplitude shifts from asymptotic travelling wave coupling

We explore the sensitivity of finite-frequency body-wave traveltimes and amplitudes to perturbations in 3-D seismic velocity structure relative to a spherically symmetric model. Using the approach of coupled travelling wave theory, we consider the effect of a structural perturbation on an isolated portion of the seismogram. By convolving the spectrum of the differential seismogram with the spectrum of a narrow window taper, and using a Taylor's series expansion for wavenumber as a function of frequency on a mode dispersion branch, we derive semi-analytic expressions for the sensitivity kernels. Far-field effects of wave interactions with the free surface or internal discontinuities are implicitly included, as are wave conversions upon scattering. The kernels may be computed rapidly for the purpose of structural inversions. We give examples of traveltime sensitivity kernels for regional wave propagation at 1 Hz. For the direct SV wave in a simple crustal velocity model, they are generally complicated because of interfering waves generated by interactions with the free surface and the Mohorovic??ic?? discontinuity. A large part of the interference effects may be eliminated by restricting the travelling wave basis set to those waves within a certain range of horizontal phase velocity. ?? Journal compilation ?? 2006 RAS.

Geophysical Journal International↗
USGS1997

Introduction to special section: The Trans-Alaska Crustal Transect (TACT) across Arctic Alaska

This special section of the Journal of Geophysical Research addresses the composition and structural evolution of the lithosphere in northern Alaska. Investigations reported in this section were mainly undertaken as part of the Trans‐Alaska Crustal Transect (TACT), an integrated geological and geophysical transect of the entire Alaskan lithosphere along a north‐south corridor undertaken from 1984 to 1992 (Figure 1). The onshore segment of the transect approximately follows along the route of the trans‐Alaskan pipeline; the offshore segment extends across the continental margin in the Gulf of Alaska to the Pacific plate. The TACT line is unique in that it provides a coordinated onshore/offshore geological and geophysical traverse of the North American plate in Alaska from the active convergent Pacific margin to the passive Arctic margin of the continent.

Alaska↗
USGS2010

Time-dependent seismic tomography

Of methods for measuring temporal changes in seismic-wave speeds in the Earth, seismic tomography is among those that offer the highest spatial resolution. 3-D tomographic methods are commonly applied in this context by inverting seismic wave arrival time data sets from different epochs independently and assuming that differences in the derived structures represent real temporal variations. This assumption is dangerous because the results of independent inversions would differ even if the structure in the Earth did not change, due to observational errors and differences in the seismic ray distributions. The latter effect may be especially severe when data sets include earthquake swarms or aftershock sequences, and may produce the appearance of correlation between structural changes and seismicity when the wave speeds are actually temporally invariant. A better approach, which makes it possible to assess what changes are truly required by the data, is to invert multiple data sets simultaneously, minimizing the difference between models for different epochs as well as the rms arrival-time residuals. This problem leads, in the case of two epochs, to a system of normal equations whose order is twice as great as for a single epoch. The direct solution of this system would require twice as much memory and four times as much computational effort as would independent inversions. We present an algorithm, tomo4d, that takes advantage of the structure and sparseness of the system to obtain the solution with essentially no more effort than independent inversions require. No claim to original US government works Journal compilation ?? 2010 RAS.

Geophysical Journal International↗
USGS2006

Regional intensity attenuation models for France and the estimation of magnitude and location of historical earthquakes

Intensity assignments for 33 calibration earthquakes were used to develop intensity attenuation models for the Alps, Armorican, Provence, Pyrenees and Rhine regions of France. Intensity decreases with ?? most rapidly in the French Alps, Provence and Pyrenees regions, and least rapidly in the Armorican and Rhine regions. The comparable Armorican and Rhine region attenuation models are aggregated into a French stable continental region model and the comparable Provence and Pyrenees region models are aggregated into a Southern France model. We analyse MSK intensity assignments using the technique of Bakun & Wentworth, which provides an objective method for estimating epicentral location and intensity magnitude MI. MI for the 1356 October 18 earthquake in the French stable continental region is 6.6 for a location near Basle, Switzerland, and moment magnitude M is 5.9-7.2 at the 95 per cent (??2??) confidence level. MI for the 1909 June 11 Trevaresse (Lambesc) earthquake near Marseilles in the Southern France region is 5.5, and M is 4.9-6.0 at the 95 per cent confidence level. Bootstrap resampling techniques are used to calculate objective, reproducible 67 per cent and 95 per cent confidence regions for the locations of historical earthquakes. These confidence regions for location provide an attractive alternative to the macroseismic epicentre and qualitative location uncertainties used heretofore. ?? 2006 The Authors Journal compilation ?? 2006 RAS.

Geophysical Journal International↗
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↗
USGS2012

Introduction to special section on the California-Arizona crustal transect: CACTIS, Part 3

The CACTIS (California‐Arizona Crustal Transect Interim Synthesis) workshop in May 1988 brought researchers together in Flagstaff, Arizona, to discuss the geologic evolution and crustal structure of the southern Cordillera between the San Andreas fault in southeastern California and the Colorado Plateau in Arizona [ Sass et al ., 1988]. The first set of papers resulting from the workshop appeared in the Journal of Geophysical Research (JGR) as the special CACTIS 1 section in January 1990 [ Howard et al ., 1990], and a second as CACTIS 2 in November 1990 [ Haxel et al ., 1990]. Twelve papers in this issue form part 3 of the CACTIS series. A final collection of papers will appear in the future as a joint special section with the California Consortium for Crustal Studies (CALCRUST).

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Post-seismic relaxation following the great 2004 Sumatra-Andaman earthquake on a compressible self-gravitating Earth

he Mw ??? 9.0 2004 December 26 Sumatra-Andaman and Mw =8.7 2005 March 28 Nias earthquakes, which collectively ruptured approximately 1800 km of the Andaman and Sunda subduction zones, are expected to be followed by vigorous viscoelastic relaxation involving both the upper and lower mantle. Because of these large spatial dimensions it is desirable to fully account for gravitational coupling effects in the relaxation process. We present a stable method of computing relaxation of a spherically-stratified, compressible and self-gravitating viscoelastic Earth following an impulsive moment release event. The solution is cast in terms of a spherical harmonic expansion of viscoelastic normal modes. For simple layered viscoelastic models, which include a low-viscosity oceanic asthenosphere, we predict substantial post-seismic effects over a region several 100s of km wide surrounding the eastern Indian Ocean. We compare observed GPS time-series from ten regional sites (mostly in Thailand and Indonesia), beginning in 2004 December, with synthetic time-series that include the coseismic and post-seismic effects of the 2004 December 26 and 2005 March 28 earthquakes. A viscosity structure involving a biviscous (Burgers body) rheology in the asthenosphere explains the pattern and amplitude of post-seismic offsets remarkably well. ?? 2006 The Authors Journal compilation ?? 2006 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↗
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↗