Search NASA⌕ Search

SEARCH · Search NASA

Results for “initial conditions”

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.

At least 181 records · Page 10

On the Correspondence between Seasonal Forecast Biases and Long-Term Climate Biases in Sea Surface Temperature

In this paper, the correspondence between mean sea surface temperature (SST) biases in retrospective seasonal forecasts (hindcasts) and long-term climate simulations from five global climate models is examined to diagnose the degree to which systematic SST biases develop on seasonal time scales. The hindcasts are from the North American Multi-Model Ensemble and the climate simulations are from the Coupled Model Intercomparison Project. The analysis suggests that most robust climatological SST biases begin to form within 6 months of a realistically initialized integration, although the growth rate varies with location, time, and model. In regions with large biases, interannual variability and ensemble spread is much smaller than the climatological bias. Additional ensemble hindcasts of the Community Earth System Model with a different initialization method suggest that initial conditions do matter for the initial bias growth, but the overall global bias patterns are similar after 6 months. A hindcast approach is more suitable to study biases over the tropics and sub-tropics than over the extra-tropics because of smaller initial biases and faster bias growth. The rapid emergence of SST biases makes it likely that fast processes with times scales shorter than the seasonal time scales in the atmosphere and upper ocean are responsible for a substantial part of the climatological SST biases. Studying the growth of biases may provide important clues to the causes and ultimately the amelioration of these biases. Further, initialized seasonal hindcasts can profitably be used in the development of high-resolution coupled ocean-atmosphere models.

54 ENVIRONMENTAL SCIENCES↗

Assimilating Scanning Radar Data into Cloud-Scale Models

The aim of this project is to improve the initial conditions of large eddy simulations (LES) of deep convective storms by assimilating scanning radar reflectivity and radial velocity data as well as local environmental observations. Accomplishing this goal often results in more accurate numerical forecasts of clouds and precipitation than if the background initial meteorological conditions are solely relied upon for convection initiation and intensification, as is typical of most operational and research models. Conducted in a storm-scale LES framework, these more accurate numerical analyses can serve as tool to better understand and parameterize land-atmosphere coupling, boundary layer turbulence, and cloud processes, each used for model parameterization development of convection. Our work examined the sensitivity of the numerical representation of cloud-scale wind and microphysical structures to: i) the assimilation of radiosonde observations characterizing the near-storm inflow, ii) the assimilation of radar data of varying resolution for model simulations with grid spacings of 1-km and 0.25-km, and iii) to the assimilation of three-dimensional plan-projection-indicator (PPI) and range-height indicator (RHI) radar volumes. Few studies have explored deep convective storm analyses and forecasts produced by assimilation of radar observations to cloud-scale LES. Furthermore, to our knowledge, no storm-scale DA study has explored the value of assimilating high-resolution sector RHI volumes to determine if they can improve the retrieval of the vertical kinematic and microphysical storm structure aloft over the use of traditionally utilized PPI scans.

54 ENVIRONMENTAL SCIENCES↗

Post-irradiation examination of legacy high burnup fuel to support safety testing

Safety/transient testing to evaluate performance under off-normal conditions is an essential pillar for both the development of Accident Tolerant Fuels (ATF) and the optimization of fuel operation economics beyond current discharge burnups. Among other factors, the successful interpretation of the transient testing results relies upon the knowledge of the initial conditions of the test, including the characteristics of the fuel system under scrutiny. When testing pre-irradiated material, the assumptions that the fuel and the cladding still have the same properties as in the pre-irradiation stage is obviously wrong and could affect the results of the test. This is particularly true the more burnup accumulates in the fuel rod and irradiation progresses. The knowledge of the initial microstructure of both fuel and cladding allows a clearer interpretation of the subsequent transient testing results, provides validation of the physical phenomena underlying the model predictions and eliminates the uncertainties related to the limited knowledge of the sample status before the test. One example is the phenomenon of fine fragmentation that occurs in Light Water reactor (LWR) fuel. During a Loss of Coolant Accident (LOCA) or Reactivity Initiated Transient (RIA) the fuel can severely fragment. During LOCA, high burnup fuel tend to finely fragment, which has raised safety concerns due to the increased likelihood of dispersal of such small particles once the cladding has burst and due to the increased fission gas release. Therefore, efforts have been devoted to the assessment of a pulverization threshold that could determine the conditions under which fine fragmentation is predominant. However, the lack of information regarding the initial conditions of the fuel, and the connections between those conditions and the pre-transient irradiation history, have hindered the development of a fully mechanistic fragmentation and pulverization criterion. The empirical relationships rely on conservative estimations, due to the lack of information on critical material properties and characteristics. More generally, experimental evidence of the irradiation-induced modifications at microstructural scale are necessary to determine the behavior of the material at the macroscopic scale, with the latter being the one of technological interest. Significant progress has been made in the last two decades in the developments of analytical materials science techniques that can be applied to highly radioactive materials, such as high burnup fuels. The availability of new techniques and the improvement of existing ones has enabled investigations previously not possible that can deepen the understanding of the fuel characteristics and properties at high burnup. The better knowledge of material behavior and irradiation-induced phenomena could help the prediction of its performance. In this context, the scope of the present work is to apply a wide portfolio of advanced characterization techniques to determine properties that are relevant for safety and performance. The results are interpreted in the context of engineering scale post-irradiation examinations and available information on the irradiation conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Process development and scale-up of value added molecules (CRADA Final Report)

In this project, a genetically modified Saccharomyces cerevisiae was used to produce Histamine via fermentation of dextrose. Zymergen provided a tech transfer document based on their 1 L studies, where they observed a titer of 0.9 g/L with this strain. The ABPDU staff performed two bench scale campaigns, with four 2 L reactors in each campaign, to carry out tech transfer and process optimization studies. Optimization was pursued by varying initial glucose concentration, feed initiation conditions, feed strategies, vitamin supplementation, inoculum size, etc. We chose a 2 L process condition, with 4% inoculum volume, continuous feed initiated after complete ethanol depletion, that yielded the most promising titers (2.41 g/L) for scale up to 300 L.

60 APPLIED LIFE SCIENCES↗

Utah FORGE: Data for 3-D Model Development - Lithology, Temperature, Pressure, and Stress

This submission comprises two downloadable .zip archives. Each archive contains spreadsheets with 3-D data and a .txt document describing the data. The Mesh Files .zip download contains data regarding the lithologic contacts of the granitoid and overlying sedimentary basin fill. The Initial Conditions .zip download contains data regarding temperature, pressure, and stress. These data were generated during Utah FORGE Phase 2 and exported from Seequent Leapfrog Geothermal. They are intended for use in 3-D earth modeling.

15 GEOTHERMAL ENERGY↗

Prioritizing the selection of CMIP6 model ensemble members for downscaling projections of CONUS temperature and precipitation

Given the mismatch between the large volume of data archived for the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and limited personnel and computational resources for downscaling, only a small fraction of the CMIP6 archive can be downscaled. In this work, we develop an approach to robustly sample projected hydroclimate states in CMIP6 for downscaling to test whether the selection of a single initial condition (IC) ensemble member from each CMIP6 model is sufficient to span the range of modeled hydroclimate over the conterminous United States (CONUS) and CONUS sub-regions. We calculate the pattern-centered root mean square difference of IC ensemble member anomalies relative to each model’s historical climatology for shared socioeconomic pathway (SSP) projections over 30-year time periods and compare the ratio of inter-model to intra-model variability for this metric. Regardless of SSP, inter-model variability is generally much greater than intra-model variability at the scales of the CONUS as a whole, as well as for most CONUS sub-regions. However, for some variables and scenarios, inter- and intra-model variability are similar at sub-CONUS scales, indicating that selecting a single IC ensemble member per model may be sufficient to sample the range of projected hydroclimate states in the 21st Century across CONUS, but for specific regions and variables, more careful selection of ensemble members may be necessary. Regionally resolved Taylor diagrams identify where more IC ensemble member downscaling efforts should be focused if resources are available to do so. Our results suggest that, with parsimonious sampling, the requisite computational expense of downscaling temperature and precipitation fields over the CONUS for subsequent CMIP activities may increase only marginally despite the great increase in data volumes with each successive CMIP phase.

54 ENVIRONMENTAL SCIENCES↗

Local-thermal-gradient and large-scale-circulation impacts on turbine-height wind speed forecasting over the Columbia River Basin

Abstract. We investigate the sensitivity of turbine-height wind speed forecast to initial condition (IC) uncertainties over the Columbia River Gorge (CRG) and Columbia River Basin (CRB) for two typical weather phenomena, i.e., local-thermal-gradient-induced marine air intrusion and a cold frontal passage. Four types of turbine-height wind forecast anomalies and their associated IC uncertainties related to local thermal gradients and large-scale circulations are identified using the self-organizing map (SOM) technique. The four SOM types are categorized into two patterns, each accounting for half of the ensemble members. The first pattern corresponds to IC uncertainties that alter the wind forecast through a modulating weather system, which produces the strongest wind anomalies in the CRG and CRB. In the second pattern, the moderate uncertainties in local thermal gradient and large-scale circulation jointly contribute to wind forecast anomaly. We analyze the cross section of wind and temperature anomalies through the gorge to explore the evolution of vertical features of each SOM type. The turbine-height wind anomalies induced by large-scale IC uncertainties are more concentrated near the front. In contrast, turbine-height wind anomalies induced by the local IC thermal uncertainties are found above the surface thermal anomalies. Moreover, the wind forecast accuracy in the CRG and CRB is limited by IC uncertainties in a few specific regions, e.g., the 2 m temperature within the basin and large-scale circulation over the northeast Pacific around 140∘ W.

17 WIND ENERGY↗

Quantification of Aerosol Transmission through Stress Corrosion Crack-Like Geometries

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the commercial practice of using more significant backfill pressures in the canister, up to approximately 800 kPa. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the rates of gas flow and aerosol transmission of a spent fuel surrogate through an engineered microchannel with dimensions representative of an SCC were evaluated experimentally using coupled mass flow and aerosol analyzers. The microchannel was formed by mating two gage blocks with a linearly tapering slot orifice nominally 13 μm (0.005 in.) tall on the upstream side and 25 μm (0.0010 in.) tall on the downstream side. The orifice is 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, CeO 2 , were seeded and mixed with either helium or air inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to add to previous testing that characterized SCCs under well-controlled boundary conditions through the inclusion of testing improvements that establish initial conditions in a more consistent way. These ongoing testing efforts are focused on understanding the evolution in both size and quantity of a hypothetical release of aerosolized spent fuel particles from failed fuel to the canister interior and ultimately through an SCC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Preshot Report for Carbonara: Dynamics of a Boron Carbide Sandwich

The main purpose of the Carbonara experiment is to provide much-needed data to inform physics modeling of boron carbide (B 4 C) subjected to a sweeping high explosive (HE) detonation. The proposed experiment will include three separate shots, each employing a rectangular slab of HE (4” × 3” × 1”) sandwiched between two wedges of B 4 C, see Fig.1. In each case the top B 4 C wedge tapers from a height of 2.5 cm to 0.5 cm, while the bottom wedge tapers from 1.5 cm to 0.5 cm. The wedges are 4” in the direction of the proton beam, and 3” in the direction transverse to the beam. Three shots are currently planned, each one using a different HE: PBX 9501, PBX 9502, and PBX 9701. The HE will be initiated using a line wave generator (LWG) to preserve planar symmetry, in conjunction with a Comp. B booster to smooth out the profile for the initial condition on the main charge and to aid in initiation of the PBX 9502 shot. In each case the booster spans the entire length and height of the main HE slab and extends 0.5” in toward the LWG.

36 MATERIALS SCIENCE↗

Modeling the Effects of Microwave Heating on Densification in Chemical Vapor Infiltration

Microwave heating has great potential to accelerate the synthesis of ceramic matrix composites (CMCs) by a process called Chemical Vapor Infiltration (CVI). In CVI, reactive gases ingress a porous preform and undergo chemical transformation to deposit solid ceramic phase within the pores at high temperature thus, densifying the preform. However, the competing effects of chemical kinetics and gas transport are known to result in non-uniform densification as the outer surfaces of the preform experience faster depositions compared to the core. Achieving spatial temperature control plays a key role improving the densification quality. Microwave heating can potentially create temperature inversion such that the core of the preform is hotter than the outer surface and subsequently, lead to improved densification while keeping the manufacturing times and costs low. In the present work, a computational modeling strategy has been developed that accounts for the key physical phenomena responsible for densification of porous preforms using microwave heating. Specifically, a chemical kinetics model has been formulated for Silicon Carbide (SiC) deposition from MTS/H2 precursor. The model is implemented in a pore-resolved reactive transport solver, Quilt. The CVI simulations are performed for several conditions. Initially, parameterized temperature control is imposed to identify optimum conditions for good densification quality at a fraction of processing time. Further, simulations of microwave heating of porous SiC preforms are performed using OpenFOAM. A strategy to achieve and enhance temperature inversion is identified. The resulting temperature profiles are used in the pore-resolved densification simulations to analyze the densification behavior. It is observed that the temperature inversion achieved by microwave heating leads to significant improvements in densification quality and at the same time, keeps the manufacturing time low.

36 MATERIALS SCIENCE↗

Uncertainties in Local Intense Precipitation Flood Modeling

In the Task 1 report for the U.S. Nuclear Regulatory Commission (NRC) funded Local Intense Precipitation (LIP) Probabilistic Flood Hazard Assessment (PFHA) Pilot Study, the Pacific Northwest National Laboratory (PNNL) reviewed hydrologic and hydraulic modeling approaches for estimation of LIP flooding and their implementation in readily available simulation software packages (Prasad and Yuan 2020). NRC’s interest is in developing guidance for performing LIP PFHA at nuclear power plant (NPP) sites. This report describes the work done under Task 2 of the study that focused on describing uncertainties associated with LIP flood modeling. Two main types of uncertainties exist in LIP flood modeling: (1) aleatory uncertainties that arise from the inherent natural variability of the hydrometeorological system and (2) epistemic uncertainties that arise from the analysts’ incomplete knowledge of the hydrometeorological and site configuration. Principal sources for aleatory variability in LIP flood simulations at NPP sites include precipitation (i.e., magnitude, duration, and temporal distribution of LIP events), initial conditions (e.g., soil moisture content, stormwater drainage discharge, surface storage), boundary conditions (e.g., upstream discharge, downstream water surface elevations), and long-term temporal trends (e.g., climate change). This report describes approaches used for estimating the aleatory variability in precipitation including precipitation-frequency analysis, numerical weather prediction, and stochastic weather generation. Data sources that can be used to estimate aleatory variability in initial and boundary conditions are also described. Approaches to include effects of long-term trends like those from climate change into estimation of aleatory variability are summarized. Sources of epistemic uncertainty in LIP flood simulations at NPP sites include process representation (e.g., multiple approaches to represent runoff generation, stormwater drainage, and hydraulic routing), site configurations (e.g., site layout, flow features, status of temporary flood protection, blockage of drains), model resolution, and long-term temporal trends (e.g., known/planned site alterations, land-use changes at and in the vicinity of the site). The report describes alternative process representations (methods and models), particularly those implemented in the LIP flood simulation software packages reviewed in Task 1 report, and lists the associated model parameters. Approaches for estimating model parameters when surface and subsurface water exchanges occur are also described.

42 ENGINEERING↗

Ensemble flow reconstruction in the atmospheric boundary layer from spatially limited measurements through latent diffusion models

Due to costs and practical constraints, field campaigns in the atmospheric boundary layer typically only measure a fraction of the atmospheric volume of interest. Machine learning techniques have previously successfully reconstructed unobserved regions of flow in canonical fluid mechanics problems and two-dimensional geophysical flows, but these techniques have not yet been demonstrated in the three-dimensional atmospheric boundary layer. Here, we conduct a numerical analogue of a field campaign with spatially limited measurements using large-eddy simulation. We pose flow reconstruction as an inpainting problem, and reconstruct realistic samples of turbulent, three-dimensional flow with the use of a latent diffusion model. The diffusion model generates physically plausible turbulent structures on larger spatial scales, even when input observations cover less than 1% of the volume. Through a combination of qualitative visualization and quantitative assessment, we demonstrate that the diffusion model generates meaningfully diverse samples when conditioned on just one observation. These samples successfully serve as initial conditions for a large-eddy simulation code. We find that diffusion models show promise and potential for other applications for other turbulent flow reconstruction problems.

17 WIND ENERGY↗

A spatially-resolved model of neutron-irradiated tungsten coupling stochastic cluster dynamics and finite deformation plasticity

Structural materials used in nuclear reactors face severe degradation in mechanical properties, such as hardening and embrittlement. At the microscopic scale, this occurs due to creation and accumulation of irradiation-induced defects and their interaction with system dislocations. Although techniques exist which can model evolution of irradiation defects, for instance kinetic transport theory-based models, their interaction with mechanical deformation of the bulk material has not been investigated extensively. In this work, we demonstrate a novel spatially-resolved multiscale coupling between microscopic irradiation defect evolution, modeled using Stochastic Cluster Dynamics (SCD) and macroscopic mechanical deformation modeled using a finite-deformation plasticity model. SCD is used to determine the statistically averaged defect cluster spacing, dependent on operating conditions such as irradiation dose and temperature. This acts as an initial condition that governs the critical resolved shear stress of dislocation glide in the macroscopic plasticity model. This framework is used to predict mechanical behavior in post-mortem test of irradiated Tungsten samples, which has found its importance as structural material used in nuclear reactors. The results obtained using the coupled approach are in good agreement with experimental data of uniaxial tension tests. The model is able to capture the effect of temperature and irradiation dose on the material hardening. Two methods are proposed to estimate hardness – using Tabor's Law relating uniaxial yield stress to hardness and from flat-punch simulations. The results are in reasonable agreement with hardness data from micro-indentation experiments of irradiated Tungsten samples. Finally, the model is also able to reveal microstructural details such as spatial variation in defect density and local stress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

MHD simulation on magnetic compression of field reversed configurations with NIMROD

Abstract Magnetic compression has long been proposed a promising method for plasma heating in a field reversed configuration (FRC). However, it remains a challenge to fully understand the physical mechanisms underlying the compression process, due to its highly dynamic nature beyond the one-dimensional (1D) adiabatic theory model (Spencer et al 1983 Phys. Fluids 26 1564). In this work, magnetohydrodynamics simulations on the magnetic compression of FRCs using the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355) and their comparisons with the 1D theory have been performed. The effects of the assumptions of the theory on the compression process have been explored, and the detailed profiles of the FRC during compression have been investigated. The pressure evolution agrees with the theoretical prediction under various initial conditions. The axial contraction of the FRC can be affected by the initial density profile and the ramping rate of the compression magnetic field, but the theoretical predictions on the FRC’s length in general and the relation r s = 2 r o in particular hold approximately well during the whole compression process, where r s is the major radius of FRC separatrix and r o is that of the magnetic axis. The evolutions of the density and temperature can be affected significantly by the initial equilibrium profile and the ramping rate of the compression magnetic field. During the compression, the major radius of the FRC is another parameter that is susceptible to the ramping rate of the compression field. Basically, for the same magnetic compression ratio, the peak density is higher and the FRC’s radius r s is smaller than the theoretical predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fuzzy Euclidean wormholes in the inflationary universe

In this paper, we investigate complex-valued Euclidean wormholes in the Starobinsky inflation. Due to the properties of the concave inflation potential, the classicality condition at both ends of the wormhole can be satisfied, as long as the initial condition of the inflation field is such that it is located sufficiently close to the hilltop. Furthermore, we compare the probabilities of classicalized wormholes with the Hartle–Hawking compact instantons and conclude that the Euclidean wormholes are probabilistically preferred than compact instantons, if the inflation lasts more than 50 e-foldings. Our result assumes that the Euclidean path integral is the correct effective description of quantum gravity. This opens a new window for various future investigations that can be either confirmed or refuted by future experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct Numerical Simulation of Partial Fuel Stratification Assisted Lean Premixed Combustion for Assessment of Hybrid G-Equation/Well-Stirred Reactor Model

Partial fuel stratification (PFS) is a promising fuel injection strategy to stabilize lean premixed combustion in spark-ignition (SI) engines. PFS creates a locally stratified mixture by injecting a fraction of the fuel, just before spark timing, into the engine cylinder containing homogeneous lean fuel/air mixture. Further, this locally stratified mixture, when ignited, results in complex flame structure and propagation modes similar to partially premixed flames and allows for faster and more stable flame propagation than a homogeneous lean mixture. This study focuses on understanding the detailed flame structures associated with PFS-assisted lean premixed combustion. First, a two-dimensional direct numerical simulation (DNS) is performed using detailed fuel chemistry, experimental pressure trace, and realistic initial conditions mapped from a prior engine large-eddy simulation (LES), replicating practical lean SI operating conditions. DNS results suggest that the conventional triple flame structure is prevalent during the initial stage of flame kernel growth. Both premixed and nonpremixed combustion modes are present with the premixed mode contributing dominantly to the total heat release. Detailed analysis further reveals the effects of flame stretch and fuel pyrolysis on flame displacement speed. Based on the DNS findings, the accuracy of a hybrid G-equation/well-stirred reactor (WSR) combustion model is assessed for the PFS-assisted lean operation in the LES context. It is found that the G-equation model qualitatively captures the premixed branches of the triple flame, while the WSR model predicts the nonpremixed branch of the triple flame. Finally, potential needs for improvements to the hybrid G-equation/WSR modeling approach are discussed.

33 ADVANCED PROPULSION SYSTEMS↗

The Importance of an Ensemble Approach for Modeling Aerosol‐Convection Interactions

Aerosol-convection interactions modulate cloud microphysics, thermodynamics, and updraft intensity, contributing to climate-scale aerosol-radiative forcing. However, quantifying aerosol indirect effects in mixed-phase deep convection remains challenging due to uncertainties in parameterized physics and initial conditions driving nonlinear evolution of convective processes. This study investigates the convective updraft sensitivity to random initial temperature perturbations using an idealized ensemble modeling framework informed by in situ thermodynamic and aerosol observations from the DOE TRACER field campaign. We analyze the impact of small-scale initial perturbations on updraft velocity and supersaturation, and determine the ensemble size required to minimize stochastic internal variability. Results show that minor thermodynamic perturbations can produce updraft variability comparable to aerosol-induced changes reported in prior work. An ensemble of 10 members sufficiently reduces variability, enabling robust investigation of aerosol-related updraft invigoration and informing the design of future ensemble-based aerosol-convection interaction studies in environments with significant spatiotemporal mesoscale thermodynamic and aerosol heterogeneity.

aerosol-convection interactions↗

Impact of volcanic eruptions on CMIP6 decadal predictions: a multi-model analysis

Abstract. In recent decades, three major volcanic eruptions of different intensity have occurred (Mount Agung in 1963, El Chichón in 1982 and Mount Pinatubo in 1991), with reported climate impacts on seasonal to decadal timescales that could have been potentially predicted with accurate and timely estimates of the associated stratospheric aerosol loads. The Decadal Climate Prediction Project component C (DCPP-C) includes a protocol to investigate the impact of volcanic aerosols on the climate experienced during the years that followed those eruptions through the use of decadal predictions. The interest of conducting this exercise with climate predictions is that, thanks to the initialisation, they start from the observed climate conditions at the time of the eruptions, which helps to disentangle the climatic changes due to the initial conditions and internal variability from the volcanic forcing. The protocol consists of repeating the retrospective predictions that are initialised just before the last three major volcanic eruptions but without the inclusion of their volcanic forcing, which are then compared with the baseline predictions to disentangle the simulated volcanic effects upon climate. We present the results from six Coupled Model Intercomparison Project Phase 6 (CMIP6) decadal prediction systems. These systems show strong agreement in predicting the well-known post-volcanic radiative effects following the three eruptions, which induce a long-lasting cooling in the ocean. Furthermore, the multi-model multi-eruption composite is consistent with previous work reporting an acceleration of the Northern Hemisphere polar vortex and the development of El Niño conditions the first year after the eruption, followed by a strengthening of the Atlantic Meridional Overturning Circulation the subsequent years. Our analysis reveals that all these dynamical responses are both model- and eruption-dependent. A novel aspect of this study is that we also assess whether the volcanic forcing improves the realism of the predictions. Comparing the predicted surface temperature anomalies in the two sets of hindcasts (with and without volcanic forcing) with observations we show that, overall, including the volcanic forcing results in better predictions. The volcanic forcing is found to be particularly relevant for reproducing the observed sea surface temperature (SST) variability in the North Atlantic Ocean following the 1991 eruption of Pinatubo.

Bilbao, Roberto (ORCID:0000000307294980)↗