Parametric analysis of modeled ion escape from Mars     

Curtis V. Manning  Yingjuan Ma  Christopher P. McKay 《Icarus》2011,212(1):131-137

We develop a parametric fit to the results of a detailed magnetohydrodynamic (MHD) study of the response of ion escape rates (O⁺, and ) to strongly varied solar forcing factors, as a way to efficiently extend the MHD results to different conditions. We then use this to develop a second, evolutionary model of solar forced ion escape. We treat the escape fluxes of ion species at Mars as proportional to the product of power laws of four factors - that of the EUV flux R_euv, the solar wind particle density R_ρ, its velocity (squared) R_v², and the interplanetary magnetic field pressure R_B², where forcing factors are expressed in units of the current epoch-averaged values. Our parametric model is: , where ?(i) is the escape flux of ion i. We base our study on the results of just six provided MHD model runs employing large forcing factor variations, and thus construct a successful, first-order parametric model of the MHD program. We perform a five-dimensional least squares fit of this power law model to the MHD results to derive the flux normalizations and the indices of the solar forcing factors. For O⁺, we obtain the values, 1.73 × 10²⁴ s⁻¹, 0.782, 0.251, 0.382, and 0.214, for ?₀, α, β, γ, and δ, respectively. For , the corresponding values are 1.68 × 10²⁴ s⁻¹, −0.393, 0.798, 0.967, and 0.533. For , they are 8.66 × 10²² s⁻¹, −0.427, 1.083, 1.214, and 0.690. The fit reproduces the MHD results to an average error of about 5%, suggesting that the power laws are broadly representative of the MHD model results. Our analysis of the MHD model shows that by itself an increase in R_EUV enhances O⁺ loss, but suppresses the escape of and , whereas increases in solar wind (i.e., in , and R_B², with R_euv constant) favors the escape of heavier ions more than light ions. The ratios of escaping ions detectable at Mars today can be predicted by this parametric fit as a function of the solar forcing factors. We also use the parametric model to compute escape rates over martian history. This second parametric model expresses ion escape functions of one variable (per ion), ?(i) = ?₀(i)(t/t₀)^−ξ(i). The ξ(i) are linear combinations of the epoch-averaged ion escape sensitivities, which are seen to increase with ion mass. We integrate the and oxygen ion escape rates over time, and find that in the last 3.85 Gyr, Mars would have lost about mbars of , and of water (from O⁺ and ) from ion escape.  相似文献   

Mars high resolution gravity fields from MRO, Mars seasonal gravity, and other dynamical parameters     

Alex S. Konopliv  Sami W. Asmar  Özgür Karatekin  Suzanne E. Smrekar  Maria T. Zuber 《Icarus》2011,211(1):401-428

With 2 years of tracking data collection from the MRO spacecraft, there is noticeable improvement in the high frequency portion of the spherical harmonic Mars gravity field. The new JPL Mars gravity fields, MRO110B and MRO110B2, show resolution near degree 90. Additional years of MGS and Mars Odyssey tracking data result in improvement for the seasonal gravity changes which compares well to global circulation models and Odyssey neutron data and Mars rotation and precession (). Once atmospheric dust is accounted for in the spacecraft solar pressure model, solutions for Mars solar tide are consistent between data sets and show slightly larger values (k₂ = 0.164 ± 0.009, after correction for atmospheric tide) compared to previous results, further constraining core models. An additional 4 years of Mars range data improves the Mars ephemeris, determines 21 asteroid masses and bounds solar mass loss (dGM_Sun/dt < 1.6 × 10⁻¹³ GM_Sun year⁻¹).  相似文献   

Oxygen ion escape at Mars in a hybrid model: High energy and low energy ions     

Esa Kallio  Kaijun Liu  Valter Pohjola 《Icarus》2010,206(1):152-163

We have studied the escape and energization of several O⁺ populations and an population at Mars by using a hybrid model. The quasi-neutral hybrid model, HYB-Mars model, included five oxygen ion populations making it possible to distinguish photoions from oxygen ions originating from charge exchange processes and from the ionosphere.We have identified two high-energy ion components and one low-energy ion component of oxygen. They have different spatial and energy distributions near Mars. The two high-energy oxygen ion components, consisting of a high-energy “beam” and a high-energy “halo”, have different origins. (1) The high-energy (>∼100 eV) “beam” of O⁺ and ions are originating from the ionosphere. These ions form a highly asymmetric spatial distribution of escaping oxygen ions with respect to the direction of the convective electric field in the solar wind. (2) The high-energy (>∼100 eV) “halo” component contains O⁺ ions which are formed from the oxygen neutral exosphere by extreme ultraviolet radiation (EUV) and by charge exchange processes. These energetic halo ions can be found all around Mars. (3) The low energy O⁺ and ions (<∼100 eV) form a relatively symmetric spatial distribution around the Mars-Sun line. They originate from the ionosphere and from charge exchange processes between protons and exospheric oxygen atoms.The existence of the low- and the high-energy oxygen components is in agreement with recent in situ plasma measurements made by the ASPERA-3 instrument on the Mars Express mission. The analysis of the escaping oxygen ions suggests that the global energization of escaping planetary ions in the martian tail is controlled by the convective electric field.  相似文献   

New Horizons Alice ultraviolet observations of a stellar occultation by Jupiter’s atmosphere     

Thomas K. Greathouse  G.R. Gladstone  S.A. Stern  R.J. Vervack Jr.  M.H. Versteeg  L.A. Young  H. Throop 《Icarus》2010,208(1):293-305

The Alice ultraviolet spectrograph onboard the New Horizons spacecraft observed two occultations of the bright star χ Ophiucus by Jupiter’s atmosphere on February 22 and 23, 2007 during the approach phase of the Jupiter flyby. The ingress occultation probed the atmosphere at 32°N latitude near the dawn terminator, while egress probed 18°N latitude near the dusk terminator. A detailed analysis of both the ingress and egress occultations, including the effects of molecular hydrogen, methane, acetylene, ethylene, and ethane absorptions in the far ultraviolet (FUV), constrains the eddy diffusion coefficient at the homopause level to be  cm² s⁻¹, consistent with Voyager measurements and other analyses (Festou, M.C., Atreya, S.K., Donahue, T.M., Sandel, B.R., Shemansky, D.E., Broadfoot, A.L. [1981]. J. Geophys. Res. 86, 5717-5725; Vervack Jr., R.J., Sandel, B.R., Gladstone, G.R., McConnell, J.C., Parkinson, C.D. [1995]. Icarus 114, 163-173; Yelle, R.V., Young, L.A., Vervack Jr., R.J., Young, R., Pfister, L., Sandel, B.R. [1996]. J. Geophys. Res. 101 (E1), 2149-2162). However, the actual derived pressure level of the methane homopause for both occultations differs from that derived by [Festou et al., 1981] and [Yelle et al., 1996] from the Voyager ultraviolet occultations, suggesting possible changes in the strength of atmospheric mixing with time. We find that at 32°N latitude, the methane concentration is  cm⁻³ at 70,397 km, the methane concentration is  cm⁻³ at 70,383 km, the acetylene concentration is  cm⁻³ at 70,364 km, and the ethane concentration is  cm⁻³ at 70,360 km. At 18°N latitude, the methane concentration is  cm⁻³ at 71,345 km, the methane concentration is  cm⁻³ at 71,332 km, the acetylene concentration is cm⁻³ at 71,318 km, and the ethane concentration is  cm⁻³ at 71,315 km. We also find that the H₂ occultation light curve is best reproduced if the atmosphere remains cold in the microbar region such that the base of the thermosphere is located at a lower pressure level than that determined by in situ instruments aboard the Galileo probe (Seiff, A., Kirk, D.B., Knight, T.C.D., Young, R.E., Mihalov, J.D., Young, L.A., Milos, F.S., Schubert, G., Blanchard, R.C., Atkinson, D. [1998]. J. Geophys. Res. 103 (E10), 22857-22889) - the Sieff et al. temperature profile leads to too much absorption from H₂ at high altitudes. However, this result is highly model dependent and non-unique. The observations and analysis help constrain photochemical models of Jupiter’s atmosphere.  相似文献   

Impact of atmospheric gravity waves on the jovian ionosphere     

Daniel Barrow  Katia I. Matcheva 《Icarus》2011,211(1):609-622

We investigate the effects of atmospheric gravity waves on the vertical and horizontal structure of the ionosphere of Jupiter. The presented non-linear, two-dimensional model of the jovian ionosphere allows for spatially and temporally varying neutral wind and temperature fields and tracks the time evolution of six ionospheric species, , and . An analytical approach is used to validate the model results for linear, small-amplitude waves and to elucidate the mechanisms that leads to perturbations in the density of the main ion species, H⁺ and . We demonstrate that the long-lived H⁺ ions are perturbed directly by wave dynamics whereas short-lived ions such as are perturbed by chemical interactions with other perturbed ion species. The model is then applied using larger gravity wave amplitudes consistent with observations. Atmospheric gravity waves propagating at high altitudes create layers of enhanced electron density similar to the system of layers observed during the J0-ingress radio occultation of the Galileo spacecraft. Our best fit to the J0-ingress observation is achieved using an 82 min period forcing wave with horizontal and vertical wavelengths of 500 km and 60 km respectively, and peaks at 510 km above the 1 bar pressure level. We further investigate the effects of the wave-induced ion flux on the background ionospheric structure and demonstrate that in the presence of a gravity wave the background density profiles of the H⁺ and ions are significantly modified. We also find that the column density of has variations that can exceed 10% as the wave propagates.  相似文献   

Extension of atmospheric dust loading to high altitudes during the 2001 Mars dust storm: MGS TES limb observations     

R. Todd Clancy  Michael J. Wolff  Bruce A. Cantor  Timothy H. McConnochie 《Icarus》2010,207(1):98-109

Mars Global Surveyor (MGS) visible (solarband bolometer) and thermal infrared (IR) spectral limb observations from the Thermal Emission Spectrometer (TES) support quantitative profile retrievals for dust opacity and particle sizes during the 2001 global dust event on Mars. The current analysis considers the behavior of dust lifted to altitudes above 30 km during the course of this storm; in terms of dust vertical mixing, particle sizes, and global distribution. TES global maps of visible (solarband) limb brightness at 60 km altitude indicate a global-scale, seasonally evolving (over 190-240° solar longitudes, L_S) longitudinal corridor of vertically extended dust loading (which may be associated with a retrograde propagating, wavenumber 1 Rossby wave). Spherical radiative transfer analysis of selected limb profiles for TES visible and thermal IR radiances provide quantitative vertical profiles of dust opacity, indicating regional conditions of altitude-increasing dust mixing ratios. Observed infrared spectral dependences and visible-to-infrared opacity ratios of dust scattering over 30-60 km altitudes indicate particle sizes characteristic of lower altitudes (cross-section weighted effective radius, ), during conditions of significant dust transport to these altitudes. Conditions of reduced dust loading at 30-60 km altitudes present smaller dust particle sizes . These observations suggest rapid meridional transport at 30-80 km altitudes, with substantial longitudinal variation, of dust lifted to these altitudes over southern hemisphere atmospheric regions characterized by extraordinary (m/s) vertical advection velocities. By L_S=230° dust loading above 50 km altitudes decreased markedly at southern latitudes, with a high altitude (60-80 km) haze of fine (likely) water ice particles appearing over 10°S-40°N latitudes.  相似文献   

Rosetta-Alice observations of exospheric hydrogen and oxygen on Mars     

Paul D. Feldman  Andrew J. Steffl  Michael F. A’Hearn  S. Alan Stern  David C. Slater  Henry B. Throop  Lori M. Feaga 《Icarus》2011,214(2):394-399

The European Space Agency’s Rosetta spacecraft, en route to a 2014 encounter with comet 67P/Churyumov-Gerasimenko, made a gravity assist swing-by of Mars on 25 February 2007, closest approach being at 01:54 UT. The Alice instrument on board Rosetta, a lightweight far-ultraviolet imaging spectrograph optimized for in situ cometary spectroscopy in the 750-2000 Å spectral band, was used to study the daytime Mars upper atmosphere including emissions from exospheric hydrogen and oxygen. Offset pointing, obtained five hours before closest approach, enabled us to detect and map the H i Lyman-α and Lyman-β emissions from exospheric hydrogen out beyond 30,000 km from the planet’s center. These data are fit with a Chamberlain exospheric model from which we derive the hydrogen density at the 200 km exobase and the H escape flux. The results are comparable to those found from the Ultraviolet Spectrometer experiment on the Mariner 6 and 7 fly-bys of Mars in 1969. Atomic oxygen emission at 1304 Å is detected at altitudes of 400-1000 km above the limb during limb scans shortly after closest approach. However, the derived oxygen scale height is not consistent with recent models of oxygen escape based on the production of suprathermal oxygen atoms by the dissociative recombination of .  相似文献   

An investigation of Pluto’s troposphere using stellar occultation light curves and an atmospheric radiative-conductive-convective model     

Angela M. Zalucha  Xun Zhu  Darrell F. Strobel  J.L. Elliot 《Icarus》2011,214(2):685-700

We use a radiative-conductive-convective model to assess the height of Pluto’s troposphere, as well as surface pressure and surface radius, from stellar occultation data from the years 1988, 2002, and 2006. The height of the troposphere, if it exists, is less than 1 km for all years analyzed. Pluto has at most a planetary boundary layer and not a troposphere. As in previous analyses of Pluto occultation light curves, we find that the surface pressure is increasing with time, assuming that latitude and longitude variations in Pluto’s atmosphere are negligible. The surface pressure is found to be slightly higher ( μbar in 1988,  μbar in 2002, and 18.5 ± 4.7 μbar in 2006) than in our previous analyses with the troposphere excluded. The surface radius is determined to be . Comparison of the minimum reduced chi-squared values between the best-fit radiative-conductive-convective (i.e., troposphere-included) model and best-fit radiative-conductive (i.e., troposphere-excluded) shows that the troposphere-included model is only a slightly better fit to the data for all 3 years. Uncertainties in the small-scale physical processes of Pluto’s lower atmosphere and consequently the functional form of the model troposphere lend more confidence to the troposphere-excluded results.  相似文献   

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris   总被引：2，自引：0，他引：2  

Alex S. Konopliv  Charles F. Yoder  E. Myles Standish  Dah-Ning Yuan  William L. Sjogren 《Icarus》2006,182(1):23-50

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.  相似文献   

The Planetary Fourier Spectrometer (PFS) onboard the European Mars Express mission     

V. Formisano  F. Angrilli  S. Atreya  D. Biondi  M.I. Blecka  L. Colangeli  F. Esposito  M. Giuranna  V. Gnedykh  G. Hansen  I. Khatuntsev  N. Ignatiev  A. Jurewicz  J. Lopez Moreno  A. Mattana  E. Mencarelli  V. Moroz  F. Nespoli  R. Orfei  V. Orofino  D. Patsaev  M. Rataj  J. Rodriguez  B. Saggin  L. Zasova 《Planetary and Space Science》2005,53(10):963-974

The Planetary Fourier Spectrometer (PFS) for the Mars Express mission is an infrared spectrometer optimised for atmospheric studies. This instrument has a short wave (SW) channel that covers the spectral range from 1700 to (1.2-) and a long-wave (LW) channel that covers 250- (5.5-). Both channels have a uniform spectral resolution of . The instrument field of view FOV is about 1.6^° (FWHM) for the Short Wavelength channel (SW) and 2.8^° (FWHM) for the Long Wavelength channel (LW) which corresponds to a spatial resolution of 7 and 12 km when Mars is observed from an height of 250  km. PFS can provide unique data necessary to improve our knowledge not only of the atmosphere properties but also about mineralogical composition of the surface and the surface-atmosphere interaction.The SW channel uses a PbSe detector cooled to 200-220 K while the LW channel is based on a pyroelectric (L_iT_aO₃) detector working at room temperature. The intensity of the interferogram is measured every 150 nm of physical mirrors displacement, corresponding to 600 nm optical path difference, by using a laser diode monochromatic light interferogram (a sine wave), whose zero crossings control the double pendulum motion. PFS works primarily around the pericentre of the orbit, only occasionally observing Mars from large distances. Each measurements take 4 s, with a repetition time of 8.5 s. By working roughly 0.6 h around pericentre, a total of 330 measurements per orbit will be acquired 270 looking at Mars and 60 for calibrations. PFS is able to take measurements at all local times, facilitating the retrieval of surface temperatures and atmospheric vertical temperature profiles on both the day and the night side.  相似文献   

Low-speed impacts between rubble piles modeled as collections of polyhedra, 2     

D.G. Korycansky  Erik Asphaug 《Icarus》2009,204(1):316-329

We present the results of additional calculations involving the collisions of km-scale rubble piles. In new work, we used the Open Dynamics Engine (ODE), an open-source library for the simulation of rigid-body dynamics that incorporates a sophisticated collision-detection and resolution routine. We found that using ODE resulted in a speed-up of approximately a factor of 30 compared with previous code. In this paper we report on the results of almost 1200 separate runs, the bulk of which were carried out with 1000-2000 elements. We carried out calculations with three different combinations of the coefficients of friction η and (normal) restitution ?: low (η=0,?=0.8), medium (η=0,?=0.5), and high (η=0.5,?=0.5) dissipation.For target objects of ∼1 km in radius, we found reduced critical disruption energy values in head-on collisions from 2 to 100 J kg⁻¹ depending on dissipation and impactor/target mass ratio. Monodisperse objects disrupted somewhat more easily than power-law objects in general. For oblique collisions of equal-mass objects, mildly off-center collisions (b/b₀=0.5) seemed to be as efficient or possibly more efficient at collisional disruption as head-on collisions. More oblique collisions were less efficient and the most oblique collisions we tried (b/b₀=0.866) required up to ∼200 J kg⁻¹ for high-dissipation power-law objects. For calculations with smaller numbers of elements (total impactor or 200 elements) we found that collisions were more efficient for smaller numbers of more massive elements, with values as low as for low-dissipation cases. We also analyzed our results in terms of the relations proposed by Stewart and Leinhardt [Stewart, S.T., Leinhardt, Z.M., 2009. Astrophys. J. 691, L133-L137] where where Q_R is the impact kinetic energy per unit total mass m_i+m_T. Although there is a significant amount of scatter, our results generally bear out the suggested relation.  相似文献   

Orbital effects of non-isotropic mass depletion of the atmospheres of evaporating hot Jupiters in extrasolar systems     

Lorenzo Iorio 《New Astronomy》2012,17(3):356-361

  相似文献   

Photochemical escape of oxygen from Mars: A comparison of the exobase approximation to a Monte Carlo method     

Jane L. Fox  Aleksander B. Ha? 《Icarus》2009,204(2):527-544

The non-thermal escape of neutral O atoms from Mars at the current epoch is largely due to dissociative recombination of :