Earth Sciences Commons^™

Brief Communication: Is Vertical Shear In An Ice Shelf (Still) Negligible?, Chris Miele, Timothy C. Bartholomaus, Ellyn M. Enderlin

Geosciences Faculty Publications and Presentations

Vertical shear is recognized today as a key component of the stress balance of ice shelves. However, the first ice shelf models were built on the neglect of vertical shear. Partly due to its historical treatment, it remains common to discuss vertical shear as though it were still considered negligible in ice shelf models. Here, we offer a historical perspective on the changing treatment of vertical shear over time, and we emphasize the term's non-negligibility in current ice shelf modeling. We illustrate our discussion in the simplest context of an analytic, isothermal, shallow-ice-shelf model.

Hurricane Impacts On Land In The Central And Eastern Caribbean Since 1494 Ce From Written Records, Michael Chenoweth, Ian Howard

Geosciences Faculty Publications and Presentations

Written accounts of hurricanes and their impacts in the Caribbean region date back to 1494 CE We report a new compilation of hurricanes in the longest settled regions of the central and eastern Caribbean region (CECR) that is the most complete yet produced and present basic statistics. We assess likely undercounts of hurricanes due to incomplete documentary records in the earliest part of the record by using the most complete reporting from Puerto Rico relative to the remainder of the CECR. We compare our results with other documentary and proxy data from the region and demonstrate using wavelet analysis significant …

Forecasting Sea Level Rise-Driven Inundation In Diked And Tidally Restricted Coastal Lowlands, K. M. Befus, A. P.D. Kurnizki, K. D. Kroeger, M. J. Eagle, T. P. Smith

Geosciences Faculty Publications and Presentations

Diked and drained coastal lowlands rely on hydraulic and protective infrastructure that may not function as designed in areas with relative sea-level rise. The slow and incremental loss of the hydraulic conditions required for a well-drained system make it difficult to identify if and when the flow structures no longer discharge enough water, especially in tidal settings where two-way flows occur through the dike. We developed and applied a hydraulic mass-balance model to quantify how water levels in the diked and tidally restricted coastal wetlands and water bodies dynamically respond to sea-level rise, specifically applied to the Herring River Estuary …

Early Permian Zircon Ages From The P. Confluens And P. Pseudoreticulata Spore-Pollen Zones In The Southern Bonaparte And Canning Basins, Northwestern Australia, A. J. Mory, J. Crowley, J. Backhouse, R. S. Nicoll, J. D. Gorter

Geosciences Faculty Publications and Presentations

The Pseudoreticulatispora confluens–P. pseudoreticulata spore-pollen zonal datum typically coincides with the end of widespread Permian glacial deposits in Western Australia. Although previously attributed to the mid-Sakmarian, chemical abrasion isotope dilution thermal ionisation mass spectrometry (TIMS) dating of zircons from volcanic tuffs in the Ditji Formation of the Bonaparte Basin and the Grant Group in the Canning Basin point to an Asselian age of about 295.25 Ma for this datum. All dated zircons from the Ditji Formation came from petroleum well cuttings but the accompanying palynology was mostly from sidewall cores; however, all Grant Group samples were from conventional core. TIMS …

Apatite Reference Materials For Sims Microanalysis Of Isotopes And Trace Elements, James L. Crowley, Mark Schmitz

Geosciences Faculty Publications and Presentations

Twelve apatite samples have been tested as secondary ion mass spectrometry (SIMS) reference materials. Laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) analysis shows that the SLAP, NUAN and GR40 apatite gems are internally homogeneous, with most trace element mass fractions having 2 standard deviations (2s) ≤ 2.0%. BR2, BR5, OL2, AFG2 and AFB1, which have U > 63 μg g^-1, ²⁰⁶Pb/²⁰⁴Pb > 283, and homogeneous SIMS U-Pb data, have respective isotope dilution thermal ionisation mass spectrometry (ID-TIMS) ages of 2053.83 ± 0.21 Ma, 2040.34 ± 0.09 Ma, 868.87 ± 0.25 Ma, 478.71 ± 0.22 Ma and …

Diverged Landscape Of Restaurant Recovery From The Covid-19 Pandemic In The United States, Siqin Wang, Xiao Huang, Bing She, Zhenlong Li

Geosciences Faculty Publications and Presentations

The COVID-19 pandemic has imposed catastrophic impacts on the restaurant industry as a crucial socioeconomic sector that contributes to the global economy. However, the understanding of how the restaurant industry was recovered from COVID-19 remains underexplored. This study constructs a spatially explicit evaluation of the effect of COVID-19 on the restaurant industry in the US, drawing on the attributes of +200,000 restaurants from Yelp and +600 million individual-level restaurant visitations provided by SafeGraph from 1st January 2019 to 31st December 2021. We produce quantitative evidence of lost restaurant visitations and revenue amid the pandemic, the changes in the customers’ origins, …

Estimating Snow Accumulation And Ablation With L-Band Interferometric Synthetic Aperture Radar (Insar), Jack Tarricone, Ryan W. Webb, Hans-Peter Marshall, Anne W. Nolin, Franz J. Meyer

Geosciences Faculty Publications and Presentations

Snow is a critical water resource for the western United States and many regions across the globe. However, our ability to accurately measure and monitor changes in snow mass from satellite remote sensing, specifically its water equivalent, remains a challenge. To confront these challenges, NASA initiated the SnowEx program, a multiyear effort to address knowledge gaps in snow remote sensing. During SnowEx 2020, the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) team acquired an L-band interferometric synthetic aperture radar (InSAR) data time series to evaluate the capabilities and limitations of repeat-pass L-band InSAR for tracking changes in snow water equivalent …

Volcano Opto-Acoustics: Mapping The Infrasound Wavefield At Yasur Volcano (Vanuatu), J. B. Johnson, T. Boyer, L. M. Watson, J. F. Anderson

Geosciences Faculty Publications and Presentations

We explore the capabilities of volcano opto-acoustics, a promising technique for measuring explosion and infrasound resonance phenomena at open-vent volcanoes. Joint visual and infrasound study at Yasur Volcano (Vanuatu) demonstrate that even consumer-grade cameras are capable of recording infrasound with high fidelity. Passage of infrasonic waves, ranging from as low as 5 Pa to hundreds of Pa, from both explosions and persistent tremor, pressurizes and depressurizes ambient plumes inducing visible vaporization and condensation respectively. Optical tracking of these pressure wavefields can be used to identify spectral characteristics, which vary within Yasur's two deep craters and are distinct for explosion and …

Infrasound Detection Of Approaching Lahars, J. B. Johnson, J. F. Anderson, J. Mock, A. Bosa

Geosciences Faculty Publications and Presentations

Infrasound may be used to detect the approach of hazardous volcanic mudflows, known as lahars, tens of minutes before their flow fronts arrive. We have analyzed signals from more than 20 secondary lahars caused by precipitation events at Fuego Volcano during Guatemala’s rainy season in May through October of 2022. We are able to quantify the capabilities of infrasound monitoring through comparison with seismic data, time lapse camera imagery, and high-resolution video of a well-recorded event on August 17. We determine that infrasound sensors, deployed adjacent to the lahar path and in small-aperture (10 s of meters) arrays, are particularly …

New U–Pb Geochronology For The Central Atlantic Magmatic Province, Critical Reevaluation Of High-Precision Ages And Their Impact On The End-Triassic Extinction Event, Alisson L. Oliveira, Mark D. Schmitz, Corey J. Wall, James L. Crowley, Antomat A. Macêdo Filho, Maria Helena B. M. Hollanda

Geosciences Faculty Publications and Presentations

The end-Triassic extinction (ETE) event represents one of the ‘big five’ episodes of mass extinction. The leading hypothesis for the cause of the ETE is the intrusion of voluminous magmas of the Central Atlantic Magmatic Province (CAMP) into carbon-rich sediments of two South American sedimentary basins, around 201.5 Ma. The timing of dikes and sills emplacement, however, must be considered in light of age models from CAMP rocks occurring in North America. In this work, we present new high-precision ages for critical samples in NE Brazil (201.579 ± 0.057 Ma) and Canada (201.464 ± 0.017 Ma), in order to evaluate …

Critical Workforce Skills For Bachelor-Level Geoscientists: An Analysis Of Geoscience Job Advertisements, G. W. Shafer, K. Viskupic, A. E. Egger

Geosciences Faculty Publications and Presentations

Understanding the skills bachelor-level geoscientists need to enter the workforce is critical to their success. The goal of this study was to identify the workforce skills that are most requested from a broad range of geoscience employers. We collected 3668 job advertisements for bachelor-level geoscientists and used a case-insensitive, code-matching function in Matlab to determine the skills geoscience employers seek. Written communication (67%), field skills (63%), planning (53%), and driving (51%) were most frequently requested. Field skills and data collection were frequently found together in the ads. Written communication skills were common regardless of occupation. Quantitative skills were requested less …

Tectonomagmatic Evolution Of Southwestern Laurentia: Insights From Zircon U-Pb Geochronology And Hafnium Isotopic Composition Of The Red Bluff Granite Suite, West Texas, Usa, Munazzam Ali Mahar, Philip C. Goodell, Jason W. Ricketts, Eric J. Kappus, James L. Crowley, Anthony M. Alvarez

Geosciences Faculty Publications and Presentations

We provide laser ablation–multicollector–inductively coupled plasma–mass spectrometry (LA-MC-ICP-MS) and high-precision chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb ages and Hf isotopic compositions of zircons from the Red Bluff Granite Suite and mafic dikes in the Franklin Mountains of El Paso, Texas, USA. Granitoids exposed in the Franklin Mountains were previously divided into five magmatic stages based on cross-cutting relationships. Major and trace element compositions showed that these granitoids are ferroan, alkaline, and A2 type. Homogeneity in the whole-rock geochemistry suggests that the granite stages are genetically related and share similar petrogenetic histories. Weighted mean zircon ²⁰⁶Pb/²³⁸U …

Advances In Our Understanding Of Pyroclastic Current Behavior From The 1980 Eruption Sequence Of Mount St. Helens Volcano (Washington), Usa, Brittany D. Brand, Aaron A. Marshall

Geosciences Faculty Publications and Presentations

This review summarizes what the volcanology community has learned thus far from studying the deposits of pyroclastic currents (PCs) from the 1980 eruption sequence at Mount St. Helens. The review includes mass flow events during the May 18 eruption, including the lateral blast, the afternoon column collapse and boil-over PC activity, and some aspects of the debris avalanche. We also include a summary of PCs generated in the smaller eruptions following the climactic May 18 event. Our objective is to summarize the state of our understanding of PC transport and emplacement mechanisms from the combination of field and laboratory observations, …

New Insights Into The Age And Origin Of Two Small Cretaceous Seamount Chains Proximal To The Northwestern Hawaiian Ridge, Arturo Sotomayor, Andrea Balbas, Kevin Konrad, Anthony A. P. Koppers, Jasper G. Konter, V. Dorsey Wanless, Thomas F. Hourigan, Christopher Kelley, Nicole Raineault

Geosciences Faculty Publications and Presentations

The Northwestern Hawaiian Ridge is an age-progressive volcanic chain sourced from the Hawaiian mantle plume. Proximal to the Northwestern Hawaiian Ridge are several clusters of smaller seamounts and ridges with limited age constraints and unknown geodynamic origins. This study presents new bathymetric data and ⁴⁰Ar/³⁹Ar age determinations from lava flow samples recovered by remotely operated vehicle (ROV) from two east–west-trending chains of seamounts that lie north of the Pūhāhonu and Mokumanamana volcanoes. The previously unexplored Naifeh Chain (28°48′N,167°48′W) and Plumeria Chain (25°36′N, 164°35′W) contain five volcanic structures each, including three guyots in the Naifeh Chain. New ^{40 …}

Dynamic Mass Loss From Greenland's Marine-Terminating Peripheral Glaciers (1985–2018), Katherine E. Bollen, Ellyn M. Enderlin, Rebecca Muhlheim

Geosciences Faculty Publications and Presentations

Global glacier mass balance decreased rapidly over the last two decades, exceeding mass loss from the Greenland and Antarctic Ice Sheets. In Greenland, peripheral glaciers and ice caps (GICs) cover only ~5% of Greenland's area but contributed ~20% of the island's ice mass loss between 2000 and 2018. Although Greenland GIC mass loss due to surface meltwater runoff has been estimated using atmospheric models, mass lost to changes in ice discharge into oceans (i.e., dynamic mass loss) remains unquantified. We use the flux gate method to estimate discharge from Greenland's 585 marine-terminating peripheral glaciers between 1985 and 2018, and compute …

Lidar And Deep Learning Reveal Forest Structural Controls On Snowpack, Ahmad Hojatimalekshah, Joel Gongora, Josh Enterkine, Nancy F. Glenn, T. Trevor Caughlin, Hans-Peter Marshall, Christopher A. Hiemstra

Geosciences Faculty Publications and Presentations

Forest structure has a strong relationship with abiotic components of the environment. For example, canopy morphology controls snow depth through interception and modifies incoming thermal radiation. In turn, snow water availability affects forest growth, carbon sequestration, and nutrient cycling. We investigated how structural diversity and topography affect snow depth patterns across scales. The study site, Grand Mesa, Colorado, is representative of many areas worldwide where declining snowpack and its consequences for forest ecosystems are increasingly an environmental concern. On the basis of a convolution neural network model (R² of 0.64; root mean squared error of 0.13 m), we …

Operational Water Forecast Ability Of The Hrrr-Isnobal Combination: An Evaluation To Adapt Into Production Environments, Joachim Meyer, John Horel, Patrick Kormos, Andrew Hedrick, Ernesto Trujillo, S. Mckenzie Skiles

Geosciences Faculty Publications and Presentations

Operational water-resource forecasters, such as the Colorado Basin River Forecast Center (CBRFC) in the Western United States, currently rely on historical records to calibrate the temperature-index models used for snowmelt runoff predictions. This data dependence is increasingly challenged, with global and regional climatological factors changing the seasonal snowpack dynamics in mountain watersheds. To evaluate and improve the CBRFC modeling options, this work ran the physically based snow energy balance iSnobal model, forced with outputs from the High-Resolution Rapid Refresh (HRRR) numerical weather prediction model across 4 years in a Colorado River Basin forecast region. Compared to in situ, remotely sensed, …

Precise Zircon U-Pb Dating Of The Mesoproterozoic Gawler Large Igneous Province, South Australia, Elizabeth A. Jagodzinski, Anthony J. Reid, James L. Crowley, Claire E. Wade, Stacey Curtis

Geosciences Faculty Publications and Presentations

The Mesoproterozoic Gawler Range Volcanics and Benagerie Volcanic Suite of the Gawler Craton and Curnamona Province, South Australia, together with associated intrusive magmatism, define an intracontinental, subaerial large igneous province (LIP) preserving an estimated 110 000 km³ of volcanic rock, which hosts one of the world's largest orebodies, the Fe oxide-Cu-Au-U deposit at Olympic Dam, and numerous other related Cu-Au deposits. New high-precision Chemical Abrasion Isotope Dilution Thermal Ionization Mass Spectrometry (CA-TIMS) U-Pb dates on volcanic zircons allow for regional correlations between stratigraphic units of the GRV and BVS, and an understanding of how magmatic styles, temperatures, composition and …

Drivers Of Spatiotemporal Patterns Of Surface Water Inputs In A Catchment At The Rain-Snow Transition Zone Of The Water-Limited Western United States, E. Trujillo, J. Mcnamara

Geosciences Faculty Publications and Presentations

Spatial and temporal dynamics of rainfall and snowmelt (i.e., surface water inputs, SWI) control soil moisture, groundwater recharge, and streamflow at annual, seasonal, and event scales. In the rain-snow transition zone, comprising a large portion of the mountainous western United States, there is limited understanding of the sensitivity of spatiotemporal SWI dynamics across hydrologically variable water years (WYs). We modeled rainfall and snowpack dynamics in a small headwater catchment (1.8 km²) spanning the rain-snow transition in southwestern Idaho, USA, for two hydrologically distinct WYs (2011 and 2014). In wet WY 2011 and dry WY 2014, total precipitation drove …

Elastic Thermobarometry, Matthew J. Kohn, Mattia L. Mazzucchelli, Matteo Alvaro

Geosciences Faculty Publications and Presentations

Upon exhumation and cooling, contrasting compressibilities and thermal expansivities induce differential strains (volume mismatches) between a host crystal and its inclusions. These strains can be quantified in situ using Raman spectroscopy or X-ray diffraction. Knowing equations of state and elastic properties of minerals, elastic thermobarometry inverts measured strains to calculate the pressure-temperature conditions under which the stress state was uniform in the host and inclusion. These are commonly interpreted to represent the conditions of inclusion entrapment. Modeling and experiments quantify corrections for inclusion shape, proximity to surfaces, and (most importantly) crystal-axis anisotropy, and they permit accurate application of the more …

Pacific Northwest Streamflow Data Landscape: A Report From The 2022 Streamflow Data Roundtables, Kendra E. Kaiser

Geosciences Faculty Publications and Presentations

This project was funded by the U.S. Geological Survey Northwest Climate Adaptation Science Center to catalog the location, temporal extent and purpose of non-USGS streamflow datasets. As part of this project, roundtable meetings convened local, state and federal agencies, and nonprofits to explore the complexity of gathering and integrating the identified datasets and identify issues surrounding data-sharing across organizations.

This report synthesizes discussions from each of the state roundtable discussions convened in the spring of 2022, and highlights common challenges and needs across the region. Additional information from organizations not able to be present at the meetings were added after …

New Age Constraints On The Break-Up Of Rodinia And Amalgamation Of Southwestern Gondwana From The Choquequirao Formation In Southwestern Peru, Eben Blake Hodgin, Victor Carlotto, Francis A. Macdonald, Mark D. Schmitz, James L. Crowley

Geosciences Faculty Publications and Presentations

The Choquequirao Formation is a >3 km-thick amphibolite-grade succession that outcrops in the Central Andes of southern Peru. To constrain its age and tectonostratigraphic setting, detrital zircon and metamorphic zircon, titanite, and rutile U–Pb isotopic analyses were conducted. Mantle-derived c. 640 Ma detrital zircons constrain the maximum age of the lower part of the succession and 550–490 Ma metamorphic zircon domains constrain its minimum age. The absence of early Paleozoic detrital zircons suggests that deposition predated early Paleozoic orogenesis in southwestern Gondwana. The close similarity of detrital zircon age spectra to those from sediments deposited on the Arequipa basement suggests …

Causes, Responses, And Implications Of Anthropogenic Versus Natural Flow Intermittence In River Networks, Kendra E. Kaiser

Geosciences Faculty Publications and Presentations

Rivers that do not flow year-round are the predominant type of running waters on Earth. Despite a burgeoning literature on natural flow intermittence (NFI), knowledge about the hydrological causes and ecological effects of human-induced, anthropogenic flow intermittence (AFI) remains limited. NFI and AFI could generate contrasting hydrological and biological responses in rivers because of distinct underlying causes of drying and evolutionary adaptations of their biota. We first review the causes of AFI and show how different anthropogenic drivers alter the timing, frequency and duration of drying, compared with NFI. Second, we evaluate the possible differences in biodiversity responses, ecological functions, …

Carbonates In The Critical Zone, M. D. Covington, J. B. Martin, L. E. Toran, J. L. Macalady, N. Sekhon, P. L. Sullivan, Á. A. García Jr., J. B. Heffernan, W. D. Graham

Geosciences Faculty Publications and Presentations

Earth's Critical Zone (CZ), the near-surface layer where rock is weathered and landscapes co-evolve with life, is profoundly influenced by the type of underlying bedrock. Previous studies employing the CZ framework have focused primarily on landscapes dominated by silicate rocks. However, carbonate rocks crop out on approximately 15% of Earth's ice-free continental surface and provide important water resources and ecosystem services to ∼1.2 billion people. Unlike silicates, carbonate minerals weather congruently and have high solubilities and rapid dissolution kinetics, enabling the development of large, interconnected pore spaces and preferential flow paths that restructure the CZ. Here we review the state …

Uncertainty Of Icesat-2 Atl06- And Atl08-Derived Snow Depths For Glacierized And Vegetated Mountain Regions, Ellyn M. Enderlin, Colten M. Elkin, Madeline Gendreau, H.P. Marshall, Shad O'Neel, Christopher Mcneil, Caitlyn Florentine, Louis Sass

Geosciences Faculty Publications and Presentations

Seasonal snow melt dominates the hydrologic budget across a large portion of the globe. Snow accumulation and melt vary over a broad range of spatial scales, preventing accurate extrapolation of sparse in situ observations to watershed scales. The lidar onboard the Ice, Cloud, and land Elevation, Satellite (ICESat-2) was designed for precise mapping of ice sheets and sea ice, and here we assess the feasibility of snow depth-mapping using ICESat-2 data in more complex and rugged mountain landscapes. We explore the utility of ATL08 Land and Vegetation Height and ATL06 Land Ice Height differencing from reference elevation datasets in two …

Predictions And Drivers Of Sub-Reach-Scale Annual Streamflow Permanence For The Upper Missouri River Basin: 1989–2018, Kendra E. Kaiser

Geosciences Faculty Publications and Presentations

The presence of year-round surface water in streams (i.e., streamflow permanence) is an important factor for identifying aquatic habitat availability, determining the regulatory status of streams, managing land use change, allocating water resources, and designing scientific studies. However, accurate, high resolution, and dynamic prediction of streamflow permanence that accounts for year-to-year variability at a regional extent is a major gap in modeling capability. Herein, we expand and adapt the U.S. Geological Survey (USGS) PRObability of Streamflow PERmanence (PROSPER) model from its original implementation in the Pacific Northwest (PROSPER_PNW) to the upper Missouri River basin (PROSPER_UM), a …

The Formation Of The 8˚20’ N Seamount Chain, East Pacific Rise, Dorsey Wanless

Geosciences Faculty Publications and Presentations

Near-axis seamounts provide a unique setting to investigate three-dimensional mantle processes associated with the formation of new oceanic crust and lithosphere. Here, we investigate the characteristics and evolution of the 8˚20’N Seamount Chain, a lineament of seamounts that extends ~ 175 km west of the East Pacific Rise (EPR) axis, just north of the fracture zone of the Siqueiros Transform Fault. Shipboard gravity, magnetic, and bathymetric data acquired in 2016 are utilized to constrain models of seamount emplacement and evolution. Geophysical observations indicate that these seamounts formed during four distinct episodes of volcanism coinciding with changes in regional plate motion …

Statistical Analysis Of Aquifer Hydraulic Properties By A Continuous Pumping Tomography Test: Application To The Boise Hydrogeophysical Research Site, Kan Bun Cheng, Gedeon Dagan, Warren Barrash, Michael Cardiff, Avinoam Rabinovich

Geosciences Faculty Publications and Presentations

Characterizing aquifer heterogeneity is paramount for accurate flow and transport modeling. In this work, we present a new approach for statistical analysis of hydraulic properties in continuous pumping tomography tests of a phreatic aquifer. The method entails determining equivalent hydraulic conductivity (K_eq), specific storage (S_s,eq), and specific yield (S_y,eq) at many locations in the field and then calculating statistical moments of the equivalent properties, assuming they are random space variables. Equivalent properties are defined as the ones pertinent to a homogeneous aquifer for which the head time …

Emplacement Of The Franklin Large Igneous Province And Initiation Of The Sturtian Snowball Earth, Mark D. Schmitz

Geosciences Faculty Publications and Presentations

During the Cryogenian (720 to 635 Ma ago) Snowball Earth glaciations, ice extended to sea level near the equator. The cause of this catastrophic failure of Earth’s thermostat has been unclear, but previous geochronology has suggested a rough coincidence of glacial onset with one of the largest magmatic episodes in the geological record, the Franklin large igneous province. U-Pb geochronology on zircon and baddeleyite from sills associated with the paleo-equatorial Franklin large igneous province in Arctic Canada record rapid emplacement between 719.86 ± 0.21 and 718.61 ± 0.30 Ma ago, 0.9 to 1.6 Ma before the onset of widespread glaciation. …

Infrasonic Gliding Reflects A Rising Magma Column At Mount Etna (Italy), Mariangela Sciotto, Leighton M. Watson, Andrea Cannata, Massimo Cantarero, Emanuela De Beni, Jeffrey B. Johnson

Geosciences Faculty Publications and Presentations

Infrasound is increasing applied as a tool to investigate magma dynamics at active volcanoes, especially at open-vent volcanoes, such as Mt. Etna (Italy), which are prodigious sources of infrasound. Harmonic infrasound signals have been used to constrain crater dimensions and track the movement of magma within the shallow plumbing system. This study interprets the remarkable systematic change in monotonic infrasound signals preceding a lava fountaining episode at Mt. Etna on 20 February 2021. We model the changing tones (0.7 to 3 Hz fundamental frequency) as a rise in the magma column from 172 ± 25 m below the crater rim …