Search NASA⌕ Search

SEARCH · Search NASA

Results for “Arctic, Earth System 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 307 records · Page 17

Arctic shrub expansion generates regional variation in litter decomposition by altering litter quality and the decomposition environment

Abstract The expansion of deciduous shrubs into the graminoid‐dominated arctic tundra is expected to alter litter decomposition by changing litter quality and local abiotic and biotic conditions. However, it remains unclear how shrub expansion will affect litter decomposition at regional scales, where macroclimate is expected to be the dominant regulator of decomposition. To determine the relative influence of macroclimate and local controls on regional patterns of litter mass loss and nitrogen release, we conducted two hierarchical litter decomposition experiments across spatial scales. We decomposed leaf and root litter from a prominent graminoid ( Eriophorum vaginatum ) and three genera of deciduous shrubs ( Betula , Alnus and Salix ) for 1 year within replicated plots at five sites spanning a 160 km latitudinal gradient in northern Alaska. Using Eriophorum litter as a substrate, we found that macroclimate was the primary regulator of mass loss but had opposing effects on leaf and root litter. As summer temperature increased along the latitudinal gradient (11.9 to 13.9°C), leaf litter mass loss increased by 20% whereas root litter mass loss decreased by 33%. Leaf nitrogen release also increased with summer temperature. Conversely, root nitrogen release was controlled by the vegetation type of the decomposition environment. Using different shrub litters as substrates, we found that litter quality and its interaction with soil microclimate and macroclimate controlled decomposition. Overall, shrub root litters decomposed faster than Eriophorum root litter, losing 53% more mass and 190% more nitrogen across all sites and decomposition environments. For leaf litter, however, patterns varied by litter genus, with Salix losing more mass and Betula and Alnus losing less mass than Eriophorum . Our findings demonstrate that shrub expansion in the Arctic can regulate leaf and root litter decomposition at the regional scale through its effects on local controls, primarily litter quality. Ongoing increases in shrub cover are likely to accelerate the turnover of root litter carbon and nitrogen pools in tundra ecosystems. Therefore, including shrub‐related processes in Earth system models will improve our ability to predict regional‐scale litter decomposition and its effects on carbon and nutrient cycling in a warming Arctic. Read the free Plain Language Summary for this article on the Journal blog.

Vozzo, Justin T. [Department of Natural Resources ↗

Quantifying the Known Unknown: Including Marine Sources of Greenhouse Gases in Climate Modeling

Researchers have recently estimated that Arctic submarine permafrost currently traps 60 billion tons of methane and contains 560 billion tons of organic carbon in seafloor sediments and soil, a giant pool of carbon with potentially large feedbacks on the climate system. Unlike terrestrial permafrost, the submarine permafrost system has remained a “known unknown” because of the difficulty in acquiring samples and measurements. Consequently, this potentially large carbon stock never yet considered in global climate models or policy discussions, represents a real wildcard in our understanding of Earth’s climate. This report summarizes our group’s effort at developing a numerical modeling framework designed to produce a first-of-its-kind estimate of Arctic methane gas releases from the marine sediments to the water column, and potentially to the atmosphere, where positive climate feedback may occur. Newly developed modeling capability supported by the Laboratory Directed Research and Development (LDRD) program at Sandia National Laboratories now gives us the ability to probabilistically map gas distribution and quantity in the seabed by using a hybrid approach of geospatial machine learning, and predictive numerical thermodynamic ensemble modeling. The novelty in this approach is its ability to produce maps of useful data in regions that are only sparsely sampled, a common challenge in the Arctic, and a major obstacle to progress in the past. By applying this model to the circum-Arctic continental shelves and integrating the flux of free gas from in situ methanogenesis and dissociating gas hydrates from the sediment column under climate forcing, we can provide the most reliable estimate of a spatially and temporally varying source term for greenhouse gas flux that can be used by global oceanographic circulation and Earth system models (such as DOE’s E3SM). The result will allow us to finally tackle the wildcard of the submarine permafrost carbon system, and better inform us about the severity of future national security threats that sustained climate change poses.

54 ENVIRONMENTAL SCIENCES↗

High Latitude Atmospheric Rivers: Teleconnections and Impacts

Atmospheric rivers (ARs) are long, narrow synoptic scale features that act as moisture transport vehicles by moving water vapor poleward. They are typically found in storm tracks and often have profound implications for regional hydroclimate. For high latitude locales such as Antarctica and the Arctic, ARs are fewer in number but often have an outsized impact by either producing significant snow accumulation events, or, depending on their thermal characteristics, initiating melt events on ice sheets, ice shelves, and sea ice. In this study, we apply both reanalysis and climate model data to evaluate high latitude ARs for both the Arctic and Antarctic in the context of teleconnections and modes of variability. Specifically, we fold in both MERRA-2 and ERA5 reanalysis as well as historical and climate change simulations from the CESM2 (Community Earth System Model, Version 2) and E3SMv2 (Energy Exascale Earth System Model, v2) with multiple ensemble members to better understand natural variability versus climate change signals. An overview of AR climatology and impacts will be presented, including a baseline uncertainty quantification analysis using data from ARTMIP (the Atmospheric River Tracking Method Intercomparison Project). AR impacts during specific modes of variability will be characterized by diagnosing precipitation and boundary layer temperature associated with Ars influencing the cryosphere. Modes of variability investigated include the SAM (Southern Annular Mode), PDO (Pacific Decadal Oscillation), PSA2 (Pacific South American Mode 2), the IOD (Indian Ocean Dipole), and the AO (Arctic Oscillation).

Atmospheric Rivers↗

Improved ELMv1-ECA simulations of zero-curtain periods and cold-season CH 4 and CO 2 emissions at Alaskan Arctic tundra sites

Field measurements have shown that cold-season methane (CH 4 ) and carbon dioxide (CO 2 ) emissions contribute a substantial portion to the annual net carbon emissions in permafrost regions. However, most earth system land models do not accurately reproduce cold-season CH 4 and CO 2 emissions, especially over the shoulder (i.e., thawing and freezing) seasons. Here we use the Energy Exascale Earth System Model (E3SM) land model version 1 (ELMv1-ECA) to tackle this challenge and fill the knowledge gap of how cold-season CH 4 and CO 2 emissions contribute to the annual totals at Alaska Arctic tundra sites. Specifically, we improved the ELMv1-ECA soil water phase-change scheme, environmental controls on microbial activity, and the methane module. Results demonstrate that both soil temperature and the duration of zero-curtain periods (i.e., the fall period when soil temperatures linger around 0°C) simulated by the updated ELMv1-ECA were greatly improved; e.g., the mean absolute error (MAE) in zero-curtain durations at 12 cm depth was reduced by 62 % on average. Furthermore, the MAEs of simulated cold-season carbon emissions at three tundra sites were improved by 72 % and 70 % on average for CH 4 and CO 2 , respectively. Overall, CH 4 emitted during the early cold season (September and October), which often includes most of the zero-curtain period in Arctic tundra, accounted for more than 50 % of the total emissions throughout the entire cold season (September to May) in the model, compared with around 49.4 % (43 %–58 %) in observations. From 1950 to 2017, both CO 2 emissions during the zero-curtain period and during the entire cold season showed increasing trends, for example, of 0.17 and 0.36 gC m -2 yr -1 at Atqasuk. This study highlights the importance of zero-curtain periods in facilitating cold-season CH 4 and CO 2 emissions from tundra ecosystems.

58 GEOSCIENCES↗

Daytime Cloud Property Retrievals Over the Arctic from Multispectral MODIS Data

Improving climate model predictions over Earth's polar regions requires a complete understanding of polar clouds properties. Passive satellite remote sensing techniques can be used to retrieve macro and microphysical properties of polar cloud systems. However, over the Arctic, there is minimal contrast between clouds and the background snow surface observed in satellite data, especially for visible wavelengths. This makes it difficult to identify clouds and retrieve their properties from space. Variable snow and ice cover, temperature inversions, and the predominance of mixed-phase clouds further complicate cloud property identification. For this study, the operational Clouds and the Earth s Radiant Energy System (CERES) cloud mask is first used to discriminate clouds from the background surface in Terra Moderate Resolution Imaging Spectroradiometer (MODIS) data. A solar-infrared infrared nearinfrared technique (SINT) first used by Platnick et al. (2001) is used here to retrieve cloud properties over snow and ice covered regions.

Spangenberg, Douglas A.↗

NGEE Arctic Field-to-Model

This repository contains workshop documentation and shell scripts/computational infrastructure for building, running, and analyzing output from the DOE Energy Exascale Earth System Model (E3SM) Land Model (ELM).

Fiorella, Rich [@lanl]↗

Significant impact of forcing uncertainty in a large ensemble of climate model simulations

Forcing due to solar and volcanic variability, on the natural side, and greenhouse gas and aerosol emissions, on the anthropogenic side, are the main inputs to climate models. Reliable climate model simulations of past and future climate change depend crucially upon them. Here we analyze large ensembles of simulations using a comprehensive Earth System Model to quantify uncertainties in global climate change attributable to differences in prescribed forcings. The different forcings considered here are those used in the two most recent phases of the Coupled Model Intercomparison Project (CMIP), namely CMIP5 and CMIP6. We show significant differences in simulated global surface air temperature due to volcanic aerosol forcing in the second half of the 19th century and in the early 21st century. The latter arise from small-to-moderate eruptions incorporated in CMIP6 simulations but not in CMIP5 simulations. We also find significant differences in global surface air temperature and Arctic sea ice area due to anthropogenic aerosol forcing in the second half of the 20th century and early 21st century. These differences are as large as those obtained in different versions of an Earth System Model employing identical forcings. In simulations from 2015 to 2100, we find significant differences in the rates of projected global warming arising from CMIP5 and CMIP6 concentration pathways that differ slightly but are equivalent in terms of their nominal radiative forcing levels in 2100. Our results highlight the influence of assumptions about natural and anthropogenic aerosol loadings on carbon budgets, the likelihood of meeting Paris targets, and the equivalence of future forcing scenarios.

54 ENVIRONMENTAL SCIENCES↗

FUTURE LEVELS OF AEROSOLS & RADIATIVE FORCING OVER THE ARCTIC

The Arctic Monitoring and Assessment Programme (AMAP) is currently assessing the impacts of Short-Lived Climate Forcers (SLCF) on Arctic climate and air quality. In support of the assessment, we used the NASA Goddard Institute of Space Sciences (GISS) Earth System Model (modelE2.1), with prescribed sea surface temperature and sea-ice fraction, to simulate SLCF concentrations globally between 1995 and 2015. Two simulations were conducted, using the One-Moment Aerosol (OMA) and the Multiconfiguration Aerosol TRacker of mIXing state (MATRIX) aerosol modules. OMA is a mass-based scheme in which aerosols are assumed to remain externally mixed and have a prescribed and constant size distribution, while MATRIX is an aerosol microphysics scheme based on the quadrature method of moments, which is able to explicitly simulate the mixing state of aerosols. Anthropogenic emissions from the ECLIPSE v6b emissions database were used, along with emissions from aircrafts and open biomass burning from the Coupled Model Intercomparison Project Phase 6 (CMIP6), while the natural emissions of sea salt, DMS, isoprene and dust are calculated interactively. The simulated monthly surface concentrations of sulfate (SO4), black carbon (BC), organic carbon (OA), and ozone (O3) are compared with observations from a set of Arctic stations, extracted from the EBAS and IMPROVE databases, as well as a few additional locations. Simulated aerosol optical depths (AOD) are also compared with Advanced Very-High Resolution Radiometer (AVHRR). The study will present the evaluation of the modelE2.1 in simulating SLCF levels over the Arctic using different aerosol schemes, along with observed and simulated trends of SLCFs over the Arctic between 1995 and 2015.

Aersosols↗

Time irreversibility as an indicator of approaching tipping points in Earth subsystems

With shifting environmental trends, many Earth system elements may be poised to undergo critical transitions or ‘tipping’. Reliable anticipation of these tipping elements is vital to inform policy decisions. Many of the current methods for tipping point detection are based on loss of resilience or ‘critical slowdown’ of the system as it approaches a tipping point. However, these methods are prone to false alarms; the detected slowdown may be an artifact of nonstationary noise unrelated to tipping behavior. Here, we explore the efficacy of early warning signs based on a nonequilibrium thermodynamics framework. The model-free detection method relies on the increased intrinsic time-irreversibility due to detailed balance breaking, preceding the onset of tipping or instabilities. We demonstrate that these EWSs are effective for tipping point detection and robust against false alarms due to nonstationary noise, using idealized models for two key elements of the Earth system that are prone to tipping: the Atlantic Meridional Overturning Circulation and Arctic sea-ice loss.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Flow Routing on the Unstructured Voronoi Meshes in Earth System Modeling

Flow routing is a fundamental process of Earth System Models' (ESMs) river component. Traditional flow routing models rely on Cartesian rectangular meshes, which exhibit limitations, particularly when coupled with unstructured mesh-based ocean components. They also lack the support for regionally refined models. While previous studies have highlighted the potential benefits of unstructured meshes for flow routing, their widespread application and comprehensive evaluation within ESMs remain limited. This study extends the river component of the Energy Exascale Earth System Model to unstructured Voronoi meshes. We evaluated the model's performance in simulating river discharge and water depth across three watersheds spanning the Arctic, temperate, and tropical regions. The results show that while providing several benefits, unstructured mesh-based flow routing can achieve comparable performance to structured mesh-based routing, and their difference is often less than 10%. Although the unstructured mesh-based method could address several existing limitations, this research also shows that additional improvements in the numerical method are needed to fully exploit the advantages of unstructured mesh for hydrologic and ESMs.

54 ENVIRONMENTAL SCIENCES↗

Extrapolar Cloud Feedbacks as a Driver of Arctic Amplification

The role of cloud feedbacks in Arctic amplification (AA) of anthropogenic warming remains unclear. Traditional feedback analysis diagnoses the net cloud feedback as strongly positive in the tropics but either weak or negative in the Arctic, suggesting that AA would be amplified if cloud feedbacks were suppressed. However, in cloud-locking experiments using the slab ocean version of the Energy Exascale Earth System Model (E3SM), we find that suppressing cloud feedbacks results in a substantial decrease in AA under greenhouse gas forcing. We show that the increase in AA from cloud feedbacks arises from two main mechanisms: 1) the additional energy contributed by positive cloud feedbacks in the tropics leads to increased poleward moist atmospheric heat transport (AHT) which then amplifies Arctic warming; and 2) the additional Arctic warming is amplified by positive noncloud feedbacks in the region, together making extrapolar cloud feedbacks amplify AA. We also find that cloud changes can modify the strength of noncloud feedback, but that modification has a small effect on Arctic warming. We further examine the role of cloud feedbacks in AA using a moist energy balance model, which demonstrates that interactions of cloud feedbacks with moist AHT and other positive feedbacks dominate the influence of clouds on the pattern of surface warming. However, the contribution of cloud-induced changes in noncloud feedbacks on AA is relatively minor. Furthermore, these results demonstrate that traditional attributions of AA, that are based on local feedback analysis, overlook key interactions between extrapolar cloud changes, poleward AHT, and noncloud feedbacks in the Arctic.

58 GEOSCIENCES↗

Summertime low clouds mediate the impact of the large-scale circulation on Arctic sea ice

The rapid Arctic sea ice retreat in the early 21 st century is believed to be driven by several dynamic and thermodynamic feedbacks, such as ice-albedo feedback and water vapor feedback. However, the role of clouds in these feedbacks remains unclear since the causality between clouds and these processes is complex. Here, we use NASA CERES satellite products and NCAR CESM model simulations to suggest that summertime low clouds have played an important role in driving sea ice melt by amplifying the adiabatic warming induced by a stronger anticyclonic circulation aloft. The upper-level high pressure regulates low clouds through stronger downward motion and increasing lower troposphere relative humidity. The increased low clouds favor more sea ice melt via emitting stronger longwave radiation. Then decreased surface albedo triggers a positive ice-albedo feedback, which further enhances sea ice melt. Considering the importance of summertime low clouds, accurate simulation of this process is a prerequisite for climate models to produce reliable future projections of Arctic sea ice.

54 ENVIRONMENTAL SCIENCES↗

Water balance response of permafrost-affected watersheds to changes in air temperatures

Abstract Observations show increases in river discharge to the Arctic Ocean especially in winter over the last decades but the physical mechanisms driving these changes are not yet fully understood. We hypothesize that even in the absence of a precipitation increase, permafrost degradation alone can lead to increased annual river runoff. To test this hypothesis we perform 12 millennium-long simulations over an idealized hypothetical watershed (1 km 2 ) using a distributed, physically based water balance model (Water flow and Balance Simulation Model, WaSiM). The model is forced by both a hypothetical warming defined by an air temperature increase of 7.5 ∘ C over 100 years, and a corresponding cooling scenario. To assess model sensitivity we vary soil saturated hydraulic conductivity and lateral subsurface flow configuration. Under the warming scenario, changes in subsurface water transport due to ground temperature changes result in a 7%–14% increase in annual runoff accompanied by a 6%–20% decrease in evapotranspiration. The increase in runoff is most pronounced in winter. Hence, the simulations demonstrate that changes in permafrost characteristics due to climate warming and associated changes in evapotranspiration provide a plausible mechanism for the observed runoff increases in Arctic watersheds. In addition, our experiments show that when lateral subsurface moisture transport is not included, as commonly done in global-scale Earth System Models, the equilibrium water balance in response to the warming or cooling is similar to the water balance in simulations where lateral subsurface transport is included. However, the transient changes in water balance components prior to reaching equilibrium differ greatly between the two. For example, for high saturated hydraulic conductivity only when lateral subsurface transport is considered, a period of decreased runoff occurs immediately after the warming. This period is characterized by a positive change in soil moisture storage caused by the soil moisture deficit developed during prior cooling.

Debolskiy, Matvey V. (ORCID:0000000296343627)↗

Near‐Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub‐Arctic Watershed

Abstract Discontinuous permafrost environments exhibit strong spatial heterogeneity at scales too small to be driven by weather forcing or captured by Earth System Models. Here we analyze effects of observed spatial heterogeneity in soil and vegetation properties, hydrology, and thermal dynamics on ecosystem carbon dynamics in a watershed on the Seward Peninsula in Alaska. We apply a Morris global sensitivity analysis to a process‐rich, successfully tested terrestrial ecosystem model (TEM), ecosys , varying soil properties, boundary conditions, and weather forcing. We show that landscape heterogeneity strongly impacts soil temperatures and vegetation composition. Snow depth, O‐horizon thickness, and near‐surface water content, which vary at scales of O(m), control the soil thermal regime more than an air temperature gradient corresponding to a 140 km north–south distance. High shrub productivity is simulated only in talik (perennially unfrozen) soils with high nitrogen availability. Through these effects on plant and permafrost dynamics, landscape heterogeneity impacts ecosystem productivity. Simulations with near‐surface taliks have higher microbial respiration (by 78.0 gC m −2 yr −1 ) and higher net primary productivity (by 104.9 gC m −2 yr −1 ) compared to runs with near‐surface permafrost, and simulations with high shrub productivity have outlying values of net carbon uptake. We explored the prediction uncertainty associated with ignoring observed landscape heterogeneity, and found that watershed net carbon uptake is 60% larger when heterogeneity is accounted for. Our results highlight the complexity inherent in discontinuous permafrost environments and demonstrate that missing representation of subgrid heterogeneity in TEMs could bias predictions of high‐latitude carbon budget.

54 ENVIRONMENTAL SCIENCES↗

Seasonal Dependent Impact of Ice Cloud Longwave Scattering on the Polar Climate

Most climate models neglect cloud longwave (LW) scattering because scattering is considered negligible compared to strong LW absorption by clouds and greenhouse gases. While this rationale is valid for simulating extrapolar regions, it is questionable for the polar regions, where the atmosphere is dry and hence has weak absorption, and ice clouds that have strong scattering capability frequently occur. Using the slab-ocean Community Earth System Model, we show that ice cloud LW scattering can warm winter surface air temperature by 0.8–1.8 K in the Arctic and 1.3–1.9 K in the Antarctic, while this warming becomes much weaker in polar summer. Such scattering effect cannot be correctly assessed when sea surface temperature and sea ice are prescribed as this effect is manifested through a surface-atmosphere coupling. Furthermore, cloud LW scattering is a necessity for the correct simulation of polar climate and surface radiation budget, especially in the winter.

54 ENVIRONMENTAL SCIENCES↗

Expanding the design space of stratospheric aerosol geoengineering to include precipitation-based objectives and explore trade-offs

Previous climate modeling studies demonstrate the ability of feedback-regulated, stratospheric aerosol geoengineering with injection at multiple independent latitudes to meet multiple simultaneous temperature-based objectives in the presence of anthropogenic climate change. However, the impacts of climate change are not limited to rising temperatures but also include changes in precipitation, loss of sea ice, and many more; knowing how a given geoengineering strategy will affect each of these climate metrics is vital to understanding the limits and trade-offs of geoengineering. Here, we first introduce a new method of visualizing the design space in which desired climate outcomes are represented by 2-D surfaces on a 3-D graph. Surface orientations represent how different injection choices influence that objective, and intersecting surfaces represent objectives which can be met simultaneously. Using this representation as a guide, we present simulations of two new strategies for feedback-regulated aerosol injection, using the Community Earth System Model with the Whole Atmosphere Community Climate Model – CESM1(WACCM). The first simultaneously manages global mean temperature, tropical precipitation centroid, and Arctic sea ice extent, while the second manages global mean precipitation, tropical precipitation centroid, and Arctic sea ice extent. Both simulations control the tropical precipitation centroid to within 5 % of the goal, and the latter controls global mean precipitation to within 1% of the goal. Additionally, the first simulation overcompensates sea ice, while the second undercompensates sea ice; all of these results are consistent with the expectations of our design space model. In addition to showing that precipitation-based climate metrics can be managed using feedback alongside other goals, our simulations validate the utility of our design space visualization in predicting our climate model behavior under a given geoengineering strategy, and together they help illustrate the fundamental limits and trade-offs of stratospheric aerosol geoengineering.

58 GEOSCIENCES↗

Snow Depth from AMSR-2 Using Multispectral Satellite Data in an Artificial Neural Network

By using diffusion theory and Monte Carlo lidar radiative transfer simulations, Hu et al. (2022b) has derived snow depth from the first-, second- and third-order moments of the lidar backscattering pathlength distribution. Lu et al. (2022) calculated the snow depth by applying the methods to the satellite ICESat-2 lidar measurements over the Arctic sea ice, as well as land surfaces of Northern Hemisphere. In this paper, an artificial neural network (ANN) algorithm, employing several channels from Advanced Microwave Scanning Radiometer 2 (AMSR-2) and the humidity vertical profiles from Global Modeling and Assimilation Office (GMAO) Goddard Earth Observing System for Instrument Teams (GEOS-IT) product, is trained to determine snow depth identified by time and geolocation matched 2019 ICESat-2 snow-depth data during winter months over the Arctic sea ice. The trained ANN snow-depth was applied to 2018 AMSR-2 clear pixel data, although the algorithms perform reasonably well in thinner clouds. The validation data (different from the training set) of ANN snow depth from AMSR-2 showed a good agreement with time matched and co-located snow-depth values from ICESat-2. The bias was near zero, with mean absolute error (MAE) 0.05 cm and a root-mean-square-error (RMSE) 0.08 cm. Prior applying the trained ANN snow depth to AMSR-2 data, a cloud screening algorithm was developed with a similar approach. A separate ANN cloud mask was trained to determine an AMSR-2 pixel is clear or cloudy with time and geolocation matched 2015 CALIOP Vertical Feature Mask (VFM) over Arctic sea ice. The ANN cloud mask from AMSR-2 under-estimated cloud fraction by 3-6% compared to CALIOP . The additional research is needed to conclusively evaluate the ANN cloud mask accuracy. Finally, this paper will lay the foundation for a sustained long-term snowfall and snow-storm monitoring system. The future Cloud Aerosol LIdar for Global scale Observations of the ocean-Land Atmosphere system (CALIGOLA) mission will provide a means to calculate snow depth from the lidar backscattering pathlength distribution, benefiting from the UV, visible and infrared pulses. With the calculated snow depth as the truth one could develop a machine learning algorithm, as it was done in this paper, using a passive microwave instrument available at that time to generate a wide range of snow depth data, covering extensive spatial areas in the cross-orbit direction.

Neural Network↗

The Salience of Nonlinearities in the Boreal Winter Response to ENSO: Arctic Stratosphere and Europe

The salience of nonlinearities in the Arctic stratospheric response to El Niño - Southern Oscillation (ENSO) is evaluated in a 41 member ensemble of the period 1980 to 2009 in the NASA Goddard Earth Observing System Chemistry-Climate Model. Central Pacific El Niño (CP EN) events preferentially increase zonal wavenumber 2, while Eastern Pacific El Niño (EP EN) events preferentially increase zonal wavenumber 1, though these nonlinearities become statistically robust only after more than 45 events are considered. This leads to a difference in the morphology of SSW events preferentially forced by each ENSO flavor. While the zonal-mean Arctic stratospheric response is little different between EP and CP EN once one accounts for the relative weakness of CP events, the tropospheric subpolar negative North Atlantic Oscillation response is stronger for CP events. Other nonlinearities are even less salient: there is no indication of any nonlinearities in the wave driving between moderate EP EN and extreme EP EN, nor between EN and La Niña (LN). However in approximately 15% of the integrations within our ensemble, there is no difference in SSW frequency between EN and LN, suggesting that a similar SSW frequency for both EN and LN, as has occurred over the past 60 years, can occur by chance. The intra-ensemble variability in LN SSW frequency can be related back to intra-ensemble variability in the tropospheric response to LN.

Arctic↗