This paper describes the GISS‐E2.1 contribution to the Coupled Model Intercomparison Project, Phase 6 (CMIP6). This model version differs from the predecessor model (GISS‐E2) chiefly due to parameterization improvements to the atmospheric and ocean model components, while keeping atmospheric resolution the same. Model skill when compared to modern era climatologies is significantly higher than in previous versions. Additionally, updates in forcings have a material impact on the results. In particular, there have been specific improvements in representations of modes of variability (such as the Madden‐Julian Oscillation and other modes in the Pacific) and significant improvements in the simulation of the climate of the Southern Oceans, including sea ice. The effective climate sensitivity to 2 × CO₂ is slightly higher than previously at 2.7–3.1°C (depending on version) and is a result of lower CO₂ radiative forcing and stronger positive feedbacks.

The Earth's climate is rapidly changing. Over the past centuries, aerosols, via their ability to absorb or scatter solar radiation and alter clouds, played an important role in counterbalancing some of the greenhouse gas (GHG) caused global warming. The multicentury anthropogenic aerosol cooling effect prevented present‐day climate from reaching even higher surface air temperatures and subsequent more dramatic climate impacts. Trends in aerosol concentrations and optical depth show that in many polluted regions such as Europe and the United States, aerosol precursor emissions decreased back to levels of the 1950s. More recent polluting countries such as China may have reached a turning point in recent years as well, while India still follows an upward trend. Here we study aerosol trends in the Coupled Model Intercomparison Project Phase 6 (CMIP6) simulations of the GISS ModelE2.1 climate model using a fully coupled atmosphere composition configuration, including interactive gas‐phase chemistry and either an aerosol microphysical (MATRIX) or a mass‐based (One‐Moment Aerosol, OMA) aerosol module. Results show that whether global aerosol radiative forcing is already declining depends on the aerosol scheme used. Using the aerosol microphysical scheme, where the aerosol system reacts more strongly to the trend in sulfur dioxide (SO₂) emissions, global peak direct aerosol forcing was reached in the 1980s, whereas the mass‐based scheme simulates peak direct aerosol forcing around 2010.

Simulations of the CMIP6 historical period 1850–2014, characterized by the emergence of anthropogenic climate drivers like greenhouse gases, are presented for different configurations of the NASA Goddard Institute for Space Studies (GISS) Earth System ModelE2.1. The GISS‐E2.1 ensembles are more sensitive to greenhouse gas forcing than their CMIP5 predecessors (GISS‐E2) but warm less during recent decades due to a forcing reduction that is attributed to greater longwave opacity in the GISS‐E2.1 pre‐industrial simulations. This results in an atmosphere less sensitive to increases in opacity from rising greenhouse gas concentrations, demonstrating the importance of the base climatology to forcing and forced climate trends. Most model versions match observed temperature trends since 1979 from the ocean to the stratosphere. The choice of ocean model is important to the transient climate response, as found previously in CMIP5 GISS‐E2: the model that more efficiently exports heat to the deep ocean shows a smaller rise in tropospheric temperature. Model sea level rise over the historical period is traced to excessive drawdown of aquifers to meet irrigation demand with a smaller contribution from thermal expansion. This shows how fully coupled models can provide indirect observational constraints upon forcing, in this case, constraining irrigation rates with observed sea level changes. The overall agreement of GISS‐E2.1 with observed trends is familiar from evaluation of its predecessors, as is the conclusion that these trends are almost entirely anthropogenic in origin.

We examine the anthropogenically forced climate response for the 21st century representative concentration pathway (RCP) emission scenarios and their extensions for the period 2101–2500. The experiments were performed with ModelE2, a new version of the NASA Goddard Institute for Space Sciences (GISS) coupled general circulation model that includes three different versions for the atmospheric composition components: a noninteractive version (NINT) with prescribed composition and a tuned aerosol indirect effect (AIE), the TCAD version with fully interactive aerosols, whole‐atmosphere chemistry, and the tuned AIE, and the TCADI version which further includes a parameterized first indirect aerosol effect on clouds. Each atmospheric version is coupled to two different ocean general circulation models: the Russell ocean model (GISS‐E2‐R) and HYCOM (GISS‐E2‐H). By 2100, global mean warming in the RCP scenarios ranges from 1.0 to 4.5°C relative to 1850–1860 mean temperature in the historical simulations. In the RCP2.6 scenario, the surface warming in all simulations stays below a 2°C threshold at the end of the 21st century. For RCP8.5, the range is 3.5–4.5°C at 2100. Decadally averaged sea ice area changes are highly correlated to global mean surface air temperature anomalies and show steep declines in both hemispheres, with a larger sensitivity during winter months. By the year 2500, there are complete recoveries of the globally averaged surface air temperature for all versions of the GISS climate model in the low‐forcing scenario RCP2.6. TCADI simulations show enhanced warming due to greater sensitivity to CO₂, aerosol effects, and greater methane feedbacks, and recovery is much slower in RCP2.6 than with the NINT and TCAD versions. All coupled models have decreases in the Atlantic overturning stream function by 2100. In RCP2.6, there is a complete recovery of the Atlantic overturning stream function by the year 2500 while with scenario RCP8.5, the E2‐R climate model produces a complete shutdown of deep water formation in the North Atlantic.

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 × 10⁶ km²/°C) and E2.1‐H models (−0.94 × 10⁶ km²/°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).