Effects of thermal tides on the Venus atmospheric superrotation

Similar documents
Dynamical effects of the thermal tides in the Venus atmosphere

Recent Development of CCSR/NIES AGCM for Venus Atmospheric Sciences

Simulations of superrotation using a GCM for Venus middle atmosphere

The Fully Developed Superrotation Simulated by a General Circulation Model of a Venus-like Atmosphere

no eddies eddies Figure 3. Simulated surface winds. Surface winds no eddies u, v m/s φ0 =12 φ0 =0

A Further Study of Dynamics of the Four-Day Circulation. By Yoshihisa Matsuda

Anomalous solar heating dependence of Venus s cloud-level convection

Traveling planetary-scale Rossby waves in the winter stratosphere: The role of tropospheric baroclinic instability

CHAPTER 4. THE HADLEY CIRCULATION 59 smaller than that in midlatitudes. This is illustrated in Fig. 4.2 which shows the departures from zonal symmetry

Superrotation in a Venus general circulation model

Lecture #1 Tidal Models. Charles McLandress (Banff Summer School 7-13 May 2005)

Wave-driven equatorial annual oscillation induced and modulated by the solar cycle

Icarus. Decadal variations in a Venus general circulation model

Traveling planetary-scale Rossby waves in the winter stratosphere: The role of tropospheric baroclinic instability

Super rotation in a Venus General Circulation Model

A mechanistic model study of quasi-stationary wave reflection. D.A. Ortland T.J. Dunkerton NorthWest Research Associates Bellevue WA

High initial time sensitivity of medium range forecasting observed for a stratospheric sudden warming

PUBLICATIONS. Journal of Geophysical Research: Planets

What kind of stratospheric sudden warming propagates to the troposphere?

arxiv: v1 [astro-ph.ep] 21 Sep 2016

Equatorial Superrotation on Tidally Locked Exoplanets

Lecture 10a: The Hadley Cell

Studies of Venus using a General Circulation Model

What We Know Today About the Venus Middle Atmosphere

Responses of mesosphere and lower thermosphere temperatures to gravity wave forcing during stratospheric sudden warming

Eliassen-Palm Theory

Numerical simulation of the equatorial wind jet in the thermosphere

7 The General Circulation

Physical Processes in Acoustic Wave Heating of the Thermosphere

Dynamical and Thermal Effects of Gravity Waves in the Terrestrial Thermosphere-Ionosphere

2. Meridional atmospheric structure; heat and water transport. Recall that the most primitive equilibrium climate model can be written

Solar Semidiurnal Tide in the Dusty Atmosphere of Mars

Tidal Coupling in the Earth s Atmosphere. Maura Hagan NCAR High Altitude Observatory

EART164: PLANETARY ATMOSPHERES

3. Midlatitude Storm Tracks and the North Atlantic Oscillation

Modeling of Jupiter s stratosphere: new radiation code and impacts on the dynamics

Introduction to Atmospheric Circulation

Analysis of Ultra-fast Kelvin Waves Simulated by the Kyushu University GCM

Day-to-day variations of migrating semidiurnal tide in the mesosphere and thermosphere

A Numerical Study on Appearance of the Runaway Greenhouse State of a Three-Dimensional Gray Atmosphere

Meridional Transport in the Stratosphere of Jupiter

Dynamics of Giant Planet Atmospheres. Tapio Schneider (with Junjun Liu) California Institute of Technology

Venus Upper Atmosphere Circulation Models and Coupling from Above and Below

Alan Z. Liu Embry Riddle Aeronautical University - Daytona Beach, Chester S. Gardner

Stratosphere Troposphere Coupling in a Relatively Simple AGCM: Impact of the Seasonal Cycle

The Transition to Superrotation in Terrestrial Atmospheres

The stratospheric response to extratropical torques and its relationship with the annular mode

Imperial College London

First detection of wave interactions in the middle atmosphere of Mars

Model description of AGCM5 of GFD-Dennou-Club edition. SWAMP project, GFD-Dennou-Club

Scales of Linear Baroclinic Instability and Macroturbulence in Dry Atmospheres

7 The Quasi-Biennial Oscillation (QBO)

Radiative-Convective Models. The Hydrological Cycle Hadley Circulation. Manabe and Strickler (1964) Course Notes chapter 5.1

A new pathway for communicating the 11-year solar cycle signal to the QBO

Introduction to Climate ~ Part I ~

The Tropospheric Jet Response to Prescribed Zonal Forcing in an Idealized Atmospheric Model

An Examination of Anomalously Low Column Ozone in the Southern Hemisphere Midlatitudes During 1997

Lecture #2 Planetary Wave Models. Charles McLandress (Banff Summer School 7-13 May 2005)

Introduction to Isentropic Coordinates: a new view of mean meridional & eddy circulations. Cristiana Stan

The effect of varying forcing on the transport of heat by transient eddies.

Hadley Circulation as a Modulator of the Extratropical Climate

Is Antarctic climate most sensitive to ozone depletion in the middle or lower stratosphere?

A Mechanism for the Effect of Tropospheric Jet Structure on the Annular Mode Like Response to Stratospheric Forcing

1 Climatological balances of heat, mass, and angular momentum (and the role of eddies)

Downward migration of extratropical zonal wind anomalies

P4.2 THE THREE DIMENSIONAL STRUCTURE AND TIME EVOLUTION OF THE DECADAL VARIABILITY REVEALED IN ECMWF REANALYSES

Eliassen-Palm Cross Sections Edmon et al. (1980)

JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 102, NO. E3, PAGES , MARCH 25, balance of the Venusian midlatitude

The Matsuno Gill model and equatorial superrotation

Influence of Doubled CO 2 on Ozone via Changes in the Brewer Dobson Circulation

Zonal Momentum Balance in the Tropical Atmospheric Circulation during the Global Monsoon Mature Months

Effect of zonal asymmetries in stratospheric ozone on simulated Southern Hemisphere climate trends

Four ways of inferring the MMC. 1. direct measurement of [v] 2. vorticity balance. 3. total energy balance

Dynamical formation of an extra-tropical tropopause inversion layer in a relatively simple general circulation model

The Effect of a Hadley Circulation on the Propagation and Reflection of Planetary Waves in a Simple One-Layer Model

A non-hydrostatic and compressible 2-D model simulation of Internal Gravity Waves generated by convection

Dynamics in Planetary Atmospheric Physics: Comparative Studies of Equatorial Superrotation for Venus, Titan, and

General Circulation. Nili Harnik DEES, Lamont-Doherty Earth Observatory

The Impact of the State of the Troposphere on the Response to Stratospheric Heating in a Simplified GCM

Dynamics and Kinematics

Upper atmosphere response to stratosphere sudden warming: Local time and height dependence simulated by GAIA model

Geophysics Fluid Dynamics (ESS228)

Parity of solar global magnetic field determined by turbulent diffusivity

Goal: Understand the dynamics and thermodynamics of the Hadley circulation

Sensitivity of Precipitation in Aqua-Planet Experiments with an AGCM

Characteristics of Wave Induced Oscillations in Mesospheric O2 Emission Intensity and Temperature

Influence of Sudden Stratosphere Warmings on the Ionosphere and Thermosphere

1. The vertical structure of the atmosphere. Temperature profile.

The role of synoptic eddies in the tropospheric response to stratospheric variability

Idealized numerical experiments on microscale eddies in the Venusian cloud layer

Extremely cold and persistent stratospheric Arctic vortex in the winter of

Modeling the Downward Influence of Stratospheric Final Warming events

The general circulation of the atmosphere

A ground-to-exosphere Martian general circulation model: 2. Atmosphere during solstice conditions Thermospheric polar warming

By STEVEN B. FELDSTEINI and WALTER A. ROBINSON* University of Colorado, USA 2University of Illinois at Urbana-Champaign, USA. (Received 27 July 1993)

Forced Annular Mode Patterns in a Simple Atmospheric General Circulation Model

A new perspective on gravity waves in the Martian atmosphere: Sources and features

Lecture #3: Gravity Waves in GCMs. Charles McLandress (Banff Summer School 7-13 May 2005)

A Truncated Model for Finite Amplitude Baroclinic Waves in a Channel

Mechanisms for influence of the stratosphere on the troposphere

Transcription:

JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jd007901, 2007 Effects of thermal tides on the Venus atmospheric superrotation M. Takagi 1 and Y. Matsuda 2 Received 10 August 2006; revised 9 December 2006; accepted 5 January 2007; published 9 May 2007. [1] A nonlinear dynamical model on the sphere has been numerically integrated to investigate a generation mechanism of the Venus atmospheric superrotation by the thermal tides. By using the solar heating exciting the diurnal and semidiurnal tides, the atmospheric superrotation extending from the ground to 80 km is generated. The vertical distributions of the mean zonal flow obtained in our experiments are similar to the observations. Velocity of the mean zonal wind on the equator reaches about 60 100 m s 1 near the cloud top level. A linear theory suggests that the atmospheric superrotation obtained in the present study is generated and maintained by the momentum transport associated with the thermal tides. Namely, the downward transport of zonal momentum that is associated with the downward propagating semidiurnal tide excited in the cloud layer induces the mean zonal flow opposite to the Venus rotation in the lowest layer adjacent to the ground. Surface friction acting on this counter flow provides the atmosphere with the net angular momentum from the solid part of Venus. It is examined how the atmospheric superrotation depends on vertical eddy viscosity and Newtonian cooling. The result shows that magnitude of the atmospheric superrotation is not so sensitive to vertical eddy viscosity but is strongly influenced by Newtonian cooling. Citation: Takagi, M., and Y. Matsuda (2007), Effects of thermal tides on the Venus atmospheric superrotation, J. Geophys. Res., 112,, doi:10.1029/2006jd007901. 1. Introduction [2] The atmospheric superrotation of Venus is one of the most remarkable phenomena in the planetary meteorology. The Venus rotation is so slow that the atmospheric motion on Venus was expected to be approximately axisymmetric about a line between the subsolar and antisolar points (for stability of this axisymmetric circulation, see the works of Takagi and Matsuda [1999, 2000]). However, a number of observations show that the fast zonal flow (superrotation) is predominant in the Venus atmosphere [Schubert, 1983]. It is well known that the almost entire atmosphere of Venus rotates much faster than the Venus itself. The mean zonal wind increases monotonically from the ground and reaches about 100 m s 1 near the cloud top level. The corresponding angular velocity is about 60 times larger than that of the solid part of Venus. Although many studies have been made so far, the mechanism generating the superrotation of Venus atmosphere is not still clearly understood. [3] Fels and Lindzen [1974] and Plumb [1975] showed a possibility that the flow retrograde to the solar motion (i.e., atmospheric superrotation) can be generated by momentum transport associated with vertical propagation of the thermal 1 Department of Earth and Planetary Science, University of Tokyo, Tokyo, Japan. 2 Department of Astronomy and Earth Science, Tokyo Gakugei University, Tokyo, Japan. Copyright 2007 by the American Geophysical Union. 0148-0227/07/2006JD007901 tides excited in the heating layer. Covey et al. [1986] also argued that the diurnal tide behaves like a thermally direct subsolar-to-antisolar circulation in a diffusive lower Venus atmosphere and the superrotation may be generated there. In the Venus atmosphere, the thick cloud covers the whole planet at 45 70 km altitudes and the solar heating is mainly absorbed there. It is likely that the thermal tides strongly excited in the cloud layer have important roles on the atmospheric circulation of Venus. The dynamical effects of the thermal tides on the Venus atmospheric superrotation have been examined further by Baker and Leovy [1987], Hou et al. [1990], and Newman and Leovy [1992]. By using a numerical model that spans 30 110 km altitudes, Newman and Leovy showed that the atmospheric superrotation near the cloud top can be maintained by the momentum transport of thermal tides. However, it should be noted that increase of the net angular momentum that enables the superrotation to form at the cloud levels is produced by Rayleigh friction acting near the upper boundary in their model. Moreover, since the lower boundary is located at 30 km, the downward propagation of the thermal tides (i.e., downward momentum transport) cannot be correctly represented in their model. A vertically uniform background mean zonal flow of solid body rotation with equatorial speed of 75 m s 1 is also assumed in their study. [4] By using a Venus-like atmospheric general circulation model (GCM) with the zonally uniform solar heating, Yamamoto and Takahashi [2003a, 2003b] reproduced atmospheric superrotation faster than 100 m s 1 near 60 km altitude. In their results, a strong meridional circulation with 1of8

a large single cell extending from the ground to about 80 km altitude appears and seems to play an important role on the angular momentum transport. They attempt to interpret the superrotation appearing in their model by the generation mechanism proposed by Gierasch [1975] (see also the works of Matsuda [1980, 1982]). It should be noted, however, that the solar heating used in their model is too strong at the lower levels. From observational data reported by Tomasko et al. [1980], it can be estimated that the averaged solar heating rate is about 3.0 10 3 K day 1 in a layer of 0 10 km altitudes (assuming that the solar flux absorbed at the ground is redistributed in this layer), whereas that used by Yamamoto and Takahashi is 0.52 K day 1 there. Hence it is possible that the meridional circulation induced in their model is overemphasized by this heating. Yamamoto and Takahashi [2004] extended these works by adopting a solar heating with diurnal variation. It is shown that the atmospheric superrotation is generated by a similar mechanism except that the angular momentum advected poleward by the meridional circulation is returned back to the lower latitudes by the thermal tides. The strong meridional circulation reproduced by the strong heating in the lower atmosphere is comparable with that of the preceding studies. Their studies are very interesting since the superrotation comparable with the observations is reproduced in the GCM by the strong meridional circulation. However, the reproduction of superrotation by the overestimated heating below the cloud layer does not necessarily demonstrate that the generating mechanism specified in their model is working in the real Venus atmosphere. Lee et al. [2005] also succeeded in reproducing a significant global atmospheric superrotation and discussed the atmospheric waves and the polar vortex which appeared in their result. It seems, however, that the generating mechanism of superrotation in their model is not fully clarified. [5] Recently, Takagi and Matsuda [2005, 2006] examined the thermal tides in the Venus atmosphere in detail by using a linear model based on the observed zonal flows. It is clearly shown that the thermal (semidiurnal) tide excited in the cloud layer propagates downward to the ground and accelerates the atmosphere in the direction opposite to the Venus rotation at altitudes of 0 10 km. They argued from these results that the surface friction acts on this counter flow, and as a result the net momentum is supplied from the solid part of Venus to maintain the atmospheric superrotation. It should be noticed here that no mean zonal flow retrograde to the Venus rotation has been observed by any of several missions. However, this fact is not necessarily contradictory to the above hypothesis (for details, see section 3.1). [6] In the present study, by extending these linear studies, we investigate whether the atmospheric superrotation of Venus can be generated by the nonlinear effect of the thermal tides, integrating a nonlinear dynamical model in spherical geometry. 2. Model [7] In order to investigate the generation mechanism of the Venus atmospheric superrotation by examining the interaction among the thermal tides, the mean zonal flows and the surface friction, a full nonlinear dynamical model on the sphere is constructed. The basic equations are primitive ones in the spherical and s coordinates [Hoskins and Simmons, 1975]. Time integration is performed using the semi-implicit method. Model performance of the dynamical core is checked by comparison with the results of Held and Suarez [1994]. Mean zonal jets and temperature distributions reproduced by the present model are very similar to those of Held and Suarez. The model atmosphere extends from the ground to 120 km, and 60 levels are taken at a regular spacing of 2 km. The horizontal resolution is T10 (triangular truncation at wave number 10). This horizontal resolution is so coarse that numerical convergence and/or truncation errors may well be questioned. It has been confirmed for a few cases, however, that the present model with T21 resolution gives similar results to those described in the following sections. The values of physical parameters for this Venus atmospheric model are as follows: the radius of planet is 6052 km, the period of planetary rotation 243 Earth days, the solar day 117 Earth days, the gravity 8.9 m s 2,the standard surface pressure 9.2 10 6 Pa, and the gas constant 191.4 J kg 1 K 1 for the CO 2 atmosphere. [8] In order to represent the radiative processes simply, the solar heating, Q, is prescribed in the present study. Its vertical structure is determined by the works of Tomasko et al. [1980] and Crisp [1986]. Vertical distribution of the specific heat at constant pressure, C p, is taken from the Venus international reference atmosphere (VIRA) data [Seiff et al., 1985]. Figure 1 shows the vertical profiles of C p and Q at the subsolar point. The solar flux absorbed at the ground is assumed to be redistributed equally into the lowest atmosphere with 10 km thickness just above the ground. The solar heating above 80 km is artificially attenuated in order to focus on the dynamical effects of the thermal tides excited below the cloud top level. The heating rates represented in K s 1 become extremely small below 50 km due to the dense lower atmosphere of Venus. It is assumed after the work of Crisp [1986] that the horizontal structure of the solar heating on the day side is proportional to cos 0.6 q S (cos 1.5 q S ) above (below) 71 km altitude where q S is the solar zenith angle. In the following calculations, the solar heating is decomposed into the mean zonal component, Q, and the deviation from it, Q 0 : Qðf; q; zþ ¼Qðq; zþþq 0 ðf; q; zþ; where f is longitude, q latitude, and z altitude. It is expected that the thermal tides and the mean meridional flow (Hadley circulation) are excited by Q 0 and Q, respectively. In the following, Q 0 is referred to as tidal component of the solar heating. [9] The radiative forcing in the infrared region is replaced by a linear Newtonian cooling scheme with coefficients based on the work of Crisp [1986]. The temperature field is relaxed to the prescribed horizontally uniform temperature field, T 0 (z), which is taken from the VIRA data. The thermal relaxation to a horizontally uniform temperature field is consistent with the exclusion of the mean zonal component of solar heating. The vertical profiles of T 0 and the relaxation time of Newtonian cooling are shown in Figure 2. [10] Coefficients of vertical eddy viscosity and heat diffusion are set to be 2.5 10 3,2.5 10 2,or2.5 10 1 m 2 s 1, ð1þ 2of8

Figure 1. Vertical profiles of the solar heating at the subsolar point, Q s (z) (solid line), and the specific heat at constant pressure, C p (z) (dotted line). and horizontal eddy viscosity is represented by the fourthorder hyperviscosity with relaxation time 1.0 Earth day for the maximum wave number component. In addition, we adopt molecular viscosity and heat diffusion that are dependent on temperature and density of the basic state field. Their effects are important above 90 km altitude. Rayleigh friction is not adopted in the present model except at the lowest level, where the surface friction acts on horizontal winds. The relaxation time of Rayleigh friction at the lowest level is assumed to be 1.0 Earth day. [11] The initial state for time integration is an atmosphere at rest for all experiments. It should be emphasized that the background atmospheric superrotation nor the meridional gradient of temperature is not assumed at all throughout the present study. 3. Results [12] In order to isolate specific effects of the thermal tides on the generation and maintenance of the Venus atmospheric superrotation, only the tidal component of solar heating, Q 0, is used for time integration in the present study. The model is integrated for 300 Earth years for all experiments. Parameter values chosen for five experiments conducted in the present study (cases of M1, M2, M3, W2, and S2) are summarized in Table 1. 3.1. Standard Case [13] Figure 3 shows time evolution of the mean zonal flow on the equator at 65, 45, and 5 km altitudes obtained for the case of M2 (see Table 1). In the present study, zonal velocity in the same direction as the planetary rotation is defined to be positive. Though the relaxation time of the vertical eddy viscosity is estimated to be larger than 1000 Earth years, it seems that the wind field is in a nearly equilibrium state. The zonal velocity at 65 km shows temporal variations with amplitude of a few to 10 m s 1. It is interesting to note that temporal variation has been also observed for the real atmospheric superrotation of Venus [Gierasch et al., 1997]. [14] Meridional-height distribution of the mean zonal flow and the mean temperature deviation from the prescribed equilibrium state obtained at 300 Earth years is illustrated in Figure 4. [15] The atmospheric superrotation has been generated in a wide range of the atmospheric layers below 80 km altitude. The mean zonal flow takes a maximum velocity of about 65 75 m s 1 near 63 km altitude on the equator. The second local maximum of the mean zonal flow exists just below 80 km. It is likely that this profile of the mean zonal flow with two maxima results from the prescribed vertical profile of solar heating with two peaks at the corresponding altitudes, as shown in Figure 1. The mean zonal flow increases with altitude almost linearly from the ground to 63 km and decreases sharply above 75 km. Above 80 km, the mean zonal flow remains very weak in all latitudes. It is also shown in Figure 4b that there is a reversal of the latitudinal temperature gradient near 60 km. In the light of the gradient wind balance, this reversed temperature gradient is well consistent with the (negative) vertical shear of mean zonal flow. Since the mean zonal component of solar heating, Q, is excluded from our experiments, the distribution of mean temperature deviation is not forced by meridional differential heating but induced by the formation of atmospheric superrotation. [16] These features of the mean zonal flow are similar to the observational results [Schubert, 1983], though the zonal velocity produced in this experiment is slower by 20 30 m s 1 than that observed in the cloud layer. Moreover, observed temperature fields indicate that midlatitude jets exist at cloud levels [Newman et al., 1984], which are not reproduced in the present result. This should be ascribed to the fact that the exclusion of the mean zonal component of the solar heating results in the absence of poleward advection of zonal momentum by the meridional circulation. [17] It is also shown in Figures 3 and 4a that the mean zonal flow retrograde to the planetary rotation (u < 0)is generated at the lowest levels in lower latitudes than 30. The time evolution of mean zonal flow at 5 km altitude Figure 2. Vertical profiles of the prescribed temperature field, T 0 (z) (solid line), and the relaxation time of Newtonian cooling, t N (z) (dotted line). 3of8

Table 1. Parameter Values Adopted in Five Numerical Experiments Case Vertical eddy viscosity (m 2 s 1 ) Factor of Newtonian cooling rates a M1 2.5 10 3 1 M2 2.5 10 2 1 M3 2.5 10 1 1 W2 2.5 10 2 1/3 S2 2.5 10 2 3 a Standard values of Newtonian cooling rate are taken from the work of Crisp [1986]. shown in Figure 3 indicates that this retrograde mean zonal flow is stronger in the early stage (about 0 100 Earth years) of the time integration. This distribution of the retrograde mean zonal flow seems very consistent with the numerical result and suggestion of Takagi and Matsuda [2006]. Figure 5 shows a vertical profile of vertical momentum flux, ru 0 w 0, averaged over the region equatorward of 30 degrees and 200 250 Earth years. It is indicated by this profile that the mean zonal flows tend to be accelerated at about 50 70 km altitudes and decelerated at about 0 10 and 70 80 km. It is also noted that this result is very consistent with linear calculations for the semidiurnal tide [Takagi, 2002; Takagi and Matsuda, 2006]. These results enable us to infer that in the present model the thermal tides excited in the cloud layer propagate to the ground and the momentum transport associated with the downward propagation of the semidiurnal tide induces the mean zonal flow retrograde to the planetary rotation in the lower latitudes at 0 10 km altitudes. Since Rayleigh friction is not adopted except at the lowest level in the present model, the net angular momentum supporting the atmospheric superrotation must be provided from the solid part by the surface friction acting on the mean zonal flow adjacent to the ground. As mentioned above, this retrograde mean zonal flow appears only in the early stage of the time integration and nearly vanishes after the superrotation has reached a (quasi) steady state. This result seems consistent with the fact that no mean zonal flow retrograde to the Venus rotation has been observed yet. [18] Figure 6 depicts the time evolution of the vertical profiles of mean zonal flow on the equator at 10, 50, and 300 Earth years. It is seen from these profiles that the atmospheric superrotation is generated near 60 km first and then extends above and below. The extension of the superrotation in the vertical direction results from the effect of vertical eddy and molecular viscosity. Before 10 Earth year, the superrotation is accompanied with the counter flows above and below it; this fact is consistent with the prediction by Fels and Lindzen [1974]. However, the counter flow just below the superrotation is very weak and instead another faster counter flow is found in the lowest layer of 0 10 km altitudes as shown in Figure 6b. It should be noted that the equatorial velocity of the mean zonal flow above 80 km appears to converge at 4 ms 1, which is the zonal phase velocity of the solar heating on the equator. This result means that the critical level is formed just above 80 km and the upward propagating thermal tides are absorbed there by thermal damping due to Newtonian cooling. (Note that Figure 3. Time evolution of mean zonal flow on the equator at 65, 45, and 5 km altitudes obtained for the case of M2. Figure 4. Meridional-height distribution of (a) mean zonal flow (m s 1 ) and (b) mean temperature deviation (K) at 300 Earth years obtained for the case of M2. Contour intervals are (a) 10 m s 1 and (b) 0.3 K, respectively. Contour lines for 0 K are omitted in Figure 4b. 4of8

Figure 5. A vertical profile of vertical momentum flux, ru 0 w 0, averaged over the region equatorward of 30 degrees and 200 250 Earth years for the M2 case. The unit is kg m 1 s 2. thermal damping due to Newtonian cooling is very weak below the cloud bottom.) As a result, the mean zonal flow is decelerated below the critical level and the vertical shear of the mean zonal flow between the cloud top and the critical level is maintained in this experiment. Since the strong solar heating above the cloud top levels, which is excluded from our experiments, would induce strong thermal tides there, this vertical profile of mean zonal flow with a critical level in this layer cannot appear in the real atmosphere of Venus. 3.2. Dependency on Vertical Eddy Viscosity [19] In order to examine how the atmospheric superrotation generated by the thermal tides depends on the vertical eddy viscosity, numerical experiments are conducted for the cases of small vertical eddy viscosity (M1) and large one (M3). See Table 1 for their values. [20] Figure 7 shows the vertical profiles of the mean zonal flow obtained for M1, M2, and M3. In all the cases, peaks of the mean zonal flow are located near 63 km. It is found that their maximum velocity is not so sensitive to the vertical eddy viscosity, though the velocity obtained for M3 is slightly smaller than those for M1 and M2. This result implies that the magnitude of the mean zonal flow is influenced by the mean meridional temperature gradient that is in the gradient wind balance with the zonal flow rather than the vertical eddy viscosity (see also section 3.3). It should be noted that this is not contradictory to the fact that the mean temperature deviation is induced by the formation of superrotation because the mean temperature deviation induced in this way is strongly influenced by the radiative process. [21] In the case of M1, the mean zonal flow retrograde to the Venus rotation remains in the lowest levels of 0 10 km, as shown in Figure 7. In other cases, such a retrograde flow vanishes in an early stage of the time integration due to the stronger vertical eddy viscosity. Figure 6. Vertical profiles of mean zonal flow on the equator at 10 (dotted line), 50 (dashed line), and 300 (solid line) Earth years obtained for the case of M2. Figure 6b enlarges those at 0 20 km altitudes. Figure 7. Vertical profiles of mean zonal flow on the equator obtained for M1 (solid line), M2 (dashed line), and M3 (dotted line). Note that the profiles are averaged over 250 300 Earth years. 5of8

Figure 8. As in Figure 7 but for W2 (solid line), M2 (dashed line), and S2 (dotted line). 3.3. Dependency on Newtonian Cooling [22] In this subsection, in order to examine how Newtonian cooling influences the generation of atmospheric superrotation, numerical experiments are conducted for the cases with coefficients of Newtonian cooling multiplied by factors of 1/3 (W2) and 3 (S2). The adopted value of vertical eddy viscosity is the same as M2. [23] Figure 8 shows the vertical profiles of mean zonal flow on the equator obtained for the cases of W2, M2, and S2. The mean zonal flow is more sensitive to Newtonian cooling than to the vertical eddy viscosity. Velocities of the mean zonal flow obtained in the three cases are considerably different in a layer of 40 80 km altitudes. The maximum velocity of mean zonal flow increases (decreases) to about 95 (60) m s 1 for weak (strong) Newtonian cooling. It is also found from time series of the mean zonal flow that the maximum velocity fluctuates between about 80 and 110 m s 1 for the W2 case (not shown). In the case of S2, the mean zonal flow decreases sharply with altitude above 63 km. This is because the thermal tides are damped there by strong Newtonian cooling. The corresponding meridional distributions of zonal mean temperature deviation at altitudes of 43 47 km obtained for W2, M2, and S2 are shown in Figure 9. It is understood that the larger (smaller) vertical shear of mean zonal flow is associated with the larger (smaller) meridional temperature gradient in the weak (strong) Newtonian cooling case. The magnitude of the atmospheric superrotation is strongly influenced by Newtonian cooling in the present model. 4. Discussion 4.1. Pumping Mechanism of Momentum [24] Though the solar heating used in the present study does not involve its zonal mean component, the present results shown in section 3 suggest that the atmospheric superrotation of Venus can be generated and maintained by the dynamical effects of the thermal tides. As suggested by Takagi and Matsuda [2006], the surface friction acting on the mean zonal flow in the lowest layer, which is created by the downward transport of the zonal momentum associated with the semidiurnal tide excited in the cloud layer, can provide the atmosphere with the net angular momentum required for the maintenance of the global superrotation. It might be speculated, however, that this pumping mechanism of the momentum is inconsistent with the work of Fels and Lindzen [1974]. They predicted that the critical levels were formed just above and below the heating layer and the thermal tides could not propagate across these critical levels. Since the angular momentum cannot be transported across these levels, the generation of the strong superrotation should be prevented. However, only the critical level above the cloud layer is formed in the present experiments. As shown by Takagi and Matsuda [2006], the thermal tides propagating downward from the cloud layer are not necessarily damped below the cloud layer, but reach the ground. This is because Newtonian cooling is much weaker in the lower atmosphere. [25] Figure 10 is a schematic illustration of acceleration mechanism of the mean zonal flow induced by the thermal tides in the work of Fels and Lindzen [1974] and the present one. Since the critical level below the cloud layer is not formed in the present experiments, the angular momentum associated with the downward propagating thermal tides is transported from the cloud layer to the lowest layer adjacent to the ground. This view is also supported by the vertical momentum flux shown in Figure 5. [26] It should be noted here that Hou and Farrel [1987] examined the effect of critical levels in the Venus atmosphere. They proposed that the superrotation may be induced by the critical level absorption of the upward propagating gravity waves excited in the lower atmosphere. This mechanism requires gravity waves with zonal phase velocities much faster than the mean zonal flow at which the waves are excited. It is difficult, however, to suppose such waves in the real atmosphere. Figure 9. Meridional distributions of zonal mean temperature deviation at altitudes of 43 47 km obtained for W2 (solid line), M2 (dashed line), and S2 (dotted line). 6of8

Figure 10. Schematic illustration of acceleration mechanism of mean zonal flow by the thermal tides (a) in the work of Fels and Lindzen [1974] and (b) in the present model. The dashed lines denote zonal phase velocity of the solar heating. The short arrow under the ground represents the surface friction force exerted on the mean flow adjacent to the surface. 4.2. Initial Development of Mean Zonal Flow [27] Since the basic state used in the linear calculation of Takagi and Matsuda [2006] is a mean zonal flow based on the observations, their result cannot be applied directly to the initial stages of our time integration, which is started from a resting atmosphere. It is confirmed by linear calculations with a basic atmosphere at rest that both the diurnal and semidiurnal tides excited in the cloud layer cannot propagate to the ground (not shown), as expected from the dispersion relation of the internal gravity wave. This result implies that the initial development of the mean zonal flow must be caused by other processes except the thermal tides, which are stationary observed from the sun. [28] Figure 11 shows time height distributions of the temperature deviations and the mean zonal flow on the equator in the initial evolution obtained in the case of M2. It is seen that two wave packets propagate downward from 50 km altitude to the ground in the period of 30 100 Earth days. As a result, the mean zonal flow retrograde to the planetary rotation is induced at altitudes of 0 5 km at about 100 Earth days and the net angular momentum held by the atmosphere increases by the action of surface friction. After 100 Earth days, wave packets propagate to the ground one after another and the mean zonal flow retrograde to the planetary rotation is reinforced at altitudes of 0 10 km. This result indicates that the downward propagating transient waves have important roles on the generation of the mean zonal flow in the initial stage. 5. Conclusion and Remarks [29] A full nonlinear dynamical model on the sphere is constructed to investigate the generation and maintenance of the Venus atmospheric superrotation by the thermal tides. In order to focus on the generation mechanism proposed by Takagi and Matsuda [2006], numerical experiments are conducted by using the solar heating without the mean zonal component. The results show that the atmospheric superrotation similar to the observations is reproduced at altitudes of 0 80 km. The superrotation is generated and maintained by the mechanism summarized as follows. The thermal tides excited in the cloud layer propagate downward to the ground. The momentum transport associated with the downward propagating thermal tides induces the mean zonal flow retrograde to the planetary rotation at altitudes of 0 10 km, on which the surface friction acts. As a result, the net momentum required for the generation of the superrotation is pumped up from the solid part to the atmosphere. It is confirmed that the theory proposed by Takagi and Matsuda [2006] may work and the atmospheric superrotation is generated in the nonlinear model. It should be emphasized that this mechanism does not require effects of mean meridional circulation at all. Above the cloud layer, the vertical shear of the mean zonal flow is maintained by the decelerating effect of the thermal tides which are damped by Newtonian cooling. [30] The vertical structure of the atmospheric superrotation generated by the thermal tides is not so sensitive to the vertical eddy viscosity. The maximum velocity in the equatorial cloud layer is about 60 70 m s 1 for the cases of standard Newtonian cooling. It is also confirmed that velocity of the mean zonal flow at about 40 80 km altitudes is strongly influenced by Newtonian cooling. In the weak (strong) Newtonian cooling case, the faster (slower) mean zonal flow about 80 110 (50 60) m s 1 is produced in the equatorial cloud layer. These results may be ascribed to the following points: (1) mean meridional temperature gradient is reduced by stronger Newtonian Figure 11. Time height distributions of (a) temperature deviations (K) at 0 50 km and (b) mean zonal flow (m s 1 ) at 0 10 km on the equator. Contour intervals are 1 K and 0.05 m s 1, respectively. The period is 0 200 Earth days of the initial stage. 7of8

cooling, and (2) stronger Newtonian cooling tends to inhibit the vertical propagation of thermal tides, as shown by Takagi and Matsuda [2005]. [31] The zonal mean component of solar heating is excluded in the present study. In order to demonstrate that the mechanism proposed in this study can work to generate the superrotation in the real Venus atmosphere, the present work should be extended to include the zonal mean component of the solar heating and the mean meridional circulation excited by it. It should be also noted that the radiative process in the Venus atmosphere with enormous optical depth cannot be correctly represented by a linear Newtonian cooling scheme. In order to simulate and understand the Venus atmospheric superrotation better, it is necessary to incorporate a radiation scheme suitable for the Venus atmosphere into the present model in future. [32] Acknowledgment. Figures were produced by GrADS, GNUPLOT, Tgif, and GFD-DENNOU Library. References Baker, N., and C. B. Leovy (1987), Zonal winds near Venus cloud top level: A model study of the interaction between the zonal mean circulation and the semidiurnal tide, Icarus, 69, 202 220. Covey, C., R. L. Walterscheid, and G. Schubert (1986), Dissipative tides: Application to Venus lower atmosphere, J. Atmos. Sci., 43, 3273 3278. Crisp, D. (1986), Radiative forcing of the Venus mesosphere: I. solar fluxes and heating rates, Icarus, 67, 484 514. Fels, S. B., and R. S. Lindzen (1974), The interaction of thermally excited gravity waves with mean flows, Geophys. Fluid Dyn., 6, 149 191. Gierasch, P. J. (1975), Meridional circulation and the maintenance of the Venus atmospheric rotation, J. Atmos. Sci., 32, 1038 1044. Gierasch, P. J., et al. (1997), The general circulation of the venus atmosphere: An assessment, in Venus II, edited by S. W. Bougher, D. M. Hunten, and R. J. Phillips, pp. 459 500, Univ. of Arizona Press, Tucson. Held, I. M., and M. J. Suarez (1994), A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models, Bull. Amer. Meteorol. Soc., 75, 1825 1830. Hoskins, B. J., and A. J. Simmons (1975), A multi-layer spectral model and the semi-implicit method, Quart. J. R. Meteorol. Soc., 101, 637 655. Hou, A. Y., and B. F. Farrel (1987), Superrotation induced by critical-level absorption of gravity waves on Venus, J. Atmos. Sci., 44, 1049 1061. Hou, A. Y., S. B. Fels, and R. M. Goody (1990), Zonal superrotation above Venus cloud base induced by the semidiurnal tide and the mean meridional circulation, J. Atmos. Sci., 47, 1894 1901. Lee, C., S. R. Lewis, and P. L. Read (2005), A numerical model of the atmosphere of Venus, Adv. Space Res., 36, 2142 2145. Matsuda, Y. (1980), Dynamics of the four-day circulation in the Venus atmosphere, J. Meteorol. Soc. Japan, 58, 443 470. Matsuda, Y. (1982), A further study of dynamics of the four-day circulation in the Venus atmosphere, J. Meteorol. Soc. Japan, 60, 245 254. Newman, M., and C. Leovy (1992), Maintenance of strong rotational winds in Venus middle atmosphere by thermal tides, Science, 257, 647 650. Newman, M., G. Schubert, A. J. Kliore, and I. R. Patel (1984), Zonal winds in the middle atmosphere of Venus from pioneer Venus radio occultation data, J. Atmos. Sci., 41, 1901 1913. Plumb, R. A. (1975), Momentum transport by the thermal tide in the stratosphere of Venus, Quart. J. R. Meteorol. Soc., 101, 763 776. Schubert, G. (1983), General circulation and the dynamical state of the Venus atmosphere, in Venus, edited by D. M. Hunten et al., pp. 681 765, Univ. of Arizona Press, Tucson. Seiff, A., et al. (1985), Models of the structure of the atmosphere of Venus from the surface to 100 kilometers altitude, Adv. Space Res., 5(11), 3 58, doi:10.1016/0273-1177(85)90197-8. Takagi, M. (2002), Thermal tides and topographic waves in the atmosphere of Venus, Ph.D. thesis, Univ. of Tokyo. Takagi, M., and Y. Matsuda (1999), Interpretation of Thompson s mechanism for Venus atmospheric super-rotation and its extension to the spherical geometry, J. Meteorol. Soc. Japan, 77, 971 983. Takagi, M., and Y. Matsuda (2000), Stability of convection between the day and night sides on the rotating sphere, J. Meteorol. Soc. Japan, 78, 181 186. Takagi, M., and Y. Matsuda (2005), Sensitivity of thermal tides in the Venus atmosphere to basic zonal flow and Newtonian cooling, Geophys. Res. Lett., 32, L02203, doi:10.1029/2004gl022060. Takagi, M., and Y. Matsuda (2006), Dynamical effect of thermal tides in the lower Venus atmosphere, Geophys. Res. Lett., 33, L13102, doi:10.1029/ 2006GL026168. Tomasko, M. G., L. R. Doose, P. H. Smith, and A. P. Odell (1980), Measurement of the flux of sunlight in the atmosphere of Venus, J. Geophys. Res., 85, 8167 8186. Yamamoto, M., and M. Takahashi (2003a), The fully developed superrotation simulated by a general circulation model of a Venus-like atmosphere, J. Atmos. Sci., 60(3), 561 574. Yamamoto, M., and M. Takahashi (2003b), Superrotation and equatorial waves in a T21 Venus-like AGCM, Geophys. Res. Lett., 30(9), 1449, doi:10.1029/2003gl016924. Yamamoto, M., and M. Takahashi (2004), Dynamics of Venus superrotation: The eddy momentum transport processes newly found in a GCM, Geophys. Res. Lett., 31, L09701, doi:10.1029/2004gl019518. Y. Matsuda, Department of Astronomy and Earth Science, Tokyo Gakugei University, 4-4-1 Nukuikita-machi, Koganei, Tokyo 184-8501, Japan. (ymatsuda@u-gakugei.ac.jp) M. Takagi, Department of Earth and Planetary Science, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan. (takagi@eps. s.u-tokyo.ac.jp) 8of8