Atmospheric tides on Venus: III. The planetary boundary layer

Diurnal solar heating of Venus' surface produces variable temperatures, winds, and pressure gradients within a shallow layer at the botton of the atmosphere. The corresponding asymmetric mass distribution experiences a tidal torque tending to maintain Venus' slow retrograde rotation. It is shown that including viscosity in the boundary layer does not materially affect the balance of torques. On the other hand, friction between the air and ground can reduce the predicted wind speeds from ∽5 to ∽1 m/sec in the lower atmosphere, more consistent with the observations from Venus landers and descent probes. Implications for aeolian activity on Venus' surface and for future missions are discussed.     

Television imaging of fast auroral waves

Fast auroral waves have been recorded with an Image Intensifier-TV system and a new characteristic of this phenomenon has been revealed, namely, an abrupt transition from wave motion equatorwards to an E to W motion associated with an auroral band. For equatorward motion speeds ranged from 15 to 150 km/sec which is significantly lower than values previously reported. For the E to W motion speeds ranged from 20 to 90 km/sec.     

Shear Instabilities as a Probe of Jupiter's Atmosphere

Linear wave patterns in Jupiter's clouds with wavelengths strongly clustered around 300 km are commonly observed in the planet's equatorial atmosphere (F. M. Flasar and P. J. Gierasch, 1986, J. Atmos. Sci.43, 2683-2707). We propose that the preferred wavelength is related to the thickness of an unstable shear layer within the clouds (A. P. Ingersoll and D. W. Koerner 1989, Bull. Am. Astron. Soc.21, 943). We numerically analyze the linear stability of wavelike disturbances that have nonzero horizontal phase speeds in Jupiter's atmosphere and find that, if the static stability in the shear layer is very low (but still nonnegative), a deep vertical shear layer like the one measured by the Galileo probe (D. H. Atkinson et al. 1998, J. Geophys. Res.103, 22911-22928) can generate the instabilities. The fastest growing waves grow exponentially within an hour, and their wavelengths match the observations. Close to zero values of static stability that permit the growth of instabilities are within the range of values measured by the Galileo probe in a hot spot (A. Seiff et al. 1998, J. Geophys. Res.103, 22857-22889). Our model probes Jupiter's equatorial atmosphere below the cloud deck and suggests that thick regions of wind shear and low static stability exist outside hot spots.     

Global patterns in cloud forms on Mars

Mariner 9 images and all Viking orbiter images through July 1979 were searched for cloud forms. A computer-accessible catalog was assembled, consisting of a classification of cloud type (lee wav wave, for example) and properties (directionality, wavelength, for example). Lee wave directionality shows a pattern and seasonal variation at high latitudes which is consistent with predictions of theoretical modeling. Fog and haze occurrence shows no obvious correlation with water abundance or any other simple causal factor. Lee waves are rare at equatorial latitudes. Plumes (probably dust) occur preferentially at locations where strong boundary layer convection is expected.     

Acceleration of the solar wind by whistler waves from coronal holes

We investigate the possibility of an additional acceleration of the high speed solar wind by whistler waves propagating outward from a coronal hole. We consider a stationary, spherically symmetric model and assume a radial wind flow as well as a radial magnetic field. The energy equation consists of (a) energy transfer of the electron beam which excites the whistler waves, and (b) energy transfer of the whistler waves described by conservation of wave action density. The momentum conservation equation includes the momentum transfer of two gases (a thermal gas and an electron beam). The variation of the temperature is described by a polytropic law. The variation of solar wind velocity with the radial distance is calculated for different values of energy density of the whistler waves. It is shown that the acceleration of high speed solar wind in the coronal hole due to the whistler waves is very important. We have calculated that the solar wind velocity at the earth's orbit is about equal to 670 km/sec (for wave energy density about 10^?4 erg cm^?3 at 1.1R⊙). It is in approximate agreement with the observed values.     

Low energy ions in co-rotating interaction regions at 1 AU: Evidence for statistical ion acceleration

The possibility of the statistical acceleration of solar wind ions to energies above 10 keV in the vicinity of co-rotating high speed solar wind streams by scattering from hydromagnetic waves is considered. We find that this process may occur only in the compressed fast stream plasma within the interaction region between the stream interface and the trailing edge, and may account for the energetic ion enhancements observed in this region by Richardson and Zwickl (Planet. Space Sci. 32, 1179, 1984). When statistical acceleration occurs in the outer heliosphere, the accelerated ions may provide a source population for acceleration at the co-rotating reverse shock.     

The acceleration of minor ion species in the solar wind

This paper provides a comprehensive analysis of the dynamics of the flow of minor ion species in the solar wind under the combined influences of gravity, Coulomb friction (with protons), rotational forces (arising from the Sun's rotation and the interplanetary spiral magnetic field) and wave forces (induced in the minor ion flow by Alfvén waves propagating in the solar wind). It is assumed that the solar wind can be considered as a proton-electron plasma which is, to a first approximation, unaffected by the presence of minor ions. In the dense hot region near the Sun Coulomb friction accelerates minor ions outwards against the gravitational force, part of which is cancelled by the charge-separation electric field. Once the initial acceleration has been achieved, wave and rotational forces assist Coulomb friction in further increasing the minor ion speed so that it becomes comparable with, or perhaps even exceeds, the solar wind speed. A characteristic feature of the non-resonant wave force is that it tends to bring the minor ion flow into an equilibrium where the radial speed matches the Alfvén speed relative to the solar wind speed, whereas Coulomb friction and rotational forces tend to bring the flow into an equilibrium where the radial speed of the minor ions equals the solar wind speed. Therefore, provided that there is sufficient wave energy and Coulomb friction is weak, the minor ion speed can be trapped between these two speeds. This inteststing result is in qualitative agreement with observational findings to the effect that the differential flow speed between helium ions and protons is controlled by the ratio of the solar wind expansion time to the ion-proton collision time. If the thermal speeds of the protons and minor ions are small compared to the Alfvén speed, two stable equilibrium speeds can exist because the rapid decrease in the Coulomb cross-section with increasing differential flow speed allows the non-resonant wave force to balance Coulomb friction at more than one ion speed. However, it must be emphasized that resonant wave acceleration and/or strong ion partial pressure gradients are required to achieve radial speeds of minor ions in excess of the proton speed, since, as is shown in Section 4, the non-resonant wave acceleration on protons and minor ions are identical when their radial speeds are the same, with the result that, in the solar wind, non-resonant wave acceleration tends (asymptotically) to equalize minor ion and proton speeds.     

Systematic winds at heights between 350 and 675 km from analysis of the orbits of four balloon satellites

Six values of the rate of rotation of the Earth's upper atmosphere have been obtained by analysis of the orbital inclinations of four balloon satellites in the intervals just before the final decay of their orbits. The effective heights of these results range from about 350 to about 675 km. The values themselves range from 0·8 to 1·4 times the Earth's rotation and correspond to zonal wind speeds between 100 m/sec westward and 200 m/sec eastward. All the results correspond to fairly specific local times and are consistent with a diurnal wind pattern in low latitudes having a strong eastward maximum near local midnight and a lesser westward maximum near 10:00 LT. They argue against the contention of a sharp decrease in the rate with respect to that of the Earth, which is supposed to begin at about 360 km. The factors involved in the determination of these values and the method used are discussed in considerable detail.     

Hydromagnetic waves driven by velocity shear instability in the magnetospheric boundary layer

The bulk flow of the solar wind plasma in the flank-side of the magnetospheric boundary layer, where the magnetic field lines are closed, has a component transverse to the ambient field. There is quite a strong velocity shear. The theoretical model ignores inhomogeneities in the ambient field and the mass density which occur at the magnetopause on about the same length scale as that of the velocity shear.Consideration is restricted to hydromagnetic waves which have a k-vector nearly normal to the B_o-V_o plane, i.e., approximately the magnetopause surface (k_x >k_z ～ k_y′k_xL_B > 1 and L_B = 0.1 ～ 1.0 R_E where L_B is a characteristic length of the boundary layer). It is found that a long-period (T ? 40 sec) hydromagnetic wave [the Alfvén-like wave ($\text{Ω}{}_{\mathrm{<mtext>A</mtext>}}$)] driven by velocity shear instability can be excited in the shear plasma. It is also found that the group velocity of the HM-wave is directed almost along the magnetic field line and that the magnetic variance in the shear plasma tends to be parallel to the B_o-V_o plane. The velocity shear instability in the magnetospheric boundary layer is judged to be a likely source of long-period magnetic pulsations.     

Some problems of Venus' atmospheric dynamics

A brief survey is given of the observational data on wind speeds in the atmosphere of Venus, as well as results of the theory of similitude and of a scale analysis for estimation of the wind speeds and temperature contrasts. It is shown that, in the lower portion of the atmosphere, containing roughly half of the mass, the first method produces results which are in somewhat better agreement with the measurements. Measurements of the wind distribution in the atmosphere are discussed. It is shown that, in the slowly rotating atmosphere of Venus, we should expect the Solberg mechanism of inertial instability of the circulation to be effective with respect to axisymmetrical perturbations. The numerical experiments of G.P. Williams (1968, J. Atmos. Sci., 25, 34–1045; 1970, Geophys. Fluid Dyn., 1, 357–369) indicate that in this case the circulation in the meridional plane can be broken down into a series of forward and reverse cells. The existence of such cells can serve to preserve the angular momentum of the planet with its atmosphere—the total momentum of the atmospheric frictional forces against the surface should on the average equal zero. This supports the hypothesis of G. Schubert et al. (1980, J. Geophys. Res., 85, 8007–8025) concerning the multicellular structure of the meridional circulation. Data are analyzed with regard to the time variability of the circulation. If the angular momentum of Venus′ atmosphere can change by 30% (which is not excluded by the presently available data; in Earth's atmosphere seasonal variations of the momentum reach 50%), then the relative variations in the length of a Venusian day will attain 10^?3, i.e., several hours. The surface boundary layer is considered. It is shown that, due to the small transparency of the atmosphere to thermal radiation, heat transfer between the surface and the atmosphere should basically take place by turbulent heat exchange. The basic parameters of the dynamic and thermal regimes of this layer are estimated. Questions of light refraction in the boundary layer are discussed. A strict theory of refraction, developed for these conditions, confirms the preliminary rough estimates of V.I. Moroz (1976, Cosmic Res., 14, No. 5, 691–692; Space Sci. Rev., 25, 3–127), viz, that the horizon is visible on the panorama at a distance of order 100m, due to a relatively sharp negative gradient near the surface.     

Solar Wind and Chromospheric Network

A physical model of the transition region, including upflow of the plasma in magnetic field funnels that are open to the overlying corona, is presented. A numerical study of the effects of Alfvén waves on the heating and acceleration of the nascent solar wind originating in the chromospheric network is carried out within the framework of a two-fluid model for the plasma. It is shown that waves with reasonable amplitudes can, through their pressure gradient together with the thermal pressure gradient, cause a substantial initial acceleration of the wind (on scales of a few Mm) to locally supersonic flows in the rapidly expanding magnetic field trunks of the transition region network. The concurrent proton heating is due to the energy supplied by cyclotron damping of the high-frequency Alfvén waves, which are assumed to be created through small-scale magnetic activity. The wave energy flux of the model is given as a condition at the upper chromosphere boundary, located above the thin layer where the first ionization of hydrogen takes place.Among the new numerical results are the following: Alfvén waves with an assumed f ^-1 power spectrum in the frequency range from 1 to ⁴ Hz, and with an integrated mean amplitude ranging between 25 and 75 km s⁴, can produce very fast acceleration and also heating through wave dissipation. This can heat the lower corona to a temperature of 5× 10⁵ K at a height of h=12,000 km, starting from 5× 10⁴ K at h=3000 km. The resulting thermal and wave pressure gradients can accelerate the wind to speeds of up to 150 km s^-1 at h=12,000 km, starting from 20 km s^-1 at h=3000 km in a rapidly diverging flux tube. Thus the nascent solar wind becomes supersonic at heights well below the classical Parker-Type sonic point. This is a consequence of the fact that any large wave-energy flux, if it is to be conducted through the expanding funnel to the corona, implies the building-up of an associated wave-pressure gradient. Because of the diverging field geometry, this might lead to a strong initial acceleration of the flow. There is a multiplicity of solutions, depending mainly on the coronal pressure. Here we discuss two new (as compared with a static transition region model) possibilities, namely that either the flow remains supersonic or slows down abruptly by shock formation, which then yields substantial coronal heating up to the canonical 10⁶ K for the proton temperature.     

Titan global climate model: A new 3-dimensional version of the IPSL Titan GCM

We have developed a new 3-dimensional climate model for Titan’s atmosphere, using the physics of the IPSL Titan 2-dimensional climate model with the current version of the LMDZ General Circulation Model dynamical core. Microphysics and photochemistry are still computed as zonal averages. This GCM covers altitudes from surface to 500 km altitude, with barotropic waves now being resolved and the diurnal cycle included. The boundary layer scheme has been changed, yielding a strong improvement in the tropospheric zonal wind profile modeled at Huygens descent position and season. The potential temperature profile is fairly consistent with Huygens observations in the lowest 10 km. The latitudinal profile of the near-surface temperature is close to observed values. The minimum of zonal wind observed by the Huygens probe just above the tropopause is also present in these simulations, and its origin is discussed by comparing solar heating and dynamical transport of energy. The stratospheric temperature and wind fields are consistent with our previous works. Compared to observations, the zonal wind peak is too weak (around 120 m/s) and too low (around 200 km). The temperature structures appear to be compressed in altitude, and depart strongly from observations in the upper stratosphere. These discrepancies are correlated, and most probably related to the altitude of the haze production. The model produces a detached haze layer located more than 150 km lower than observed by the Cassini instruments. This low production altitude is due to the current position of the GCM upper boundary. However, the temporal behaviour of the detached haze layer in the model may explain the seasonal differences observed between Cassini and Voyager 1. The waves present in the GCM are analyzed, together with their respective roles in the angular momentum budget. Though the role of the mean meridional circulation in momentum transport is similar to previous work, and the transport by barotropic waves is clearly seen in the stratosphere, a significant part of the transport at high latitudes is done all year long through low-frequency tropospheric waves that may be baroclinic waves.     

Estimates of the wind speeds required for particle motion on Mars

We have obtained estimates of the threshold wind speed V_gt near the top of the atmospheric boundary layer on Mars and of the rotation angle α between this wind velocity and the direction of the surface stress. this calculation has been accomplished by combining wind tunnel determinations of the friction velocity with semi-empirical theories of the Earth's atmospheric boundary layer. Calculations have been performed for a variety of values of the surface pressure, ground temperature, roughness height, boundary layer height, atmospheric composition atmospheric stability, particle density, particle diameter, and strength of the cohesive force between the particles.The curve of threshold wind speed as a function of particle diameter monotonically decreases with decreasing particle diameter for a cohesionless soil but has the classical $\text{U}$ shape for a soil with cohesion. Observational data indicate that the latter condition holds on Mars. Under “favorable” conditions minimum threshold wind speeds between about 50 and 100m/sec are required to cause particle motion. These minimum values lie close to the highest wind speeds predicted by general circulation models. Hence, particle motion should be an infrequent occurence and should be strongly correlated with nearness to small topographic features. The latter prediction is in accord with the correlation found between albedo markings and topographic obstacles such as craters. For equal wind speeds at the midpoint of the boundary layer, particle movement occurs more readily in general at night than during the day, more readily in the winter polar areas than the equatorial areas noon, and more readily for ice particles than for silicate particles.The boundary between saltating and suspendable particles is located at a particle diameter of about 100 μm. This value is close to the diameter at which the V_gt curve has its minimum. Hence, the wind can set directly into motion both saltating and larger-sized suspendable particles, but dust-storm-sized particles usually require impact by a saltating particle for motion to be initiated. Albedo changes occur most often in regions containing a mixture of dust-stoorm-sized particles and saltating particles. The threshold wind speed for surfaces containing large, nonerodible roughness elements can either be larger or smaller than the value for surfaces with only erodible material. The former condition for V_gt holds when the roughness height z₀ is less than about 1 cm and may be illustrated by craters that have experienced less erosion than their environs. The latter condition for V_gt may be partly responsible for albedo changes detected on the elevated shield volcano, Pavonis Mons. Values of the angle α generally lie between 10 and 30°. These figures place a modest limitation on the utility of surface albedo streaks as wind direction indicators.     

Ionospheric winds in the F-region and their effects on the limiting periods of gravity waves

Spectrum analyses of ionospheric electron density and content fluctuations show periods with a lower limit near 5 min. Interpretation of this cut off in terms of gravity waves in a windless atmosphere leads to unacceptably low thermospheric temperatures near 180°K. It is concluded that neutral winds reduce the apparent cut-off period in the ionosphere. The maximum horizontal wind speed obtained from cut-off data is about 100 m/sec.     

Eddy mixing and circulation at ionospheric levels

The rates of heat input into the mesosphere and lower thermosphere are calculated and compared with the heat losses. The worldwide average eddy diffusion coefficient required to maintain continuity in the heat budget is calculated and found to vary from about 10⁷ cm²/sec at 120 km down to about 10⁵ cm²/sec at 60 km. From the global asymmetry in heating at the solstice, it is concluded that a systematic pattern of vertical velocities prevails ranging from less than 1 cm/sec in the mesosphere up to 10 cm/sec near 120 km, upward over the summer polar region and downward over the winter polar region. This can be balanced by a wind system towards the winter polar region with velocities near 1 m/sec at 60 km increasing to 30 m/sec at 120 km. Such a wind system provides an explanation for the helium bulge in the upper thermosphere over the winter polar region.     

Analysis of the orbit of cosmos 268 rocket (1969-20B)

The orbit of Cosmos 268 rocket (1969-20B) has been determined at 28 epochs during its 342-day life, with the aid of the PROP5 orbit refinement program. All available observations were used, including 16 from the Hewitt camera at Malvern, 28 from the 200-mm camera at Meudon, 56 from the kinetheodolite at the Cape Observatory, 700 visual observations from volunteer observers, 500 US Navy observations and 200 British radar observations. The orbits are of very good accuracy for such a high-drag satellite, most of the values of inclination having standard deviations less than 0.002°. The most accurate orbits are those utilizing photographic observations, and the best of these has standard deviations of 0.00001 in eccentricity and 0.0001° in inclination.

The values of inclination obtained, after correction to allow for the effects of other perturbing forces, have been analysed to determine zonal wind speeds in the upper atmosphere at heights a little above perigee (230–250 km) averaged over latitudes up to about 25°. The results show a clear distinction between the wind at night (21 to 03 hr local time), which is west-to-east with an average speed of 140 ± 50 m/sec, and the wind by day (08 to 17 hr), which is east-to-west with an average speed of 110 ± 50 m/sec.     

Propagation characteristics of hydromagnetic waves in various layers of the ionosphere

This paper derives the basic propagation characteristics of hydromagnetic waves in various layers of the ionosphere. It is shown that propagation in the upper ionosphere and the F₂ layer is largely isotropic. In the lower region of the ionosphere there are two possible modes of propagation, both being anisotropic. Propagation characteristics of waves in this lower region, however, are relatively independent of the direction of horizontal propagation. Calculations of intrinsic wave attenuation show that ducted propagation of Pc 1 signals over appreciable horizontal distances may only take place in the upper layers of the ionosphere.     

Impacts into oceans and seas

Impacts of cosmic bodies into oceans and seas lead to the formation of very high waves. Numerical simulations of 3-km and 1-km comets impacting into a 4 km depth ocean with a velocity of 20 km/sec have been conducted. For a 1-km body, depth of the interim crater in the sea bed is about 8 km below ocean level, and the height of the water wave is 10 m at a distance of 2000 km from the impact point. As the water wave runs into shallows, a huge tsunami hits the coast. The height of the wave strongly depends on the coastal and sea bed topography.If the impact occurred near the shore, the huge mass of water strikes the cliffs and the near shore mountain ridges and can cause displacement of the rocks, initiate landslides, and change the relief. Thus, impact into oceans and seas is an important geological factor.Cosmic bodies of small sizes are disrupted by aerodynamic forces. Fragments of a 100-m radius comet striking the water surface create an unstable cavity in the water of about 1 km radius. Its collapse also creates tsunami.A simple estimate has been made using the light curves from recent atmosphere explosions detected by satellites. The results of our assessment of the characteristics of meteoroids which caused these intense light flashes suggests that fragments of a 25-m stony body with initial impact velocity 15 to 20 km/sec will hit the surface. For a 75-m iron body striking the sea with a depth of 600 m, the height of the wave is 10 m at 200–300 km distance from the impact.     

Effects of Steady Flow on Magnetoacoustic-Gravity Surface Waves: I. The Weak Field Case

Magnetoacoustic gravity (MAG) waves have been studied for some time. In this article, we investigate the effect that a shear flow at a tangential discontinuity embedded in a gravitationally stratified and magnetised plasma has on MAG surface waves. The dispersion relation found is algebraically analogous to the relation of the non-flow cases obtained by Miles and Roberts (Solar Phys. 141, 205, 1992), except for the introduction of a Doppler-shifted frequency for the eigenvalue. This feature, however, introduces rather interesting physics, including the asymmetric presence of forward- and backward-propagating surface waves. We find that increasing the equilibrium flow speed leads to a shift in the permitted regions of propagation for surface waves. For most wave number combinations this leads to the fast mode being completely removed, as well as more limited phase speed regimes for slow-mode propagation. We also find that upon increasing the flow, the phase speeds of the backward propagating waves are increased. Eventually, at high enough flow speeds, the wave’s direction of propagation is reversed and is in the positive direction. However, the phase speed of the forward-propagating wave remains mainly the same. For strong enough flows we find that the Kelvin–Helmholtz instability can also occur when the forward- and backward-propagating modes couple.     

Small-scale temperature fluctuations seen by the VeRa Radio Science Experiment on Venus Express

The Venus Express Radio Science Experiment VeRa retrieves atmospheric profiles in the mesosphere and troposphere of Venus in the approximate altitude range of 40–90 km. A data set of more than 500 profiles was retrieved between the orbit insertion of Venus Express in 2006 and the end of occultation season No. 11 in July 2011. The atmospheric profiles cover a wide range of latitudes and local times, enabling us to study the dependence of vertical small-scale temperature perturbations on local time and latitude.Temperature fluctuations with vertical wavelengths of 4 km or less are extracted from the measured temperature profiles in order to study small-scale gravity waves. Significant wave amplitudes are found in the stable atmosphere above the tropopause at roughly 60 km as compared with the only shallow temperature perturbations in the nearly adiabatic region of the adjacent middle cloud layer, below.Gravity wave activity shows a strong latitudinal dependence with the smallest wave amplitudes located in the low-latitude range, and an increase of wave activity with increasing latitude in both hemispheres; the greatest wave activity is found in the high-northern latitude range in the vicinity of Ishtar Terra, the highest topographical feature on Venus.We find evidence for a local time dependence of gravity wave activity in the low latitude range within ±30° of the equator. Gravity wave amplitudes are at their maximum beginning at noon and continuing into the early afternoon, indicating that convection in the lower atmosphere is a possible wave source.The comparison of the measured vertical wave structures with standard linear-wave theory allows us to derive rough estimates of the wave intrinsic frequency and horizontal wavelengths, assuming that the observed wave structures are the result of pure internal gravity waves. Horizontal wavelengths of the waves at 65 km altitude are on the order of ≈300–450 km with horizontal phase speeds of roughly 5–10 m/s.