Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “The Earth Observer”

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 1,117 records · Page 62Linked to original sources
USGS2005

Improvement in absolute calibration accuracy of Landsat-5 TM with Landsat-7 ETM+ data

The ability to detect and quantify changes in the Earth's environment depends on satellites sensors that can provide calibrated, consistent measurements of Earth's surface features through time. A critical step in this process is to put image data from subsequent generations of sensors onto a common radiometric scale. To evaluate Landsat-5 (L5) Thematic Mapper's (TM) utility in this role, image pairs from the L5 TM and Landsat-7 (L7) Enhanced Thematic Mapper Plus (ETM+) sensors were compared. This approach involves comparison of surface observations based on image statistics from large common areas observed eight days apart by the two sensors. The results indicate a significant improvement in the consistency of L5 TM data with respect to L7 ETM+ data, achieved using a revised Look-Up-Table (LUT) procedure as opposed to the historical Internal Calibrator (IC) procedure previously used in the L5 TM product generation system. The average percent difference in reflectance estimates obtained from the L5 TM agree with those from the L7 ETM+ in the Visible and Near Infrared (VNIR) bands to within four percent and in the Short Wave Infrared (SWIR) bands to within six percent.

Conference Paper↗
USGS1982

A Proposed Model for the International Geomagnetic Reference Field-1965

A best current model of the main geomagnetic field is presented as a response to a need for an “International Geomagnetic Reference Field”. This model is described by a series of 120 spherical harmonic coefficients and their first and second time derivatives from an epoch 1960.0. It was derived from a sample of all magnetic survey data available from the interval 1900-1964 plus a recent global distribution of preliminary total field observations from the OGO-2 (1965-81A) spacecraft for epoch 1965.8. A duplicate data selection was made and the resulting field model compared with the first to help evaluate the minimum error. It was noted that the root —mean—square difference between the two models was about 30γ in the force components, 0.04 degrees in dip and 0.3 degrees in declination at the earth's surface for 1965.0.

Journal of Geomagnetism & Geoelectricity↗
USGS2015

Near-surface versus fault zone damage following the 1999 Chi-Chi earthquake: Observation and simulation of repeating earthquakes

We observe crustal damage and its subsequent recovery caused by the 1999  M 7.6 Chi-Chi earthquake in central Taiwan. Analysis of repeating earthquakes in Hualien region, ~70 km east of the Chi-Chi earthquake, shows a remarkable change in wave propagation beginning in the year 2000, revealing damage within the fault zone and distributed across the near surface. We use moving window cross correlation to identify a dramatic decrease in the waveform similarity and delays in the S wave coda. The maximum delay is up to 59 ms, corresponding to a 7.6% velocity decrease averaged over the wave propagation path. The waveform changes on either side of the fault are distinct. They occur in different parts of the waveforms, affect different frequencies, and the size of the velocity reductions is different. Using a finite difference method, we simulate the effect of postseismic changes in the wavefield by introducing S wave velocity anomaly in the fault zone and near the surface. The models that best fit the observations point to pervasive damage in the near surface and deep, along-fault damage at the time of the Chi-Chi earthquake. The footwall stations show the combined effect of near-surface and the fault zone damage, where the velocity reduction (2–7%) is twofold to threefold greater than the fault zone damage observed in the hanging wall stations. The physical models obtained here allow us to monitor the temporal evolution and recovering process of the Chi-Chi fault zone damage.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Refining fault slip rates using multiple displaced terrace risers-An example from the Honey Lake fault, NE California, USA

Faulted terrace risers are semi-planar features commonly used to constrain Quaternary slip rates along strike-slip faults . These landforms are difficult to date directly and therefore their ages are commonly bracketed by age estimates of the adjacent upper and lower terrace surfaces. However, substantial differences in the ages of the upper and lower terrace surfaces (a factor of 2.4 difference observed globally) produce large uncertainties in the slip-rate estimate. In this investigation, we explore how the full range of displacements and bounding ages from multiple faulted terrace risers can be combined to yield a more accurate fault slip rate. We use 0.25-m cell size digital terrain models derived from airborne lidar data to analyze three sites where terrace risers are offset right-laterally by the Honey Lake fault in NE California, USA. We use ages for locally extensive subhorizontal surfaces to bracket the time of riser formation: an upper surface is the bed of abandoned Lake Lahontan having an age of 15.8 ± 0.6 ka and a lower surface is a fluvial terrace abandoned at 4.7 ± 0.1 ka. We estimate lateral offsets of the risers ranging between 6.6 and 28.3 m (median values), a greater than fourfold difference in values. The amount of offset corresponds to the riser's position relative to modern stream meanders: the smallest offset is in a meander cutbank position, whereas the larger offsets are in straight channel or meander point-bar positions. Taken in isolation, the individual terrace-riser offsets yield slip rates ranging from 0.3 to 7.1 mm/a. However, when the offset values are collectively assessed in a probabilistic framework, we find that a uniform (linear) slip rate of 1.6 mm/a (1.4–1.9 mm/a at 95% confidence) can satisfy the data, within their respective uncertainties. This investigation demonstrates that integrating observations of multiple offset elements (crest, midpoint, and base) from numerous faulted and dated terrace risers at closely spaced sites can refine slip-rate estimates on strike-slip faults.

California↗
USGS2000

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

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

Arthur's Pass↗
USGS2023

A large sediment accretion wave along a northern California littoral cell

The northern California littoral cell of the Klamath River, which is a mixed rocky and sandy system with significant shoreline curvature, was investigated by examining ∼40 yr of satellite-derived shoreline positions and historical records. We find that an accretion wave of sediment was initiated near the Klamath River mouth in the late 1980s and translated downcoast over the subsequent decades. The wave passed rapidly (∼2,500 m/yr) through a rocky coastal reach with more oblique wave directions and slowly through a sandy reach (∼200 m/yr) where wave crests approach at more normal angles. Within the sandy reach, the accretion wave extended over 200 m offshore, was ∼10 km long, incorporated 20 ± 6 million m 3 of sediment, and averaged 1.3 ± 0.4 million m 3 /yr of longshore sediment transport over a 20-yr interval. Diffusion of the accretion wave was observed, but the diffusivity coefficient ( ε obs ∼0.01 m 2 /s) was lower than values predicted by theory, which we attribute to net sediment transport convergence in the study area caused by the curvature of the shoreline. Examining historical records, we find that increased sediment discharge in the Klamath River occurred during the 20th century from industrial-scale logging and climatic extremes. Thus, we hypothesize that increased river sediment discharge introduced new sediment to the littoral cell that initiated the observed accretion wave. These hypotheses can be tested with stratigraphic and mineralogic investigations of the broad study area beach that has formed during the past 150 years.

California↗
USGS2020

Holocene paleofloods and their climatological context, Upper Colorado River Basin, USA

Given its singular importance for water resources in the southwestern U.S., the Upper Colorado River Basin (UCRB) is remarkable for the paucity of its conventional hydrological record of extreme flooding. This study uses paleoflood hydrology to examine a small portion the underutilized, but very extensive natural record of Holocene extreme floods in the UCRB. We perform a meta-analysis of 77 extreme paleofloods from seven slackwater deposit sites in the UCRB to show linkages between Holocene climate patterns and extreme floods. The analysis demonstrates several clusters of extreme flood activity: 8040-7790, 3600-3460, 2880-2740, 2330-700, and 620-0 years BP. The extreme paleofloods were found to occur during both dry and wet periods in the paleoclimate record. When compared with independent paleoclimatic records across the Rocky Mountains and the southwestern U.S., the observed temporal clustering pattern of UCRB extreme paleofloods shows associations with periods of abruptly intensified North Pacific-derived storms connected with enhanced El Niño variability.

Upper Colorado River Basin↗
USGS2022

Using a grid-search approach to validate the Graves-Pitarka broadband simulation method

This work assesses the ability of the Graves–Pitarka simulation approach to reproduce observed ground motions for 12 California and Baja California earthquakes. A total of 240 realizations are computed for each earthquake and compared with recorded strong motions from near-fault sites. In addition to spatial variability in slip, each realization samples from discrete combinations of average rupture speed, rise time, and down-dip fault width. Ground motions for each realization of a given earthquake are compared to the observations using a pseudo-spectral acceleration goodness-of-fit (GoF) metric. Parameters for the lowest misfit cases are then tabulated to develop relations for estimating median values and ranges for future applications. The results are generally consistent with published scaling relations for estimating rupture dimensions as a function of magnitude. Additionally, I find the average rise time scales as 2.3 × α T × 10 − 9 × M o 1 / 3 "> 2.3 × α T × 10 - 9 × M o 1 / 3 ( M o is seismic moment in dyne-cm) and average rupture speed as ( 0.765 ± 0.075 ) × α T − 1 × V s "> ( 0.765 ± 0.075 ) × α T - 1 × V s ( V s is local shear wave velocity), where 2.3 × α T × 10 − 9 × M o 1 / 3 "> α T is a mechanism adjustment ranging from a value of 1.0 for pure strike-slip to 0.9 for pure reverse-slip cases. There are cases where slightly different combinations of parameters produce equally good fits to the observations, demonstrating the non-uniqueness of using a single GoF metric in this approach. Nonetheless, the results reinforce the importance of adequately sampling ranges of rupture parameters when performing validations, as well as when simulating ground motions for future events.

Earth, Planets and Space↗
USGS2004

Mineral dissolution in the Cape Cod aquifer, Massachusetts, USA: I . Reaction stoichiometry and impact of accessory feldspar and glauconite on strontium isotopes, solute concentrations, and REY distribution

To compare relative reaction rates of mineral dissolution in a mineralogically simple groundwater aquifer, we studied the controls on solute concentrations, Sr isotopes, and rare earth element and yttrium (REY) systematics in the Cape Cod aquifer. This aquifer comprises mostly carbonate-free Pleistocene sediments that are about 90% quartz with minor K-feldspar, plagioclase, glauconite, and Fe-oxides. Silica concentrations and pH in the groundwater increase systematically with increasing depth, while Sr isotopic ratios decrease. No clear relationship between 87 Sr/ 86 Sr and Sr concentration is observed. At all depths, the 87 Sr/ 86 Sr ratio of the groundwater is considerably lower than the Sr isotopic ratio of the bulk sediment or its K-feldspar component, but similar to that of a plagioclase-rich accessory separate obtained from the sediment. The Si- 87 Sr/ 86 Sr-depth relationships are consistent with dissolution of accessory plagioclase. In addition, solutes such as Sr, Ca, and particularly K show concentration spikes superimposed on their respective general trends. The K-Sr- 87 Sr/ 86 Sr systematics suggests that accessory glauconite is another major solute source to Cape Cod groundwater. Although the authigenic glauconite in the Cape Cod sediment is rich in Rb, it is low in in-grown radiogenic 87 Sr because of its young Pleistocene age. The low 87 Sr/ 86 Sr ratios are consistent with equilibration of glauconite with seawater. The impact of glauconite is inferred to vary due to its variable abundance in the sediments. In the Cape Cod groundwater, the variation of REY concentrations with sampling depth resembles that of K and Rb, but differs from that of Ca and Sr. Shale-normalized REY patterns are light REY depleted, show negative Ce anomalies and super-chondritic Y/Ho ratios, but no Eu anomalies. REY input from feldspar, therefore, is insignificant compared to input from a K-Rb-bearing phase, inferred to be glauconite. These results emphasize that interpretation of groundwater chemistry, even in relatively simple aquifers, may be complicated by solute contributions from “exotic” accessory minerals such as glauconite. To detect such peculiarities, groundwater studies should combine the study of elemental concentration and isotopic composition of several solutes that show different geochemical behavior.

Massachusetts↗
USGS2011

Volcanic plume height measured by seismic waves based on a mechanical model

In August 2008 an unmonitored, largely unstudied Aleutian volcano, Kasatochi, erupted catastrophically. Here we use seismic data to infer the height of large eruptive columns such as those of Kasatochi based on a combination of existing fluid and solid mechanical models. In so doing, we propose a connection between a common, observable, short-period seismic wave amplitude to the physics of an eruptive column. To construct a combined model, we estimate the mass ejection rate of material from the vent on the basis of the plume height, assuming that the height is controlled by thermal buoyancy for a continuous plume. Using the estimated mass ejection rate, we then derive the equivalent vertical force on the Earth through a momentum balance. Finally, we calculate the far-field surface waves resulting from the vertical force. The model performs well for recent eruptions of Kasatochi and Augustine volcanoes if v , the velocity of material exiting the vent, is 120-230 m s -1 . The consistency between the seismically inferred and measured plume heights indicates that in these cases the far-field ~1 s seismic energy radiated by fluctuating flow in the volcanic jet during the eruption is a useful indicator of overall mass ejection rates. Thus, use of the model holds promise for characterizing eruptions and evaluating ash hazards to aircraft in real time on the basis of far-field short-period seismic data. This study emphasizes the need for better measurements of eruptive plume heights and a more detailed understanding of the full spectrum of seismic energy radiated coeruptively.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

Hydrogeology and water quality of the Floridan aquifer system and effect of Lower Floridan aquifer withdrawals on the Upper Floridan aquifer at Barbour Pointe Community, Chatham County, Georgia, 2013

Two test wells were completed at the Barbour Pointe community in western Chatham County, near Savannah, Georgia, in 2013 to investigate the potential of using the Lower Floridan aquifer as a source of municipal water supply. One well was completed in the Lower Floridan aquifer at a depth of 1,080 feet (ft) below land surface; the other well was completed in the Upper Floridan aquifer at a depth of 440 ft below land surface. At the Barbour Pointe test site, the U.S. Geological Survey completed electromagnetic (EM) flowmeter surveys, collected and analyzed water samples from discrete depths, and completed a 72-hour aquifer test of the Floridan aquifer system withdrawing from the Lower Floridan aquifer. Based on drill cuttings, geophysical logs, and borehole EM flowmeter surveys collected at the Barbour Pointe test site, the Upper Floridan aquifer extends 369 to 567 ft below land surface, the middle semiconfining unit, separating the two aquifers, extends 567 to 714 ft below land surface, and the Lower Floridan aquifer extends 714 to 1,056 ft below land surface. A borehole EM flowmeter survey indicates that the Upper Floridan and Lower Floridan aquifers each contain four water-bearing zones. The EM flowmeter logs of the test hole open to the entire Floridan aquifer system indicated that the Upper Floridan aquifer contributed 91 percent of the total flow rate of 1,000 gallons per minute; the Lower Floridan aquifer contributed about 8 percent. Based on the transmissivity of the middle semiconfining unit and the Floridan aquifer system, the middle semiconfining unit probably contributed on the order of 1 percent of the total flow. Hydraulic properties of the Upper Floridan and Lower Floridan aquifers were estimated based on results of the EM flowmeter survey and a 72-hour aquifer test completed in Lower Floridan aquifer well 36Q398. The EM flowmeter data were analyzed using an AnalyzeHOLE-generated model to simulate upward borehole flow and determine the transmissivity of water-bearing zones. Aquifer-test data were analyzed with a two-dimensional, axisymmetric, radial, transient, groundwater-flow model using MODFLOW–2005. The flowmeter-survey and aquifer-test simulations provided an estimated transmissivity of about 60,000 square feet per day for the Upper Floridan aquifer and about 5,000 square feet per day for the Lower Floridan aquifer. Water in discrete-depth samples collected from the Upper Floridan aquifer, middle semiconfining unit, and Lower Floridan aquifer during the EM flowmeter survey in August 2013 was low in dissolved solids. Tested constituents were in concentrations within established U.S. Environmental Protection Agency drinking water-quality criteria. Concentrations of measured constituents in water samples from Lower Floridan aquifer well 36Q398 collected at the end of the 72-hour aquifer test in November 2013 were generally higher than in the discrete-depth samples collected during EM flowmeter testing in August 2013 but remained within established drinking water-quality criteria. Water-level data for the aquifer test were filtered for external influences such as barometric pressure, earth-tide effects, and long-term trends to enable detection of small (less than 1 ft) water-level responses to aquifer-test withdrawal. During the 72-hour aquifer test, the Lower Floridan aquifer was pumped at a rate of 750 gallons per minute resulting in a drawdown response of 35.5 ft in the pumped well; 1.6 ft in the Lower Floridan aquifer observation well located about 6,000 ft west of the pumped well; and responses of 0.7, 0.6, and 0.4 ft in the Upper Floridan aquifer observation wells located about 36 ft, 6,000 ft, and 6,800 ft from the pumped well, respectively

Georgia↗
USGS2026

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter↗
USGS2012

Morphological and chemical evidence of stromatolitic deposits in the 2.75 Ga Carajás banded iron formation, Brazil

We describe evidence of biogenicity in the morphology and carbon content of well-preserved, Neoarchean samples of banded iron formation (BIF) from Carajás, Brazil. Silica-rich BIF layers contain translucent ellipsoidal or trapezoidal structures (∼5–10 μm diameter) composed of silica, hematite, and kerogen, which are arranged in larger ring-like forms (rosettes). Stable carbon isotope analysis yields a δ 13 C value of −24.5‰ indicating that the contained carbon is likely biogenic. Raman and SEM analyses, as well as wavelength-dispersive X-ray elemental maps, show kerogen inside the rosette forms. Within the iron-rich BIF layers, tubular structures (0.5–5 μm) were observed between hematite granules and blades. Kerogen and kaolinite are present in these structures. Both the rosettes and the tubular structures resemble morphologies that are characteristic of some bacterial species. We hypothesize that the Carajás BIFs originated as biomats formed by one or more species that over time produced large stromatolitic structures. The rosettes and the tubular structures, associated with chert-rich and iron-rich BIF layers, respectively, may represent two different species, or perhaps, two phases of a bacterium life cycle. For example, some modern myxobacteria exhibit similar morphologies in their resting and vegetative stages. Fe(III) precipitation may have occurred by contact of Fe(II) with bacterial slime, leading to oxidation by chemical reactions with exposed polysaccharide hydroxyl and carboxyl groups. The Fe(III) would then have been available for use as a source of energy in a dissimilatory iron reduction type of metabolism. Organic carbon input presumably came from primary producers (not necessarily aerobic) within the local water column, perhaps in shallow-water communities. Alternatively, the carbon may have originated by Fischer–Tropsch synthesis at ocean hydrothermal vents. The observed lateral continuity of BIF layers may perhaps be explained by chemical signaling by the bacteria of favorable or unfavorable environmental conditions, leading to nearly synchronous cell morphogenesis from a vegetative to resting phase and vice versa.

Earth and Planetary Science Letters↗
USGS2007

Cosmogenic 10Be and 36Cl geochronology of offset alluvial fans along the northern Death Valley fault zone: Implications for transient strain in the eastern California shear zone

The northern Death Valley fault zone (NDVFZ) has long been recognized as a major right-lateral strike-slip fault in the eastern California shear zone (ECSZ). However, its geologic slip rate has been difficult to determine. Using high-resolution digital topographic imagery and terrestrial cosmogenic nuclide dating, we present the first geochronologically determined slip rate for the NDVFZ. Our study focuses on the Red Wall Canyon alluvial fan, which exposes clean dextral offsets of seven channels. Analysis of airborne laser swath mapping data indicates ???297 ?? 9 m of right-lateral displacement on the fault system since the late Pleistocene. In situ terrestrial cosmogenic 10Be and 36C1 geochronology was used to date the Red Wall Canyon fan and a second, correlative fan also cut by the fault. Beryllium 10 dates from large cobbles and boulders provide a maximum age of 70 +22/-20 ka for the offset landforms. The minimum age of the alluvial fan deposits based on 36Cl depth profiles is 63 ?? 8 ka. Combining the offset measurement with the cosmogenic 10Be date yields a geologic fault slip rate of 4.2 +1.9/-1.1 mm yr-1, whereas the 36Cl data indicate 4.7 +0.9/-0.6 mm yr-1 of slip. Summing these slip rates with known rates on the Owens Valley, Hunter Mountain, and Stateline faults at similar latitudes suggests a total geologic slip rate across the northern ECSZ of ???8.5 to 10 mm yr-1. This rate is commensurate with the overall geodetic rate and implies that the apparent discrepancy between geologic and geodetic data observed in the Mojave section of the ECSZ does not extend north of the Garlock fault. Although the overall geodetic rates are similar, the best estimates based on geology predict higher strain rates in the eastern part of the ECSZ than to the west, whereas the observed geodetic strain is relatively constant. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Three-dimensional geophysical mapping of rock alteration and water content at Mount Adams, Washington: Implications for lahar hazards

[1] Hydrothermally altered rocks, particularly if water saturated, can weaken stratovolcanoes, thereby increasing the potential for catastrophic sector collapses that can lead to far-traveled, destructive debris flows. Evaluating the hazards associated with such alteration is difficult because alteration has been mapped on few active volcanoes and the distribution and intensity of subsurface alteration are largely unknown on any active volcano. At Mount Adams, some Holocene debris flows contain abundant hydrothermal minerals derived from collapse of the altered edifice. Intense hydrothermal alteration significantly reduces the resistivity and magnetization of volcanic rock, and therefore hydrothermally altered rocks can be identified with helicopter electromagnetic and magnetic measurements. Electromagnetic and magnetic data, combined with geological mapping and rock property measurements, indicate the presence of appreciable thicknesses of hydrothermally altered rock in the central core of Mount Adams north of the summit. We identify steep cliffs at the western edge of this zone as the likely source for future large debris flows. In addition, the electromagnetic data identified water in the brecciated core of the upper 100–200 m of the volcano. Water helps alter the rocks, reduces the effective stress, thereby increasing the potential for slope failure, and acts, with entrained melting ice, as a lubricant to transform debris avalanches into lahars. Therefore knowing the distribution of water is also important for hazard assessments. Our results demonstrate that high-resolution geophysical and geological observations can yield unprecedented views of the three-dimensional distribution of altered rock and shallow pore water aiding evaluation of the debris avalanche hazard.

Journal of Geophysical Research B: Solid Earth↗
USGS2021

A unified model of crustal stress heterogeneity from borehole breakouts and earthquake focal mechanisms

Observations of crustal stress orientation from the regional inversion of earthquake focal mechanisms often conflict with those from borehole breakouts, possibly indicating local stress heterogeneity, either laterally or with depth. To investigate this heterogeneity, we compiled SHmax estimates from previous studies for 57 near‐vertical boreholes with measured breakout azimuths across the Los Angeles region. We identified subsets of earthquake focal mechanisms from established earthquake catalogs centered around each borehole with various criteria for maximum depth and maximum lateral distance from the borehole. Each subset was independently inverted for 3‐D stress orientation and corresponding SHmax probability distributions, then compared with the corresponding borehole breakout‐derived estimate. We find good agreement when both methods sample the basement stress (breakouts are close to the sediment‐basement interface), or when both methods sample the mid‐basin stress (sufficient earthquakes are present within a sedimentary basin). Along sedimentary basin margins, in contrast, we find acceptable agreement only when focal mechanisms are limited to shallow and close earthquakes, implying short‐length‐scale heterogeneity of <20 km. While the region as a whole shows evidence of both lateral and vertical stress orientation heterogeneity, we find a more homogeneous stress state within basement rock, over length scales of 1–35 km. These results reconcile the apparently conflicting observations of short‐length‐scale heterogeneity observed in boreholes, which sample primarily the basins, with the relative homogeneity of stress inferred from focal mechanisms, which sample primarily the basement, and imply distinct regimes for the appropriate use of each type of stress indicator.

JGR Solid Earth↗
USGS2022

Microgravity change during the 2008-2018 Kı̄lauea summit eruption: Nearly a decade of subsurface mass accumulation

Results from nine microgravity campaigns from Kı̄lauea, Hawaiʻi, spanning most of the volcano's 2008–2018 summit eruption, indicate persistent mass accumulation at shallow levels. A weighted least squares approach is used to recover microgravity results from a network of benchmarks around Kı̄lauea's summit, eliminate instrumental drift, and restore suspected data tares. A total mass of 1.9 × 10 11 kg was determined from these microgravity campaigns to have accumulated below Kı̄lauea Caldera during 2009–2015 at an estimated depth of 1.3 km below sea level. Only a fraction of this mass is reflected in surface deformation, and this is consistent with previously reported discrepancies between subsurface mass accumulation and observed surface deformation. The discrepancy, amongst other independent evidence from gas emissions, seismicity, and continuous gravimetry, indicate densification of magma in the reservoirs below the volcano summit. This densification may have been driven by degassing through the summit vent. It is hypothesized that during the final years of the summit eruption, magma densification resulted in a buildup of pressure in the reservoirs that may have contributed to the lower East Rift Zone outbreak of 2018. The observed mass accumulation beneath Kı̄lauea could not have been detected through other techniques and illustrates the importance of microgravity measurements in volcano monitoring.

Hawai'i↗
USGS2023

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

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

Agung Volcano, Bali↗