C. 1D radiative-convective model up previous next
C.b. Finite difference equations of the 1D model, simulation setup and results

The space differencing is evaluated by the second order centered scheme. The time integration is performed by the forward scheme for the radiative heating term and the Crank-Nicolson scheme for the diffusion term. Detail description of finite difference form of the radiative transfer equation and the vertical diffusion term are same as those shown in appendix B.d.i, appendix B.d.ii and appendix B.e

The arrangement of vertical grid and its interval are same as those of the 2D model (section 2.b). The time step is 30 seconds. The solar flux at the top of model atmosphere diurnally changes under the condition of Ls = 100° at 20°N. The initial atmospheric and ground surface temperature is constant (220 K and 210 K respectively). The numerical integration time is 30 days.

The vertical profiles of the atmospheric and the ground surface temperature at LT = 6:00 are shown in Figure C1 and Figure C2 . These profiles are used for the initial condition of the 2D model simulation (section 2.b). Even if the initial atmospheric temperature is 160 K, the result is similar to that obtained in the case that the initial atmospheric temperature is 220 K (Figure C3, Figure C4).

The diurnal changes of the atmospheric and the ground surface temperature at 30 days are shown in Figure C5 and Figure C6, respectively.

Figure C1: Vertical profile of the atmospheric temperature at LT = 6:00 on 30 days calculated by the 1D model. Figure C2: Vertical profile of the ground surface temperature at LT = 6:00 on 30th day calculated by the 1D model. The unit of vertical axis is the diurnal skin depth.

Figure C3: Time development of the vertical profile of atmospheric temperature at LT = 6:00 every 10 days calculated by the 1D model. with 220 K nitial atmospheric temperature. Number on each line indicates calculation day. Figure C4: Time development of the vertical profile of atmospheric temperature at LT = 6:00 every 10 days calculated by the 1D model with 220 K nitial atmospheric temperature. Number on each line indicates calculation day.

Figure C5: Time development of the vertical profile of atmospheric temperature every 4 hours on 30th day calculated by the 1D model. Number on each line indicates local time.

Figure C6: (left panel) Time development of vertical profile of the ground temperature from LT = 8:00 to 4:00 every 4 hours on 30th day. Number on each line indcats local time. The unit of vertical axis is the diurnal skin depth. (right panel) Time evolution of horizontal mean surface temperature (red line) and atmospheric temperature at about 1.5 m height (blue line).


A numerical simulation of thermal convection in the Martian lower atmosphere.
Odaka, Nakajima, Ishiwatari, Hayashi,   Nagare Multimedia 2001
up previous next