GISS-E2.1: Configurations and Climatology.

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

Habitable Climate Scenarios for Proxima Centauri b with a Dynamic Ocean.

The nearby exoplanet Proxima Centauri b will be a prime future target for characterization, despite questions about its retention of water. Climate models with static oceans suggest that Proxima b could harbor a small dayside surface ocean despite its weak instellation. We present the first climate simulations of Proxima b with a dynamic ocean. We find that an ocean-covered Proxima b could have a much broader area of surface liquid water but at much colder temperatures than previously suggested, due to ocean heat transport and/or depression of the freezing point by salinity. Elevated greenhouse gas concentrations do not necessarily produce more open ocean because of dynamical regime transitions between a state with an equatorial Rossby-Kelvin wave pattern and a state with a day-night circulation. For an evolutionary path leading to a highly saline ocean, Proxima b could be an inhabited, mostly open ocean planet with halophilic life. A freshwater ocean produces a smaller liquid region than does an Earth salinity ocean. An ocean planet in 3:2 spin-orbit resonance has a permanent tropical waterbelt for moderate eccentricity. A larger versus smaller area of surface liquid water for similar equilibrium temperature may be distinguishable by using the amplitude of the thermal phase curve. Simulations of Proxima Centauri b may be a model for the habitability of weakly irradiated planets orbiting slightly cooler or warmer stars, for example, in the TRAPPIST-1, LHS 1140, GJ 273, and GJ 3293 systems.

Albedos, Equilibrium Temperatures, and Surface Temperatures of Habitable Planets.

The potential habitability of known exoplanets is often categorized by a nominal equilibrium temperature assuming a Bond albedo of either ∼0.3, similar to Earth, or 0. As an indicator of habitability, this leaves much to be desired, because albedos of other planets can be very different, and because surface temperature exceeds equilibrium temperature due to the atmospheric greenhouse effect. We use an ensemble of general circulation model simulations to show that for a range of habitable planets, much of the variability of Bond albedo, equilibrium temperature and even surface temperature can be predicted with useful accuracy from incident stellar flux and stellar temperature, two known parameters for every confirmed exoplanet. Earth's Bond albedo is near the minimum possible for habitable planets orbiting G stars, because of increasing contributions from clouds and sea ice/snow at higher and lower instellations, respectively. For habitable M star planets, Bond albedo is usually lower than Earth's because of near-IR H2O absorption, except at high instellation where clouds are important. We apply relationships derived from this behavior to several known exoplanets to derive zeroth-order estimates of their potential habitability. More expansive multivariate statistical models that include currently non-observable parameters show that greenhouse gas variations produce significant variance in albedo and surface temperature, while increasing length of day and land fraction decrease surface temperature; insights for other parameters are limited by our sampling. We discuss how emerging information from global climate models might resolve some degeneracies and help focus scarce observing resources on the most promising planets.

Climates of Warm Earth-like Planets I: 3-D Model Simulations.

We present a large ensemble of simulations of an Earth-like world with increasing insolation and rotation rate. Unlike previous work utilizing idealized aquaplanet configurations we focus our simulations on modern Earth-like topography. The orbital period is the same as modern Earth, but with zero obliquity and eccentricity. The atmosphere is 1 bar N2-dominated with CO2=400 ppmv and CH4=1 ppmv. The simulations include two types of oceans; one without ocean heat transport (OHT) between grid cells as has been commonly used in the exoplanet literature, while the other is a fully coupled dynamic bathtub type ocean. The dynamical regime transitions that occur as day length increases induce climate feedbacks producing cooler temperatures, first via the reduction of water vapor with increasing rotation period despite decreasing shortwave cooling by clouds, and then via decreasing water vapor and increasing shortwave cloud cooling, except at the highest insolations. Simulations without OHT are more sensitive to insolation changes for fast rotations while slower rotations are relatively insensitive to ocean choice. OHT runs with faster rotations tend to be similar with gyres transporting heat poleward making them warmer than those without OHT. For slower rotations OHT is directed equator-ward and no high latitude gyres are apparent. Uncertainties in cloud parameterization preclude a precise determination of habitability but do not affect robust aspects of exoplanet climate sensitivity. This is the first paper in a series that will investigate aspects of habitability in the simulations presented herein. The datasets from this study are opensource and publicly available.

An Improved Convective Ice Parameterization for the NASA GISS Global Climate Model and Impacts on Cloud Ice Simulation.

Partitioning of convective ice into precipitating and detrained condensate presents a challenge for GCMs since partitioning depends on the strength and microphysics of the convective updraft. It is an important issue because detrainment of ice from updrafts influences the development of stratiform anvils, impacts radiation, and can affect GCM climate sensitivity. Recent studies have shown that the CMIP5 configurations of the Goddard Institute for Space Studies (GISS) GCM simulated upper-tropospheric ice water content (IWC) that exceeded an estimated upper bound by a factor of 2. Partly in response to this bias, a new GCM parameterization of convective cloud ice has been developed that incorporates new ice particle fall speeds and convective outflow particle size distributions (PSDs) from the NASA African Monsoon Multidisciplinary Analyses (NAMMA), NASA Tropical Composition, Cloud and Climate Coupling (TC4), DOE ARM-NASA Midlatitude Continental Convective Clouds Experiment (MC3E), and DOE ARM Small Particles in Cirrus (SPARTICUS) field campaigns. The new parameterization assumes a normalized gamma PSD with two novel developments: no explicit assumption for particle habit in the calculation of mass distributions, and a formulation for translating ice particle fall speeds as a function of maximum diameter into fall speeds as a function of melted-equivalent diameter. Two parameters (particle volume- and projected area-weighted equivalent diameter) are diagnosed as a function of temperature and IWC in the convective plume, and these parameters constrain the shape and scale of the normalized gamma PSD. The diagnosed fall speeds and PSDs are combined with the GCM's parameterized convective updraft vertical velocity to partition convective updraft condensate into precipitating and detrained components. A 5-yr prescribed sea surface temperature GCM simulation shows a 30%-50% decrease in upper-tropospheric deep convective IWC, bringing the tropical and global mean ice water path into closer agreement with CloudSat best estimates.

The relationship between boundary layer stability and cloud cover in the post-cold frontal region.

Using NASA-Aqua MODIS and AIRS data, the relationship between low-level cloud cover (cloud top below the 700 hPa level) and boundary layer stability is explored in post-cold frontal conditions. A linear relationship is found between seasonal cloud cover and two separate measures of inversion strength, the lower tropospheric stability (LTS) and the estimated inversion strength (EIS), for two specific regions in the north Atlantic and Pacific in quiescent and weakly subsiding conditions. The relationship barely changes when considering dynamically active and subsiding post-cold frontal conditions for the same regions. To explore the generality of this result and increase sample size, cold front centered composites of cloud cover and stability are constructed. The northern and southern hemisphere seasonal cloud cover and stability distributions in the post-cold frontal regions are then compared. A fairly good correlation between cloud cover and EIS is found in both hemispheres across all seasons, suggesting that a linear relationship between cloud cover and inversion strength proposed for quiescent conditions exists also in more dynamically active subsiding post-cold frontal conditions. However, for a given season and hemisphere, the correlation between cloud cover and EIS degrades in post-cold frontal regions, especially in the northern hemisphere. At these scales, other large scale factors tend to correlate better with cloud cover.

Responses of Tropical Ocean Clouds and Precipitation to the Large-Scale Circulation:Atmospheric-Water-Budget-Related Phase Space and Dynamical Regimes.

An atmospheric-water-budget-related phase space is constructed with the tendency terms related to dynamical convergence (QCON ≡ -Q∇ · V) and moisture advection (QADV ≡ -V · ∇Q) in the water budget equation. Over the tropical oceans, QCON accounts for large-scale dynamical conditions related to conditional instability, and QADV accounts for conditions related to lower-tropospheric moisture gradient. Two reanalysis products [MERRA and ERA-Interim (ERAi)] are used to calculate QCON and QADV. Using the phase space as a reference frame, the Moderate Resolution Imaging Spectroradiometer (MODIS) cloud-top pressure (CTP) and cloud optical depth (COD) are used to evaluate simulated clouds in the GISS-E2 general circulation model. In regimes of divergence over the tropical oceans, moist advection yields frequent high- to midlevel medium-thickness to thick clouds associated with moderate stratiform precipitation, while dry advection yields low-level thin clouds associated with shallow convection with lowered cloud tops. In regimes with convergence, moist and dry advection modulate the relative abundance of high-level thick clouds and low-level thin to medium-thickness clouds. GISS-E2 qualitatively reproduces the cloud property dependence on moisture budget tendencies in regimes of convergence but with larger COD compared to MODIS. Low-level thick clouds in GISS-E2 are the most frequent in regimes of near-zero convergence and moist advection instead of those of large-scale divergence. Compared to the Global Precipitation Climatology Project product, MERRA, ERAi, and GISS-E2 have more rain in regimes with deep convection and less rain in regimes with shallow convection.

Was Venus the First Habitable World of our Solar System?

Present-day Venus is an inhospitable place with surface temperatures approaching 750K and an atmosphere 90 times as thick as Earth's. Billions of years ago the picture may have been very different. We have created a suite of 3-D climate simulations using topographic data from the Magellan mission, solar spectral irradiance estimates for 2.9 and 0.715 Gya, present-day Venus orbital parameters, an ocean volume consistent with current theory, and an atmospheric composition estimated for early Venus. Using these parameters we find that such a world could have had moderate temperatures if Venus had a rotation period slower than ~16 Earth days, despite an incident solar flux 46-70% higher than Earth receives. At its current rotation period, Venus's climate could have remained habitable until at least 715 million years ago. These results demonstrate the role rotation and topography play in understanding the climatic history of Venus-like exoplanets discovered in the present epoch.

Strong temporal variation over one Saturnian year: from Voyager to Cassini.

Here we report the combined spacecraft observations of Saturn acquired over one Saturnian year (~29.5 Earth years), from the Voyager encounters (1980-81) to the new Cassini reconnaissance (2009-10). The combined observations reveal a strong temporal increase of tropic temperature (~10 Kelvins) around the tropopause of Saturn (i.e., 50 mbar), which is stronger than the seasonal variability (~a few Kelvins). We also provide the first estimate of the zonal winds at 750 mbar, which is close to the zonal winds at 2000 mbar. The quasi-consistency of zonal winds between these two levels provides observational support to a numerical suggestion inferring that the zonal winds at pressures greater than 500 mbar do not vary significantly with depth. Furthermore, the temporal variation of zonal winds decreases its magnitude with depth, implying that the relatively deep zonal winds are stable with time.

Dynamics of Saturn's south polar vortex.

The camera onboard the Cassini spacecraft has allowed us to observe many of Saturn's cloud features. We present observations of Saturn's south polar vortex (SPV) showing that it shares some properties with terrestrial hurricanes: cyclonic circulation, warm central region (the eye) surrounded by a ring of high clouds (the eye wall), and convective clouds outside the eye. The polar location and the absence of an ocean are major differences. It also shares properties with the polar vortices on Venus, such as polar location, cyclonic circulation, warm center, and long lifetime, but the Venus vortices have cold collars and are not associated with convective clouds. The SPV's combination of properties is unique among vortices in the solar system.

Imaging of Titan from the Cassini spacecraft.

Titan, the largest moon of Saturn, is the only satellite in the Solar System with a substantial atmosphere. The atmosphere is poorly understood and obscures the surface, leading to intense speculation about Titan's nature. Here we present observations of Titan from the imaging science experiment onboard the Cassini spacecraft that address some of these issues. The images reveal intricate surface albedo features that suggest aeolian, tectonic and fluvial processes; they also show a few circular features that could be impact structures. These observations imply that substantial surface modification has occurred over Titan's history. We have not directly detected liquids on the surface to date. Convective clouds are found to be common near the south pole, and the motion of mid-latitude clouds consistently indicates eastward winds, from which we infer that the troposphere is rotating faster than the surface. A detached haze at an altitude of 500 km is 150-200 km higher than that observed by Voyager, and more tenuous haze layers are also resolved.

Cassini imaging of Jupiter's atmosphere, satellites, and rings.

The Cassini Imaging Science Subsystem acquired about 26,000 images of the Jupiter system as the spacecraft encountered the giant planet en route to Saturn. We report findings on Jupiter's zonal winds, convective storms, low-latitude upper troposphere, polar stratosphere, and northern aurora. We also describe previously unseen emissions arising from Io and Europa in eclipse, a giant volcanic plume over Io's north pole, disk-resolved images of the satellite Himalia, circumstantial evidence for a causal relation between the satellites Metis and Adrastea and the main jovian ring, and information on the nature of the ring particles.

Evidence for strengthening of the tropical general circulation in the 1990s.

Satellite observations suggest that the thermal radiation emitted by Earth to space increased by more than 5 watts per square meter, while reflected sunlight decreased by less than 2 watts per square meter, in the tropics over the period 1985-2000, with most of the increase occurring after 1990. By analyzing temporal changes in the frequency of occurrence of emitted thermal and reflected solar fluxes, the effects of El Niño-Southern Oscillation are minimized, and an independent longer-time-scale variation of the radiation budget is identified. Similar analyses of upper tropospheric humidity, cloud amount, surface air temperature, and vertical velocity confirm that these flux changes are associated with a decadal-time-scale strengthening of the tropical Hadley and Walker circulations. Equatorial convective regions have intensified in upward motion and moistened, while both the equatorial and subtropical subsidence regions have become drier and less cloudy.