SEARCH · Geology Reports

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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24Linked to original sources

USGS2012

Effects of compression direction on the plasticity and rheology of hydrolytically weakened synthetic quartz crystals at atmospheric pressure

A hydrothermally grown synthetic quartz crystal with 370±60 ppm hydroxyl impurity was cut into right rectangular prisms in eight crystallographic orientations. We compressed the prisms under constant axial force corresponding to a uniaxial stress of 140.0±0.5 MPa, and temperatures of 510° and 750°C. All but one of the samples sustained permanent axial strains of 2–3%. We established the operating slip systems from specimen shape change, slip bands and dislocation etch pits on polished surfaces, crystallographic orientation changes, stress optical features in thin sections, and transmission electron microscopy. The observed creep behavior and plasticity divided the samples into three groups: (1) Crystals compressed at 45° to [0001] and [ ] and those compressed ⊥( ) and ⊥( ) deformed principally by slip parallel to [0001]. Creep rates were relatively high and were not strongly sensitive to test temperature. Dislocation arrays approximately parallel to ( ) are common. Dislocation loops are elongate parallel to [0001], indicating that the edge segments were more mobile than the screw segments. (2) The second group of samples were loaded normal to [0001] in three orientations: ⊥( ), ⊥( ), and at 45° to ( ). These samples deformed primarily by { } 〈a〉 slip with some evidence for secondary slip on the other systems. They were more creep resistant than the first group and displayed a much higher sensitivity of creep rate to test temperature. Dislocation loops are very elongate parallel to [0001], indicating that the screw dislocation segments were much more mobile than the edge segments. (3) A sample compressed parallel to [0001] at 750°C crept at a barely detectable rate (∼10 s) and no optical scale slip features were observed. These results confirm our earlier work on one orientation each from groups 1 and 2, which indicated a strong creep anisotropy for this same crystal. This creep anisotropy parallels a remarkably similar anisotropy in the diffusivity of impurities in quartz, suggesting a causal relationship between impurity diffusion and creep associated with hydrolytic weakening. Appendix Tables A1‐A5 are available with entire article on microfiche. Order from American Geophysical Union, 2000 Florida Avenue, N.W., Washington, DC 20009. Document B84‐004; $2.50. Payment must accompany order.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Continuous monitoring of surface deformation at Long Valley Caldera, California, with GPS

Continuous Global Positioning System (GPS) measurements at Long Valley Caldera, an active volcanic region in east central California, have been made on the south side of the resurgent dome since early 1993. A site on the north side of the dome was added in late 1994. Special adaptations for autonomous operation in remote regions and enhanced vertical precision were made. The data record ongoing volcanic deformation consistent with uplift and expansion of the surface above a shallow magma chamber. Measurement precisions (1 standard error) for “absolute” position coordinates, i.e., relative to a global reference frame, are 3–4 mm (north), 5–6 mm (east), and 10–12 mm (vertical) using 24 hour solutions. Corresponding velocity uncertainties for a 12 month period are about 2 mm/yr in the horizontal components and 3–4 mm/yr in the vertical component. High precision can also be achieved for relative position coordinates on short (less than 10 km) baselines using broadcast ephemerides and observing times as short as 3 hours, even when data are processed rapidly on site. Comparison of baseline length changes across the resurgent dome between the two GPS sites and corresponding two-color electronic distance measurements indicates similar extension rates within error (∼2 mm/yr) once we account for a random walk noise component in both systems that may reflect spurious monument motion. Both data sets suggest a pause in deformation for a 3.5 month period in mid-1995, when the extension rate across the dome decreased essentially to zero. Three dimensional positioning data from the two GPS stations suggest a depth (5.8±1.6 km) and location (west side of the resurgent dome) of a major inflation center, in agreement with other geodetic techniques, near the top of a magma chamber inferred from seismic data. GPS systems similar to those installed at Long Valley can provide a practical method for near real-time monitoring and hazard assessment on many active volcanoes.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Spatiotemporal evolution of dike opening and décollement slip at Kīlauea Volcano, Hawai'i

Rapid changes in ground tilt and GPS positions on Kīlauea Volcano, Hawai'i, are interpreted as resulting from a shallow, two‐segment dike intrusion into the east rift zone that began at 1217 UTC (0217 HST) on 17 June 2007 and lasted almost 3 days. As a result of the intrusion, a very small volume of basalt (about 1500 m 3 ) erupted on 19 June. Northward tilt at a coastal tiltmeter, subsidence of south flank GPS sites, southeastward displacements at southwestern flank GPS sites, and a swarm of flank earthquakes suggest that a slow slip event occurred on the décollement beneath Kīlauea's south flank concurrent with the rift intrusion. We use 4 min GPS positions that include estimates of time‐dependent tropospheric gradients and ground tilt data to study the spatial and temporal relationships between the two inferred shallow, steeply dipping dike segments extending from the surface to about 2 km depth and décollement slip at 8 km depth. We invert for the temporal evolution of distributed dike opening and décollement slip in independent inversions at each time step using a nonnegative least squares algorithm. On the basis of these inversions, the intrusion occurred in two stages that correspond spatially and temporally with concentrated rift zone seismicity. The dike opening began on the western of the two segments before jumping to the eastern segment, where the majority of opening accumulated. Dike opening preceded the start of décollement slip at an 84% confidence level; the latter is indicated by the onset of northward tilt of a coastal tiltmeter. Displacements at southwest flank GPS sites began about 18 h later and are interpreted as resulting from slow slip on the southwestern flank. Additional constraints on the evolution of the intrusion and décollement slip come from inversion of an Envisat interferogram that spans the intrusion until 0822 UTC on 18 June 2007, combined with GPS and tilt data. This inversion shows that up to 0822 UTC on 18 June, décollement slip is only required in a limited region offshore of Ka'ena Point. A similar inversion of the complete event, which includes GPS and tilt data up to 21 June and a second Envisat interferogram spanning the complete intrusion until 21 June, shows décollement slip spread westward across the south flank. This may suggest westward migration of the décollement slip as the event progressed.

Hawaii↗

USGS2019

Locating shallow seismic sources with waves scattered by surface topography: Validation of the method at the Nevada Test Site

Accurate and robust source location is fundamental in seismology. Previously, we proposed a new full waveform location method using waves scattered by topography near the source, and we demonstrated its potential in obtaining accurate source location using synthetic data. In the work presented here, we validate this method with observed seismic data from the 1993 Non-Proliferation Experiment, a chemical explosion with a precisely known location in a region with moderate topography at the Nevada Test Site. We select the first arrivals (Pn/Pg) and their immediate codas to determine the source location and compare our solution with the known location. We use a collocated-grid finite-difference method to calculate the strain Green’s tensor in a grid-search volume containing the source, and obtain three-component synthetic waveforms at 12 broadband seismic stations at regional distances using source-receiver reciprocity. We assess the solution using a least-squares misfit between the observed and synthetic waveforms. When scattered coda waves are incorporated, the best solution is within a few hundreds of meters of the exact source location, and the estimated uncertainty of the solution is reduced compared to the waveform based solution using only the P waves. The solution is robust to the choice of the frequency content and to the addition of random velocity heterogeneity. We conclude that the full waveform source location method is effective for real seismic data, though more validation tests like this one are needed to further understand its efficacy for source-station geometry, roughness of topography, signal quality and other factors.

Nevada↗

USGS2023

High-resolution InSAR reveals localized pre-eruptive deformation inside the crater of Agung Volcano, Indonesia

During a volcanic crisis, high-rate, localized deformation can indicate magma close to the surface, with important implications for eruption forecasting. However, only a few such examples have been reported, because frequent, dense monitoring is needed. High-resolution Synthetic Aperture Radar (SAR) is capable of achieving <1 m spatial resolution and sub-weekly revisit times, but is under-used. Here we use high-resolution satellite SAR imagery from COSMO-SkyMed, TerraSAR-X, and Sentinel-1 to detect intra-crater uplift preceding the November 2017 onset of eruptive activity at Agung, Indonesia. Processing the SAR imagery with an up-to-date, accurate, high-resolution digital elevation model was crucial for preventing aliasing of the deformation signal and for accurate georeferencing. We show that >15 cm of line-of-sight shortening occurred over a 400-by-400 m area on the crater floor in September-October 2017, accompanying a deep seismic swarm and flank dyke intrusion. We attribute the deformation to the pressurization of a shallow (<200 m deep) hydrothermal system by the injection of magmatic gases and fluids. We also observe a second pulse of intra-crater deformation of 3–5 cm within 4 days to 11 hr prior to the first phreatomagmatic eruption, which is consistent with interaction between the hydrothermal system and the ascending magma. This phreatomagmatic eruption created the central pathway used during the final stages of magma ascent. Our observations have important implications for understanding unrest and eruption forecasting, and demonstrate the potential of monitoring with high-resolution SAR.

Agung Volcano, Bali↗

USGS1999

Abrupt along-strike change in tectonic style: San Andreas fault zone, San Francisco Peninsula

Seismicity and high-resolution aeromagnetic data are used to define an abrupt change from compressional to extensional tectonism within a 10- to 15-km-wide zone along the San Andreas fault on the San Francisco Peninsula and offshore from the Golden Gate. This 100-km-long section of the San Andreas fault includes the hypocenter of the M w = 7.8 1906 San Francisco earthquake as well as the highest level of persistent microseismicity along that ∼470-km-long rupture. We define two distinct zones of deformation along this stretch of the fault using well-constrained relocations of all post-1969 earthquakes based a joint one-dimensional velocity/hypocenter inversion and a redetermination of focal mechanisms. The southern zone is characterized by thrust- and reverse-faulting focal mechanisms with NE trending P axes that indicate “fault-normal” compression in 7- to 10-km-wide zones of deformation on both sides of the San Andreas fault. A 1- to 2-km-wide vertical zone beneath the surface trace of the San Andreas is characterized by its almost complete lack of seismicity. The compressional deformation is consistent with the young, high topography of the Santa Cruz Mountains/Coast Ranges as the San Andreas fault makes a broad restraining left bend (∼10°) through the southernmost peninsula. A zone of seismic quiescence ∼15 km long separates this compressional zone to the south from a zone of combined normal-faulting and strike-slip-faulting focal mechanisms (including a M L = 5.3 earthquake in 1957) on the northernmost peninsula and offshore on the Golden Gate platform. Both linear pseudogravity gradients, calculated from the aeromagnetic data, and seismic reflection data indicate that the San Andreas fault makes an abrupt ∼3-km right step less than 5 km offshore in this northern zone. A similar right-stepping (dilatational) geometry is also observed for the subparallel San Gregorio fault offshore. Persistent seismicity and extensional tectonism occur within the San Andreas right stepover region and at least 15 km along-strike both to the SE and NW. The 1906 San Francisco earthquake may have nucleated within the San Andreas right stepover, which may help explain the bilateral nature of rupture of this event. Our analysis suggests two seismic hazards for the San Francisco Peninsula in addition to the hazard associated with a M = 7 to 8 strike-slip earthquake along the San Andreas fault: the potential for a M ≈ 6 normal-faulting earthquake just 5–8 km west of San Francisco and a M = 6+ thrust faulting event in the southern peninsula.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

Groundwater flow, heat transport, and water table position within volcanic edifices: Implications for volcanic processes in the Cascade Range

The position of the water table within a volcanic edifice has significant implications for volcano hazards, geothermal energy, and epithermal mineralization. We have modified the HYDROTHERM numerical simulator to allow for a free-surface (water table) upper boundary condition and a wide range of recharge rates, heat input rates, and thermodynamic conditions representative of continental volcano-hydrothermal systems. An extensive set of simulations was performed on a hypothetical stratovolcano system with unconfined groundwater flow. Simulation results suggest that the permeability structure of the volcanic edifice and underlying material is the dominant control on water table elevation and the distribution of pressures, temperatures, and fluid phases at depth. When permeabilities are isotropic, water table elevation decreases with increasing heat flux and increases with increasing recharge, but when permeabilities are anisotropic, these effects can be much less pronounced. Several conditions facilitate the ascent of a hydrothermal plume into a volcanic edifice: a sufficient source of heat and magmatic volatiles at depth, strong buoyancy forces, and a relatively weak topography-driven flow system. Further, the plume must be connected to a deep heat source through a pathway with a time-averaged effective permeability ???1 ?? 10-16 m2, which may be maintained by frequent seismicity. Topography-driven flow may be retarded by low permeability in the edifice and/or the lack of precipitation recharge; in the latter case, the water table may be relatively deep. Simulation results were compared with observations from the Quaternary stratovolcanoes along the Cascade Range of the western United States to infer hydrothermal processes within the edifices. Extensive ice caps on many Cascade Range stratovolcanoes may restrict recharge on the summits and uppermost flanks. Both the simulation results and limited observational data allow for the possibility that the water table beneath the stratovolcanoes is relatively deep.

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Rapid fluid disruption: A source for self-potential anomalies on volcanoes

Self-potential (SP) anomalies observed above suspected magma reservoirs, dikes, etc., on various volcanoes (Kilauea, Hawaii; Mount Unzen, Japan; Piton de la Fournaise, Reunion Island, Miyake Jima, Japan) result from transient surface electric fields of tens of millivolts per kilometer and generally have a positive polarity. These SP anomalies are usually attributed to electrokinetic effects where properties controlling this process are poorly constrained. We propose an alternate explanation that contributions to electric fields of correct polarity should be expected from charge generation by fluid vaporization/disruption. As liquids are vaporized or removed as droplets by gas transport away from hot dike intrusions, both charge generation and local increase in electrical resistivity by removal of fluids should occur. We report laboratory observations of electric fields in hot rock samples generated by pulses of fluid (water) through the rock at atmospheric pressure. These indicate the relative amplitudes of rapid fluid disruption (RFD) potentials and electrokinetic potentials to be dramatically different and the signals are opposite in sign. Above vaporization temperatures, RFD effects of positive sign in the direction of gas flow dominate, whereas below these temperatures, effects of negative sign dominate. This suggests that the primary contribution to observed self-potential anomalies arises from gas-related charge transport processes at temperatures high enough to produce vigorous boiling and vapor transport. At lower temperatures, the primary contribution is from electrokinetic effects modulated perhaps by changing electrical resistivity and RFD effects from high-pressure but low-temperature CO2 and SO2 gas flow ripping water molecules from saturated crustal rocks. If charge generation is continuous, as could well occur above a newly emplaced dike, positive static potentials will be set up that could be sustained for many years, and the simplest method for identifying these hot, active regions would be to identify the SP anomalies they generate.

Hawai'i↗

USGS2022

Orbital and in-situ investigation of periodic bedrock ridges in Glen Torridon, Gale Crater, Mars

Wind has been the dominant agent of landscape modification on Mars for the past ~3 billion years. Among the diversity of features formed by aeolian abrasion on the surface of Mars are periodic bedrock ridges (PBRs), a relatively recently recognized class of erosional bedforms on Mars for which Earth analogues are rare. Gale crater, the field site for NASA’s Mars Science Laboratory Curiosity rover since it landed there in 2012, contains a diverse and extensive record of aeolian deposition and erosion. This study focuses on a series of periodic, linear bedrock ridges that occur within the Fe/Mg-smectite clay-bearing Glen Torridon region of Aeolis Mons (informally Mount Sharp). During Curiosity’s exploration of the Glen Torridon region between sols ~2300-3080, the rover drove through this field of ridges, providing the first opportunity for the in situ observation of these enigmatic erosional features. This study characterizes the Glen Torridon ridges using orbiter and rover data to determine their morphology, spatial distribution, compositional and material properties, and association with other aeolian features in the area. Based on these observations, the Glen Torridon ridges are interpreted to be consistent with an origin as wind-eroded periodic bedrock ridges carved during the most recent exhumation of Mount Sharp into the present-day mound. Although there is evidence for multidirectional winds in the Glen Torridon region based on the orientation of modern ripples, megaripples, TARs and other bedrock indicators, the consistent orientation of the Glen Torridon ridges, coupled with morphologic asymmetries within the ridges, support formation and elongation of the Glen Torridon PBRs forms parallel to a net regional northerly wind direction in and around Gale crater.

Journal of Geophysical Research E: Planets↗

USGS1997

Correlated errors in geodetic time series: Implications for time-dependent deformation

Analysis of frequent trilateration observations from the two-color electronic distance measuring networks in California demonstrate that the noise power spectra are dominated by white noise at higher frequencies and power law behavior at lower frequencies. In contrast, Earth scientists typically have assumed that only white noise is present in a geodetic time series, since a combination of infrequent measurements and low precision usually preclude identifying the time-correlated signature in such data. After removing a linear trend from the two-color data, it becomes evident that there are primarily two recognizable types of time-correlated noise present in the residuals. The first type is a seasonal variation in displacement which is probably a result of measuring to shallow surface monuments installed in clayey soil which responds to seasonally occurring rainfall; this noise is significant only for a small fraction of the sites analyzed. The second type of correlated noise becomes evident only after spectral analysis of line length changes and shows a functional relation at long periods between power and frequency of 1/ƒ α , where ƒ is frequency and α≈2. With α=2, this type of correlated noise is termed random-walk noise, and its source is mainly thought to be small random motions of geodetic monuments with respect to the Earth's crust, though other sources are possible. Because the line length changes in the two-color networks are measured at irregular intervals, power spectral techniques cannot reliably estimate the level of 1/ƒ α noise. Rather, we also use here a maximum likelihood estimation technique which assumes that there are only two sources of noise in the residual time series (white noise and random-walk noise) and estimates the amount of each. From this analysis we find that the random-walk noise level averages about 1.3 mm/√yr and that our estimates of the white noise component confirm theoretical limitations of the measurement technique. In addition, the seasonal noise can be as large as 3 mm in amplitude but typically is less than 0.5 mm. Because of the presence of random-walk noise in these time series, modeling and interpretation of the geodetic data must account for this source of error. By way of example we show that estimating the time-varying strain tensor (a form of spatial averaging) from geodetic data having both random-walk and white noise error components results in seemingly significant variations in the rate of strain accumulation; spatial averaging does reduce the size of both noise components but not their relative influence on the resulting strain accumulation model.

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Sources of global warming in upper ocean temperature during El Niño

Global average sea surface temperature (SST) from 40°S to 60°N fluctuates ±0.3°C on interannual period scales, with global warming (cooling) during El Niño (La Niña). About 90% of the global warming during El Niño occurs in the tropical global ocean from 20°S to 20°N, half because of large SST anomalies in the tropical Pacific associated with El Niño and the other half because of warm SST anomalies occurring over &sim;80% of the tropical global ocean. From examination of National Centers for Environmental Prediction [ Kalnay et al. , 1996] and Comprehensive Ocean-Atmosphere Data Set [ Woodruff et al. , 1993] reanalyses, tropical global warming during El Niño is associated with higher troposphere moisture content and cloud cover, with reduced trade wind intensity occurring during the onset phase of El Niño. During this onset phase the tropical global average diabatic heat storage tendency in the layer above the main pycnocline is 1–3 W m −2 above normal. Its principal source is a reduction in the poleward Ekman heat flux out of the tropical ocean of 2–5 W m −2 . Subsequently, peak tropical global warming during El Niño is dissipated by an increase in the flux of latent heat to the troposphere of 2–5 W m −2 , with reduced shortwave and longwave radiative fluxes in response to increased cloud cover tending to cancel each other. In the extratropical global ocean the reduction in poleward Ekman heat flux out of the tropics during the onset of El Niño tends to be balanced by reduction in the flux of latent heat to the troposphere. Thus global warming and cooling during Earth's internal mode of interannual climate variability arise from fluctuations in the global hydrological balance, not the global radiation balance. Since it occurs in the absence of extraterrestrial and anthropogenic forcing, global warming on decadal, interdecadal, and centennial period scales may also occur in association with Earth's internal modes of climate variability on those scales.

Journal of Geophysical Research C: Oceans↗

USGS2003

Imaging spectroscopy: Earth and planetary remote sensing with the USGS Tetracorder and expert systems

Imaging spectroscopy is a tool that can be used to spectrally identify and spatially map materials based on their specific chemical bonds. Spectroscopic analysis requires significantly more sophistication than has been employed in conventional broadband remote sensing analysis. We describe a new system that is effective at material identification and mapping: a set of algorithms within an expert system decision‐making framework that we call Tetracorder. The expertise in the system has been derived from scientific knowledge of spectral identification. The expert system rules are implemented in a decision tree where multiple algorithms are applied to spectral analysis, additional expert rules and algorithms can be applied based on initial results, and more decisions are made until spectral analysis is complete. Because certain spectral features are indicative of specific chemical bonds in materials, the system can accurately identify and map those materials. In this paper we describe the framework of the decision making process used for spectral identification, describe specific spectral feature analysis algorithms, and give examples of what analyses and types of maps are possible with imaging spectroscopy data. We also present the expert system rules that describe which diagnostic spectral features are used in the decision making process for a set of spectra of minerals and other common materials. We demonstrate the applications of Tetracorder to identify and map surface minerals, to detect sources of acid rock drainage, and to map vegetation species, ice, melting snow, water, and water pollution, all with one set of expert system rules. Mineral mapping can aid in geologic mapping and fault detection and can provide a better understanding of weathering, mineralization, hydrothermal alteration, and other geologic processes. Environmental site assessment, such as mapping source areas of acid mine drainage, has resulted in the acceleration of site cleanup, saving millions of dollars and years in cleanup time. Imaging spectroscopy data and Tetracorder analysis can be used to study both terrestrial and planetary science problems. Imaging spectroscopy can be used to probe planetary systems, including their atmospheres, oceans, and land surfaces.

Journal of Geophysical Research E: Planets↗

USGS2012

Variations in magma supply rate at Kilauea Volcano, Hawaii

When an eruption of Kilauea lasts more than 4 months, so that a well-defined conduit has time to develop, magma moves freely through the volcano from a deep source to the eruptive site at a constant rate of 0.09 km 3 /yr. At other times, the magma supply rate to Kilauea, estimated from geodetic measurements of surface displacements, may be different. For example, after a large withdrawal of magma from the summit reservoir, such as during a rift zone eruption, the magma supply rate is high initially but then lessens and exponentially decays as the reservoir refills. Different episodes of refilling may have different average rates of magma supply. During four year-long episodes in the 1960s, the annual rate of refilling varied from 0.02 to 0.18 km 3 /yr, bracketing the sustained eruptive rate of 0.09 km 3 /yr. For decade-long or longer periods, our estimate of magma supply rate is based on long-term changes in eruptive rate. We use eruptive rate because after a few dozen eruptions the volume of magma that passes through the summit reservoir is much larger than the net change of volume of magma stored within Kilauea. The low eruptive rate of 0.009 km 3 /yr between 1840 and 1950, compared to an average eruptive rate of 0.05 km 3 /yr since 1950, suggests that the magma supply rate was lower between 1840 and 1950 than it has been since 1950. An obvious difference in activity before and since 1950 was the frequency of rift zone eruptions: eight rift zone eruptions occurred between 1840 and 1950, but more than 20 rift zone eruptions have occurred since 1950. The frequency of rift zone eruptions influences magma supply rate by suddenly lowering pressure of the summit magma reservoir, which feeds magma to rift zone eruptions. A temporary drop of reservoir pressure means a larger-than-normal pressure difference between the reservoir and a deeper source, so magma is forced to move upward into Kilauea at a faster rate.

Journal of Geophysical Research B: Solid Earth↗

USGS2000

The enigma of the Arthur's Pass, New Zealand, earthquake: 1. Reconciling a variety of data for an unusual earthquake sequence

The 1994 Arthur's Pass earthquake ( M W 6.7) is the largest in a recent sequence of earthquakes in the central South Island, New Zealand. No surface rupture was observed, the aftershock distribution was complex, and routine methods of obtaining the faulting orientation of this earthquake proved contradictory. We use a range of data and techniques to obtain our preferred solution, which has a centroid depth of 5 km, M 0 =1.3 × 10 19 N m, and a strike, dip, and rake of 221°, 47°, 112°, respectively. Discrepancies between this solution and the Harvard centroid moment tensor, together with the Global Positioning System (GPS) observations and unusual aftershock distribution, suggest that the rupture may not have occurred on a planar fault. A second, strike slip, subevent on a more northerly striking plane is suggested by these data but neither the body wave modeling nor regional broadband recordings show any complexity or late subevents. We relocate the aftershocks using both one-dimensional and three-dimensional velocity inversions. The depth range of the aftershocks (1–10 km) agrees well with the preferred mainshock centroid depth. The aftershocks near the hypocenter suggest a structure dipping toward the NW, which we interpret to be the mainshock fault plane. This structure and the Harper fault, ∼15 km to the south, appear to have acted as boundaries to the extensive aftershock zone trending NNW-SSE. Most of the M L ≥ 5 aftershocks, including the two largest ( M L 6.1 and M L 5.7), clustered near the Harper fault and have strike slip mechanisms consistent with motion on this fault and its conjugates. Forward modeling of the GPS data suggests that a reverse slip mainshock, combined with strike slip aftershock faulting in the south, is able to match the observed displacements. The occurrence of this earthquake sequence implies that the level of seismic hazard in the central South Island is greater than previous estimates.

Arthur's Pass↗

USGS2006

Distribution of stress drop, stiffness, and fracture energy over earthquake rupture zones

Using information provided by slip models and the methodology of McGarr and Fletcher (2002), we map static stress drop, stiffness (k = ????/u, where ???? is static stress drop and u is slip), and fracture energy over the slip surface to investigate the earthquake rupture process and energy budget. For the 1994 M6.7 Northridge, 1992 M7.3 Landers, and 1995 M6.9 Kobe earthquakes, the distributions of static stress drop show strong heterogeneity, emphasizing the importance of asperities in the rupture process. Average values of static stress drop are 17, 11, and 4 Mpa for Northridge, Landers, and Kobe, respectively. These values are substantially higher than estimates based on simple crack models, suggesting that the failure process involves the rupture of asperities within the larger fault zone. Stress drop as a function of depth for the Northridge and Landers earthquakes suggests that stress drops are limited by crustal strength. For these two earthquakes, regions of high slip are surrounded by high values of stiffness. Particularly for the Northridge earthquake, the prominent patch of high slip in the central part of the fault is bordered by a ring of high stiffness and is consistent with expectations based on the failure of an asperity loaded at its edge due to exterior slip. Stiffness within an asperity is inversely related to its dimensions. Estimates of fracture energy, based on static stress drop, slip, and rupture speed, were used to investigate the nature of slip weakening at four locations near the hypocenter of the Kobe earthquake for comparison with independent results based on a dynamic model of this earthquake. One subfault updip and to the NE of the hypocenter has a fracture energy of 1.1 MJ/m2 and a slip-weakening distance, Dc, of 0.66 m. Right triangles, whose base and height are Dc and the dynamic stress drop, respectively, approximately overlie the slip-dependent stress given by Ide and Takeo (1997) for the same locations near the hypocenter. The total fracture energy for the Kobe earthquake, 3.7 ?? 1014 J, is about the same as the seismic energy (Ea = 3.2 ?? 1014 J.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Basin and range crustal and upper mantle structure, northwest to central Nevada

We present an interpretation of the crustal and uppermost mantle structure of the Basin and Range of northwestern Nevada based on seismic refraction/wide‐angle reflection, near‐vertical reflection, and gravity data. In comparison to most previous estimates, we find that the crust is somewhat thicker (32–36 km versus 22–30 km), and the uppermost mantle velocity is somewhat higher (8.0 km/s versus 7.3–7.9 km/s). Along our transects, the crust is thinnest (32 km) in the Carson Sink‐Buena Vista Valley region and increases by 2–4 km to the west and east, respectively. There is considerable complexity throughout the crust where velocities range from of 2.5 km/s at the surface to 7.4 km/s in the lowermost crust. Variations in velocity and structure of the upper crustal layers reveal apparent basement velocity depressions (areas of lower velocities extending up to 10 km in depth) that underlie some surface ranges as well as the basins. The middle crust rises from about 20 km beneath central Nevada to within 12 km of the surface beneath the area of thinnest crust and is characterized by a modest (∼0.1 km/s) change in velocity and low‐velocity gradients. These midcrustal layers mark the onset of high crustal reflectivity and the apparent limiting depth to which Basin and Range faults can be traced in near‐vertical reflection profiles, suggesting that these midcrustal layers represent the transition between the brittle and ductile zones of the crust. The lower crust is more structurally complex, with layers thickening and thinning in a systematic manner with the upper crustal layers; generally, where there are velocity depressions in the upper crust, the lower crust is thickest and shallowest. The geometry of these lower crustal layers (derived from refraction modeling) coincides with changes in the crustal reflectivity, determined from the Consortium of Continental Reflection Profiling reflection data. The lower crustal layer is unusually high in velocity (7.4 km/s) and is likely the layer identified as mantle in some previous studies. We do not identify the 7.4 km/s layer as mantle because (1) there is an underlying layer with a velocity (8.0 km/s) that is more consistent with the worldwide average velocity for the upper mantle, and (2) the 7.4 km/s layer does not correspond to the “reflection” Moho. Gravity modeling and comparison to existing seismic models show a general consensus in many aspects with respect to crustal structure. This new model forms the basis for speculation on some of the processes associated with rifting of the Basin and Range Province. One such process, lithospheric magmatism, is inferred from the strong attenuation of transmitted seismic waves, which occurs at the same interface at which high‐amplitude, bright spot reflections originate. Unlike previous models, the overall structure and velocity of the crust and uppermost mantle of our new model are similar to other regions worldwide which have undergone high degrees of extension.

Nevada↗

USGS1996

Crustal structure of the Colorado Plateau, Arizona: Application of new long-offset seismic data analysis techniques

The Colorado Plateau is a large crustal block in the southwestern United States that has been raised intact nearly 2 km above sea level since Cretaceous marine sediments were deposited on its surface. Controversy exists concerning the thickness of the plateau crust and the source of its buoyancy. Interpretations of seismic data collected on the plateau vary as to whether the crust is closer to 40 or 50 km thick. A thick crust could support the observed topography of the Colorado Plateau isostatically, while a thinner crust would indicate the presence of an underlying low-density mantle. This paper reports results on long-offset seismic data collected during the 1989 segment of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that extended from the Transition Zone into the Colorado Plateau in northwest Arizona. We apply two new methods to analyze long-offset data that employ finite difference travel time calculations: (1) a first-arrival time inverter to find upper crustal velocity structure and (2) a forward-modeling technique that allows the direct use of the inverted upper crustal solution in modeling secondary reflected arrivals. We find that the crustal thickness increases from 30 km beneath the metamorphic core complexes in the southern Basin and Range province to about 42 km beneath the northern Transition Zone and southern Colorado Plateau margin. We observe some crustal thinning (to ∼37 km thick) and slightly higher lower crustal velocities farther inboard; beneath the Kaibab uplift on the north rim of the Grand Canyon the crust thickens to a maximum of 48 km. We observe a nonuniform crustal thickness beneath the Colorado Plateau that varies by ∼15% and corresponds approximately to variations in topography with the thickest crust underlying the highest elevations. Crustal compositions (as inferred from seismic velocities) appear to be the same beneath the Colorado Plateau as those in the Basin and Range province to the southwest, implying that the plateau crust represents an unextended version of the Basin and Range. Some of the variability in crustal structure appears to correspond to preserved lithospheric discontinuities that date back to the Proterozoic Era.

Journal of Geophysical Research B: Solid Earth↗

USGS2002

Plume-driven plumbing and crustal formation in Iceland

Through combination of surface wave and body wave constraints we derive a three-dimensional (3-D) crustal S velocity model and Moho map for Iceland. It reveals a vast plumbing system feeding mantle plume melt into upper crustal magma chambers where crustal formation takes place. The method is based on the partitioned waveform inversion to which we add additional observations. Love waves from six local events recorded on the HOTSPOT- SIL networks are fitted, S n travel times from the same events measured, previous observations of crustal thickness are added, and all three sets of constraints simultaneously inverted for our 3-D model. In the upper crust (0–15 km) an elongated low-velocity region extends along the length of the Northern, Eastern and Western Neovolcanic Zones. The lowest velocities (−7%) are found at 5–10 km below the two most active volcanic complexes: Hekla and Bárdarbunga-Grímsvötn. In the lower crust (>15 km) the low-velocity region can be represented as a vertical cylinder beneath central Iceland. The low-velocity structure is interpreted as the thermal halo of pipe work which connects the region of melt generation in the uppermost mantle beneath central Iceland to active volcanoes along the neovolcanic zones. Crustal thickness in Iceland varies from 15–20 km beneath the Reykjanes Peninsula, Krafla and the extinct Snæfellsnes rift zone, to 46 km beneath central Iceland. The average crustal thickness is 29 km. The variations in thickness can be explained in terms of the temporal variation in plume productivity over the last ∼20 Myr, the Snæfellsnes rift zone being active during a minimum in plume productivity. Variations in crustal thickness do not depart significantly from an isostatically predicted crustal thickness. The best fit linear isostatic relation implies an average density jump of 4% across the Moho. Rare earth element inversions of basalt compositions on Iceland suggest a melt thickness (i.e., crustal thickness) of 15–20 km, given passive upwelling. The observed crustal thickness of up to 46 km implies active fluxing of source material through the melt zone by the mantle plume at up to 3 times the passive rate.

Journal of Geophysical Research B: Solid Earth↗