USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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 253 records · Page 14Linked to original sources
USGS2017

Magnetotelluric imaging of lower crustal melt and lithospheric hydration in the Rocky Mountain Front transition zone, Colorado, USA

We present an electrical resistivity model of the crust and upper mantle from two‐dimensional (2‐D) anisotropic inversion of magnetotelluric data collected along a 450 km transect of the Rio Grande rift, southern Rocky Mountains, and High Plains in Colorado, USA. Our model provides a window into the modern‐day lithosphere beneath the Rocky Mountain Front to depths in excess of 150 km. Two key features of the 2‐D resistivity model are (1) a broad zone (~200 km wide) of enhanced electrical conductivity (<20 Ωm) in the midcrust to lower crust that is centered beneath the highest elevations of the southern Rocky Mountains and (2) hydrated lithospheric mantle beneath the Great Plains with water content in excess of 100 ppm. We interpret the high conductivity region of the lower crust as a zone of partially molten basalt and associated deep‐crustal fluids that is the result of recent (less than 10 Ma) tectonic activity in the region. The recent supply of volatiles and/or heat to the base of the crust in the late Cenozoic implies that modern‐day tectonic activity in the western United States extends to at least the western margin of the Great Plains. The transition from conductive to resistive upper mantle is caused by a gradient in lithospheric modification, likely including hydration of nominally anhydrous minerals, with maximum hydration occurring beneath the Rocky Mountain Front. This lithospheric “hydration front” has implications for the tectonic evolution of the continental interior and the mechanisms by which water infiltrates the lithosphere.

Colorado↗
USGS2017

Decaying lava extrusion rate at El Reventador Volcano, Ecuador measured using high-resolution satellite radar

Lava extrusion at erupting volcanoes causes rapid changes in topography and morphology on the order of tens or even hundreds of meters. Satellite radar provides a method for measuring changes in topographic height over a given time period to an accuracy of meters, either by measuring the width of radar shadow cast by steep sided features, or by measuring the difference in radar phase between two sensors separated in space. We measure height changes, and hence estimate extruded lava volume flux, at El Reventador, Ecuador, between 2011 and 2016, using data from the RADARSAT-2 and TanDEM-X satellite missions. We find that 39 new lava flows were extruded between 9 February 2012 and 24 August 2016, with a cumulative volume of 44.8M m 3 dense rock equivalent and a gradually decreasing eruption rate. The average dense rock rate of lava extrusion during this time is 0.31 ± 0.02 m 3 s −1 , which is similar to the long-term average from 1972 to 2016. Apart from a volumetrically small dyke opening event between 9 March and 10 June 2012, lava extrusion at El Reventador is not accompanied by any significant magmatic ground deformation. We use a simple physics-based model to estimate that the volume of the magma reservoir under El Reventador is greater than 3 km 3 . Our lava extrusion data can be equally well fit by models representing a closed reservoir depressurising during the eruption with no magma recharge, or an open reservoir with a time-constant magma recharge rate of up to 0.35 ± 0.01 m 3 s −1 .

El Reventador Volcano↗
USGS2018

Temporal stress changes caused by earthquakes: A review

Earthquakes can change the stress field in the Earth’s lithosphere as they relieve and redistribute stress. Earthquake-induced stress changes have been observed as temporal rotations of the principal stress axes following major earthquakes in a variety of tectonic settings. The stress changes due to the 2011 Mw9.0 Tohoku-Oki, Japan, earthquake were particularly well documented. Earthquake stress rotations can inform our understanding of earthquake physics, most notably addressing the long-standing problem of whether the Earth’s crust at plate boundaries is “strong” or “weak.” Many of the observed stress rotations, including that due to the Tohoku-Oki earthquake, indicate near-complete stress drop in the mainshock. This implies low background differential stress, on the order of earthquake stress drop, supporting the weak crust model. Earthquake stress rotations can also be used to address other important geophysical questions, such as the level of crustal stress heterogeneity and the mechanisms of postseismic stress reloading. The quantitative interpretation of stress rotations is evolving from those based on simple analytical methods to those based on more sophisticated numerical modeling that can capture the spatial-temporal complexity of the earthquake stress changes.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Source spectral properties of small-to-moderate earthquakes in southern Kansas

The source spectral properties of injection-induced earthquakes give insight into their nucleation, rupture processes, and influence on ground motion. Here we apply a spectral decomposition approach to analyze P -wave spectra and estimate Brune-type stress drop for more than 2000 M L 1.5–5.2 earthquakes occurring in southern Kansas from 2014 to 2016. We find that these earthquakes are characterized by low stress drop values (median ∼0.4MPa) compared to natural seismicity in California. We observe a significant increase in stress drop as a function of depth, but the shallow depth distribution of these events is not by itself sufficient to explain their lower stress drop. Stress drop increases with magnitude from M1.5–M3.5, but this scaling trend may weaken above M4 and also depends on the assumed source model. Although we observe a nonstationary, sequence-specific temporal evolution in stress drop, we find no clear systematic relation with the activity of nearby injection wells.

Kansas↗
USGS2018

Deformation of the Pacific/North America plate boundary at Queen Charlotte Fault: The possible role of rheology

The Pacific/North America (PA/NA) plate boundary between Vancouver Island and Alaska is similar to the PA/NA boundary in California in its kinematic history and the rate and azimuth of current relative motion, yet their deformation styles are distinct. The California plate boundary shows a broad zone of parallel strike slip and thrust faults and folds, whereas the 49‐mm/yr PA/NA relative plate motion in Canada and Alaska is centered on a single, narrow, continuous ~900‐km‐long fault, the Queen Charlotte Fault (QCF). Using gravity analysis, we propose that this plate boundary is centered on the continent/ocean boundary (COB), an unusual location for continental transform faults because plate boundaries typically localize within the continental lithosphere, which is weaker. Because the COB is a boundary between materials of contrasting elastic properties, once a fault is established there, it will probably remain stable. We propose that deformation progressively shifted to the COB in the wake of Yakutat terrane's northward motion along the margin. Minor convergence across the plate boundary is probably accommodated by fault reactivation on Pacific crust and by an eastward dipping QCF. Underthrusting of Pacific slab under Haida Gwaii occurs at convergence angles >14°–15° and may have been responsible for the emergence of the archipelago. The calculated slab entry dip (5°–8°) suggests that the slab probably does not extend into the asthenosphere. The PA/NA plate boundary at the QCF can serve as a structurally simple site to investigate the impact of rheology and composition on crustal deformation and the initiation of slab underthrusting.

Journal of Geophysical Research B: Solid Earth↗
USGS2018

Mass addition at Mount St. Helens, Washington, inferred from repeated gravity surveys

Measurements of subtle changes in the Earth’s gravityfield can provideinformation on the addition/loss of mass (e.g., magma or aqueousfluids) beneath a volcano. In this study, wemeasured gravity at Mount St. Helens from 2010 to 2016 to investigate possible mass changes followingthe 2004–2008 dome-forming eruption. The raw gravity measurements were corrected for changes in themass and shape of Crater Glacier and for deformation of the volcanic edifice to obtain residual gravity values.We found positive residual gravity changes that we interpret as evidence for partial recharge of the magmareservoir that fed the 2004–2008 eruption and/or for accumulation of groundwater in one or more shallowaquifers. Most of the gravity signal can be explained by groundwater accumulation; magma recharge is notrequired. However, there is additional evidence for at least partial recharge from seismic, deformation, andgeochemical observations

Washington↗
USGS2018

Constraints on friction, dilatancy, diffusivity, and effective stress from low-frequency earthquake rates on the deep San Andreas Fault

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor on the San Andreas Fault in central California are sensitive to tidal stresses. LFEs occur at all levels of the tides, are strongly correlated and in phase with the ~200 Pa shear stresses, and weakly and not systematically correlated with the ~2 kPa tidal normal stresses. We assume that LFEs are small sources that repeatedly fail during shear within a much larger scale, aseismically slipping fault zone and consider two different models of the fault slip: (1) modulation of the fault slip rate by the tidal stresses or (2) episodic slip, triggered by the tides. LFEs are strongly clustered with duration much shorter than the semidiurnal tide; they cannot be significantly modulated on that time scale. The recurrence times of clusters, however, are many times longer than the semidiurnal, leading to an appearance of tidal triggering. In this context we examine the predictions of laboratory-observed triggered frictional (dilatant) fault slip. The undrained end-member model produces no sensitivity to the tidal normal stress, and slip onsets are in phase with the tidal shear stress. The tidal correlation constrains the diffusivity to be less than ~1 × 10 −6 /s and the product of the friction and dilatancy coefficients to be at most 5 × 10 −7 , orders of magnitude smaller than observed at room temperature. In the absence of dilatancy the effective normal stress at failure would be about ~55 kPa. For this model the observations require intrinsic weakness, low dilatancy, and lithostatic pore fluid.

Journal of Geophysical Research↗
USGS2018

Combining multiphase groundwater flow and slope stability models to assess stratovolcano flank collapse in the Cascade Range

Hydrothermal alteration can create low‐permeability zones, potentially resulting in elevated pore‐fluid pressures, within a volcanic edifice. Strength reduction by rock alteration and high pore‐fluid pressures have been suggested as a mechanism for edifice flank instability. Here we combine numerical models of multiphase heat transport and groundwater flow with a slope‐stability code that incorporates three‐dimensional distributions of strength and pore‐water pressure to address the following questions: (1) What permeability distributions and contrasts produce elevated pore‐fluid pressures in a stratovolcano? (2) What are the effects of these elevated pressures on flank stability? (3) Finally, what are the effects of magma intrusion on potential flank failure in an edifice? Simulation results show that under a range of plausible parameters, water tables in a stratovolcano can be elevated or perched. These elevated water tables result in universally lower stability (lower factor of safety) compared with equivalent dry edifices, indicating a higher likelihood of flank collapse. Low‐permeability (<1 × 10 −17 m 2 ) layers such as altered pyroclastic deposits or breccias can result in locally saturated regions (perched water) and lower factors of safety near the ground surface but may actually reduce liquid water saturation and pore pressures in the core of the edifice and thus may favor small, shallow collapses over larger, deeper collapses. Magma intrusion into the base of the edifice increases pore‐fluid pressures and decreases the factor of safety. However, the shear strength of edifice rocks also exerts a significant control on stability, so both mechanical properties and pore‐fluid pressures are important for stability assessments.

Washington↗
USGS2018

Seismicity induced by wastewater injection in Washington County, Ohio: Influence of preexisting structure, regional stress regime, and well operations

Recent seismicity in Washington County, Ohio, has been suggested to be induced by wastewater disposal operations despite injection ~2 km above the Precambrian basement. We investigated the relationships between disposal well locations and operational histories, spatiotemporal patterns of seismicity enhanced by waveform correlation, and mapped subsurface structures. We also analyzed proxies for in situ stress conditions to evaluate existing principal stress orientations and magnitudes. Double‐difference relocations of small ( M L ≤ 2.1) seismicity revealed linear clusters defining a NE‐SW orientation that parallels the trend of small‐amplitude folds in Upper Paleozoic rocks within 1–2 km of the injection well. Hypocenter depths of ~3.1 to 4.3 km indicate seismicity occurred along faults below the Silurian injection interval in lower Paleozoic strata and Precambrian crystalline basement. Similar folds/faults imaged in nearby seismic reflection lines were associated with basement‐involved fault systems that cut the injection interval targeted by the injection well, providing permeability pathways for fluid pressure increases that initiated slip. Interpreted faults were also optimally oriented in the regional stress field for reactivation due to injection‐related decreases in effective stress. Well completion records suggest that during the recorded seismicity, the Long Run‐1 well was injecting disposal fluids close to or above the fracture pressure of the targeted Silurian injection interval. A comparison of observed seismicity rates and monthly injection records from the disposal well also indicate a close correspondence between earthquake frequency and injection rate. Taken together, these results provide a physical framework for understanding why wastewater disposal operations induced the Washington County seismic sequence.

Ohio↗
USGS2017

Interactions of estuarine shoreline infrastructure with multiscale sea level variability

Sea level rise increases the risk of storms and other short‐term water‐rise events, because it sets a higher water level such that coastal surges become more likely to overtop protections and cause floods. To protect coastal communities, it is necessary to understand the interaction among multiday and tidal sea level variabilities, coastal infrastructure, and sea level rise. We performed a series of numerical simulations for San Francisco Bay to examine two shoreline scenarios and a series of short‐term and long‐term sea level variations. The two shoreline configurations include the existing topography and a coherent full‐bay containment that follows the existing land boundary with an impermeable wall. The sea level variability consists of a half‐meter perturbation, with duration ranging from 2 days to permanent (i.e., sea level rise). The extent of coastal flooding was found to increase with the duration of the high‐water‐level event. The nonlinear interaction between these intermediate scale events and astronomical tidal forcing only contributes ∼1% of the tidal heights; at the same time, the tides are found to be a dominant factor in establishing the evolution and diffusion of multiday high water events. Establishing containment at existing shorelines can change the tidal height spectrum up to 5%, and the impact of this shoreline structure appears stronger in the low‐frequency range. To interpret the spatial and temporal variability at a wide range of frequencies, Optimal Dynamic Mode Decomposition is introduced to analyze the coastal processes and an inverse method is applied to determine the coefficients of a 1‐D diffusion wave model that quantify the impact of bottom roughness, tidal basin geometry, and shoreline configuration on the high water events

Journal of Geophysical Research: Oceans↗
USGS2017

Persistent shoreline shape induced from offshore geologic framework: Effects of shoreface connected ridges

Mechanisms relating offshore geologic framework to shoreline evolution are determined through geologic investigations, oceanographic deployments, and numerical modeling. Analysis of shoreline positions from the past 50 years along Fire Island, New York, a 50 km long barrier island, demonstrates a persistent undulating shape along the western half of the island. The shelf offshore of these persistent undulations is characterized with shoreface-connected sand ridges (SFCR) of a similar alongshore length scale, leading to a hypothesis that the ridges control the shoreline shape through the modification of flow. To evaluate this, a hydrodynamic model was configured to start with the US East Coast and scale down to resolve the Fire Island nearshore. The model was validated using observations along western Fire Island and buoy data, and used to compute waves, currents and sediment fluxes. To isolate the influence of the SFCR on the generation of the persistent shoreline shape, simulations were performed with a linearized nearshore bathymetry to remove alongshore transport gradients associated with shoreline shape. The model accurately predicts the scale and variation of the alongshore transport that would generate the persistent shoreline undulations. In one location, however, the ridge crest connects to the nearshore and leads to an offshore-directed transport that produces a difference in the shoreline shape. This qualitatively supports the hypothesized effect of cross-shore fluxes on coastal evolution. Alongshore flows in the nearshore during a representative storm are driven by wave breaking, vortex force, advection and pressure gradient, all of which are affected by the SFCR.

New York↗
USGS2017

A Bayesian-based system to assess wave-driven flooding hazards on coral reef-lined coasts

Many low-elevation, coral reef-lined, tropical coasts are vulnerable to the effects of climate change, sea level rise, and wave-induced flooding. The considerable morphological diversity of these coasts and the variability of the hydrodynamic forcing that they are exposed to make predicting wave-induced flooding a challenge. A process-based wave-resolving hydrodynamic model (XBeach Non-Hydrostatic, “XBNH”) was used to create a large synthetic database for use in a “Bayesian Estimator for Wave Attack in Reef Environments” (BEWARE), relating incident hydrodynamics and coral reef geomorphology to coastal flooding hazards on reef-lined coasts. Building on previous work, BEWARE improves system understanding of reef hydrodynamics by examining the intrinsic reef and extrinsic forcing factors controlling runup and flooding on reef-lined coasts. The Bayesian estimator has high predictive skill for the XBNH model outputs that are flooding indicators, and was validated for a number of available field cases. It was found that, in order to accurately predict flooding hazards, water depth over the reef flat, incident wave conditions, and reef flat width are the most essential factors, whereas other factors such as beach slope and bed friction due to the presence or absence of corals are less important. BEWARE is a potentially powerful tool for use in early warning systems or risk assessment studies, and can be used to make projections about how wave-induced flooding on coral reef-lined coasts may change due to climate change.

Journal of Geophysical Research C: Oceans↗
USGS2018

Direct measurements of mean Reynolds stress and ripple roughness in the presence of energetic forcing by surface waves

Direct covariance observations of the mean flow Reynolds stress and sonar images of the seafloor collected on a wave‐exposed inner continental shelf demonstrate that the drag exerted by the seabed on the overlying flow is consistent with boundary layer models for wave‐current interaction, provided that the orientation and anisotropy of the bed roughness are appropriately quantified. Large spatial and temporal variations in drag result from nonequilibrium ripple dynamics, ripple anisotropy, and the orientation of the ripples relative to the current. At a location in coarse sand characterized by large two‐dimensional orbital ripples, the observed drag shows a strong dependence on the relative orientation of the mean current to the ripple crests. At a contrasting location in fine sand, where more isotropic sub‐orbital ripples are observed, the sensitivity of the current to the orientation of the ripples is reduced. Further, at the coarse site under conditions when the currents are parallel to the ripple crests and the wave orbital diameter is smaller than the wavelength of the relic orbital ripples, the flow becomes hydraulically smooth. This transition is not observed at the fine site, where the observed wave orbital diameter is always greater than the wavelength of the observed sub‐orbital ripples. Paradoxically, the dominant along‐shelf flows often experience lower drag at the coarse site than at the fine site, despite the larger ripples, highlighting the complex dynamics controlling drag in wave‐exposed environments with heterogeneous roughness.

Journal of Geophysical Research C: Oceans↗
USGS2017

Presentation and analysis of a worldwide database of earthquake-induced landslide inventories

Earthquake-induced landslide (EQIL) inventories are essential tools to extend our knowledge of the relationship between earthquakes and the landslides they can trigger. Regrettably, such inventories are difficult to generate and therefore scarce, and the available ones differ in terms of their quality and level of completeness. Moreover, access to existing EQIL inventories is currently difficult because there is no centralized database. To address these issues, we compiled EQIL inventories from around the globe based on an extensive literature study. The database contains information on 363 landslide-triggering earthquakes and includes 66 digital landslide inventories. To make these data openly available, we created a repository to host the digital inventories that we have permission to redistribute through the U.S. Geological Survey ScienceBase platform. It can grow over time as more authors contribute their inventories. We analyze the distribution of EQIL events by time period and location, more specifically breaking down the distribution by continent, country, and mountain region. Additionally, we analyze frequency distributions of EQIL characteristics, such as the approximate area affected by landslides, total number of landslides, maximum distance from fault rupture zone, and distance from epicenter when the fault plane location is unknown. For the available digital EQIL inventories, we examine the underlying characteristics of landslide size, topographic slope, roughness, local relief, distance to streams, peak ground acceleration, peak ground velocity, and Modified Mercalli Intensity. Also, we present an evaluation system to help users assess the suitability of the available inventories for different types of EQIL studies and model development.

Journal of Geophysical Research F: Earth Surface↗
USGS2017

Compositional signatures in acoustic backscatter over vegetated and unvegetated mixed sand-gravel riverbeds

Multibeam acoustic backscatter has considerable utility for remote characterization of spatially heterogeneous bed sediment composition over vegetated and unvegetated riverbeds of mixed sand and gravel. However, the use of high-frequency, decimeter-resolution acoustic backscatter for sediment classification in shallow water is hampered by significant topographic contamination of the signal. In mixed sand-gravel riverbeds, changes in the abiotic composition of sediment (such as homogeneous sand to homogeneous gravel) tend to occur over larger spatial scales than is characteristic of small-scale bedform topography (ripples, dunes, and bars) or biota (such as vascular plants and periphyton). A two-stage method is proposed to filter out the morphological contributions to acoustic backscatter. First, the residual supragrain-scale topographic effects in acoustic backscatter with small instantaneous insonified areas, caused by ambiguity in the local (beam-to-beam) bed-sonar geometry, are removed. Then, coherent scales between high-resolution topography and backscatter are identified using cospectra, which are used to design a frequency domain filter that decomposes backscatter into the (unwanted) high-pass component associated with bedform topography (ripples, dunes, and sand waves) and vegetation, and the (desired) low-frequency component associated with the composition of sediment patches superimposed on the topography. This process strengthens relationships between backscatter and sediment composition. A probabilistic framework is presented for classifying vegetated and unvegetated substrates based on acoustic backscatter at decimeter resolution. This capability is demonstrated using data collected from diverse settings within a 386 km reach of a canyon river whose bed varies among sand, gravel, cobbles, boulders, and submerged vegetation.

Colorado River, Grand Canyon, Glen Canyon, Marble ↗
USGS2017

Can beaches survive climate change?

Anthropogenic climate change is driving sea level rise, leading to numerous impacts on the coastal zone, such as increased coastal flooding, beach erosion, cliff failure, saltwater intrusion in aquifers, and groundwater inundation. Many beaches around the world are currently experiencing chronic erosion as a result of gradual, present-day rates of sea level rise (about 3 mm/year) and human-driven restrictions in sand supply (e.g., harbor dredging and river damming). Accelerated sea level rise threatens to worsen coastal erosion and challenge the very existence of natural beaches throughout the world. Understanding and predicting the rates of sea level rise and coastal erosion depends on integrating data on natural systems with computer simulations. Although many computer modeling approaches are available to simulate shoreline change, few are capable of making reliable long-term predictions needed for full adaption or to enhance resilience. Recent advancements have allowed convincing decadal to centennial-scale predictions of shoreline evolution. For example, along 500 km of the Southern California coast, a new model featuring data assimilation predicts that up to 67% of beaches may completely erode by 2100 without large-scale human interventions. In spite of recent advancements, coastal evolution models must continue to improve in their theoretical framework, quantification of accuracy and uncertainty, computational efficiency, predictive capability, and integration with observed data, in order to meet the scientific and engineering challenges produced by a changing climate.

Journal of Geophysical Research F: Earth Surface↗
USGS2017

A framework for modeling connections between hydraulics, water surface roughness, and surface reflectance in open channel flows

This paper introduces a framework for examining connections between the flow field, the texture of the air-water interface, and the reflectance of the water surface and thus evaluating the potential to infer hydraulic information from remotely sensed observations of surface reflectance. We used a spatial correlation model describing water surface topography to illustrate the application of our framework. Nondimensional relations between model parameters and flow intensity were established based on a prior flume study. Expressing the model in the spatial frequency domain allowed us to use an efficient Fourier transform-based algorithm for simulating water surfaces. Realizations for both flume and field settings had water surface slope distributions positively correlated with velocity and water surface roughness. However, most surface facets were gently sloped and thus unlikely to yield strong specular reflections; the model exaggerated the extent of water surface features, leading to underestimation of facet slopes. A ray tracing algorithm indicated that reflectance was greatest when solar and view zenith angles were equal and the sensor scanned toward the Sun to capture specular reflections of the solar beam. Reflected energy was concentrated in a small portion of the sky, but rougher water surfaces reflected rays into a broader range of directions. Our framework facilitates flight planning to avoid surface-reflected radiance while mapping other river attributes, or to maximize this component to exploit relationships between hydraulics and surface reflectance. This initial analysis also highlighted the need for improved models of water surface topography in natural rivers.

Journal of Geophysical Research F: Earth Surface↗
USGS2018

Using stereo satellite imagery to account for ablation, entrainment, and compaction in volume calculations for rock avalanches on Glaciers: Application to the 2016 Lamplugh Rock Avalanche in Glacier Bay National Park, Alaska

The use of preevent and postevent digital elevation models (DEMs) to estimate the volume of rock avalanches on glaciers is complicated by ablation of ice before and after the rock avalanche, scour of material during rock avalanche emplacement, and postevent ablation and compaction of the rock avalanche deposit. We present a model to account for these processes in volume estimates of rock avalanches on glaciers. We applied our model by calculating the volume of the 28 June 2016 Lamplugh rock avalanche in Glacier Bay National Park, Alaska. We derived preevent and postevent 2‐m resolution DEMs from WorldView satellite stereo imagery. Using data from DEM differencing, we reconstructed the rock avalanche and adjacent surfaces at the time of occurrence by accounting for elevation changes due to ablation and scour of the ice surface, and postevent deposit changes. We accounted for uncertainties in our DEMs through precise coregistration and an assessment of relative elevation accuracy in bedrock control areas. The rock avalanche initially displaced 51.7 ± 1.5 Mm3 of intact rock and then scoured and entrained 13.2 ± 2.2 Mm3 of snow and ice during emplacement. We calculated the total deposit volume to be 69.9 ± 7.9 Mm3. Volume estimates that did not account for topographic changes due to ablation, scour, and compaction underestimated the deposit volume by 31.0–46.8 Mm3. Our model provides an improved framework for estimating uncertainties affecting rock avalanche volume measurements in glacial environments. These improvements can contribute to advances in the understanding of rock avalanche hazards and dynamics.

Alaska↗