Impact of atmospheric gravity waves on the jovian ionosphere

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.     

On the population, physical decay and orbital distribution of Jupiter family comets: Numerical simulations

We study the Jupiter family comet (JFC) population assumed to come from the Scattered Disk and transferred to the Jupiter’s zone through gravitational interactions with the Jovian planets. We shall define as JFCs those with orbital periods and Tisserand parameters in the range 2<T?3.1, while those comets coming from the same source, but that do not fulfill the previous criteria (mainly because they have periods ) will be called ‘non-JFCs’. We performed a series of numerical simulations of fictitious comets with a purely dynamical model and also with a more complete dynamical-physical model that includes besides nongravitational forces, sublimation and splitting mechanisms. With the dynamical model, we obtain a poor match between the computed distributions of orbital elements and the observed ones. However with the inclusion of physical effects in the complete model we are able to obtain good fits to observations. The best fits are attained with four splitting models with a relative weak dependence on q, and a mass loss in every splitting event that is less when the frequency is high and vice versa. The mean lifetime of JFCs with radii and is found to be of about 150-200 revolutions (∼. The total population of JFCs with radii within Jupiter’s zone is found to be of 450±50. Yet, the population of non-JFCs with radii in Jupiter-crossing orbits may be ∼4 times greater, thus leading to a whole population of JFCs + non-JFCs of ∼2250±250. Most of these comets have perihelia close to Jupiter’s orbit. On the other hand, very few non-JFCs reach the Earth’s vicinity (perihelion distances ) which gives additional support to the idea that JFCs and Halley-type comets have different dynamical origins. Our model allows us to define the zones of the orbital element space in which we would expect to find a large number of JFCs. This is the first time, to our knowledge, that a physico-dynamical model is presented that includes sublimation and different splitting laws. Our work helps to understand the role played by these erosion effects in the distribution of the orbital elements and lifetimes of JFCs.     

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.     

Calculation of the enrichment of the giant planet envelopes during the “late heavy bombardment”

The giant planets of our solar system possess envelopes consisting mainly of hydrogen and helium but are also significantly enriched in heavier elements relatively to our Sun. In order to better constrain how these heavy elements have been delivered, we quantify the amount accreted during the so-called “late heavy bombardment”, at a time when planets were fully formed and planetesimals could not sink deep into the planets. On the basis of the “Nice model”, we obtain accreted masses (in terrestrial units) equal to for Jupiter, and for Saturn. For the two other giant planets, the results are found to depend mostly on whether they switched position during the instability phase. For Uranus, the accreted mass is with an inversion and without an inversion. Neptune accretes in models in which it is initially closer to the Sun than Uranus, and otherwise. With well-mixed envelopes, this corresponds to an increase in the enrichment over the solar value of 0.033±0.001 and 0.074±0.007 for Jupiter and Saturn, respectively. For the two other planets, we find the enrichments to be 2.1±1.4 (w/ inversion) or 1.2±0.7 (w/o inversion) for Uranus, and 2.0±1.2 (w/ inversion) or 2.7±1.6 (w/o inversion) for Neptune. This is clearly insufficient to explain the inferred enrichments of ∼4 for Jupiter, ∼7 for Saturn and ∼45 for Uranus and Neptune.     

Saturn eddy momentum fluxes and convection: First estimates from Cassini images

We apply an automated cloud feature tracking algorithm to estimate eddy momentum fluxes in Saturn's southern hemisphere from Cassini Imaging Science Subsystem near-infrared continuum image sequences. Voyager Saturn manually tracked images had suggested no conversion of eddy to mean flow kinetic energy, but this was based on a small sample of <1000 wind vectors. The automated procedure we use for the Cassini data produces an order of magnitude more usable wind vectors with relatively unbiased sampling. Automated tracking is successful in and around the westward jet latitudes on Saturn but not in the vicinity of most eastward jets, where the linearity and non-discrete nature of cloud features produces ambiguous results. For the regions we are able to track, we find peak eddy fluxes and a clear positive correlation between eddy momentum fluxes and meridional shear of the mean zonal wind, implying that eddies supply momentum to eastward jets and remove momentum from westward jets at a rate . The behavior we observe is similar to that seen on Jupiter, though with smaller eddy-mean kinetic energy conversion rates per unit mass of atmosphere (). We also use the appearance and rapid evolution of small bright features at continuum wavelengths, in combination with evidence from weak methane band images where possible, to diagnose the occurrence of moist convective storms on Saturn. Areal expansion rates imply updraft speeds of over the convective anvil cloud area. As on Jupiter, convection preferentially occurs in cyclonic shear regions on Saturn, but unlike Jupiter, convection is also observed in eastward jet regions. With one possible exception, the large eddy fluxes seen in the cyclonic shear latitudes do not seem to be associated with convective events.     

The luminosity function of the hot and cold Kuiper belt populations

We have performed an ecliptic survey of the Kuiper belt, with an areal coverage of 8.9 square degrees to a 50% limiting magnitude of , and have detected 88 Kuiper belt objects, roughly half of which received follow-up 1–2 months after detection. Using this survey data alone, we have measured the luminosity function of the Kuiper belt, thus avoiding any biases that might come from the inclusion of other observations. We have found that the Cold population defined as having inclinations less than 5° has a luminosity function slope α_Cold = 0.82 ± 0.23, and is different from the Hot population, which has inclinations greater than 5° and a luminosity function slope α_Hot = 0.35 ± 0.21. As well, we have found that those objects closer than 38 AU have virtually the same luminosity function slope as the Hot population. This result, along with similar findings of past surveys demonstrates that the dynamically Cold Kuiper belt objects likely have a steep size distribution, and are unique from all of the excited populations which have much shallower distributions. This suggests that the dynamically excited population underwent a different accretion history and achieved a more evolved state of accretion than the Cold population. As well, we discuss the similarities of the Cold and Hot populations with the size distributions of other planetesimal populations. We find that while the Jupiter family comets and the scattered disk exhibit similar size distributions, a power-law extrapolation to small sizes for the scattered disk cannot account for the observed influx of comets. As well, we have found that the Jupiter Trojan and Hot populations cannot have originated from the same parent population, a result that is difficult to reconcile with scattering models similar to the NICE model. We conclude that the similarity between the size distributions of the Cold population and the Jupiter Trojan population is a striking coincidence.     

Persistent evidence of a jovian mass solar companion in the Oort cloud

We present updated dynamical and statistical analyses of outer Oort cloud cometary evidence suggesting that the Sun has a wide-binary jovian mass companion. The results support a conjecture that there exists a companion of mass ≈ orbiting in the innermost region of the outer Oort cloud. Our most restrictive prediction is that the orientation angles of the orbit plane in galactic coordinates are centered on Ω, the galactic longitude of the ascending node = 319° and i, the galactic inclination = 103° (or the opposite direction) with an uncertainty in the orbit normal direction subtending <2% of the sky. Such a companion could also have produced the detached Kuiper Belt object Sedna. If the object exists, the absence of similar evidence in the inner Oort cloud implies that common beliefs about the origin of observed inner Oort cloud comets must be reconsidered. Evidence of the putative companion would have been recorded by the Wide-field Infrared Survey Explorer (WISE) which has completed its primary mission and is continuing on secondary objectives.     

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

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.     

Fragment properties at the catastrophic disruption threshold: The effect of the parent body’s internal structure

Numerical simulations of asteroid breakups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by , which results in the escape of half of the target’s mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of , and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of on the impact speed of the projectile.In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. Indeed, decreasing the strength of non-porous targets make them easier to disrupt in this regime, while increasing the strength of porous targets has much less influence on the value of . Therefore, one cannot say that non-porous targets are systematically easier or more difficult to disrupt than porous ones, as the outcome highly depends on the assumed strength values. In the gravity regime, we also confirm that the process of gravitational reaccumulation is at the origin of the largest remnant’s mass in both cases. We then propose some power-law relationships between and both target’s size and impact speed that can be used in collisional evolution models. The resulting fragment size distributions can also be reasonably fitted by a power-law whose exponent ranges between −2.2 and −2.7 for all target diameters in both cases and independently on the impact velocity (at least in the small range investigated between 3 and 5 km/s). Then, although ejection velocities in the gravity regime tend to be higher from porous targets, they remain on the same order as the ones from non-porous targets.     

Eclipsing binary Trojan asteroid Patroclus: Thermal inertia from Spitzer observations

We present mid-infrared observations of the binary L5-Trojan system (617) Patroclus-Menoetius before, during, and after two shadowing events, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. For the first time, we effectively observe changes in asteroid surface temperature in real time, allowing the thermal inertia to be determined very directly. A new detailed binary thermophysical model is presented which accounts for the system’s known mutual orbit, arbitrary component shapes, and thermal conduction in the presence of eclipses.We obtain two local thermal-inertia values, representative of the respective shadowed areas: and . The average thermal inertia is estimated to be , potentially with significant surface heterogeneity. This first thermal-inertia measurement for a Trojan asteroid indicates a surface covered in fine regolith. Independently, we establish the presence of fine-grained (<a few μm) silicates on the surface, based on emissivity features near 10 and similar to those previously found on other Trojans.We also report V-band observations and report a lightcurve with complete rotational coverage. The lightcurve has a low amplitude of peak-to-peak, implying a roughly spherical shape for both components, and is single-periodic with a period equal to the period of the mutual orbit, indicating that the system is fully synchronized.The diameters of Patroclus and Menoetius are 106±11 and , respectively, in agreement with previous findings. Taken together with the system’s known total mass, this implies a bulk mass density of , significantly below the mass density of L4-Trojan asteroid (624) Hektor and suggesting a bulk composition dominated by water ice.All known physical properties of Patroclus, arguably the best studied Trojan asteroid, are consistent with those expected in icy objects with devolatilized surface (extinct comets), consistent with what might be implied by recent dynamical modeling in the framework of the Nice Model.     

Interferometric millimeter observations of water vapor on Mars and comparison with Mars Express measurements

We present interferometric mapping of the 225.9-GHz HDO and 203.4-GHz lines on Mars obtained with the IRAM Plateau de Bure facility (PdBI). The observations were performed during martian year 28 (MY28), at L_s=320.3° for the HDO line, and at L_s=324.3° for the line. The HDO line is detected at the eastern (morning) and western (evening) limbs of the northern hemisphere, corresponding to a water column density in the range 3-6 pr.-μm. The line is not detected, which is compatible with the column densities derived from the HDO line. Quasi-simultaneous far infrared measurements obtained by the Planetary Fourier Spectrometer (PFS) onboard the Mars Express spacecraft confirm our PdBI results, yielding a 5±1 pr.-μm meridionally constant water column abundance.Such a low water abundance during the southern mid-autumn of MY28 does not correspond to the standard martian climatology as observed during the previous years. It was however already retrieved from near-infrared observations performed by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) onboard the Mars Reconnaissance Orbiter spacecraft [Smith, M.D., Wolff, M.J., Clancy, R.T., Murchie, S.L. 2009. CRISM observations of water vapor and carbon monoxide. J. Geophys. Res. 114, doi: 10.1029/2008JE003288]. Our observations thus confirm that the planet-encircling dust storm that occurred during MY28 significantly affected the martian water cycle. Our observations also demonstrate the usefulness of interferometric submillimeter observations to survey the martian water cycle from ground-based facilities.     

Comet nucleus size distributions from HST and Keck telescopes

The Wide Field Camera (WFC) on the Hubble Space Telescope and the Low Resolution Imaging Spectrograph (LRIS) on the Keck II telescope have been used to image 21 distant dynamically new, long-period (LP) and short-period (SP) Jupiter-family (JF) comet nuclei (near aphelion), as part of a long-term program to search for physical differences between short-period comets and Oort cloud comets. WFC data were obtained on Comets C/1987 H1 (Shoemaker) and C/1984 K1 (Shoemaker) during Cycle 5 (1995 December) and on C/1988 B1 (Shoemaker), C/1987 F1 (Torres), and C/1983 O1 (?ernis) during Cycle 6 (1997 April, May, and June). The HST comets were at heliocentric distances 20.4 < r[AU] < 29.5. Each comet observation was allocated 7 orbits, for ≈3.6 hrs of integration. The most difficult part of the image reduction was the removal of cosmic rays. We present our scheme for cosmic ray removal. None of the HST comet nuclei was detected to the 3-σ level at m_R∼27. The inferred upper limits to the nucleus radii are . The SP comets range in radius between , with a median value of R_N∼1.61 km. The LP comets ranged in size between <4.0-56 km. Over a range of radii between 1-10 km, the nuclei can be fit with a cumulative distribution N(>R_N)∝R_N^−α with α=1.45±0.05, and for nuclei in the range 2-5 km, α=1.91±0.06. Statistical analysis and modeling shows that the slopes of the observed TNO and JF comet distributions are not compatible, suggesting that the intrinsic distribution of JF comet nuclei is a differential a^−3.5 power law truncated at small nucleus radii between 0.3 and 2.0 km.     

A high-amplitude thermal inertia anomaly of probable magnetospheric origin on Saturn’s moon Mimas

Spectral maps of Mimas’ daytime thermal emission show a previously unobserved thermal anomaly on Mimas’ surface. A sharp V-shaped boundary, centered at 0°N and 180°W, separates relatively warm daytime temperatures from a cooler anomalous region occupying low- to mid-latitudes on the leading hemisphere. Subsequent observations show the anomalous region is also warmer than its surroundings at night, indicating high thermal inertia. Thermal inertia in the anomalous region is , compared to < outside the anomaly. Bolometric Bond albedos are similar between the two regions, in the range 0.49-0.70. The mapped portion of the thermally anomalous region coincides in shape and location to a region of high-energy electron deposition from Saturn’s magnetosphere, which also has unusually high near-UV reflectance. It is therefore likely that high-energy electrons, which penetrate Mimas’ surface to the centimeter depths probed by diurnal temperature variations, also alter the surface texture, dramatically increasing its thermal inertia.     

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

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⁻¹).     

Heavy ion escape from Mars, influence from solar wind conditions and crustal magnetic fields

We have used more than 4 years of Mars Express ion data to estimate the escape of heavy ions ( and ) from Mars. To take the limited field of view of the instrument into account, the data has been binned into spatial bins and angular bins to create average distribution functions for different positions in the near Mars space. The net escape flux for the studied low solar activity period, between May 2007 and May 2011, is 2.0 ± 0.2 × 10²⁴ s⁻¹. The escape has been calculated independently for four different quadrants in the Y_MSO − Z_MSO plane, south, dusk, north and dawn. Escape is highest from the northern and dusk quadrants, 0.6 ± 0.1 × 10²⁴ s⁻¹, and smallest from the south and dawn quadrants, 0.4 ± 0.1 × 10²⁴ s⁻¹. The flux ratio of molecular ( and ) to O⁺ ions is 0.9 ± 0.1, averaged over all quadrants. The flux difference between the north and south quadrants is statistically significant, and is presumed to be due to the presence of significant crustal magnetic fields in the southern hemisphere, reducing the outflow. The difference between the dawn and dusk quadrants is likely due to the magnetic tension associated with the nominal Parker angle spiral, which should lead to higher average magnetic tension on the dusk side. The escape increases during periods of high solar wind flux and during times when co-rotating interaction regions (CIR) affect Mars. In the latter case the increase is a factor 2.4-2.9 as compared to average conditions.     

UV and IR auroral emission model for the outer planets: Jupiter and Saturn comparison

Planetary aurora display the dynamic behavior of the plasma gas surrounding a planet. The outer planetary aurora are most often observed in the ultraviolet (UV) and the infrared (IR) wavelengths. How the emissions in these different wavelengths are connected with the background physical conditions are not yet well understood. Here we investigate the sensitivity of UV and IR emissions to the incident precipitating auroral electrons and the background atmospheric temperature, and compare the results obtained for Jupiter and Saturn. We develop a model which estimates UV and IR emission rates accounting for UV absorption by hydrocarbons, ion chemistry, and non-LTE effects. Parameterization equations are applied to estimate the ionization and excitation profiles in the H₂ atmosphere caused by auroral electron precipitation. The dependences of UV and IR emissions on electron flux are found to be similar at Jupiter and Saturn. However, the dependences of the emissions on electron energy are different at the two planets, especially for low energy (<10 keV) electrons; the UV and IR emissions both decrease with decreasing electron energy, but this effect in the IR is less at Saturn than at Jupiter. The temperature sensitivity of the IR emission is also greater at Saturn than at Jupiter. These dependences are interpreted as results of non-LTE effects on the atmospheric temperature and density profiles. The different dependences of the UV and IR emissions on temperature and electron energy at Saturn may explain the different appearance of polar emissions observed at UV and IR wavelengths, and the differences from those observed at Jupiter. These results lead to the prediction that the differences between the IR and UV aurora at Saturn may be more significant than those at Jupiter. We consider in particular the occurrence of bright polar infrared emissions at Saturn and quantitatively estimate the conditions for such IR-only emissions to appear.     

Modelling meteor ablation in the venusian atmosphere

A comparative study of meteor ablation in the atmospheres of the Earth and Venus is presented. The classical single body meteor ablation model is extended to incorporate a heat penetration depth estimate allowing the simulation of larger meteoroids, than would an isothermal model. The ablation of icy and rocky meteoroids, with densities of 1.0 and 3.4 g cm⁻³, respectively, and initial radii of up to for rock and for ice (equivalent to an initial mass of in both cases), was simulated in both atmospheres. In general venusian meteors are brighter than terrestrial equivalents. Large, slow, rocky objects may be up to 0.7 mag brighter on Venus, while small, icy particles with entry speeds in the range 30-60 km s⁻¹, are found to be upwards of 2.7 mag brighter than at the Earth. Venusian meteors reach maximum brightness at greater altitudes than would similar particles at the Earth. Rocky meteoroids have their points of maximum brightness some 15-35 km higher up at Venus, between 90 and 120 km, whereas, for icy particles this altitude difference is about 5-25 km higher up than at the Earth, in the range 100-125 km. These findings agree, for the most part, with recent analytical studies. Venusian meteors, which last from 100 ms to , tend to be shorter-lived than terrestrial meteors, with correspondingly shorter visible trails. Large (), slow () icy particles reach a maximum magnitude of ∼−2 at Venus and remain visible for about one second, with a large section of the smaller faster meteoroids simulated here remaining visible for several hundred milliseconds. In light of recent space-based meteor observations at the Earth [Jenniskens, P., Tedesco, E., Muthry, J., Laux, C.O., Price, S., 2002. Meteorit. Planet. Sci. 37, 1071-1078], such brightness, height and duration estimates as suggested in this work, may be used in developing future observational campaigns to be carried out from Venus orbit.