Open Access
Volume 642, October 2020
Article Number L11
Number of page(s) 6
Section Letters to the Editor
Published online 09 October 2020

© H. Chen et al. 2020

Licence Creative Commons
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Meteoroids are remnants of the early solar system. When a meteoroid enters Earth’s atmosphere, it generates a luminous trace known as a meteor. Determining meteor trajectories allows us to determine the origin of the meteoroid (e.g., parent body and age, Drummond 1981; Olsson-Steel 1986; Jenniskens 2006; Šegon et al. 2017; Guennoun et al. 2019; Ryabova et al. 2019). In addition, knowledge of meteoroid flux density can constrain environment models, such as the IMEM2, which helps us to plan safe space activities in the entire solar system (Divine 1993; Lemcke 1997; McBride 1997; Soja et al. 2019). Knowledge of the flux density near Earth is especially important for planetary defense. Many ground-based tracking networks (e.g., FRIPON, CAMS, DFN, EDMOND, and AMOS) have been developed with the purpose of determining the flux density (Colas et al. 2014, 2018; Rudawska et al. 2015; Tóth et al. 2015; Day et al. 2016; Jenniskens et al. 2016). Moreover, space-borne observation of meteors has recently been investigated. The space-proof camera SPOSH was developed for meteor detection (Oberst et al. 2011). Two projects have been launched with the purpose of observing meteors from space, namely the SCUBE CubeSat and the METEOR experiment on the International Space Station (ISS; Ishimaru et al. 2014; Arai et al. 2019). The CubeSat project Meteorix is planned for detecting meteors and re-entering space debris (Rambaux et al. 2018, 2019).

As a meteor is a transient event, the reconstruction of its orbit generally requires observations from distributed stations. General methods for meteor triangulation have been well documented (e.g., Wray 1967). Scenarios and algorithms that employ multiple images taken from space only or from the ground only have been intensively studied (Öpik & Singer 1959; Bronshtehn 1981; Ceplecha 1987; Borovicka 1990; Gural 2012; Hajduková 2017; Egal et al. 2017; Vida et al. 2020; Sansom et al. 2019; Chen et al. 2020). However, a combination of space- and ground-based observations has not yet been explored. With the establishment of country- to continental-scale meteor networks and the emergence of space-borne meteor observations, the chance to capture meteors from both the ground and space is increasing with time. Moreover, any Earth observation satellite can participate in meteor observations in conjunction with ground stations; the number of Earth observation satellite is increasing thanks to the establishment of new satellite constellations (e.g., the Starlink project). In particular, the future Meteorix can work with any station of ground networks (e.g., FRIPON, CAMS, and DFN) to enable and enhance meteor orbit determination. The long baseline between the satellite and the ground station seems favorable in terms of stereoscopic measurements. Nevertheless, the positioning accuracy is subject to the geometry formed by the moving target and observers; the setup of the systems, such as the attitude and orbit determination accuracy of the satellite; and image resolution.

This Letter presents an error analysis against the above-mentioned factors, and shows the requirements for system setups and optimal geometry. After briefly recalling the space- and ground-based meteor observations (Sect. 2), we present the principle of stereoscopic measurement and methods of deriving position errors for the scenario of interest (Sect. 3), the results showing the influences of different error sources as well as the necessary system setup (Sect. 4), and a discussion on the impact of the observation geometry (Sect. 5).

2. Basic observation settings

Concerning the space-borne observation, we refer to the design of the Meteorix project. The CubeSat is equipped with a camera and a detection chain that features near real-time detection at low power consumption (Lacassagne et al. 2016). The camera is placed at the end of the 3U CubeSat and pointed toward the nadir. The science orbit is set to a sun-synchronous orbit, which is the most common orbit that maximizes launch opportunities. It is assumed that the orbit is nearly circular with an altitude around 500 km. The period is 94.7 min. The local time of the ascending node is set to 10:30. Consequently, the night time lasts around 32 min per orbit. It has been estimated that this setup allows a detection of around 100 meteors in one year (Rambaux et al. 2019). The selected image sensor has 2048 × 2048 pixels, and the camera has a field of view of 40°. For the ground-based observation, we adopt the specification of a current FRIPON station: field of view of 180° (like many other all-sky networks, such as CAMS and AMOS), resolution of roughly 1000 × 1000 pixels, and axis of the camera pointed toward the zenith.

The usual observable altitude of meteors lies between 80 and 120 km (Ceplecha et al. 1998). The speed of a meteor with respect to the ground is in the range from 11 to 72 km s−1. Considering a frame rate of 10 fps, to capture the motion of the slowest meteors the position accuracy should be better than 1.1 km.

3. Space and ground stereoscopic measurement

3.1. Principle

With only one camera the direction from the camera to the target can be measured, while the distance cannot be. With simultaneous observations by stations distributed in space (e.g., Meteroix CubeSat) and on the ground (e.g., FRIPON), the distance between the target and the two stations as well as the three-dimensional coordinates can be determined, as depicted in Fig. 1. Without loss of generality, the coordinate system is chosen such that the origin, O, coincides with the location of the ground station, the x-axis points to the CubeSat, the y-axis is perpendicular to the x-axis and aligned with the direction to the meteor, and the z-axis completes the right-handed coordinate system. Let [x1,y1,z1]T denote the coordinate of the CubeSat. The true coordinate will be [x1,0,0]T. Let [xm,ym,zm]T denote the coordinate of the targeted meteor, and i, j, and k denote the unit vector along the x-, y-, and z-axis, respectively. The position vector of the meteor is


thumbnail Fig. 1.

Stereoscopic measurement by a CubeSat and a ground station, and the defined coordinate system.

Open with DEXTER

Let φ1 and φ2 represent angles between the directions from the observers to the meteor and the baseline between the two observers, as depicted in Fig. 1. The direction from an observer to the target can be derived from the two-dimensional coordinates on the image. The coordinates of the meteor can be computed from




3.2. Measurement errors

The attitude and orbit determination errors of the CubeSat and image resolutions will translate into errors of coordinates of the meteor. Differentiating Eq. (1) yields


Considering only in-plane (i.e., the x − y plane) errors, orbit determination error δy1 along j results in a change of the reference frame, which is expressed as


where B = x1, representing the length of the baseline between the two cameras. For out-of-plane errors the last term of Eq. (4) is zero as zm is zero. Let δr also be expressed as




Differentiating xm and ym in Eq. (2) with respect to the variables yields


The direction error δφ1 is a result of in-plane image resolution, ; the in-plane attitude determination error, ; and the position error of the CubeSat along the y-axis, which is expressed as


Similarly, δφ2 is extended to


where represents the in-plane direction error due to the image resolution of the ground station.

It is tricky to determine the offset δzm along k when the measured directions from the cameras to the target do not intersect in three-dimensional space due to errors. Let δzI and δzII denote the out-of-plane errors attributed to the CubeSat and the ground station, respectively. The z coordinate error δzm can be approximated by


Here δzI is a result of the out-of-plane image resolution, ; attitude determination error, ; and orbit determination error, δz1, which is expressed as


and δzII, due to the out-of-plane image resolution of the ground station, , is expressed as


It can be reasonably assumed that in-plane and out-of-plane direction errors vary independently. Uncertainties of errors have the relationship,


where σ indicates the 1σ uncertainty of the variable following it, θ1 and θ2 respectively denote the fields of views, and px1 and px2 respectively denote the pixels (in a line) of the CubeSat and the ground station. It is also assumed that distribution of orbit determination error in space is homogeneous:


The position uncertainty σr is defined as the root sum of variances of the three position components, which is expressed as


where σ2 indicates the variance of the variable following it. According to Eqs. (7)–(13) and given that the error sources are uncorrelated, σ2ri, σ2rj, and σ2rk can be computed from




3.3. Reduction of geometries

Figure 2 depicts the geometry of space- and ground-based observations. Altitudes of the CubeSat, hsat, and meteors, hmet, vary within small ranges (around 500 km and between 80 and 120 km, respectively). In the following analyses, hsat is set to 500 km and hmet 100 km to represent a general situation. The observation geometry can then be defined by the directions from the ground station to the CubeSat and the meteor. For convenience, the direction of the CubeSat with respect to the zenith of the ground station, α, and the direction of the observed meteor from the CubeSat with respect to the camera axis, β, are chosen to generalize the observing geometry of interest. As α and β are limited by the fields of view of the two observers, possible observation geometries can thus be represented by an α-β map. Taking into account the round surface of the Earth (where the radius of the Earth R = 6371 km is used), B, φ1, and φ2 entering in Eqs. (17)–(19) are computed as follows. First,


thumbnail Fig. 2.

Simplified geometry for meteor observation.

Open with DEXTER

The triangle formed by the observers and target does not necessarily lie on the vertical plane. The angle φ1 is actually within an interval expressed as

As β is to be varied in the range of the field of view (e.g., −20° to 20° in this work), which can cover situations at two ends of the intervals, the problem can be reasonably simplified by fixing φ1 for a given β to




where the distances from the meteor to the CubeSat, L1, and the ground station, L2, are computed from


4. Results

In order to provide a guideline for system setups, we look at the influences of the orbit determination error (e.g., σx1), attitude determination error σφatt, image error of the CubeSat θ1/px1, and image error of the ground station θ2/px2. The baseline setup assumes σx1 = 10 km, corresponding to coarse orbit determination based on TLE data only (Smith 2019), σφatt = 4°, corresponding to a coarse attitude determination based on magnetometer and sun sensors only (KU Leuven 2018), θ1 = 40°, px1 = 2048, θ2 = 180°, and px2 = 1000, corresponding to the baseline design of Meteorix and FRIPON. Then α ranges from 0 to 90°, and β from −20° to 20°.

The position accuracy σr computed from Eqs. (16)–(19) is presented in Fig. 3 as a function of 90 ° −α (i.e., elevation of the CubeSat) and β. The four panels correspond to the cases considering orbit determination error only (upper left), attitude determination error only (upper right), image error of the CubeSat only (lower left), and image error of the ground station only (lower right). The lower left panel shows that the CubeSat image resolution is good enough to not cause σr >  0.5 km, except in the singular-geometry region (see next section for discussion). The upper right panel shows that coarse attitude estimation results in a σr around 40 km for most observation geometries. Therefore, the CubeSat should perform fine attitude estimation where a star tracker is necessary. In this case, σφatt = 0.04° is expected (KU Leuven 2018), and the contribution of σφatt to σr decreases to 0.4 km. Even with this improvement the resulting σr can still go beyond 15 km, due to the orbit determination error (see the upper left panel). It is recommended to perform fine orbit determination as well. A GPS receiver can enable a fine orbit determination at an accuracy of 0.1 km (Smith 2019), which will reduce the error contribution to 0.15 km for most geometries. The resulting σr can go up to 2 km for some situations, given the effect of the image resolution of the ground station. To widely achieve a desired σr of 1.1 km, it is recommended to increase the image resolution of the ground station, for instance to 2048 × 2048 pixels. With an upgraded resolution, the minimum σr will go down to 1 km for most geometries. The position uncertainty σr following the recommended setups is shown in Fig. 4. The pattern is similar to the lower left panel of Fig. 3, which means the contribution of σφ2 due to the image resolution of ground station is dominant.

thumbnail Fig. 3.

Positioning uncertainty (color levels) contributed by baseline orbit and attitude determination errors (upper panels) of the CubeSat, and baseline image resolution of the CubeSat and the ground station (lower panels).

Open with DEXTER

5. Discussion

Figure 4 also displays the elevation of the meteor with respect to the ground station, which is derived from α and β. The shaded area where the meteor is below the horizon is not a valid area for observation, and thus is excluded from our discussion. In the figure, an elevation > 90° (which is not a common usage) indicates that the meteor and the CubeSat are on different sides of the zenith of the ground station.

thumbnail Fig. 4.

Positioning uncertainty (color levels) for the recommended system setup, and elevation (gray lines) of the meteor with respect to the ground station.

Open with DEXTER

Figure 4 shows that σr strongly depends on the elevation of the meteor. Following the recommended setup, σr <  2 km is widely achieved except for the elevation interval [65°,90°] (labeled case 2 in the figure). The corresponding geometry is referred to as the singular geometry. When the meteor is at an elevation of around 30° (labeled case 1) and also captured by the CubeSat camera, the geometry is optimal with σr ranging from 0.5 to 0.8 km, which depends on the observing directions from the CubeSat. Another optimal region appears in the elevation interval [110°,140°] (labeled case 3), which corresponds to the commonly used elevation [40°,70°] with the CubeSat observing from the other side of the zenith.

Figure 5 schematically visualizes the singular and optimal geometries revealed in Fig. 4. With only one camera, the error (green band) resulting from the angular size of a pixel is infinite along the observing direction (green dotted line), which should be constrained by the observation from another direction. If two directions are aligned, the error is still infinite. Because of the limited field of view of a CubeSat, it will be roughly on top of the target when capturing the target. When the meteor is near the zenith of the ground station (case 2 in Fig. 5) and the CubeSat observes it from the top, the error along the observing direction is not effectively constrained. Therefore, the resulting radius of the error zone (i.e., the intersecting green zone) is large. This explains the singularity associated with the meteor elevation interval [65°,90°]. When the elevation of the meteor is low (30° in case 1 and 40° in case 3), the included angle between the two observing directions is close to 90°, leading to effective triangulation and small errors. Even though the error across the observing direction is proportional to the square of the distance, and thus large for low-elevation observations, the resulting range of errors, which is mainly contributed by the error along the observing direction, is still minimized for low elevations.

thumbnail Fig. 5.

Intersecting error zones depending on the geometry.

Open with DEXTER

6. Conclusions

To achieve a desirable positioning accuracy of 1 km, the CubeSat should bring high-accuracy sensors (e.g., star tracker and GPS receiver) to perform fine attitude (i.e., accuracy ∼0.05°) and orbit (i.e., accuracy ∼100 m) determination. In addition, the image resolution of FRIPON is recommended to upgrade to 2048 × 2048 pixels. The performance of meteor triangulation is best (i.e., with smallest position error) when the meteor elevation is around 30°. With the recommended system setup, as long as a meteor is captured by both the CubeSat and the ground station, and at an elevation < 65° with respect to the ground station, the position and the orbit of the meteor can be determined.


We acknowledge the ESEP (Exploration Spatiale des Environnements Planétaires) postdoctoral fellowship, DIM ACAV+ Région Île-de-France, Janus CNES, and IDEX Sorbonne Universités for funding this research.


  1. Arai, T., Kobayashi, M., Yamada, M., et al. 2019, 82nd Annual Meeting of the Meteoritical Society, 82, 6493 [Google Scholar]
  2. Borovicka, J. 1990, Bull. Astron. Inst. Czechoslov., 41, 391 [Google Scholar]
  3. Bronshtehn, V. A. 1981, Physics of Meteor Phenomena (Moskva: Nauka) [Google Scholar]
  4. Ceplecha, Z. 1987, Bull. Astron. Inst. Czechoslov., 38, 222 [Google Scholar]
  5. Ceplecha, Z., Borovička, J., Elford, W. G., et al. 1998, Space Sci. Rev., 84, 327 [NASA ADS] [CrossRef] [Google Scholar]
  6. Chen, H., Rambaux, N., Matha, R., & Chelil, R. 2020, 5th IAA Conference on University Satellite Missions, Rome, Italy [Google Scholar]
  7. Colas, F., Zanda, B., Bouley, S., et al. 2014, in Asteroids, Comets, Meteors 2014, eds. K. Muinonen, A. Penttilä, M. Granvik, et al., 110 [Google Scholar]
  8. Colas, F., Baillié, K., Bouley, S., et al. 2018, European Planetary Science Congress, EPSC2018-1105 [Google Scholar]
  9. Day, B., Bland, P., Al-Rajeh, Y., & Schmidt, G. 2016, Soc. Astron. Sci. Annu. Symp., 35, 165 [Google Scholar]
  10. Divine, N. 1993, J. Geophys. Res., 98, 17029 [NASA ADS] [CrossRef] [Google Scholar]
  11. Drummond, J. D. 1981, Icarus, 45, 545 [NASA ADS] [CrossRef] [Google Scholar]
  12. Egal, A., Gural, P. S., Vaubaillon, J., Colas, F., & Thuillot, W. 2017, Icarus, 294, 43 [NASA ADS] [CrossRef] [Google Scholar]
  13. Guennoun, M., Vaubaillon, J., Čapek, D., Koten, P., & Benkhaldoun, Z. 2019, A&A, 622, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Gural, P. S. 2012, Meteorit. Planet. Sci., 47, 1405 [NASA ADS] [CrossRef] [Google Scholar]
  15. Hajduková, M. Jr., Koten, P., Kornoš, L., & Tóth, J. 2017, Planet. Space Sci., 143, 89 [CrossRef] [Google Scholar]
  16. Ishimaru, R., Sakamoto, Y., Kobayashi, M., et al. 2014, Lunar Planet. Sci. Conf., 45, 1846 [Google Scholar]
  17. Jenniskens, P. 2006, Meteor Showers and their Parent Comets (Cambridge, UK: Cambridge University Press) [CrossRef] [Google Scholar]
  18. Jenniskens, P., Nénon, Q., Gural, P. S., et al. 2016, Icarus, 266, 384 [NASA ADS] [CrossRef] [Google Scholar]
  19. KU Leuven 2018, Accurate and Agile Attitude Determination and Control for CubeSats [Google Scholar]
  20. Lacassagne, L., Cabaret, L., Etiemble, D., Hebache, F., & Petreto, A. 2016, Proceedings of the 3rd Workshop on Programming Models for SIMD/Vector Processing, WPMVP’16 (New York, NY, USA: Association for Computing Machinery) [Google Scholar]
  21. Lemcke, C. 1997, Second European Conference on Space Debris, ESA SP-393 (Darmstadt, Germany: ESA, ESOC), 649 [Google Scholar]
  22. McBride, N. 1997, Adv. Space Res., 20, 1513 [NASA ADS] [CrossRef] [Google Scholar]
  23. Oberst, J., Flohrer, J., Elgner, S., et al. 2011, Planet. Space Sci., 59, 1 [NASA ADS] [CrossRef] [Google Scholar]
  24. Olsson-Steel, D. 1986, MNRAS, 219, 47 [CrossRef] [Google Scholar]
  25. Öpik, E. J., & Singer, S. F. 1959, Phys. Today, 12, 54 [Google Scholar]
  26. Rambaux, N., Galayko, D., Guignan, G., et al. 2018, COSPAR 2018, 42nd Assembly, Pasadena, United States [Google Scholar]
  27. Rambaux, N., Vaubaillon, J., Lacassagne, L., et al. 2019, 1st ESA NEO and Debris Detection Conference (Darmstadt, Germany: ESA Space Safety Programme Office), 9 [Google Scholar]
  28. Rudawska, R., Matlovič, P., Tóth, J., & Kornoš, L. 2015, Planet. Space Sci., 118, 38 [NASA ADS] [CrossRef] [Google Scholar]
  29. Ryabova, G. O., Asher, D. J., & Campbell-Brown, M. J. 2019, Meteoroids: Sources of Meteors on Earth and Beyond (Cambridge: Cambridge University Press) [CrossRef] [Google Scholar]
  30. Sansom, E. K., Jansen-Sturgeon, T., Rutten, M. G., et al. 2019, Icarus, 321, 388 [NASA ADS] [CrossRef] [Google Scholar]
  31. Smith, A. 2019, SpaceOps Workshop, Santa Clara, United States [Google Scholar]
  32. Soja, R. H., Grün, E., Strub, P., et al. 2019, A&A, 628, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Tóth, J., Kornoš, L., Zigo, P., et al. 2015, Planet. Space Sci., 118, 102 [NASA ADS] [CrossRef] [Google Scholar]
  34. Vida, D., Gural, P. S., Brown, P. G., Campbell-Brown, M., & Wiegert, P. 2020, MNRAS, 491, 2688 [NASA ADS] [CrossRef] [Google Scholar]
  35. Šegon, D., Vaubaillon, J., Gural, P. S., et al. 2017, A&A, 598, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Wray, J. D. 1967, The Computation of Orbits of Doubly Photographed Meteors (Albquequerque: University of New Mexico Press) [Google Scholar]

All Figures

thumbnail Fig. 1.

Stereoscopic measurement by a CubeSat and a ground station, and the defined coordinate system.

Open with DEXTER
In the text
thumbnail Fig. 2.

Simplified geometry for meteor observation.

Open with DEXTER
In the text
thumbnail Fig. 3.

Positioning uncertainty (color levels) contributed by baseline orbit and attitude determination errors (upper panels) of the CubeSat, and baseline image resolution of the CubeSat and the ground station (lower panels).

Open with DEXTER
In the text
thumbnail Fig. 4.

Positioning uncertainty (color levels) for the recommended system setup, and elevation (gray lines) of the meteor with respect to the ground station.

Open with DEXTER
In the text
thumbnail Fig. 5.

Intersecting error zones depending on the geometry.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.