FRIPON: A worldwide network to track incoming meteoroids

F. Colas1, 9 , B. Zanda2, 1, 9, S. Bouley3, 1, 9, S. Jeanne1, 9, A. Malgoyre7, 9, M. Birlan1, 9, C. Blanpain7, 9, J. Gattacceca5, 9, L. Jorda4, 9, J. Lecubin7, 9, C. Marmo3, J.L. Rault1, 9, 53, J. Vaubaillon1, 9, P. Vernazza4, 9, C. Yohia7, 9, D. Gardiol10, A. Nedelcu157, 37, B. Poppe18, 40, J. Rowe45, M. Forcier16, 17, D. Koschny35, 51, 200, J.M. Trigo-Rodriguez34, 231, 232, H. Lamy33, 137, R. Behrend66, 41, L. Ferrière32, 114, D. Barghini10, 11, A. Buzzoni12, A. Carbognani12, M. Di Carlo26, M. Di Martino10, C. Knapic13, E. Londero13, G. Pratesi14, 22, S. Rasetti10, W. Riva15, G.M. Stirpe12, G.B. Valsecchi22, 23, C.A. Volpicelli10, S. Zorba13 D. Coward271, 272, E. Drolshagen18, 40, G. Drolshagen18, 40, O. Hernandez16, 17, E. Jehin33, 134, M. Jobin16, 17, A. King191, 45, 46, C. Nitschelm31, 156, T. Ott18, 40, A. Sanchez-Lavega19, 20, A. Toni35, 51, P. Abraham54, F. Affaticati187, M. Albani187, A. Andreis188, T. Andrieu219, S. Anghel37, 74, 157, E. Antaluca55, K. Antier9, 44, 53, T. Appéré56, A. Armand117, G. Ascione199, Y. Audureau3, G. Auxepaules57, T. Avoscan198, D. Baba Aissa60, 206, P. Bacci245, O. Bǎdescu37, 157, R. Baldini246, R. Baldo58, A. Balestrero15, D. Baratoux59, E. Barbotin265, M. Bardy61, S. Basso30, O. Bautista62, L. D. Bayle63, P. Beck64, 65, R. Bellitto253, R. Belluso27, C. Benna10, M. Benammi67, 68, E. Beneteau145, Z. Benkhaldoun38, 69, P. Bergamini70, F. Bernardi247, M.E. Bertaina11, P. Bessin154, L. Betti234, F. Bettonvil50, 35, D. Bihel71, C. Birnbaum9, 43, O. Blagoi157, 37, E. Blouri9, 214, I. Boacă157, 37, R. Boatǎ211, 37, B. Bobiet72, R. Bonino11, K. Boros235, E. Bouchet196, 41, V. Borgeot131, E. Bouchez73, D. Boust75, V. Boudon76, T. Bouman77, P. Bourget78, 31, S. Brandenburg49, 35, Ph. Bramond79, E. Braun80, A. Bussi198, P. Cacault81, B. Caillier82, A. Calegaro137, 33, J. Camargo83, 39, S. Caminade8, A.P.C. Campana84, P. Campbell-Burns45, R. Canal-Domingo168, 34, O. Carell71, S. Carreau85, E. Cascone237, C. Cattaneo248, P. Cauhape129, P. Cavier86, S. Celestin87, A. Cellino10, M. Champenois89, H. Chennaoui Aoudjehane92, 69, S. Chevrier87, P. Cholvy139, L. Chomier90, A. Christou91, 45, D. Cricchio238, P. Coadou103, J.Y. Cocaign95, 223, F. Cochard93, S. Cointin94, E.


Abstract
Context: Until recently, camera networks designed for monitoring fireballs worldwide were not fully automated, implying that in case of a meteorite fall, the recovery campaign was rarely immediate. This was an important limiting factor as the most fragile -hence precious -meteorites must be recovered rapidly to avoid their alteration.

Aims: The Fireball Recovery and InterPlanetary Observation Network (FRIPON) scientific
project was designed to overcome this limitation. This network comprises a fully automated camera and radio network deployed over a significant fraction of western Europe and a small fraction of Canada. As of today, it consists of 150 cameras and 25 European radio receivers and covers an area of about 1.5 × 10 6 km 2 .
Methods: The FRIPON network, fully operational since 2018, has been monitoring meteoroid entries since 2016, thereby allowing the characterization of their dynamical and physical properties.
In addition, the level of automation of the network makes it possible to trigger a meteorite recovery campaign only a few hours after it reaches the surface of the Earth. Recovery campaigns are only organized for meteorites with final masses estimated of at least 500 g, which is about one event per year in France. No recovery campaign is organized in the case of smaller final masses on the order of 50 g to 100 g, which happens about three times a year; instead, the information is delivered to the local media so that it can reach the inhabitants living in the vicinity of the fall.
Results: Nearly 4,000 meteoroids have been detected so far and characterized by FRIPON.
The distribution of their orbits appears to be bimodal, with a cometary population and a main belt population. Sporadic meteors amount to about 55% of all meteors. A first estimate of the absolute meteoroid flux (mag < -5; meteoroid size ≥∼1 cm) amounts to 1,250/year/10 6 km 2 . This value is compatible with previous estimates. Finally, the first meteorite was recovered in Italy (Cavezzo, January 2020) thanks to the PRISMA network, a component of the FRIPON science project.

Introduction
The study of the physical and dynamical properties of interplanetary matter, such as interplanetary dust particles (IDPs), meteoroids, asteroids, comets, is crucial to our understanding of the formation and evolution of the solar system. This matter exists in many sizes, from micron-sized dust grains to several hundred kilometer-sized bodies. Whereas the largest bodies are routinely studied via Article number, page 12 of 54 Fig. 1. In nineteenth century France, 45 meteorites were recovered after their fall was observed, a number that fell by a factor of 5 in the twentieth century. Even in the nineteenth century, witnessed falls were not randomly distributed. They were mostly located in the great river plains (Seine and Loire in the northwest, Garonne in the southwest, and Rhône valley in the southeast). In these regions, the population was denser, the view is free of obstacles (such as mountains), and the skies are often clear. The striking difference between the two centuries illustrates the need for distributed observers for meteorite recovery. Rural populations have declined because of urbanization in the twentieth century. A camera network such as FRIPON can monitor atmospheric entries and take over that role that was previously played by human observers. However, trained human eyes are still required to recover the meteorites; this is the aim of the Vigie-Ciel citizen science program (Colas et al. 2015).
Earth-based telescopic observations as well as less frequent interplanetary missions, the smallest bodies (diameter ≤10 m) are for the most part only observed and characterized when they enter the Earth's atmosphere as their entry generates enough light to be recorded by even the simplest types of cameras; the smaller particles are called meteors and the larger bodies are fireballs.
We know that ∼100 tons of extraterrestrial material collide with the Earth daily, mostly as small particles less than 0.2 mm in size (Zolensky et al. 2006, Rojas et al. 2019. At present, these small particles, called IDPs, are actively being collected in the stratosphere, from polar ices (Duprat et al. 2007), and within impact features on spacecraft (Moorhead et al. 2020). For such particles, the stratospheric collections provide the least contaminated and heated samples. At the other end of the size distribution of extraterrestrial material colliding with the Earth, meteorites are fragments that have survived the passage through the atmosphere without internal chemical alteration, which have been recovered at the surface of the Earth. To date, all known meteorites are pieces of either asteroids, the Moon, or Mars, with asteroidal fragments dominating the flux of material, whereas IDPs originate mostly from comets and possibly from asteroids (Bradley et al. 1996;Vernazza et al. 2015). The most detailed information on the processes, conditions, timescales, and chronology of the early history of the solar system (e.g., Neveu & Vernazza 2019; Kruijer & Kleine 2019 and references therein), including the nature and evolution of the particles in the pre-planetary solar nebula, has so far come from the study of all these extraterrestrial materials. Recovering intact samples of such materials is therefore a critical goal of planetary studies.
However, we are not very efficient at recovering the meteorites that hit the Earth. Estimates based on previous surveys (Bland et al. 1996) and on collected falls [Meteoritical Bulletin database 1 ] indicate that, for meteorites with masses greater than 100 g, probably less than 1 in 500 that fall on Earth are currently recovered. In addition, taking France as an example, recovery rates were significantly higher in the nineteenth century than they are now: 45 meteorites were observed to fall and found on the ground in the nineteenth century, whereas they were 5 times fewer in the twentieth century (Fig. 1), showing that there is at present a large potential for improvement. Hot and cold deserts are privileged dense collection areas, but most meteorites are found hundreds and up to millions of years after their fall (Hutzler et al. 2016;Drouard et al. 2019). They have thus been exposed to terrestrial alteration, which has partly obliterated the scientific information they contain. Also, the critical information regarding their pre-atmospheric orbit is no longer available.
The most efficient approach for recovering freshly fallen meteorites is to witness their bright atmospheric entry via dense (60-120 km spacing) camera and radio networks. These networks make it possible to accurately calculate their trajectory from which both their pre-atmospheric orbit and their fall location (with an accuracy on the order of a few hundred meters) can be constrained.
Records of incoming meteorites started with the appearance of photographic plates at the end of the twentieth century. A first attempt to observe incoming bolides was made in the United States and consisted of a small camera network that was operated between 1936 and 1951 (Whipple 1938), but it was only in the middle of the twentieth century that the first fireball observation networks were developed with the aim of recovering meteorites. Two such networks were established in the 1960s. The first was the Prairie Network (McCrosky & Boeschenstein 1965) in the center of the United States, which remained operational from 1964 to 1975. This network comprised 16 stations located 250 km apart. Only one meteorite was recovered thanks to this network (Lost City, 1970;McCrosky et al. 1971). The low efficiency of the Prairie Network, despite the large area it covered (750, 000 km 2 ) mainly resulted from the low efficiency of the photographic plates, the large distance between the stations, and the slow pace of the data reduction process.
The European Fireball Network (EFN) was also developed in the 1960s, under the guidance of the Ondrejov Observatory, following the recovery of the Příbram meteorite in 1959 (Ceplecha 1960). It is still active, currently covers 1 × 10 6 km 2 with about 40 cameras, (Oberst et al. 1998) and benefits from modern equipment. So far, this network has enabled the recovery of nine meteorites (Table 1).
In 1971, the Meteorite Observation and Recovery Project (MORP) project was established over part of Canada and led to the recovery of the Innisfree meteorite (Halliday et al. 1978). The modern digital camera extension of this network, called the Southern Ontario Meteor Network, led to the recovery of the Grimsby meteorite (Brown et al. 2011). The MORP project comprises 16 cameras and covers a surface area of 700, 000 km 2 . Other networks using photographic techniques have also been developed, such as the Tajikistan Fireball Network (Kokhirova et al. 2015), which consists of 5 cameras and covers 11 000 km 2 . However, none of these other networks have made it possible to recover meteorites so far. We note the existence of other networks such as the SPMN network, which facilitated the recovery of the Villalbeto de la Peña (Trigo-Rodríguez et al. 2006) and Puerto Article number, page 14 of 54 Lápice ) meteorites, as well as the Finnish Fireball Network, which facilitated recovering the Annama meteorite Trigo-Rodríguez et al. 2015). Last, the Desert Fireball Network ) was implemented in Australia in 2007. This network is based on high-resolution digital cameras and has made it possible to recover four meteorites: Bunburra Rockhole in 2007 , Mason Gully in 2010 (Dyl et al. 2016), Murrili in 2015 , and Dingle Dell in 2016 (Devillepoix et al. 2018). The success of this network results from the efficiency of the cameras and the size of the network as well as an efficient data reduction and analysis process (Sansom et al. 2019a). A method to construct a successful fireball network is discussed in Howie et al. (2017).
As of today, there are 38 meteorites with reliable reconstructed orbits, 22 of which were detected by camera networks (see Table 1). Among the remaining 16 meteorites, 14 are the result of random visual observations such as the Chelyabinsk event (data from security cameras were used for orbit computation; Borovička et al. 2013a) and two meteorites were detected as asteroids before their fall (Almahata Sitta and 2018LA). During the same time interval , 397 meteorites were recovered after their falls were witnessed by eye (Meteoritical Bulletin Database).
The main limitation of current networks is their size. Most of these networks consist of a fairly small number of cameras spread over a comparatively small territory. Altogether, they cover only 2% of the total surface of the Earth (Devillepoix et al. 2020). This implies that the number of bright events per year witnessed by these networks is small and that decades would be necessary to yield a significant number (≥100) of samples.
The Fireball Recovery and InterPlanetary Observation Network (FRIPON) scientific project was designed to contribute to this global effort to recover fresh meteorites. It comprises a network deployed over a large fraction of western Europe and a small fraction of Canada (see Fig. 2). As of today, this network consists of 150 cameras and 25 receivers for radio detection and covers an area of 1.5 × 10 6 km 2 (section 3). The FRIPON network is coupled in France with the Vigie-Ciel citizen science program, the aim of which is to involve the general public in the search for meteorites in order to improve their recovery rate. In the present paper, we first describe the technology of the FRIPON network and its architecture, and finally we give the first results obtained after four years of observations and report on the first meteorite recovery in Italy 2 ).

General description of the network
The FRIPON science project was originally designed by a core team of six French scientists from the Paris Observatory (IMCCE), the French National Museum of Natural History (MNHN-IMPMC), Université Paris-Saclay (GEOPS), and Aix-Marseille University (LAM / CEREGE / OSU Pythéas) to: i) monitor the atmospheric entry of fireballs, that is, interplanetary matter with typical sizes greater than ∼1 centimeter; ii) characterize their orbital properties to constrain both Table 1. Thirty-eight known meteorites with reliable orbit reference discovered by networks ("N"), visual observations ("V") or telescopic observations ("T"). their origin and fall location; and iii) recover freshly fallen meteorites. This project benefited from a grant from the French National research agency (Agence Nationale de la Recherche: ANR) in 2013 to install a network of charged coupled device (CCD) cameras and radio receivers to cover the entire French territory. Specifically, the grant was used to design the hardware (section 3.2), building on experience gained from previous networks; develop an efficient and automatic detection and Article number, page 16 of 54 Tunisia. As of today, 150 cameras, using FRIPON technology, and 25 radio receivers are operational around the world (see Fig. 2).
The FRIPON science project regroups all the above-mentioned national networks, with all the cameras monitored and remotely controlled by the Service Informatique Pythéas (SIP; Aix-Marseille University, France), which maintains the whole network with the support of the scientific team. All the data from the FRIPON network are stored and processed in Marseille. The data processing consists of monthly astrometric and photometric reduction of the calibration images and daily processing of multi-detections. Two databases host the data. One stores the raw data and the other stores higher-level, processed data, such as orbits and trajectories. These data are available to all coinvestigators of the network 4 . On request, national data can be sent to a different reduction pipeline for alternate processing and storage 5 .

Optical cameras
Since the early 2000s, digital cameras have been used by all networks that are deployed to monitor fireballs. Two alternate technical solutions are adopted. The first is based on a low-resolution detector (e.g., Southern Ontario Meteor Network; Brown et al. 2011), while the second relies on a high-resolution detector (e.g., Desert Fireball Network; Bland et al. 2012). The measurements acquired by low-resolution cameras can be accurate enough to compute orbits and strewn fields as long as the network is dense, with numerous cameras. For example, the Southern Ontario Meteor Network, which has been operating in Canada since 2004, led to the recovery of the Grimsby meteorite (Brown et al. 2011). In the case of the FRIPON network, we followed the philosophy of the Canadian Fireball Network (Brown et al. 2011) as detailed hereafter.
We used a CCD Sony ICX445 chip with 1296x964 pixels and a pixel size of 3.75 x 3.75 microns. For the optical design, we used a 1.25 mm focal length F/2 fish-eye camera lens, which leads to a pixel scale of 10 arcmin. Given that fireballs are typically observed at an altitude between 100 km and 40 km, we designed a network with a median distance of 80 km between cameras to perform an optimal triangulation. Jeanne et al. (2019) showed that the astrometric accuracy is on the order of 1 arcmin, equivalent to 30 m at a distance of 100 km. In section 4, we show that the final accuracy on the trajectory is on the order of 20 m for the position and of 100 m/s for the velocity; this value is required for the identification of meteorite source regions in the solar system as shown by Granvik & Brown (2018).
The optical device and the CCD were embedded into a special case (Fig. 3) sealed with a transparent dome, thereby allowing us to record full-sky images. Moreover, these cases are equipped with a passive radiator, which serves to release the heat produced by the electronics during the warm periods of the year to minimize CCD dark current.
Each camera is controlled by an Intel NUCi3 computer on which the data are temporarily stored. A single power over ethernet (PoE) cable is used for data transfer and for powering and remotely managing the camera through a TPLINK (TL-SG22110P or 1500G-10PS) switch. Such a solution makes it easy to install the optical station and operate it remotely and to use cables up to 100 meters long between the camera and the computer. Fig. 3 shows the design 6 of the camera as well as its installation at the Pic du Midi Observatory.

Radio receivers
In addition to optical observations, we used the powerful signal of the GRAVES radar of the French Air Force. This radar is particularly well adapted for the detection, identification, and tracking of space targets including incoming meteoroids (Michal et al. 2005). Located near Dijon (Burgundy, central eastern France), its four main beams transmit nominally on a half-volume located south of a line between Austria and western France. However, the secondary radiation lobes of the radar make it possible to also detect meteors that disintegrate in the northern part of France. For such observations we do not need as tight a mesh as we do for the optical network. We have 25 stations with an average distance of 200 km, mainly in France, but also in Belgium, United Kingdom, Italy, Switzerland, Spain, and Austria. The GRAVES radar system transmits on 143.050 MHz in a continuous wave (CW) mode 24 hours a day. A meteoroid entering the E and D layers of the Earth ionosphere produces ions and free electrons generated by the ionization of air and of meteoroid molecules. The free electrons have the property of scattering radio waves according to "back or forward meteor scatter" modes when they are illuminated by a radio transmitter. The FRIPON radio setup is presented in the Appendix.

Data storage and access
The FRIPON stations are composed of a Linux minicomputer, a wide-angle camera, and a manageable switch guaranteeing the isolation of the network of the host institute. The installation is done with an automated deployment system based on a USB key.  (Audureau et al. 2014) and a set of scripts. The data, which include astrometric long exposures images, single detection (stacked images), and multiple detections (both optical and radio raw data) are subsequently transferred to the central server.
The data collected on the server are then indexed in a database. During this operation, visuals are generated. When an optical event groups at least two stations, the FRIPON pipeline is executed to generate the dynamical and physical properties of the incoming meteoroid such as its orbit, its mass and its impact zone.
All the data are made available through a web interface that is accessible to the worldwide community in real time 7 . This interface makes it possible to display and download data in the form of an archive that complies with the data policy of the project by means of access right management.

Detection strategy
The acquisition and detection software FreeTure was specifically developed by the FRIPON team and runs permanently on the minicomputers (see Audureau et al. 2014 for a full description). The images corresponding to single detections by FreeTure are stored locally and a warning (time and location) is sent to the central server in Marseille. If at least one other station detects an event within +/-3 seconds, it is then treated as a "multiple detection". We note that we implemented a distance criterion of less than 190 km to avoid false detections. This value was determined empirically by manually checking one year of double detections. This strategy works well during the night, but leads to 30% of false detections mainly during twilight.
Radio data corresponding to the last week of acquisition are only stored locally. Only radio data acquired at the time of an optical multi-detection are uploaded from the radio stations to the Marseille data center for processing.

Optical data
Scientific optical data are CCD observations recorded at a rate of 30 frames per second (fps). This acquisition rate is necessary to avoid excessive elongation of the meteor in the images in the case of high speed fireballs. For example, a typical bolide with an average speed of 40 km/s at 100 km altitude at the zenith leads to a 20 • /s apparent speed on the sky and to a four pixel elongated trail on the CCD. It is larger than the average width of the point spread function (PSF; typically 1.8 pixels), but still easy to process for centroid determination. No dark and flatfield corrections are made.
However, almost no reference star is measurable on a single frame with such an acquisition speed, as the limiting magnitude is about zero. It is thus necessary to record images with a longer exposure time for calibration. We therefore recorded five second exposure images every ten minutes; the goal is to have a decent signal-to-noise ratio (S/N) up to a magnitude of 4.5 and to only marginally affect detection efficiency. Such a calibration strategy allows the detection of a few thousand calibration stars for a given camera on a clear night. To mitigate the effect of cloudy nights and breakdowns, we computed an astrometric calibration once per month for each station.
This works for most cameras as their mounts are rigid. However, we occasionally detected flexible mounts based on the repeated calibrations, which led us to shorten the masts of such stations.
Calibration procedure uses the ICRF2 8 reference frame. The distortion function of the optical system is computed in the topocentric horizontal reference system. This allows for an astrometric solution for stars above 10 degrees of elevation with an accuracy of 1 arcmin. Our procedure leads to the calculation of the azimuth and the elevation of the bolides in the J2000 reference frame.
More details regarding our astrometric calibration procedure can be found in Jeanne et al. (2019).
For the photometric reduction, we used the same frames as for the astrometric calibration, namely the long exposure frames. We then established a correspondence between the observed stars and those present in the Hipparcos catalog (Bessell 2000). The following steps are subsequently applied to calculate the absolute magnitude light curve of a meteor, namely: i) determination of the flux of an equivalent magnitude 0 star at zenith and the linear extinction function of the air mass for one-month cumulative observation; ii) measurement of the bolide flux on individual frames and conversion in magnitude; and iii) conversion of the meteor magnitude Mag into an absolute magnitude AMag, defined as its magnitude at a distance of 100 km, .
(1)  Bottom: Absolute magnitude light curves of the event as seen by 15 cameras; the red curve is Beaumont-lès-Valence. It is clear that the saturation limit is around magnitude -8 (all the other light curves fall above this limit). Cameras located further away may be able to measure more non-saturated data, but all the cameras become heavily saturated as the bolide reaches its maximum luminosity.
applied in the future. At this point, we point out several limitations of our data reduction procedure as follows: -Clouds may partly cover the night sky, which may bias the measure of instrumental magnitudes.
Article number, page 22 of 54 -Meteors are mainly detected at small elevations (typically below 30 • ). These records are therefore affected by nonlinearities of the atmospheric extinction.
-A uniform cloud layer can be the source of an under estimation of bolide magnitude.
The first photometric measurements of the FRIPON network are reflected in the histogram of all detections in subsection 4.1.2. Routines to merge all light curves into one are now under development. As our data reduction is based on dynamics, the photometric curves are only used at present to detect major events.
To summarize, the astrometric reduction allows us to obtain an accuracy of one-tenth pixel or 1 arcmin for meteor measurements. Photometry is at that time only usable for events with an absolute magnitude lower than −8 with an accuracy of 0.5 magnitude.

Trajectory determination
Most of our method is described in Jeanne et al. (2019) and in Jeanne (2020)  Therefore, we developed a modified least-squares regression to fit the data taking into account the internal and external or systematic error on each camera.
We first estimate the internal error of each camera by fitting a plane passing through the observation station and all the measured points. The average internal error of the cameras amounts to 0.75 arcmin, which corresponds to 0.07 pixel. We also compute a first estimation of the external error by averaging distances between the observed positions of stars and those calculated from the Hipparcos catalog (Bessell 2000) in a neighborhood of 100 pixels around the meteor. We then compute a global solution using the modified least-squares estimator of the trajectory T χ 2 given by the minimization of the following sum: Article number, page 23 of 54 where i j (T ) is the residual between the j th measure taken by the i th camera and the trajectory T , σ i is the internal error of the i th camera, s i is the systematic error of the i th camera, and n i is the number of images taken by the i th camera.
This method allows us to characterize the systematic errors of our cameras (e.g., a misaligned lens), but not errors such as the location of the camera. To tackle these errors, we compute a first estimate of the trajectory and we compare the residuals with the expected random and systematic errors. If they are larger than expected for a specific camera, we iteratively decrease its weight during the calculation of the trajectory. The final systematic error is usually on the order of 0.3 arcmin, which ends the iterative process.
Two geometric configurations lead to important errors or degeneracies in the trajectory determination: stations located too far from the fireball and stations aligned with the trajectory of the fireball. However, most of the time, the final bright flight straight line trajectory is known with a precision of a few tens of meters. In a second step, all individual data points with time stamps are projected on the straight line to be used afterward for dynamical purposes.

Orbit, drag, and ablation model
To compute the orbit of the bolide parent body, we need to measure its velocity before it has experienced significant interaction with the upper atmosphere. This interaction starts well before the bright flight. Therefore, we need a deceleration model to estimate the infinite velocity, even if the deceleration is not measurable, which happens to be the case for many events (especially the high speed events). This problem is complex because physical parameters evolve during atmospheric entry and moreover several parameters are unknown such as drag coefficient, object size, shape, density and strength. Like other teams (Lyytinen & Gritsevich 2016, Bouquet et al. 2014, Sansom et al. 2019b..) we use a simple physical model to fit the bright flight data.
We used a dynamic model from Bronshten (1983), equations (3) and (4). This model describes the deceleration and ablation of a meteoroid in an atmosphere based on the following three equations : where c d is the drag coefficient, c h the heat-transfer coefficient, H is the enthalpy of destruction, ρ atm is the gas density, m is the normalized meteoroid mass, M e is the pre-entry mass, s is the normalized cross-section area, S e is the pre-entry cross-section area, µ is the so-called shape change coefficient.
Article number, page 24 of 54 The atmospheric gas density ρ atm is taken from the empirical model NRLMSISE-00 (Lyytinen & Gritsevich 2016).
These three equations can be rewritten into two independent equations (Turchak & Gritsevich 2014). The equation of motion is written as and the equation of mass is written as where A is a deceleration parameter (in square meters per kilogram) and B is an ablation parameter (in square meters per joule) as follows: We

Dark flight
At the end of the bright flight, a meteoroid is subject only to aerodynamic drag (including winds) and gravity. At this stage, the meteoroid speed is too low to cause ablation (hence dark flight).
The equation of motion during dark flight is as follows: where A f (V w ) is the deceleration parameter of the fragment, which depends on the wind velocity (relative to the fragment) V w . We used a local atmospheric model of wind retrieved from meteorological offices.
Article number, page 25 of 54 The end of the bright flight simulation gives us the initial conditions of the dark flight motion, namely the initial position, speed, and acceleration of the fragment. The initial condition of acceleration gives us a definition of A f 0 , the limit of A f when wind velocity is huge in front of sound velocity c s as follows: The evolution of A f as a function of wind velocity can be retrieved in Ceplecha (1987). Finally, we performed several computations using the Monte Carlo method to take into account the measurement errors of all the initial parameters to obtain a ground map (strewn field) as a function of the final mass of the bolide.
Of course, owing to the various simplifying assumptions made, we can only underestimate the size of the strewn field. However, we can see that varying unknowns, such as the object density or the drag parameter, only cause the strewn field to slide along its center line. In the end, the main unknown is the width of that strewn field in the direction perpendicular to its center line, which can be several hundred meters up to 1 km.
Taking the example of the 1 January 2020 fall in Italy ), our determined strewn field with a 99% confidence level consisted of a thin strip 5.6 km long and 100 meters wide. The actual meteorite was found only 200 meters from the central line of this strip.
This demonstrates the accuracy of our method, the offset being mainly due to our approximating the meteorite shape as a sphere.

Statistics and network efficiency
One of the main objectives of the FRIPON network is to measure the unbiased incoming flux of extraterrestrial matter. In this section, we first present the raw statistics of detected falls. Next, we attempt to constrain the absolute flux of incoming material.    Table 2 gives the number of detections for each shower per year. This figure presents that the main showers are detected and that the sporadic meteors are uniformly distributed over the celestial sphere except for part of the southern hemisphere, which is not at present within the reach of FRIPON. Overall, sporadic meteors represent 55% of the data.  Of the objects, 55% are sporadic: their radiants cover the sky uniformly except for its southern part, which is invisible from European latitudes. The north toroidal sporadic source is visible in the top left corner and low speed objects are shown along the ecliptic plane coming from the anti-helion source. Fig. 8. Histogram of the absolute magnitude of all the events detected by the network showing that the exhaustive detection regime is only reached around mag −5. The slope is compatible with that obtained by previous studies such as Brown et al. (2002), as shown on Fig. 11, which describes the distribution of interplanetary matter from 1 cm to 1 km. The global shape of the histogram is similar to that in Ott et al. (2014), which is shifted as CILBO cameras are more sensitive than FRIPON cameras.

Raw meteoroid detections
To estimate that flux, we need an estimation of the cloud coverage, the percentage of operational stations, and the sensitivity of our network as a function of meteor brightness. Regarding that last Article number, page 29 of 54 point, Fig. 8 shows the absolute magnitude histogram after three years of observations. Assuming a power-law size distribution for interplanetary matter (Brown et al. 2002), it appears that FRIPON is clearly not fully efficient for events fainter than -5 in magnitude. This detection threshold is similar to that of the Prairie network (Halliday et al. 1996) and implies, as for other networks (Devillepoix et al. 2020), a minimum detection size of ∼1 cm for incoming meteoroids. We note that smaller objects can nevertheless be detected if their entry speed is high enough.
To calculate the efficiency of FRIPON, we only used the French stations as these were the first to be installed and France was fully covered in 2017. We considered its area, with a 120 km band added around it (Fig. C.1 in Appendix) for a total of 10 6 km 2 , which was the basis for the calculation. For ≥1 cm meteoroids (i.e., for magnitude < -5 fireballs), we obtained an average rate of 250 events/year/10 6 km 2 . Last, to estimate the incoming meteoroid flux for ≥1 cm bodies, we needed to correct for dead time (day time: 0.5 and average cloud cover: 0.4). The dead timecorrected meteoroid flux for ≥1 cm meteoroids is 1,250/year/10 6 km 2 , which is comparable to the 1,500/year/10 6 km 2 value given by Halliday et al. (1996). Our determination is raw and requires that we carry out a more detailed analysis in the future with more data. Our analysis shows that the network has reached a complete efficiency for the French territory for meteoroids larger than 1 cm.

Orbit precision
A precise determination of the orbit requires the extraction of a realistic initial velocity for the object. This can only be achieved by taking into account its deceleration in the upper atmosphere before the bright flight. Therefore our model of drag and ablation depends on three parameters (see section 3.3.3): the initial velocity V, a drag coefficient A, and an ablation coefficient B. Depending on the quality of the data-for example, the number of cameras, weather conditions, and distance of the camera to the bolide-these three parameters do not have the same influence on the trajectory calculation and cannot be determined with the same accuracy. We classified the meteors in three categories: 1. Those whose deceleration is hardly noticeable (A/σ A < 2), which represent 65% of all meteors.
2. Those for which only the deceleration is noticeable (A/σ A > 2 and B/σ B < 2), which represent 21% of all meteors. In those cases, the ablation is not observed.
3. Those for which both the deceleration and the ablation are noticeable (A/σ A > 2 and B/σ B > 2), which represent 14% of all meteors.
For dynamical studies, only the detections that fall in one of the last two categories (35% of all detections) can be used. The typical velocity accuracy is then 100 m/s, which is required both for the identification of meteorite source regions in the solar system (Granvik & Brown 2018) and for the search for interstellar meteoroids (Hajduková et al. 2019).
Article number, page 30 of 54 4.2. Dynamical properties of the observed meteoroids Fig. 9. Histogram of sporadic fireball entry velocities. Two populations can be observed: 1) low speed objects corresponding mostly to asteroidal orbits and 2) fast objects corresponding to TNOs or comet-like objects. This dichotomy has also been observed by Drolshagen et al. (2014) with the CILBO network for smaller objects.
In the following, we restrict our analysis to sporadic meteors. The histogram of initial velocities is shown in Fig. 9. It reveals two populations of meteoroids whose entry velocities differ by about 50 km/s, suggesting an asteroidal (55%) and a cometary (45%) population. This result can also be inferred from the histogram of meteoroid detections as a function of the inverse of the semimajor axis of their orbit (Fig. 10). This figure clearly shows a main belt population with semimajor axes between that of Mars and that of Jupiter, as well as a cometary population, possibly including Oort cloud material, with semimajor axes greater than that of Jupiter. Last, we note the presence of a few meteoroids with negative semimajor axes. However, rather than concluding that interstellar matter was detected, we attribute these events to large errors associated with the calculation of their initial velocity. As a matter of fact, these events have semimajor axes that differ significantly from that of the interstellar object 1I/Oumuamua.
It is clear that in more than three years of observation, FRIPON has not detected any interstellar object so far. This compares to results obtained by other networks such as CMOR (Weryk & Brown Article number, page 31 of 54 2004), who found that only 0.0008 % of the objects detected might be of interstellar origin; while a more recent work (Moorhead 2018) did not find interstellar candidate in CMOR data. In the case of the FRIPON network, only an upper limit of 0.1% can be given, but we expect the real value to be much lower. Hajduková et al. (2019) showed that no network so far has ever experienced a conclusive detection of an interstellar meteoroid. Most false detections are likely to stem from a bad error estimation, especially that of the initial speed, which requires an estimation of the drag coefficient.

Meteorite falls and first field search
Based on Halliday et al. (1989), about ten meteorites weighing more than 100 g must fall each year over the area covered by the FRIPON network. Table 3 lists the events that produced a computed significant initial and/or final mass. The fall rate that we observe for final masses equal or greater than 100 g is 2.7 per year. This value is compatible with that of Halliday et al. (1996), once corrected to take into account the 20% overall efficiency of the FRIPON network (see above), as this yields a corrected rate of 14 falls per year. Among these events, only 1 led to the recovery of meteorite fragments. This event occurred near Cavezzo in Italy  and was detected by PRISMA Article number, page 32 of 54 cameras. Further details regarding the meteorite and its recovery will be presented in a forthcoming paper. This recovery is particularly important in showing that it is possible to find a 3 g stone thanks to the mobilization of the public with the help of various media (e.g., internet and newspapers). This strategy has worked well and can be reproduced for all comparatively small falls (typically a few dozen grams). In such cases, it is clear that the chances of finding the stone are low and do not warrant the organization of large searches, while an appeal to the general public may be fruitful. In the Cavezzo case, the meteorite was found on a path by a walker and his dog.
We chose to compute the flux of objects with final masses greater than 100 g for which the accuracy is moderate to high (m/σ m > 2). This flux is 14 meteorites/year/10 6 km 2 (see above). We extrapolated it down to a mass of 10 g, assuming a power-law distribution of the final masses of the meteorites (Huss 1990), and obtained a value of 94 meteorites/year/10 6 km 2 , close to the value from Halliday et al. (1989); this is also based on fireball data. This value is, however, lower than the other estimates (Bland et al. 1996 andDrouard et al. 2019), which are based on field searches.
The Evatt et al. 2020 estimate based on the study of meteorites found in Antarctic blue ice gives a mid-range value that is consistent with all previous estimates.
Article number, page 33 of 54

Extension of the network
Significantly increasing the area covered by the network (by at least an order of magnitude) will be fundamental in increasing the recovery rate of meteorites, as this will lead to the detection, over a reasonable period, of a statistically significant number of very bright meteors that might be recovered on the ground as meteorites. Hence, there is a major interest in extending the FRIPON network over all of Europe and to other parts of the world. Such an extension has already begun (see Fig. 2) and will be pursued over the coming years. The development plan includes, as a priority, the densification of the European coverage as well as its extension to southern countries such as Morocco, Algeria, and Tunisia. For Spain, FRIPON is complemented by the SPMN network (Trigo-Rodríguez et al. 2004), with which we already collaborate for trans-national events and with whom we organized a search for a possible meteorite fall in January 2019. Such a southern extension would be sufficient to generate a network area about ten times larger than that of metropolitan France. In addition, the network is currently also being developed in Canada in North America and in Chile in South America. Fig. 11 shows that 30 objects larger than one meter fall on Earth (510 × 10 6 km 2 ) every year. Taking into account the current surface area of the FRIPON network, the average expected detection rate of such objects is limited to an average of one in ten years. Extending the area of the network is thus a necessity to reach an acceptable detection rate for 1 meter objects. An extension to Europe and North Africa would make it reach a surface area of 6 × 10 6 km 2 , which is comparable to that of the Australian DFN network (Devillepoix et al. 2016), leading to a probability of a one-meter event approximately every years.

Software
The reduction pipeline is operational and only requires minor improvements. The acquisition software FreeTure still shows a surprisingly high false detection rate, which requires that daylight observations are turned off at the moment. A new version using deep learning techniques is being developed so that daytime observations will become possible. The development of a tool to compute the light curve of heavily saturated events (Anghel et al. 2019b) is also planned.

Hardware
The hardware currently in use in the network corresponds to pre-2014 technology. A complete hardware update after five years of utilization is thus desirable to improve the temporal resolution of the light curves and the performance and flexibility of the acquisition computers. A non-exhaustive list of improvements includes upgrading from CCD to CMOS detectors and switching the current PCs to Raspberry Pi4 single board computers (SBCs).
Article number, page 34 of 54 Fig. 11. Flux of small near-Earth objects colliding with the Earth (Brown et al. 2002). Data are shown over a range of 14 magnitudes in energy. The statistical model is based on near-Earth population for big sizes and, for the smaller objects, it is derived from a decade-long survey of ground-based observations of meteor and fireballs. The FRIPON network lies exactly between minor planets (detected by telescopes and planetary impacts) and interplanetary dust (detected by meteor networks). The solid arrow corresponds to FRIPON nominal mode; the dashed line is for rare events, observable by FRIPON but with a very low probability.
In addition, a prototype of an all-sky radiometer is presently under development (Rault & Colas 2019), to resolve the saturation issue and improve on the bandwidth of the cameras. This radiometer covers the visible and near-infrared wavelengths. It is based on a 16 PIN photodiode matrix, followed by a trans-impedance amplification chain and a 14 bit industrial USB data acquisition module, which samples at a rate of 20 kHz. As an example, we superimposed on Fig. 12 the FRIPON camera light curve for an event of magnitude -9.5, which occurred on 14 August 2019 at 03:07:02 UTC and the corresponding high data rate radiometer light curve.

Radio
The aim of FRIPON radio receivers is an accurate measurement of meteor velocities through the Doppler effect, allowing a much better determination of the orbital data (especially semimajor axes). In Table 4, we present the value of the initial velocity and effective surface-to-mass ra- achieved with the radio data leads to errors one order of magnitude lower compared to that achieved with only the visible images. However, it seems at present that only about 30 % of the optical detections lead to a detectable radio signal and that several bright radio events do not have any visible counterpart. For this reason, radio data have not been widely used yet, and further work is needed to improve our understanding of the complex phenomena associated with the generation of radio echoes by the plasma surrounding the meteors. Over time, we came to the conclusion that detailed information on the fragmentation and final destruction of bolides might also be obtained thanks to the head echoes produced by the GRAVES HPLA radar. Last, we sometimes detected unexpected oscillations on the usually smooth Doppler shift curves (Rault et al. 2018), which indicates cyclic fluctuations on the radial positions of the radar cross section (RCS) of the plasma envelope surrounding the meteor bodies (see Fig. 13).  (Cecconi et al. 2018). Since the Summer of 2017, the software FreeTure contains an experimental real-time algorithm for the detection of TLEs. This algorithm runs along with the meteor detection part on selected stations with a view to help localize TLEs observed by the future CNES space mission TARANIS (Blanc et al. 2017). The FRIPON network infrastructure can also be used to conduct large-scale light pollution monitoring campaigns using the all-sky calibration images collected over time (Jechow et al. 2018).

Observation from space
The network can also be extended vertically by combining space measurements with ground measurements. Space-borne observations have several advantages, such as providing a wide geographical coverage with one camera, longer recording times, and no weather constraints. The small satellite sector is evolving very quickly, opening up new opportunities for scientific missions (Millan et al. 2019). In particular, relatively inexpensive missions make it possible to design swarms of satellites or even constellations dedicated to monitoring the Earth and therefore meteors. In this framework, a Universitary Cubesat demonstrator called Meteorix is under study (Rambaux et al. 2019). The Meteorix mission is dedicated to the observation and characterization of meteors and space debris entering the Earth's atmosphere. The orbit chosen for Meteorix is a low Earth sunsynchronous orbit at an altitude of 500 km. Such configuration will make it possible to detect on average a sporadic meteoroid entry per day and about 20 meteors during a major meteor shower.
The nominal mission lifetime is one year. Three-dimensional astrometry and photometry would become possible in case of a detection over the FRIPON network.

Conclusion
The FRIPON scientific network, originally developed to cover the French territory, is now a fully automated network monitoring fireballs above part of western Europe and a small fraction of Canada. As of today, it consists of 150 cameras and 25 radio receivers covering an area of about 1.5 × 10 6 km 2 . The level of automation of the network is such that a recovery campaign can be triggered only a few hours after a meteorite reached the surface of the Earth.
The FRIPON scientific project has been monitoring meteoroid entries in western Europe since 2016, thereby allowing the characterization of the dynamical and physical properties of nearly 4,000 meteoroids. It has thus allowed us to significantly enhance the statistics of orbital parameters of meteoroids, while also searching for possible interstellar meteoroids. The FRIPON observations show that the distribution of the orbits of incoming bolides appears bimodal, comprising a cometary population and a main belt population. Sporadic meteors amount to about 55% of all meteoroids. In addition, we found no evidence for the presence of interstellar meteoroids in our sample. Overall, Article number, page 38 of 54 it appears that the range of sensitivity of the FRIPON network encompasses particles originating both from comets and asteroids. A first estimate of the absolute flux of meteoroids bigger than 1 cm amounts to 1,250 /year/10 6 km 2 , which is a value compatible with previous reports. We also estimate the flux of meteorites heavier than 100g to 14/year/10 6 km 2 , which is a value compatible with data from other fireball networks but lower than those obtained from collecting meteorites.
Finally, the first meteorite has been recovered in Italy following observations by the PRISMA network, a component of the FRIPON network.
Further extension of the FRIPON network is under way. In the coming years, it will be extended to North and West Africa as well as Canada and to the southern hemisphere in South America and South Africa. The goal is to reach a size large enough to allow the recovery of at least one fresh meteorite per year. In addition to the geographical extension of the network, technical developments will be conducted to improve the photometry of saturated images. Moreover, we plan to implement new algorithms in the detection software, so that daytime observations become possible and useful.
Finally, we plan to fully exploit the radio network, both to improve current orbits and to reach a better understanding of the physical mechanism of meteoroid entries.