Mapping of the cusp plasma precipitation on the surface of Mercury

The presence of a magnetosphere around Mercury plays a fundamental role on the way the solar wind plasma interacts with the planet. Since the observations suggest that Mercury should occupy a large fraction of its magnetosphere and because of lack of an atmosphere, significant differences in solar wind-magnetosphere coupling are expected to exist with respect to the Earth case. On the basis of a modified Tsyganenko T96 model we describe the geometry of the magnetic field that could characterize Mercury, and its response to the variations of the impinging solar wind and of the interplanetary magnetic field. The investigation is focused on the shape and dimension of the open magnetic field regions (cusps) that allow the direct penetration of magnetosheath plasma through the exosphere of Mercury, down to its surface. The precipitating particle flux and energy are evaluated as a function of the open field line position, according to different solar wind conditions. A target of this study is the evaluation of the sputtered particles from the crust of the planet, and their contribution to the exospheric neutral particle populations. Such estimates are valuable in the frame of a neutral particle analyser to be proposed on board of the ESA/BepiColombo mission.     

Solar wind-magnetosphere interaction as simulated by a 3-D em particle code

We have studied the solar wind-magnetosphere interaction using a 3-D electromagnetic particle code. The results for an unmagnetized solar wind plasma streaming past a dipole magnetic field show the formation of a magnetopause and a magnetotail, the penetration of energetic particles into cusps and radiation belt and dawn-dusk asymmetries. The effects of interplanetary magnetic field (IMF) have been investigated in a similar way as done by MHD simulations. The simulation results with a southward IMF show the shrunk magnetosphere with great particle entry into the cusps and nightside magnetosphere. This is a signature of a magnetic reconnection at the dayside magnetopause. After a quasi-stable state is established with an unmagnetized solar wind we switched on a solar wind with an northward IMF. In this case the significant changes take place in the magnetotail. The waving motion was seen in the magnetotail and its length was shortened. This phenomena are consistent with the reconnections which occur at the high latitude magnetopause. In our simulations kinetic effects will determine the self-consistent anomalous resistivity in the magnetopause that causes reconnections.Deceased January 24, 1993; R. Bunemanet al. 1993.     

太阳圆面活动体光谱分析方法的比较

在研究暗条、耀斑等日面活动体的光谱时，观测谱线上都叠加上一定强度的背景光影响。如何消去这种影响对分析结果影响较大。目前，经典的方法有两大类：第一类是减法消去法，这种方法广泛为国内学者采用；第二类方法是除法消去法，即国外学者提出的（微）方模型方法。本文通过误差传递方法，比较了上述二大类方法，结果发现，所有方法的误差传递形式都是一样的，但（微）方模型的背景光误差可能要小些，至少从理论上讲要处理得合理些，因此（微）方模型要优越些。然而，从根本上讲，可以说，到目前为止，太阳圆面背景光处理问题还没有彻底解决。这方面的工作有待进一步探索。     

Strong interplanetary field enhancements at Ulysses—evidence of dust trails' interaction with the solar wind?

Interplanetary field enhancements were first discovered in the vicinity of Venus. These events are characterised by an increase in the magnitude of the heliospheric magnetic field with a near-symmetrical, sometimes thorn-shaped profile, and last from minutes to hours. Surveys of the events near Venus and Earth indicated clustering of the events in inertial space, which suggested that their sources were Solar System objects other than the Sun. A survey is presented of strong events of this type detected by the Ulysses spacecraft from 1990 to late 2001. Most of the events are accompanied by a discontinuity in the field direction near the events' centres. Other discontinuities are often symmetrical about the enhancement. The majority of events last less than two hours. When examined as a whole, the events tend to be accompanied by subtle changes in some plasma parameters. The majority of the enhancements are accompanied by magnetic holes on their fringes. The enhancements' occurrence rate increases with decreasing heliocentric distance. Possible formation mechanisms are discussed. No link was found with solar, or solar wind sources. Several aspects of the survey results are consistent with an origin related to cometary dust trails. Possible processes associated with a dust-solar wind interaction are discussed.     

On the propagation of hydromagnetic waves in a plasma of thermal and suprathermal components

The propagation of MHD waves is studied when two ideal fluids, thermal and suprathermal gases, coupled by magnetic field are moving with the steady flow velocity. The fluids move independently in a direction perpendicular to the magnetic field but gets coupled along the field. Due to the presence of flow in suprathermal and thermal fluids there appears forward and backward waves. All the forward and backward modes propagate in such a way that their rate of change of phase speed with the thermal Mach number is same. It is also found that besides the usual hydromagnetic modes there appears a suprathermal mode which propagates with faster speed. Surface waves are also examined on an interface formed with composite plasma (suprathermal and thermal gases) on one side and the other is a non-magnetized plasma. In this case, the modes obtained are two or three depending on whether the sound velocity in thermal gas is equal to or greater than the sound velocity in suprathermal gas. The results lead to the conclusion that the interaction of thermal and suprathermal components may lead to the occurrence of an additional mode called suprathermal mode whose phase velocity is higher than all the other modes.     

High kinetic energy density jets in the Earth’s magnetosheath: A case study

We study several high kinetic energy density jets observed during a traversal of the dayside magnetosheath by the Cluster spacecraft on March 17, 2001, at various distances from the magnetopause, generally characterised by anomalously high values of the local magnetosonic Mach number. We concentrate on two jets observed just outside the magnetopause, the first almost parallel to the GSM x axis and the second directed northward-tailward along the nominal magnetopause surface. We present evidence that none of them can be ascribed to magnetic reconnection at the magnetopause and show that the magnetopause is severely deformed by the jets, so that its local normal forms an angle of 97° with the quiet time magnetopause normal. On these grounds, we suggest that the indentation of the magnetopause is caused by an anti-sunward jet ramming into the magnetopause slightly equatorward of the northern cusp and that the northward-tailward jet is the result of its reflection at the deformed magnetopause. Finally, we briefly discuss our results by comparing them with past studies of events which in some way recall the one analysed herein.     

Radar absorption due to a corotating interaction region encounter with Mars detected by MARSIS

Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) is a subsurface and topside ionosphere radar sounder aboard the European Space Agency spacecraft Mars Express, in orbit at Mars since 25 December 2003, and in operation since 17 June 2005. The ionospheric sounding mode of MARSIS is capable of detecting the reflection of the sounding wave from the martian surface. This ability has been used in previous work to show that the surface reflection is absorbed and disappears during periods when high fluxes of energetic particles are incident on the ionosphere of Mars. These absorption events are believed to be the result of increased collisional damping of the sounding wave, caused by increased electron density below the spacecraft, in turn caused by impact ionization from the impinging particles. In this work we identify two absorption events that were isolated during periods when the surface reflection is consistently visible and when Mars is nearly at opposition. The visibility of the surface reflection is viewed in conjunction with particle and photon measurements taken at both Mars and Earth. Both absorption events are found to coincide with Earth passing through solar wind speed and ion flux signatures indicative of a corotating interaction region (CIR). The two events are separated by an interval of approximately 27 days, corresponding to one solar rotation. The first of the two events coincides with abruptly enhanced particle fluxes seen in situ at Mars. Simultaneous with the particle enhancement there are an abrupt decrease in the intensity of electron oscillations, typically seen by the Mars Express particle instrument ASPERA-3 between the magnetic pileup boundary and the martian bow shock, and a sharp drop in the solar wind pressure, seen in the proxy quantity based on MGS magnetometer observations. The decrease in oscillation intensity is therefore the probable effect of a relaxation of the martian bow shock. The second absorption event does not show a particle enhancement and complete ASPERA-3 data during that time are unavailable. Other absorption events are the apparent result of solar X-ray and XUV enhancements. We conclude that surface reflection absorption events are sometimes caused by enhanced ionospheric ionization from high energy particles accelerated by the shocks associated with a CIR. A full statistical analysis of CIRs in relation to observed absorption events in conjunction with a quantitative analysis of the deposition of ionization during space weather events is needed for a complete understanding of this phenomenon. If such analyses can be carried out, radar sensing of the martian ionosphere might be useful as a space weather probe.     

Solar wind variation with the cycle

The cyclic evolution of the heliospheric plasma parameters is related to the time-dependent boundary conditions in the solar corona. “Minimal” coronal configurations correspond to the regular appearance of the tenuous, but hot and fast plasma streams from the large polar coronal holes. The denser, but cooler and slower solar wind is adjacent to coronal streamers. Irregular dynamic manifestations are present in the corona and the solar wind everywhere and always. They follow the solar activity cycle rather well. Because of this, the direct and indirect solar wind measurements demonstrate clear variations in space and time according to the minimal, intermediate and maximal conditions of the cycles. The average solar wind density, velocity and temperature measured at the Earth’s orbit show specific decadal variations and trends, which are of the order of the first tens per cent during the last three solar cycles. Statistical, spectral and correlation characteristics of the solar wind are reviewed with the emphasis on the cycles.     

Location of the bow shock and ion composition boundaries at Venus—initial determinations from Venus Express ASPERA-4

For the first time since 1992 when the Pioneer Venus Orbiter (PVO) ceased to operate, there is again a plasma instrument in orbit around Venus, namely the ASPERA-4 flown on Venus Express (inserted into an elliptical polar orbit about the planet on April 11, 2006). In this paper we report on measurements made by the ion and electron sensors of ASPERA-4 during their first five months of operation and, thereby, determine the locations of both the Venus bow shock (BS) and the ion composition boundary (ICB) under solar minimum conditions. In contrast to previous studies based on PVO data, we employ a 3-parameter fit to achieve a realistic shape for the BS. We use a different technique to fit the ICB because this latter boundary cannot be represented by a conic section. Additionally we investigate the dependence of the location of the BS on solar wind ram pressure (based on ASPERA-4 solar wind data) and solar EUV flux (using a proxy from Earth).     

The sodium tail of the Moon

During the few days centered about new Moon, the lunar surface is optically hidden from Earth-based observers. However, the Moon still offers an observable: an extended sodium tail. The lunar sodium tail is the escaping “hot” component of a coma-like exosphere of sodium generated by photon-stimulated desorption, solar wind sputtering and meteoroid impact. Neutral sodium atoms escaping lunar gravity experience solar radiation pressure that drives them into the anti-solar direction forming a comet-like tail. During new Moon time, the geometry of the Sun, Moon and Earth is such that the anti-sunward sodium flux is perturbed by the terrestrial gravitational field resulting in its focusing into a dense core that extends beyond the Earth. An all-sky camera situated at the El Leoncito Observatory (CASLEO) in Argentina has been successfully imaging this tail through a sodium filter at each lunation since April 2006. This paper reports on the results of the brightness of the lunar sodium tail spanning 31 lunations between April 2006 and September 2008. Brightness variability trends are compared with both sporadic and shower meteor activity, solar wind proton energy flux and solar near ultra violet (NUV) patterns for possible correlations. Results suggest minimal variability in the brightness of the observed lunar sodium tail, generally uncorrelated with any single source, yet consistent with a multi-year period of minimal solar activity and non-intense meteoric fluxes.     

Comparison of X-rays from Comets LINEAR (C/1999 S4) and McNaught-Hartley (C/1999 T1)

The Chandra X-ray Observatory (CXO) observations of Comets McNaught-Hartley (MH) and LINEAR S4 (S4) have been processed in the same way to compare X-rays from those comets. The X-ray isophotes are crescent-like in S4 and more circular in MH because of the different phase angles (98° and 44°, respectively). The peak X-ray brightness is greater in S4 than that in MH by a factor of 1.5 and smaller by a factor of 1.7 after the correction for heliocentric distance. The X-ray luminosities of MH and S4 are equal to 8.6 and 1.4×10¹⁵ erg s⁻¹ inside the apertures of ρ=1.5 and 0.5×10⁴ km, respectively. (Brightness is 20% of the peak value at these ρ.) Efficiencies of X-ray excitation corrected to the solar wind flow are similar and equal to 4.3×10⁻¹⁴ erg AU^3/2 in both comets. This confirms the solar wind excitation of X-rays in comets. Spectra of the comets were extracted with a special care of the background correction and using an energy-dependent spectral resolution code. The MH spectrum consists of ten emissions instead of nine emissions in the previously published spectrum. The new emission at 307 eV fills in a strong minimum in the previous spectrum and removes the major difference between that spectrum and the synthetic spectrum. This emission is assigned to the C⁺⁴ and Mg⁺⁹ lines. The positions of the other emissions and their identification are similar to those in the previous spectrum. The S4 spectrum consists of eight emissions, and four emissions are the same as in MH. The line identification is given. Ion ratios in the solar wind have been extracted from the spectra. O⁺⁸/O⁺⁷ is equal to 0.29±0.04 and 0.14±0.02 in MH and S4, and this difference correlates with the higher solar wind speed in S4. Ne⁺⁹/O⁺⁷ is (15±6)×10⁻³ and (19±7)×10⁻³, and these are the first data on Ne⁺⁹ in the solar wind. C⁺⁶/O⁺⁷ is 0.7±0.2 in both MH and S4. X-ray spectroscopy of comets may be used as a diagnostic tool to study the solar wind composition.     

On the effect of the martian crustal magnetic field on atmospheric erosion

Without the shielding of a strong intrinsic magnetic field, the martian atmosphere directly interacts with the impacting solar wind. The neutral constituents of the atmospheric corona can be ionized, and then picked up and accelerated by the magnetic field and convection electric field in the solar wind. A significant fraction of pickup ions escape Mars’ gravitational pull and are lost to space. This non-thermal escape process of heavy species is an important mechanism responsible for atmospheric erosion. While there is a perception that the martian magnetic anomalies are significant for the ionospheric density distribution and the bow shock standoff location, little is known about the quantitative influence of the martian crustal magnetic field on the global distribution of escaping pickup ions. In this paper, we apply a newly developed Monte Carlo ion transport model to resolve the crustal field effect on the pickup oxygen ion distribution around Mars. The background magnetic and electric fields, in which test particles are followed, are calculated using an independent three-dimensional multispecies MHD model. The effects of the crustal magnetic field on particle escape are quantified by varying the crustal field orientation in the model setup and comparing the corresponding test particle simulation results. The comparison is made by turning on or off the crustal field or changing the local time of the strongest field from the dayside to the dawnside. It is found that without the protection of the crustal magnetic field, the total amount of atmospheric escape through the tail region would be enhanced by more than a factor of two. It is shown that the crustal magnetic field not only regionally deflects the solar wind around the martian atmosphere, but also has an important global effect on atmospheric erosion and thus on long-term atmospheric evolution.     

Ion escape from Mars as a function of solar wind conditions: A statistical study

The influence of solar EUV and solar wind conditions on ion escape at Mars is investigated using ion data from the Aspera-3 instrument on Mars Express, combined with solar wind proxy data obtained from the Mars Global Surveyor (MGS) spacecraft. A solar EUV flux proxy based on data from the Earth position, scaled and shifted in time for Mars, is used to study relatively long time scale changes related to solar EUV variability. Data from May 2004 until November 2005 has been used. A clear dependence on the strength of the subsolar magnetic field as inferred from MGS measurements is seen in the ion data. The region of significant heavy ion flows is compressed and the heavy ion flux density is higher for high subsolar magnetic field strength. Because of the difference in outflow area, the difference in estimated total outflow is somewhat less than the difference in average flux density. We confirm previous findings that escaping planetary ions are mainly seen in the hemisphere into which the solar wind electric field is pointed. The effect is more pronounced for the high subsolar magnetic field case.The average ion motion has a consistent bias towards the direction of the solar wind electric field, but the main motion is in the antisunward direction. The antisunward flow velocity increases with tailward distance, reaching above at 2 to 3 martian radii downtail from Mars for O⁺ ions. Different ion species reach approximately the same bulk flow energy. We did not find any clear correlation between the solar EUV flux and the ion escape distribution or rate, probably because the variation of the solar EUV flux over our study interval was too small. The results indicate that the solar wind and its magnetic field directly interacts with the ionosphere of Mars, removing more ions for high subsolar magnetic field strength. The interaction region and the tail heavy ion flow region are not perfectly shielded from the solar wind electric field, which accelerates particles over relatively large tail distances.     

A plasma flow velocity boundary at Mars from the disappearance of electron plasma oscillations

The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) on the Mars Express (MEX) spacecraft is capable of measuring ionospheric electron density by the use of two main methods: remote radar sounding and from the excitation of local plasma oscillations. The frequency of the locally excited electron plasma oscillations is used to measure the local electron density. However, plasma oscillations are not observed when the plasma flow velocity is higher than about 160 km/s, which occurs mainly in the solar wind and magnetosheath. As a consequence, in many passes, there is a sudden disappearance of the plasma oscillations as the spacecraft enters into the magnetosheath. This fact allows us to identify a flow velocity boundary on the dayside, between the ionosphere of Mars and the shocked solar wind. This paper summarizes the results of the measurement of 552 orbits mostly over a period from August 4, 2005 to August 17, 2007. The boundary points found using MARSIS have been verified by measurements from the Analyzer of Space Plasma and Energetic Atoms (ASPERA-3) Electron Spectrometer (ELS) instrument on Mars Express. The average position of the flow velocity boundary is compared to flow velocity simulations computed using hybrid model and other boundaries. The boundary altitude is slightly lower than the magnetic pile-up boundary determined using Phobos 2 and Mars Global Surveyor (MGS) crossings, but it is in good agreement with the induced magnetospheric boundary determined by ASPERA-3. Investigation of the effect of the crustal magnetic field revealed that the flow velocity boundary is raised at the locations with strong crustal magnetic fields.     

The 27-Day Signal in Sunspot Number Series and the Solar Dynamo

We analyze the Wolf number daily series WN (1849 to present) as well as two other related series characterizing solar activity. Our analysis consists in computing the amplitude of a given Fourier component in a sliding time window and examining its long-term evolution. We start with the well-known 27.03- and 27.6-day periods and observe strong decadal variations of this amplitude as well as a sharp increase of the average value starting around 1905. We then consider a packet of 31 lines with periods from 25.743 to 28.453 days, which is shown to be a better representation of the synodic solar rotation. We first examine the temporal evolution of individual lines, then the energy of the packet. The energy of the packet increases sharply at the beginning of the 20th century, leading by more than two decades the well-known increase of the Wolf number. The nonaxisymmetry of sunspots increases before the total increase of activity and may be considered as a precursor. We discuss briefly and tentatively this observation in terms of solar dynamo theory.     

Mercury’s magnetosphere-solar wind interaction for northward and southward interplanetary magnetic field: Hybrid simulation results

Analysis of global hybrid simulations of Mercury’s magnetosphere-solar wind interaction is presented for northward and southward interplanetary magnetic field (IMF) orientations in the context of MESSENGER’s first two encounters with Mercury. The global kinetic simulations reveal the basic structure of this interaction, including a bow shock, ion foreshock, magnetosheath, cusp regions, magnetopause, and a closed ion ring belt formed around the planet within the magnetosphere. The two different IMF orientations induce different locations of ion foreshock and different magnetospheric properties: the dayside magnetosphere is smaller and cusps are at lower latitudes for southward IMF compared to northward IMF whereas for southward IMF the nightside magnetosphere is larger and exhibits a thin current sheet with signatures of magnetic reconnection and plasmoid formation. For the two IMF orientations the ion foreshock and quasi-parallel magnetosheath manifest ion-beam-driven large-amplitude oscillations, whereas the quasi-perpendicular magnetosheath shows ion-temperature-anisotropy-driven wave activity. The ions in Mercury’s belt remain quasi-trapped for a limited time before they are either absorbed by Mercury’s surface or escape from the magnetosphere. The simulation results are compared with MESSENGER’s observations.     

冕洞演化与太阳活动周的关系及在日地预报中的应用

冕洞的研究在近二十多年里取得了丰硕的成果。本文回顾了冕洞的发现及观测历史，系统阐述了冕洞的结构特征、形成及演化规律，讨论了冕洞对日地空间产生的影响，冕洞与超级活动区的关系以及冕洞在太阳活动预报中所起的作用，在此基础上利用１９７０—１９９５年的冕洞资料对冕洞的时空分布和磁极性演化规律与太阳活动周的关系，以及冕洞与太阳风速度、地磁扰动等方面进行分析研究，得出以下结论：（１）冕洞在南北半球的分布在形态上基本是对称的，但在冕洞数量上北半球稍占优势；（２）冕洞的盛衰演化呈周期性，表现为赤道冕洞周期与黑子周期是完全一致的，极冕洞周期与黑子周期相位相差１８０°；（３）赤道冕洞的纬度分布随太阳活动周上升而上升，当太阳活动周达到极大值时，它也达到极大，然后再随太阳活动周下降而下降，极冕洞的纬度延伸方向演化与赤道冕洞相反；（４）极冕洞的极场呈１１年周期性，并且极场反转出现在太阳活动峰年期间；（５）太阳风和地磁扰动与冕洞的演化有着密切的关系     

Periodic Appearance of Coronal Holes and the Related Variation of Solar Wind Parameters

We compared the variability of coronal hole (CH) areas (determined from daily GOES/SXI images) with solar wind (daily ACE data) and geomagnetic parameters for the time span 25 January 2005 until 11 September 2005 (late declining phase of solar cycle 23). Applying wavelet spectral analysis, a clear 9-day period is found in the CH time series. The GOES/SXI image sequence suggests that this periodic variation is caused by a mutual triangular distribution of CHs ∼120° apart in longitude. From solar wind parameters a 9-day periodicity was obtained as well, simultaneously with the 9-day period in the CH area time series. These findings provide strong evidence that the 9-day period in solar wind parameters, showing up as higher harmonic of the solar rotation frequency, is caused by the “periodic” longitudinal distribution of CHs on the Sun recurring for several solar rotations. The shape of the wavelet spectrum from the Dst index matches only weakly with that from the CH areas and is more similar to the wavelet spectrum of the solar wind magnetic field magnitude. The distinct 9-day period does not show up in sunspot group areas which gives further evidence that the solar wind modulation is strongly related to CH areas but not to active region complexes. The wavelet power spectra for the whole ACE data range (∼1998 – 2006) suggest that the 9-day period is not a singular phenomenon occurring only during a specific time range close to solar minimum but is occasionally also present during the maximum and decay phase of solar cycle 23. The main periods correspond to the solar rotation (27^d) as well as to the second (13.5^d) and third (9^d) harmonic. Electronic Supplementary Material The online version of this article () contains supplementary material, which is available to authorized users.     

Quasi-Periodic Releases of Streamer Blobs and Velocity Variability of the Slow Solar Wind near the Sun

We search for persistent and quasi-periodic release events of streamer blobs during 2007 with the Large Angle Spectrometric Coronagraph on the Solar and Heliospheric Observatory and assess the velocity of the slow solar wind along the plasma sheet above the corresponding streamer by measuring the dynamic parameters of blobs. We find ten quasi-periodic release events of streamer blobs lasting for three to four days. In each day of these events, we observe three – five blobs. The results are in line with previous studies using data observed near the last solar minimum. Using the measured blob velocity as a proxy for that of the mean flow, we suggest that the velocity of the background slow solar wind near the Sun can vary significantly within a few hours. This provides an observational manifestation of the large velocity variability of the slow solar wind near the Sun.     

Parametric analysis of modeled ion escape from Mars

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.