Global temperature leaped more than 0.4 degrees C (0.7 degrees F) during the past two years, the 12-month average peaking in August 2024 at +1.6 degrees C relative to the temperature at the beginning of last century (the 1880-1920 average). This temperature jump was spurred by one of the periodic tropical El Ni & ntilde;o warming events, but many Earth scientists were baffled by the magnitude of the global warming, which was twice as large as expected for the weak 2023-2024 El Ni & ntilde;o. We find that most of the other half of the warming was caused by a restriction on aerosol emissions by ships, which was imposed in 2020 by the International Maritime Organization to combat the effect of aerosol pollutants on human health. Aerosols are small particles that serve as cloud formation nuclei. Their most important effect is to increase the extent and brightness of clouds, which reflect sunlight and have a cooling effect on Earth. When aerosols - and thus clouds - are reduced, Earth is darker and absorbs more sunlight, thus enhancing global warming. Ships are the main aerosol source in the North Pacific and North Atlantic Oceans. We quantify the aerosol effect from the geographical distribution of sunlight reflected by Earth as measured by satellites, with the largest expected and observed effects in the North Pacific and North Atlantic Oceans. We find that aerosol cooling, and thus climate sensitivity, are understated in the best estimate of the United Nations Intergovernmental Panel on Climate Change (IPCC). Global warming caused by reduced ship aerosols will not go away as tropical climate moves into its cool La Ni & ntilde;a phase. Therefore, we expect that global temperature will not fall much below +1.5 degrees C level, instead oscillating near or above that level for the next few years, which will help confirm our interpretation of the sudden global warming. High sea surface temperatures and increasing ocean hotspots will continue, with harmful effects on coral reefs and other ocean life. The largest practical effect on humans today is increase of the frequency and severity of climate extremes. More powerful tropical storms, tornadoes, and thunderstorms, and thus more extreme floods, are driven by high sea surface temperature and a warmer atmosphere that holds more water vapor. Higher global temperature also increases the intensity of heat waves and - at the times and places of dry weather - high temperature increases drought intensity, including "flash droughts" that develop rapidly, even in regions with adequate average rainfall. Polar climate change has the greatest long-term effect on humanity, with impacts accelerated by the jump in global temperature. We find that polar ice melt and freshwater injection onto the North Atlantic Ocean exceed prior estimates and, because of accelerated global warming, the melt will increase. As a result, shutdown of the Atlantic Meridional Overturning Circulation (AMOC) is likely within the next 20-30 years, unless actions are taken to reduce global warming - in contradiction to conclusions of IPCC. If AMOC is allowed to shut down, it will lock in major problems including sea level rise of several meters - thus, we describe AMOC shutdown as the "point of no return." We suggest that an alternative perspective - a complement to the IPCC approach - is needed to assess these issues and actions that are needed to avoid handing young people a dire situation that is out of their control. This alternative approach will make more use of ongoing observations to drive modeling and more use of paleoclimate to test modeling and test our understanding. As of today, the threats of AMOC shutdown and sea level rise are poorly understood, but better observations of polar ocean and ice changes in response to the present accelerated global warming have the potential to greatly improve our understanding.

We present the second generation of ROCKE-3D (Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics), a generalized three-dimensional general circulation model (GCM) for use in Solar System and exoplanetary simulations of rocky planet climates. ROCKE-3D version 2.0 is a descendant of ModelE2.1, the flagship Earth system model of the NASA Goddard Institute for Space Studies (GISS) used in the most recent Intergovernmental Panel on Climate Change (IPCC) assessments. ROCKE-3D is a continuous effort to expand the capabilities of GISS ModelE to handle a broader range of planetary conditions, including different atmospheric planet sizes, gravities, pressures, and rotation rates; more diverse chemistry schemes and atmospheric compositions; diverse ocean and land distributions and topographies; and potential basic biosphere functions. In this release we present updated physics and many more supported configurations which can serve as starting points to simulate the atmospheres of rocky terrestrial planets of interest. Two different radiation schemes are supported, the GISS radiation, valid only for atmospheres similar to that of modern Earth, and SOCRATES, which is more generalized but more computationally expensive. While ROCKE-3D can simulate a very wide range of planetary and atmospheric configurations, we describe here a small subset of them, with the goal of demonstrating the structural capabilities, rather than the scientific breadth, of the model. Three different atmospheric composition options (preindustrial Earth, the aerosol-free and ozone-free atmosphere used in ROCKE-3D 1.0, and an anoxic atmosphere with no aerosols), three ocean configurations (prescribed, Q-flux, and dynamic), and two resolutions are described: the medium resolution (4x5 degrees in latitude and longitude, previously used in ROCKE-3D 1.0) and the fine resolution, which has double the resolution in the atmosphere and 4 times the horizontal and 3 times the vertical resolution in the ocean. Finally, for the land surface hydrology, we have introduced generalized physics for arbitrary topography in the pooling and evaporation of water and river transport of water between grid cells, as well as for the vertical stratification of temperature in dynamic lakes. We quantify how the different component choices affect model results and discuss the strengths and limitations of using each component, together with how one can select which component to use. ROCKE-3D is publicly available, and tutorial sessions are available for the community, greatly facilitating its use by any interested group.

The historical global temperature record is an essential data product for quantifying the variability and change of the Earth system. In recent years, better characterization of observational uncertainty in global and hemispheric trends has become available, but the methodologies are not necessarily applicable to analyses at smaller regional areas, or monthly or seasonal means, where station sparsity and other systematic issues contribute to greater uncertainty. This study presents a gridded uncertainty ensemble of historical surface temperature anomalies from the Goddard Institute for Space Studies (GISS) Surface Temperature (GISTEMP) product. This ensemble characterizes the complex spatial and temporal correlation structure of uncertainty, enabling better uncertainty propagation for climate and applied science in applications of historical temperature products at spatial scales from global to regional and temporal scales from centennial to monthly. This work details the methodology for generating the uncertainty ensemble, presents key statistics of the uncertainty evolution over space and time, and provides best practices for using the uncertainty ensemble in future studies. Summary statistics from the uncertainty ensemble agree well with the previous GISTEMP global uncertainty assessment, providing confidence in both.

Marine oxygen plays a fundamental role in regulating the transfer of organic carbon and nutrients to their dissolved inorganic forms, serving as the terminal electron acceptor for heterotrophic respiration. Oxygen can become limiting to these processes in coastal and open-ocean oxygen deficient zones (ODZs). The maintenance of ODZs depends on the balance between physical processes such as ventilation and biogeochemical processes such as remineralization. These two processes act in opposing directions on ODZs, with ventilation responsible for transporting oxygen from the surface, where it is near saturation with the partial pressure of O-2 in the atmosphere, to deeper oxygen-depleted layers, and remineralization acting to consume oxygen through biogeochemical processes such as microbial degradation throughout the water column. Remineralization is represented in all CMIP6 models, but the actual parameterizations, as well as their magnitude, are widely varying. In this study, we examine the sensitivity to remineralization of the equatorial Pacific ODZ (O-2 <60 mu mol/kg) in a model; the NASA GISS Ocean Biogeochemical Model (NOBM-G) embedded in the NASA-GISS coupled atmosphere-ocean model. We find that increasing the remineralization rate shoals the ODZ onset depth (i.e. the regionally averaged upper bound of the ODZ), and decreases ODZ volume and regional-averaged thickness. Changes in biological consumption and vertical convergence of oxygen are identified as the primary contributors to changes in ODZ onset depth. ODZ thickness is mainly influenced by the shoaling of the bottom boundary of the ODZ, which decreases with maximum remineralization rate due to decreasing oxygen consumption and increasing horizontal oxygen convergence in deep waters. On the other hand, vertically-averaged ODZ area has a more complex, non-monotonic relationship with maximum remineralization rate. While these findings suggest an important role for remineralization in determining the shape of the ODZ in our model, the relative importance of remineralization vs. other physical parameterizations remains to be established. While our results reflect the structure of our ocean biogeochemical model, the relationship of remineralization to ODZ characteristics in other models should be examined to better inform model inter-comparisons.

Abstract. While paleoclimate simulations have been a priority for Earth system modelers over the past three decades, little attention has been paid to the period between the mid-Holocene and the Last Millennium, although this is an important period for the emergence of complex societies. Here, we consider the climate of 2500 BP (550 BCE), a period when compared to late preindustrial time, greenhouse gas concentrations were slightly lower, and orbital forcing led to a stronger seasonal cycle in high latitude insolation. To capture the influence of land cover on climate, we asynchronously coupled the NASA GISS ModelE Earth system model with the LPJ-LMfire dynamic global vegetation model. We simulated global climate and assessed our results in the context of independent paleoclimate reconstructions. We also explored a set of combinations of model performance parameters (bias and variability) and demonstrated their importance for the asynchronous coupling framework. The coupled model system shows substantial vegetation albedo feedback to climate. In the absence of a bias correction, while driving LPJ-LMfire in the coupling process, ModelE drifts towards colder conditions in the high latitudes of the Northern Hemisphere in response to land cover simulated by LPJ-LMfire. A regional precipitation response is also prominent in the various combinations of the coupled model system, with a substantial intensification of the Summer Indian Monsoon and a drying pattern over Europe. Evaluation of the simulated climate against reconstructions of temperature from multiple proxies and the isotopic composition of precipitation (δ18Op) from speleothems demonstrated the skill of ModelE in simulating past climate. A regional analysis of the simulated vegetation-climate response further confirmed the validity of this approach. The coupled model system is sensitive to the representation of shrubs and this land cover type requires particular attention as a potentially important driver of climate in regions where shrubs are abundant. Our results further demonstrate the importance of bias correction in coupled paleoclimate simulations.

Improved knowledge of glacial-to-interglacial global temperature change yields Charney (fast-feedback) equilibrium climate sensitivity 1.2 ± 0.3°C (2σ) per W/m2, which is 4.8°C ± 1.2°C for doubled CO2. Consistent analysis of temperature over the full Cenozoic era—including ‘slow’ feedbacks by ice sheets and trace gases—supports this sensitivity and implies that CO2 was 300–350 ppm in the Pliocene and about 450 ppm at transition to a nearly ice-free planet, exposing unrealistic lethargy of ice sheet models. Equilibrium global warming for today’s GHG amount is 10°C, which is reduced to 8°C by today’s human-made aerosols. Equilibrium warming is not ‘committed’ warming; rapid phaseout of GHG emissions would prevent most equilibrium warming from occurring. However, decline of aerosol emissions since 2010 should increase the 1970–2010 global warming rate of 0.18°C per decade to a post-2010 rate of at least 0.27°C per decade. Thus, under the present geopolitical approach to GHG emissions, global warming will exceed 1.5°C in the 2020s and 2°C before 2050. Impacts on people and nature will accelerate as global warming increases hydrologic (weather) extremes. The enormity of consequences demands a return to Holocene-level global temperature. Required actions include: (1) a global increasing price on GHG emissions accompanied by development of abundant, affordable, dispatchable clean energy, (2) East-West cooperation in a way that accommodates developing world needs, and (3) intervention with Earth’s radiation imbalance to phase down today’s massive human-made ‘geo-transformation’ of Earth’s climate. Current political crises present an opportunity for reset, especially if young people can grasp their situation.

The historical global temperature record is an essential data product for quantifying the variability and change of the Earth system. In recent years, better characterization of observational uncertainty in global and hemispheric trends has become available, but the methodologies are not necessarily applicable to analyses at smaller regional areas, or monthly means, where station sparsity and other systematic issues contribute to greater uncertainty.This work details a gridded uncertainty ensemble of historical temperature anomalies from the Goddard Institute for Space Studies (GISS) Surface Temperature product (GISTEMP) product. This ensemble characterizes the complex spatial and temporal correlation structure of uncertainty in gridded historical temperature, enabling proper uncertainty propagation for climate and social science at regional and monthly scales. This work details the methodology for generating the uncertainty ensemble, key statistics of the uncertainty evolution over space and time, and provides best practices for using the uncertainty ensemble in future studies. Summary statistics from the uncertainty ensemble are in good agreement with production GISTEMP.Two applications of the uncertainty ensemble are also presented. First, the warmest year on record is shown to most likely be 2016 with a 53.2% chance and 2020 as the second most likely with a 44.4% chance. Second, it is shown that the arctic is warming 2.5 - 5 times faster than the globe, significantly faster than the regularly quoted twice as fast.

Abstract 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).

Abstract Climate change will affect both the mean state and seasonality of marine physical and biogeochemical properties, with important implications for the oceanic sink of atmospheric CO2. Here, we investigate the seasonal cycle of the air‐sea exchange of CO2 and pCO2,sw (surface seawater pCO2) and their long term changes using the CMIP6 submission of the NASA‐GISS modelE (GISS‐E2.1‐G). In comparison to the CMIP5 submission (GISS‐E2‐R), we find that on the global scale, the seasonal cycles of the CO2 flux and NPP have improved, while the seasonal cycles of dissolved inorganic carbon (DIC), alkalinity, and macronutrients have deteriorated. Moreover, for all ocean biogeochemistry fields, changes in skill between E2.1‐G and E2‐R display large regional variability. For E2.1‐G, we find similar modeled and observed CO2 flux seasonal cycles in the subtropical gyres, where seasonal anomalies of pCO2,sw and the flux are temperature‐driven, and the Southern Ocean, where anomalies are DIC‐driven. Biases in these seasonal cycles are largest in the subpolar and equatorial regions, driven by a combination of biases in temperature, DIC, alkalinity, and wind speed. When comparing the historical simulation to a simulation with an idealized increase in atmospheric pCO2, we find that the seasonal amplitudes of the CO2 flux and pCO2,sw generally increase. These changes are produced by increases in the sensitivity of pCO2,sw to its respective drivers. These findings are consistent with the notion that the seasonality of pCO2,sw is expected to increase due to the increase of atmospheric pCO2, with changes in the seasonality of temperature, DIC, and alkalinity having secondary influences.

AbstractSimulations 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.

Much attention has been given to the climate dynamics and habitable boundaries of synchronously rotating planets around low mass stars. However, other rotational states are possible, particularly when higher eccentricity orbits can be maintained in a system, including spin-orbit resonant configurations. Additionally, the oscillating strain as a planet moves from periastron to apoastron results in friction and tidal heating, which can be an important energy source. Here, we simulate the climate of ocean-covered planets near the inner edge of the habitable zone around M to solar stars with ROCKE-3D, and leverage the planetary evolution software package, VPLanet, to calculate tidal heating rates for Earth-sized planets orbiting 2600 K and 3000 K stars. This study is the first to use a 3-D General Circulation Model that implements tidal heating to investigate habitability for multiple resonant states. We find that in the absence of tidal heating, the resonant state has little impact on the inner edge, because for a given stellar flux, higher-order states tend to be warmer than synchronous rotators, but for a given temperature, have drier upper atmospheres. However, when strong tidal heating is present, the rotational component implies a strong dependence of habitable conditions on the system evolution and rotational state. Since tidal and stellar heating both decrease with orbital distance, this results in a compact orbital width separating temperate and uninhabitable climates. We summarize these results and also compare ROCKE-3D to previously published simulations of the inner edge that used a modified version of the NCAR CAM4 model.

Understanding the climate of terrestrial planetary atmospheres has been increasingly the focus of research worldwide, in light of the increasing amount of rocky planet discoveries orbiting other stars in or near their habitable zone. Here we present simulations with the new version of the 3D climate model ROCKE-3D, whose version 2.0 will soon become publicly available. A wide range of configurations will be supported, compared to a handful ones in its predecessor, version 1.0 (Way et al., 2017). These include two model resolutions (4x5 and 2x2.5), two radiation schemes (GISS and SOCRATES), three atmospheric configurations (Earth-like, Earth-like without O3 and aerosols, and N2-dominated), and three ocean setups (prescribed sea-surface temperatures and ice cover, q-flux, and dynamic). Simulations of all different configuration combinations have been performed and will become available for use by the community. Key results will be presented across those configurations, together of the role of the structural uncertainty in model setup in the resulting climate calculated by the model.

Abstract We present results from the NASA GISS ModelE2.1‐G‐CC Earth System Model with coupled climate‐carbon cycle simulations that were submitted to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) Coupled Climate‐Carbon Cycle MIP (C4MIP). Atmospheric CO2 concentration and carbon budgets for the land and ocean in the historical simulations were generally consistent with observations. Low simulated atmospheric CO2 concentrations during 1850–1950 were due to excess uptake from prescribed land cover change, which erroneously replaced arid shrublands with higher biomass crops, and assumed high 2004 LAI values in vegetated lands throughout the historical simulation. At the end of the historical period, slightly higher simulated CO2 than observed resulted from the land being an insufficient net carbon sink, despite the net effect of CO2 fertilization and warming‐induced increases to leaf photosynthetic capacity. The global ocean carbon uptake agreed well with the observations with the largest discrepancies in the low latitudes. Future climate projection at 2091–2100 agreed with CMIP5 models in the northward shift, of temperate deciduous forest climate and expansion across Eurasia along 60°N latitude, and dramatic regional biome shifts from drying and warming in continental Europe. Carbon feedback parameters were largely similar to the CMIP5 model ensemble. For our model, the variation of land feedback parameters within the uncertainty arises from the fertilization feedback being less sensitive due to lack of increased vegetation growth, and the comparably more negative ocean carbon‐climate feedback is due to the large slowdown of the Atlantic overturning circulation.

AbstractThis 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 × CO2 is slightly higher than previously at 2.7–3.1°C (depending on version) and is a result of lower CO2 radiative forcing and stronger positive feedbacks.

This paper is dedicated to Wally Broecker, the "father of global warming", whose inquisitive mind has stimulated much of the world"s research aimed at understanding global climate.