Probabilistic tephra hazard modeling for the Jemez Mountains Volcanic Field and Valles Caldera, NM USA
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Stratospheric aerosols, most of which originate from explosive volcanic sulfur emissions into the stratosphere, are a key natural driver of climate variability. They are thus a forcing provided by the Coupled Model Intercomparison Project (CMIP) Climate Forcings Task Team to climate modelling groups participating in phase 7 of CMIP. For the historical period, we provide two datasets covering 1750–2023: (i) a volcanic upper tropospheric-stratospheric sulfur emission dataset, documented in a companion paper; and (ii) a stratospheric sulfate aerosol optical property dataset, documented in this paper. For the satellite era (from 1979 onwards), stratospheric aerosol optical properties are derived from the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) dataset. For the pre-satellite era (1750–1978), optical properties are derived from our volcanic SO 2 emission dataset using a new version of the reduced-complexity volcanic aerosol model Easy Volcanic Aerosol (Height) (EVA_H v2). A background, non-volcanic stratospheric aerosol climatology is derived from the 1998–2001 period with a trend over 1850–1978 accounting for increasing anthropogenic aerosols. A monthly stratospheric aerosol climatology is derived from the 1850–2021 average for both pre-industrial and scenario (future) simulations, with a 9-year ramp over 2022–2030 for scenario simulations to ensure a smooth transition from the historical period. CMIP7 simulations should use version 2.2.1 of our pre-industrial and historical datasets, and version 2.2.2 of our scenario dataset. Our methodology to produce historical aerosol optical properties significantly differs from CMIP6 for the pre-satellite era, and the resulting forcings in turn largely differ. In particular, the CMIP6 dataset was mostly based on the sparse and uncertain pyrheliometer record, which resulted in strongly underrepresented emissions from small-to-moderate magnitude eruptions. The resulting bias is addressed in CMIP7, which is entirely emission-derived in the pre-satellite era and uses more recent ice-core-based volcanic sulfur emission inventories than CMIP6. Our approach results in an overall larger volcanic aerosol forcing for CMIP7, with the 1850–2014 mean mid-visible global mean stratospheric aerosol optical depth (SAOD) in CMIP7 (0.0138) being 29 % higher than in CMIP6 (0.0107). The pre-industrial mean of the same variable is 26 % higher in CMIP7 (0.0135, derived from the historical 1850–2021) than CMIP6 (0.0107, derived from the historical 1850–2014 mean). Using a reduced-complexity climate model, we simulate a global mean surface temperature that is 0.07 °C colder for 1850–1900 when using the CMIP7 dataset instead of CMIP6, whereas 2000–2014 is 0.03 °C warmer in CMIP7. Our dataset also exhibits lower forcing for 1960–1980, resulting in temperatures 0.06 °C warmer when averaged over 1960–1990, a period for which CMIP6 climate models exhibit a cold bias. Given the large uncertainties characterizing the dataset, in particular for the pre-satellite era, we advise against treating the CMIP7 or CMIP6 dataset as uniquely superior for any specific year and highlight the need for further evaluation. We conclude the study by discussing sources of uncertainty for the dataset, future research avenues to improve it, as well as requirements to operationalize the production of the dataset, i.e. extend it and update it on an annual basis instead of every 5–7 years following CMIP cycles.
Wavelength-dependent mid-infrared (400–1,800 cm -1 ; 5.56–25 μm) optical constants (real and imaginary indices of refraction; n and k) are determined using reflectance spectra at a spectral sampling of 4 cm -1 for several silicate glasses of varying SiO 2 wt% which include (a) basaltic volcanic glass from Kilauea, Hawaii, (b) synthetic andesite, (c) two synthetic dacites, (d) obsidian volcanic glass from Mount Lassen, California, (e) synthetic rhyodacite, and (f) rhyolitic volcanic glass from Mexico. Because glasses are optically isotropic, no specific orientation was required for spectral measurements, and polished glass samples were measured at random orientations using micro-FTIR spectrometer. Lorentz-Lorenz dispersion theory and Fresnel reflectance model for high symmetry materials were used to model the optical constants by optimizing oscillator parameters in modeled spectra to match laboratory spectra. A genetic algorithm (GA) approach automatically finds the natural oscillators and their parameters for each spectrum, and then uses these parameters in a non-linear least squares optimization routine. The study compared spectral parameters such as Christiansen feature, reststrahlen bands (RBs), and the peaks centered around 860–1,100 (peak 1) and 400–480 cm -1 (peak 2) of n and k to their respective SiO 2 wt% of the glasses. CF, RBs, peak 1 of n and k, and peak 2 of n shift to higher wavenumbers with increased SiO 2 wt%, whereas peak 2 of k shifts to lower wavenumbers with increased SiO 2 wt%. Derived optical constants of these glasses will improve quantitative abundance mapping of volcanic materials on the surfaces of silicate targets in the Solar System.
The angrites are a class of achondrites that encompass a wide range of igneous textures from quenched, volcanic, and subvolcanic stones to slowly cooled, plutonic rocks. The compositions of the various geochemical reservoirs generating this variety of rocks have not been investigated fully because historically the numbers and masses of angrites available for study have been quite small. However, the rapid increase in meteorites from Northwest Africa (NWA) has enabled a renewed opportunity for such an investigation. In particular, three samples, NWA 2999, 4590, and 4801, have provided a new window into our understanding of the origin and petrogenesis of the coarse-grained (plutonic) angrites. We report here the trace element abundances in individual mineral grains of pyroxene, plagioclase, olivine/kirschsteinite, and phosphate and in the whole-rock samples. We utilize these data to constrain the petrogenetic history of each of these samples on the angrite parent body (e.g., parental melt compositions and oxygen fugacity conditions) and assess genetic relationships to previously investigated angrites. The trace element abundances in each of the three coarse-grained angrites studied here indicate a unique history for each. The observed trace element abundances and patterns in NWA 2999 show similarities with previously studied fine-grained, volcanic angrites and potentially indicate a common geochemical source reservoir, even though NWA 2999 is temporally distinct from the volcanic angrites. In contrast, NWA 4590 has trace element characteristics analogous to other coarse-grained angrites (e.g., Lewis Cliff [LEW] 86010), suggesting that these samples originated from geochemically similar source reservoirs. The third angrite, NWA 4801, exhibits geochemical characteristics most similar to the plutonic, coarse-grained angrites, but also appears to have some affinities in its trace element characteristics to the volcanic, fine-grained angrites. This suggests that NWA 4801 may represent a petrogenetic link between two distinct geochemical reservoirs on the angrite parent body. In aggregate, the trace element distributions in these three plutonic angrites suggest that while they may have originated up to several million years after the fine-grained angrites, they sampled a range of source reservoirs on the angrite parent body. Here, some of these source reservoirs were likely similar to those of the fine-grained angrites, but others had distinct geochemical characteristics.
Hydrothermally active submarine volcanoes are mineral-rich biological oases contributing significantly to chemical fluxes in the deep sea, yet little is known about the microbial communities inhabiting these systems. Here we investigate the diversity of microbial life in hydrothermal deposits and their metagenomics-inferred physiology in light of the geological history and resulting hydrothermal fluid paths in the subsurface of Brothers submarine volcano north of New Zealand on the southern Kermadec arc. From metagenome-assembled genomes we identified over 90 putative bacterial and archaeal genomic families and nearly 300 previously unknown genera, many potentially endemic to this submarine volcanic environment. While magmatically influenced hydrothermal systems on the volcanic resurgent cones of Brothers volcano harbor communities of thermoacidophiles and diverse members of the superphylum “DPANN,” two distinct communities are associated with the caldera wall, likely shaped by two different types of hydrothermal circulation. The communities whose phylogenetic diversity primarily aligns with that of the cone sites and magmatically influenced hydrothermal systems elsewhere are characterized predominately by anaerobic metabolisms. These populations are probably maintained by fluids with greater magmatic inputs that have interacted with different (deeper) previously altered mineral assemblages. However, proximal (a few meters distant) communities with gene-inferred aerobic, microaerophilic, and anaerobic metabolisms are likely supported by shallower seawater-dominated circulation. Furthermore, mixing of fluids from these two distinct hydrothermal circulation systems may have an underlying imprint on the high microbial phylogenomic diversity. Collectively our results highlight the importance of considering geologic evolution and history of subsurface processes in studying microbial colonization and community dynamics in volcanic environments.
Lava caves, tubes, and fumaroles in Hawai‘i present a range of volcanic, oligotrophic environments from different lava flows and host unexpectedly high levels of bacterial diversity. These features provide an opportunity to study the ecological drivers that structure bacterial community diversity and assemblies in volcanic ecosystems and compare the older, more stable environments of lava tubes, to the more variable and extreme conditions of younger, geothermally active caves and fumaroles. Using 16S rRNA amplicon-based sequencing methods, we investigated the phylogenetic distinctness and diversity and identified microbial interactions and consortia through co-occurrence networks in 70 samples from lava tubes, geothermal lava caves, and fumaroles on the island of Hawai‘i. Our data illustrate that lava caves and geothermal sites harbor unique microbial communities, with very little overlap between caves or sites. We also found that older lava tubes (500–800 yrs old) hosted greater phylogenetic diversity (Faith's PD) than sites that were either geothermally active or younger (<400 yrs old). Geothermally active sites had a greater number of interactions and complexity than lava tubes. Average phylogenetic distinctness, a measure of the phylogenetic relatedness of a community, was higher than would be expected if communities were structured at random. This suggests that bacterial communities of Hawaiian volcanic environments are phylogenetically over-dispersed and that competitive exclusion is the main driver in structuring these communities. This was supported by network analyses that found that taxa (Class level) co-occurred with more distantly related organisms than close relatives, particularly in geothermal sites. Network “hubs” (taxa of potentially higher ecological importance) were not the most abundant taxa in either geothermal sites or lava tubes and were identified as unknown families or genera of the phyla, Chloroflexi and Acidobacteria. These results highlight the need for further study on the ecological role of microbes in caves through targeted culturing methods, metagenomics, and long-read sequence technologies.
Climate instability driven by emission of volatiles during emplacement of large igneous provinces (LIPs) is frequently invoked as a potential cause of mass extinctions. However, documenting this process in the geologic record requires a holistic understanding of eruption rates, the location of eruptive centers, and potential sources of climate-changing volatiles. We present new chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb zircon geochronology from Malwa Plateau basalts on the northern margin of the Deccan LIP, India. These basalts have been previously interpreted as either an extension of the province's main volcanic stratigraphy or as an independent eruptive center active up to millions of years prior to the main eruptive phase. Our data instead demonstrate that the lower Malwa Plateau basalts are temporally correlative with the first pulse of Deccan volcanism and provide new constraints on its initiation and duration. Paleomagnetic data further indicate that upper Malwa Plateau basalts may be age-equivalent to the second, third, and fourth pulses of Deccan volcanism. The relative thicknesses of age-equivalent packages of basalt are consistent with eruption of the Deccan LIP from a southward-migrating eruptive center. The first eruptive pulse is coeval with a ~200 kyr Late Maastrichtian warming event preserved globally in contemporaneous stratigraphic sections. We propose that the first pulse of Deccan magmatism was more voluminous in the north, where it erupted through organic-rich sedimentary rocks of the Narmada-Tapti rift basin. Thermal metamorphism of these sediments could have been a source of sufficient CO 2 to drive the Late Maastrichtian warming event, which if true would reconcile the apparent dampened warming signals associated with later Deccan eruptive pulses.
Plume volcanism may sample mantle sources deeper than mid-ocean ridge and arc volcanism. Ocean island basalts (OIBs) are commonly related to plume volcanism, and their diverse isotopic and elemental compositions can be described using a limited number of mantle endmembers. However, the origins and depths of these mantle endmembers are highly debated. Here we show that the HIMU (high μ, μ = 238 U/ 204 Pb) endmember may reside in the transition zone. Specifically, we report the geochemical signature of a high-pressure multiphase diamond inclusion, entrapped at 420–440 km depth and 1450 ±50 K, which matches exactly the geochemical patterns of the HIMU-rich OIBs. Since the HIMU component is variably sampled by almost all OIBs, our finding implies that the transition zone causes a major overprint of the geochemical features of mantle plumes. Furthermore, some mantle plumes, like those feeding Bermuda, St Helena, Tubuai and Mangaia, appear to be dominated by this source. Furthermore, our finding highlights the importance of the transition zone in highly incompatible element budget of the mantle.
The Snake River Plain (SRP) volcanic province overlies the track of the Yellowstone hotspot, a thermal anomaly that extends deep into the mantle. Most of the area is underlain by a basaltic volcanic province that overlies a mid-crustal intrusive complex, which in turn provides the long-term heat flux needed to sustain geothermal systems. Previous studies have identified several known geothermal resource areas within the SRP. For the geothermal study presented herein, our goals were to: (1) adapt the methodology of Play Fairway Analysis (PFA) for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal PFA model for the SRP and identify the most promising plays, using GIS-based software tools that are standard in the petroleum industry. The study focused on identifying three critical resource parameters for exploitable hydrothermal systems in the SRP: heat source, reservoir and recharge permeability, and cap or seal. Data included in the compilation for heat source were heat flow, distribution and ages of volcanic vents, groundwater temperatures, thermal springs and wells, helium isotope anomalies, and reservoir temperatures estimated using geothermometry. Reservoir and recharge permeability was inferred from the analysis of stress orientations and magnitudes, post-Miocene faults, and subsurface structural lineaments based on magnetics and gravity data. Data for cap or seal included the distribution of impermeable lake sediments and clay-seal associated with hydrothermal alteration below the regional aquifer. These data were used to compile Common Risk Segment maps for heat, permeability, and seal, which were combined to create a Composite Common Risk Segment map for all southern Idaho that reflects the risk associated with geothermal resource exploration and identifies favorable resource tracks. Our regional assessment indicated that undiscovered geothermal resources may be located in several areas of the SRP. Two of these areas, the western SRP and Camas Prairie, were selected for more detailed assessment, during which heat, permeability, and seal were evaluated using newly collected field data and smaller grid parameters to refine the location of potential resources. These higher resolution assessments illustrate the flexibility of our approach over a range of scales.
Magma and related hydrothermal fluid movement, and their interaction with solid rock, in active volcanic regions, generate a wide variety of seismic waves whose characterization can mitigate the risk of a potential eruption. Located in the western region of the Trans-Mexican Volcanic Belt, Ceboruco Volcano, whose last eruptive period was 1870–1875, is considered to be one of the most hazardous volcanoes in Mexico. We have conducted a detailed study of the seismicity in the surroundings of Ceboruco's volcanic edifice to assess the current state of this volcano. A dense temporary seismic network with 25 seismic stations in an area of 16 km × 16 km was deployed between November 2016 and July 2017, as part of the P-24 project of the CeMIEGeo consortium; this effort has allowed the detection of 81 earthquakes concentrated beneath the crater with depths between 4 and 8 km. In this study, we observe that the recorded seismicity occurs in swarms, and we specifically identify four sequences that we characterize in detail via the first focal mechanisms available for this volcano. Our results suggest a change in the local seismicity distribution compared to earlier observations, which reported seismic activity near the volcano edifice associated with fluid migration along zones of weakness related to the extensional stresses of the Tepic-Zacoalco rift. The changes in seismic patterns and obtained focal mechanisms are consistent with observed fluid effects at many geothermal sites worldwide, but also could suggest resumption of activity at this currently dormant volcano.
Abstract Volcanic eruptions are shallow phenomena that represent the final stage of density- and viscosity- driven processes of melt migration from source rocks at upper mantle depths. In this experimental study, we investigated the effect of pressure (0.7–7.0 GPa) and temperature (1335–2000 °C) on the viscosity and the atomic melt structure of a synthetic anhydrous primitive alkaline basalt, an analogue of the pre-eruptive magma that likely feeds the Campi Flegrei Volcanic District at present day. Obtained viscosities (0.5–3.0 Pa s), mobility (0.1–0.4 g cm 3 Pa −1 s −1 ) and ascent velocity (1.5–6.0 m yr −1 ) are presented to support geochemical and geophysical observations of Campi Flegrei as a critical volcanic district currently undergoing gradual magma recharge at depth.
Venus and Earth are rocky planets of roughly the same size and bulk density, yet their surface volcanic and tectonic features appear substantially different. On Venus, the coexistence of large volcanic highlands—interpreted as the surface expression of long-lived mantle plumes—alongside coronae, smaller features thought to be caused by transient thermal diapirs, remains enigmatic. Using two-dimensional numerical models of mantle convection with sharp and broad mineral phase transitions for pyrolite, we show that both scales of upwellings can be generated in a stagnant lid planet with an interior temperature 250 to 400 K warmer than Earth’s. The smaller plumes originate from a ~600 km deep internal layer that exists as a consequence of the different sequence of mineral phase transitions that occur in warmer mantles less processed and differentiated by partial melting and volcanism. Future models that include melting will provide further tests of our hypothesis.
Plant-associated microbial communities have diverse phenotypic effects on their hosts that are only beginning to be revealed. We hypothesized that morpho-physiological variations in the tropical tree Tabebuia heterophylla, observed on different geological substrates, arise in part due to microbial processes in the rhizosphere. We characterized the microbiota of the rhizosphere and soil communities associated with T. heterophylla trees in high and low altitude sites (with varying temperature and precipitation) of volcanic, karst and serpentine geologies across Puerto Rico. We sampled 6 areas across the island in three geological materials including volcanic, serpentine and karst soils. Collection was done in 2 elevations (>450m and 0-300m high), that included 3 trees for each site and 4 replicate soil samples per tree of both bulk and rhizosphere. Genomic DNA was extracted from 144 samples, and 16S rRNA V4 sequencing was performed on the Illumina MiSeq platform. Proteobacteria, Actinobacteria, and Verrucomicrobia were the most dominant phyla, and microbiomes clustered by geological substrate and elevation. Volcanic samples were enriched in Verrucomicrobia; karst was dominated by nitrogen-fixing Proteobacteria, and serpentine sites harbored the most diverse communities, with dominant Cyanobacteria. Sites with similar climates but differing geologies showed significant differences on rhizobiota diversity and composition demonstrating the importance of geology in shaping the rhizosphere microbiota, with implications for the plant's phenotype. Our study sheds light on the combined role of geology and climate in the rhizosphere microbial consortia, likely contributing to the phenotypic plasticity of the trees.
The Fe 3+ /Fe T ratios (Fe 3+ /[Fe 2+ +Fe 3+ ]) in minerals can be used to understand their crystallization and post-crystallization conditions. However, as natural minerals are often zoned and contain inclusions, bulk techniques, e.g., wet chemistry, may not provide accurate Fe 3+ /Fe T values for a single phase of interest. We determined Fe 3+ /Fe T ratios of amphiboles in different crystallographic orientations by single-crystal synchrotron Mössbauer spectroscopy (SMS) in energy and time domain modes from four volcanic localities (Long Valley Caldera, Mount St. Helens, Lassen Volcanic Center, U.S.A., and Mt. Pinatubo, Philippines). The high spatial resolution (as low as 12 × 12 μm spot size) and standard-free nature of SMS allow the detection of intra-grain compositional heterogeneities in Fe 3+ /Fe T with relatively low uncertainties. Here, we combine SMS with major element compositions, water contents, and hydrogen isotope compositions to document the Fe 3+ /Fe T ratios as a function of mineral composition and post-crystallization dehydrogenation. Spectra were fitted with up to five distinct sites: ferrous iron on M(1), M(2), M(3), and ferric iron on M(2) and M(3), consistent with X-ray diffraction studies on single crystals of amphibole. The Fe 3+ /Fe T ratios range from 0.14 ± 0.03 (Long Valley Caldera), 0.51 to 0.63 ± 0.02 (representing intra-grain heterogeneities, Mount St. Helens) to 0.86 ± 0.03 (Lassen Volcanic Center). The latter grain experienced post-crystallization dehydrogenation, shown by its low water content (0.6 ± 0.05 wt%) and its elevated hydrogen isotope composition (δD = +25 ± 3‰ relative to SMOW). The Fe 3+ /Fe T ratios of 0.62 ± 0.01 and 0.20 ± 0.01 of two Mt. Pinatubo grains correlate with high-Al 2 O 3 cores and low-Al 2 O 3 rims and smaller phenocrysts in the sample, respectively. This study shows that SMS is capable of distinguishing two different domains with dissimilar Fe 3+ /Fe T values formed under different crystallization conditions, demonstrating that SMS in combination with major element, water, and hydrogen isotope compositions allows the interpretation of amphibole Fe 3+ /Fe T ratios in the context of crystallization and post-crystallization processes.
The June 15, 1991 Mt. Pinatubo eruption is simulated in E3SM by injecting 10 Tg of SO 2 gas in the stratosphere, turning off prescribed volcanic aerosols, and enabling E3SM to treat stratospheric volcanic aerosols prognostically. This experimental prognostic treatment of volcanic aerosols in the stratosphere results in some realistic behaviors (SO 2 evolves into H 2 SO 4 which heats the lower stratosphere), and some expected biases (H 2 SO 4 aerosols sediment out of the stratosphere too quickly). Climate fingerprinting techniques are used to establish a Mt. Pinatubo fingerprint based on the vertical profile of temperature from the E3SMv1 DECK ensemble. By projecting reanalysis data and preindustrial simulations onto the fingerprint, the Mt. Pinatubo stratospheric heating anomaly is detected. Projecting the experimental prognostic aerosol simulation onto the fingerprint also results in a detectable heating anomaly, but, as expected, the duration is too short relative to reanalysis data.
The Snake River Plain (SRP) volcanic province overlies a thermal anomaly that extends deep into the mantle; it represents one of the highest heat flow provinces in North America. The Yellowstone hotspot continues to feed a magma system that underlies much of southern Idaho and has produced basaltic volcanism as young as 2000 years old. It has been estimated to host up to 855 MW of potential geothermal power production, most of which is associated with the Snake River Plain volcanic province in Idaho, which lies outside the area of Yellowstone National Park (Neely, K.W. and Galinato, G., 2007, Geothermal power generation in Idaho: an overview of current developments and future potential, Open File Report, Idaho Office of Energy Resources).
Abstract. This paper describes the addition of a stratospheric prognostic aerosol (SPA) capability – developed with the goal of accurately simulating sulfate aerosol formation and evolution in the stratosphere – in the Department of Energy (DOE) Energy Exascale Earth System Model, version 2 (E3SMv2). The implementation includes changes to the four-mode Modal Aerosol Module microphysics in the stratosphere to allow for larger particle growth and more accurate stratospheric aerosol lifetime following the Pinatubo eruption. E3SMv2-SPA reasonably reproduces stratospheric aerosol lifetime, burden, aerosol optical depth, and top-of-atmosphere flux when compared to remote sensing observations. E3SMv2-SPA also has close agreement with the interactive chemistry–climate model CESM2-WACCM (Community Earth System Model version 2–Whole Atmosphere Community Climate Model) – which has a more complete chemical treatment – and the observationally constrained, prescribed volcanic aerosol treatment in E3SMv2. Global stratospheric aerosol size distributions identify the nucleation and growth of sulfate aerosol from volcanically injected SO2 from both major and minor volcanic eruptions from 1991 to 1993. The modeled aerosol effective radius is consistently lower than satellite and in situ measurements (max differences of ∼ 30 %). Comparisons with in situ size distribution samples indicate that this simulated underestimation in both E3SMv2-SPA and CESM2-WACCM is due to overly small accumulation and coarse-mode aerosols 6–18 months post-eruption, with E3SMv2-SPA simulating ∼ 50 % of the coarse-mode geometric mean diameters of observations 11 months post-eruption. Effective radii from the models and observations are used to calculate offline scattering and absorption efficiencies to explore the implications of smaller simulated aerosol size for the Pinatubo climate impacts. Scattering efficiencies at wavelengths of peak solar irradiance (∼ 0.5 µm) are 10 %–80 % higher for daily samples in models relative to observations through 1993, suggesting higher diffuse radiation at the surface and a larger cooling effect in the models due to the smaller simulated aerosol; absorption efficiencies at the peak wavelengths of outgoing terrestrial radiation (∼ 10 µm) are 15 %–40 % lower for daily samples in models relative to observations, suggesting an underestimation in stratospheric heating in the models due to the smaller simulated aerosol. These potential biases are based on aerosol size alone and do not take into account differences in the aerosol number. The overall agreement of E3SMv2-SPA with observations and its similar performance to the well-validated CESM2-WACCM makes E3SMv2-SPA a viable alternative to simulating climate impacts from stratospheric sulfate aerosols.