CO2 Ocean Bistability on Terrestrial Exoplanets.

Cycling of carbon dioxide between the atmosphere and interior of rocky planets can stabilize global climate and enable planetary surface temperatures above freezing over geologic time. However, variations in global carbon budget and unstable feedback cycles between planetary sub-systems may destabilize the climate of rocky exoplanets toward regimes unknown in the Solar System. Here, we perform clear-sky atmospheric radiative transfer and surface weathering simulations to probe the stability of climate equilibria for rocky, ocean-bearing exoplanets at instellations relevant for planetary systems in the outer regions of the circumstellar habitable zone. Our simulations suggest that planets orbiting G- and F-type stars (but not M-type stars) may display bistability between an Earth-like climate state with efficient carbon sequestration and an alternative stable climate equilibrium where CO2 condenses at the surface and forms a blanket of either clathrate hydrate or liquid CO2. At increasing instellation and with ineffective weathering, the latter state oscillates between cool, surface CO2-condensing and hot, non-condensing climates. CO2 bistable climates may emerge early in planetary history and remain stable for billions of years. The carbon dioxide-condensing climates follow an opposite trend in pCO2 versus instellation compared to the weathering-stabilized planet population, suggesting the possibility of observational discrimination between these distinct climate categories.

Vertically Resolved Magma Ocean-Protoatmosphere Evolution: H2, H2O, CO2, CH4, CO, O2, and N2 as Primary Absorbers.

The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositional and thermal scenarios of early planetary evolution, even though these are of prime importance to aid astronomical inferences on the environmental context and geological history of extrasolar planets. Here, we present a coupled numerical framework that links an evolutionary, vertically resolved model of the planetary silicate mantle with a radiative-convective model of the atmosphere. Using this method, we investigate the early evolution of idealized Earth-sized rocky planets with end-member, clear-sky atmospheres dominated by either H2, H2O, CO2, CH4, CO, O2, or N2. We find central metrics of early planetary evolution, such as energy gradient, sequence of mantle solidification, surface pressure, or vertical stratification of the atmosphere, to be intimately controlled by the dominant volatile and outgassing history of the planet. Thermal sequences fall into three general classes with increasing cooling timescale: CO, N2, and O2 with minimal effect, H2O, CO2, and CH4 with intermediate influence, and H2 with several orders of magnitude increase in solidification time and atmosphere vertical stratification. Our numerical experiments exemplify the capabilities of the presented modeling framework and link the interior and atmospheric evolution of rocky exoplanets with multiwavelength astronomical observations.

Dynamics of atmospheres with a non-dilute condensible component.

The diversity of characteristics for the host of recently discovered exoplanets opens up a great deal of fertile new territory for geophysical fluid dynamics, particularly when the fluid flow is coupled to novel thermodynamics, radiative transfer or chemistry. In this paper, we survey one of these new areas-the climate dynamics of atmospheres with a non-dilute condensible component, defined as the situation in which a condensible component of the atmosphere makes up a substantial fraction of the atmospheric mass within some layer. Non-dilute dynamics can occur for a wide range of condensibles, generically applying near both the inner and the outer edges of the conventional habitable zone and in connection with runaway greenhouse phenomena. It also applies in a wide variety of other planetary circumstances. We first present a number of analytical results developing some key features of non-dilute atmospheres, and then show how some of these features are manifest in simulations with a general circulation model adapted to handle non-dilute atmospheres. We find that non-dilute atmospheres have weak horizontal temperature gradients even for rapidly rotating planets, and that their circulations are largely barotropic. The relative humidity of the condensible component tends towards 100% as the atmosphere becomes more non-dilute, which has important implications for runaway greenhouse thresholds. Non-dilute atmospheres exhibit a number of interesting organized convection features, for which there is not yet any adequate theoretical understanding.

The effect of host star spectral energy distribution and ice-albedo feedback on the climate of extrasolar planets.

Planetary climate can be affected by the interaction of the host star spectral energy distribution with the wavelength-dependent reflectivity of ice and snow. In this study, we explored this effect with a one-dimensional (1-D), line-by-line, radiative transfer model to calculate broadband planetary albedos as input to a seasonally varying, 1-D energy balance climate model. A three-dimensional (3-D) general circulation model was also used to explore the atmosphere's response to changes in incoming stellar radiation, or instellation, and surface albedo. Using this hierarchy of models, we simulated planets covered by ocean, land, and water-ice of varying grain size, with incident radiation from stars of different spectral types. Terrestrial planets orbiting stars with higher near-UV radiation exhibited a stronger ice-albedo feedback. We found that ice extent was much greater on a planet orbiting an F-dwarf star than on a planet orbiting a G-dwarf star at an equivalent flux distance, and that ice-covered conditions occurred on an F-dwarf planet with only a 2% reduction in instellation relative to the present instellation on Earth, assuming fixed CO(2) (present atmospheric level on Earth). A similar planet orbiting the Sun at an equivalent flux distance required an 8% reduction in instellation, while a planet orbiting an M-dwarf star required an additional 19% reduction in instellation to become ice-covered, equivalent to 73% of the modern solar constant. The reduction in instellation must be larger for planets orbiting cooler stars due in large part to the stronger absorption of longer-wavelength radiation by icy surfaces on these planets in addition to stronger absorption by water vapor and CO(2) in their atmospheres, which provides increased downwelling longwave radiation. Lowering the IR and visible-band surface ice and snow albedos for an M-dwarf planet increased the planet's climate stability against changes in instellation and slowed the descent into global ice coverage. The surface ice-albedo feedback effect becomes less important at the outer edge of the habitable zone, where atmospheric CO(2) could be expected to be high such that it maintains clement conditions for surface liquid water. We showed that ∼3-10 bar of CO(2) will entirely mask the climatic effect of ice and snow, leaving the outer limits of the habitable zone unaffected by the spectral dependence of water ice and snow albedo. However, less CO(2) is needed to maintain open water for a planet orbiting an M-dwarf star than would be the case for hotter main-sequence stars.

The dynamics behind Titan's methane clouds.

We present results of an axisymmetric global circulation model of Titan with a simplified suite of atmospheric physics forced by seasonally varying insolation. The recent discovery of midlatitude tropospheric clouds on Titan has caused much excitement about the roles of surface sources of methane and the global circulation in forming clouds. Although localized surface sources, such as methane geysers or "cryovolcanoes," have been invoked to explain these clouds, we find in this work that clouds appear in regions of convergence by the mean meridional circulation and over the poles during solstices, where the solar forcing reaches its seasonal maximum. Other regions are inhibited from forming clouds because of dynamical transports of methane and strong subsidence. We find that for a variety of moist regimes, i.e., with the effect of methane thermodynamics included, the observed cloud features can be explained by the large-scale dynamics of the atmosphere. Clouds at the solsticial pole are found to be a robust feature of Titan's dynamics, whereas isolated midlatitude clouds are present exclusively in a variety of moist dynamical regimes. In all cases, even without including methane thermodynamics, our model ceases to produce polar clouds approximately 4-6 terrestrial years after solstices.

High levels of atmospheric carbon dioxide necessary for the termination of global glaciation.

The possibility that the Earth suffered episodes of global glaciation as recently as the Neoproterozoic period, between about 900 and 543 million years ago, has been widely discussed. Termination of such 'hard snowball Earth' climate states has been proposed to proceed from accumulation of carbon dioxide in the atmosphere. Many salient aspects of the snowball scenario depend critically on the threshold of atmospheric carbon dioxide concentrations needed to trigger deglaciation. Here I present simulations with a general circulation model, using elevated carbon dioxide levels to estimate this deglaciation threshold. The model simulates several phenomena that are expected to be significant in a 'snowball Earth' scenario, but which have not been considered in previous studies with less sophisticated models, such as a reduction of vertical temperature gradients in winter, a reduction in summer tropopause height, the effect of snow cover and a reduction in cloud greenhouse effects. In my simulations, the system remains far short of deglaciation even at atmospheric carbon dioxide concentrations of 550 times the present levels (0.2 bar of CO2). I find that at much higher carbon dioxide levels, deglaciation is unlikely unless unknown feedback cycles that are not captured in the model come into effect.

The hydrologic cycle in deep-time climate problems.

Hydrology refers to the whole panoply of effects the water molecule has on climate and on the land surface during its journey there and back again between ocean and atmosphere. On its way, it is cycled through vapour, cloud water, snow, sea ice and glacier ice, as well as acting as a catalyst for silicate-carbonate weathering reactions governing atmospheric carbon dioxide. Because carbon dioxide affects the hydrologic cycle through temperature, climate is a pas des deux between carbon dioxide and water, with important guest appearances by surface ice cover.

Decay of passive scalars under the action of single scale smooth velocity fields in bounded two-dimensional domains: from non-self-similar probability distribution functions to self-similar eigenmodes.

We examine the decay of passive scalars with small, but nonzero, diffusivity in bounded two-dimensional (2D) domains. The velocity fields responsible for advection are smooth (i.e., they have bounded gradients) and of a single large scale. Moreover, the scale of the velocity field is taken to be similar to the size of the entire domain. The importance of the initial scale of variation of the scalar field with respect to that of the velocity field is strongly emphasized. If these scales are comparable and the velocity field is time periodic, we see the formation of a periodic scalar eigenmode. The eigenmode is numerically realized by means of a deterministic 2D map on a lattice. Analytical justification for the eigenmode is available from theorems in the dynamo literature. Weakening the notion of an eigenmode to mean statistical stationarity, we provide numerical evidence that the eigenmode solution also holds for aperiodic flows (represented by random maps). Turning to the evolution of an initially small scale scalar field, we demonstrate the transition from an evolving (i.e., non-self-similar) probability distribution function (pdf) to a stationary (self-similar) pdf as the scale of variation of the scalar field progresses from being small to being comparable to that of the velocity field (and of the domain). Furthermore, the non-self-similar regime itself consists of two stages. Both stages are examined and the coupling between diffusion and the distribution of the finite time Lyapunov exponents is shown to be responsible for the pdf evolution.

Surface quasigeostrophic turbulence: The study of an active scalar.

We study the statistical and geometrical properties of the potential temperature (PT) field in the surface quasigeostrophic (SQG) system of equations. In addition to extracting information in a global sense via tools such as the power spectrum, the g-beta spectrum, and the structure functions we explore the local nature of the PT field by means of the wavelet transform method. The primary indication is that an initially smooth PT field becomes rough (within specified scales), though in a qualitatively sparse fashion. Similarly, initially one-dimensional iso-PT contours (i.e., PT level sets) are seen to acquire a fractal nature. Moreover, the dimensions of the iso-PT contours satisfy existing analytical bounds. The expectation that the roughness will manifest itself in the singular nature of the gradient fields is confirmed via the multifractal nature of the dissipation field. Following earlier work on the subject, the singular and oscillatory nature of the gradient field is investigated by examining the scaling of a probability measure and a sign singular measure, respectively. A physically motivated derivation of the relations between the variety of scaling exponents is presented, the aim being to bring out some of the underlying assumptions which seem to have gone unnoticed in previous presentations. Apart from concentrating on specific properties of the SQG system, a broader theme of the paper is a comparison of the diagnostic inertial range properties of the SQG system with both the two- and three-dimensional Euler equations. (c) 2002 American Institute of Physics.

Lattice models of advection-diffusion.

We present a synthesis of theoretical results concerning the probability distribution of the concentration of a passive tracer subject to both diffusion and to advection by a spatially smooth time-dependent flow. The freely decaying case is contrasted with the equilibrium case. A computationally efficient model of advection-diffusion on a lattice is introduced, and used to test and probe the limits of the theoretical ideas. It is shown that the probability distribution for the freely decaying case has fat tails, which have slower than exponential decay. The additively forced case has a Gaussian core and exponential tails, in full conformance with prior theoretical expectations. An analysis of the magnitude and implications of temporal fluctuations of the conditional diffusion and dissipation is presented, showing the importance of these fluctuations in governing the shape of the tails. Some results concerning the probability distribution of dissipation, and concerning the spatial scaling properties of concentration fluctuation, are also presented. Though the lattice model is applied only to smooth flow in the present work, it is readily applicable to problems involving rough flow, and to chemically reacting tracers. (c) 2000 American Institute of Physics.