Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research, Solid Earth”

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

Magma flux at Okmok Volcano, Alaska, from a joint inversion of continuous GPS, campaign GPS, and interferometric synthetic aperture radar

Volcano deformation is usually measured using satellite geodetic techniques including interferometric synthetic aperture radar (InSAR), campaign GPS, and continuous GPS. Differences in the spatial and temporal sampling of each system mean that most appropriate inversion scheme to determine the source parameters from each data set is different. Most studies either compare results from independent inversions or subsample the data sets to the lowest common factor. It is unclear whether differences in the solution reflect differences in source behavior, differences in measurement bias, or differences in inversion technique. Here we develop a single inversion procedure that captures the benefits of each system, especially the daily sampling of continuous GPS and the high spatial resolution of InSAR. Okmok Volcano, Alaska, is an ideal target for such a test because a long series (<15 years) of InSAR and continuous GPS measurement exists and the source is almost continuously active and in a stable location.

Alaska↗
USGS2010

Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly-heterogeneous crustal stress

It has been proposed that the crustal stress field contains small-length-scale heterogeneity of much larger amplitude than the uniform background stress. This model predicts that earthquake focal mechanisms should reflect the loading stress rather than the uniform background stress. So, if the heterogeneous stress hypothesis is correct, focal mechanisms before and after a large earthquake should align with the tectonic loading and the earthquake-induced static stress perturbation, respectively. However, I show that the off-fault triggered aftershocks of the 1992 M7.3 Landers, California, earthquake align with the same stress field as the pre-Landers mechanisms. The aftershocks occurred on faults that were well oriented for failure in the pre-Landers stress field and then loaded by the Landers-induced static stress change. Aftershocks in regions experiencing a 0.05 to 5 MPa coseismic differential stress change align with the modeled Landers-induced static stress change, implying that they were triggered by the stress perturbation. Contrary to the heterogeneous stress hypothesis, these triggered aftershocks are also well aligned with the pre-Landers stress field obtained from inverting the pre-Landers focal mechanisms. Therefore, the inverted pre-Landers stress must represent the persistent background stress field. Earthquake focal mechanisms provide an unbiased sample of the spatially coherent background stress field, which is large relative to any small-scale stress heterogeneity. The counterexample provided by the Landers earthquake is strong evidence that the heterogeneous stress model is not widely applicable.

California↗
USGS2011

Stress transfer among en echelon and opposing thrusts and tear faults: Triggering caused by the 2003 M w = 6.9 Zemmouri, Algeria, earthquake

The essential features of stress interaction among earthquakes on en echelon thrusts and tear faults were investigated, first through idealized examples and then by study of thrust faulting in Algeria. We calculated coseismic stress changes caused by the 2003 M w = 6.9 Zemmouri earthquake, finding that a large majority of the Zemmouri afterslip sites were brought several bars closer to Coulomb failure by the coseismic stresses, while the majority of aftershock nodal planes were brought closer to failure by an average of ~2 bars. Further, we calculated that the shallow portions of the adjacent Thenia tear fault, which sustained ~0.25 m slip, were brought >2 bars closer to failure. We calculated that the Coulomb stress increased by 1.5 bars on the deeper portions of the adjacent Boumerdes thrust, which lies just 10–20 km from the city of Algiers; both the Boumerdes and Thenia faults were illuminated by aftershocks. Over the next 6 years, the entire south dipping thrust system extending 80 km to the southwest experienced an increased rate of seismicity. The stress also increased by 0.4 bar on the east Sahel thrust fault west of the Zemmouri rupture. Algiers suffered large damaging earthquakes in A.D. 1365 and 1716 and is today home to 3 million people. If these shocks occurred on the east Sahel fault and if it has a ~2 mm/yr tectonic loading rate, then enough loading has accumulated to produce a M w = 6.6–6.9 shock today. Thus, these potentially lethal faults need better understanding of their slip rate and earthquake history.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Postseismic relaxation following the 1994 M w 6.7 Northridge earthquake, southern California

We have reexamined the postearthquake deformation of a 65 km long linear array of 11 geodetic monuments extending north–south across the rupture (reverse slip on a blind thrust dipping 40°S–20°W) associated with the 1994 M w 6.7 Northridge earthquake. That array was surveyed frequently in the interval from 4 to 2650 days after the earthquake. The velocity of each of the monuments over the interval 100–2650 days postearthquake appears to be constant. Moreover, the profile of those velocities along the length of the array is very similar to a preearthquake velocity profile for a nearby, similarly oriented array. We take this to indicate that significant postseismic relaxation is evident only in the first 100 days postseismic and that the subsequent linear trend is typical of the interseismic interval. The postseismic relaxation (postseismic displacement less displacement that would have occurred at the preseismic velocity) is found to be almost wholly parallel (N70°W) to the nearby (40 km) San Andreas Fault with only negligible relaxation in the direction of coseismic slip (N20°E) on the Northridge rupture. We suggest that the N70°W relaxation is caused by aseismic, right-lateral slip at depth on the San Andreas Fault, excess slip presumably triggered by the Northridge rupture. Finally, using the Dieterich (1994) stress-seismicity relation, we show that return to the preseismic deformation rate within 100 days following the earthquake could be consistent with the cumulative number of M > 2.5 earthquakes observed following the main shock.

California↗
USGS2011

Toward a consistent model for strain accrual and release for the New Madrid Seismic Zone, central United States

At the heart of the conundrum of seismogenesis in the New Madrid Seismic Zone is the apparently substantial discrepancy between low strain rate and high recent seismic moment release. In this study we revisit the magnitudes of the four principal 1811–1812 earthquakes using intensity values determined from individual assessments from four experts. Using these values and the grid search method of Bakun and Wentworth (1997), we estimate magnitudes around 7.0 for all four events, values that are significantly lower than previously published magnitude estimates based on macroseismic intensities. We further show that the strain rate predicted from postglacial rebound is sufficient to produce a sequence with the moment release of one M max 6.8 every 500 years, a rate that is much lower than previous estimates of late Holocene moment release. However, M w 6.8 is at the low end of the uncertainty range inferred from analysis of intensities for the largest 1811–1812 event. We show that M w 6.8 is also a reasonable value for the largest main shock given a plausible rupture scenario. One can also construct a range of consistent models that permit a somewhat higher M max , with a longer average recurrence rate. It is thus possible to reconcile predicted strain and seismic moment release rates with alternative models: one in which 1811–1812 sequences occur every 500 years, with the largest events being M max ∼6.8, or one in which sequences occur, on average, less frequently, with Mmax of ∼7.0. Both models predict that the late Holocene rate of activity will continue for the next few to 10 thousand years.

New Madrid Seismic Zone↗
USGS2011

Frictional strengths of talc-serpentine and talc-quartz mixtures

Talc is a constituent of faults in a variety of settings, and it may be an effective weakening agent depending on its abundance and distribution within a fault. We conducted frictional strength experiments under hydrothermal conditions to determine the effect of talc on the strengths of synthetic gouges of lizardite and antigorite serpentinites and of quartz. Small amounts of talc weaken serpentinite gouges substantially more than predicted by simple weight averaging. In comparison, mixtures of quartz and talc show a linear trend of strength reduction at talc concentrations 15 wt % and enhanced weakening at higher concentrations. All of the strength data are fit by a modified version of the Reuss mixing law that allows for the dominance of one mineral over the other. The difference in the behavior of serpentinite-talc and quartz-talc mixtures at low talc concentrations is a reflection of their different textures. Lizardite, antigorite, and talc all have platy habits, and displacement within gouges composed of these minerals is localized to narrow shears along which the platy grains have rotated into alignment with the shear surfaces. The shears in the mixed phyllosilicate gouges maximize the proportion of the weaker mineral within them. When mixed with a strong, rounded mineral such as quartz, some minimum concentration of talc is needed to form connected pathways that enhance strength reductions. The typical development of talc by the reaction of Si-rich fluids with serpentinite or dolomite would tend to localize its occurrence in a natural fault and result in enhanced weakening.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Multiplets: Their behavior and utility at dacitic and andesitic volcanic centers

Multiplets, or groups of earthquakes with similar waveforms, are commonly observed at volcanoes, particularly those exhibiting unrest. Using triggered seismic data from the 1980-1986 Mount St. Helens (MSH) eruption, we have constructed a catalog of multiplet occurrence. Our analysis reveals that the occurrence of multiplets is related, at least in part, to the viscosity of the magma. We also constructed catalogs of multiplet occurrence using continuous seismic data from the 2004 eruption at MSH and 2007 eruption at Bezymianny Volcano, Russia. Prior to explosions at MSH in 2004 and Bezymianny in 2007, the multiplet proportion of total seismicity (MPTS) declined, while the average amplitudes and standard deviations of the average amplitude increased. The life spans of multiplets (time between the first and last event) were also shorter prior to explosions than during passive lava extrusion. Dome-forming eruptions that include a partially solidified plug, like MSH (1983-1986 and 2004-2008), often possess multiplets with longer life spans and MPTS values exceeding 50%. Conceptually, the relatively unstable environment prior to explosions is characterized by large and variable stress gradients brought about by rapidly changing overpressures within the conduit. We infer that such complex stress fields affect the number of concurrent families, MPTS, average amplitude, and standard deviation of the amplitude of the multiplets. We also argue that multiplet detection may be an important new monitoring tool for determining the timing of explosions and in forecasting the type of eruption.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

The magnitude distribution of earthquakes near Southern California faults

We investigate seismicity near faults in the Southern California Earthquake Center Community Fault Model. We search for anomalously large events that might be signs of a characteristic earthquake distribution. We find that seismicity near major fault zones in Southern California is well modeled by a Gutenberg-Richter distribution, with no evidence of characteristic earthquakes within the resolution limits of the modern instrumental catalog. However, the b value of the locally observed magnitude distribution is found to depend on distance to the nearest mapped fault segment, which suggests that earthquakes nucleating near major faults are likely to have larger magnitudes relative to earthquakes nucleating far from major faults. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Elevated gas hydrate saturation within silt and silty clay sediments in the Shenhu area, South China Sea

Gas hydrate saturations were estimated using five different methods in silt and silty clay foraminiferous sediments from drill hole SH2 in the South China Sea. Gas hydrate saturations derived from observed pore water chloride values in core samples range from 10 to 45% of the pore space at 190–221 m below seafloor (mbsf). Gas hydrate saturations estimated from resistivity (R t ) using wireline logging results are similar and range from 10 to 40.5% in the pore space. Gas hydrate saturations were also estimated by P wave velocity obtained during wireline logging by using a simplified three-phase equation (STPE) and effective medium theory (EMT) models. Gas hydrate saturations obtained from the STPE velocity model (41.0% maximum) are slightly higher than those calculated with the EMT velocity model (38.5% maximum). Methane analysis from a 69 cm long depressurized core from the hydrate-bearing sediment zone indicates that gas hydrate saturation is about 27.08% of the pore space at 197.5 mbsf. Results from the five methods show similar values and nearly identical trends in gas hydrate saturations above the base of the gas hydrate stability zone at depths of 190 to 221 mbsf. Gas hydrate occurs within units of clayey slit and silt containing abundant calcareous nannofossils and foraminifer, which increase the porosities of the fine-grained sediments and provide space for enhanced gas hydrate formation. In addition, gas chimneys, faults, and fractures identified from three-dimensional (3-D) and high-resolution two-dimensional (2-D) seismic data provide pathways for fluids migrating into the gas hydrate stability zone which transport methane for the formation of gas hydrate. Sedimentation and local canyon migration may contribute to higher gas hydrate saturations near the base of the stability zone.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

A reevaluation of the Pallett Creek earthquake chronology based on new AMS radiocarbon dates, San Andreas fault, California

The Pallett Creek paleoseismic record occupies a keystone position in most attempts to develop rupture histories for the southern San Andreas fault. Previous estimates of earthquake ages at Pallett Creek were determined by decay counting radiocarbon methods. That method requires large samples which can lead to unaccounted sources of uncertainty in radiocarbon ages because of the heterogeneous composition of organic layers. In contrast, accelerator mass spectrometry (AMS) radiocarbon dates may be obtained from small samples that have known carbon sources and also allow for a more complete sampling of the section. We present 65 new AMS radiocarbon dates that span nine ground‐rupturing earthquakes at Pallett Creek. Overall, the AMS dates are similar to and reveal no dramatic bias in the conventional dates. For many layers, however, individual charcoal samples were younger than the conventional dates, leading to earthquake ages that are overall slightly younger than previously reported. New earthquake ages are determined by Bayesian refinement of the layer ages based on stratigraphic ordering and sedimentological constraints. The new chronology is more regular than previously published records in large part due to new samples constraining the age of event R. The closed interval from event C to 1857 has a mean recurrence of 135 years ( σ = 83.2 years) and a quasiperiodic coefficient of variation (COV) of 0.61. We show that the new dates and resultant earthquake chronology have a stronger effect on COV than the specific membership of this long series and dating precision improvements from sedimentation rates.

California↗
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↗
USGS2011

A model for seasonal changes in GPS positions and seismic wave speeds due to thermoelastic and hydrologic variations

It is known that GPS time series contain a seasonal variation that is not due to tectonic motions, and it has recently been shown that crustal seismic velocities may also vary seasonally. In order to explain these changes, a number of hypotheses have been given, among which thermoelastic and hydrology-induced stresses and strains are leading candidates. Unfortunately, though, since a general framework does not exist for understanding such seasonal variations, it is currently not possible to quickly evaluate the plausibility of these hypotheses. To fill this gap in the literature, I generalize a two-dimensional thermoelastic strain model to provide an analytic solution for the displacements and wave speed changes due to either thermoelastic stresses or hydrologic loading, which consists of poroelastic stresses and purely elastic stresses. The thermoelastic model assumes a periodic surface temperature, and the hydrologic models similarly assume a periodic near-surface water load. Since all three models are two-dimensional and periodic, they are expected to only approximate any realistic scenario; but the models nonetheless provide a quantitative framework for estimating the effects of thermoelastic and hydrologic variations. Quantitative comparison between the models and observations is further complicated by the large uncertainty in some of the relevant parameters. Despite this uncertainty, though, I find that maximum realistic thermoelastic effects are unlikely to explain a large fraction of the observed annual variation in a typical GPS displacement time series or of the observed annual variations in seismic wave speeds in southern California. Hydrologic loading, on the other hand, may be able to explain a larger fraction of both the annual variations in displacements and seismic wave speeds. Neither model is likely to explain all of the seismic wave speed variations inferred from observations. However, more definitive conclusions cannot be made until the model parameters are better constrained. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Seismic and geodetic signatures of fault slip at the Slumgullion Landslide Natural Laboratory

We tested the hypothesis that the Slumgullion landslide is a useful natural laboratory for observing fault slip, specifically that slip along its basal surface and side-bounding strike-slip faults occurs with comparable richness of aseismic and seismic modes as along crustal- and plate-scale boundaries. Our study provides new constraints on models governing landslide motion. We monitored landslide deformation with temporary deployments of a 29-element prism array surveyed by a robotic theodolite and an 88-station seismic network that complemented permanent extensometers and environmental instrumentation. Aseismic deformation observations show that large blocks of the landslide move steadily at approximately centimeters per day, possibly punctuated by variations of a few millimeters, while localized transient slip episodes of blocks less than a few tens of meters across occur frequently. We recorded a rich variety of seismic signals, nearly all of which originated outside the monitoring network boundaries or from the side-bounding strike-slip faults. The landslide basal surface beneath our seismic network likely slipped almost completely aseismically. Our results provide independent corroboration of previous inferences that dilatant strengthening along sections of the side-bounding strike-slip faults controls the overall landslide motion, acting as seismically radiating brakes that limit acceleration of the aseismically slipping basal surface. Dilatant strengthening has also been invoked in recent models of transient slip and tremor sources along crustal- and plate-scale faults suggesting that the landslide may indeed be a useful natural laboratory for testing predictions of specific mechanisms that control fault slip at all scales.

Colorado↗
USGS2011

InSAR observations of aseismic slip associated with an earthquake swarm in the Columbia River flood basalts

In 2009 a swarm of small shallow earthquakes occurred within the basalt flows of the Columbia River Basalt Group (CRBG). The swarm occurred within a dense seismic network in the U.S. Department of Energys Hanford Site. Data from the seismic network along with interferometric synthetic aperture radar (InSAR) data from the European Space Agencys (ESA) ENVISAT satellite provide insight into the nature of the swarm. By modeling the InSAR deformation data we constructed a model that consists of a shallow thrust fault and a near horizontal fault. We suggest that the near horizontal lying fault is a bedding-plane fault located between basalt flows. The geodetic moment of the modeled fault system is about eight times the cumulative seismic moment of the swarm. Precise location estimates of the swarm earthquakes indicate that the area of highest slip on the thrust fault, ???70mm of slip less than ???0.5km depth, was not located within the swarm cluster. Most of the slip on the faults appears to have progressed aseismically and we suggest that interbed sediments play a central role in the slip process. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Historical perspective on seismic hazard to Hispaniola and the northeast Caribbean region

We evaluate the long-term seismic activity of the North-American/Caribbean plate boundary from 500 years of historical earthquake damage reports. The 2010 Haiti earthquakes and other earthquakes were used to derive regional attenuation relationships between earthquake intensity, magnitude, and distance from the reported damage to the epicenter, for Hispaniola and for Puerto Rico and the Virgin Islands. The attenuation relationship for Hispaniola earthquakes and northern Lesser Antilles earthquakes is similar to that for California earthquakes, indicating a relatively rapid attenuation of damage intensity with distance. Intensities in Puerto Rico and the Virgin Islands decrease less rapidly with distance. We use the intensity-magnitude relationships to systematically search for the location and intensity magnitude MI which best fit all the reported damage for historical earthquakes. Many events occurred in the 20th-century along the plate-boundary segment from central Hispaniola to the NW tip of Puerto Rico, but earlier events from this segment were not identified. The remaining plate boundary to the east to Guadeloupe is probably not associated with M > 8 historical subduction-zone earthquakes. The May 2, 1787 earthquake, previously assigned an M 8–8.25, is probably only M I 6.9 and could be located north, west or SW of Puerto Rico. An M I 6.9 earthquake on July 11, 1785 was probably located north or east of the Virgin Islands. We located M I < 8 historical earthquakes on April 5, 1690, February 8, 1843, and October 8, 1974 in the northern Lesser Antilles within the arc. We speculate that the December 2, 1562 (M I 7.7) and May 7, 1842 (M I 7.6) earthquakes ruptured the Septentrional Fault in northern Hispaniola. If so, the recurrence interval on the central Septentrional Fault is ~300 years, and only 170 years has elapsed since the last event. The recurrence interval of large earthquakes along the Hispaniola subduction segment is likely longer than the historical record. Intra-arc M ≥ 7.0 earthquakes may occur every 75–100 years in the 410-km-long segment between the Virgin Islands and Guadeloupe.

Caribbean Sea↗
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↗
USGS2012

Slab1.0: A three-dimensional model of global subduction zone geometries

We describe and present a new model of global subduction zone geometries, called Slab1.0. An extension of previous efforts to constrain the two-dimensional non-planar geometry of subduction zones around the focus of large earthquakes, Slab1.0 describes the detailed, non-planar, three-dimensional geometry of approximately 85% of subduction zones worldwide. While the model focuses on the detailed form of each slab from their trenches through the seismogenic zone, where it combines data sets from active source and passive seismology, it also continues to the limits of their seismic extent in the upper-mid mantle, providing a uniform approach to the definition of the entire seismically active slab geometry. Examples are shown for two well-constrained global locations; models for many other regions are available and can be freely downloaded in several formats from our new Slab1.0 website, http://on.doi.gov/d9ARbS . We describe improvements in our two-dimensional geometry constraint inversion, including the use of ‘average’ active source seismic data profiles in the shallow trench regions where data are otherwise lacking, derived from the interpolation between other active source seismic data along-strike in the same subduction zone. We include several analyses of the uncertainty and robustness of our three-dimensional interpolation methods. In addition, we use the filtered, subduction-related earthquake data sets compiled to build Slab1.0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far, concluding that the width of these seismogenic zones is on average 30% larger than previous studies have suggested.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Shallow conduit system at Kilauea Volcano, Hawaii, revealed by seismic signals associated with degassing bursts

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in March, 2008 and continuing to the present time is characterized by episodic explosive bursts of gas and ash from a vent within Halemaumau Pit Crater. These bursts are accompanied by seismic signals that are well recorded by a broadband network deployed in the summit caldera. We investigate in detail the dimensions and oscillation modes of the source of a representative burst in the 1−10 s band. An extended source is realized by a set of point sources distributed on a grid surrounding the source centroid, where the centroid position and source geometry are fixed from previous modeling of very-long-period (VLP) data in the 10–50 s band. The source time histories of all point sources are obtained simultaneously through waveform inversion carried out in the frequency domain. Short-scale noisy fluctuations of the source time histories between adjacent sources are suppressed with a smoothing constraint, whose strength is determined through a minimization of the Akaike Bayesian Information Criterion (ABIC). Waveform inversions carried out for homogeneous and heterogeneous velocity structures both image a dominant source component in the form of an east trending dike with dimensions of 2.9 × 2.9 km. The dike extends ∼2 km west and ∼0.9 km east of the VLP centroid and spans the depth range 0.2–3.1 km. The source model for a homogeneous velocity structure suggests the dike is hinged at the source centroid where it bends from a strike E 27°N with northern dip of 85° west of the centroid, to a strike E 7°N with northern dip of 80° east of the centroid. The oscillating behavior of the dike is dominated by simple harmonic modes with frequencies ∼0.2 Hz and ∼0.5 Hz, representing the fundamental mode ν 11 and first degenerate mode ν 12 = ν 21 of the dike. Although not strongly supported by data in the 1–10 s band, a north striking dike segment is required for enhanced compatibility with the model elaborated in the 10–50 s band. This dike provides connectivity between the east trending dike and the new vent within Halemaumau Pit Crater. Waveform inversions with a dual-dike model suggest dimensions of 0.7 × 0.7 km to 2.6 × 2.6 km for this segment. Further elaboration of the complex dike system under Halemaumau does not appear to be feasible with presently available data.

Hawai'i↗