SGP4/SDP4模型精度分析

本文基于最新发布的SGP4/SDP4(Simplified General Perturbation Version 4/Simplified Deep-space Perturbation Version 4)模型设计了一套定轨方案,从空间目标库中挑选出不同类型和轨道参数的1120个目标进行计算,定量给出了SGP4/SDP4模型处理不同类型空间目标的定轨预报精度.结果表明:近地目标定轨精度为百米量级;半同步和同步轨道定轨精度平均为0.7和1.9km.椭圆轨道目标的定轨精度与偏心率有关,除少数e>0.8的椭圆轨道目标,绝大多数椭圆轨道目标定轨误差均小于10km.用SGP4/SDP4模型对近地目标预报3天,半同步轨道预报30天,同步轨道预报15天,椭圆轨道预报1天,预报误差一般不超过40km.     

Study on the Orbit Prediction Errors of Space Objects Based on Historical TLE Data

Two line element (TLE) released by the North American Aerospace Defense Command (NORAD) is widely used by aerospace workers, and the matched SGP4/SDP4 (Simplified General Perturbation Version 4/Simplified Deep-space Perturbation Version 4) model is used to propagate it. Nevertheless, no corresponding information about its accuracy and covariance is clearly given, thus the application of the TLE data is greatly restricted. In this paper, the determined and predicted orbits are compared to generate the orbit error data, based on the historical TLE data obtained from the Space-Track website and the SGP4/SDP4 model. By dividing different time bins, the fitting coefficients of the variation of orbit prediction error with time are given for each space object, and the characteristics of the error evolution are further discussed for the different types of orbits. The mean analytic model of the orbit prediction error evolution with time is given respectively for the four orbit types of space objects, which provides a valuable reference for extending the application of the TLE data.     

基于历史TLE的空间目标轨道预报误差演化规律研究

北美防空司令部(North American Aerospace Defense Command, NORAD)发布的双行根数(Two Line Element, TLE)是广大航天工作者最常用的轨道根数,与其对应的轨道模型是SGP4/SDP4 (Simplified General Perturbation Version 4/Simplified Deep-space Perturbation Version 4)解析模型.由于TLE中并没有包含相应的轨道精度信息,编目轨道的应用范围受到很大的限制.基于Space-Track网站发布的历史TLE数据和配套的SGP4/SDP4动力学模型,采用定轨标预报的方法统计并生成了大量目标轨道的预报误差,通过对预报轨道的时间区间划分给出了每个目标的预报误差随预报时间变化的拟合系数,并进一步对不同类型轨道预报误差的演化规律和特征进行了分类讨论,给出了4种轨道类型目标的轨道预报误差随时间演化的平均解析模型,为拓展双行根数的应用提供有价值的参考.     

Venus Monitoring Camera for Venus Express

The Venus Express mission will focus on a global investigation of the Venus atmosphere and plasma environment, while additionally measuring some surface properties from orbit. The instruments PFS and SPICAV inherited from the Mars Express mission and VIRTIS from Rosetta form a powerful spectrometric and spectro-imaging payload suite. Venus Monitoring Camera (VMC)—a miniature wide-angle camera with 17.5° field of view—was specifically designed and built to complement these experiments and provide imaging context for the whole mission. VMC will take images of Venus in four narrow band filters (365, 513, 965, and 1000 nm) all sharing one CCD. Spatial resolution on the cloud tops will range from 0.2 km/px at pericentre to 45 km/px at apocentre when the full Venus disc will be in the field of view. VMC will fulfill the following science goals: (1) study of the distribution and nature of the unknown UV absorber; (2) determination of the wind field at the cloud tops (70 km) by tracking the UV features; (3) thermal mapping of the surface in the 1 μm transparency “window” on the night side; (4) determination of the global wind field in the main cloud deck (50 km) by tracking near-IR features; (5) study of the lapse rate and H₂O content in the lower 6–10 km; (6) mapping O₂ night-glow and its variability.     

Cassini RPWS observations of dust in Saturn's E Ring

The Cassini radio and plasma wave science (RPWS) instrument is sensitive to few-micron dust grains impacting on the spacecraft at relative speeds of order 10 km/s. Through the first year or so of operations in orbit at Saturn, the RPWS has made a number of both inclined and equatorial crossings of the E ring, particularly near the orbit of Enceladus. Assuming water ice grains, the typical size particle detected by the RPWS has a radius of a few microns. Peak impact rates of about 50 s⁻¹ are found near the orbit of Enceladus corresponding to densities of order 5×10⁻⁴ m⁻³. The variation of dust fluxes as a function of height above or below the equator is well described by a Gaussian distribution with a scale height of about 2800 km although there is usually some non-Gaussian variation near the peak fluxes suggesting some structure in the core of the ring. Offsets of the peak number densities are typically of the order of a few hundred km from the geometric equator. A near-equatorial radial profile through the orbit of Enceladus shows a sharply peaked distribution at the orbit of the moon. A size distribution averaged over several passes through the orbit of Enceladus is determined which varies as m^−2.80. The peak in dust number density at the orbit of Enceladus is consistent with previous optical measurements and strongly supports the suggestion that Enceladus is a primary source for E ring particles.     

日月轨道计算对TLE深空编目目标预报精度的影响

双行根数(Two Line Element, TLE)作为一类广泛使用的空间物体编目数据, 其预报精度和误差特性是TLE编目 在空间碎片研究中所要关注的问题之一. TLE编目需要配合SGP4/SDP4 (Simplified General Perturbations 4/Simplified Deep Space 4)模型进行轨道预报, 对深空物体来说, 主要考虑带谐项$J_2$、$J_3$、$J_4$摄动、 第三体日月摄动和特殊轨道共振问题修正等. 其中, SGP4/SDP4模型第三体摄动计算时, 对日月轨道近似采用了长期进动根数和 简单平运动的方式, 在外推10d时存在约2$^\circ$--3${^\circ     

The outcome of the protoplanetary disk of very massive stars

We suggest that planets, brown dwarfs, and even low mass stars can be formed by fragmentation of protoplanetary disks around very massive stars (M ? 100 M_⊙). We discuss how fragmentation conditions make the formation of very massive planetary systems around very massive stars favorable. Such planetary systems are likely to be composed of brown dwarfs and low mass stars of ～0.1–0.3 M_⊙, at orbital separations of ～ few × 100–10⁴ AU. In particular, scaling from solar-like stars suggests that hundreds of Mercury-like planets might orbit very massive stars at ～10³ AU where conditions might favor liquid water. Such fragmentation objects can be excellent targets for the James Webb Space Telescope and other large telescopes working in the IR bands. We predict that deep observations of very massive stars would reveal these fragmentation objects, orbiting in the same orbital plane in cases where there are more than one object.     

Characteristics of planetary-scale waves simulated by a new venusian mesosphere and thermosphere general circulation model

We have developed a new general circulation model (GCM) for the venusian mesosphere and thermosphere (80-about 180 km). Our GCM simulations show that winds in the subsolar-to-antisolar direction (SS–AS) are predominant above about 90 km. A weak return flow of the SS–AS is seen below about 90 km. We performed GCM simulations imposing the planetary-scale waves (thermal tides, Rossby wave, and Kelvin wave) at the lower boundary. Although the diurnal and semidiurnal tides are damped below 95 km, the Rossby wave propagates up to around 130 km. However, the amplitude of the Rossby wave is too small (<1 m/s) to affect the general circulation. On the other hand, the Kelvin wave propagates up to about 130 km with a maximum zonal wind fluctuation of approximately 5.9 m/s on average. The amplitude of the Kelvin wave sometimes exceeds 10 m/s around the terminator. The Kelvin wave causes a temporal variation in the wind velocity at the altitude of the O₂-1.27 μm nightglow emission (about 95 km). Using a newly developed 1-D nightglow model and the composition distribution calculated from our GCM, we investigated the impact of the Kelvin wave on the nightglow distribution. Our results suggest that the Kelvin wave would cause temporal variations in the nightglow emission in the 23:50–00:20 LT region with an intensity of 1.1–1.3 MR and a period of approximately 4 days.     

First ever in situ observations of Venus’ polar upper atmosphere density using the tracking data of the Venus Express Atmospheric Drag Experiment (VExADE)

On its highly elliptical 24 h orbit around Venus, the Venus Express (VEX) spacecraft briefly reaches a periapsis altitude of nominally 250 km. Recently, however, dedicated and intense radio tracking campaigns have taken place in August 2008, October 2009, February and April 2010, for which the periapsis altitude was lowered to the 186–176 km altitude range in order to be able to probe the upper atmosphere of Venus above the North Pole for the first time ever in situ. As the spacecraft experiences atmospheric drag, its trajectory is measurably perturbed during the periapsis pass, allowing us to infer total atmospheric mass density at the periapsis altitude. A Precise Orbit Determination (POD) of the VEX motion is performed through an iterative least-squares fitting process to the Doppler tracking data, acquired by the VEX radioscience experiment (VeRa). The drag acceleration is modelled using an initial atmospheric density model (VTS3 model, Hedin, A.E., Niemann, H.B., Kasprzak, W.T., Seiff, A. [1983]. J. Geophys. Res. 88, 73–83). A scale factor of the drag acceleration is estimated for each periapsis pass, which scales Hedin’s density model in order to best fit the radio tracking data. Reliable density scale factors have been obtained for 10 passes mainly from the second (October 2009) and third (April 2010) VExADE campaigns, which indicate a lower density by a factor of about 1.8 than Hedin’s model predicts. These first ever in situ polar density measurements at solar minimum have allowed us to construct a diffusive equilibrium density model for Venus’ thermosphere, constrained in the lower thermosphere primarily by SPICAV-SOIR measurements and above 175 km by the VExADE drag measurements (Müller-Wodarg et al., in preparation). The preliminary results of the VExADE campaigns show that it is possible to obtain with the POD technique reliable estimates of Venus’ upper atmosphere densities at an altitude of around 175 km. Future VExADE campaigns will benefit from the planned further lowering of VEX pericenter altitude to below 170 km.     

Latitudinal distribution of HDO abundance above Venus’ clouds by ground-based 2.3 μm spectroscopy

The abundance of HDO above the clouds in the dayside atmosphere of Venus was measured by ground-based 2.3 μm spectroscopy over 4 days. This is the first HDO observation above the clouds in this wavelength region corresponding to a new height region. The latitudinal distributions found show no clearly defined structure. The disk-averaged mixing ratio is 0.22 ± 0.03 ppm for a representative height region of 62–67 km. This is consistent with measurements found in previous studies. Based on previous H₂O measurements, the HDO/H₂O ratio is found to be 140 ± 20 times larger than the telluric ratio. This lies between the ratios of 120 ± 40 and 240 ± 25, respectively, reported for the 30–40 km region by ground-based nightside spectroscopy and for the 80–100 km region by solar occultation measurement on board the Venus Express.     

The planetary fourier spectrometer (PFS) onboard the European Venus Express mission

The planetary fourier spectrometer (PFS) for the Venus Express mission is an infrared spectrometer optimized for atmospheric studies. This instrument has a short wavelength (SW) channel that covers the spectral range from 1700 to 11400 cm⁻¹ (0.9–5.5 μm) and a long wavelength (LW) channel that covers 250–1700 cm⁻¹ (5.5–45 μm). Both channels have a uniform spectral resolution of 1.3 cm⁻¹. The instrument field of view FOV is about 1.6 ° (FWHM) for the short wavelength channel and 2.8 ° for the LW channel which corresponds to a spatial resolution of 7 and 12 km when Venus is observed from an altitude of 250 km. PFS can provide unique data necessary to improve our knowledge not only of the atmospheric properties but also surface properties (temperature) and the surface-atmosphere interaction (volcanic activity).PFS works primarily around the pericentre of the orbit, only occasionally observing Venus from larger distances. Each measurements takes 4.5 s, with a repetition time of 11.5 s. By working roughly 1.5 h around pericentre, a total of 460 measurements per orbit will be acquired plus 60 for calibrations. PFS is able to take measurements at all local times, enabling the retrieval of atmospheric vertical temperature profiles on both the day and the night side.The PFS measures a host of atmospheric and surface phenomena on Venus. These include the:(1) thermal surface flux at several wavelengths near 1 μm, with concurrent constraints on surface temperature and emissivity (indicative of composition); (2) the abundances of several highly-diagnostic trace molecular species; (3) atmospheric temperatures from 55 to 100 km altitude; (4) cloud opacities and cloud-tracked winds in the lower-level cloud layers near 50-km altitudes; (5) cloud top pressures of the uppermost haze/cloud region near 70–80 km altitude; and (6) oxygen airglow near the 100 km level. All of these will be observed repeatedly during the 500-day nominal mission of Venus Express to yield an increased understanding of meteorological, dynamical, photochemical, and thermo-chemical processes in the Venus atmosphere. Additionally, PFS will search for and characterize current volcanic activity through spatial and temporal anomalies in both the surface thermal flux and the abundances of volcanic trace species in the lower atmosphere.Measurement of the 15 μm CO₂ band is very important. Its profile gives, by means of a complex temperature profile retrieval technique, the vertical pressure-temperature relation, basis of the global atmospheric study.PFS is made of four modules called O, E, P and S being, respectively, the interferometer and proximity electronics, the digital control unit, the power supply and the pointing device.     

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.     

The mass,orbit, and tidal evolution of the Quaoar–Weywot system

Here we present new adaptive optics observations of the Quaoar–Weywot system. With these new observations we determine an improved system orbit. Due to a 0.39 day alias that exists in available observations, four possible orbital solutions are available with periods of ～11.6, ～12.0, ～12.4, and ～12.8 days. From the possible orbital solutions, system masses of 1.3–1.5 ± 0.1 × 10²¹ kg are found. These observations provide an updated density for Quaoar of 2.7–5.0 g cm^?3. In all cases, Weywot’s orbit is eccentric, with possible values ～0.13–0.16. We present a reanalysis of the tidal orbital evolution of the Quaoar–Weywot system. We have found that Weywot has probably evolved to a state of synchronous rotation, and has likely preserved its initial inclination over the age of the Solar System. We find that for plausible values of the effective tidal dissipation factor tides produce a very slow evolution of Weywot’s eccentricity and semi-major axis. Accordingly, it appears that Weywot’s eccentricity likely did not tidally evolve to its current value from an initially circular orbit. Rather, it seems that some other mechanism has raised its eccentricity post-formation, or Weywot formed with a non-negligible eccentricity.     

Comparisons of Cassini flybys of the Titan magnetospheric interaction with an MHD model: Evidence for organized behavior at high altitudes

Recent papers suggest the significant variability of conditions in Saturn’s magnetosphere at the orbit of Titan. Because of this variability, it was expected that models would generally have a difficult time regularly comparing to data from the Titan flybys. However, we find that in contrast to this expectation, it appears that there is underlying organization of the interaction features roughly above ～1800 km (1.7 Rt) altitude by the average external field due to Saturn’s dipole moment. In this study, we analyze Cassini’s plasma and magnetic field data collected at 9 Titan encounters during which the external field is close to the ideal southward direction and compare these observations to the results from a 2-fluid (1 ion, 1 electron) 7-species MHD model simulations obtained under noon SLT conditions. Our comparative analysis shows that under noon SLT conditions the Titan plasma interaction can be viewed in two layers: an outer layer between 6400 and 1800 km where interaction features observed in the magnetic field are in basic agreement with a purely southward external field interaction and an inner layer below 1800 km where the magnetic field measurements show strong variations and deviate from the model predictions. Thus the basic features inferred from the Voyager 1 flyby seem to be generally present above ～1800 km in spite of the ongoing external variations from SLT excursions, time variability and magnetospheric current systems as long as a significant southward external field component is present. At around ～1800 km kinetic effects (such as mass loading and heavy ion pickup) and below 1800 km ionospheric effects (such as drag of ionospheric plasma due to coupling with neutral winds and/or magnetic memory of Titan’s ionosphere) complicate what is observed.     

An updated period analysis for ER Orionis: A definitive solution

An updated period analysis for the overcontact eclipsing binary ER Orionis is presented. Featured is an improved derivation of parameters for the light time effect (LTE) due to the third star (in actuality, a pair of stars) utilising the latest set of eclipse timings. The very good fit between the eclipse timing differences (ETD) plot (otherwise known as an O–C diagram) and the theoretical ETD curve makes possible an improved determination of the rate of mass interchange between the binary pair, dm₁/dt = +1.83(6) × 10⁻⁷ M_ʘ/year. In addition, the mass of the companion system (in actuality, m₃ sin i) and the elements of its orbit were computed. A suggestion is made for a method of future determination of the inclination of the orbit of the companion system.     

Assessing the long-term variability of Venus winds at cloud level from VIRTIS–Venus Express

The Venus Express (VEX) mission has been in orbit to Venus for more than 4 years now. The Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) instrument onboard VEX observes Venus in two channels (visible and infrared) obtaining spectra and multi-wavelength images of the planet that can be used to sample the atmosphere at different altitudes. Day-side images in the ultraviolet range (380 nm) are used to study the dynamics of the upper cloud at 66–72 km while night-side images in the near infrared (1.74 μm) map the opacity of the lower cloud deck at 44–48 km. Here we present a long-term analysis of the global atmospheric dynamics at these levels using a large selection of orbits from the VIRTIS-M dataset covering 860 Earth days that extends our previous work (Sánchez-Lavega, A. et al. [2008]. Geophys. Res. Lett. 35, L13204) and allows studying the variability of the global circulation at the two altitude levels. The atmospheric superrotation is evident with equatorial to mid-latitudes westward velocities of 100 and 60 m s^?1 in the upper and lower cloud layers. These zonal velocities are almost constant in latitude from the equator to 50°S. From 50°S to 90°S the zonal winds at both cloud layers decrease steadily to zero at the pole. Individual cloud tracked winds have errors of 3–10 m s^?1 with a mean of 5 m s^?1 and the standard deviations for a given latitude of our zonal and meridional winds are 9 m s^?1. The zonal winds in the upper cloud change with the local time in a way that can be interpreted in terms of a solar tide. The zonal winds in the lower cloud are stable at mid-latitudes to the tropics and present variability at subpolar latitudes apparently linked to the activity of the South polar vortex. While the upper cloud presents a net meridional motion consistent with the upper branch of a Hadley cell with peak velocity v = 10 m s^?1 at 50°S, the lower cloud meridional motions are less organized with some cloud features moving with intense northwards and southwards motions up to v = ±15 m s^?1 but, on average, with almost null global meridional motions at all latitudes. We also examine the long-term behavior of the winds at these two vertical layers by comparing our extended wind tracked data with results from previous missions.     

Vertical structure of the Venus cloud top from the VeRa and VIRTIS observations onboard Venus Express

We investigate the Venus cloud top structure by joint analysis of the data from Visual and Thermal Infrared Imaging Spectrometer (VIRTIS) and the atmospheric temperature sounding by the Radio Science experiment (VeRa) onboard Venus Express. The cloud top altitude and aerosol scale height are derived by fitting VIRTIS spectra at 4–5 μm with temperature profiles taken from the VeRa radio occultation. Our study shows gradual descent of the cloud top from 67.2 ± 1.9 km in low latitudes to 62.8 ± 4.1 km at the pole and decrease of the aerosol scale height from 3.8 ± 1.6 km to 1.7 ± 2.4 km. These changes correlate with the mesospheric temperature field. In the cold collar and high latitudes the cloud top position remarkably coincides with the sharp minima in temperature inversions suggesting importance of radiative cooling in their maintenance. This behaviour is consistent with the earlier observations. Spectral trend of the cloud top altitude derived from a comparison with the earlier observations in 1.6–27 μm wavelength range is qualitatively consistent with sulphuric acid composition of the upper cloud and suggests that particle size increases from equator to pole.     

Vertical profiling of SO2 and SO above Venus’ clouds by SPICAV/SOIR solar occultations

New measurements of sulfur dioxide (SO₂) and monoxide (SO) in the atmosphere of Venus by SPICAV/SOIR instrument onboard Venus Express orbiter provide ample statistics to study the behavior of these gases above Venus’ clouds. The instrument (a set of three spectrometers) is capable to sound atmospheric structure above the clouds in several observation modes (nadir, solar and stellar occultations) either in the UV or in the near IR spectral ranges. We present the results from solar occultations in the absorption ranges of SO₂ (190–230 nm, and at 4 μm) and SO (190–230 nm). The dioxide was detected by the SOIR spectrometer at the altitudes of 65–80 km in the IR and by the SPICAV spectrometer at 85–105 km in the UV. The monoxide’s absorption was measured only by SPICAV at 85–105 km. We analyzed 39 sessions of solar occultation, where boresights of both spectrometers are oriented identically, to provide complete vertical profiling of SO₂ of the Venus’ mesosphere (65–105 km). Here we report the first firm detection and measurements of two SO₂ layers. In the lower layer SO₂ mixing ratio is within 0.02–0.5 ppmv. The upper layer, also conceivable from microwave measurements by Sandor et al. (Sandor, B.J., Todd Clancy, R., Moriarty-Schieven, G., Mills, F.P. [2010]. Icarus 208, 49–60) is characterized by SO₂ increasing with the altitude from 0.05 to 2 ppmv, and the [SO₂]/[SO] ratio varying from 1 to 5. The presence of the high-altitude SO_x species could be explained by H₂SO₄ photodissociation under somewhat warmer temperature conditions in Venus mesosphere. At 90–100 km the content of the sulfur dioxide correlates with temperature increasing from 0.1 ppmv at 165–170 K to 0.5–1 ppmv at 190–192 K. It supports the hypothesis of SO₂ production by the evaporation of H₂SO₄ from droplets and its subsequent photolysis at around 100 km.     

Upper limits for H2SO4 in the mesosphere of Venus

Rapid temporal variability of SO₂ and SO in the Venus 85–100 km mesosphere (Sandor, B.J., Clancy, R.T., Moriarty-Schieven G.H. [2007]. Bull. Am. Astron. Soc. 39, 503; Sandor, B.J., Clancy, R.T., Moriarty-Schieven, G.H., Mills, F.P. [2010]. Icarus 208, 49–60) requires in situ sources and sinks for these molecules. While many loss mechanisms are recognized, no process for in situ production is known. Observational investigations to find, or constrain other potential sulfur reservoirs offer one method toward understanding the applicable photochemistry. Here, we report upper limits for gas-phase H₂SO₄ (sulfuric acid) abundances in Venus’ 85–100 km upper mesosphere, derived from 16 ground-based sub-mm spectroscopic observations in the period 2004–2008. Unlike the ubiquitous sulfuric acid solid/liquid aerosol, the gas phase would be photochemically active, potentially both source and sink for SO and SO₂. H₂SO₄ is retrieved from sub-mm lines located in the same bandpass as the SO₂ and SO lines described by Sandor et al. (Sandor, B.J., Clancy, R.T., Moriarty-Schieven, G.H., Mills, F.P. [2010]. Icarus 208, 49–60). H₂SO₄ upper limits reported here are thus simultaneous and spatially coincident with measurements of SO₂ and SO, providing for analysis of the three sulfur species collectively. The average H₂SO₄ abundance over 16 observations is 1 ± 2 ppb (i.e. <3 ppb). Upper limits for individual observations range from 3 to 44 ppb, where quality of the observing weather is the dominant constraint on measurement precision. The sum of H₂SO₄, SO₂ and SO varies widely. In one comparison, the sum [H₂SO₄ + SO₂ + SO] measured on one date differs by 10-σ from the sum measured 2 months later. We conclude that upper mesospheric sulfur atoms are not conserved among the three molecules, that H₂SO₄ is not a significant sulfur reservoir for balancing the observed variations of [SO₂ + SO], and is not relevant to the (still unknown) photochemistry responsible for observed behavior of SO₂ and SO. Having ruled out H₂SO₄, we infer that elemental sulfur is the most probable candidate for the needed third reservoir.     

Small-scale variability in Saturn’s lower ionosphere

We perform and present a wavelet analysis on all 31 Cassini electron density profiles published to date (Nagy, A.F. et al. [2006]. J. Geophys. Res. 111 (A6), CiteID A06310; Kliore, A.J. et al. [2009]. J. Geophys. Res. 114 (A4), CiteID A04315). We detect several discrete scales of variability present in the observations. Small-scale variability (S < 700 km) is observed in almost all data sets at different latitudes, both at dawn and dusk conditions. The most typical scale of variability is 300 km with scales between 200 km and 450 km being commonly present in the vast majority of the profiles. A low latitude dawn/dusk asymmetry is noted in the prevalent scales with the spectrum peaking sharply at the 300 km scale at dusk conditions and being broader at dawn conditions. Compared to dawn conditions the dusk ionosphere also shows more significant variability at the 100 km scale. The 300 km vertical scale is also present in the few available profiles from the northern hemisphere. Early observations from 2005 show a dominant scale at 350 km whereas later in 2007–2008 the spectrum shifts to the shorter scales with the most prominent scale being 300 km. The performed wavelet analysis and the obtained results are independent of assumptions about the nature of the layers and do not require a definition for a “background” electron density profile.In the second part of the paper we present a gravity wave propagation/dissipation model for Saturn’s upper atmosphere and compare the wave properties to the characteristics of the observed electron density variability at different scales. The general features observed in the data are consistent with gravity waves being present in the lower ionosphere and causing layering of the ions and the electrons. The wave-driving mechanism provides a simultaneous explanation for several of the properties of the observed variability: (i) lack of variability in the electron density above the predicted region of wave dissipation; (ii) in most cases the peak amplitude of variability occurs within the altitude range for dissipation of gravity waves or below; (iii) shorter scales have smaller amplitudes than the longer scales; (iv) shorter scales are present at lower altitudes whereas longer scales persist to higher altitudes; and (v) several layers often form a system of equally spaced maxima and minima that can be traced over a large altitude range.