Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ocean Modeling Climate Modeling”

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 217 records · Page 12

Evaluation of Ice Nucleating Particles and Their Sources in the Central Arctic during MOSAiC

The Arctic is warming faster than any other region on Earth, causing glaciers to melt, frozen ground to thaw, and the Arctic Ocean's ice to shrink. These consequential changes induce feedback loops that exacerbate warming and affect weather and climate worldwide. Microscopic airborne particles called "aerosols" and clouds in the sky are crucial for regulating heat and light reaching the Arctic surface. However, the magnitude of their effects is not adequately quantified, especially in the central Arctic, where they impact temperatures directly over the sea ice. Unique aerosols called "ice nucleating particles" (INPs), which play a significant role in cloud ice production, remain understudied. Understanding how ice forms in clouds is critical in the Arctic, as it affects cloud lifespan, interactions with heat and light, and precipitation. In this project, we conducted the first-ever observations of INPs in the middle of the Arctic over a whole year, covering the entire period when sea ice grows and melts. Furthermore, these are the first observations of INPs in different size ranges throughout the year, anywhere in the world. We use DNA sequencing to evaluate the presence of various types of microorganisms, while INP measurements on seawater, sea ice, snow, and meltwater samples help assess potential local Arctic sources of INPs. The results from this work are currently being used to improve the accuracy of Arctic cloud formation in various models.

54 ENVIRONMENTAL SCIENCES↗

Global Precipitation Experiment—A New World Climate Research Programme Lighthouse Activity

The future state of the global water cycle and the prediction of freshwater availability for humans around the world remain among the challenges of climate research and are relevant to several United Nations Sustainable Development Goals. The Global Precipitation Experiment (GPEX) takes on the challenge of improving the prediction of precipitation quantity, phase, timing, and intensity, characteristics that are products of a complex integrated system. It will achieve this by leveraging existing World Climate Research Programme (WCRP) activities and community capabilities in satellite, surface-based, and airborne observations, modeling, and experimental research and by conducting new and focused activities. It was launched in October 2023 as a WCRP Lighthouse Activity. Here, we present an overview of the GPEX Science Plan that articulates the primary science questions related to precipitation measurements, process understanding, model performance and improvements, and plans for capacity development. The central phase of GPEX is the WCRP Years of Precipitation for 2–3 years with coordinated global field campaigns focusing on different storm types (atmospheric rivers, mesoscale convective systems, monsoons, and tropical cyclones, among others) over different regions and seasons. Activities are planned over the three phases (before, during, and after the Years of Precipitation) spanning a decade. These include gridded data evaluation and development, advanced modeling, enhanced understanding of processes critical to precipitation, multiscale prediction of precipitation events across scales, and capacity development. These activities will be further developed as part of the GPEX Implementation Plan.

Climate change↗

Comparative Evaluation of the Ability of the MYNN‐EDMF PBL Scheme in WRF Model to Reproduce Near Surface Wind Speed Over Different Topographical Types

Abstract This study systematically evaluates the performance of the Mellor‐Yamada‐Nakanishi‐Niino‐Eddy‐Diffusion‐Mass‐Flux planetary boundary layer (PBL) scheme within the Weather Research and Forecasting (WRF) model in simulating near‐surface wind speeds across various topographies in New York State (NYS). Simulated wind speeds are compared with in‐situ measurements from 22 surface sites, grouped into six topographic categories: continental plain (CT), lakeside (LS), river valley (RV), Long Island (LI), Block Island (BI), and offshore ocean (OO). A quantitative evaluation based on Relative Euclidean Distance shows that wind speeds at the OO site are the most accurately reproduced, followed by those at LI sites, while the model performs less accurately for the remaining topographic groups. Wind speeds over CT sites tend to be overestimated by approximately 1 m/s, although their diurnal variability (DV) is well captured. In contrast, the model underestimates wind DV at LS, RV, LI, and BI sites, with the largest biases occurring at LI and BI, resulting in underestimated daytime wind speed and/or overestimated nighttime wind speed. The OO winds exhibit minimal diurnal variation, accurately captured by our WRF model. The surface wind diurnal variation is closely linked to PBL development. Among the indicators of PBL development, surface potential temperature biases most strongly correlate with wind speed biases. Our WRF model faces challenges in capturing the distinctions between winds influenced by local circulations and those over continental plains, and the significantly stronger winds at OO compared to BI. Potential causes for these biases are discussed, offering pathways for improving surface wind simulations in future.

54 ENVIRONMENTAL SCIENCES↗

Subsurface Cooling and Sea Surface Temperature Pattern Formation Over the Equatorial Pacific

Abstract The equatorial cold tongue region has not warmed up in response to historical radiative forcing in the real world, contrary to the strong warming often simulated by climate models. Here we demonstrate that climate models fail to represent one or both of the key processes driving observed sea surface temperature (SST) pattern formation: a realistic surface wind stress pattern shaping subsurface cooling through wind‐driven circulation changes, and effective connectivity between subsurface and surface temperatures via upwelling and mixing. Consequently, none of the models approximate the observed lack of cold tongue SST warming and strengthening of zonal SST gradient across the equatorial Pacific. Furthermore, those that come closest achieve this due to interhemispheric warming differences rather than equatorial dynamics as observed. Addressing different origins of subsurface cooling in observations and simulations, and how they connect to SST, will lead to improved understanding of tropical Pacific SST changes to date and how they will evolve in the future.

Jiang, Feng [Lamont‐Doherty Earth Observatory Colu↗

Air‐Ice‐Ocean Coupling During a Strong Mid‐Winter Cyclone: Observing Coupled Dynamic Interactions Across Scales

Abstract Arctic cyclones are key drivers of sea ice and ocean variability. During the 2019–2020 Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, joint observations of the coupled air‐ice‐ocean system were collected at multiple spatial scales. Here, we present observations of a strong mid‐winter cyclone that impacted the MOSAiC site as it drifted in the central Arctic pack ice. The sea ice dynamical response showed spatial structure at the scale of the evolving and translating cyclonic wind field. Internal ice stress and ocean stress play significant roles, resulting in timing offsets between the atmospheric forcing and the ice response and post‐cyclone inertial ringing in the ice and ocean. Ice motion in response to the wind field then forces the upper ocean currents through frictional drag. The strongest impacts to the sea ice and ocean from the passing cyclone occur as a result of the surface impacts of a strong atmospheric low‐level jet (LLJ) behind the trailing cold front and changing wind directions between the warm‐sector LLJ and post cold‐frontal LLJ. Impacts of the cyclone are prolonged through the coupled ice‐ocean inertial response. Local impacts of the approximately 120 km wide LLJ occur over a 12 hr period or less and at scales of a kilometer to a few tens of kilometers, meaning that these impacts occur at combined smaller spatial scales and faster time scales than most satellite observations and coupled Earth system models can resolve.

Meteorology & Atmospheric Sciences↗

Eastward-shifting boreal summer intraseasonal oscillation amplifies North America heatwave and wildfire risks in warming climate

Boreal summer intraseasonal oscillation (BSISO) is a key component of tropical climate variability, characterized by eastward and northeastward propagation of organized convection across the Indo-Pacific region. BSISO influences weather and climate extremes through atmospheric teleconnections, but its future changes and related extremes remain unclear. Here, based on Coupled Model Intercomparison Project Phase 6, we find that under a high-emission scenario, BSISO convection will shift eastward by ~3° (~300 km) in 2065–2099, driven by moisture profile changes associated with sea surface temperature warming over central-to-eastern Pacific. This eastward-shifted BSISO convection, along with a northward expansion of the westerly jet, strengthens the BSISO teleconnections that extend into broader areas of North America. These changes will increase the BSISO-related heatwave risks by 23% and wildfire risks by 3.5 times over northern North America compared to present-day conditions. The findings underscore the amplification of BSISO’s influence on extreme risks of extreme weather, emphasizing the need for improved adaptation and mitigation strategies for future climates.

BSISO teleconnection↗

The Influence of Anomalous Biomass Emissions on ENSO in CESM2

Abstract The influence of biomass burning (BB) aerosols arising from wildfires and agricultural fires on the transient coupled evolution of El Niño–Southern Oscillation (ENSO) is explored in Community Earth System Model, version 2 (CESM2). For both El Niño and La Niña, two 20-member ensembles are generated from initial states that are predisposed to evolve into ENSO events. For each ENSO phase, one ensemble is forced with the observed BB emissions during satellite-era ENSO events while the other is forced with a climatological annual cycle, with the responses to anomalous BB emissions estimated from interensemble differences. It is found that the regional responses to anomalous BB emissions occur mainly during boreal fall, which is also the time of the climatological seasonal maximum in emissions. Transient responses are identified in precipitation, clouds, and radiation in both the tropics and extratropics. At the onset of El Niño, these include increased precipitation in the northern branch of the intertropical convergence zone (ITCZ) and an enhancement of cloud albedo and amount across the Maritime Continent and eastern subtropical Pacific Ocean. Additional responses are identified through the course of El Niño and successive La Niña events, the net effect of which is to strengthen sea surface temperature (SST) anomalies in the eastern Pacific Ocean during El Niño and warm the tropical Pacific Ocean during La Niña. These responses improve the simulation of ENSO power, diversity, and asymmetry in CESM2. Significance Statement Biomass burning emissions from both wildfires and agricultural fires during ENSO act as a climate feedback, as they are both driven by fire weather anomalies tied to ENSO and influence clouds, radiation, and precipitation on a global scale. While the largest biomass burning (BB) emissions emanate mainly from Indonesia, other regions also exhibit systematic responses. When simulated in CESM2, the emissions collectively impact the flows of energy across the tropical Pacific Ocean in a spatially and temporally complex fashion, with the net effect of improving CESM2’s simulation of ENSO events.

Fasullo, John T. [National Center for Atmospheric ↗

CMIP7 Data Request: atmosphere priorities and opportunities

This paper presents a comprehensive overview of the Coupled Model Intercomparison Project Phase 7 (CMIP7) request for data unlocking key research avenues in atmospheric science and provides justification for the resources needed to produce this data. Topics within the CMIP7 Atmosphere Theme centre around processes and feedbacks in atmospheric science such as clouds, aerosols and atmospheric chemistry, atmospheric circulation, temperature variability and extremes, radiative forcings, and Earth system model evaluation. These topics are summarised in this paper as scientific “opportunities” which will be realised through CMIP7 experiments and Earth system model outputs. These opportunities were submitted by a thematic group of atmospheric science community representatives combined with an extended consultation process. The production of these variables will close key gaps and uncertainties identified during previous rounds of CMIP, and will be broadly used by scientific, policy, governmental, industry, and other communities that rely on climate model projections for research and decision making, including supporting the 7th Intergovernmental Panel on Climate Change Assessment Report (AR7). As an author group, we also reflect on the process used to collate this data request and make recommendations to future CMIP governance on implementing a consultation on this scale in the future.

58 GEOSCIENCES↗

Oceanic fronts shape hemispheric contrasts in polar stratospheric extremes

Sudden Stratospheric Warmings (SSWs) and Polar Stratospheric Clouds (PSCs) exhibit striking inter-hemispheric asymmetries: SSWs are frequent in the Arctic but rare in the Antarctic, while PSCs are more persistent in the Antarctic. Although land–sea thermal contrast and orography (LSCO) have been traditionally invoked to explain these asymmetries, here we show using semi-idealized model experiments that even though LSCO strongly impacts the mean state of SSWs and PSCs, it alone cannot fully account for the observed differences. Using model experiments, we reveal that midlatitude oceanic sea surface temperature (SST) fronts represent a crucial supporting additional driver of the hemispheric stratospheric differences. Like LSCO, SST fronts enhance stratospheric convergence of resolved waves, strengthening the Brewer–Dobson Circulation and inducing high-latitude adiabatic warming. This warming significantly enhances Arctic SSW frequency and strongly suppresses PSC formation. Sub-grid-scale non–resolved wave forcing modulates the stratospheric resolved waves effect. The oceanic impact is dominated by North Pacific SST fronts. Our results highlight the indispensable role of SST fronts in shaping Arctic–Antarctic asymmetries in stratospheric dynamics and associated extremes.

Omrani, Nour-Eddine [Univ. of Bergen (Norway)] (OR↗

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗

Biophysical model of eelgrass and water quality in Coos Bay, OR shows greater mitigation potential for ocean acidification than hypoxia

Seagrass beds provide important ecosystem services and are valued, in part, for their potential to mediate stressors such as ocean acidification and hypoxia (OAH) for sensitive species. However, the susceptibility of seagrasses to anthropogenic impacts and recent declines motivate the need to better understand the drivers of seagrass and the water quality consequences that occur with variation in seagrass abundance. To meet this need, we leveraged existing monitoring data (water quality and seagrass), hydrodynamic circulation model, and biogeochemical model framework with seagrass submodel, to produce a biophysical model of Coos Bay estuary, Oregon, U.S. The model includes biogeochemical processes involving water quality, plankton, seagrass, and sediment-water interactions. Ecosystem models like this are useful for evaluating complex estuarine systems because they allow us to extend our understanding of system dynamics beyond existing observations and perform experiments to identify the processes driving observed patterns. We used the biophysical model of Coos Bay to evaluate the dynamics of water quality and native eelgrass (Zostera marina) under three eelgrass abundance scenarios (zero eelgrass, current extent, and maximum observed extent) to elucidate the relationship between eelgrass and OAH. Including eelgrass in the Coos Bay model produced results that more closely resembled water quality observations - dissolved oxygen (DO) and pH were more dynamic in simulations with eelgrass, often having both higher highs and lower lows. While there were some areas of the estuary where DO improved with the addition of eelgrass to the model there was overall a small net increase in harmful DO conditions (based on a salmon physiological threshold). In contrast, ocean acidification conditions, pH and calcium carbonate saturation state for aragonite (Ω), were improved (based on oyster requirements) with the addition of eelgrass - although the magnitude of improvement differed seasonally and spatially. Our new model represents a useful tool - one which accounts for and controls the relevant physical and biogeochemical processes - to evaluate conditions that confer resilience or enhance vulnerability to OAH in an important Pacific Northwest coastal estuary and results can inform the OAH-related dynamics occurring in other eastern boundary current estuaries.

FVCOM-ICM↗

Detecting tropospheric composition and climate responses to US air pollution controls in the context of internally-arising variability

Since the 1970s, air pollutant emissions controls in the United States (US) have lowered concentrations of ozone (O 3 ) and aerosols, which have opposing radiative effects on surface temperature. Using a pair of initial-condition ensembles generated by a fully-coupled chemistry-climate model, we simulate the “world avoided” by US air pollution controls. In this counterfactual world, we find tropospheric column O 3 increases, robust to natural internal variability, extending across the Northern Hemisphere. Robust aerosol increases, dominated by sulfate, remain localized near the US. Wintertime Northwest Atlantic cloud droplet number concentration is particularly sensitive to US aerosol. While an ensemble mean US surface cooling signal (−0.4 °C) implies that aerosol-driven cooling prevails over any O 3 -induced warming, we find that large regional internal variability will confound its detection in any single transient realization. Larger signal-to-noise ratios for composition versus climate variables underscore the greater detectability of emissions-driven changes in tropospheric composition compared to their associated climate impacts.

54 ENVIRONMENTAL SCIENCES↗

A Comparative Study of Cloud Properties Between Northern and Southern High Latitudes Based on ARM Observations and EAMv2 Simulations

Earth's high-latitude regions show a considerable amount of cloud coverage, which significantly affects global radiation budgets. Here, in this study, we examine ground- and ship-based remote sensing observations from two Atmospheric Radiation Measurement (ARM) field campaigns conducted in the high-latitude regions of Northern and Southern Hemispheres (NH and SH). A comparison of macrophysical properties for low, mid-level, high, and deep clouds shows that the cloud top heights and thickness values are higher in the Arctic for each cloud type compared with the Southern Ocean and Antarctic regions, possibly due to higher sea surface temperatures in the north Atlantic. The Arctic also shows lower ice phase occurrence frequencies compared with the SH high latitudes for low clouds. For high, mid-level and deep clouds, the cloud phase distributions as a function of temperature show small geographical variabilities. The observations are then compared with the Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) to evaluate the representations of these cloud properties. Here, the results show that the model captures cloud base and top heights well and reproduces the differences in cloud macrophysical properties between the two hemispheres. But the model underestimates ice phase frequencies for all cloud types. Such model cloud phase bias is strongly correlated with positive relative humidity biases and negative aerosol number concentration biases, indicating an important role of relative humidity and aerosols in modulating cloud phase partitioning in the model.

54 ENVIRONMENTAL SCIENCES↗

Thermodynamics-guided machine learning model for predicting convective boundary layer height and its multi-site applicability

Accurate estimation of convective boundary layer height (CBLH) is vital for weather, climate, and air quality modeling. Machine learning (ML) shows promise in CBLH prediction, but input parameter selection often lacks physical grounding, limiting generalizability. This study introduces a novel ML framework for CBLH prediction, integrating thermodynamic constraints and the diurnal CBLH cycle as an implicit physical guide. Boundary layer growth is modeled as driven by surface heat fluxes and atmospheric heat absorption represented with the low tropospheric stability, using the diurnal cycle as input and output. TPOT and AutoKeras are employed to select optimal models, validated against Doppler lidar-derived CBLH data, achieving an R 2 of 0.84 across untrained years. Comparisons of eddy covariance (ECOR) and energy balance Bowen ratio (EBBR) flux measurements show the same prediction capability. Models trained on the ARM SGP C1 site with ECOR data and tested at E37 and E39 yield R 2 values of 0.79 and 0.81, respectively, demonstrating their adaptability. The ML model trained with all sites' data slightly enhances the performance compared with ML models trained over single-site data. The interquartile range for predicted CBLH is consistently narrower than that for DL-derived CBLH, reflecting lower variability in predicted CBLH compared to DL-derived CBLH, which is influenced by additional factors, which are not well represented with the model inputs. The model's generalizability across multiple sites at the ARM SGP site demonstrates its potential for transfer to greater distances, offering a scalable approach for enhancing boundary layer parameterization in atmospheric models.

Chu, Yufei [Stony Brook Univ., NY (United States)]↗

WEC-Sim Modeling of Laminar Scientific's Patented Seesaw Wave Energy Converter: Preprint

Laminar Scientific's patented seesaw wave energy converter was modeled in WEC-Sim to predict performance. The device operates by utilizing ocean surface waves to rotate a truss in pitch about a pivot. The pivot is located at the top of two pylons, which are embedded in the seafloor. The seesaw has a float on either end, and the buoyancy forces from each float cause the system to rise or fall with passing waves. Device performance relies upon seesaw length and ocean wavelength creating an antiphase effect. The seesaw truss has an adjustable length intended to achieve this effect. The operation method enforces a narrow band of wavelengths which induce the largest rotational motion from the device. The hydrodynamic analysis of the device was performed using Capytaine, and the results confirmed that the device operates best in a narrow frequency band. Four float-to-float spacing cases and three pylon radii were examined. The hydrodynamic results indicate a match between the model and the physical expectations for the device, and that varying the pylon radii by 0.1-m increments for three instances creates minimal changes in hydrodynamic properties. Power matrices for three float spacing cases of the device were simulated with Joint North Sea Wave Project spectra waves and optimal power take-off damping in WEC-Sim. The maximum average power production for the 15-m spacing case was 14.1 kW with a 5.0-s peak wave period and 4-m significant wave heights. Plots of capture-width ratios indicated that the device performance was linear and confirmed that the device is optimal in a narrow frequency band. The maximum percentage of the available wave power produced by the 15-m device was approximately 16%. Simulations of the device in regular waves were used to produce plots of average power compared to a ratio of float spacing to wavelength. These plots indicate that the power production is maximized at a ratio of 0.5, and further confirm that the device has a narrow frequency response. The device was simulated at an example field location, where the device produced an annual average power rating of 1.6 kW given an average omnidirectional wave climate of 10.3 kW m-1 and an optimal, linearized power take-off model. While the maximum predicted device performance is reliant upon a narrow band of wave frequencies, the conducted analysis provides an opportunity to improve device design prior to prototyping and testing. Modifying the design to respond to a broader frequency range would improve device performance.

numerical modeling↗

Impact of Turbulence Representation on the Relationship Between Cloud Feedback and Aerosol‐Cloud Interaction in an E3SMv2 Perturbed Parameter Ensemble

Recent studies reveal an anti‐correlation between global cloud feedback (CF) and effective radiative forcing due to aerosol‐cloud interaction (ERFaci) in Earth system models, but the physical mechanisms underlying it remain uncertain. Here we investigate how different turbulence representations contribute to this relationship over the global ocean using an ensemble of Energy Exascale Earth System Model version 2 simulations with perturbed turbulence parameters. The anti‐correlation appears only in the tropical ascent regime. In the Northern Hemisphere midlatitude and high latitude regimes, there is no significant correlation, and in the tropical marine low cloud and Southern Ocean regimes, the correlation is positive. These opposite correlations are primarily driven by opposing CF responses to perturbed parameters. We find that the mean‐state turbulent mixing strength affects both CF and ERFaci, enabling strong correlations in certain regimes. This study highlights the complex linkages between CF and ERFaci through turbulent processes across diverse cloud regimes.

54 ENVIRONMENTAL SCIENCES↗

Multi-model analysis of the impact of water vapor on the radiative forcing of volcanic aerosols after the 2022 Hunga Eruption

On 15 January 2022, the Hunga volcano eruption released unprecedented amounts of water vapor into the atmosphere alongside a modest amount of SO 2 . In this work we analyse results from multiple Earth system models as part of the Hunga Tonga-Hunga Ha'apai Volcano Impact Model Observation Comparison Project. Our results show a good model agreement over the climatic outcomes of the eruption, overall indicating a significant negative radiative forcing from the Hunga eruption. The multi-model mean of global instantaneous radiative forcing averaged over 2022–2023 is estimated at −0.19 ± 0.06 W m −2 at the top-of-atmosphere (TOA), and −0.16 ± 0.06 W m −2 at the surface. Simulations with free-running meteorology and climatological sea surface temperatures and sea ice yield a global mean TOA forcing of −0.14 ± 0.10 W m −2 across two models for the first 2 years, decreasing to −0.09 ± 0.10 W m −2 on average between 2022 and 2027. However, these global values may be underestimated by about 50 %, considering that recent SO 2 injection retrievals suggest nearly twice the amount than the 0.5 Tg-SO 2 used in the protocol. We also find that the contribution from added stratospheric water vapor is minimal and that the injected SO 2 and the resulting formation of stratospheric sulfate dominate the radiative forcing. However, water vapor played a key role in the initial aerosol growth, leading to a stronger negative radiative forcing during the first 6 months after the eruption compared to simulations without water vapor co-injection.

Geosciences↗

The Role of Internal Variability in Springtime Arctic Amplification from 1980 to 2022

Arctic amplification (AA) refers to the enhanced warming of the Arctic relative to the global average due to rising greenhouse gases, measured as the ratio of Arctic-mean to global-mean surface air temperature (SAT) trends. From 1980 to 2022, annual-mean AA reached 4.2 (Arctic defined as north of 70°N). Climate models simulate AA but fail to reproduce its magnitude. Sweeney et al. attributed much of this model–observation discrepancy to internal variability. AA shows seasonality and so does the discrepancy. Spring (March–May) shows the largest gap: Observed AA is 4.2, while the multimodel mean is 2.7. This raises several questions: 1) What role does internal variability play in observed spring AA? 2) How does simulated spring AA compare to observations when internal variability is removed? 3) If internal variability is significant, what mechanisms drive it? To address these, we adapted the machine learning algorithm from Sweeney et al., training on simulated multidecadal spring SAT and sea level pressure (SLP) trend maps. Our results show that internal variability enhanced spring Arctic warming by 37% and reduced global warming by 10%. Removing internal variability reconciles the spring AA discrepancy. The estimated internal contribution to Arctic spring warming is supported by an independent dynamical adjustment approach. We identify an atmospheric circulation pattern in observations associated with this internal warming. Observed internal Siberian SAT and SLP trends follow the simulated SAT–SLP relationship but lie at the distribution’s extreme, suggesting models generally underestimate internal variability unless the observed configuration reflects a rare real-world realization.

Arctic↗