SEARCH · USGS Science

Results for “Earth Observer”

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 1,549 records · Page 86Linked to original sources

USGS2025

Multi-scale intense geoelectric and geomagnetic field perturbations observed after an interplanetary magnetic field turning

Intense geoelectric fields during geomagnetic storms generate geomagnetically induced currents in power grids and other infrastructure, necessitating an understanding of their causes, for example, through coordinated space and ground observations. This study investigates localized intense geoelectric ( E ) and geomagnetic ( B ) field perturbations following an Interplanetary Magnetic Field (IMF) turning during a geomagnetic storm on 25 October 2011. Observations from EarthScope magnetotelluric sites in the upper Midwest United States revealed shorter period ( ∼ 1 min) ultra-low-frequency (ULF) waves superimposed on longer period ( ∼ 10 min) perturbations in both E and B fields. These sites, located at ∼ 19 hr magnetic local time and 56 −57 ° magnetic latitude, recorded large amplitude E and B perturbations. Ground-based all-sky imagers showed auroral brightening with sunward and poleward propagation, while upstream spacecraft linked the perturbations to an IMF turning and solar wind dynamic pressure impulse. The longer-period E and B field perturbations likely stem from localized ionospheric currents tied to substorm auroral activity post-IMF turning. The combination of ionospheric currents, ULF waves, and the Earth's varying conductivity produces intense geoelectric fields of ≥ 2 V/km in the upper Midwest. A comparison using input data and software compatible with the NOAA/USGS geoelectric field nowcast model revealed its limitations in capturing such events due to the temporal and spatial resolution of the underlying data. Using 1-s geomagnetic field data can improve geoelectric field models by capturing short-period and large spatial scale waves, although localized magnetic perturbations remain underestimated due to insufficient ground magnetometer density.

Space Weather↗

USGS2024

A conserved interdomain microbial network underpins cadaver decomposition despite environmental variables

Microbial breakdown of organic matter is one of the most important processes on Earth, yet the controls of decomposition are poorly understood. Here we track 36 terrestrial human cadavers in three locations and show that a phylogenetically distinct, interdomain microbial network assembles during decomposition despite selection effects of location, climate and season. We generated a metagenome-assembled genome library from cadaver-associated soils and integrated it with metabolomics data to identify links between taxonomy and function. This universal network of microbial decomposers is characterized by cross-feeding to metabolize labile decomposition products. The key bacterial and fungal decomposers are rare across non-decomposition environments and appear unique to the breakdown of terrestrial decaying flesh, including humans, swine, mice and cattle, with insects as likely important vectors for dispersal. The observed lockstep of microbial interactions further underlies a robust microbial forensic tool with the potential to aid predictions of the time since death.

Nature Microbiology↗

USGS2017

Using mineral geochemistry to decipher slab, mantle, and crustal input in the generation of high-Mg andesites and basaltic andesites from the northern Cascade Arc

To better understand the role of slab melt in the petrogenesis of North Cascades magmas, this study focuses on petrogenesis of high-Mg lavas from the two northernmost active volcanoes in Washington. High-Mg andesites (HMA) and basaltic andesites (HMBA) in the Cascade Arc have high Mg# [molar Mg/(Mg+Fe 2+ )] relative to their SiO 2 contents, elevated Nd/Yb, and are Ni- and Cr-enriched. The rock units examined here include the Tarn Plateau HMBA (51.8–54.0 wt% SiO 2 , Mg# 68–70) and Glacier Creek HMA (58.3–58.7 wt% SiO 2 , Mg# 63–64) from the Mount Baker Volcanic Field, and the Lightning Creek HMBA (54.8–54.6 SiO 2 , Mg# 69–73) from Glacier Peak. This study combines major and trace element compositions of minerals and whole rocks to test several petrogenetic hypotheses and to determine which, if any, are applicable to North Cascades HMA and HMBA. In the Tarn Plateau HMBA, rare earth element (REE) equilibrium liquids calculated from clinopyroxene compositions have high Nd/Yb that positively correlates with Mg#. This correlation suggests an origin similar to that proposed for Aleutian adakites, where intermediate, high Nd/Yb slab-derived melts interact with the overlying mantle to become Mg-rich, and subsequently mix with low Nd/Yb, mantle-derived mafic magmas with lower Mg#. In the Glacier Creek HMA, elevated whole-rock MgO and SiO 2 contents resulted from accumulation of xenocrystic olivine and differentiation processes, respectively, but the cause of high Nd/Yb is less clear. However, high whole-rock Sr/P (fluid mobile/fluid immobile) values indicate a mantle source that was fluxed by an enriched, hydrous slab component, likely producing the observed high Nd/Yb REE signature. The Lightning Creek HMBA is a hybridized rock unit with at least three identifiable magmatic components, but only one of which has HMA characteristics. Cr and Mg contents in Cr-spinel and olivine pairs in this HMA component suggest that its source is a strongly depleted mantle, and high whole-rock Sr/P values indicate mantle melting that was induced through hydration, likely adding the component responsible for the observed high Nd/Yb REE pattern. The elevated SiO 2 contents (54.6 wt%) of the HMA component resulted from differentiation or high degrees of partial melting of ultramafic material through the addition of H 2 O. Therefore the Lightning Creek HMBA is interpreted to have originated from a refractory mantle source that underwent melting through interaction with an enriched slab component. Our results indicate that in addition to slab-derived fluids, slab-derived melts also have an important role in the production of HMA and HMBA in the north Cascade Arc.

American Mineralogist↗

USGS2005

The seasonal cycle of diabatic heat storage in the Pacific Ocean

This study quantifies uncertainties in closing the seasonal cycle of diabatic heat storage (DHS) over the Pacific Ocean from 20??S to 60??N through the synthesis of World Ocean Circulation Experiment (WOCE) reanalysis products from 1993 to 1999. These products are DHS from Scripps Institution of Oceanography (SIO); near-surface geostrophic and Ekman currents from Earth and Space Research (ESR); and air-sea heat fluxes from Comprehensive Ocean-Atmosphere Data Set (COADS), National Centers for Environmental Prediction (NCEP), and European Center for Mid-Range Weather Forecasts (ECMWF). With these products, we compute residual heat budget components by differencing long-term monthly means from the long-term annual mean. This allows the seasonal cycle of the DHS tendency to be modeled. Everywhere latent heat flux residuals dominate sensible heat flux residuals, shortwave heat flux residuals dominate longwave heat flux residuals, and residual Ekman heat advection dominates residual geostrophic heat advection, with residual dissipation significant only in the Kuroshio-Oyashio current extension. The root-mean-square (RMS) of the differences between observed and model residual DHS tendencies (averaged over 10??latitude-by-20??longitude boxes) is <20 W m-2 in the interior ocean and <100 W m-2 in the Kuroshio-Oyashio current extension. This reveals that the residual DHS tendency is driven everywhere by some mix of residual latent heat flux, shortwave heat flux, and Ekman heat advection. Suppressing bias errors in residual air-sea turbulent heat fluxes and Ekman heat advection through minimization of the RMS differences reduces the latter to <10 W m-2 over the interior ocean and <25 W m -2 in the Kuroshio-Oyashio current extension. This reveals air-sea temperature and specific humidity differences from in situ surface marine weather observations to be a principal source of bias error, overestimated over most of ocean but underestimated near the Intertropical Convergence Zone. ?? 2005 Elsevier Ltd. All rights reserved.

Progress in Oceanography↗

USGS2023

A 1.2 billion pixel human-labeled dataset for data-driven classification of coastal environments

The world’s coastlines are spatially highly variable, coupled-human-natural systems that comprise a nested hierarchy of component landforms, ecosystems, and human interventions, each interacting over a range of space and time scales. Understanding and predicting coastline dynamics necessitates frequent observation from imaging sensors on remote sensing platforms. Machine Learning models that carry out supervised (i.e., human-guided) pixel-based classification, or image segmentation, have transformative applications in spatio-temporal mapping of dynamic environments, including transient coastal landforms, sediments, habitats, waterbodies, and water flows. However, these models require large and well-documented training and testing datasets consisting of labeled imagery. We describe “Coast Train,” a multi-labeler dataset of orthomosaic and satellite images of coastal environments and corresponding labels. These data include imagery that are diverse in space and time, and contain 1.2 billion labeled pixels, representing over 3.6 million hectares. We use a human-in-the-loop tool especially designed for rapid and reproducible Earth surface image segmentation. Our approach permits image labeling by multiple labelers, in turn enabling quantification of pixel-level agreement over individual and collections of images.

Scientific Data↗

USGS1989

Velocities of antarctic outlet glaciers determined from sequential Landsat images

Approximately 91.0 percent of the volume of present-day glacier ice on Earth is in Antarctica; Greenland contains about another 8.3 percent of the volume. Thus, together, these two great ice sheets account for an estimated 99.3 percent of the total. Long-term changes in the volume of glacier ice on our planet are the result of global climate change. Because of the relationship of global ice volume to sea level (± 330 cubic kilometers of glacier ice equals ± 1 millimeter sea level), changes in the mass balance of the antarctic ice sheet are of particular importance. Whether the mass balance of the east and west antarctic ice sheets is positive or negative is not known. Estimates of mass input by total annual precipitation for the continent have been made from scattered meteorological observations (Swithinbank 1985). The magnitude of annual ablation of the ice sheet from calving of outlet glaciers and ice shelves is also not well known. Although the velocities of outlet glaciers can be determined from field measurements during the austral summer,the technique is costly, does not cover a complete annual cycle,and has been applied to just a few glaciers. To increase the number of outlet glaciers in Antarctica for which velocities have been determined and to provide additional data for under-standing the dynamics of the antarctic ice sheets and their response to global climate change, sequential Landsat image of several outlet glaciers were measured.

Antarctic Journal of the United States↗

USGS2014

Seismicity, faulting, and structure of the Koyna-Warna seismic region, Western India from local earthquake tomography and hypocenter locations

Although seismicity near Koyna Reservoir (India) has persisted for ~50 years and includes the largest induced earthquake ( M 6.3) reported worldwide, the seismotectonic framework of the area is not well understood. We recorded ~1800 earthquakes from 6 January 2010 to 28 May 2010 and located a subset of 343 of the highest-quality earthquakes using the tomoDD code of Zhang and Thurber (2003) to better understand the framework. We also inverted first arrivals for 3-D Vp , Vs , and Vp / Vs and Poisson's ratio tomography models of the upper 12 km of the crust. Epicenters for the recorded earthquakes are located south of the Koyna River, including a high-density cluster that coincides with a shallow depth (<1.5 km) zone of relatively high Vp and low Vs (also high Vp / Vs and Poisson's ratios) near Warna Reservoir. This anomalous zone, which extends near vertically to at least 8 km depth and laterally northward at least 15 km, is likely a water-saturated zone of faults under high pore pressures. Because many of the earthquakes occur on the periphery of the fault zone, rather than near its center, the observed seismicity-velocity correlations are consistent with the concept that many of the earthquakes nucleate in fractures adjacent to the main fault zone due to high pore pressure. We interpret our velocity images as showing a series of northwest trending faults locally near the central part of Warna Reservoir and a major northward trending fault zone north of Warna Reservoir.

Maharashtra↗

USGS2005

Magmatic effects of the Cobb hot spot on the Juan de Fuca Ridge

The interaction of the Juan de Fuca Ridge with the Cobb hot spot has had a considerable influence on the magmatism of the Axial Segment of the ridge, the second-order segment that overlies the hot spot. In addition to the construction of the large volcanic edifice of Axial Seamount, the Axial Segment has shallow bathymetry and a prevalence of constructional volcanic features along its 100-km length, suggesting that hot spot-derived magmas supplement and oversupply the ridge. Lavas are generally more primitive at Axial Seamount and more evolved in the Axial Segment rift zones, suggesting that fractional crystallization is enhanced with increasing distance from the hot spot because of a reduced magma supply and more rapid cooling. Although the Cobb hot spot is not an isotopically enriched plume, it produces lavas with some distinct geochemical characteristics relative to normal mid-ocean ridge basalt, such as enrichments in alkalis and highly incompatible trace elements, that can be used as tracers to identify the presence and prevalence of the hot spot influence along the ridge. These characteristics are most prominent at Axial Seamount and decline in gradients along the Axial Segment. The physical model that can best explain the geochemical observations is a scenario in which hot spot and mid-ocean ridge basalt (MORB) magmas mix to varying degrees, with the proportions controlled by the depth to the MORB source. Modeling of two-component mixing suggests that MORB is the dominant component in most Axial Segment basalts. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Transient rheology of the upper mantle beneath central Alaska inferred from the crustal velocity field following the 2002 Denali earthquake

The M7.9 2002 Denali earthquake, Alaska, is one of the largest strike-slip earthquakes ever recorded. The postseismic GPS velocity field around the 300-km-long rupture is characterized by very rapid horizontal velocity up to ???300 mm/yr for the first 0.1 years and slower but still elevated horizontal velocity up to ???100 mm/yr for the succeeding 1.5 years. I find that the spatial and temporal pattern of the displacement field may be explained by a transient mantle rheology. Representing the regional upper mantle as a Burghers body, I infer steady state and transient viscosities of ??1 = 2.8 ?? 1018 Pa s and ??2 = 1.0 ?? 1017 Pa s, respectively, corresponding to material relaxation times of 1.3 and 0.05 years. The lower crustal viscosity is poorly constrained by the considered horizontal velocity field, and the quoted mantle viscosities assume a steady state lower crust viscosity that is 7??1. Systematic bias in predicted versus observed velocity vectors with respect to a fixed North America during the first 3-6 months following the earthquake is reduced when all velocity vectors are referred to a fixed site. This suggests that the post-Denali GPS time series for the first 1.63 years are shaped by a combination of a common mode noise source during the first 3-6 months plus viscoelastic relaxation controlled by a transient mantle rheology.

Journal of Geophysical Research B: Solid Earth↗

USGS1975

Compositional evidence regarding the influx of interplanetary materials onto the lunar surface

Siderophilic element/Ir ratios are higher in mature lunar soils from highlands sites than in those from mare sites. We infer that the population of materials responsible for the early intense bombardment of the Moon had high ratios, and that the population responsible for the essentially constant flux has low ratios. No group of chondrites has siderophile/Ir ratios identical to those in the mare or highlands soils; CM chondrites are the most similar, and CM-like materials may account for a major fraction of Earth-crossing materials during the past 3.7 b.y. Siderophile/Ir ratios may be used to determine the amount of highlands regolith in soils or breccias from the mare-highlands interface areas (Apollo 15 and 17), and to infer the time of formation of highlands breccias whose sideropbiles originated in mature soils. Arguments are summarized against the viewpoint that the siderophiles in most highlands breccias originated in basin-forming projectiles. Differences in mature soil siderophile concentrations at Apollo 14 and 16 indicate a substantially greater concentration at the latter site immediately following the Imbrium event. Siderophile concentrations are used to estimate mean regolith depths at the landing sites; as relative values these are more precise than estimates based on seismic or crater observations. The longlived flux is calculated to be 2.9 g cm −2 b.y. −1 averaged over the past 3.7 b.y. A consideration of the relationship between mass fluence and time indicates that the mass flux decreased with a half-life of about 40 m.y. immediately following the Imbrium event.

The Moon↗

USGS2011

Global climate change implications for coastal and offshore oil and gas development

The discussion and debate about climate change and oil and gas resource development has generally focused on how fossil fuel use affects the Earth's climate. This paper explores how the changing climate is likely to affect oil and gas operations in low-lying coastal areas and the outer continental shelf. Oil and gas production in these regions comprises a large sector of the economies of many energy producing nations. Six key climate change drivers in coastal and marine regions are characterized with respect to oil and gas development: changes in carbon dioxide levels and ocean acidity, air and water temperature, precipitation patterns, the rate of sea level rise, storm intensity, and wave regime. These key drivers have the potential to independently and cumulatively affect coastal and offshore oil and gas exploration, production, and transportation, and several impacts of climate change have already been observed in North America.

Energy Policy↗

USGS2016

Geologic and geochemical results from boreholes drilled in Yellowstone National Park, Wyoming, 2007 and 2008

Between 2007 and 2008, seven Earthscope Plate Boundary Observatory (PBO) boreholes ranging in depth from about 200 to 800 feet deep were drilled in and adjacent to the Yellowstone caldera in Yellowstone National Park, for the purpose of installing volcano monitoring instrumentation. Five of the seven boreholes were equipped with strainmeters, downhole seismometers, and tiltmeters. Data collected during drilling included field observations of drill cuttings, stratigraphy within the boreholes, water temperature, and water and drill cuttings samples from selected depths. Six of the seven boreholes encountered rhyolite lavas and tuffs. The rhyolite lavas compose the Canyon flow, the Gardner River flow, the Gibbon River flow, the Hayden Valley flow, the Nez Perce Creek flow, and the West Thumb flow. Boreholes also penetrated a vertical sequence through the Lava Creek Tuff and the Tuff of Bluff Point. In addition, one borehole drilled through a Swan Lake Flat Basalt sequence and terminated in a rhyolite lava flow. After drilling the seven PBO boreholes, cuttings were examined and selected for preparation of grain mounts, thin sections, and geochemical analysis. Major ions and trace elements (including rare earth elements) of selected cuttings were determined by x-ray fluorescence (XRF) and inductively coupled plasma-mass spectrometry (ICP-MS); the ICP-MS provided more precise trace-element analysis than XRF. A preliminary interpretation of the results of geochemical analyses generally shows a correlation between borehole cuttings and previously mapped geology. The geochemical data and borehole stratigraphy presented in this report provide a foundation for future petrologic, geochemical, and geophysical studies.

Wyoming↗

USGS2008

In situ Raman spectroscopic investigation of the structure of subduction-zone fluids

In situ Raman spectra of synthetic subduction-zone fluids (KAlSi3O8-H2O system) were measured to 900?? and 2.3 GPa using a hydrothermal diamond-anvil cell. The structures of aqueous fluid and hydrous melt become closer when conditions approach the second critical endpoint. Almost no three-dimensional network was observed in the supercritical fluid above 2 GPa although a large amount of silicate component is dissolved, suggesting that the physical and chemical properties of these phases change drastically at around the second critical endpoint. Our experimental results indicate that the fluids released from a subducting slab change from aqueous fluid to supercritical fluid with increasing depth under the volcanic arcs. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Tracking magma volume recovery at Okmok Volcano using GPS and an unscented kalman filter

Changes beneath a volcano can be observed through position changes in a GPS network, but distinguishing the source of site motion is not always straightforward. The records of continuous GPS sites provide a favorable data set for tracking magma migration. Dense campaign observations usually provide a better spatial picture of the overall deformation field, at the expense of an episodic temporal record. Combining these observations provides the best of both worlds. A Kalman filter provides a means for integrating discrete and continuous measurements and for interpreting subtle signals. The unscented Kalman filter (UKF) is a nonlinear method for time‐dependent observations. We demonstrate the application of this technique to deformation data by applying it to GPS data collected at Okmok volcano. Seven years of GPS observations at Okmok are analyzed using a Mogi source model and the UKF. The deformation source at Okmok is relatively stable at 2.5 km depth below sea level, located beneath the center of the caldera, which means the surface deformation is caused by changes in the strength of the source. During the 7 years of GPS observations more than 0.5 m of uplift has occurred, a majority of that during the time period January 2003 to July 2004. The total volume recovery at Okmok since the last eruption in 1997 is ∼60–80%. The UKF allows us to solve simultaneously for the time‐dependence of the source strength and for the location without a priori information about the source. 2009 by the American Geophysical Union.

Alaska↗

USGS2015

Indium: bringing liquid-crystal displays into focus

Introduction Indium is rare in the Earth’s crust. The continental crust contains an average of about 50 parts per billion of indium, whereas the oceanic crust contains about 72 parts per billion, which is similar to meteoritic abundances and comparable to the crustal abundance of silver. Indium minerals are rare in nature and only 12 indium minerals are known. In its elemental form, indium is a soft, lustrous, silver-white metal with a low melting point relative to other metals. It is ductile and malleable, even at temperatures approaching absolute zero, making it ideal for cryogenic applications. Indium was discovered in the mid-1800s by two German chemists who were investigating zinc ores from Freiberg, Saxony. They named it after the distinctive indigo-blue color observed in its emission spectrum. For years indium remained only a scientific curiosity and early applications of indium were few, but included manufacturing of light-emitting diodes and coatings for bearings used in aircraft engines. Indium-bearing nuclear control rods became more widely used in the 1970s, and today the major application of indium is in manufacturing liquid-crystal displays. Compared to more abundant industrial metals such as lead and zinc, information about the behavior and toxicity of indium in the environment is limited. However, many indium compounds have been proven to be toxic to animals.

Fact Sheet↗

USGS2022

Thermophysical and compositional properties of paleobedforms on Mars

Bedforms on Earth and Mars are often preserved in the rock record in the form of sedimentary rock with distinct cross-bedding. On rare occasions, the full-surface geometry of a bedform can be preserved through burial and lithification. These features, known as paleobedforms, are found in a variety of geographic locations on Mars. Evidence in the morphology of paleobedforms, such as the retention of impact craters and steep erosional scarps, suggests that these features are well-lithified and capable of withstanding prolonged weathering and erosion. Here, we present results from thermophysical and compositional analyses on a subset of the best preserved paleobedform candidate fields on Mars. Thermophysical modeling elucidates the changes these bedforms underwent from their unconsolidated, particulate nature to their currently observed properties. Certain paleobedforms have elevated thermal inertias (e.g., ∼300–500 J·m −2 ·s −1/2 ·K −1 ) when compared with modern bedforms (∼250 J·m −2 ·s −1/2 ·K −1 ), and modeling indicates that they have cement volumes of 0.8%–1.5% even as high as 30%. However, most paleobedform candidates have unexpectedly low thermal inertia when compared with modern dunes. Additionally, compositional analyses reveal a range of spectral characteristics within paleobedforms (e.g., primary and secondary alteration products). These features add to the already existing class of Martian surfaces in which thermal inertia does not seem to correspond to erodibility, cohesion, or mechanical strength. Studying paleobedforms with both raised and nonraised thermal inertia has provided new insights into lithification on Mars and constrained the environmental conditions leading to the formation of these enigmatic features.

Journal of Geophysical Research - Planets↗

USGS2023

Forecasts of polar bear (Ursus maritimus) land use in the southern Beaufort and Chukchi Seas, 2040–65

This report provides analysis to extend the 2040 forecasts of polar bear ( Ursus maritimus ) land use for the southern Beaufort and Chukchi Sea populations presented in a recent publication (Rode and others, 2022) through the year 2065. To inform long-term polar bear management considerations, we provide point-estimate forecasts and 95-percent prediction intervals of the proportion of polar bear populations summering onshore for 21 days or more (≥) and their duration onshore every 5 years from 2040 to 2065. Because sea-ice projections based on earth system models show greater divergence with emission scenarios after 2040, we have provided forecasts for three greenhouse gas emission scenarios, SSP1-2.6, SSP2-4.5 and SSP5-8.5, compared to the two emission scenarios used in Rode and others (2022). Our forecasting methods estimated that 61–97 percent of polar bears in the southern Beaufort Sea and 80–100 percent of polar bears in the Chukchi Sea populations may summer onshore for ≥21 days by 2065. Forecasts of mean duration onshore were 105–158 days for polar bears in the southern Beaufort Sea and 111–178 days in the Chukchi Sea populations by 2065. Sea ice conditions projected to occur by 2065 could alter the current patterns of bear behavior from what has been observed over the past 30 years. As a result, these extended projections are associated with a higher degree of uncertainty than estimates through 2040, especially under the SSP5-8.5 scenario.

Open-File Report↗

USGS2010

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

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

Hawai'i↗