It has been recognized that
radiative heating caused by atmospheric dust
is one of the most important factors in determining
the thermal and circulation structure
of the Martian atmosphere.
The vertical profiles of temperature observed by spacecrafts
frequently give stable lapse rate compared to dry adiabat
(e.g., Lindal et
al., 1979).
Those stable profiles are considered to be caused by radiative
heating associated with dust as is demonstrated by
one-dimensional (1D) radiative convective models
(Gierasch and Goody,
1972; Pollack et
al., 1979).
Dust is also considered to intensify the magnitude of
the Martian large scale atmospheric circulation,
as is suggested by the comparison between
the simulation results
under dusty and clear sky (i.e., dust-free) conditions
obtained with general circulation models (GCMs)
(e.g.,
Pollack et al., 1990).
However, GCMs have not succeeded in
the prognostic calculation of the amount of dust,
the important element of the Martian atmosphere.
The global dust storm, which is one of the most striking
phenomena in the Martian atmosphere
(e.g., Briggs
et al., 1979),
has not been self consistently
simulated in the current GCMs.
Starting from a dusty atmosphere as the initial condition,
the GCMs produce intense large scale wind
to inject dust into the atmosphere to maintain the
dust amount in the atmosphere.
However, starting from a dust-free or
a small amount of dust atmosphere,
the intensity of the large scale winds calculated in the GCMs
is so weak that the model surface stress is not large
enough to raise dust from the surface.
Current GCMs can not describe
a spontaneous transition from the dust-free Mars to
the dusty Mars
(Joshi et al., 1997;
Wilson and Hamilton, 1996).
It is argued by
Wilson and Hamilton (1996)
that the amount of surface stress
which is necessary for dust injection into the atmosphere
could be supplemented by small-scale wind fluctuations
which are not resolved in GCMs.
However, they presented neither the nature nor the origin of
small-scale wind motions.
As a candidate of the small-scale motions which are not
represented in GCM, one can think of vertical convection
driven by radiative forcing and sensible heat from the
surface.
Actually, it has been pointed out that
the data obtained at the site of Viking Lander 1
demonstrates the existence of vertical convection
(Hess et al., 1977;
Ryan and Lucich, 1983).
According to the studies by the use of 1D radiative convective
models,
the depth of convection layer is estimated to be
about from 9 to 10 km under the dust-free condition
( Flasar and Goody, 1976;
Pollack et al., 1979),
and
is about from 3 to 4 km under the dusty condition
when the visible optical depth of dust is from 0.3 to 0.5
(Savijärvi, 1991b;
Haberle et al., 1993).
However, since the fluid dynamical aspects of vertical
convection in the Martian atmosphere has not been
intensely examined so far,
the fundamental characteristics such as cell patterns and
wind intensity associated with vertical convection
have not been well revealed.
Vertical convection in the Martian atmosphere can be regarded
as dry convection in general.
Since the amount of water vapor is quite small,
condensation heating of water vapor can be neglected
compared to radiative heating in the Martian atmosphere
(e.g.,
Zurek et al., 1992).
Condensation of CO2,
the major component of the Martian atmosphere,
does not have to be considered except in the poler
region.
Dry convection occurs also in the terrestrial atmosphere
in the planetary boundary layer near the surface.
However,
dry convection in the Martian atmosphere
covers the whole height of Martian troposphere.
The characteristics of such "deep" dry convection
has been little examined so far.
In this study, we construct a numerical model which can
explicitly represent convective fluid motion, and
investigate possible characteristics of the vertical
convection in the Martian atmosphere.
We perform the following two cases of numerical simulations;
-
Simulation of vertical convection without dust
(the dust-free case):
Focus is placed on the patterns of circulation and intensity
of wind, and the amount of surface friction due to the
wind.
By using the calculated values of surface friction,
we discuss the possibility of dust injection from the
surface, which has not been realized in GCMs under the
dust free condition.
-
Simulation of vertical convection without dust
(the dusty case):
Assuming that the convective wind injects dust into
the atmosphere,
we investigate the characteristics of dust mixing by thermal
convection and the effects of radiative heating due to
dust on the circulation structure of vertical convection.
The numerical model utilized here is spatially two-dimensional
(2D).
The advantage of restricting the computational domain
to two-dimensional space is that
we can employ both a large computational domain and
a sufficiently high spatial resolution that
resolves convective plumes explicitly.
Moreover,
by using a 2D model,
it is expected to be easier to recognize some
characteristic features of convection, if any, than by
using a three-dimensional (3D) model.
|