USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ocean Modelling”

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 919 records · Page 51Linked to original sources
USGS2004

Suspended sediment and hydrodynamics above mildly sloped long wave ripples

We investigate the spatial and temporal distribution of suspended sediment and the associated hydrodynamics over mildly sloped long wave ripples on the inner shelf. These bedforms had wavelengths of approximately 1 m and heights of approximately 5 cm, in a mean water depth of 4 m. The vertical and temporal structures of the suspended sediment concentration (SSC) are consistent with the entrainment of sediment on the offshore flank of the ripple, and rapid vertical mixing at the time of flow reversal, followed by advection onshore by the onshore fluid motion. This work confirms that the mechanism for sediment suspension above low-amplitude, long wave ripples is similar to the vortex formation process expected over steeper vortex ripples. Numerical simulations of the flow using the Dune2d model indicate that a separated rotational flow structure is generated at the flank of the ripple on the seaward side of the ripple crest, near the time of flow reversal. The simulations indicate that only one vortex is formed during each wave period, in agreement with the field observations. This asymmetry is due mainly to the presence of an offshore mean near-bed current of approximately 6 to 8 cm/s. The SSC is calculated by the model and compared to the field observations. A hydraulic bed roughness of 10 to 15 median grain diameters ( d 50 ) was used in order to match the model prediction to the observed SSC approximately 1 cm above the seabed (cab). However, the modeled SSC and turbulent kinetic energy were significantly lower than the field observations at elevations exceeding approximately 2 cab.

North Carolina↗
USGS2012

Validation of a coupled wave-flow model in a high-energy setting: the mouth of the Columbia River

A monthlong time series of wave, current, salinity, and suspended-sediment measurements was made at five sites on a transect across the Mouth of Columbia River (MCR). These data were used to calibrate and evaluate the performance of a coupled hydrodynamic and wave model for the MCR based on the Delft3D modeling system. The MCR is a dynamic estuary inlet in which tidal currents, river discharge, and wave-driven currents are all important. Model tuning consisted primarily of spatial adjustments to bottom drag coefficients. In combination with (near-) default parameter settings, the MCR model application is able to simulate the dominant features in the tidal flow, salinity and wavefields observed in field measurements. The wave-orbital averaged method for representing the current velocity profile in the wave model is considered the most realistic for the MCR. The hydrodynamic model is particularly effective in reproducing the observed vertical residual and temporal variations in current structure. Density gradients introduce the observed and modeled reversal of the mean flow at the bed and augment mean and peak flow in the upper half of the water column. This implies that sediment transport during calmer summer conditions is controlled by density stratification and is likely net landward due to the reversal of flow near the bed. The correspondence between observed and modeled hydrodynamics makes this application a tool to investigate hydrodynamics and associated sediment transport.

Oregon, Washington↗
USGS2021

Velocity-porosity relations in carbonate and siliciclastic subduction zone input materials

The mechanical, physical, and frictional properties of incoming materials play an important role in subduction zone structure and slip behavior because these properties influence the strength of the accretionary wedge and megathrust plate boundary faults. Incoming sediment sections often show an increase in compressional wave speed (Vp) and a decrease in porosity with depth due to consolidation. These relations allow seismic-velocity models to be used to elucidate properties and conditions at depth. However, variations in these properties are controlled by lithology and composition as well as cementation and diagenesis. We present an analysis of shipboard measurements of Vp and porosity on incoming sediment cores from International Ocean Discovery Program (IODP) expeditions at the Hikurangi Margin, Nankai Trough, Aleutian Trench, Middle America Trench, and Sunda Trench. Porosity for these samples ranges from 5% to 85% and Vp ranges from 1.5 to 6 km/s. Vp-porosity relations developed by Erikson & Jarrad (1998), https://doi.org/10.1029/98JB02128 and Hoffman & Tobin (2004) https://10.2973/odp.proc.sr.190196.355.2004, with a critical porosity of ∼30%, can represent carbonate-poor (<50 wt% CaCO3), mainly hemipelagic, incoming sediment regardless of the margin. But these relations tend to underestimate porosity in incoming sediments with carbonate content greater than 50 wt%, which appear to have a critical porosity of between 45% and 50%. This discrepancy will lead to inaccuracy in estimates of fluid budget and overpressure in subduction zones. The velocity-porosity relation in carbonate sediments is non-unique due to the complexity that results from the greater susceptibility of carbonate rocks to diagenetic processes.

Geochemistry, Geophysics, Geosystems↗
USGS2012

Methane hydrates and global climate

Methane hydrates are globally widespread in permafrost regions and beneath the sea in sediment of outer continental margins. The amount of methane sequestered in gas hydrates is probably enormous, but estimates of the amounts are speculative and range over three orders of magnitude (about 10 3 to 10 6 GT (gigatons = 10 15 g)). A question of current interest concerns the possible consequences of an addition of methane to the atmosphere from destabilized methane hydrates due to global warming. Models of greenhouse warming predict that climatic change will be greatest in the Arctic. Thus, if methane from destabilized gas hydrates contributes to greenhouse warming, this destabilization will most likely take place first in the Arctic, particularly in the shallow nearshore regions of the Arctic Ocean where offshore permafrost is found. The process of permafrost warming and release of methane from gas hydrates may already be in progress, but the amount being released now and to be released in the 21st century is probably small. The positive feedback of this atmospheric methane on global climates will likely be minimal.

Global Biogeochemical Cycles↗
USGS2014

A previously unrecognized path of early Holocene base flow and elevated discharge from Lake Minong to Lake Chippewa across eastern Upper Michigan

It has long been hypothesized that flux of fresh meltwater from glacial Lake Minong in North America's Superior Basin to the North Atlantic Ocean triggered rapid climatic shifts during the early Holocene. The spatial context of recent support for this idea demands a reevaluation of the exit point of meltwater from the Superior Basin. We used ground penetrating radar (GPR), foundation borings from six highway bridges, a GIS model of surface topography, geologic maps, U.S. Department of Agriculture–Natural Resources Conservation Service soils maps, and well logs to investigate the possible linkage of Lake Minong with Lake Chippewa in the Lake Michigan Basin across eastern Upper Michigan. GPR suggests that a connecting channel lies buried beneath the present interlake divide at Danaher. A single optical age hints that the channel aggraded to 225 m as elevated receipt of Lake Agassiz meltwater in the Superior Basin began to wane <10.6 ka. The large supply of sediment required to accommodate aggradation was immediately available at the channel's edge in the littoral shelves of abandoned Lake Algonquin and in distal parts of post-Algonquin fans. As discharge decreased further, the aggraded channel floor was quickly breached and interbasin flow to Lake Chippewa was restored. Basal radiocarbon ages on wood from small lakes along the discharge path and a GIS model of Minong's shoreline are consistent with another transgression of Minong after ca. 9.5 ka. At the peak of the latter transgression, the southeastern rim of the Superior Basin (Nadoway Drift Barrier) failed, ending Lake Minong. Upon Minong's final drop, aggradational sediments were deposited at Danaher, infilling the prior breach.

Michigan↗
USGS2026

Development of projected depth-duration-frequency curves for precipitation in Florida, 2020–59 and 2050–89

The planning, permitting, and design of stormwater-management projects require estimates of the depths of extreme precipitation for current and future events with specified durations and return periods. In this project, precipitation data from six downscaled climate datasets were used to determine changes in precipitation depth-duration-frequency curves from the period 1966–2005 to the periods 2020–59 and 2050–89. The downscaled climate datasets are from the Coupled Model Intercomparison Project Phases 5 and 6 and include (1) Coordinated Regional Downscaling Experiment (CORDEX), (2) Localized Constructed Analogs (LOCA), (3) Multivariate Adaptive Constructed Analogs (MACA), (4) Jupiter Intelligence Weather Research and Forecasting model (JupiterWRF), (5) LOCA version 2 (LOCA2), and (6) National Aeronautics and Space Administration Earth Exchange Global Daily Downscaled Projections (NEX-GDDP). Change factors—multiplicative changes in expected extreme precipitation magnitude from a historical to future period—were computed for grid cells containing National Oceanic and Atmospheric Administration Atlas 14 stations in Florida. Change factors for specific durations and return periods were developed to scale the National Oceanic and Atmospheric Administration Atlas 14 historical depth-duration-frequency values to the periods 2020–59 and 2050–89 on the basis of changes in extreme precipitation derived from six downscaled climate datasets. Overall, a large variation in change factors across downscaled climate datasets was found, with change factors generally being greater than 1 and increasing with return period. In general, median change factors were found to range within 1.01–1.58 for 2020–59 and 1.01–1.63 for 2050–89, depending on the downscaled climate dataset, region, duration, and return period, indicating a projected overall increase in future extreme-precipitation events. When data from all datasets are considered together, median change factors range within 1.04–1.18 for the period 2020–59 and within 1.04–1.23 for the period 2050–89, depending on the region, duration, and return period. Spatial patterns in median change factors were found to vary by dataset.

Florida↗
USGS2020

Using video survey to examine the effect of habitat on gag grouper encounter

Gag is a reef fish that was declared overfished in the Gulf of Mexico (GOM) in 2009. Although Gag are no longer listed as overfished, fisheries managers are concerned that stocks may not be recovering. Our objective was to identify habitat characteristics important to Gag, and their effect on the probability of Gag occurrence. We obtained data from three separate fisheries-independent video surveys that sampled in the eastern GOM from 2010-2017: the National Atmospheric and Oceanic Administration (NOAA) Panama City, FL Office, the NOAA Southeast Area Monitoring and Assessment Program, and the Florida Fish and Wildlife Research Institute. We ran a separate mixed effects logistic regression for each survey, and used Akaike’s Information Criteria to determine the best fitting models. Some variables - percent rock coverage, vertical relief, latitude, and depth - were present in all confidence models. Depth did not have the same relationship with Gag across all surveys, suggesting that shallower habitats (<50 m) might be more suitable for juveniles, whereas deeper habitats (>50 m) might be more suitable for adults. Managers may be able to help Gag and encourage their recovery by using these data to establish or expand protected areas throughout shallower waters.

Florida↗
USGS2015

Seismic velocities within the sedimentary succession of the Canada Basin and southern Alpha-Mendeleev Ridge, Arctic Ocean: evidence for accelerated porosity reduction?

The Canada Basin and the southern Alpha-Mendeleev ridge complex underlie a significant proportion of the Arctic Ocean, but the geology of this undrilled and mostly ice-covered frontier is poorly known. New information is encoded in seismic wide-angle reflections and refractions recorded with expendable sonobuoys between 2007 and 2011. Velocity–depth samples within the sedimentary succession are extracted from published analyses for 142 of these records obtained at irregularly spaced stations across an area of 1.9E + 06 km 2 . The samples are modelled at regional, subregional and station-specific scales using an exponential function of inverse velocity versus depth with regionally representative parameters determined through numerical regression. With this approach, smooth, non-oscillatory velocity–depth profiles can be generated for any desired location in the study area, even where the measurement density is low. Practical application is demonstrated with a map of sedimentary thickness, derived from seismic reflection horizons interpreted in the time domain and depth converted using the velocity–depth profiles for each seismic trace. A thickness of 12–13 km is present beneath both the upper Mackenzie fan and the middle slope off of Alaska, but the sedimentary prism thins more gradually outboard of the latter region. Mapping of the observed-to-predicted velocities reveals coherent geospatial trends associated with five subregions: the Mackenzie fan; the continental slopes beyond the Mackenzie fan; the abyssal plain; the southwestern Canada Basin; and, the Alpha-Mendeleev magnetic domain. Comparison of the subregional velocity–depth models with published borehole data, and interpretation of the station-specific best-fitting model parameters, suggests that sandstone is not a predominant lithology in any of the five subregions. However, the bulk sand-to-shale ratio likely increases towards the Mackenzie fan, and the model for this subregion compares favourably with borehole data for Miocene turbidites in the eastern Gulf of Mexico. The station-specific results also indicate that Quaternary sediments coarsen towards the Beaufort-Mackenzie and Banks Island margins in a manner that is consistent with the variable history of Laurentide Ice Sheet advance documented for these margins. Lithological factors do not fully account for the elevated velocity–depth trends that are associated with the southwestern Canada Basin and the Alpha-Mendeleev magnetic domain. Accelerated porosity reduction due to elevated palaeo-heat flow is inferred for these regions, which may be related to the underlying crustal types or possibly volcanic intrusion of the sedimentary succession. Beyond exploring the variation of an important physical property in the Arctic Ocean basin, this study provides comparative reference for global studies of seismic velocity, burial history, sedimentary compaction, seismic inversion and overpressure prediction, particularly in mudrock-dominated successions.

Geophysical Journal International↗
USGS2017

A physical model for extreme drought over southwest Asia

The socioeconomic difficulties of southwest Asia, defined as the area bound by the domain 25°N–40°N and 40°E–70°E, are exacerbated by extreme precipitation deficits during the November–April rainy season. The precipitation deficits during many southwest Asia droughts have been examined in terms of the forcing by climate variability originating over the Pacific Ocean as a result of the El Niño–Southern Oscillation (ENSO), Pacific decadal variability (PDV), and the long-term warming of Pacific (LT) sea surface temperatures (SST). Here we examine how the most extreme November–April southwest Asia droughts relate to global SSTs and the associated large-scale atmospheric circulation anomalies and analyze the specific atmospheric forcing mechanisms responsible for changes in regional southwest Asian precipitation. The driest November–April seasons during 1948–2012 over southwest Asia are forced by subsidence and reductions of moisture fluxes as a result of the interaction of the mean flow with anomalous zonally symmetric high pressure throughout the Northern Hemisphere. The anomalous zonally symmetric high pressure throughout the Northern Hemisphere occurs simultaneously with cool central and eastern Pacific SST anomalies associated with La Niña and the negative phase of PDV and a warm west Pacific Ocean caused in part by the long-term warming of the west Pacific Ocean.

Book chapter↗
USGS2010

Poroelastic stress-triggering of the 2005 M8.7 Nias earthquake by the 2004 M9.2 Sumatra-Andaman earthquake

The M9.2 Sumatra-Andaman earthquake (SAE) occurred three months prior to the M8.7 Nias earthquake (NE). We propose that the NE was mechanically triggered by the SAE, and that poroelastic effects were a major component of this triggering. This study uses 3D finite element models (FEMs) of the Sumatra-Andaman subduction zone (SASZ) to predict the deformation, stress, and pore pressure fields of the SAE. The coseismic slip distribution for the SAE is calibrated to near-field GPS data using FEM-generated Green's Functions and linear inverse methods. The calibrated FEM is then used to predict the postseismic poroelastic contribution to stress-triggering along the rupture surface of the NE, which is adjacent to the southern margin of the SAE. The coseismic deformation of the SAE, combined with the rheologic configuration of the SASZ produces two transient fluid flow regimes having separate time constants. SAE coseismic pore pressures in the relatively shallow forearc and volcanic arc regions (within a few km depth) dissipate within one month after the SAE. However, pore pressures in the oceanic crust of the down-going slab persist several months after the SAE. Predictions suggest that the SAE initially induced MPa-scale negative pore pressure near the hypocenter of the NE. This pore pressure slowly recovered (increased) during the three-month interval separating the SAE and NE due to lateral migration of pore fluids, driven by coseismic pressure gradients, within the subducting oceanic crust. Because pore pressure is a fundamental component of Coulomb stress, the MPa-scale increase in pore pressure significantly decreased stability of the NE fault during the three-month interval after the SAE and prior to rupture of the NE. A complete analysis of stress-triggering due to the SAE must include a poroelastic component. Failure to include poroelastic mechanics will lead to an incomplete model that cannot account for the time interval between the SAE and NE. Our transient poroelastic model explains both the spatial and temporal characteristics of triggering of the NE by the SAE.

Earth and Planetary Science Letters↗
USGS2019

Lithosphere and shallow asthenosphere rheology from observations of post-earthquake relaxation

In tectonically active regions, post-earthquake motions are generally shaped by a combination of continued fault slippage (afterslip) on a timescale of days to months and viscoelastic relaxation of the lower crust and upper mantle on a timescale of days to years. Transient crustal motions have been observed following numerous magnitude >~7 earthquakes in various tectonic settings: continental rift zones (Basin and Range), continental plate boundary zones (San Andreas fault corridor; Alaska; Turkey), subduction zones (Japan, Chile, Sumatra), ongoing continental collision zones (Arabia; Tibet), and mid-ocean rifting zones (Iceland). When afterslip can be discriminated from viscoelastic relaxation and when temporal coverage of the postseismic measurements is broad (i.e., geodetic surveys of at least several years duration are available), a wide spectrum of relaxation timescales are usually identified. Current temporal resolution and modeling approaches (e.g., Burgers body analog) allow identification of transient (Kelvin) and steady-state (Maxwell) viscosities that are operable in the short-term and long-term, respectively. I compile results from 40 studies of post-earthquake motions, augmented by ten studies of contemporary surface loading or unloading, that illuminate current estimates of transient and steady-state viscosity of the lower crust and/or uppermost mantle. Lower crust viscosity estimates range from ~ 10 18 to 10 21 Pa s, with most estimates near the upper end except in areas of overthickened crust. Mantle lithosphere and asthenosphere viscosity estimates are particularly abundant and yield a picture of transient viscosity ranging from ~10 16 to 10 19 Pa s and steady-state viscosity ranging from ~10 18 to 10 21 Pa s. To first order, both transient and steady-state viscosities are well correlated with regional heat flow, and steady-state viscosities are comparable with temperature and strain-rate dependent rock viscosities from laboratory-based flow laws.

Physics of Earth and Planetary Interiors↗
USGS2021

Twenty-first-century projections of shoreline change along inlet-interrupted coastlines

Sandy coastlines adjacent to tidal inlets are highly dynamic and widespread landforms, where large changes are expected due to climatic and anthropogenic influences. To adequately assess these important changes, both oceanic (e.g., sea-level rise) and terrestrial (e.g., fluvial sediment supply) processes that govern the local sediment budget must be considered. Here, we present novel projections of shoreline change adjacent to 41 tidal inlets around the world, using a probabilistic, reduced complexity, system-based model that considers catchment-estuary-coastal systems in a holistic way. Under the RCP 8.5 scenario, retreat dominates (90% of cases) over the twenty-first century, with projections exceeding 100 m of retreat in two-thirds of cases. However, the remaining systems are projected to accrete under the same scenario, reflecting fluvial influence. This diverse range of response compared to earlier methods implies that erosion hazards at inlet-interrupted coasts have been inadequately characterised to date. The methods used here need to be applied widely to support evidence-based coastal adaptation.

Nature--Scientific Reports↗
USGS1993

The global carbon dioxide budget

The increase in atmospheric CO 2 levels during the last deglaciation was comparable in magnitude to the recent historical increase. However, global CO 2 budgets for these changes reflect fundamental differences in rates and in sources and sinks. The modern oceans are a rapid net CO 2 sink, whereas the oceans were a gradual source during the deglaciation. Unidentified terrestrial CO 2 sinks are important uncertainties in both the deglacial and recent CO 2 budgets. The deglacial CO 2 budget represents a complexity of long-term dynamic behavior that is not adequately addressed by current models used to forecast future atmospheric CO 2 levels.

Science↗
USGS1991

Volcano spacing and plate rigidity

In-plane stresses, which accompany the flexural deformation of the lithosphere under the load of adjacent volcanoes, may govern the spacing of volcanoes in hotspot provinces. Specifically, compressive stresses in the vicinity of a volcano prevent new upwelling in this area, forcing a new volcano to develop at a minimum distance that is equal to the distance in which the radial stresses change from compressional to tensile (the inflection point). If a volcano is modeled as a point load on a thin elastic plate, then the distance to the inflection point is proportional to the thickness of the plate to the power of 3/4. Compilation of volcano spacing in seven volcanic groups in East Africa and seven volcanic groups of oceanic hotspots shows significant correlation with the elastic thickness of the plate and matches the calculated distance to the inflection point. In contrast, volcano spacing in island arcs and over subduction zones is fairly uniform and is much larger than predicted by the distance to the inflection point, reflecting differences in the geometry of the source and the upwelling areas.

Geology↗
USGS1999

Aeromagnetic legacy of early Paleozoic subduction along the Pacific margin of Gondwana

Comparison of the aeromagnetic signatures and geology of southeastern Australia and northern Victoria Land, Antarctica, with similar data from ancient subduction zones in California and Japan, provides a framework for reinterpretation of the plate tectonic setting of the Pacific margin of early Paleozoic Gondwana. In our model, the plutons in the Glenelg (south-eastern Australia) and Wilson (northern Victoria Land) zones formed the roots of continental-margin magmatic arcs. Eastward shifting of arc magmatism resulted in the Stavely (south-eastern Australia) and Bowers (northern Victoria Land) volcanic eruptions onto oceanic forearc crust. The turbidites in the Stawell (southeastern Australia) and Robertson Bay (northern Victoria Land zones) shed from the Glenelg and Wilson zones, respectively, were deposited along the trench and onto the subducting oceanic plate. The margin was subsequently truncated by thrust faults and uplifted during the Delamerian and Ross orogenies, leading to the present-day aeromagnetic signatures.

Geology↗
USGS2010

High-resolution seismic-reflection data offshore of Dana Point, southern California borderland

The U.S. Geological Survey collected high-resolution shallow seismic-reflection profiles in September 2006 in the offshore area between Dana Point and San Mateo Point in southern Orange and northern San Diego Counties, California. Reflection profiles were located to image folds and reverse faults associated with the San Mateo fault zone and high-angle strike-slip faults near the shelf break (the Newport-Inglewood fault zone) and at the base of the slope. Interpretations of these data were used to update the USGS Quaternary fault database and in shaking hazard models for the State of California developed by the Working Group for California Earthquake Probabilities. This cruise was funded by the U.S. Geological Survey Coastal and Marine Catastrophic Hazards project. Seismic-reflection data were acquired aboard the R/V Sea Explorer, which is operated by the Ocean Institute at Dana Point. A SIG ELC820 minisparker seismic source and a SIG single-channel streamer were used. More than 420 km of seismic-reflection data were collected. This report includes maps of the seismic-survey sections, linked to Google Earth? software, and digital data files showing images of each transect in SEG-Y, JPEG, and TIFF formats.

Open-File Report↗
USGS2020

Ideas and perspectives: A strategic assessment of methane and nitrous oxide measurements in the marine environment

>In the current era of rapid climate change, accurate characterization of climate-relevant gas dynamics – namely production, consumption, and net emissions – is required for all biomes, especially those ecosystems most susceptible to the impact of change. Marine environments include regions that act as net sources or sinks for numerous climate-active trace gases including methane (CH4) and nitrous oxide (N2O). The temporal and spatial distributions of CH4 and N2O are controlled by the interaction of complex biogeochemical and physical processes. To evaluate and quantify how these mechanisms affect marine CH4 and N2O cycling requires a combination of traditional scientific disciplines including oceanography, microbiology, and numerical modeling. Fundamental to these efforts is ensuring that the datasets produced by independent scientists are comparable and interoperable. Equally critical is transparent communication within the research community about the technical improvements required to increase our collective understanding of marine CH4 and N2O. A workshop sponsored by Ocean Carbon and Biogeochemistry (OCB) was organized to enhance dialogue and collaborations pertaining to marine CH4 and N2O. Here, we summarize the outcomes from the workshop to describe the challenges and opportunities for near-future CH4 and N2O research in the marine environment.

Biogeosciences↗
USGS2020

Satellite observations of the 2016–2017 eruption of Bogoslof volcano: Aviation and ash fallout hazard implications from a water-rich eruption

Satellite data (GOES, VIIRS, MODIS, AVHRR) were used in near real-time during the 2016–2017 eruption of Bogoslof volcano, Alaska, to detect explosive eruptive activity and to characterize the resulting volcanic clouds. This study examines satellite data to estimate volcanic cloud heights and mass eruption rates for 47 of the 70 explosive events. Eighteen of the volcanic clouds reached an altitude in excess of 8.5 km asl, where they posed a potential hazard to aviation. We estimate eruption rates were in the range of 10 4 –10 7 kg/s. Eruption durations were available for 28 events, and the total mass of these events was 5.7 × 10 10 kg. Most of the explosions occurred from submarine vents, producing volcanic clouds with water-rich characteristics in satellite data. We infer that these water-rich (phreatomagmatic) events contained ice-coated ash particles, which changed their visible and multispectral characteristics. Only two of the explosions produced clouds with satellite characteristics that would be considered ash-rich. We conclude that these events were relatively dry eruptions with limited access to ocean water. Although some of the explosive events transitioned from submarine to subaerial vents, we observed no change in the character of the volcanic clouds during these transitions. We speculate that enhanced removal of fine-grained volcanic ash likely occurred due to aggregation, with implications for modeling ash transport and fallout. We conclude that the majority of Bogoslof’s water-rich volcanic clouds did contain volcanic ash despite the lack of conventional “ash signature” in satellite data. This has implications for satellite monitoring of future water-rich/shallow submarine eruptions.

Alaska↗