This article discussesl a few of the problems that arise in geophysical fluid dynamics and climate that are associated with the presence of moisture in the air, its condensation and release of latent heat. Our main focus is Earth's atmosphere but we also discuss how these problems might manifest themselves on other planetary bodies, with particular attention to Titan where methane takes on the role of water. GFD has traditionally been concerned with understanding the very basic problems that lie at the foundation of dynamical meteorology and ocean\-ography. Conventionally, and a little ironically, the subject mainly considers `dry' fluids, meaning it does not concern itself overly much with phase changes. The subject is often regarded as dry in another way, because it does not consider problems perceived as relevant to the real world, such as clouds or rainfall, which have typically been the province of complicated numerical models. Those models often rely on parameterizations of unresolved processes, parameterizations that may work very well but that often have a semi-empirical basis. The consequent dichotomy between the foundations and the applications prevents progress being made that has both a secure basis in scientific understanding and a relevance to the Earth's climate, especially where moisture is concerned. The dichotomy also inhibits progress in understanding the climate of other planets, where observations are insufficient to tune the parameterizations that weather and climate models for Earth rely upon, and a more fundamental approach is called for. Here we discuss four diverse examples of the problems with moisture: the determination of relative humidity and cloudiness; the transport of water vapor and its possible change under global warming; the moist shallow water equations and the Madden-Julian Oscillation; and the hydrology cycle on other planetary bodies.

The climate and circulation of a terrestrial planet are governed by, among other things, the distance to its host star, its size, rotation rate, obliquity, atmospheric composition and gravity. Here we explore the effects of the last of these, the Newtonian gravitational acceleration, on its atmosphere and climate. We first demonstrate that if the atmosphere obeys the hydrostatic primitive equations, which are a very good approximation for most terrestrial atmospheres, and if the radiative forcing is unaltered, changes in gravity have no effect at all on the circulation except for a vertical rescaling. That is to say, the effects of gravity may be completely scaled away and the circulation is unaltered. However, if the atmosphere contains a dilute condensible that is radiatively active, such as water or methane, then an increase in gravity will generally lead to a cooling of the planet because the total path length of the condensible will be reduced as gravity increases, leading to a reduction in the greenhouse effect. Furthermore, the specific humidity will decrease, leading to changes in the moist adiabatic lapse rate, in the equator-to-pole heat transport, and in the surface energy balance because of changes in the sensible and latent fluxes. These effects are all demonstrated both by theoretical arguments and by numerical simulations with moist and dry general circulation models.

Using an idealised general circulation model, we investigate the atmospheric circulation of Earth-like terrestrial planets in a variety of orbital configurations. We relax the common assumption of the planet being tidally-locked, and look at the role atmospheric dynamics can have in the observed thermal phase curve when the substellar point is non-stationary. In slowly rotating planets, a moving forcing can induce strong jets in the upper troposphere, both prograde and retrograde, sensitive to the speed and direction of the diurnal forcing. We find that, consistent with previous shallow water model experiments, the thermal phase curve offset is sensitive to the velocity of the substellar point moving across the surface of the planet. For a planet with a known orbital period, the results show that the observed hotspot on the planet could be either east or west of the substellar point, depending on whether the planet is tidally-locked or not.

The zonal-mean atmospheric flow of an idealized terrestrial planet is analyzed using both numerical simulations and zonally symmetric theories, focusing largely on the limit of low planetary rotation rate. Two versions of a zonally symmetric theory are considered, the standard Held-Hou model, which features a discontinuous zonal wind at the edge of the Hadley cell, and a variant with continuous zonal wind but discontinuous temperature. The two models have different scalings for the boundary latitude and zonal wind. Numerical simulations are found to have smoother temperature profiles than either model, with no temperature or velocity discontinuities even in zonally symmetric simulations. Continuity is achieved because of the presence of an overturning circulation poleward of the point of maximum zonal wind, which allows the zonal velocity profile to be smoother than the original theory without the temperature discontinuities of the variant theory. Zonally symmetric simulations generally fall between the two sets of theoretical scalings, and have a faster polar zonal flow than either. Three-dimensional simulations that allow for eddy motion fall closer to the scalings of the variant model. At very low rotation rates the maximum zonal wind falls with falling planetary rotation rate, even in the three-dimensional simulations, and collapses completely at zero rotation. Nevertheless, the low-rotation limit of the overturning circulation is strong enough to drive the temperature profile close to a state of nearly constant potential temperature.

Using a shallow water model with time-dependent forcing we show that the peak of an exoplanet thermal phase curve is, in general, offset from secondary eclipse when the planet is rotating. That is, the planetary hot-spot is offset from the point of maximal heating (the substellar point) and may lead or lag the forcing; the extent and sign of the offset is a function of both the rotation rate and orbital period of the planet. We also find that the system reaches a steady-state in the reference frame of the moving forcing. The model is an extension of the well studied Matsuno-Gill model into a full spherical geometry and with a planetary-scale translating forcing representing the insolation received on an exoplanet from a host star. The speed of the gravity waves in the model is shown to be a key metric in evaluating the phase curve offset. If the velocity of the substellar point (relative to the planet's surface) exceeds that of the gravity waves then the hotspot will lag the substellar point, as might be expected by consideration of forced gravity wave dynamics. However, when the substellar point is moving slower than the internal wavespeed of the system the hottest point can lead the passage of the forcing. We provide an interpretation of this result by consideration of the Rossby and Kelvin wave dynamics as well as, in the very slowly rotating case, a one-dimensional model that yields an analytic solution. Finally, we consider the inverse problem of constraining planetary rotation rate from an observed phase curve.

We show that by changing a single non-dimensional number, the thermal Rossby number, global atmospheric simulations with only axisymmetric forcing pass from an Earth-like atmosphere to a superrotating atmosphere that more resembles the atmospheres of Venus or Titan. The transition to superrotation occurs under conditions in which equatorward-propagating Rossby waves generated by baroclinic instability at intermediate and high latitudes are suppressed, which will occur when the deformation radius exceeds the planetary radius. At large thermal Rossby numbers following an initial, nearly axisymmetric phase, a global baroclinic wave of zonal wavenumber one generated by mixed barotropic-baroclinic instability dominates the eddy flux of zonal momentum. The global wave converges eastward zonal momentum to the equator and deposits westward momentum at intermediate latitudes during spinup and before superrotation emerges, and the baroclinic instability ceases once superrotation is established. A global barotropic mode of zonal wavenumber one generated by a mix of high- and low-latitude barotropic instability is responsible for maintaining superrotation in the statistically steady state. At intermediate thermal Rossby numbers, momentum flux by the global baroclinic mode is subdominant relative to smaller baroclinic modes, and thus strong superrotation does not develop.