USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research, Earth Surface”

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 379 records · Page 21Linked to original sources
USGS2010

Geodetic evidence for en echelon dike emplacement and concurrent slow slip during the June 2007 intrusion and eruption at Kilauea volcano, Hawaii

A series of complex events at Kīlauea Volcano, Hawaii, 17 June to 19 June 2007, began with an intrusion in the upper east rift zone (ERZ) and culminated with a small eruption (1500 m 3 ). Surface deformation due to the intrusion was recorded in unprecedented detail by Global Positioning System (GPS) and tilt networks as well as interferometric synthetic aperture radar (InSAR) data acquired by the ENVISAT and ALOS satellites. A joint nonlinear inversion of GPS, tilt, and InSAR data yields a deflationary source beneath the summit caldera and an ENE-striking uniform-opening dislocation with ~2 m opening, a dip of ∼80° to the south, and extending from the surface to ~2 km depth. This simple model reasonably fits the overall pattern of deformation but significantly misfits data near the western end of an inferred dike-like source. Three more complex dike models are tested that allow for distributed opening including (1) a dike that follows the surface trace of the active rift zone, (2) a dike that follows the symmetry axis of InSAR deformation, and (3) two en echelon dike segments beneath mapped surface cracks and newly formed steaming areas. The en echelon dike model best fits near-field GPS and tilt data. Maximum opening of 2.4 m occurred on the eastern segment beneath the eruptive vent. Although this model represents the best fit to the ERZ data, it still fails to explain data from a coastal tiltmeter and GPS sites on Kīlauea's southwestern flank. The southwest flank GPS sites and the coastal tiltmeter exhibit deformation consistent with observations of previous slow slip events beneath Kīlauea's south flank, but inconsistent with observations of previous intrusions. Slow slip events at Kīlauea and elsewhere are thought to occur in a transition zone between locked and stably sliding zones of a fault. An inversion including slip on a basal decollement improves fit to these data and suggests a maximum of ~15 cm of seaward fault motion, comparable to previous slow-slip events.

Hawai'i↗
USGS2020

Methods for rapidly estimating velocity precision from GNSS time series in the presence of temporal correlation: A new method and comparison of existing methods

Time series of position estimates from Global Navigational Satellite System (GNSS) are used to measure the velocities of points on the surface of the Earth. Along with the velocity estimates, a measure of the precision is needed to assess the quality of the velocity measurement. Here, I evaluate rate uncertainties provided by four different methods that have been applied to geodetic time series. The most rigorous approach uses a data covariance that incorporates a variety of noise processes relevant to geodetic time series but is computationally demanding. Two other approaches are efficient algorithms and are used widely, but both can provide less rigorous estimates of the rate uncertainty. I propose and evaluate a fourth method, which provides estimates of rate uncertainty closer to the rigorous approach but is significantly less computationally demanding. I have evaluated all three methods against the more rigorous method using both simulations and time series from 190 GNSS sites. For data best characterized as having a flicker type noise process, one of the widely used methods overestimates the uncertainty by up to a factor of 2, while the other widely used method underestimates the uncertainty by less than a factor of 2. For a random-walk process, both methods underestimate the rate uncertainty by a factor of 3 to 5.

Journal of Geophysical Research↗
USGS2012

Cooling rate and thermal structure determined from progressive magnetization of the dacite dome at Mount St. Helens, Washington

Our study of a magnetic anomaly associated with the recently active dacite dome at Mount St. Helens suggests that the dome consists of a hot, nonmagnetized core surrounded by a cool, magnetized carapace and flanking talus. The talus does not contribute to the anomaly because its constituent blocks are randomly oriented. Temporal changes in the magnetic anomaly indicate that the magnetized carapace thickened at an average rate of 0.03±0.01 m/d from 1984 to 1986. Petrographic and rock magnetic properties of dome samples indicate that the dominant process responsible for these changes is magnetization of extensively oxidized rock at progressively deeper levels within the dome as the rock cools through its blocking temperature, rather than subsequent changes in magnetization caused by further oxidation. Newly extruded material cools rapidly for a short period as heat is conducted outward in response to convective heat loss from its surface. The cooling rate gradually declines for several weeks, and thereafter the material cools at a relatively constant rate by convective heat loss from its interior along fractures that propagate inward. The rate of internal convective heat loss through fractures varies with rainfall, snowmelt, and large-scale fracturing during subsequent eruptive episodes. In accordance with a model for solidification of the 1959 lava lake at Kilauea Iki, Hawaii, we picture the dome's magnetized carapace as being a two-phase, porous, convective zone separated from the nonmagnetized core of the dome by a thin, single-phase conductive zone. As a consequence of the heat balance between the conductive and convective zones, the blocking-temperature isotherm migrates inward at a relatively constant rate. If the dome remains inactive, the time scale for its complete magnetization is estimated to be 18–36 years, a forecast which can be refined by shallow drilling into the dome and by continuing studies of its growing magnetic anomaly.

Journal of Geophysical Research Solid Earth↗
USGS2012

Comparison of ground motion from tremors and explosions in deep gold mines

Seismic body waves, from tamped chemical explosions, two with yields of 50 and one of 150 kg, were compared with corresponding data from three mining-induced tremors with a view to testing methods of discriminating between the two types of events. Detonated at depths of about 2 km, all three explosions generated P waves for which the low-frequency spectral asymptotes agree well with corresponding results scaled down from nuclear shots at the Nevada Test Site. For the two smaller explosions, recorded underground at hypocentral distances ranging from 234 to 871 m, the P wave corner frequencies, defined by the intersections of the low- and high-frequency asymptotes, are also in fine agreement with scaled-down results from the Nevada Test Site; for the 150-kg explosion, which was recorded only at the surface, it appears that 2.39 km of upward propagation caused the corner frequency to be at least a factor of 4 lower than anticipated due to attenuation. All three explosions generated S waves that appear to be a consequence of deviatoric stress release in the immediate environs of the explosions. The three tremors analyzed here generated P and S waves whose spectra and source parameters agree well with standard earthquake source models. With regard to discrimination based on P wave spectra, we conclude that for events of fixed low-frequency spectral asymptotes, the explosions typically have higher corner frequencies than tremors or earthquakes, although counterexamples certainly exist. Interestingly, the 150-kg explosion was identified as such on the basis of P and S wave polarities that are incompatible with the normally expected double-couple source model; instead, these initial motions are consistent with an explosion in conjunction with normal faulting. The body wave spectra of this explosion and those of a nearby tremor, however, were indistinguishable.

Journal of Geophysical Research Solid Earth↗
USGS2026

Role of volcano-tectonic interactions during early-phase magma-assisted continental rifting

We study volcanic-tectonic interactions during the early phase of rifting in the Natron Rift, a magmatic segment of the East African Rift System. We use the open-source crustal deformation modeling code, PyLith, to model deflation during the 2007–2008 explosive eruptions and inflation during 2016–2023. We compare predictions of surface displacements resulting from magma reservoir inflation from 3D models with and without topography and homogeneous and heterogeneous material properties against Global Navigation Satellite System displacements from 2016 to 2023. We prefer the model with heterogeneous material properties and topography for fitting the observed displacements. Including topography improves the fit to GNSS data and modulates stress changes from shallow magma sources. Our numerical modeling of Coulomb stress changes suggests that deflation of a magma reservoir under the active Ol Doinyo Lengai volcano could trigger oblique normal slip on the Natron Fault with Coulomb stress increases exceeding ∼0.1 MPa on the volcano-adjacent fault segment. Fault slip on the Natron Fault is less likely during inflation of the magma reservoir because normal compressive stress increases on the fault. We also demonstrate that Coulomb stress changes are sensitive to the inclusion of topography in our model, particularly because the depth of the magma reservoir is similar to the amplitude of the topographic relief.

Ol Doinyo Lenga↗
USGS2003

Oceans on Mars: An assessment of the observational evidence and possible fate

If the large Late Hesperian outflow channels were eroded by extensive floods, as appears likely, then large bodies of water must have once occupied the northern plains during that period. Previous estimates of the sizes of bodies of water in the northern lowlands range up to 3 ?? 108 km3. Several contacts have been previously mapped around the edges of the northern plains and interpreted to be shorelines remaining from these former standing bodies of water. We examine the elevations and geologic relations along these contacts in detail and find little support for their interpretation as shorelines. Some contacts are clearly of volcanic origin, and all have significant variations in elevation. Better support for the former presence of water over large parts of the northern plains is provided by the Vastitas Borealis Formation (VBF). Most of the post-Noachian fill within the northern basin is ridged plains of Lower Hesperian age, interpreted to be volcanic in origin. Overlying the ridged plains is the VBF, a thin veneer of material of Upper Hesperian age. The VBF may have been deposited from large floods. Support for this interpretation is the similarity in age between the outflow channels and the VBF, the presence of the VBF at the lower ends of the outflow channels, and identification of numerous features in the outcrop areas of the VBF that are suggestive of basal melting of an ice sheet. To cover all the area over which the VBF is exposed would require ???2.3 ?? 107 km3 of water. Spread over the entire surface of Mars, this volume is equal to a global layer (global equivalent layer, or GEL) ???156 m deep. We find no support for the larger estimates of ocean volumes that range up to 3 ?? 108 km3 and which imply comparable amounts of water per unit area as are currently present on the surface of the Earth. Under present climatic conditions on Mars an ocean would freeze in a geologically short time period (???104 years), then would sublimate away at rates strongly dependent on the presence or absence of debris on the ice surface. The present VBF is interpreted as a sublimation residue from the ponded outflow channel effluents. The fate of a volume of water thought to have been emplaced by the outflow channels (???2.3 ?? 107 km3) is largely accounted for by the presence of other existing reservoirs on the planet. An approximately 20-30 m GEL of water is estimated to be in the present polar caps, and a 50 m GEL may have escaped to space since the Hesperian, leaving ???80 m GEL unaccounted for. This amount may be partly trapped in other volatile-rich deposits on the surface, and a significant amount could have reentered the groundwater system by south polar basal melting and been progressively cold-trapped at the base of a growing cryosphere. On the basis of our assessment of the Hesperian-aged deposits, we predict that testing of the Clifford and Parker [2001] hypothesis that a Noachian-aged ocean covered up to one third of the surface of Mars will be made very difficult by the enhanced degradation rates in the Noachian and subsequent geological events in the northern lowlands.

Journal of Geophysical Research E: Planets↗
USGS2024

Limited preservation of strike-slip surface displacement in the geomorphic record

Offset geomorphic markers are commonly used to interpret slip history of strike-slip faults and have played an important role in forming earthquake recurrence models. These data sets are typically analyzed using cumulative probability methods to interpret average amounts of slip in past earthquakes. However, interpretation of the geomorphic record to infer surface slip history is complicated by slip variability, measurement uncertainty, and modification of offset features in the landscape. To investigate how well geomorphic data record surface slip, we use offset measurements from recent strike-slip surface ruptures ( n = 39), faults with geomorphic evidence of multiple strike-slip earthquakes ( n = 29), and synthetic slip distributions with added noise ( n > 10,000) to examine the constraints of the geomorphic record and the underlying assumptions of the cumulative offset probability distribution analysis method. We find that the geomorphic record is unlikely to resolve more than two paleo-slip distributions, except in specific cases with low slip variability, high slip-per-event, and semiarid climate. In cases where site-specific conditions allow for interpretation of more than two earthquakes, lateral extrapolation along a fault is not straightforward because on-fault displacement and distributed deformation may be spatially variable in each earthquake. We also find that average slip in modern earthquakes is adequately recovered by probability methods, but the reported prevalence of strike-slip faults with characteristic slip history is not supported by geomorphic data. We also propose updated methods to interpret slip history and construct uncertainty bounds for paleo-slip distributions.

Journal of Geophysical Research: Solid Earth↗
USGS2006

Erosion rates at the Mars Exploration Rover landing sites and long-term climate change on Mars

Erosion rates derived from the Gusev cratered plains and the erosion of weak sulfates by saltating sand at Meridiani Planum are so slow that they argue that the present dry and desiccating environment has persisted since the Early Hesperian. In contrast, sedimentary rocks at Meridiani formed in the presence of groundwater and occasional surface water, and many Columbia Hills rocks at Gusev underwent aqueous alteration during the Late Noachian, approximately coeval with a wide variety of geomorphic indicators that indicate a wetter and likely warmer environment. Two-toned rocks, elevated ventifacts, and perched and undercut rocks indicate localized deflation of the Gusev plains and deposition of an equivalent amount of sediment into craters to form hollows, suggesting average erosion rates of ???0.03 nm/yr. Erosion of Hesperian craters, modification of Late Amazonian craters, and the concentration of hematite concretions in the soils of Meridiani yield slightly higher average erosion rates of 1-10 nm/yr in the Amazonian. These erosion rates are 2-5 orders of magnitude lower than the slowest continental denudation rates on Earth, indicating that liquid water was not an active erosional agent. Erosion rates for Meridiani just before deposition of the sulfate-rich sediments and other eroded Noachian areas are comparable with slow denudation rates on Earth that are dominated by liquid water. Available data suggest the climate change at the landing sites from wet and likely warm to dry and desiccating occurred sometime between the Late Noachian and the beginning of the Late Hesperian (3.7-3.5 Ga). Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2022

Three Mw ≥ 4.7 earthquakes within the Changning (China) shale gas field ruptured shallow faults intersecting with hydraulic fracturing wells

From 2017 to 2019, three destructive earthquakes (27 January 2017 Mw 4.7, 16 December 2018 Mw 5.2, and 3 January 2019 Mw 4.8) occurred in the Changning shale gas field in the southwest Sichuan Basin, China. Previous seismological studies attributed these events to hydraulic fracturing (HF), but were unable to identify the causative seismogenic faults and their slip behaviors. Here, we use Sentinel-1 synthetic aperture radar data to measure surface deformation triggered by the three events and conduct geodetic inversions to characterize their rupture models. The resulting coseismic interferograms show prominent surface deformation with the maximum line-of-sight displacements of up to 4 cm. The inversion results show that all three earthquakes mainly ruptured sedimentary formations above the shale gas bed, in the upper 3 km of the crust, with slip magnitudes ranging from 8.5 to 15 cm, and stress drops ranging from ∼1.8 to ∼3.3 MPa. Their source faults intersect with horizontal HF wells, but do not root in the crystalline basement. Combined with the reported difficulty of increasing HF operation pressures prior to the three events, we argue that they were most likely induced by direct injection of pressurized fluids into the fault zones. Crustal deformation patterns inferred from regional topography and GPS velocities highlight that the Changning field is located within a triple junction region near the southeastern margin of the Tibetan Plateau with large deformation gradients; such conditions are not only favorable to the development of critically stressed faults, but also facilitate the occurrence of at least moderate magnitude earthquakes.

Changning shale gas field, Sichuan Basin↗
USGS2016

Origin of the pulse-like signature of shallow long-period volcano seismicity

Short-duration, pulse-like long-period (LP) events are a characteristic type of seismicity accompanying eruptive activity at Mount Etna in Italy in 2004 and 2008 and at Turrialba Volcano in Costa Rica and Ubinas Volcano in Peru in 2009. We use the discrete wave number method to compute the free surface response in the near field of a rectangular tensile crack embedded in a homogeneous elastic half space and to gain insights into the origin of the LP pulses. Two source models are considered, including (1) a vertical fluid-driven crack and (2) a unilateral tensile rupture growing at a fixed sub-Rayleigh velocity with constant opening on a vertical crack. We apply cross correlation to the synthetics and data to demonstrate that a fluid-driven crack provides a natural explanation for these data with realistic source sizes and fluid properties. Our modeling points to shallow sources (<1 km depth), whose signatures are representative of the Rayleigh pulse sampled at epicentral distances >∼1 km. While a slow-rupture failure provides another potential model for these events, the synthetics and resulting fits to the data are not optimal in this model compared to a fluid-driven source. We infer that pulse-like LP signatures are parts of the continuum of responses produced by shallow fluid-driven sources in volcanoes.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

The competition between thermal contraction and differentiation in the stress history of the Moon

The scarcity of both extension and compression features on the Moon strongly constrains the history of the lunar radius—to variations of less than ±1 km over the past 3.8 Gyr. This limit has traditionally been interpreted as requiring a delicate balance between thermal contraction of the near‐surface and expansion of a substantial cold interior region. Recent theories of lunar origin (e.g., giant impact), in contrast, favor a “hot” initial state. We propose that a reconciliation may be possible by taking account of the volume change Δ V / V | d due to differentiation. We calculate STP densities based on simplified normative mineralogies for a suite of estimates of the bulk lunar composition, of primary lunar basalt, and of the residuum left when the maximum amount of the latter is extracted from the former. Typically Δ V / V | d ≃ 2 to 5%—an expansion equivalent to heating by ∼10 3 K. Provided the timing of differentiation is correct, one might offset the cooling of a magma ocean as much as 630 km deep by differentiation of the remainder of the Moon (which need not start much below the solidus temperature). A large but not impossible amount of gabbroic melt production is implied: ∼100 times the volume of mare basalts known to have been extruded. We do not address the detailed genetic relationship of this melt to the basalts observed on the lunar surface but point out that it need not have reached the surface directly or even have entered the crust in order for the expansion to have occurred. To assess the timing of melt formation, we investigate a simple conductive lunar thermal model which takes account of both Δ V / V | d and thermal contraction. Our initial state is characterized by a central temperature T c and a depth Z 0 above which the material (derived from the magma ocean) is already at the solidus and is not suceptible to volume changes upon further differentiation. We find a range of models satisfying the limits on radius increase and decrease. The hottest has T c = 1210 K, Z 0 = 400 km; without Δ V / V | d , we would need a larger or colder (or both) core, e.g., T c ≲ 700 K for Z 0 = 200–400 km, in agreement with previous investigators. Our modeling thus lends credence to the idea that the Moon could have been initially ≳50% molten (with the remainder relatively close to the solidus) and yet experienced little volume change over the last 3.8 Gyr.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

Synthetic aperture radar interferometry of Okmok volcano, Alaska: radar observations

ERS-1/ERS-2 synthetic aperture radar interferometry was used to study the 1997 eruption of Okmok volcano in Alaska. First, we derived an accurate digital elevation model (DEM) using a tandem ERS-1/ERS-2 image pair and the preexisting DEM. Second, by studying changes in interferometric coherence we found that the newly erupted lava lost radar coherence for 5-17 months after the eruption. This suggests changes in the surface backscattering characteristics and was probably related to cooling and compaction processes. Third, the atmospheric delay anomalies in the deformation interferograms were quantitatively assessed. Atmospheric delay anomalies in some of the interferograms were significant and consistently smaller than one to two fringes in magnitude. For this reason, repeat observations are important to confidently interpret small geophysical signals related to volcanic activities. Finally, using two-pass differential interferometry, we analyzed the preemptive inflation, coeruptive deflation, and posteruptive inflation and confirmed the observations using independent image pairs. We observed more than 140 cm of subsidence associated with the 1997 eruption. This subsidence occurred between 16 months before the eruption and 5 months after the eruption, was preceded by ∼18 cm of uplift between 1992 and 1995 centered in the same location, and was followed by ∼10 cm of uplift between September 1997 and 1998. The best fitting model suggests the magma reservoir resided at 2.7 km depth beneath the center of the caldera, which was ∼5 km from the eruptive vent. We estimated the volume of the erupted material to be 0.055 km 3 and the average thickness of the erupted lava to be ∼7.4 m. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

California↗
USGS2004

Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults

We argue that key features of thrust earthquake triggering, inhibition, and clustering can be explained by Coulomb stress changes, which we illustrate by a suite of representative models and by detailed examples. Whereas slip on surface-cutting thrust faults drops the stress in most of the adjacent crust, slip on blind thrust faults increases the stress on some nearby zones, particularly above the source fault. Blind thrusts can thus trigger slip on secondary faults at shallow depth and typically produce broadly distributed aftershocks. Short thrust ruptures are particularly efficient at triggering earthquakes of similar size on adjacent thrust faults. We calculate that during a progressive thrust sequence in central California the 1983 M w = 6.7 Coalinga earthquake brought the subsequent 1983 M w = 6.0 Nuñez and 1985 M w = 6.0 Kettleman Hills ruptures 10 bars and 1 bar closer to Coulomb failure. The idealized stress change calculations also reconcile the distribution of seismicity accompanying large subduction events, in agreement with findings of prior investigations. Subduction zone ruptures are calculated to promote normal faulting events in the outer rise and to promote thrust-faulting events on the periphery of the seismic rupture and its downdip extension. These features are evident in aftershocks of the 1957 M w = 9.1 Aleutian and other large subduction earthquakes. We further examine stress changes on the rupture surface imparted by the 1960 M w = 9.5 and 1995 M w = 8.1 Chile earthquakes, for which detailed slip models are available. Calculated Coulomb stress increases of 2–20 bars correspond closely to sites of aftershocks and postseismic slip, whereas aftershocks are absent where the stress drops by more than 10 bars. We also argue that slip on major strike-slip systems modulates the stress acting on nearby thrust and strike-slip faults. We calculate that the 1857 M w = 7.9 Fort Tejon earthquake on the San Andreas fault and subsequent interseismic slip brought the Coalinga fault ∼1 bar closer to failure but inhibited failure elsewhere on the Coast Ranges thrust faults. The 1857 earthquake also promoted failure on the White Wolf reverse fault by 8 bars, which ruptured in the 1952 M w = 7.3 Kern County shock but inhibited slip on the left-lateral Garlock fault, which has not ruptured since 1857. We thus contend that stress transfer exerts a control on the seismicity of thrust faults across a broad spectrum of spatial and temporal scales.

California↗
USGS2002

Imaging the mantle beneath Iceland using integrated seismological techniques

Using a combination of body wave and surface wave data sets to reveal the mantle plume and plume head, this study presents a tomographic image of the mantle structure beneath Iceland to 400 km depth. Data comes primarily from the PASSCAL-HOTSPOT deployment of 30 broadband instruments over a period of 2 years, and is supplemented by data from the SIL and ICEMELT networks. Three sets of relative teleseismic body wave arrival times are generated through cross correlation: S and SKS arrivals at 0.03–0.1 Hz, and P and PKIKP arrivals at 0.03–0.1 and 0.8–2.0 Hz. Prior to inversion the crustal portion of the travel time anomalies is removed using the crustal model ICECRTb. This step has a significant effect on the mantle velocity variations imaged down to a depth of ∼250 km. Inversion of relative arrival times only provides information on lateral velocity variations. Surface waves are therefore used to provide absolute velocity information for the uppermost mantle beneath Iceland. The average wave number for the Love wave fundamental mode at 0.020 and 0.024 Hz is measured and used to invert for the average S velocity. Combination of the body wave and surface wave information reveals a predominantly horizontal low-velocity anomaly extending from the Moho down to ∼250 km depth, interpreted as a plume head. Below the plume head a near-cylindrical low-velocity anomaly with a radius of ∼100 km and peak V P and V S anomalies of −2% and −4%, respectively, extends down to the maximum depth of resolution at 400 km. Within the plume head, in the uppermost mantle above the core of the plume, there is a relatively high velocity with a maximum V P and V S anomaly of +2%. This high-velocity anomaly may be the result of the extreme degree of melt extraction necessary to generate the thick (46 km) crust in central Iceland. Comparison of the plume volumetric flux implied by our images, the crustal generation rate, and the degree of melting suggested by rare earth element inversions, suggests that (1) mantle material must be flowing horizontally away from the plume core faster than the overlying lithosphere and (2) the bulk of the plume material does not participate in melting beneath Iceland.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Three-dimensional velocity structure of Siletzia and other accreted terranes in the Cascadia forearc of Washington

Eocene mafic crust with high seismic velocities underlies much of the Oregon and Washington forearc and acts as a backstop for accretion of marine sedimentary rocks from the obliquely subducting Juan de Fuca slab. Arc-parallel migration of relatively strong blocks of this terrane, known as Siletzia, focuses upper crustal deformation along block boundaries, which are potential sources of earthquakes. In a three-dimensional velocity model of coastal Washington, we have combined surface geology, well data, and travel times from earthquakes and controlled source seismic experiments to resolve the major boundaries of the Siletz terrane with the adjacent accreted sedimentary prism and volcanic arc. In southern Washington and northern Oregon the Siletz terrane appears to be a thick block (???20 km) that extends west of the coastline and makes a high-angle contact with the offshore accreted sedimentary prism. On its east flank the high-velocity Siletz terrane boundary coincides with an en echelon zone of seismicity in the arc. In northern Washington the western edge of Siletzia makes a lower-angled, fault-bound contact with the accretionary prism. In addition, alternating, east-west trending uplifts and downwarps of the Siletz terrane centered on the antiformal Olympic Mountains may reflect focusing of north-south compression in the northern part of the Siletz terrane. This compressional strain may result from northward transport and clockwise rotation of the Siletz terrane into the relatively fixed Canadian Coast Mountains restraining bend along the coast.

Oregon, Washington↗
USGS2023

Quantifying site effects and their influence on earthquake source parameter estimations using a dense array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

Oklahoma↗
USGS2012

Limitations of quantitative analysis of deep crustal seismic reflection data: Examples from GLIMPCE

Amplitude preservation in seismic reflection data can be obtained by a relative true amplitude (RTA) processing technique in which the relative strength of reflection amplitudes is preserved vertically as well as horizontally, after compensating for amplitude distortion by near-surface effects and propagation effects. Quantitative analysis of relative true amplitudes of the Great Lakes International Multidisciplinary Program on Crustal Evolution seismic data is hampered by large uncertainties in estimates of the water bottom reflection coefficient and the vertical amplitude correction and by inadequate noise suppression. Processing techniques such as deconvolution, F-K filtering, and migration significantly change the overall shape of amplitude curves and hence calculation of reflection coefficients and average reflectance. Thus lithological interpretation of deep crustal seismic data based on the absolute value of estimated reflection strength alone is meaningless. The relative strength of individual events, however, is preserved on curves generated at different stages in the processing. We suggest that qualitative comparisons of relative strength, if used carefully, provide a meaningful measure of variations in reflectivity. Simple theoretical models indicate that peg-leg multiples rather than water bottom multiples are the most severe source of noise contamination. These multiples are extremely difficult to remove when the water bottom reflection coefficient is large (>0.6), a condition that exists beneath parts of Lake Superior and most of Lake Huron.

Great Lakes↗