Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atmospheric Model Parameterizations”

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 199 records · Page 11

Numerical Modeling of the Global Atmosphere

Under this grant, we continued development and evaluation of the updraft downdraft model for cumulus parameterization. The model includes the mass, rainwater and vertical momentum budget equations for both updrafts and downdrafts. The rainwater generated in an updraft falls partly inside and partly outside the updraft. Two types of stationary solutions are identified for the coupled rainwater budget and vertical momentum equations: (1) solutions for small tilting angles, which are unstable; (2) solutions for large tilting angles, which are stable. In practical applications, we select the smallest stable tilting angle as an optimum value. The model has been incorporated into the Arakawa-Schubert (A-S) cumulus parameterization. The results of semi-prognostic and single-column prognostic tests of the revised A-S parameterization show drastic improvement in predicting the humidity field. Cheng and Arakawa presents the rationale and basic design of the updraft-downdraft model, together with these test results. Cheng and Arakawa, on the other hand gives technical details of the model as implemented in current version of the UCLA GCM.

Arakawa, Akio↗

Interseasonal Variations in the Middle Atmosphere Forced by Gravity Waves

In our Numerical Spectral Model (NSM), which incorporates Hines' Doppler Spread Parameterization, gravity waves (GW) propagating in the east/west direction can generate the essential features of the observed equatorial oscillations in the zonal circulation and in particular the QBO (quasi-biennial oscillation) extending from the stratosphere into the upper mesosphere. We report here that the NSM also produces inter-seasonal variations in the zonally symmetric (m = 0) meridional circulation. A distinct but variable meridional wind oscillation (MWO) is generated, which appears to be the counterpart to the QBO. With a vertical grid-point resolution of about 0.5 km, the NSM produces the MWO through momentum deposition of GWs propagating in the north/south direction. The resulting momentum source represents a third (generally odd) order non-linear function of the meridional winds, and this enables the oscillation, as in the case of the QBO for the zonal winds. Since the meridional winds are relatively small compared to the zonal winds, however, the vertical wavelength that maintains the MWO is much smaller, i.e., only about 10 km instead of 40 km for the QBO. Consistent with the associated increase of the viscous stress, the period of the MWO is then short compared with that of the QBO, i.e., only about two to four months. Depending on the strength of the GW forcing, the computed amplitudes of the MWO are typically 4 m/s in the upper stratosphere and mesosphere, and the associated temperature amplitudes are between about 2 and 3 K. These amplitudes may be observable with the instruments on the TIMED spacecraft. Extended computer simulations with the NSM in 2D (two-dimensional) and 3D (three-dimensional) reveal that the MWO is modulated by and in turn influences the QBO.

Mayr, H. G.↗

Interseasonal Variations in the Middle Atmosphere Forced by Gravity Waves

In our Numerical Spectral Model (NSM), which incorporates Hines' Doppler Spread Parameterization, gravity waves (GW) propagating in the east/west direction can generate the essential features of the observed equatorial oscillations of the zonal circulation and in particular the QBO (quasi-biennial oscillation) extending from the stratosphere into the upper mesosphere. We report here that the NSM also produces inter-seasonal variations in the zonally symmetric meridional circulation. A distinct meridional oscillation (MO) is generated, which appears to be the counterpart to the QBO. With a vertical grid-point resolution of about 0.5 km, the NSM produces the MO through momentum deposition of GW's propagating in the north/south direction. This process is inherently non-linear, of third (odd) order, which enables the oscillation. Since the meridional winds are relatively small compared to the zonal winds, the vertical wavelength required to maintain the MO is also smaller, i.e., only about 10 km instead of the 30 km for the QBO. The corresponding viscous stress is then larger, and the period of the MO is thus short compared with that of the QBO, i.e., only about 3 to 4 months. Depending on the strength of the GW forcing, the computed amplitudes of the meridional wind oscillation are typically 5 m/s in the upper stratosphere and mesosphere, and the associated temperature amplitudes are between about 2 and 3 K. These amplitudes may be observable with the instruments on the TIMED spacecraft. Extended computer simulations with the NSM in 2D and 3D reveal that the MO at low latitudes is modulated by the QBO and in turn can influence it to produce a hemispherically asymmetric component. The annual circulation from the summer to the winter hemisphere is likely to play an important role.

Mayr, H. G.↗

Improvement of the GEOS-5 AGCM upon Updating the Air-Sea Roughness Parameterization

The impact of an air-sea roughness parameterization over the ocean that more closely matches recent observations of air-sea exchange is examined in the NASA Goddard Earth Observing System, version 5 (GEOS-5) atmospheric general circulation model. Surface wind biases in the GEOS-5 AGCM are decreased by up to 1.2m/s. The new parameterization also has implications aloft as improvements extend into the stratosphere. Many other GCMs (both for operational weather forecasting and climate) use a similar class of parameterization for their air-sea roughness scheme. We therefore expect that results from GEOS-5 are relevant to other models as well.

Garfinkel, C. I.↗

Maximum entropy modeling to identify physical drivers of shallow snowpack heterogeneity using unpiloted aerial system (UAS) lidar

Understanding the spatial variability of the snowpack is valuable for hydrologists and ecologists seeking to predict hydrological processes in a cold region. Snow distribution is a function of interactions among static variables, such as terrain, vegetation, and soil properties, and dynamic meteorological variables, such as solar radiation, wind speed and direction, and soil moisture. However, identifying the dominant physical drivers responsible for spatial patterns of the snowpack, particularly for ephemeral, shallow snowpacks, has been challenging due to the lack of the high-resolution snowpack and physical variables with high vertical accuracy as well as inherent limitations in traditional approaches. This study uses an Unpiloted Aerial System (UAS) lidar-based snow depth and static variables (1-m spatial resolution) to analyze field-scale spatial structures of snow depth and apply the Maximum Entropy (MaxEnt) model to identify primary controls over open terrain and forests at the University of New Hampshire Thompson Farm Research Observatory, New Hampshire, United States. We found that, among nine topographic and soil variables, plant functional type and terrain roughness contribute up to 80% and 76% of relative importance in the MaxEnt framework to predict locations of deeper or shallower snowpacks, respectively, across a mixed temperate forested and field landscape. Soil variables, such as organic matter and saturated hydraulic conductivity, were also important controls (up to 70% and 81%) on snow depth spatial variations for both open and forested landscapes suggesting spatial variations in soil variables under snow can control thermal transfer among soil, snowpack, and surface-atmosphere. This work contributes to improving land surface and snow models by informing parameterization of the sub-grid scale snow depths, down-scaling remotely sensed snow products, and understanding field scale snow states.

Eunsang Cho↗

Surface Hydrology in Global River Basins in the Off-Line Land-Surface GEOS Assimilation (OLGA) System

Land surface hydrology for the Off-line Land-surface GEOS Analysis (OLGA) system and Goddard Earth Observing System (GEOS-1) Data Assimilation System (DAS) has been examined using a river routing model. The GEOS-1 DAS land-surface parameterization is very simple, using an energy balance prediction of surface temperature and prescribed soil water. OLGA uses near-surface atmospheric data from the GEOS-1 DAS to drive a more comprehensive parameterization of the land-surface physics. The two global systems are evaluated using a global river routing model. The river routing model uses climatologic surface runoff from each system to simulate the river discharge from global river basins, which can be compared to climatologic river discharge. Due to the soil hydrology, the OLGA system shows a general improvement in the simulation of river discharge compared to the GEOS-1 DAS. Snowmelt processes included in OLGA also have a positive effect on the annual cycle of river discharge and source runoff. Preliminary tests of a coupled land-atmosphere model indicate improvements to the hydrologic cycle compared to the uncoupled system. The river routing model has provided a useful tool in the evaluation of the GCM hydrologic cycle, and has helped quantify the influence of the more advanced land surface model.

Bosilovich, Michael G.↗

On the Development of Parameterized Linear Analytical Longitudinal Airship Models

In order to explore Titan, a moon of Saturn, airships must be able to traverse the atmosphere autonomously. To achieve this, an accurate model and accurate control of the vehicle must be developed so that it is understood how the airship will react to specific sets of control inputs. This paper explains how longitudinal aircraft stability derivatives can be used with airship parameters to create a linear model of the airship solely by combining geometric and aerodynamic airship data. This method does not require system identification of the vehicle. All of the required data can be derived from computational fluid dynamics and wind tunnel testing. This alternate method of developing dynamic airship models will reduce time and cost. Results are compared to other stable airship dynamic models to validate the methods. Future work will address a lateral airship model using the same methods.

aerobot↗

Sensitivity of June Near-Surface Temperatures and Precipitation in the Eastern United States to Historical Land Cover Changes Since European Settlement

Land cover changes alter the near surface weather and climate. Changes in land surface properties such as albedo, roughness length, stomatal resistance, and leaf area index alter the surface energy balance, leading to differences in near surface temperatures. This study utilized a newly developed land cover data set for the eastern United States to examine the influence of historical land cover change on June temperatures and precipitation. The new data set contains representations of the land cover and associated biophysical parameters for 1650, 1850, 1920, and 1992, capturing the clearing of the forest and the expansion of agriculture over the eastern United States from 1650 to the early twentieth century and the subsequent forest regrowth. The data set also includes the inferred distribution of potentially water-saturated soils at each time slice for use in the sensitivity tests. The Regional Atmospheric Modeling System, equipped with the Land Ecosystem-Atmosphere Feedback (LEAF-2) land surface parameterization, was used to simulate the weather of June 1996 using the 1992, 1920, 1850, and 1650 land cover representations. The results suggest that changes in surface roughness and stomatal resistance have caused present-day maximum and minimum temperatures in the eastern United States to warm by about 0.3 C and 0.4 C, respectively, when compared to values in 1650. In contrast, the maximum temperatures have remained about the same, while the minimums have cooled by about 0.1 C when compared to 1920. Little change in precipitation was found.

Strack, John E.↗

Parameterizing Gravity Waves and Understanding Their Impacts on Venus' Upper Atmosphere

The complexity of Venus’ upper atmospheric circulation is still being investigated. Simulations of Venus’ upper atmosphere largely depend on the utility of Rayleigh Friction (RF) as a driver and necessary process to reproduce observations (i.e. temperature, density, nightglow emission). Currently, there are additional observations which provide more constraints to help characterize the driver(s) of the circulation. This work will largely focus on the impact parameterized gravity waves have on Venus’ upper atmosphere circulation within a three dimensional hydrodynamic model (Venus Thermospheric General Circulation Model).

Brecht, A. S.↗

Understanding Cirrus Ice Crystal Number Variability for Different Heterogeneous Ice Nucleation Spectra

Along with minimizing parameter uncertainty, understanding the cause of temporal and spatial variability of the nucleated ice crystal number, Ni, is key to improving the representation of cirrus clouds in climate models. To this end, sensitivities of Ni to input variables like aerosol number and diameter provide valuable information about nucleation regime and efficiency for a given model formulation. Here we use the adjoint model of the adjoint of a cirrus formation parameterization (Barahona and Nenes, 2009b) to understand Ni variability for various ice-nucleating particle (INP) spectra. Inputs are generated with the Community Atmosphere Model version 5, and simulations are done with a theoretically derived spectrum, an empirical lab-based spectrum and two field-based empirical spectra that differ in the nucleation threshold for black carbon particles and in the active site density for dust. The magnitude and sign of Ni sensitivity to insoluble aerosol number can be directly linked to nucleation regime and efficiency of various INP. The lab-based spectrum calculates much higher INP efficiencies than field-based ones, which reveals a disparity in aerosol surface properties. Ni sensitivity to temperature tends to be low, due to the compensating effects of temperature on INP spectrum parameters; this low temperature sensitivity regime has been experimentally reported before but never deconstructed as done here.

Cirrus↗

Sensitivity of Age-of-Air Calculations to the Choice of Advection Scheme

The age of air has recently emerged as a diagnostic of atmospheric transport unaffected by chemical parameterizations, and the features in the age distributions computed in models have been interpreted in terms of the models' large-scale circulation field. This study shows, however, that in addition to the simulated large-scale circulation, three-dimensional age calculations can also be affected by the choice of advection scheme employed in solving the tracer continuity equation, Specifically, using the 3.0deg latitude X 3.6deg longitude and 40 vertical level version of the Geophysical Fluid Dynamics Laboratory SKYHI GCM and six online transport schemes ranging from Eulerian through semi-Lagrangian to fully Lagrangian, it will be demonstrated that the oldest ages are obtained using the nondiffusive centered-difference schemes while the youngest ages are computed with a semi-Lagrangian transport (SLT) scheme. The centered- difference schemes are capable of producing ages older than 10 years in the mesosphere, thus eliminating the "young bias" found in previous age-of-air calculations. At this stage, only limited intuitive explanations can be advanced for this sensitivity of age-of-air calculations to the choice of advection scheme, In particular, age distributions computed online with the National Center for Atmospheric Research Community Climate Model (MACCM3) using different varieties of the SLT scheme are substantially older than the SKYHI SLT distribution. The different varieties, including a noninterpolating-in-the-vertical version (which is essentially centered-difference in the vertical), also produce a narrower range of age distributions than the suite of advection schemes employed in the SKYHI model. While additional MACCM3 experiments with a wider range of schemes would be necessary to provide more definitive insights, the older and less variable MACCM3 age distributions can plausibly be interpreted as being due to the semi-implicit semi-Lagrangian dynamics employed in the MACCM3. This type of dynamical core (employed with a 60-min time step) is likely to reduce SLT's interpolation errors that are compounded by the short-term variability characteristic of the explicit centered-difference dynamics employed in the SKYHI model (time step of 3 min). In the extreme case of a very slowly varying circulation, the choice of advection scheme has no effect on two-dimensional (latitude-height) age-of-air calculations, owing to the smooth nature of the transport circulation in 2D models. These results suggest that nondiffusive schemes may be the preferred choice for multiyear simulations of tracers not overly sensitive to the requirement of monotonicity (this category includes many greenhouse gases). At the same time, age-of-air calculations offer a simple quantitative diagnostic of a scheme's long-term diffusive properties and may help in the evaluation of dynamical cores in multiyear integrations. On the other hand, the sensitivity of the computed ages to the model numerics calls for caution in using age of air as a diagnostic of a GCM's large-scale circulation field.

Eluszkiewicz, Janusz↗

A Solar Radiation Parameterization for Atmospheric Studies

The solar radiation parameterization (CLIRAD-SW) developed at the Goddard Climate and Radiation Branch for application to atmospheric models are described. It includes the absorption by water vapor, O3, O2, CO2, clouds, and aerosols and the scattering by clouds, aerosols, and gases. Depending upon the nature of absorption, different approaches are applied to different absorbers. In the ultraviolet and visible regions, the spectrum is divided into 8 bands, and single O3 absorption coefficient and Rayleigh scattering coefficient are used for each band. In the infrared, the spectrum is divided into 3 bands, and the k-distribution method is applied for water vapor absorption. The flux reduction due to O2 is derived from a simple function, while the flux reduction due to CO2 is derived from precomputed tables. Cloud single-scattering properties are parameterized, separately for liquid drops and ice, as functions of water amount and effective particle size. A maximum-random approximation is adopted for the overlapping of clouds at different heights. Fluxes are computed using the Delta-Eddington approximation.

Chou, Ming-Dah↗

Integration for Airborne Dust Prediction Systems and Vegetation Phenology to Track Pollen for Asthma Alerts in Public Health Decision Support Systems

The objective of the program is to assess the feasibility of combining a dust transport model with MODIS derived phenology to study pollen transport for integration with a public health decision support system. The use of pollen information has specifically be identified as a critical need by the New Mexico State Health department for inclusion in the Environmental Public Health Tracking (EPHT) program. Material and methods: Pollen can be transported great distances. Local observations of plan phenology may be consistent with the timing and source of pollen collected by pollen sampling instruments. The Dust REgional Atmospheric Model (DREAM) is an integrated modeling system designed to accurately describe the dust cycle in the atmosphere. The dust modules of the entire system incorporate the state of the art parameterization of all the major phases of the atmospheric dust life such as production, diffusion, advection, and removal. These modules also include effects of the particles size distribution on aerosol dispersion. The model was modified to use pollen sources instead of dust. Pollen release was estimated based on satellite-derived phenology of key plan species and vegetation communities. The MODIS surface reflectance product (MOD09) provided information on the start of the plant growing season, growth stage, and pollen release. The resulting deterministic model is useful for predicting and simulating pollen emission and downwind concentration to study details of phenology and meteorology and their dependencies. The proposed linkage in this project provided critical information on the location timing and modeled transport of pollen directly to the EPHT> This information is useful to support the centers for disease control and prevention (CDC)'s National EPHT and the state of New Mexico environmental public health decision support for asthma and allergies alerts.

Luvall, J. C.↗

Top-of-Atmosphere Albedo Estimation from Angular Distribution Models using Scene Identification from Satellite Cloud Property Retrievals

The next generation of Earth radiation budget satellite instruments will routinely merge estimates of global top-of-atmosphere radiative fluxes with cloud properties. This information will offer many new opportunities for validating radiative transfer models and cloud parameterizations in climate models. In this study, five months of POLarization and Directionality of the Earth's Reflectances (POLDER) 670 nm radiance measurements are considered in order to examine how satellite cloud property retrievals can be used to define empirical Angular Distribution Models (ADMs) for estimating top-of-atmosphere (TOA) albedo. ADMs are defined for 19 scene types defined by satellite retrievals of cloud fraction and cloud optical depth. Two approaches are used to define the ADM scene types: The first assumes there are no biases in the retrieved cloud properties and defines ADMs for fixed discrete intervals of cloud fraction and cloud optical depth (fixed-tau approach). The second approach involves the same cloud fraction intervals, but uses percentile intervals of cloud optical depth instead (percentile-tau approach). Albedos generated using these methods are compared with albedos inferred directly from the mean observed reflectance field. Albedos based on ADMs that assume cloud properties are unbiased (fixed-tau approach) show a strong systematic dependence on viewing geometry. This dependence becomes more pronounced with increasing solar zenith angle, reaching approximately equals 12% (relative) between near-nadir and oblique viewing zenith angles for solar zenith angles between 60 deg and 70 deg. The cause for this bias is shown to be due to biases in the cloud optical depth retrievals. In contrast, albedos based on ADMs built using percentile intervals of cloud optical depth (percentile-tau approach) show very little viewing zenith angle dependence and are in good agreement with albedos obtained by direct integration of the mean observed reflectance field (less than 1% relative error). When the ADMs are applied separately to populations consisting of only liquid water and ice clouds, significant biases in albedo with viewing geometry are observed (particularly at low sun elevations), highlighting the need to account for cloud phase both in cloud optical depth retrievals and in defining ADM scene types. ADM-derived monthly mean albedos determined for all 5 deg x 5 deg latitude/longitude regions over ocean are in good agreement (regional RMS relative errors less than 2%) with those obtained by direct integration when ADM albedos inferred from specific angular bins are averaged together. Albedos inferred from near-nadir and oblique viewing zenith angles are the least accurate, with regional RMS errors reaching approximately 5-10% (relative). Compared to an earlier study involving ERBE ADMs, regional mean albedos based on the 19 scene types considered here show a factor of 4 reduction in bias error and a factor of 3 reduction in RMS error.

Loeb, N. G.↗

The TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI). Part II: Moist Cases - The Two Waterworlds

To identify promising exoplanets for atmospheric characterization and to make the best use of observational data, a thorough understanding of their atmospheres is needed. 3D general circulation models (GCMs) are one of the most comprehensive tools available for this task and will be used to interpret observations of temperate rocky exoplanets. Due to parameterization choices made in GCMs, they can produce different results, even for the same planet. Employing four widely-used exoplanetary GCMs - ExoCAM, LMD-G, ROCKE-3D and the UM - we continue the TRAPPIST-1 Habitable Atmosphere Intercomparison by modeling aquaplanet climates of TRAPPIST-le with a moist atmosphere dominated by either nitrogen or carbon dioxide. Although the GCMs disagree on the details of the simulated regimes, they all predict a temperate climate with neither of the two cases pushed out of the habitable state. Nevertheless, the inter-model spread in the global mean surface temperature is non-negligible: 14 K and 24 K in the nitrogen and carbon dioxide dominated case, respectively. We find substantial inter-model differences in moist variables, with the smallest amount of clouds in LMDGeneric and the largest in ROCKE-3D. ExoCAM predicts the warmest climate for both cases and thus has the highest water vapor content and the largest amount and variability of cloud condensate. The UM tends to produce colder conditions, especially in the nitrogen-dominated case due to a strong negative cloud radiative effect on the day side of TRAPPIST-le. Our study highlights various biases of GCMs and emphasizes the importance of not relying solely on one model to understand exoplanet climates.

TRAPPIST↗

Lead-induced atmospheric circulations

Leads provide a significant source of heat and moisture to the Arctic winter atmosphere, and plumes from wide leads have been observed to penetrate the Arctic inversion. We have developed a two-dimensional, high-resolution, deep anelastic numerical model to investigate the atmospheric convection from leads with widths ranging from 100 m to 10 km. A second-order turbulence closure scheme is used to parameterize the atmospheric turbulence in the horizontally inhomogeneous system. This study describes how the lead-induced circulations can enhance the vertical transport of heat into the atmospheric boundary layer. This model is compared with large-eddy simulation results and with lidar observations of a lead-induced ice crystal plume. The model is used to study the effect of varying lead widths and ambient atmospheric conditions on the resultant convection from leads, and some preliminary results are described.

Alam, Afshan↗

Evaluation of the Surface Representation of the Greenland Ice Sheet in a General Circulation Model

Simulated surface conditions of the Goddard Earth Observing System model, version 5 (GEOS 5) atmospheric general circulation model (AGCM) are examined for the contemporary Greenland Ice Sheet (GrIS). A surface parameterization that explicitly models surface processes including snow compaction, meltwater percolation and refreezing, and surface albedo is found to remedy an erroneous deficit in the annual net surface energy flux and provide an adequate representation of surface mass balance (SMB) in an evaluation using simulations at two spatial resolutions. The simulated 1980-2008 GrIS SMB average is 24.7+/-4.5 cm yr(- 1) water-equivalent (w.e.) at.5 degree model grid spacing, and 18.2+/-3.3 cm yr(- 1) w.e. for 2 degree grid spacing. The spatial variability and seasonal cycle of the simulation compare favorably to recent studies using regional climate models, while results from 2 degree integrations reproduce the primary features of the SMB field. In comparison to historical glaciological observations, the coarser resolution model overestimates accumulation in the southern areas of the GrIS, while the overall SMB is underestimated. These changes relate to the sensitivity of accumulation and melt to the resolution of topography. The GEOS-5 SMB fields contrast with available corresponding atmospheric models simulations from the Coupled Model Intercomparison Project (CMIP5). It is found that only a few of the CMIP5 AGCMs examined provide significant summertime runoff, a dominant feature of the GrIS seasonal cycle. This is a condition that will need to be remedied if potential contributions to future eustatic change from polar ice sheets are to be examined with GCMs.

AGCM↗

Future Climate Change Under SSP Emission Scenarios with GISS-E2.1

This paper presents the response to anthropogenic forcing in the GISS-E2.1 climate models for the 21st century Shared Socioeconomic Pathways emission scenarios within the Coupled Model Intercomparison Project Phase 6 (CMIP6). The experiments were performed using an updated and improved version of the NASA Goddard Institute for Space Studies (GISS) coupled general circulation model that includes two different versions for atmospheric composition: A non-interactive version (NINT) with prescribed composition and a tuned aerosol indirect effect and the One-Moment Aerosol model (OMA) version with fully interactive aerosols which includes a parameterized first indirect aerosol effect on clouds. The effective climate sensitivities are 3.0°C and 2.9°C for the NINT and OMA models, respectively. Each atmospheric version is coupled to two different ocean general circulation models: The GISS ocean model (E2.1-G) and HYCOM (E2.1-H). We describe the global mean responses for all future scenarios and spatial patterns of change for surface air temperature and precipitation for four of the marker scenarios: SSP1-2.6, SSP2-4.5, SSP4-6.0, and SSP5-8.5. By 2100, global mean warming ranges from 1.5°C to 5.2°C relative to 1,850–1,880 mean temperature. Two high-mitigation scenarios SSP1-1.9 and SSP1-2.6 limit the surface warming to below 2°C by the end of the 21st century, except for the NINT E2.1-H model that simulates 2.2°C of surface warming. For the high emission scenario SSP5-8.5, the range is 4.6–5.2°C at 2100. Due to about 15% larger effective climate sensitivity and stronger transient climate response in both NINT and OMA CMIP6 models compared to CMIP5 versions, there is a stronger warming by 2100 in the SSP emission scenarios than in the comparable Representative Concentration Pathway (RCP) scenarios in CMIP5. Changes in sea ice area are highly correlated to global mean surface air temperature anomalies and show steep declines in both hemispheres, with the largest sea ice area decreases occurring during September in the Northern Hemisphere in both E2.1-G (−1.21 × 106 km2/°C) and E2.1-H models (−0.94 × 106 km2/°C). Both coupled models project decreases in the Atlantic overturning stream function by 2100. The largest decrease of 56%–65% in the 21st century overturning stream function is produced in the warmest scenario SSP5-8.5 in the E2.1-G model, comparable to the reduction in the corresponding CMIP5 GISS-E2 RCP8.5 simulation. Both low-end scenarios SSP1-1.9 and SSP1-2.6 also simulate substantial reductions of the overturning (9%–37%) with slow recovery of about 10% by the end of the 21st century (relative to the maximum decrease at the middle of the 21st century).

climate models↗