Planetary Climate Dynamics

Abstract:

A dust transport scheme has been developed for a general circulation model of the Martian atmosphere. This enables radiatively active dust transport, with the atmospheric state responding to changes in the dust distribution via atmospheric heating, as well as dust transport being determined by atmospheric conditions. The scheme includes dust lifting, advection by model winds, atmospheric mixing, and gravitational sedimentation. Parameterizations of lifting initiated by (1) near-surface wind stress and (2) convective vortices known as dust devils are considered. Two parameterizations are defined for each mechanism and are first investigated offline using data previously output from the non-dust- transporting model. The threshold-insensitive parameterizations predict some lifting over most regions, varying smoothly in space and time. The threshold-sensitive parameterizations predict lifting only during extreme atmospheric conditions (such as exceptionally strong winds), so lifting is rarer and more confined to specific regions and times. Wind stress lifting is predicted to peak during southern summer, largely between latitudes 15° and 35°S, with maxima also in regions of strong slope winds or thermal contrast flows. These areas are consistent with observed storm onset regions and dark streak surface features. Dust devil lifting is also predicted to peak during southern summer, with a moderate peak during northern summer. The greatest dust devil lifting occurs in early afternoon, particularly in the Noachis, Arcadia/Amazonis, Sirenum, and Thaumasia regions. Radiatively active dust transport experiments reveal strong positive feedbacks on lifting by near-surface wind stress and negative feedbacks on lifting by dust devils.

Abstract:

Multiannual dust transport simulations have been performed using a Mars general circulation model containing a dust transport scheme which responds to changes in the atmospheric state. If the dust transport is "radiatively active," the atmospheric state also responds to changes in the dust distribution. This paper examines the suspended dust distribution obtained using different lifting parameterizations, including an analysis of dust storms produced spontaneously during these simulations. The lifting mechanisms selected are lifting by (1) near-surface wind stress and (2) convective vortices known as dust devils. Each mechanism is separated into two types of parameterization: threshold-sensitive and -insensitive. The latter produce largely unrealistic annual dust cycles and storms, and no significant interannual variability. The threshold-sensitive parameterizations produce more realistic annual and interannual behavior, as well as storms with similarities to observed events, thus providing insight into how real Martian dust storms may develop. Simulations for which dust devil lifting dominates are too dusty during northern summer. This suggests either that a removal mechanism (such as dust scavenging by water ice) reduces opacities at this time or that dust devils are not the primary mechanism for storm production. Simulations for which near-surface wind stress lifting dominates produce the observed low opacities during northern spring/ summer, yet appear unable to produce realistic global storms without storm decay being prevented by the occurrence of large-scale positive feedbacks on further lifting. Simulated dust levels are generally linked closely to the seasonal state of the atmosphere, and no simulation produces the observed amount of interannual variability.

Abstract:

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.