Instability of Zonal Flows in Rotating Spherical Shells: An Application to Jupiter

Measurements from the Galileo probe suggest that the zonal winds are deep rooted. Jupiter's high rotation rate makes it likely that the whole outer molecular H/He layer is involved in these long-lived jet flows. Assuming that the primary flows are geostrophic, and that the banded surface structure stretches right through the molecular H/He layer, we examine the conditions for such flows to be stable. As a first step, the linear stability of some prescribed banded zonal flows in a rotating spherical shell is explored. Incompressibility is assumed for simplicity, and the boundary condition is stress-free. We compare solutions for two aspect ratios, appropriate for the molecular H/He layers of Jupiter and Saturn, and two Ekman numbers (E=10⁻² and E=10⁻⁴). Convective and shear flow instabilities compete in our system. The convective instabilities are of the well-known columnar structure. Shear flow instabilities for the larger Ekman number are similar to the Taylor-Couette instability in rotating annuli. At the lower Ekman number, shear flow instabilities adopt a geostrophic character, assuming the form of rotating columns, similar to the convective instabilities. While the convective instability always sets in outside the tangent cylinder, shear instability can become unstable inside the tangent cylinder. If even a weak zonal flow is present inside the tangent cylinder, the flow is unstable to shear instability. This offers an explanation why the jovian zonal jet structure is much weaker at the higher latitudes that correspond to locations inside the tangent cylinder.     

Capture probability in the 3:1 mean motion resonance with Jupiter: an application to the Vesta family

We study the capture and crossing probabilities in the 3:1 mean motion resonance with Jupiter for a small asteroid that migrates from the inner to the middle Main Belt under the action of the Yarkovsky effect. We use an algebraic mapping of the averaged planar restricted three-body problem based on the symplectic mapping of Hadjidemetriou (Celest Mech Dyn Astron 56:563–599, 1993), adding the secular variations of the orbit of Jupiter and non-symplectic terms to simulate the migration. We found that, for fast migration rates, the captures occur at discrete windows of initial eccentricities whose specific locations depend on the initial resonant angles, indicating that the capture phenomenon is not probabilistic. For slow migration rates, these windows become narrower and start to accumulate at low eccentricities, generating a region of mutual overlap where the capture probability tends to 100 %, in agreement with the theoretical predictions for the adiabatic regime. Our simulations allow us to predict the capture probabilities in both the adiabatic and non-adiabatic cases, in good agreement with results of Gomes (Celest Mech Dyn Astron 61:97–113, 1995) and Quillen (Mon Not RAS 365:1367–1382, 2006). We apply our model to the case of the Vesta asteroid family in the same context as Roig et al. (Icarus 194:125–136, 2008), and found results indicating that the high capture probability of Vesta family members into the 3:1 mean motion resonance is basically governed by the eccentricity of Jupiter and its secular variations.     

Pressure-induced absorption by H2 in the atmosphere of Jupiter and Saturn

The S(1) line of the pressure-induced fundamental band of H₂ was identified and measured in the spectra of Saturn and Jupiter. This broad line at 4750 cm^?1 lies in a region free from telluric and planetary absorptions. It is about 99% absorbing in the core; the high-frequency wing extends to at least 5100 cm^?1. We compare the obseved line shape to the predictions of both a reflecting-layer model (RLM) and a homogeneous scattering model (HSM). The RLM provides a good fit to the Saturn line profile for temperatures near 150K; the derived base-level density is 0.52 (+0.26, ?0.17) amagat and the H₂ abundance is 25 (+10, ?9) km-amagat, assuming a scale height of 48 km. The Jupiter line profile is fit by both the RLM and HSM, but for widely differing temperatures, neither of which seems probable. The precise fitting of the observed S(1) line profile to computed models depends critically on the determination of the true continuum level; difficulties encountered in finding the continuum, especially for Jupiter, are discussed. Derived RLM densities and abundances for both planets are substantially lower than those derived from RLM analyses of the H₂ quadrupole lines, the 3ν₃ band of CH₄, and from other sources.     

A Kinetic Model of Coronal Heating and Acceleration by Ion-Cyclotron Waves: Preliminary Results

We present a kinetic model of the heating and acceleration of coronal protons by outward-propagating ion-cyclotron waves on open, radial magnetic flux tubes. In contrast to fluid models which typically insist on bi-Maxwellian distributions and which spread the wave energy and momentum over the entire proton population, this model follows the kinetic evolution of the collisionless proton distribution function in response to the combination of the resonant wave-particle interaction and external forces. The approximation is made that pitch-angle scattering by the waves is faster than all other processes, resulting in proton distributions which are uniform over the resonant surfaces in velocity space. We further assume, in this preliminary version, that the waves are dispersionless so these resonant surfaces are portions of spheres centered on the radial sum of the Alfvén speed and the proton bulk speed. We incorporate the fact that only those protons with radial speeds less than the bulk speed will be resonant with outward-propagating waves, so this rapid interaction acts only on the sunward half of the distribution. Despite this limitation, we find that the strong perpendicular heating of the resonant particles, coupled with the mirror force, results in substantial outward acceleration of the entire distribution. The proton distribution evolves towards an incomplete shell in velocity space, and appears vastly different from the distributions assumed in fluid models. Evidence of these distinctive distributions should be observable by instruments on Solar Probe.     

Methane absorption in the atmosphere of Jupiter from 1800 to 9500 cm and implications for vertical cloud structure

New measurements of the low-temperature near-infrared absorption of methane (Sihra, 1998, Laboratory measurements of near-infrared methane bands for remote sensing of the jovian atmosphere, Ph.D. thesis, University of Oxford) have been combined with existing, longer path-length, higher-temperature data of Strong et al. (1993, Spectral parameters of self- and hydrogen-broadened methane from 2000 to 9500 cm⁻¹ for remote sounding of the atmosphere of Jupiter, J. Quant. Spectrosc. Radiat. Trans. 50, 309-325) and fitted with band models. The combined data set is found to be more consistent with previous low-temperature methane absorption measurements than that of Strong et al. (1993, J. Quant. Spectrosc. Radiat. Trans. 50, 309-325) but covers the same wider wavelength range and accounts for both self- and hydrogen-broadening conditions. These data have been fitted with k-coefficients in the manner described by Irwin et al. (1996, Calculated k-distribution coefficients for hydrogen- and self-broadened methane in the range 2000-9500 cm⁻¹ from exponential sum fitting to band modelled spectra, J. Geophys. Res. 101, 26,137-26,154) and have been used in multiple-scattering radiative transfer models to assess their impact on our previous estimates of the jovian cloud structure obtained from Galileo Near-Infrared Mapping Spectrometer (NIMS) observations (Irwin et al., 1998, Cloud structure and atmospheric composition of Jupiter retrieved from Galileo NIMS real-time spectra, J. Geophys. Res. 103, 23,001-23,021; Irwin et al., 2001, The origin of belt/zone contrasts in the atmosphere of Jupiter and their correlation with 5-μm opacity, Icarus 149, 397-415; Irwin and Dyudina, 2002, The retrieval of cloud structure maps in the equatorial region of Jupiter using a principal component analysis of Galileo/NIMS data, Icarus 156, 52-63). Although significant differences in methane opacity are found at cooler temperatures, the difference in the optical depth of the atmosphere due to methane is found to diminish rapidly with increasing pressure and temperature and thus has negligible effect on the cloud structure inferred at deeper levels. Hence the main cloud opacity variation is still found to peak at around 1-2 bar using our previous analytical approach, and is thus still in disagreement with Galileo Solid State Imager (SSI) determinations (Banfield et al., 1998, Jupiter's cloud structure from Galileo imaging data, Icarus 135, 230-250; Simon-Miller et al., 2001, Color and the vertical structure in Jupiter's belts, zones and weather systems, Icarus 154, 459-474) which place the main cloud deck near 0.9 bar. Further analysis of our retrievals reveals that this discrepancy is probably due to the different assumptions of the two analyses. Our retrievals use a smooth vertically extended cloud profile while the SSI determinations assume a thin NH₃ cloud below an extended haze. When the main opacity in our model is similarly assumed to be due to a thin cloud below an extended haze, we find the main level of cloud opacity variation to be near the 1 bar level—close to that determined by SSI and moderately close to the expected condensation level of ammonia ice of 0.85 bar, assuming that the abundance of ammonia on Jupiter is (7±1)×¹⁰⁻⁴ (Folkner et al., 1998, Ammonia abundance in Jupiter's atmosphere derived from the attenuation of the Galileo probe's radio signal, J. Geophys. Res. 103, 22,847-22,855; Atreya et al., 1999, A comparison of the atmospheres of Jupiter and Saturn: deep atmospheric composition, cloud structure, vertical mixing, and origin, Planet. Space Sci. 47, 1243-1262). However our data in the 1-2.5 μm range have good height discrimination and our lowest estimate of the cloud base pressure of 1 bar is still too great to be consistent with the most recent estimates of the ammonia abundance of 3.5 × solar. Furthermore the observed limited spatial distribution of ammonia ice absorption features on Jupiter suggests that pure ammonia ice is only present in regions of localised vigorous uplift (Baines et al., 2002, Fresh ammonia ice clouds in Jupiter: spectroscopic identification, spatial distribution, and dynamical implications, Icarus 159, 74-94) and is subsequently rapidly modified in some way which masks its pure absorption features. Hence we conclude that the main cloud deck on Jupiter is unlikely to be composed of pure ammonia ice and instead find that it must be composed of either NH₄SH or some other unknown combination of ammonia, water, and hydrogen sulphide and exists at pressures of between 1 and 2 bar.     

从仰望星空到走向太空——纪念伽利略用天文望远镜进行天文观测400年

2009年是伽利略(Galileo Galilei,1564—1642)首次用望远镜观测天体400周年,因此被联合国确定为国际天文年,以纪念这位人类历史上第一个把望远镜对准茫茫太空的人。伽利略是近代科学的开创者之一,是科学史上的伟人。他把理论与实验相结合,形成了一套基于实验观察、数学分析、严谨实证的科学研究方法,从此人类有了现代意义上的科学。伽利略等人所开创的近现代科学,今天更加充满生机,有力地推动着人类文明的进步与发展。     

Yardangs in terrestrial ignimbrites: Synergistic remote and field observations on Earth with applications to Mars

The conditions of formation and the form of yardangs in ignimbrites in the Central Andes of Chile, Bolivia, and Argentina may be the most convincing terrestrial analog to the processes and lithology that produce the extensive yardangs of the Medusae Fossae Formation (MFF) of Mars. Through remote and field study of yardang morphologies in the Central Andes we highlight the role that variable material properties of the host lithology plays in their final form. Here, ignimbrites typically show two main facies: an indurated and jointed facies, and a weakly to poorly indurated, ash- and pumice-rich facies. Both facies are vertically arranged in large (erupted volume >100's of km³) ignimbrites resulting in a resistant capping layer, while smaller (10's of km³) ignimbrites are made predominantly of the weakly indurated facies. The two facies have quite different mechanical properties; the indurated facies behaves as strong rock, fails by block collapse and supports steep/vertical cliffs, while the non-indurated facies is more easily eroded and forms gentle slopes and manifests as more subdued erosional forms. In response to aeolian action, the presence of an upper indurated facies results in large, elongate, high aspect ratio (1:20-1:40) megayardangs that form tall (100 m), thin ridges with steep to vertical walls. These are built on a broad apron of the weakly indurated facies with abundant fallen blocks from the upper indurated facies. These terrestrial megayardangs appear to be analogous to megayardangs with associated block fields seen on Mars. Smaller-volume, weakly indurated ignimbrites are sculpted into smaller, stubbier forms with aspect ratios of 1:5-1:10 and heights rarely exceeding 10 m. Excavation of a windward basal moat suggests an erosional progression like that seen in incipient yardangs on Mars. Excavation rates of 0.007-0.003 cm/year are calculated for the weakly indurated ignimbrites. While a persistent strong unidirectional wind is the dominant parameter controlling yardang formation and orientation, a role for flow separation and vorticity is also suggested by our observations at both yardang types. While the indurated facies is commonly pervasively jointed, jointing is of secondary importance in controlling yardang orientation. Serrated margins, a common feature on Mars, result from oblique intersections of jointing with yardang flanks or scarps of ignimbrite. The processes of yardang formation we describe from ignimbrites from the Central Andes are not necessarily specific to ignimbrites, but do connote that degree and distribution of induration is a major control in yardang formation and this has implications for the lithology of the MFF on Mars.     

Airborne spectroscopy of Jupiter in the 100- to 300-cm−1 region: Global properties of ammonia gas and ice haze

A spectrum of the disk of Jupiter was obtained in January 1978 from the Kuiper Airborne Observatory, covering the 100- to 300-cm^?1 spectral range at a resolution corresponding to 1.65 cm^?1. Although taken more than a year before the Voyager 1 Jupiter encounter, this spectrum serves to extend the Voyager IRIS experiment coverage down from its lower limit of 200 cm^?1. Analysis of the spectrum provides information on global mean properties of ammonia gas and an ammonia ice haze. A vertical distribution indistinguishable from saturation equilibrium, with a sharp depletion near the temperature minimum, matches the observed shape of the rotational line absorption best. Constraints on the total optical thickness of the ammonia ice haze can be made, but other properties, such as particle size or vertical scale height, cannot be distinguished clearly from our data in this spectral region. Nevertheless, all models of the haze produce a “continuum” thermal emission between the NH₃ line manifolds which is much lower than that produced by the H₂ collision-induced dipole opacity.     

Computer Simulation of Sputtering of Lunar Regolith by Solar Wind Protons: Contribution to Change of Surface Composition and to Hydrogen Flux at the Lunar Poles

A computer simulation of the sputtering of lunar soil by solar wind protons was performed with the TRIM program. The rate of the sputtering-induced erosion of regolith particles was shown to be less than 0.2 Å per year. A preferential sputtering of Ca, Mg, and O was found along with a less intense sputtering of Fe, Si, and Ti. However, with no other selection mechanisms, surface concentrations of the atoms would differ from the volume ones by no more than 6 %. The enrichment of rims of regolith particles with iron occurs as a result of selective removal of lighter atoms from the lunar surface because of different energies of escape from the Moon's gravity. The energy distributions proved to be the same for all sorts of the sputtered atoms, except for implanted hydrogen; thus, a greater fraction of the atoms left on the lunar surface corresponds to heavier elements. According to simulation results, the concentration of reduced iron observed in the mature regolith could be attained during the time of regolith particle exposure to the present flux of solar wind (10⁵ years). Thus, sputtering can provide the concentration of Fe₀ observed in regolith. On periphery of a cloud of impact vapor the temperature is too low for an irreversible selective removal of evaporation products; thus, a meteoritic bombardment contributes to the formation of composition of the rims of regolith particles mainly through enrichment of the rims with elements from the bulk of the particles. The estimates of fluxes of backscattered solar wind protons and of sputtered protons, earlier implanted to the regolith, demonstrated that their contribution to the proton flux near the poles is only 10⁴ cm^–2 s^–1. This is by two orders of magnitude smaller than the proton flux from the Earth's magnetosphere which is, therefore, the main source of protons for permanently shaded polar craters of the Moon.     

Proton and He^2＋ Temperature Anisotropies in the Solar Wind Driven by Ion Cyclotron Waves

We carried out one-dimensional hybrid simulations of resonant scattering of protons and He^2 ions by ion cyclotron waves in an initially homogeneous, collisionless and magnetized plasma. The initial ion cyclotron waves have a power spectrum and propagate both outward and inward. Due to the resonant interaction with the protons and He^2 ions, the wave power will be depleted in the resonance region. Both the protons and He^2 ions can be resonantly heated in the direction perpendicular to the ambient magnetic field and leading to anisotropic velocity distributions, with the anisotropy higher for the He^2 ions than for the protons. At the same time, the anisotropies of the protons and He^2 ions are inversely correlated with the plasma β‖p=8πnpkBT‖p/B0^2, consistent with the prediction of the quasilinear theory (QLT).     

The dynamics of a rapidly escaping atmosphere: Applications to the evolution of Earth and Venus

A simple, idealized model for the rapid escape of a hydrogen thermosphere provides some quantitative estimates for the energy-limited flux of escaping particles. The model assumes that the atmosphere is “tightly bound” by the gravitational field at lower altitudes, that diffusion through the lower atmosphere does not limit the flux, and that the main source of heating is solar euv. Rather low thermospheric temperatures are typical of such escape and a characteristic minimum develops in the temperature profile as the escape flux approaches its maximum possible value. The flux is limited by the amount of euv energy absorbed, which is in turn controlled by the radial extent of the thermosphere. Regardless of the amount of hydrogen in the thermosphere, the low temperatures accompanying rapid escape limit its extent, and thus constrain the flux. Applied to the Earth and Venus, the results suggest that the escape of hydrogen from these planets would have been energy-limited if their primordial atmospheres contained total hydrogen mixing ratios exceeding a few percent. Hydrogen and deuterium may have been lost in bulk, but heavier elements would have remained in the atmosphere. These results place constraints on hypotheses for the origin of the planets and their subsequent evolution.     

Coronal Shock Waves,EUV Waves,and Their Relation to CMEs. III. Shock-Associated CME/EUV Wave in an Event with a Two-Component EUV Transient

On 17 January 2010, STEREO-B observed in extreme ultraviolet (EUV) and white light a large-scale dome-shaped expanding coronal transient with perfectly connected off-limb and on-disk signatures. Veronig et al. (Astrophys. J. Lett. 716, L57, 2010) concluded that the dome was formed by a weak shock wave. We have revealed two EUV components, one of which corresponded to this transient. All of its properties found from EUV, white light, and a metric type II burst match expectations for a freely expanding coronal shock wave, including correspondence with the fast-mode speed distribution, while the transient sweeping over the solar surface had a speed typical of EUV waves. The shock wave was presumably excited by an abrupt filament eruption. Both a weak shock approximation and a power-law fit match kinematics of the transient near the Sun. Moreover, the power-law fit matches the expansion of the CME leading edge up to 24 solar radii. The second, quasi-stationary EUV component near the dimming was presumably associated with a stretched CME structure; no indications of opening magnetic fields have been detected far from the eruption region.     

Diffusion and dispersion in magnetohydrodynamic turbulence: The influence of mean magnetic fields

The properties of magnetohydrodynamic (MHD) turbulence under the influence of a strong mean magnetic field are investigated from the Lagrangian viewpoint by tracking fluid particles in direct numerical simulations. The particle trajectories show characteristic bends near vortex sheets. A strong mean magnetic field leads to preferential diffusion parallel to the mean magnetic field. The two‐particle relative dispersion process shows a dependence on the orientation of the initial separation vector. The relative dispersion is slowed down for initial separation vectors aligned with the mean magnetic field. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Terrestrial planet and asteroid formation in the presence of giant planets. I. Relative velocities of planetesimals subject to Jupiter and Saturn perturbations.

We investigate the orbital evolution of 10(13)- to 10(25) -g planetesimals near 1 AU and in the asteroid belt (near 2.6 AU) prior to the stage of evolution when the mutual perturbations between the planetesimals become important. We include nebular gas drag and the effects of Jupiter and Saturn at their present masses and in their present orbits. Gas drag introduces a size-dependent phasing of the secular perturbations, which leads to a pronounced dip in encounter velocities (Venc) between bodies of similar mass. Plantesimals of identical mass have Venc approximately 1 and approximately 10 m s-1 (near 1 and 2.6 AU, respectively) while bodies differing by approximately 10 in mass have Venc approximately 10 and approximately 100 m s-1 (near 1 and 2.6 AU, respectively). Under these conditions, growth, rather than erosion, will occur only by collisions of bodies of nearly the same mass. There will be essentially no gravitational focusing between bodies less than 10(22) to 10(25) g, allowing growth of planetary embryos in the terrestrial planet region to proceed in a slower nonrunaway fashion. The environment in the asteroid belt will be even more forbidding and it is uncertain whether even the severely depleted present asteroid belt could form under these conditions. The perturbations of Jupiter and Saturn are quite sensitive to their semi-major axes and decrease when the planets' heliocentric distances are increased to allow for protoplanet migration. It is possible, though not clearly demonstrated, that this could produce a depleted asteroid belt but permit formation of a system of terrestrial planet embryos on a approximately 10(6)-year timescale, initially by nonrunaway growth and transitioning to runaway growth after approximately 10(5) years. The calculations reported here are valid under the condition that the relative velocities of the bodies are determined only by Jupiter and Saturn perturbations and by gas drag, with no mutual perturbations between planetesimals. If, while subject to these conditions, the bodies become large enough for their mutual perturbations to influence their velocity and size evolution significantly, the problem becomes much more complex. This problem is under investigation.     

星云中继假说: 星云中的原始生命和地球生命的起源

提出了关于地球生命起源的新模型---星云中继假说, 它是宇宙胚种论的修改版本. 在这个模型中, 作为宇宙``种子''的原始生命起源于太阳系的前身恒星系统中的生物化学过程, 并且在前身恒星死亡后充满整个原太阳星云. 地球生命的起源可以分为3个阶段: 太阳前身恒星的原始生命起源, 原太阳星云时期和太阳系形成与地球生命时期. 这个模型最主要的推论是原始生命(或其后裔)以及它们的化石存在于太阳系内各种天体之中.     

The Whole Heliosphere Interval in the Context of a Long and Structured Solar Minimum: An Overview from Sun to Earth

Throughout months of extremely low solar activity during the recent extended solar-cycle minimum, structural evolution continued to be observed from the Sun through the solar wind and to the Earth. In 2008, the presence of long-lived and large low-latitude coronal holes meant that geospace was periodically impacted by high-speed streams, even though solar irradiance, activity, and interplanetary magnetic fields had reached levels as low as, or lower than, observed in past minima. This time period, which includes the first Whole Heliosphere Interval (WHI 1: Carrington Rotation (CR) 2068), illustrates the effects of fast solar-wind streams on the Earth in an otherwise quiet heliosphere. By the end of 2008, sunspots and solar irradiance had reached their lowest levels for this minimum (e.g., WHI 2: CR 2078), and continued solar magnetic-flux evolution had led to a flattening of the heliospheric current sheet and the decay of the low-latitude coronal holes and associated Earth-intersecting high-speed solar-wind streams. As the new solar cycle slowly began, solar-wind and geospace observables stayed low or continued to decline, reaching very low levels by June??C?July 2009. At this point (e.g., WHI 3: CR 2085) the Sun?CEarth system, taken as a whole, was at its quietest. In this article we present an overview of observations that span the period 2008??C?2009, with highlighted discussion of CRs 2068, 2078, and 2085. We show side-by-side observables from the Sun??s interior through its surface and atmosphere, through the solar wind and heliosphere and to the Earth??s space environment and upper atmosphere, and reference detailed studies of these various regimes within this topical issue and elsewhere.     

Electron acceleration in solar-flare magnetic traps: Model properties and their observational confirmations

Using an analytical solution of the kinetic equation, we have investigated the model properties of the coronal and chromospheric hard X-ray sources in the limb flare of July 19, 2012. We calculated the emission spectrum at the flare loop footpoints in the thick-target approximation with a reverse current and showed it to be consistent with the observed one. The spectrum of the coronal source located above the flare loop was calculated in the thin-target approximation. In this case, the slope of the hard X-ray spectrum is reproduced very accurately, but the intensity of the coronal emission is lower than the observed one by several times. Previously, we showed that this contradiction is completely removed if the additional (relative to the primary acceleration in the reconnecting current layer) electron acceleration in the coronal magnetic trap that contracts in the transverse direction and decreases in length during the impulsive flare phase is taken into account. In this paper we study in detail this effect in the context of a more realistic flare scenario, where a whole ensemble of traps existed in the hard X-ray burst time, each of which was at different stages of its evolution: formation, collapse, destruction. Our results point not only to the existence of first-order Fermi acceleration and betatron electron heating in solar flares but also to their high efficiency. Highly accurate observations of a specific flare are used as an example to show that the previously predicted theoretical features of the model find convincing confirmations.     

Diffuse interstellar bands (DIB): co-planar doubly excited He and metal atoms embedded in Rydberg Matter

The interpretation of the more than 300 diffuse interstellar bands (DIBs) is one of the most long-standing problems in interstellar spectra since the two first bands were reported in 1921. We now predict the frequencies of 260 diffuse interstellar bands (DIBs) using the Rydberg Matter model we have developed previously. These transitions involve mainly He atoms, but other two-electron atoms like Ca and other metals can take part in the absorption processes. Approximately 70% of the total intensity of the DIBs is due to absorption in doubly excited states and 30% in singly excited He atoms. The doubly excited states are in inverted states while the He atoms are thermal. The possibilities to observe DIBs in the UV and NIR ranges are discussed and band positions are predicted.     

Influence of the largest asteroids on the orbital motions of terrestrial planets: application to the Earth and Mars

The level of precision of modern numerical ephemeris of the Solar System necessitates taking into account the gravitational influence of the largest asteroids on the terrestrial planets. This can be done in a straightforward manner when assuming that the mass of the asteroid is well known. Nevertheless, this is rarely the case, even for the largest asteroids. In this paper, we use recent determinations of the masses of Ceres, Pallas, and Vesta to both qualitatively and quantitatively determine the action of these asteroids on the orbital parameters of the Earth and Mars. This is done by the numerical integration by comparing the orbital motions of the perturbed planet when adding or not the perturbing asteroid to the classical 9 bodies problem (the Sun + the eight planets). Some preliminary results are discussed. Published in Russian in Astronomicheskii Vestnik, 2009, Vol. 43, No. 1, pp. 83–86. The text was submitted by the autors in English.