Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Sustainability”

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,459 records · Page 81Linked to original sources

Substantial declines in salinity observed across the Upper Colorado River Basin during the 20th century, 1929 to 2019

Salinity in the Colorado River Basin causes an estimated $300 to $400 million per year in economic damages in the U.S. To inform and improve salinity‐control efforts, this study quantifies long‐term trends in salinity (dissolved solids) across the Upper Colorado River Basin (UCRB), including time periods prior to the construction of large dams and preceding the implementation of salinity‐control projects. Weighted Regressions on Time, Discharge, and Season was used with datasets of dissolved‐solids and specific‐conductance measurements, collected as early as 1929, to evaluate long‐term trends in dissolved‐solids loads and concentrations in streams from 1929 to 2019 ( n =14). Results indicate that large, widespread, and sustained downward trends in dissolved‐solids concentrations and loads occurred over the last 50 to 90 years. For 12 of the 14 stream sites with significant downward change, median declines of ‐38% (range of ‐14 to ‐57%) and ‐40% (range of ‐9 to ‐65%) were observed for flow‐normalized concentration and load, respectively. Steepest rates of decline occurred from 1980 to 2000, coincident with the initiation of salinity‐control efforts in the 1980s. However, there was a consistent slowing or reversing of downward trends after 2000 even though salinity‐control efforts continued. Significant decreases in salinity occurred as early as the 1940s at some streams, indicating that, in addition to salinity‐control projects, changes in land cover, land use, and/or climate substantially affect salinity transport in the UCRB. Observed dissolved‐solids trends are likely the result of changes to watershed‐related processes, not due to changes in the streamflow regime.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

Fitting jet noise similarity spectra to volcano infrasound data

Infrasound (low-frequency acoustic waves) has proven useful to detect and characterize subaerial volcanic activity, but understanding the infrasonic source during sustained eruptions is still an area of active research. Preliminary comparison between acoustic eruption spectra and the jet noise similarity spectra suggests that volcanoes can produce an infrasonic form of jet noise from turbulence. The jet noise similarity spectra, empirically derived from audible laboratory jets, consist of two noise sources: large-scale turbulence (LST) and fine-scale turbulence (FST). We fit the similarity spectra quantitatively to eruptions of Mount St. Helens in 2005, Tungurahua in 2006, and Kīlauea in 2018 using nonlinear least squares fitting. By fitting over a wide infrasonic frequency band (0.05–10 Hz) and restricting the peak frequency above 0.15 Hz, we observe a better fit during times of eruption versus non-eruptive background noise. Fitting smaller overlapping frequency bands highlights changes in the fit of LST and FST spectra, which aligns with observed changes in eruption dynamics. Our results indicate that future quantitative spectral fitting of eruption data will help identify changes in eruption source parameters such as velocity, jet diameter, and ash content which are critical for effective hazard monitoring and response.

Earth and Space Science↗

How will baseflow respond to climate change in the Upper Colorado River Basin?

Baseflow is critical to sustaining streamflow in the Upper Colorado River Basin. Therefore, effective water resources management requires estimates of baseflow response to climatic changes. This study provides the first estimates of projected baseflow changes from historical (1984 – 2012) to thirty-year periods centered around 2030, 2050, and 2080 under warm/wet, median, and hot/dry climatic conditions using a hybrid statistical-deterministic baseflow model. Total baseflow supplied to the Lower Colorado River Basin may decline by up to 33%, although this value may increase in the near future by 6% under warm/wet conditions. The percentage of baseflow lost during in-stream transport is projected to increase by 1 - 5% relative to historical conditions. Results highlight that climate driven changes in high elevation hydrology have impacts on basin-wide water availability. Study results have implications for human and ecological water availability in one of the most heavily managed watersheds in the world.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

Flow development and entrainment in turbulent particle-laden jets

Explosive eruptions expel volcanic gases and particles at high pressures and velocities. Within this multiphase fluid, small ash particles affect the flow dynamics, impacting mixing, entrainment, turbulence, and aggregation. To examine the role of turbulent particle behavior, we conducted an analogue experiment using a particle-laden jet. We used compressed air as the carrier fluid, considering turbulent conditions at Reynolds numbers from approximately 5,000 to 20,000. Two different particles were examined: 14-μm diameter solid nickel spheres and 13-μm diameter hollow glass spheres. These resulted in Stokes numbers between 1 and 35 based on the convective scale. The particle mass percentage in the mixture is varied from 0.3% to more than 20%. Based on a 1-D volcanic plume model, these Stokes numbers and mass loadings corresponded to millimeter-scale particle diameters at heights of 4–8 km above the vent during large, sustained eruptions. Through particle image velocimetry, we measured the mean flow behavior and the turbulence statistics in the near-exit region, primarily focusing on the dispersed phase. We show that the flow behavior is dominated by the particle inertia, with high Stokes numbers reducing the entrainment by more than 40%. When applied to volcanic plumes, these results suggest that high-density particles can greatly increase the probability of column collapse.

Journal of Geophysical Research: Atmospheres↗

Hydration state and rheologic stratification of the lithospheric mantle beneath the North Anatolian Fault, Turkey

We present constraints on the hydration state and rheology of the lithospheric mantle beneath the North Anatolian fault zone (NAFZ). Peridotite xenoliths from the Biyikali and Çorlu volcanic centers record deformational microstructures consistent with shearing in a lithosphere-scale transcurrent fault system. Analysis by Fourier transform infrared spectroscopy indicates that nominally anhydrous phases retain some OH − , but bulk rock concentrations are generally restricted to <50 ppm H 2 O by weight. From the rock microstructure, we determined differential stress magnitude and active deformation mechanism(s); combined with estimates of hydration state, we constrained the rheology. Recrystallized grain size piezometry shows that the mantle beneath the NAFZ sustained differential stresses of 10–20 MPa, largely independent of depth. The dominant deformation mechanism(s) change with depth; xenoliths extracted from shallower depths record evidence for grain size-sensitive creep possibly in the presence of melt. At intermediate depths, both dislocation creep and grain size-sensitive mechanisms were active, and we did not observe evidence for deformation in the presence of melt. The deepest samples were dominated by dislocation creep. The strong temperature sensitivity of creep mechanisms, combined with the low variability in differential stress, contributes to a stratified viscosity profile ranging from 10 18 Pa s for the deepest samples to >10 22 Pa s at shallower depths (assuming a melt-free rheology). Although difficult to quantify from the rock record, melt likely reduced the viscosity of the shallow lithospheric mantle. Vertical stratification in viscosity beneath the NAFZ, the result of melt-present deformation and/or transitions in the dominant deformation mechanism, has important consequences for the seismic cycle of strike-slip fault systems.

North Anatolian fault zone↗

A watershed moment for western U.S. dams

The summer of 2023 is a notable time for water-resource management in the western United States: Glen Canyon Dam, on the Colorado River, turns 60 years old while the largest dam-removal project in history is beginning on the Klamath River. This commentary discusses these events in the context of a changing paradigm for dam and reservoir management in this region. Since the era of large dam building began to wane six decades ago, new challenges have arisen for dam and reservoir management owing to climate change, population increase, reservoir sedimentation, declining safety of aging dams, and more environmentally focused management objectives. Today we also better understand dams' benefits, costs, and environmental impacts, including some that were unforeseen and took decades to become apparent. Where dams have become unsafe, obsolete (e.g., due to excessive reservoir sedimentation), and uneconomical beyond saving, dam removal has become common. The science and practice of dam removal are accelerating rapidly, and some long-term physical and biological response studies are now available. Removal of four hydroelectric dams on the Klamath River will be a larger and more complex project than any previous dam removal. The imminency of this project reflects a very different situation for dam and reservoir management than 60 years ago. Looking forward, dam and reservoir management in the western United States and worldwide will require continued collaboration and innovative thinking to meet a wide range of objectives and to manage water resources sustainably for future generations.

Water Resources Research↗

Seasonal and species-level water-use strategies and groundwater dependence in dryland riparian woodlands during extreme drought

Drought-induced groundwater decline and warming associated with climate change are primary threats to dryland riparian woodlands. We used the extreme 2012–2019 drought in southern California as a natural experiment to assess how differences in water-use strategies and groundwater dependence may influence the drought susceptibility of dryland riparian tree species with overlapping distributions. We analyzed tree-ring stable carbon and oxygen isotopes collected from two cottonwood species ( Populus trichocarpa and P . fremontii ) along the semi-arid Santa Clara River. We also modeled tree source water δ 18 O composition to compare with observed source water δ 18 O within the floodplain to infer patterns of groundwater reliance. Our results suggest that both species functioned as facultative phreatophytes that used shallow soil moisture when available but ultimately relied on groundwater to maintain physiological function during drought. We also observed apparent species differences in water-use strategies and groundwater dependence related to their regional distributions. P . fremontii was constrained to more arid river segments and ostensibly used a greater proportion of groundwater to satisfy higher evaporative demand. P . fremontii maintained ∆ 13 C at pre-drought levels up until the peak of the drought, when trees experienced a precipitous decline in ∆ 13 C. This response pattern suggests that trees prioritized maintaining photosynthetic processes over hydraulic safety, until a critical point. In contrast, P . trichocarpa showed a more gradual and sustained reduction in ∆ 13 C, indicating that drought conditions induced stomatal closure and higher water use efficiency. This strategy may confer drought avoidance for P . trichocarpa while increasing its susceptibility to anticipated climate warming.

California↗

Next generation public supply water withdrawal estimation for the conterminous United States using machine learning and operational frameworks

Estimation of human water withdrawals is more important now than ever due to uncertain water supplies, population growth, and climate change. Fourteen percent of the total water withdrawal in the United States is used for public supply, typically including deliveries to domestic, commercial, and occasionally including industrial, irrigation, and thermoelectric water withdrawal. Stewards of water resources in the USA require estimates of water withdrawals to manage and plan for future demands and sustainable water supplies. This study compiled the most comprehensive conterminous United States water withdrawal data set to date and developed a machine learning framework for estimating public supply withdrawals and associated uncertainty for the period 2000–2020. The modeling approach provides service area resolution estimates to allow for annual and monthly water withdrawal estimation while incorporating a complex array of driving factors that include hydroclimatic, demographic, socioeconomic, geographic, and land use factors. Model results reveal highly variable and lognormally distributed per-capita water withdrawal, spanning from 30 to 650 gallons per capita per day (GPCD), across community, regional, and national scales, with pronounced seasonal variations. Analysis of estimated withdrawal trends indicates that the national annual average withdrawal experienced a decline at a rate of 0.58 GPCD/year during the period from 2000 to 2020. Model interpretation reveals a complex interplay between public supply withdrawal and key predictors, including population size, warm-season precipitation, counts of large buildings and houses, and areas of urban and commercial land use. The developed models can forecast future public supply driven by various climate, demographic, and socioeconomic scenarios.

conterminous United States↗

The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought

Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi‐year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage‐limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.

Colorado↗

The United States Magnetotelluric Array and the National Impedance Map

The United States Magnetotelluric Array (USMTArray) data set, collected in the years 2006–2024, consists of more than 1,700 long-period magnetotelluric stations covering the entirety of the contiguous United States on a quasi-regular 70 km grid. Funding across multiple federal agencies was critical to sustaining this effort to its completion. Important components of the project included active guidance and participation from the MT community, the open and timely availability of all data, and the application of consistent instrumentation and robust data processing. Together with parallel advancement in the development of publicly available three-dimensional (3D) inversion codes, the USMTArray has revitalized the US magnetotelluric community and increased the visibility of magnetotellurics within the Earth-science community. Taken as a whole, these data are visualized as the National Impedance Map, which, together with a 3D synthesis conductivity model of the nation, reveals the electrical architecture of the contiguous US. USMTArray data are used by researchers worldwide for fundamental and applied studies, including investigations of continental architecture and evolution, estimation of hazards to critical infrastructure due to geomagnetic storms, and assessment of the nation's undiscovered geothermal and mineral resources. We here review the history and development of the project, discuss the challenges and successes in its execution, present the National Impedance Map and synthesis conductivity model, and highlight the breadth of research stemming from this rich data set.

conterminous United States↗

Patterns and drivers of cliff erosion in Big Sur, California, USA using repeat photogrammetry, 2017–2023

Seacliff erosion in steep terrain poses major risks to transportation and critical infrastructure. In Big Sur, California, USA, seacliff erosion threatens the sustainability of the central coast stretch of California State Route 1, a transportation corridor that is critical to the region's economy. Published cliff retreat rates for the region range from 1 to 40 cm yr −1 , highlighting that high-resolution, process-based studies could enhance understanding of the causes of spatial and temporal variability. We quantified cliff erosion and investigated its drivers along ∼13 km of the Big Sur coastline at week–month timescales during the late fall to early spring wet seasons between January 2017 and June 2023 by analyzing 3D point clouds developed from aerial imagery using four-dimensional structure-from-motion (4D SfM) photogrammetry techniques. We calculated cliff face retreat rates of 2.23 ± 3.06 cm yr −1 (mean ±1 σ ), an order of magnitude lower than long-term estimated rates for the region (which included large deep-seated landslides), but in line with short-term rates reported across California. Change detection imagery comparison, cliff profiles through time, and statistical analysis reveal a cyclical cliff evolution process in which erosion by wave action at the cliff base destabilizes the cliff and primes it for subsequent failure during precipitation events. Although more erosion by volume could be attributed to precipitation-induced increases in soil moisture (784 m 3 km −1 yr −1 ) compared with erosion attributed to wave power (282 m 3 km −1 yr −1 ), our observations underscore the coupled nature of these processes in driving cliff evolution, consistent with established theory and observations.

California↗

Multi-stage crustal thickening, surface uplift, and collision in the western Himalaya-Karakoram-Tibet orogen revealed by chemical mohometry

The Himalaya-Karakoram-Tibet (HKT) orogen provides an unrivaled opportunity to study the dynamic linkages between deep and surface processes during collisional orogenesis. However, these efforts are hindered by conflicting interpretations on the number and timing of collisional events, and the timing of crustal thickening and associated surface uplift. Here, we resolve this with quantitative paleo-crustal thickness estimates in the northwestern HKT orogen. We show that: (a) the paleo-Asian margin had thick crust (50–60 km) at least 65 Ma prior to terminal collision, consistent with a continental arc setting, (b) crustal thickening to 60 km or more occurred at ca. 60–50 Ma in the Kohistan-Ladakh arc and by 40–25 Ma in the paleo-Asian margin, indicating a multi-stage Himalayan collision, and (c) modern crustal thicknesses in the northwestern HKT have been sustained since ca. 40–25 Ma suggesting an orogenic steady-state in which crustal thickening, crustal flow, and surface uplift have been balanced by erosion.

Geophysical Research Letters↗

Decadal shifts in groundwater age detected by environmental tracers across California, USA

Groundwater age offers important insight into recharge, storage, and contamination risk. Although models predict age changes can be driven by pumping and climate variability, direct observational evidence remains limited. Here, we analyzed paired environmental tracer suites (tritium, carbon-14, and tritiogenic helium-3) collected a decade apart from 268 wells across California to assess the prevalence of groundwater age transience. Travel-time distribution models and statistical tests indicated age transience at 29% of sites, occurring most often in agricultural regions, such as the San Joaquin Valley and Southern Coast Ranges, where large carbon-14 changes coincided with substantial nitrate and chloride shifts. Sites with tritiogenic helium-3 data showed more frequent age transience, underscoring the value of multi-tracer data sets. These results provide the first regional evidence of widespread groundwater age change and a method for detecting changing water balances with implications for groundwater sustainability and water quality.

California↗

Rupture continuity through intermittent pauses in Cascadia slow slip events

Cascadia slow slip events (SSEs) are often envisioned as smooth, continuous ruptures, progressively activating tremor asperities as they propagate. Macroscopically, geodetic inversions and spatiotemporal maps of tremor epicenters show steady, uniform migration. In detail tremor is more chaotic and discontinuous. Larger long-term SSEs observed in daily geodetic solutions are inferred to exhibit intermittent pauses that reflect temporary re-locking of the fault, but this temporal resolution limits tests for similar re-locking on shorter timescales. We use temporal measurements of the areal growth and radiated energy of tremor clusters to investigate SSE intermittence. We find that ruptures mirror tremor pauses. Areal growth rate, however, does not reset, and removing the pauses results in smoother and more similar growth measurements among all SSEs. The rupture similarity occurs regardless of size or location and hints at an underlying uniformity and lack of predeterminism in eventual SSE size. Epicentral uncertainty precludes quantifying early rupture stages, but for larger events areal growth follows a power-law and slows with increasing size. Temporal correlations in tremor energy with inferred SSE propagation velocities and tremor rates suggest its use as a proxy for slip velocity. We find that tremor energy is tidally modulated at daily and sub-daily frequencies, and this modulation is continuous through pauses, suggesting a memory of slip state is sustained through them. We argue these pauses reflect unsteady propagation of the slip front, marked by rapid re- and un-locking, and excluding them removes rupture complexity to reveal a diffusive-like slip process and underlying universality in growth.

California, Oregon, Washington↗

Towards affordable wetland evapotranspiration monitoring using the Variance-Bowen Ratio method: Insights from three contrasting wetlands

Accurate measurement of evapotranspiration (ET) is essential for sustainable water management. Standard methods such as eddy covariance (EC) are costly, while alternatives such as surface renewal are cheaper but require calibration and complex data processing. This study evaluates the utility of the Variance-Bowen Ratio (VBR) method for estimating ET across three California’s wetlands. Using data from 2023, half-hourly latent heat flux (λE) and daily/monthly ET from VBR were compared with EC at one non-tidal (site A) and two tidal (sites B and C) wetlands. λE VBR consistently underestimated λE EC , with root mean squared errors (RMSE) of 61.2 W m -2 at sites A, 106.1 W m -2 at site B, and 137.2 W m -2 at site C, largely due to storage fluxes across sites. Temporal integration improved VBR’s performance at tidal sites, where compensating water heat storage errors yielded low daily and monthly biases (site B: RMSE = 0.78 mm/d and 12 mm/month; r = 0.93; site C: RMSE = 0.90 mm/d and 13 mm/month; r = 0.93), with reduced major axis (RMA) regression slopes of 0.98 and ~0.91. In contrast, biomass heat storage at site A caused persistent biases (RMSEs = 0.97 mm/d and 23 mm/month; daily and monthly RMA slopes ~0.75; r = 0.85). These results highlight VBR’s limitations in environments with substantial storage fluxes. Despite this, VBR is cost-effective for estimating daily and monthly ET, with sensor costs at least tenfold lower than EC and simpler setup, making it suitable for ET monitoring in resource-limited and hard-to-access regions.

California↗

Tide gage measurements of uplift along the south coast of Alaska

Annual mean sea levels along the south coast of Alaska are used to measure uplift along the Alaska‐Aleutian subduction zone. Oceanographic effects are removed from the observed annual mean sea levels by subtracting a correction that is proportional to the sea level fluctuations observed in southeast Alaska. That correction is effective in reducing fluctuations in the observed, annual mean sea level as far west as the tip of Alaska peninsula. Additional corrections to remove the eustatic rise in sea level and the apparent fall in sea level due to postglacial isostatic rebound of the land are introduced. This corrected sea level record should provide a measure of tectonic subsidence. In the area affected by the 1964 Alaska earthquake, postseismic uplift occurs where coseismic subsidence was observed, and postseismic subsidence occurs where coseismic uplift was observed. The immediate postseismic response is damped out within the first decade, and the subsequent uplift rates appear to be steady over the 1974–1989 interval. However, some of those rates seem to be too high to be sustained over the ∼1000 year earthquake recurrence interval appropriate to this area if the interseismic deformation is only to recover the coseismic displacement. Thus a long‐term ( ∼100 years) relaxation in uplift rates is postulated. The immediate (time constant ∼5 years) postseismic relaxation is attributed to postseismic slip on the plate interface directly downdip from the coseismic rupture. The long‐term (time constant ∼100 years) relaxation is attributed to flow in the asthenosphere.

Alaska↗

Application of continuum models to deformation of the Aleutian Island Arc

Continuum models were constructed to describe large‐scale deformation of the Aleutian Island Arc over the past 5 m.y. These models consider the island arc as a continuum in the horizontal plane with the velocity boundary condition at the Pacific edge stated as a fraction of Pacific plate convergence transferred to the arc. First, a simple model of uniformly distributed strain is formulated to illustrate the mechanics of continuous deformation. Lineaments along the arc massif rotated about a vertical axis are matched by small‐element rotation calculated from the model. However, this model does not predict across‐arc variations in deformation and produces an unrealistic amount of crustal thickening after 5 m.y. A physically more meaningful model of deformation is the thin viscous sheet model based on averages of stress and rheology throughout the lithosphere. The amount of motion transferred from the Pacific plate to the arc is constrained by the rotated lineaments, while the effective stress‐strain exponent ( n ) and the ability the lithosphere has to sustain crustal thickness contrasts (the Argand number) are independent variables. Primarily, bathymetry, earthquake focal mechanisms, and styles of faulting are used to evaluate the models. The preferred model is one where the amount of motion transferred from the Pacific plate is greater in an arc‐parallel direction than in an arc‐normal direction, producing stresses consistent with strike‐slip faulting at the far western end of the arc and tensional stresses consistent with transverse normal faulting elsewhere in the arc massif. This model agrees with observations of slip vectors by Ekström and Engdahl (1989), who conclude that a portion of the arc‐parallel component of relative plate motion is taken up in the overriding plate. This model implies that compressive stress transferred to the arc is small in comparison to along‐arc shear stress and that stresses conducive to strikeslip faulting are prevalent throughout the arc.

Alaska↗

An episode of reinflation of the Long Valley Caldera, eastern California: 1989-1991

Following the episodes of inflation of the resurgent dome associated with the May 1980 earthquake sequence (four M 6 earthquakes) and the January 1983 earthquake swarm (two M 5.2 events), 7 years of frequently repeated two-color geodimeter measurements spanning the Long Valley caldera document gradually decreasing extensional strain rates from 5 ppm/yr in mid-1983, when the measurements began, to near zero in mid-1989. The corresponding seismic activity within the caldera persisted at a low rate of fewer than 10 M ≥ 1.2 earthquakes per week from 1985 through November 1989 with no events exceeding M 3.0. Early October 1989 marked a change in activity when measurements of the two-color geodimeter network showed a significant increase in extensional strain rate (9 ppm/yr) across the caldera. The seismic activity began exceeding 10 M ≥ 1.2 per week in early December 1989 and rapidly increased to a sustained level of tens of M ≥ 1.2 per week with bursts having hundreds of events per day. Many events exceeded M 3.0 and the largest event was M ≈ 4. The 1989–1991 inflation episode is the first time that we have sufficient geodetic measurements in Long Valley to define the temporal relation between onset of an inflation episode and onset of brittle failure (earthquake swarm within the caldera). Here, the onset of deformation preceded the onset of increased earthquake activity by more than 2 months. The seismicity rate began to decrease in mid-July 1990, consistent with a gradually slowing of extension across the caldera as measured by the two-color geodimeter. The recent episode of inflation can be modeled by a single Mogi point source located about 7 km beneath the center of the resurgent dome. In contrast, the deformation pattern observed between mid-1983 and mid-1989 is best reproduced by fault slip in the south moat, inflation at 6.5 km depth near Casa Diablo Hot Springs and inflation beneath the resurgent dome. It appears that the 7-km source beneath the resurgent dome that was active for the earlier episodes is the primary source for the more recent episode. The model used to satisfy the line length observations predicts 7.5 cm of uplift along leveling route along highway 395 from mid-1983 to mid-1989 and an additional 11 cm through the end of 1991. To comp are with the energy release from seismicity, the modeled inflation from late 1989 through the end of 1991 has a moment that is a factor of 40 more than the cumulative seismic moment from earthquakes located within the caldera from the same period. Thus the recent inflation episode represents a significant portion of the observed geodetic deformation with only little seismic energy release.

Journal of Geophysical Research↗