Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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 235 records · Page 13

Developing high performance RF heating scenarios on the WEST tokamak

High power experiments, up to 9.2 MW with LHCD and ICRH, have been carried out in the full tungsten tokamak WEST. Quasi non inductive discharges have been achieved allowing to extend the plasma duration to 53s with stationary conditions in particular with respect to tungsten contamination. Transitions in H mode are obtained lasting up to 4s with weak energy increment at the power crossing the separatrix is close to the threshold. Hot L mode plasmas (Te(0)>3keV) with a confinement time following the ITER L96 scaling are routinely obtained. The weak aspect ratio dependence of this scaling law is confirmed. Tungsten accumulation is generally not an operational issue on WEST. Difficulty of burning through tungsten can prevent from accessing to a hot core plasma in the ramp-up phase or can lead to rapid collapse of the central temperature when radiation is enhanced by a slight decrease of the temperature. Apart few pulses post-boronization, the plasma radiation is rather high (Prad/Ptot~50%) and is dominated by tungsten. This fraction does not vary as the RF power is ramped up and is quite similar in ICRH and/or LHCD heated plasmas. An estimate of the contribution of the RF antennas to the plasma contamination in tungsten is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Planetary and synoptic-scale interactions during the life cycle of a mid-latitude blocking anticyclone over the North Atlantic

The formation of a blocking anticyclone over the North Atlantic has been examined over its entire life-CyCle using the Zwack-Okossi (Z-O) equation as the diagnostic tool. This blocking anticyclone occurred in late October and early November of 1985. The data used were provided by the NASA Goddard Laboratory for Atmospheres on a global 2.O degree latitude by 2.5 degree longitudinal grid. The horizontal distribution of the atmospheric forcing mechanisms that were important to 500 mb block formation, maintenance and decay were examined. A scale-partitioned form of the Z-O equation was then used to examine the relative importance of forcing on the planetary and synoptic scales, and their interactions. As seen in previous studies, the results presented here show that upper tropospheric anticyclonic vorticity advection was the most important contributor to block formation and maintenance. However, adiabatic warming, and vorticity tilting were also important at various times during the block lifetime. In association with precursor surface cyclogenesis, the 300 mb jet streak in the downstream (upstream) from a long-wave trough (ridge) amplified significantly. This strengthening of the jet streak enhanced the anti-cyclonic vorticity advection field that aided the amplification of a 500 mb short-wave ridge. Tile partitioned height tendency results demonstrate that the interactions between the planetary and sn,noptic-scale through vorticity advection was the most important contributor to block formation. Planetary-scale, synoptic-scale. and their interactions contributed weakly to the maintenance of the blocking anticyclone with the advection of synoptic-scale vorticity by the planetary-scale flow playing a more important role. Planetary-scale decay ofthe long-wave ridge contributed to the demise of this blocking event.

Smith, Phillip J.↗

Revealing atomic-scale switching pathways in van der Waals ferroelectrics

Two-dimensional (2D) van der Waals (vdW) materials hold the potential for ultrascaled ferroelectric (FE) devices due to their silicon compatibility and robust polarization down to atomic scale. However, the inherently weak vdW interactions enable facile sliding between layers, introducing complexities beyond those encountered in conventional ferroelectric materials and presenting substantial challenges in uncovering intricate switching pathways. Here, we combine atomic-resolution imaging under in situ electrical biasing conditions with first-principles calculations to unravel the atomic-scale switching mechanisms in SnSe, a vdW group IV monochalcogenide. Our results uncover the coexistence of a consecutive 90° switching pathway and a direct 180° switching pathway from antiferroelectric (AFE) to FE order in this vdW system. Atomic-scale investigations and strain analysis reveal that the switching processes simultaneously induce interlayer sliding and compressive strain, while the lattice remains coherent despite the presence of multidomain structures. These findings elucidate vdW ferroelectric switching dynamics at atomic scale and lay the foundation for the rational design of 2D ferroelectric nanodevices.

Li, Xinyan [Rice University, Houston, TX (United S↗

Evidence for Suppression of Structure Growth in the Concordance Cosmological Model

We present evidence for a suppressed growth rate of large-scale structure during the dark-energy dominated era. Modeling the growth rate of perturbations with the ``growth index'' γ, we find that current cosmological data strongly prefer a higher growth index than the value γ=0.55 predicted by general relativity in a flat ΛCDM cosmology. Both the cosmic microwave background data from Planck and the large-scale structure data from weak lensing, galaxy clustering, and cosmic velocities separately favor growth suppression. When combined, they yield γ = 0.633$^{+0.025}_{–0.024}$, excluding γ = 0.55 at a statistical significance of 3.7σ. The combination of fσ8 and Planck measurements prefers an even higher growth index of γ = 0.639$^{+0.024}_{–0.025}$, corresponding to a 4.2σ-tension with the concordance model. In Planck data, the suppressed growth rate offsets the preference for nonzero curvature and fits the data equally well as the latter model. Furthermore, a higher γ leads to a higher matter fluctuation amplitude S 8 inferred from galaxy clustering and weak lensing measurements, and a lower S 8 from Planck data, effectively resolving the S 8 tension.

79 ASTRONOMY AND ASTROPHYSICS↗

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The gravitational lensing imprints of DES Y3 superstructures on the CMB: a matched filtering approach

Low-density cosmic voids gravitationally lens the cosmic microwave background (CMB), leaving a negative imprint on the CMB convergence |$\kappa$|⁠. This effect provides insight into the distribution of matter within voids, and can also be used to study the growth of structure. We measure this lensing imprint by cross-correlating the Planck CMB lensing convergence map with voids identified in the Dark Energy Survey Year 3 (DES Y3) data set, covering approximately 4200 deg|$^2$| of the sky. We use two distinct void-finding algorithms: a 2D void-finder that operates on the projected galaxy density field in thin redshift shells, and a new code, Voxel, which operates on the full 3D map of galaxy positions. We employ an optimal matched filtering method for cross-correlation, using the Marenostrum Institut de Ciències de l’Espai N-body simulation both to establish the template for the matched filter and to calibrate detection significances. Using the DES Y3 photometric luminous red galaxy sample, we measure |$A_\kappa$|⁠, the amplitude of the observed lensing signal relative to the simulation template, obtaining |$A_\kappa = 1.03 \pm 0.22$| (⁠|$4.6\sigma$| significance) for Voxel and |$A_\kappa = 1.02 \pm 0.17$| (⁠|$5.9\sigma$| significance) for 2D voids, both consistent with Lambda cold dark matter expectations. We additionally invert the 2D void-finding process to identify superclusters in the projected density field, for which we measure |$A_\kappa = 0.87 \pm 0.15$| (⁠|$5.9\sigma$| significance). The leading source of noise in our measurements is Planck noise, implying that data from the Atacama Cosmology Telescope, South Pole Telescope and CMB-S4 will increase sensitivity and allow for more precise measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

The Universe at Moderate Redshift

The report covers the work done in the past year and a wide range of fields including properties of clusters of galaxies; topological properties of galaxy distributions in terms of galaxy types; patterns of gravitational nonlinear clustering process; development of a ray tracing algorithm to study the gravitational lensing phenomenon by galaxies, clusters and large-scale structure, one of whose applications being the effects of weak gravitational lensing by large-scale structure on the determination of q(0); the origin of magnetic fields on the galactic and cluster scales; the topological properties of Ly(alpha) clouds the Ly(alpha) optical depth distribution; clustering properties of Ly(alpha) clouds; and a determination (lower bound) of Omega(b) based on the observed Ly(alpha) forest flux distribution. In the coming year, we plan to continue the investigation of Ly(alpha) clouds using larger dynamic range (about a factor of two) and better simulations (with more input physics included) than what we have now. We will study the properties of galaxies on 1 - 100h(sup -1) Mpc scales using our state-of-the-art large scale galaxy formation simulations of various cosmological models, which will have a resolution about a factor of 5 (in each dimension) better than our current, best simulations. We will plan to study the properties of X-ray clusters using unprecedented, very high dynamic range (20,000) simulations which will enable us to resolve the cores of clusters while keeping the simulation volume sufficiently large to ensure a statistically fair sample of the objects of interest. The details of the last year's works are now described.

Cen, Renyue↗

Quarkonium semiclassical transport in quark-gluon plasma: factorization and quantum correction

We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Zhang, Hai↗

Massive superclusters as a probe of the nature and amplitude of primordial density fluctuations

It is pointed out that correlation studies of galaxy positions have been widely used in the search for information about the large-scale matter distribution. The study of rare condensations on large scales provides an approach to extend the existing knowledge of large-scale structure into the weakly clustered regime. Shane (1975) provides a description of several apparent massive condensations within the Shane-Wirtanen catalog, taking into account the Serpens-Virgo cloud and the Corona cloud. In the present study, a description is given of a model for estimating the frequency of condensations which evolve from initially Gaussian fluctuations. This model is applied to the Corona cloud to estimate its 'rareness' and thereby estimate the rms density contrast on this mass scale. An attempt is made to find a conflict between the density fluctuations derived from the Corona cloud and independent constraints. A comparison is conducted of the estimate and the density fluctuations predicted to arise in a universe dominated by cold dark matter.

Kaiser, N.↗

QCD analysis of Δ S = 0 hadronic parity violation

We present a QCD analysis of the effective weak Hamiltonian at hadronic energy scales for strangeness-nonchanging (ΔS = 0) hadronic processes. Performing a leading-order renormalization group analysis in QCD from the W to the O(2 GeV) energy scale, we derive the pertinent effective Hamiltonian for hadronic parity violation, including the effects of both neutral and charged weak currents. We compute the complete renormalization group evolution of all isosectors and the evolution through heavy-flavor thresholds for the first time. We show that the additional four-quark operators that enter below the W mass scale from QCD operator mixing effects form a closed set, and they result in a 12 x 12 anomalous dimension matrix. Computing the resulting effective Hamiltonian and comparing to earlier results, we affirm the importance of operator mixing effects and find, as an example, that the parity-violating pion-nucleon coupling constant, using the factorization Ansatz and an assessment of the pertinent quark charge of the nucleon in lattice QCD at the 2 GeV scale, is in better agreement with recent experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗

Constraining the baryonic feedback with cosmic shear using the DES Year-3 small-scale measurements

ABSTRACT We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package baccoemu to accelerate the evaluation of the baryonic non-linear matter power spectrum. We design our analysis pipeline to focus on the constraints of the baryonic suppression effects, utilizing the implication given by a principal component analysis on the Fisher forecasts. Our constraint on the baryonic effects can then be used to better model and ameliorate the effects of baryons in producing cosmological constraints from the next-generation large-scale structure surveys. We detect the baryonic suppression on the cosmic shear measurements with a ∼2σ significance. The characteristic halo mass for which half of the gas is ejected by baryonic feedback is constrained to be $M_c \gt 10^{13.2} \, h^{-1} \, \mathrm{M}_{\odot }$ (95 per cent C.L.). The best-fitting baryonic suppression is $\sim 5{{\ \rm per\ cent}}$ at $k=1.0 \, {\rm Mpc}\ h^{-1}$ and $\sim 15{{\ \rm per\ cent}}$ at $k=5.0 \, {\rm Mpc} \ h^{-1}$. Our findings are robust with respect to the assumptions about the cosmological parameters, specifics of the baryonic model, and intrinsic alignments.

79 ASTRONOMY AND ASTROPHYSICS↗

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Urban Land Surface Effects on Summertime Clouds and Moist Convection in Houston Under Different Synoptic Conditions

Urban landscapes modify cloud formation and convection through complex thermodynamic and aerodynamic processes; however, their influence under different synoptic regimes remains poorly understood. This study investigates the impact of the Houston metropolitan area on summertime cloud cover and convective cell characteristics using a combination of satellite observations, radar data, and high-resolution process-based modeling. We isolate urban effects by comparing model simulations with realistic urban land cover against hypothetical scenarios where all urban areas are replaced by rural vegetation. Results reveal that Houston's land cover consistently enhances cloud fraction and convective activity relative to surrounding rural areas, altered by large-scale meteorological forcing. Under weakly forced conditions, enhanced surface heat flux primarily contributes to driving low-level convergence and vertical ascent, leading to over 8% cloud fraction increase between 2 and 6 km over the city. Under strongly forced conditions, urban influences manifest differently depending on cloud type. For non-convective clouds, the city acts as a barrier that decelerates and lifts moist southerly inflow, increasing low cloud cover over the urban core, while decreasing it downwind the city. For convective clouds, both synoptic ascent and urbanization modulate moisture redistribution and cloud structure, producing modest cloud enhancement over the city and slight suppression over the downwind area. Urbanization exerts small changes in the intensity of the convective cells; however, it significantly decreases their duration and traveling distance. This work highlights the importance of accounting for land surface heterogeneity in modeling clouds and precipitation and demonstrates that urban impacts on clouds are highly regime-dependent.

Liu, Ye [Pacific Northwest National Laboratory (PN↗

Magnetic reconnection in the era of exascale computing and multiscale experiments

Astrophysical plasmas have the remarkable ability to preserve magnetic topology, which inevitably gives rise to the accumulation of magnetic energy within stressed regions including current sheets. This stored energy is often released explosively through the process of magnetic reconnection, which produces a reconfiguration of the magnetic field, along with high-speed flows, thermal heating and nonthermal particle acceleration. Either collisional or kinetic dissipation mechanisms are required to overcome the topological constraints, both of which have been predicted by theory and validated with in situ spacecraft observations or laboratory experiments. However, major challenges remain in understanding magnetic reconnection in large systems, such as the solar corona, where the collisionality is weak and the kinetic scales are vanishingly small in comparison with macroscopic scales. The plasmoid instability or formation of multiple plasmoids in long, reconnecting current sheets is one possible multiscale solution for bridging this vast range of scales, and new laboratory experiments are poised to study these regimes. In conjunction with these efforts, we anticipate that the coming era of exascale computing, together with the next generation of observational capabilities, will enable new progress on a range of challenging problems, including the energy build-up and onset of reconnection, partially ionized regimes, the influence of magnetic turbulence and particle acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

Large ordered moment with strong easy-plane anisotropy and vortex-domain pattern in the kagome ferromagnet Fe 3 Sn

We report the magnetic anisotropy of kagome bilayer ferromagnet Fe 3 Sn probed by the bulk magnetometry and magnetic force microscopy (MFM) on high-quality single crystals. The dependence of magnetization on the orientation of the external magnetic field reveals strong easyplane magnetocrystalline anisotropy and anisotropy of the saturation magnetization. The leading magnetocrystalline anisotropy constant shows a monotonous increase from $\mathrm{K_1$≈ –1.0 \times 10^6 J/m^3}$ at 300 K to $\mathrm{–1.3 \times 10^6 J/m^3}$ at 2 K. Our ab initio electronic structure calculations yield the value of total magnetic moment of 7.1$μ_Β$/$f.u.$ and a magnetocrystalline anisotropy energy density of –0.57 $\mathrm{meV/}$$f.u.$ $\mathrm{(–1.62 \times 10^6J/m^3)}$ both being in reasonable agreement with the experimental values. The MFM imaging reveals micrometer-scale magnetic vortices with weakly pinned cores that vanish at the saturation field of ~3T applied perpendicular to the kagome plane. The observed vortex-domain structure is well reproduced by the micromagnetic simulations, using the experimentally determined value of the anisotropy and exchange stiffness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Partitioning uncertainty in projections of Arctic sea ice

Abstract Improved knowledge of the contributing sources of uncertainty in projections of Arctic sea ice over the 21st century is essential for evaluating impacts of a changing Arctic environment. Here, we consider the role of internal variability, model structure and emissions scenario in projections of Arctic sea-ice area (SIA) by using six single model initial-condition large ensembles and a suite of models participating in Phase 5 of the Coupled Model Intercomparison Project. For projections of September Arctic SIA change, internal variability accounts for as much as 40%–60% of the total uncertainty in the next decade, while emissions scenario dominates uncertainty toward the end of the century. Model structure accounts for 60%–70% of the total uncertainty by mid-century and declines to 30% at the end of the 21st century in the summer months. For projections of wintertime Arctic SIA change, internal variability contributes as much as 50%–60% of the total uncertainty in the next decade and impacts total uncertainty at longer lead times when compared to the summertime. In winter, there exists a considerable scenario dependence of model uncertainty with relatively larger model uncertainty under strong forcing compared to weak forcing. At regional scales, the contribution of internal variability can vary widely and strongly depends on the calendar month and region. For wintertime SIA change in the Greenland-Iceland-Norwegian and Barents Seas, internal variability contributes 60%–70% to the total uncertainty over the coming decades and remains important much longer than in other regions. We further find that the relative contribution of internal variability to total uncertainty is state-dependent and increases as sea ice volume declines. These results demonstrate that internal variability is a significant source of uncertainty in projections of Arctic sea ice.

Bonan, David B. (ORCID:0000000338676009)↗