Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Hazardous Materials”

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 55 records · Page 3Linked to original sources

Pore fluid pressure, apparent friction, and Coulomb failure

Many recent studies of stress-triggered seismicity rely on a fault failure model with a single free parameter, the apparent coefficient of friction, presumed to be a material constant with possible values 0 ≤ μ′ ≤ 1. These studies may present a misleading view of fault strength and the role of pore fluid pressure in earthquake failure. The parameter μ′ is intended to incorporate the effects of both friction and pore pressure, but is a material constant only if changes in pore fluid pressure induced by changes in stress are proportional to the normal stress change across the potential failure plane. Although specific models of fault zones permit such a relation, neither is it known that fault zones within the Earth behave this way, nor is this behavior expected in all cases. In contrast, for an isotropic homogeneous poroelastic model the pore pressure changes are proportional to changes in mean stress, μ′ is not a material constant, and −∞ ≤ μ′ ≤ +∞. Analysis of the change in Coulomb failure stress for tectonically loaded reverse and strike-slip faults shows considerable differences between these two pore pressure models, suggesting that such models might be distinguished from one another using observations of triggered seismicity (e.g., aftershocks). We conclude that using the constant apparent friction model exclusively in studies of Coulomb failure stress is unwise and could lead to significant errors in estimated stress change and seismic hazard.

Journal of Geophysical Research B: Solid Earth↗

Implications of silicic vent patterns for the presence of large crustal magma chambers

On the basis of the distribution of silicic vents, many volcanic fields can be grouped with (1) igneous systems that may be small and whose vent locations are controlled by regional tectonics, (2) those that include sizable crustal magma bodies which erupt at sites determined by their anomalous local stress fields, or (3) relatively small volume systems that are transitional between categories 1 and 2. Linear vent patterns that are aligned normal to the regional least principal stress (σ 3 ) commonly are associated with absence of evidence for large, shallow magma bodies. The Coso volcanic field and the Inyo‐Mono domes in California and probably the South Sister area in Oregon are examples of such systems. The 1960 dacitic fissure eruption at Cordón Caulle in southern Chile evidently is linked to tectonic stress relaxation associated with the great earthquake that occurred 48 hours before the eruption began. Large, shallow magma chambers are thought to perturb the local stress field so that areal patterns of silicic vents are diffuse, radial, or arcuate. Such systems may erupt great volumes of pyroclastic material catastrophically and produce large calderas. Well‐preserved examples of late precaldera leaks of silicic magma occur at Long Valley, California, and Mount Mazama (Crater Lake), Oregon. Other possible examples are noted. Some of these, those which formed calderas smaller than Crater Lake, apparently were preceded by silicic eruptions from aligned vents. Bearing in mind that there exists a transitional group between tectonically controlled small systems and very large magma chambers, vent distributions can be useful in evaluating potential volcanic hazards for silicic volcanic fields.

California↗

Gravitational stability of three-dimensional stratovolcano edifices

Catastrophic flank collapses have occurred at many stratovolcanoes worldwide. We present a three-dimensional (3-D) slope stability analysis for assessing and quantifying both the locations of minimum edifice stability and the expected volumes of potential failure. Our approach can search the materials underlying a topographic surface, represented as a digital elevation model (DEM), and determine the relative stability of all parts of the edifice. Our 3-D extension of Bishop's [1955] simplified limit-equilibrium analysis incorporates spherical failure surfaces, variable material properties, pore fluid pressures, and earthquake shaking. Although a variety of processes can trigger collapse, we focus here on gravitationally induced instability. Even homogeneous rock properties strongly influence the depth and volume of the least stable potential failure. For large failures in complex topography, patterns of potential instability do not mimic local ground surface slope alone. The May 18, 1980, catastrophic failure of the north flank of Mount St. Helens provides the best documented case history to test our method. Using the undeformed edifice topography of Mount St. Helens in an analysis of dry, static slope stability with homogeneous materials, as might be conducted in a precollapse hazard analysis, our method identified the northwest flank as the least stable region, although the north flank stability was within 5% of the minimum. Using estimates of the conditions that existed 2 days prior to collapse, including deformed topography with a north flank bulge and combined pore pressure and earthquake shaking effects, we obtained good estimates of the actual failure location and volume. Our method can provide estimates of initial failure volume and location to aid in assessing downslope or downstream hazards.

Washington↗

A statistical analysis of the global historical volcanic fatalities record

A new database of volcanic fatalities is presented and analysed, covering the period 1600 to 2010 AD. Data are from four sources: the Smithsonian Institution, Witham (2005), CRED EM-DAT and Munich RE. The data were combined and formatted, with a weighted average fatality figure used where more than one source reports an event; the former two databases were weighted twice as strongly as the latter two. More fatal incidents are contained within our database than similar previous works; approximately 46% of the fatal incidents are listed in only one of the four sources, and fewer than 10% are in all four. 278,880 fatalities are recorded in the database, resultant from 533 fatal incidents. The fatality count is dominated by a handful of disasters, though the majority of fatal incidents have caused fewer than ten fatalities. Number and empirical probability of fatalities are broadly correlated with VEI, but are more strongly influenced by population density around volcanoes and the occurrence and extent of lahars (mudflows) and pyroclastic density currents, which have caused 50% of fatalities. Indonesia, the Philippines, and the West Indies dominate the spatial distribution of fatalities, and there is some negative correlation between regional development and number of fatalities. With the largest disasters removed, over 90% of fatalities occurred between 5 km and 30 km from volcanoes, though the most devastating eruptions impacted far beyond these distances. A new measure, the Volcano Fatality Index, is defined to explore temporal changes in societal vulnerability to volcanic hazards. The measure incorporates population growth and recording improvements with the fatality data, and shows prima facie evidence that vulnerability to volcanic hazards has fallen during the last two centuries. Results and interpretations are limited in scope by the underlying fatalities data, which are affected by under-recording, uncertainty, and bias. Attempts have been made to estimate the extent of these issues, and to remove their effects where possible. The data analysed here are provided as supplementary material. An updated version of the Smithsonian fatality database fully integrated with this database will be publicly available in the near future and subsequently incorporate new data.

Journal of Applied Volcanology↗

Modeled coastal-ocean pathways of land-sourced contaminants in the aftermath of Hurricane Florence

Extreme precipitation during Hurricane Florence, which made landfall in North Carolina in September 2018, led to breaches of hog waste lagoons, coal ash pits, and wastewater facilities. In the weeks following the storm, freshwater discharge carried pollutants, sediment, organic matter, and debris to the coastal ocean, contributing to beach closures, algae blooms, hypoxia, and other ecosystem impacts. Here, the ocean pathways of land-sourced contaminants following Hurricane Florence are investigated using the Regional Ocean Modeling System (ROMS) with a river point source with fixed water properties from a hydrologic model (WRF-Hydro) of the Cape Fear River Basin, North Carolina's largest watershed. Patterns of contaminant transport in the coastal ocean are quantified with a finite duration tracer release based on observed flooding of agricultural and industrial facilities. A suite of synthetic events also was simulated to investigate the sensitivity of the river plume transport pathways to river discharge and wind direction. The simulated Hurricane Florence discharge event led to westward (downcoast) transport of contaminants in a coastal current, along with intermittent storage and release of material in an offshore (bulge) or eastward (upcoast) region near the river mouth, modulated by alternating upwelling and downwelling winds. The river plume patterns led to a delayed onset and long duration of contaminants affecting beaches 100 km to the west, days to weeks after the storm. Maps of the onset and duration of hypothetical water quality hazards for a range of weather conditions may provide guidance to managers on the timing of swimming/shellfishing advisories and water quality sampling.

Journal of Geophysical Research: Oceans↗

Liquefaction potential index: Field assessment

Cone penetration test (CPT) soundings at historic liquefaction sites in California were used to evaluate the predictive capability of the liquefaction potential index (LPI), which was defined by Iwasaki et al. in 1978. LPI combines depth, thickness, and factor of safety of liquefiable material inferred from a CPT sounding into a single parameter. LPI data from the Monterey Bay region indicate that the probability of surface manifestations of liquefaction is 58 and 93%, respectively, when LPI equals or exceeds 5 and 15. LPI values also generally correlate with surface effects of liquefaction: Decreasing from a median of 12 for soundings in lateral spreads to 0 for soundings where no surface effects were reported. The index is particularly promising for probabilistic liquefaction hazard mapping where it may be a useful parameter for characterizing the liquefaction potential of geologic units.

Journal of Geotechnical and Geoenvironmental Engin↗

Monogenetic origin of Ubehebe Crater maar volcano, Death Valley, California: Paleomagnetic and stratigraphic evidence

Paleomagnetic data for samples collected from outcrops of basaltic spatter at the Ubehebe Crater cluster, Death Valley National Park, California, record a single direction of remanent magnetization indicating that these materials were emplaced during a short duration, monogenetic eruption sequence ~ 2100 years ago. This conclusion is supported by geochemical data encompassing a narrow range of oxide variation, by detailed stratigraphic studies of conformable phreatomagmatic tephra deposits showing no evidence of erosion between layers, by draping of sharp rimmed craters by later tephra falls, and by oxidation of later tephra layers by the remaining heat of earlier spatter. This model is also supported through a reinterpretation and recalculation of the published age results (Sasnett et al., 2012) from an innovative and bold exposure-age study on very young materials. Their conclusion of multiple and protracted eruptions at Ubehebe Crater cluster is here modified through the understanding that some of their quartz-bearing clasts inherited from previous exposure on the fan surface (too old), and that other clasts were only exposed at the surface by wind and/or water erosion centuries after their eruption (too young). Ubehebe Crater cluster is a well preserved example of young monogenetic maar type volcanism protected within a National Park, and it represents neither a protracted eruption sequence as previously thought, nor a continuing volcanic hazard near its location.

California↗

Influence of grain size and shape on volcanic ash electrical conductivity

Few studies have examined the electrical properties of volcanic ash or considered the effects of physical characteristics, such as grain size and shape on its electrification. This study measures the resistivity of eight volcanic ash samples, three milled-samples and five natural ashfall samples from Alaska, U.S.A., using a current amplifier and examines the influence of particle size and particle shape on calculated grain conductivity. Grain conductivities are calculated by applying spherical, oblate, and prolate ellipsoidal shape parameters within a general effective media equation. Volcanic ash is an electrically insulating material with average resistivity values ranging from 1.732E 7 to 1.606E 12 Ωm. Results show that grain size has a limited impact on resistivity. Calculated grain conductivities remained within the same order of magnitude between the three different shape calculations ranging between 1.367E − 7 and 2.196E − 12 S/m using the spherical shape parameter. Electrical measurements of heterogeneous powders are complex and challenging. This study aims to expand the community's dataset and understanding of the electrical properties of volcanic ash and the impact of different physical characteristics. This has implication for understanding ash hazards to electrical infrastructure and volcanic lightning.

Journal of Volcanology and Geothermal Research↗

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↗

Tectonic relations of carbon dioxide discharges and earthquakes

CO 2 ‐rich springs occur worldwide along major zones of seismicity. They are mostly in young orogenic belts, but some are in areas of rifting continental platforms. Analyses of 13 C content indicate that much of the CO 2 is derived from the mantle and that other important sources are the metamorphism of marine carbonate‐bearing sedimentary rocks and the degradation of organic material. The presence of calc‐silicate minerals, such as pumpellyite in metagraywacke, is evidence of former conversion of carbonate‐bearing rocks into calc‐silicate minerals and release of CO 2 . The CO 2 pressure in fractured rocks of a fault zone reduces the effective normal stress and, if it is sufficiently great, allows the fault to slip. If the pressure were maintained at a sufficiently high level, the fault behavior might be characterized by frequent small earthquakes and aseismic creep such as occur along active segments of the San Andreas system. The presence of CO 2 ‐rich springs may indicate a potentially hazardous seismic region. Monitoring of CO 2 discharges could be useful in earthquake prediction.

Journal of Geophysical Research B: Solid Earth↗

Geology of El Chichon volcano, Chiapas, Mexico

The (pre-1982) 850-m-high andesitic stratovolcano El Chicho??n, active during Pleistocene and Holocene time, is located in rugged, densely forested terrain in northcentral Chiapas, Me??xico. The nearest neighboring Holocene volcanoes are 275 km and 200 km to the southeast and northwest, respectively. El Chicho??n is built on Tertiary siltstone and sandstone, underlain by Cretaceous dolomitic limestone; a 4-km-deep bore hole near the east base of the volcano penetrated this limestone and continued 770 m into a sequence of Jurassic or Cretaceous evaporitic anhydrite and halite. The basement rocks are folded into generally northwest-trending anticlines and synclines. El Chicho??n is built over a small dome-like structure superposed on a syncline, and this structure may reflect cumulative deformation related to growth of a crustal magma reservoir beneath the volcano. The cone of El Chicho??n consists almost entirely of pyroclastic rocks. The pre-1982 cone is marked by a 1200-m-diameter (explosion?) crater on the southwest flank and a 1600-m-diameter crater apparently of similar origin at the summit, a lava dome partly fills each crater. The timing of cone and dome growth is poorly known. Field evidence indicates that the flank dome is older than the summit dome, and K-Ar ages from samples high on the cone suggest that the flank dome is older than about 276,000 years. At least three pyroclastic eruptions have occurred during the past 1250 radiocarbon years. Nearly all of the pyroclastic and dome rocks are moderately to highly porphyritic andesite, with plagioclase, hornblende and clinopyroxene the most common phenocrysts. Geologists who mapped El Chicho??n in 1980 and 1981 warned that the volcano posed a substantial hazard to the surrounding region. This warning was proven to be prophetic by violent eruptions that occurred in March and April of 1982. These eruptions blasted away nearly all of the summit dome, blanketed the surrounding region with tephra, and sent pyroclastic flows down radial drainages on the flanks of the cone; about 0.3 km3 of material (density of all products normalized to 2.6 g cm-3) was erupted. More debris entered the stratosphere than from any other volcanic eruption within at least the past two decades. Halite and a calcium sulfate mineral (anhydrite?) recovered from the stratospheric cloud, and anhydrite as a common accessory mineral in 1982 juvenile erupted products may reflect contamination of El Chicho??n magma by the evaporite sequence revealed by drilling. ?? 1984.

Journal of Volcanology and Geothermal Research↗

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska↗

Selection of the Mars Pathfinder landing site

The Mars Pathfinder spacecraft will land on a depositional fan near the mouth of the catastrophic outflow channel, Ares Vallis (19.5°N, 32.8°W). This site offers the prospect of analyzing a variety of rock types from the ancient cratered highlands, intermediate-age ridged plains, and reworked channel deposits. Analyses of these rocks by Pathfinder instruments will enable first-order scientific questions to be addressed, such as differentiation of the crust, the development of weathering products, and the nature of the early environment, as well as their subsequent evolution on Mars. Constraints imposed by (1) spacecraft and rover designs (which are robust), (2) entry, descent, and landing, (3) scientific potential at various sites, and (4) safety were important considerations in site selection. Engineering constraints require a 70 km by 200 km smooth, flat (low slopes) area located between 10° and 20°N that is below 0 km elevation, with average radar reflectivity, little dust, and moderate rock abundance. Three regions on Mars are between 10° and 20°N and below 0 km elevation: Chryse, Amazonis, and Isidis-Elysium. Science considerations favor sites at the mouths of outflow channels (grab bag sites offer an assay of rock types on Mars), highland sites (early crustal differentiation and climate), and sites covered with dark (unoxidized) material. Sites are considered safe if they are clearly below 0 km elevation, appear acceptably free of hazards in high-resolution (<50 m/pixel) Viking orbiter images and have acceptable reflectivity and roughness at radar wavelengths, thermal inertia, rock abundance, red to violet ratio, and albedo. Recent 3.5-cm wavelength radar observations were used to verify elevation, reflectivity, and roughness within the landing ellipses. Three sites meet all of these criteria: Ares Vallis, Tritonis Lacus, and Isidis. Although Isidis appears to be safer than Tritonis and Ares, the greater scientific potential at Ares Vallis resulted in its selection. Comparisons of the Grand Coulee (channel) and the depositional Ephrata Fan of the Channeled Scabland in eastern Washington, with Ares Vallis and its depositional fan also suggest the Ares Vallis landing site is safe and scientifically interesting.

Journal of Geophysical Research Planets↗

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗

Multi-model comparison of computed debris flow runout for the 9 January 2018 Montecito, California post-wildfire event

Hazard assessment for post-wildfire debris flows, which are common in the steep terrain of the western United States, has focused on the susceptibility of upstream basins to generate debris flows. However, reducing public exposure to this hazard also requires an assessment of hazards in downstream areas that might be inundated during debris flow runout. Debris flow runout models are widely available, but their application to hazard assessment for post-wildfire debris flows has not been extensively tested. Necessary inputs to these models include the total volume of the mobilized flow, flow properties (either inherent material properties or calibration coefficients), and site topography. Estimates of volume are possible in post-event (“back calculation”) studies, yet before an event, volume is an uncertain quantity. We simulated debris flow runout for the well-constrained 9 January 2018 Montecito event using three models (RAMMS, FLO2D, and D-Claw) to determine the relative importance of volume and flow properties. We broke the impacted area into three domains, and for each model-domain combination, we performed a numerical sampling study in which volume and flow properties varied within a wide, but plausible range. We assessed model performance based on inundation patterns and peak flow depths. We found all models could simulate the event with comparable results. Simulation performance was most sensitive to flow volume and less sensitive to flow properties. Our results emphasize the importance of reducing uncertainty in pre-event estimates of flow volume for hazard assessment.

California↗

Insecticide contaminations in wetland habitats and their effects on fish- eating birds

An unusual mortality of fish-eating birds occurred at the Tule Lake National Wildlife Refuge in California between 1960 and 1962. Over 1100 dead birds of ten species were found during that period. Investigations of the mortality indicated that birds died as a result of their exposure to toxaphene, which was applied to adjoining agricultural areas, transported to the refuge in waste irrigation water, and accumulated in fish eaten by the birds. Toxaphene was used in agriculture for only 3 years, 1958-60, and by 1962 residues were no longer found in the marsh habitat. However, results also showed that a potentially hazardous contamination of DDT and its metabolites was present in the birds and the environment. Studies were, therefore, continued to determine the food-chain relationships of DDT in marsh ecosystems, and the effects of the insecticide on birds. Sampling for insecticide residues was undertaken at four waterfowl refuges, each representing a somewhat different set of environmental conditions. DDT residues in water at each refuge were found to be relatively low, but filtered samples showed that residues in suspended material were 10 000 to 20 000 times as great as those in the filtrate. Much of the suspended material was organic matter, which is an important source of energy for the ecosystems. In environments with invertebrates, DDT residues were relatively low in sediments and high in fish, while residues in vegetation were highest where invertebrates were absent. Results suggest that organic materials and invertebrates can be important factors in determining the fate and involvement of insecticides in aquatic environments. Studies of the movements, exposure to pesticides, productivity, and longevity of individuals in three nesting colonies of White Pelicans were begun in 1962 to depict the relationships of pesticides to a migratory population of birds dependent upon wetland habitats. To date there is no evidence that the exposure of pelicans to DDT is influencing the population dynamics of the birds.

California, Oregon↗

Leachability of uranium and other elements from freshly erupted volcanic ash

A study of leaching of freshly erupted basaltic and dacitic air-fall ash and bomb fragment samples, unaffected by rain, shows that glass dissolution is the dominant process by which uranium is initially mobilized from air-fall volcanic ash. Si, Li, and V are also preferentially mobilized by glass dissolution. Gaseous transfer followed by fixation of soluble uranium species on volcanic-ash particles is not an important process affecting uranium mobility. Gaseous transfer, however, may be important in forming water-soluble phases, adsorbed to ash surfaces, enriched in the economically and environmentally important elements Zn, Cu, Cd, Pb, B, F, and Ba. Quick removal of these adsorbed elements by the first exposure of freshly erupted ash to rain and surface water may pose short-term hazards to certain forms of aquatic and terrestrial life. Such rapid release of material may also represent the first step in transportation of economically important elements to environments favorable for precipitation into deposits of commercial interest. Ash samples collected from the active Guatemalan volcanoes Fuego and Pacaya (high-Al basalts) and Santiaguito (hornblende-hypersthene dacite); bomb fragments from Augustine volcano (andesite-dacite), Alaska, and Heimaey (basalt), Vestmann Islands, Iceland; and fragments of "rhyolitic" pumice from various historic eruptions were subjected to three successive leaches with a constant water-to-ash weight ratio of 4:1. The volcanic material was successively leached by: (1) distilled-deionized water (pH = 5.0-5.5) at room temperature for 24 h, which removes water-soluble gases and salts adsorbed on ash surfaces during eruption; (2) dilute HCl solution (pH = 3.5-4.0) at room temperature for 24 h, which continues the attack initiated by the water and also attacks acid-soluble sulfides and oxides; (3) a solution 0.05 M in both Na,CO, and NaHCO, (pH = 9.9) at 80&deg;C for one week, which preferentially dissolves volcanic glass. The first two leaches mimic interaction of ash with rain produced in the vicinity of an active eruption. The third leach accelerates the effect of prolonged contact of volcanic ash with alkaline ground water present during ash diagenesis.

Journal of Volcanology and Geothermal Research↗

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↗