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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22Linked to original sources

USGS2012

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth↗

USGS2012

Three-dimensional P velocity structure of the crust below Newberry Volcano, Oregon

Utilizing teleseismic P residuals, we have detected a column of high P velocity material extending from within 10 km of the surface below the summit of Newberry Volcano, Oregon, to midcrustal depths near 25 km. We interpret this column to be the expression of a swarm of predominantly subsolidus gabbroic sills and dikes which were intruded as the volcano was built. The high P velocities observed below the volcano severely limit the size of magma chambers which could presently exist in the crust below Newberry Volcano. Those possible include a few percent of partial melt distributed through large volumes of a mafic intrusion zone in the midcrust; a few, smaller, higher melt fraction zones in the midcrust with dimensions less than 6 km and whose aggregate volume is only a few percent of enclosing volumes of 200 km 3 ; small magma bodies with dimensions of a few kilometers located within the upper 10 km of the crust; or a mafic, crystal-rich magma of arbitrary dimensions located in the upper few km. The low P velocities detected in the upper 4 km beneath the center of the summit caldera may be partially caused by a magma chamber in the second of these catagories.

Journal of Geophysical Research Solid Earth↗

USGS2009

Thermal conductivity of hydrate-bearing sediments

A thorough understanding of the thermal conductivity of hydrate-bearing sediments is necessary for evaluating phase transformation processes that would accompany energy production from gas hydrate deposits and for estimating regional heat flow based on the observed depth to the base of the gas hydrate stability zone. The coexistence of multiple phases (gas hydrate, liquid and gas pore fill, and solid sediment grains) and their complex spatial arrangement hinder the a priori prediction of the thermal conductivity of hydrate-bearing sediments. Previous studies have been unable to capture the full parameter space covered by variations in grain size, specific surface, degree of saturation, nature of pore filling material, and effective stress for hydrate-bearing samples. Here we report on systematic measurements of the thermal conductivity of air dry, water- and tetrohydrofuran (THF)-saturated, and THF hydrate–saturated sand and clay samples at vertical effective stress of 0.05 to 1 MPa (corresponding to depths as great as 100 m below seafloor). Results reveal that the bulk thermal conductivity of the samples in every case reflects a complex interplay among particle size, effective stress, porosity, and fluid-versus-hydrate filled pore spaces. The thermal conductivity of THF hydrate–bearing soils increases upon hydrate formation although the thermal conductivities of THF solution and THF hydrate are almost the same. Several mechanisms can contribute to this effect including cryogenic suction during hydrate crystal growth and the ensuing porosity reduction in the surrounding sediment, increased mean effective stress due to hydrate formation under zero lateral strain conditions, and decreased interface thermal impedance as grain-liquid interfaces are transformed into grain-hydrate interfaces.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

The Augustine magmatic system as revealed by seismic tomography and relocated earthquake hypocenters from 1994 through 2009

We incorporate 14 years of earthquake data from the Alaska Volcano Observatory with data from a 1975 controlled‐source seismic experiment to obtain the three‐dimensional P and S wave velocity structure and the first high‐precision earthquake locations at Augustine Volcano to be calculated in a fully three‐dimensional velocity model. Velocity tomography shows two main features beneath Augustine: a narrow, high‐velocity column beneath the summit, extending from ∼2 km depth to the surface, and elevated velocities on the south flank. Our relocation results allow a thorough analysis of the spatio‐temoral patterns of seismicity and the relationship to the magmatic and eruptive activity. Background seismicity is centered beneath the summit at an average depth of 0.6 km above sea level. In the weeks leading to the January 2006 eruption of Augustine, seismicity focused on a NW‐SE line along the trend of an inflating dike. A series of drumbeat earthquakes occurred in the early weeks of the eruption, indicating further magma transport through the same dike system. During the six months following the onset of the eruption, the otherwise quiescent region 1 to 5 km below sea level centered beneath the summit became seismically active with two groups of earthquakes, differentiated by frequency content. The deep longer‐period earthquakes occurred during the eruption and are interpreted as resulting from the movement of magma toward the summit, and the post‐eruptive shorter‐period earthquakes may be due to the relaxation of an emptied magma tube. The seismicity subsequently returned to its normal background rates and patterns.

Alaska↗

USGS1999

Overview of the Mars Pathfinder Mission: Launch through landing, surface operations, data sets, and science results

Mars Pathfinder successfully landed at Ares Vallis on July 4, 1997, deployed and navigated a small rover about 100 m clockwise around the lander, and collected data from three science instruments and ten technology experiments. The mission operated for three months and returned 2.3 Gbits of data, including over 16,500 lander and 550 rover images, 16 chemical analyses of rocks and soil, and 8.5 million individual temperature, pressure and wind measurements. Path‐finder is the best known location on Mars, having been clearly identified with respect to other features on the surface by correlating five prominent horizon features and two small craters in lander images with those in high‐resolution orbiter images and in inertial space from two‐way ranging and Doppler tracking. Tracking of the lander has fixed the spin pole of Mars, determined the precession rate since Viking 20 years ago, and indicates a polar moment of inertia, which constrains a central metallic core to be between 1300 and ∼2000 km in radius. Dark rocks appear to be high in silica and geochemically similar to anorogenic andesites; lighter rocks are richer in sulfur and lower in silica, consistent with being coated with various amounts of dust. Rover and lander images show rocks with a variety of morphologies, fabrics and textures, suggesting a variety of rock types are present. Rounded pebbles and cobbles on the surface as well as rounded bumps and pits on some rocks indicate these rocks may be conglomerates (although other explanations are also possible), which almost definitely require liquid water to form and a warmer and wetter past. Air‐borne dust is composed of composite silicate particles with a small fraction of a highly magnetic mineral, interpreted to be most likely maghemite; explanations suggest iron was dissolved from crustal materials during an active hydrologic cycle with maghemite freeze dried onto silicate dust grains. Remote sensing data at a scale of a kilometer or greater and an Earth analog correctly predicted a rocky plain safe for landing and roving with a variety of rocks deposited by catstrophic floods, which are relatively dust free. The surface appears to have changed little since it formed billions of years ago, with the exception that eolian activity may have deflated the surface by ∼3–7 cm, sculpted wind tails, collected sand into dunes, and eroded ventifacts (fluted and grooved rocks). Pathfinder found a dusty lower atmosphere, early morning water ice clouds, and morning near‐surface air temperatures that changed abruptly with time and height. Small scale vortices, interpreted to be dust devils, were observed repeatedly in the afternoon by the meteorology instruments and have been imaged.

Journal of Geophysical Research E: Planets↗

USGS2009

Three‐dimensional model of Hellenic Arc deformation and origin of the Cretan uplift

[1] The Hellenic Arc of Greece is the most seismically active part of Europe, but little is know about its mechanics. We modeled deformation along the arc using a finite element model. The model was intended to capture large‐scale 3‐D structure of Nubian plate subduction beneath the Aegean block and its deformational consequences. The shape of the interface was developed using mapped traces at the surface and earthquake hypocenters at depth. Model block motions were constrained by recent compilations of GPS velocity vectors. We simulated a 10 ka period of convergence between Nubia and the Aegean and calculated the strain field in the overriding plate as well as the spatial distribution and orientation of differential stress (∣ σ 1 − σ 3 ∣). From these calculations we derived testable quantities such as the expected seismic moment rate on the interplate contact, uplift pattern, and distribution of strain modes. Our relatively simple model broadly reproduced observed uplift patterns, earthquake activity, and loci of extension and contraction. The model showed a localization of uplift near the island of Crete, where the fastest Aegean uplift rates are well documented. Comparison of calculated expected seismic moment and observed earthquake catalogs implies a nearly fully coupled interplate contact. On the basis of our modeling results, we suggest that south Aegean deformation is driven primarily by the fast moving (∼33 mm a −1 ) Aegean upper plate overriding a nearly stalled (∼5 mm a −1 ) Nubian lower plate. This tectonic setting thus more closely resembles a continental thrust than it does a typical oceanic subduction zone.

Aegean Sea, Hellenic Arc↗

USGS2012

The sea slope problem revisited

Discrepancies in sea surface topography based on comparisons between the results of steric leveling and repeated geodetic levelings have identified what is known as the ‘sea slope problem.’ This problem is actually twofold: (1) the sea surface relief based on steric leveling differs significantly from that based on geodetic leveling along several generally north-south coasts, and (2) successively propagated levelings between several widely separated tide stations indicate that the stationary sea slope seemingly has been changing with time, whereas differenced sea level means between these stations indicate that the sea surface relief has remained virtually invariant during the same intervals. Reexamination of the three reported discrepancies between geodetically and sterically determined sea slopes indicates that the Australian example is based on leveling of a quality inappropriate to the comparison. The discrepancy developed along the Atlantic coast of the United States is limited to the reach between Portsmouth–Hampton Roads, Virginia, and Key West, Florida, where the accuracy of steric leveling may be especially vulnerable owing to the dynamic effects of the Gulf Stream. Reconsideration of the example along the Pacific coast of the United States indicates that the various discrepancies are due chiefly to intrasurvey movement and resultant distortion of geodetically defined height differences between tide stations. Agreement between the results of steric and geodetic leveling along tectonically inactive north-south coasts devoid of strong boundary currents is generally good. This observation supports the conclusions (1) that any directionally dependent systematic error in geodetic leveling is measurably insignificant and (2) that where allowance is made for the possible effects of major boundary currents or intrasurvey movement during levelings between tide stations, the sea slope problem tends to vanish.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Sedimentary deposits in the northern lowland plains, Mars

The lowland plains on Mars have surfaces marked by large polygonal fracture patterns. It was recently proposed that the fracture patterns were developed on sedimentary deposits from outflow channels. We support this hypothesis because of the following observations. (1) Polygonal fracture patterns tend to occur in low areas on Mars that apparently received influx of sediments; the area of northern plains, in which the fractures occur, also forms a regional low. (2) The fracture patterns are concentrated near reentrants of the plains projecting into the southern highlands. Not only are these reentrants associated with outflow channels, but channels and polygonally fractured deposits have similar ages. (3) Sinuous ridges at the mouths of Martian outflow channels resemble ridges in Antarctic ice streams and ice shelves; the similarity suggests that the Martian material was transported from the southern highlands toward the northern plains, as the Antarctic material is transported from the continent toward the adjacent ocean. Overall, it appears that a massive transfer of material took place midway in Martian history, either through the outflow channels or elsewhere along the northern highland scarp, and that this material probably formed the deposits now characterized by polygonal fracture patterns.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Estimates of stress drop and crustal tectonic stress from the 27 February 2010 Maule, Chile, earthquake: Implications for fault strength

The great 27 February 2010 M w 8.8 earthquake off the coast of southern Chile ruptured a ∼600 km length of subduction zone. In this paper, we make two independent estimates of shear stress in the crust in the region of the Chile earthquake. First, we use a coseismic slip model constrained by geodetic observations from interferometric synthetic aperture radar (InSAR) and GPS to derive a spatially variable estimate of the change in static shear stress along the ruptured fault. Second, we use a static force balance model to constrain the crustal shear stress required to simultaneously support observed fore‐arc topography and the stress orientation indicated by the earthquake focal mechanism. This includes the derivation of a semianalytic solution for the stress field exerted by surface and Moho topography loading the crust. We find that the deviatoric stress exerted by topography is minimized in the limit when the crust is considered an incompressible elastic solid, with a Poisson ratio of 0.5, and is independent of Young's modulus. This places a strict lower bound on the critical stress state maintained by the crust supporting plastically deformed accretionary wedge topography. We estimate the coseismic shear stress change from the Maule event ranged from −6 MPa (stress increase) to 17 MPa (stress drop), with a maximum depth‐averaged crustal shear‐stress drop of 4 MPa. We separately estimate that the plate‐driving forces acting in the region, regardless of their exact mechanism, must contribute at least 27 MPa trench‐perpendicular compression and 15 MPa trench‐parallel compression. This corresponds to a depth‐averaged shear stress of at least 7 MPa. The comparable magnitude of these two independent shear stress estimates is consistent with the interpretation that the section of the megathrust fault ruptured in the Maule earthquake is weak, with the seismic cycle relieving much of the total sustained shear stress in the crust.

Maule↗

USGS2008

Meteorites on Mars observed with Mars Exploration Rovers

Reduced weathering rates due to the lack of liquid water and significantly greater typical surface ages should result in a higher density of meteorites on the surface of Mars compared to Earth. Several meteorites were identified among the rocks investigated during Opportunity's traverse across the sandy Meridiani plains. Heat Shield Rock is a IAB iron meteorite and has been officially recognized as 'Meridiani Planum.' Barberton is olivine-rich and contains metallic Fe in the form of kamacite, suggesting a meteoritic origin. It is chemically most consistent with a mesosiderite silicate clast. Santa Catarina is a brecciated rock with a chemical and mineralogical composition similar to Barberton. Barberton, Santa Catarina, and cobbles adjacent to Santa Catarina may be part of a strewn field. Spirit observed two probable iron meteorites from its Winter Haven location in the Columbia Hills in Gusev Crater. Chondrites have not been identified to date, which may be a result of their lower strengths and probability to survive impact at current atmospheric pressures. Impact craters directly associated with Heat Shield Rock, Barberton, or Santa Catarina have not been observed, but such craters could have been erased by eolian-driven erosion.

Journal of Geophysical Research E: Planets↗

USGS1999

Stress transfer by the 1988-1989 M=5.3 and 5.4 Lake Elsman foreshocks to the Loma Prieta fault: Unclamping at the site of peak mainshock slip

We study the stress transferred by the June 27, 1988, M = 5.3 and August 8, 1989, M = 5.4 Lake Elsman earthquakes, the largest events to strike within 15 km of the future Loma Prieta rupture zone during 74 years before the 1989 M = 6.9 Loma Prieta earthquake. We find that the first Lake Elsman event brought the rupture plane of the second event 0.3–1.6 bars (0.03–0.16 MPa) closer to Coulomb failure but that the Lake Elsman events did not bring the future Loma Prieta hypocentral zone closer to failure. Instead, the Lake Elsman earthquakes are calculated to have reduced the normal stress on (or “undamped”) the Loma Prieta rupture surface by 0.5–1.0 bar (0.05–0.10 MPa) at the site where the greatest slip subsequently occurred in the Loma Prieta earthquake. This association between the sites of peak unclamping and slip suggests that the Lake Elsman events did indeed influence the Loma Prieta rupture process. Unclamping the fault would have locally lowered the resistance to sliding. Such an effect could have been enhanced if the lowered normal stress permitted fluid infusion into the undamped part of the fault. Although less well recorded, the M L = 5.0 1964 and M L = 5.3 1967 Corralitos events struck within 10 km of the southwest end of the future Loma Prieta rupture. No similar relationship between the normal stress change and subsequent Loma Prieta slip is observed, although the high-slip patch southwest of the Loma Prieta epicenter corresponds roughly to the site of calculated Coulomb stress increase for a low coefficient of friction. The Lake Elsman-Loma Prieta result is similar to that for the 1987 M = 6.2 Elmore Ranch and M = 6.7 Superstition Hills earthquakes, suggesting that foreshocks might influence the distribution of mainshock slip rather than the site of mainshock nucleation.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Field Reconnaissance Geologic Mapping of the Columbia Hills, Mars: Results from MER Spirit and MRO HiRISE Observations

Chemical, mineralogic, and lithologic ground truth was acquired for the first time on Mars in terrain units mapped using orbital Mars Reconnaissance Orbiter's High Resolution Imaging Science Experiment (MRO HiRISE) image data. Examination of several dozen outcrops shows that Mars is geologically complex at meter length scales, the record of its geologic history is well exposed, stratigraphic units may be identified and correlated across significant areas on the ground, and outcrops and geologic relationships between materials may be analyzed with techniques commonly employed in terrestrial field geology. Despite their burial during the course of Martian geologic time by widespread epiclastic materials, mobile fines, and fall deposits, the selective exhumation of deep and well-preserved geologic units has exposed undisturbed outcrops, stratigraphic sections, and structural information much as they are preserved and exposed on Earth. A rich geologic record awaits skilled future field investigators on Mars. The correlation of ground observations and orbital images enables construction of a corresponding geologic reconnaissance map. Most of the outcrops visited are interpreted to be pyroclastic, impactite, and epiclastic deposits overlying an unexposed substrate, probably related to a modified Gusev crater central peak. Fluids have altered chemistry and mineralogy of these protoliths in degrees that vary substantially within the same map unit. Examination of the rocks exposed above and below the major unconformity between the plains lavas and the Columbia Hills directly confirms the general conclusion from remote sensing in previous studies over past years that the early history of Mars was a time of more intense deposition and modification of the surface. Although the availability of fluids and the chemical and mineral activity declined from this early period, significant later volcanism and fluid convection enabled additional, if localized, chemical activity.

Journal of Geophysical Research E: Planets↗

USGS2001

Seismic mapping of shallow fault zones in the San Gabriel Mountains from the Los Angeles Region Seismic Experiment, southern California

During the Los Angeles Region Seismic Experiment (LARSE), a reflection/refraction survey was conducted along a profile (line 1) extending from Seal Beach, California, northeastward to the Mojave Desert and crossing the Los Angeles and San Gabriel Valley basins and San Gabriel Mountains. In most shot gathers from the southern and central San Gabriel Mountains, clear secondary arrivals are seen that merge, or appear to merge, with first arrivals at three locations, including the location of the Vincent thrust fault, an exposed late Mesozoic/early Cenozoic megathrust. These secondary arrivals are interpretable as reflections in the shallow crust (<5 km depth) from a concave-upward interface that projects to the surface in the north near the Vincent thrust fault, is offset in its central part at the San Gabriel fault (an old branch of the San Andreas fault), and terminates in the south at 1 to 2 km depth at the southern mountain front. The velocity structure above and below this interface strongly suggests it is the Vincent thrust fault: intermediate velocities (6.2 km/s), consistent with mylonites overlying the Vincent thrust fault, are observed above it; lower velocities (5.8 km/s), consistent with the Pelona Schist underlying the Vincent thrust fault, are observed below it. Problems arise, however, in attempting to match this reflector to the exposed Vincent thrust fault, which is seen in outcrops east of line 1. The Vincent thrust fault is shallower than the reflector in most places. An unmapped structure (steep fault, monocline, or thrust fault) is required between line 1 and the outcrops that either drops the Vincent thrust fault down to the depths of the reflector or repeats the Vincent thrust fault beneath line 1 in the footwall of another thrust fault. An alternative interpretation of the reflector is a deep greenstone horizon within the Pelona Schist, although this alternative is not favored by the velocity structure. Copyright 2001 by the American Geophysical Union.

California↗

USGS2012

Hydrologic effects on water level changes associated with episodic fault creep near Parkfield, California

As part of the Parkfield, California, earthquake prediction experiment, water level is monitored in a well 460 m from the main trace of the San Andreas fault on Middle Mountain, in the preparation zone of the anticipated Parkfield earthquake. The well configuration allows water level to be monitored in two fluid reservoirs at depths of 85 and 250 m below land surface. During 1987, water level changes were recorded during 12 of the 18 episodes of accelerated fault creep detected by a creep meter spanning the fault trace 750 m northwest of the well. The creep-related water level changes in the shallow reservoir have durations of less than 1 day, whereas in the deeper reservoir the changes persist for as long as 2 months. These data suggest that the transient nature of the water level changes in the shallow interval is due to vertical flow to the water table and is not evidence that creep events propagate past the well. Phase leads of earth tidal constituents in the water level data from the shallow interval relative to the same constituents in the local volume strain tide support the interpretation of significant flow to the water table at periods of 1 day or less. The form of the water level changes in the deep interval is affected by horizontal flow to the well bore. This effect can be removed from the water level records using a theoretical response curve constrained by the phases of earth tidal constituents in the deep interval relative to the local volume strain tide. For the events where the signal in the shallow interval has been large enough to measure, the sizes of the simultaneous water level changes in the two reservoirs are consistent with the same amounts of volume strain occurring at both depths.

Journal of Geophysical Research Solid Earth↗

USGS2016

Insights into shallow magmatic processes at Kīlauea Volcano, Hawaiʻi, from a multiyear continuous gravity time series

Continuous gravity data collected near the summit eruptive vent at Kīlauea Volcano, Hawaiʻi, during 2011–2015 show a strong correlation with summit-area surface deformation and the level of the lava lake within the vent over periods of days to weeks, suggesting that changes in gravity reflect variations in volcanic activity. Joint analysis of gravity and lava level time series data indicates that over the entire time period studied, the average density of the lava within the upper tens to hundreds of meters of the summit eruptive vent remained low—approximately 1000–1500 kg/m 3 . The ratio of gravity change (adjusted for Earth tides and instrumental drift) to lava level change measured over 15 day windows rose gradually over the course of 2011–2015, probably reflecting either (1) a small increase in the density of lava within the eruptive vent or (2) an increase in the volume of lava within the vent due to gradual vent enlargement. Superimposed on the overall time series were transient spikes of mass change associated with inflation and deflation of Kīlauea's summit and coincident changes in lava level. The unexpectedly strong mass variations during these episodes suggest magma flux to and from the shallow magmatic system without commensurate deformation, perhaps indicating magma accumulation within, and withdrawal from, void space—a process that might not otherwise be apparent from lava level and deformation data alone. Continuous gravity data thus provide unique insights into magmatic processes, arguing for continued application of the method at other frequently active volcanoes.

Hawaii↗

USGS2010

Rheologic and structural controls on the deformation of Okmok volcano, Alaska: FEMs, InSAR, and ambient noise tomography

Interferometric synthetic aperture radar (InSAR) data indicate that the caldera of Okmok volcano, Alaska, subsided more than a meter during its eruption in 1997. The large deformation suggests a relatively shallow magma reservoir beneath Okmok. Seismic tomography using ambient ocean noise reveals two low‐velocity zones (LVZs). The shallow LVZ corresponds to a region of weak, fluid‐saturated materials within the caldera and extends from the caldera surface to a depth of 2 km. The deep LVZ clearly indicates the presence of the magma reservoir beneath Okmok that is significantly deeper (>4 km depth) compared to previous geodetic‐based estimates (3 km depth). The deep LVZ associated with the magma reservoir suggests magma remains in a molten state between eruptions. We construct finite element models (FEMs) to simulate deformation caused by mass extraction from a magma reservoir that is surrounded by a viscoelastic rind of country rock embedded in an elastic domain that is partitioned to account for the weak caldera materials observed with tomography. This configuration allows us to reduce the estimated magma reservoir depressurization to within lithostatic constraints, while simultaneously maintaining the magnitude of deformation required to predict the InSAR data. More precisely, the InSAR data are best predicted by an FEM simulating a rind viscosity of 7.5 × 10 16 Pa s and a mass flux of −4.2 × 10 9 kg/d from the magma reservoir. The shallow weak layer within the caldera provides a coeruption stress regime and neutral buoyancy horizon that support lateral magma propagation from the central magma reservoir to extrusion near the rim of the caldera.

Alaska↗

USGS1997

Färoe-Iceland Ridge Experiment: 1. Crustal structure of northeastern Iceland

Results from the Färoe-Iceland Ridge Experiment (FIRE) constrain the crustal thickness as 19 km under the Northern Volcanic Zone of Iceland and 35 km under older Tertiary areas of northeastern Iceland. The Moho is defined by strong P wave and S wave reflections. Synthetic seismogram modeling of the Moho reflection indicates mantle velocities of at least 8.0 km/s beneath the Tertiary areas of northeastern Iceland and at least 7.9 km/s beneath the neovolcanic zone. Crustal diving rays resolve the structure of the upper and lower crust. Surface P wave velocities are 1.1–4.0 km/s in Quaternary rocks and are rather higher, 4.4–4.7 km/s, in the Tertiary basalts that outcrop elsewhere. The highest crustal P wave velocities observed directly from diving rays are 7.1 km/s, from rays that turn at 24 km depth. Velocities of 7.35 km/s at the base of the crust are inferred from extrapolation of the lower crustal velocity gradient (0.024 s −1 ). A Poisson's ratio of approximately 0.27, equivalent to an S wave to P wave travel time ratio of 1.78, is measured throughout the crust east of the neovolcanic zone. The Poisson's ratio and the steep Moho topography (in places up to 30° from the horizontal) indicate that the entire crust outside the neovolcanic zone is cool (<800°C). Gravity data are well matched by a velocity/density conversion of our seismic crustal model and indicate a region of low mantle density beneath the neovolcanic zone, believed to be due to elevated mantle temperatures. The crustal thickness in the neovolcanic zone is consistent with geochemical estimates of the melt generation, placing constraints on the flow within the Iceland mantle plume.

Journal of Geophysical Research B: Solid Earth↗

USGS2017

Integrated geophysical characteristics of the 2015 Illapel, Chile, earthquake

On 16 September 2015, a Mw 8.3 earthquake ruptured the subduction zone offshore of Illapel, Chile, generating an aftershock sequence with 14 Mw 6.0–7.0 events. A double source W phase moment tensor inversion consists of a Mw 7.2 subevent and the main Mw 8.2 phase. We determine two slip models for the mainshock, one using teleseismic broadband waveforms and the other using static GPS and InSAR surface displacements, which indicate high slip north of the epicenter and west-northwest of the epicenter near the oceanic trench. These models and slip distributions published in other studies suggest spatial slip uncertainties of ~25 km and have peak slip values that vary by a factor of 2. We relocate aftershock hypocenters using a Bayesian multiple-event relocation algorithm, revealing a cluster of aftershocks under the Chilean coast associated with deep (20–45 km depth) mainshock slip. Less vigorous aftershock activity also occurred near the trench and along strike of the main aftershock region. Most aftershocks are thrust-faulting events, except for normal-faulting events near the trench. Coulomb failure stress change amplitudes and signs are uncertain for aftershocks collocated with deeper mainshock slip; other aftershocks are more clearly associated with loading from the mainshock. These observations reveal a frictionally heterogeneous interface that ruptured in patches at seismogenic depths (associated with many aftershocks) and with homogeneous slip (and few aftershocks) up to the trench. This event likely triggered seismicity separate from the main slip region, including along-strike events on the megathrust and intraplate extensional events.

Illapel↗