北斗系统测站钟差短期预报模型比较及其在单星定轨中的应用

北斗卫星导航系统(BDS)地面跟踪站都配置有高精度的氢原子钟,并基于精密定轨数据处理与主站的时间基准进行同步.在卫星轨道机动以及机动恢复期间,通常采用几何法定轨以及单星定轨确定卫星的轨道.而在这两种定轨模式中,需要提供精确的测站钟差作为输入.为提高定轨的实时性,需要对测站钟差进行预报处理.分析了2次多项式模型、附加周期项模型、灰色模型3种模型对北斗地面跟踪站钟差短期拟合和预报的性能,并将钟差预报结果应用于单星定轨,同时还分析了不同预报钟差用于定轨的精度.试验发现,以上3种模型对6个测站钟差的平均拟合精度分别为0.14 ns、0.05 ns、0.27 ns,预报1 h的平均精度分别为1.17 ns、0.88 ns、1.28 ns,预报2 h的平均精度分别为2.72 ns、2.09 ns、2.53 ns.采用3种模型对测站钟差进行预报并用于单星定轨,采用附加周期项的钟差预报模型轨道3维误差最小,不同模型轨道径向精度差异在3 cm以内.以上结果表明,附加周期项的站钟拟合及预报模型在北斗系统机动期间的轨道恢复数据处理具有最好的效果.     

An Accuracy Analysis of the SGP4/SDP4 Model

Based on the latest release of the SGP4/SDP4 (Simplified General Perturbation Version 4/ Simplified Deep-space Perturbation Version 4) model, in this paper we have designed an orbit determination program. Through calculations for the 1120 objects with various types and orbital elements selected from the space objects database, we have obtained the accuracies of the orbit determination prediction dealt with various types of space objects by the SGP4/SDP4 model. The results show that the accuracies of the near-earth objects are in the order of magnitude of 100 meters; the averages of the orbit determination accuracies of the semi-synchronous and geosynchronous orbits are, respectively, 0.7 and 1.9 km. The orbit determination accuracies of the elliptical orbit objects are related to their eccentricities. Except for few elliptical orbit objects with e > 0.8, the orbit determination errors of the vast majority of the elliptical orbit objects are all less than 10 km. By using the SGP4/SDP4 model to make 3 days predictions for near-earth objects, 30 days for semi-synchronous orbit objects, 15 days for geosynchronous orbit objects and 1 day for elliptical orbit objects, the errors of prediction generally don’t exceed 40 km.     

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.     

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.     

Solar modulation of galactic cosmic rays during 19–23 solar cycles

We make a detailed analysis of cross-correlation and time-lag between monthly data of galactic cosmic rays (GCRs) intensity and different solar activity indices (e.g., sunspot number, sunspot area, green coronal Fe line and 10.7 cm solar radio flux) during 19–23 solar cycles. GCRs time-series data from Kiel neutron monitor station and solar data from the last 50 years period, covering five solar cycles (19–23), and alternating solar polarity states (i.e., five A < 0 and four A > 0) have been investigated. We find a clear asymmetry in the cross-correlation between GCRs and solar activity indicators for both odd and even-numbered solar cycles. The time-lags between GCRs and solar parameters are found different in different solar cycles as well as in the opposite polarity states (A < 0 and A > 0) within the same solar cycle. Possible explanations of the observed results are discussed in light of modulation models, including drift effects.     

The first extrasolar planet detected via gravitational microlensing

In gravitational microlensing, distant planetary systems may be discovered by utilizing them as naturally occuring lenses. Efforts to find planets by this technique began in the 1990s. The first definitive detection of an extrasolar planet by microlensing was made in 2003 in the event OGLE 2003-BLG-235/MOA 2003-BLG-53, where the observed light curve was best reproduced using a binary microlensing model with a mass ratio of 0.004. Further observations with the HST revealed that the lens system comprises a 2.6 Jupiter mass planet in a 4.3 A.U. wide orbit around a 0.6 Solar mass K dwarf at a distance of 5.8 Kpc. Subsequently, the number of planets detected by microlensing is increasing.     

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.     

Prediction of Clock Errors of Atomic Clocks Based on Modified Linear Combination Model†

The prediction of the clock errors of atomic clocks plays an important role in the work on time and frequency. Each of the prediction models often used at present has its own merits and shortages. A combination of the predicted results obtained by means of these models can be used to synthesize the characteristics of various kinds of prediction models. In the light of the problem which occurs when the linear combination model is used to make the prediction of clock errors, the concept of learning weight is proposed and the modified combination prediction model is made by taking advantage of various kinds of pieces of accuracy information. For verifying the efficiency of this method the clock error sequences of the IGS (International GNSS Service) of 4 GPS satellites are selected and the predicted results of the quadratic polynomial and grey model are combined. The result shows that the modified model can further improve the stability and accuracy based on the guarantee of the reliability.     

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.     

The abundance and vertical distribution of the unknown ultraviolet absorber in the venusian atmosphere from analysis of Venus Monitoring Camera images

Observations of Venus using the ultraviolet filter of the Venus Monitoring Camera (VMC) on ESA’s Venus Express Spacecraft (VEX) provide the best opportunity for study of the spatial and temporal distribution of the venusian unknown ultraviolet absorber since the Pioneer Venus (PV) mission. We compare the results of two sets of 125 radiative transfer models of the upper atmosphere of Venus to each pixel in a subset of VMC UV channel images. We use a quantitative best fit criterion based upon the notion that the distribution of the unknown absorber should be independent of the illumination and observing geometry. We use the product of the cosines of the incidence and emission angles and search for absorber distributions that are uncorrelated with this geometric parameter, finding that two models can describe the vertical distribution of the unknown absorber. One model is a well-mixed vertical profile above a pressure level of roughly 120 mb (～63 km). This is consistent with the altitude of photochemical formation of sulfuric acid. The second model describes it as a thin layer of pure UV absorber at a pressure level roughly around 24 mb (～71 km) and this altitude is consistent with the top of upper cloud deck. We find that the average abundance of unknown absorber in the equatorial region is 0.21 ± 0.04 optical depth and it decreases in the polar region to 0.08 ± 0.05 optical depth at 365 nm.     

Thermospheric/mesospheric temperatures on Venus: Results from ground-based high-resolution spectroscopy of CO2 in 1990/1991 and comparison to results from 2009 and between other techniques

We report temperatures in Venus’ upper mesosphere/lower thermosphere, deduced from reanalyzing very high resolution infrared spectroscopy of CO₂ emission lines acquired in 1990 and 1991. Kinetic temperatures at ～110 km altitude (0.15 Pa) are derived from the Doppler width of fully-resolved single line profiles measured near 10.4 μm wavelength using the NASA GSFC Infrared Heterodyne Spectrometer (IRHS) at the NASA IRTF on Mauna Kea, HI, close to Venus inferior conjunction and two Venus solstices. Measured temperatures range from ～200 to 240 K with uncertainty typically less than 10 K. Temperatures retrieved from similar measurement in 2009 using the Cologne Tuneable Heterodyne Infrared Spectrometer (THIS) at the NOAO McMath Telescope at Kitt Peak, AZ are 10–20 K lower. Temperatures retrieved more recently from the SOIR instrument on Venus EXpress are consistent with these results when the geometry of observation is accounted for. It is difficult to compare ground-based sub-mm retrievals extrapolated to 110 km due to their much larger field of view, which includes the night side regions not accessible to infrared heterodyne observations. Temperature variability appears to be high on day-to-day as well as longer timescales. Observed short term and long term variability may be attributed to atmospheric dynamics, diurnal variability and changes over solar activity and seasons. The Venus International Reference Atmosphere (VIRA) model predicts cooler temperatures at the sampled altitudes in the lower thermosphere/upper mesosphere and is not consistent with these measurements.     

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.     

Interpretation of combined infrared,submillimeter, and millimeter thermal flux data obtained during the Rosetta fly-by of Asteroid (21) Lutetia

The European Space Agency’s Rosetta spacecraft is the first Solar System mission to include instrumentation capable of measuring planetary thermal fluxes at both near-IR (VIRTIS) and submillimeter–millimeter (smm–mm, MIRO) wavelengths. Its primary mission is a 1 year reconnaissance of Comet 67P/Churyumov–Gerasimenko beginning in 2014. During a 2010 close fly-by of Asteroid 21 Lutetia, the VIRTIS and MIRO instruments provided complementary data that have been analyzed to produce a consistent model of Lutetia’s surface layer thermal and electrical properties, including a physical model of self-heating. VIRTIS dayside measurements provided highly resolved 1 K accuracy surface temperatures that required a low thermal inertia, I < 30 J/(K m² s^0.5). MIRO smm and mm measurements of polar night thermal fluxes produced constraints on Lutetia’s subsurface thermal properties to depths comparable to the seasonal thermal wave, yielding a model of I < 20 J/(K m² s^0.5) in the upper few centimeters, increasing with depth in a manner very similar to that of Earth’s Moon. Subsequent MIRO-based model predictions of the dayside surface temperatures reveal negative offsets of ～5–30 K from the higher VIRTIS-measurements. By adding surface roughness in the form of 50% fractional coverage of hemispherical mini-craters to the MIRO-based thermal model, sufficient self-heating is produced to largely remove the offsets relative to the VIRTIS measurements and also reproduce the thermal limb brightening features (relative to a smooth surface model) seen by VIRTIS. The Lutetia physical property constraints provided by the VIRTIS and MIRO data sets demonstrate the unique diagnostic capabilities of combined infrared and submillimeter/millimeter thermal flux measurements.     

Cyclic period changes and the light-time effect in eclipsing binaries: A low-mass companion around the system VV Ursae Majoris

In this article, a period analysis of the late-type eclipsing binary VV UMa is presented. This work is based on the periodic variation of eclipse timings of the VV UMa binary. We determined the orbital properties and mass of a third orbiting body in the system by analyzing the light-travel time effect. The O−C diagram constructed for all available minima times of VV UMa exhibits a cyclic character superimposed on a linear variation. This variation includes three maxima and two minima within approximately 28,240 orbital periods of the system, which can be explained as the light-travel time effect (LITE) because of an unseen third body in a triple system that causes variations of the eclipse arrival times. New parameter values of the light-time travel effect because of the third body were computed with a period of 23.22 ± 0.17 years in the system. The cyclic-variation analysis produces a value of 0.0139 day as the semi-amplitude of the light-travel time effect and 0.35 as the orbital eccentricity of the third body. The mass of the third body that orbits the eclipsing binary stars is 0.787 ± 0.02 M_⊙, and the semi-major axis of its orbit is 10.75 AU.     

The effect of orbital evolution on the Haumea (2003 EL61) collisional family

The Haumea family is currently the only identified collisional family in the Kuiper belt. We numerically simulate the long-term dynamical evolution of the family to estimate a lower limit of the family’s age and to assess how the population of the family and its dynamical clustering are preserved over Gyr timescales. We find that the family is not younger than 100 Myr, and its age is at least 1 Gyr with 95% confidence. We find that for initial velocity dispersions of 50–400 m s^?1, approximately 20–45% of the family members are lost to close encounters with Neptune after 3.5 Gyr of orbital evolution. We apply these loss rates to two proposed models for the formation of the Haumea family, a graze-and-merge type collision between two similarly sized, differentiated KBOs or the collisional disruption of a satellite orbiting Haumea. For the graze-and-merge collision model, we calculate that >85% of the expected mass in surviving family members within 150 m s^?1 of the collision has been identified, but that one to two times the mass of the known family members remains to be identified at larger velocities. For the satellite-break-up model, we estimate that the currently identified family members account for ～50% of the expected mass of the family. Taking observational incompleteness into account, the observed number of Haumea family members is consistent with either formation scenario at the 1σ level, however both models predict more objects at larger relative velocities (>150 m s^?1) than have been identified.     

Ultraviolet spectral variations of FK Comae and V 1794 Cygni

Ultraviolet spectra of FK Comae and V1794 Cygni observed with the Hubble Space Telescope Cosmic Origins Spectrograph (HST COS) and the International Ultraviolet Explorer (IUE) satellites were analyzed for the period 1981–2011. Temporal variations of line fluxes of the O I 1306 Å, C II 1336 Å, C IV 1550 Å, He II 1640 Å and Mg II k & h 2800 Å, produced in the transition regions and chromospheres of these stars, imply variations in density and temperature changes in the line emitting regions as a result of the rapid rotation and magnetic fields responsible for stellar activity. Results are consistent with the models of Ramsey et al. (1981), Oliveira and Foing (1999), and Korhonen et al. (2000).     

A photochemical model for the Venus atmosphere at 47–112 km

The model is intended to respond to the recent findings in the Venus atmosphere from the Venus Express and ground-based submillimeter and infrared observations. It extends down to 47 km for comparison with the kinetic model for the lower atmosphere (Krasnopolsky, V.A. [2007]. Icarus 191, 25–37) and to use its results as the boundary conditions. The model numerical accuracy is significantly improved by reduction of the altitude step from 2 km in the previous models to 0.5 km. Effects of the NUV absorber are approximated using the detailed photometric observations at 365 nm from Venera 14. The H₂O profile is not fixed but calculated in the model. The model involves odd nitrogen and OCS chemistries based on the detected NO and OCS abundances. The number of the reactions is significantly reduced by removing of unimportant processes. Column rates for all reactions are given, and balances of production and loss may be analyzed in detail for each species.The calculated vertical profiles of CO, H₂O, HCl, SO₂, SO, OCS and of the O₂ dayglow at 1.27 μm generally agree with the existing observational data; some differences are briefly discussed. The OH dayglow is ～30 kR, brighter than the OH nightglow by a factor of 4. The H + O₃ process dominates in the nightglow excitation and O + HO₂ in the dayglow, because of the reduction of ozone by photolysis. A key feature of Venus’ photochemistry is the formation of sulfuric acid in a narrow layer near the cloud tops that greatly reduces abundances of SO₂ and H₂O above the clouds. Delivery of SO₂ and H₂O through this bottleneck determines the chemistry and its variations above the clouds. Small variations of eddy diffusion near 60 km result in variations of SO₂, SO, and OCS at and above 70 km within a factor of ～30. Variations of the SO₂/H₂O ratio at the lower boundary have similar but weaker effect: the variations within a factor of ～4 are induced by changes of SO₂/H₂O by ±5%. Therefore the observed variations of the mesospheric composition originate from minor variations of the atmospheric dynamics near the cloud layer and do not require volcanism. NO cycles are responsible for production of a quarter of O₂, SO₂, and Cl₂ in the atmosphere. A net effect of photochemistry in the middle atmosphere is the consumption of CO₂, SO₂, and HCl from and return of CO, H₂SO₄, and SO₂Cl₂ to the lower atmosphere. These processes may be balanced by thermochemistry in the lower atmosphere even without outgassing from the interior, though the latter is not ruled out by our models. Some differences between the model and observations and the previous models are briefly discussed.     

Sulfur chemistry in the middle atmosphere of Venus

Venus Express measurements of the vertical profiles of SO and SO₂ in the middle atmosphere of Venus provide an opportunity to revisit the sulfur chemistry above the middle cloud tops (～58 km). A one dimensional photochemistry-diffusion model is used to simulate the behavior of the whole chemical system including oxygen-, hydrogen-, chlorine-, sulfur-, and nitrogen-bearing species. A sulfur source is required to explain the SO₂ inversion layer above 80 km. The evaporation of the aerosols composed of sulfuric acid (model A) or polysulfur (model B) above 90 km could provide the sulfur source. Measurements of SO₃ and SO (a¹Δ → X³Σ^-) emission at 1.7 μm may be the key to distinguish between the two models.     

The orbital distribution of Near-Earth Objects inside Earth’s orbit

Canada’s Near-Earth Object Surveillance Satellite (NEOSSat), set to launch in early 2012, will search for and track Near-Earth Objects (NEOs), tuning its search to best detect objects with a < 1.0 AU. In order to construct an optimal pointing strategy for NEOSSat, we needed more detailed information in the a < 1.0 AU region than the best current model (Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.M., Levison, H.F., Michel, P., Metcalfe, T.S. [2002]. Icarus 156, 399–433) provides. We present here the NEOSSat-1.0 NEO orbital distribution model with larger statistics that permit finer resolution and less uncertainty, especially in the a < 1.0 AU region. We find that Amors = 30.1 ± 0.8%, Apollos = 63.3 ± 0.4%, Atens = 5.0 ± 0.3%, Atiras (0.718 < Q < 0.983 AU) = 1.38 ± 0.04%, and Vatiras (0.307 < Q < 0.718 AU) = 0.22 ± 0.03% of the steady-state NEO population. Vatiras are a previously undiscussed NEO population clearly defined in our integrations, whose orbits lie completely interior to that of Venus. Our integrations also uncovered the unexpected production of retrograde orbits from main-belt asteroid sources; this retrograde NEA population makes up ?0.1% of the steady-state NEO population. The relative NEO impact rate onto Mercury, Venus, and Earth, as well as the normalized distribution of impact speeds, was calculated from the NEOSSat-1.0 orbital model under the assumption of a steady-state. The new model predicts a slightly higher Mercury impact flux.