the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Spectral observations at the Canary Island LongBaseline Observatory (CILBO): calibration and datasets
Detlef Koschny
Regina Rudawska
Salvatore Vicinanza
Stefan Loehle
Martin Eberhart
Arne Meindl
Hans Smit
Lionel Marraffa
Rico Landman
Daphne Stam
The Canary Island LongBaseline Observatory (CILBO) is a doublestation meteor camera setup located on the Canary Islands operated by ESA's Meteor Research Group since 2010. Observations of meteors are obtained in the visual wavelength band by intensified video cameras from both stations, supplemented by an intensified video camera mounted with a spectral grating at one of the locations. The cameras observe during cloudless and precipitationfree nights, and data are transferred to a main computer located at ESA/ESTEC once a day. The image frames that contain spectral information are calibrated, corrected, and finally processed into line intensity profiles. An ablation simulation, based on Bayesian statistics using a Markov chain Monte Carlo method, allows determining a parameter space, including the ablation temperatures, chemical elements, and their corresponding line intensities, to fit against the line intensity profiles of the observed meteor spectra. The algorithm is presented in this paper and one example is discussed. Several hundred spectra have been processed and made available through the Guest Archive Facility of the Planetary Science Archive of ESA. The data format and metadata are explained.
Insights into the dust composition of our solar system have been obtained during the past decade. Groundbased observations from ultraviolet to radio wavelengths (Snodgrass et al., 2017), spacebased telescopes, such as HST (Hadamcik et al., 2016), direct spacebased remote sensing and in situ observations as obtained by the Rosetta spacecraft during its 2year nominal mission lifetime (Rubin et al., 2020; Rivilla et al., 2020; Altwegg et al., 2020), and laboratory analysis of samples returned back to Earth, e.g., by the STARDUST mission (Brownlee, 2014), have provided valuable datasets. A summary of our current knowledge of cometary dust is given in LevasseurRegourd et al. (2018). For interplanetary dust in general, see Koschny et al. (2019) and references therein. The knowledge of dust allows us to derive further information on, e.g., the origin and further evolution of comets (Rubin et al., 2020). Another indirect way to observe dust is the observation of meteors and thus the observation of the ablation process of meteoroids through Earth's atmosphere. To allow the chemical analysis of the ablation process, spectral observations are needed and have been successfully demonstrated by several teams in the past decades (Vojáček et al., 2019; Borovička et al., 2005; Abe et al., 2020; Ward, 2017; Madiedo, 2015; Jenniskens et al., 2014; Rudawska et al., 2014; Zender et al., 2002; Matlovič et al., 2019; Zender, 1994).
ESA's Meteor Research Group (MRG) maintains a number of cameras acquiring the meteor ablation phenomena in Earth's atmosphere, also known as shooting stars or fireballs for brighter events. A general overview of our activities is given in Koschny et al. (2015a).
Besides supporting campaignbased observations (Koschny et al., 2012; Vaubaillon et al., 2015), the MRG maintains a permanent double station on the Canary Islands acquiring videos from two distinct locations. The parallel acquisition of meteors allows deriving geometric information through trigonometric algorithms (Koschny et al., 2013a). In one of the CILBO stations, a dedicated camera is equipped with a spectral grating to acquire spectral information of the brightest meteors throughout the ablation process in the Earth's atmosphere (Zender et al., 2014).
In this paper, we describe the data reduction process from the acquired images to the final spectra.
A full algorithm to derive meteor spectra was discussed in the middle of the twentieth century by Ceplecha (1971). His observations were executed at the Ondr̆ojev Observatory and based on the usage of photographic cameras acquiring a long exposure in the visual wavelength regime. Spectral plates or spectral prisms diffracted the lights source, e.g., light from a star or from a meteor, into its spectral components. Once the photographic plate was developed, the light curve was measured using a device called a microdensitometer, e.g., the Mann 1140 or the Lirepho2 microdensitometer (Lirepho is an acronym for “lichtelektrisches Registrierphotometer” produced by Zeiss). The wavelength identification along the obtained light curve was directly computed due to the availability of the zero order of the meteor using an earlier algorithm of Ceplecha (1961). Given an infalling vector (the zero order of the meteor) on the grating plate, the algorithm computes the (x, y) position for each wavelength on the focal plane, taking into account rotations between the grating plate reference frame and the optical axes.
Quantitative analysis of meteor spectra was introduced by several observers (Borovicka, 1994a, b; Jenniskens, 2007). The observations and analysis of Borovicka (1994a) identified two main thermal components during the ablation process, a lowtemperature component of about 4000 K and a hot component of 10 000 K. Jenniskens (2007) provides a summary of all identified elements and confirms the two thermal components, complemented by the formalism to derive the number densities for chemical elements. The calibration of the wavelength scale is achieved by laboratory measurements using a mercury–argon calibration source applied over different source incidence angles. These allow at a later step in the algorithm confirming the correct distances between potential line intensities of different elements. The wavelength scale is described as nearly linear, which overcomes the problem of the diffraction hyperbola that is the effect of a curved spectrum due to the usage of a nongnomonic lens. The projection of a meteor spectrum onto an image is only linear under the assumption of an optimal lens (Ceplecha, 1961). The projection of an optimal lens onto an imaging plane, i.e., chargecoupled device detector, is gnomonic. In video systems used in the last 3 decades this is, however, not given, and Dubs and Schlatter (2015) propose an additional image transformation to transform an acquired image into an gnomonic projection.
This paper presents the CILBO observation locations and hardware (Sect. 2) and the data processing steps applied to the meteor spectra: the individual video frames are radiometrically processed (Sect. 3) and converted into the FITS format. In the case that the visual camera has observed the zeroorder and the grating camera has observed the firstorder spectrum of the meteor, the wavelengths are computed for each individual frame (Sect. 4) and we obtain a spectral profile. In the following, an intensity correction is applied to the spectral profile (Sect. 6). To determine the temperature regime and the chemical elements that optimally match the obtained spectral profile, we simulate the spectral response of several elements at different temperatures and use these simulated values by applying the Bayesian methodology (Sect. 7) combined with a Markov chain Monte Carlo (MCMC) sampling. As a result we obtain the line intensities, which are then further corrected to the corresponding number densities of the elements by integrating the column density over the line of sight (Sect. 6).
The data pipeline is then briefly discussed (Sect. 8) using one meteor spectrum. In Sect. 9 we describe the datasets that we make available in the Guest Archive Facility of the ESA's Planetary Science Archive Supplement (Zender and Koschny, 2023).
2.1 Equipment
The Meteor Research Group (MRG) of the European Space Agency operates the doublestation meteor camera system CILBO (Canary Island LongBaseline Observatory) (Koschny et al., 2014). Several imageintensified video cameras are operated on Tenerife at the Izaña Observatory and on La Palma at the Observatorio del Roque de los Muchachos. A photograph of one of the cameras is shown in Fig. 1. Table 1 provides details on the individual cameras and locations.
The cameras ICC7 and ICC9 acquire videos with their lines of sight pointing to a point 100 km above the Earth's surface, at a location between the two islands, at 28.5533^{∘} N and −17.1667^{∘} W. The socalled spectral camera ICC8 is offset to ICC7 by 23.0^{∘} to capture the firstorder spectrum of meteors passing through the field of view of camera ICC7. A technical overview drawing is shown in Fig. 2.
The camera systems use a lowdistortion machinevision lens from Fujinon, imaging the sky onto the entrance aperture of a secondgeneration image intensifier (DEP XX1700). The intensifier has a spectral sensitivity from about 350 to 850 nm. It is fibercoupled to a chargecoupled device (CCD), which is read out via a Toshiba Teli CS8310Bi PAL video camera. The PAL format has 576×768 pixels squared. Full frames are read out at a rate of 25 frames per second with 8bit dynamical accuracy. With the given field of view the pixel scale is 2.3^{′}. A dew remover (type Kendrick) is wrapped around the objective lens to avoid dewing.
2.2 Control software
The core of the software is the detection software MetRec (Molau, 1999) that analyses the image streams from ICC7 and ICC9, independently from each other. The ICC8 camera image stream is read into a separate detection software called SpecRec on the computer of ICC7. In the case that the MetRec detection software of camera ICC7 detects a meteor, a hardwired trigger signal (RS232 serial port) is transferred to the ICC8 detection software. The image stream, including the three images acquired before the trigger signal arrived, is stored onto a hard disk. The MetRec software stores for each individual frame of each detected meteor event the right ascension and declination of the photometric center of the meteor and the magnitude. It also stores an estimate of the cloud coverage with 1 min accuracy. The computers are synchronized to a networkprovided time signal every 2 min. A more detailed description of the CILBO setup and its cameras is available in Koschny et al. (2013b).
2.3 System performance
Since the start of operations in 2012, the cameras ICC7 and ICC8 have individually recorded over 90 000 meteor events, as indicated in Table 2. To apply the wavelength calibration to these events, the calibration routine we have developed requires that both zeroorder images recorded by ICC7 and firstorder spectra from ICC8 be available simultaneously (further explained in Sect. 4). About 14 000 were recorded by both systems simultaneously. The distribution of these events per year is listed in the column “Simultaneous events” in Table 2.
See Koschny et al. (2017) for a more detailed analysis covering 2012 to 2017. The 3D trajectory and heliocentric orbits of these meteoroids were computed using the MOTS (Meteor Orbit and Trajectory Software) code (Koschny and Diaz del Rio, 2002; Koschny et al., 2015b). The column “ICC8 spectral events” in Table 2 refers to the identified full spectra obtained. In this paper, we focus on the analysis from 2012 to 2017, as afterwards the transmissivity of the grating began to degrade. We have also excluded spectra that are not fully contained in all frames of the meteor event.
Each detected meteor event results in a number of bitmap (BMP) image files for ICC7 and ICC8, as well as an “information file” for ICC7 which contains measurement results for the meteor. The image files contain the three images acquired before the first detection of the event, as well as three images after the meteor event is lost.
As the CILBO setup did not foresee a camera door or shutter, a dark current image is only available from the time of the initial setup of the overall system in 2012. A flatfield image is produced synthetically by taking the median values for each pixel of an image sequence obtained on a clear night and then used throughout the years. First, all images, including the pre and postimages, are radiometrically corrected by dark current subtraction and the flatfield division.
To compensate for acquisitions taken during cloudy conditions or during dawn or dusk, the six radiometrically processed preevent and postevent images are stacked together into a background image: for each pixel, the median value of the six available pixels is taken.
Each radiometrically processed image in the meteor sequence is then subtracted by this background image. This operation also subtracts the stars from the images and would result in a low background image in the case that no spectrum is visible. A socalled total image is created by selecting for each pixel location the median pixel value of all images of the meteor sequence obtained so far. The total image provides a summary of the full meteor event and allows having a quick look at an event to judge the spectrum brightness.
The obtained image sequences, as well as the total image, dark current, and flat field, are stored in a single FITSformatted file referenced as the level0 dataset (see Sect. 9 for a detailed description).
Different to other algorithms described so far (Borovicka, 1994a; Dubs and Schlatter, 2015), our algorithm determines the pixel position as a function of wavelength in the ICC8 image, computed from the position in the zeroorder image from ICC7. This is done by taking the angular offset between the two cameras into account. During the operational period of the cameras, this angular offset was derived several times from measurements and confirmed to be stable. The basic idea of this step can be best compared to the twodimensional diffraction formula as we know it from the singleslit optical diffraction experiment.
Knowing the wavelength λ and the slit width d, one can compute the angle θ and thus the location of the maxima and minima.
Assuming Fraunhofer's diffraction at the grating, the light of a star or meteor is considered parallel entering light rays, and its resulting diffraction can be expressed using the threedimensional grating formula as provided in Eq. (2).
where

$\overrightarrow{{v}^{\mathrm{lens}}\left(\mathit{\lambda}\right)}$ represents the direction of the diffracted light with respect to the lens for each wavelength λ;

(x,y)^{grating} represents the coordinates on and in the plane of the grating, at which the light enters;

n is the grating order, which has +1 as a nominal value;

λ is the wavelength (in meters);

d is the groove density (in number of grooves per meter).
With the information from ICC7, especially the right ascension and declination of the meteor, the zero order, for each individual frame, this diffraction vector $\overrightarrow{{v}^{\mathrm{Lens}}\left(\mathit{\lambda}\right)}$ for any wavelength λ can be calculated. To derive the (x,y) coordinates at the CCD frame, a sequence of reference frame transformations has to be applied. We further assume that all reference frames are righthanded ones and that the reference frames' z axes are aligned with each other, in particular the CCD, the optical, and the grating z axes.
As a first step, we have to derive the (x,y)^{grating} coordinates of Eq. (2) by applying the following reference frame transformations.

ℱ_{celestial}⟶ℱ_{local horizon}
This transformation converts a meteor position (right ascension and declination in J2000, as obtained from ICC7 information) for a given time from the Earthcentered celestial frame into an observercentered frame (altitude, hour angle), with the horizon as the fundamental frame, using the Python routine eq2hor. 
ℱ_{local horizon}⟶ℱ_{grating}
This transformation converts a meteor position (altitude and hour angle) for a given time from the observerbased frame into a gratingcentered reference frame. This conversion requires knowledge of the azimuth and elevation of the +x axis and +y axis of the grating, respectively. Assuming that it is identical to the optical axis of the complete camera system, it can be computed from the right ascension and declination of the center of the field of view at any given time. The boresight angle, the angle that the z axis of the grating is tilted, allows for a nonhorizontal plate. Unfortunately, this value was not measured at the time of the CILBO installation and we use a rotation value of −11.0^{∘} as a bestfit value in the z axis. The tilt angle describes a rotation of the grating, a rotation of the x–y plane of ℱ_{grating} perpendicular to the +z axis of the grating. If these angles are known, a transformation matrix can be constructed to convert any vector in the horizontal frame into a vector in the grating frame as follows.$$\begin{array}{}\text{(3)}& \begin{array}{rl}{\mathit{x}}^{\mathrm{grating}}& ={M}_{\mathrm{azimuth}}\cdot {M}_{\mathrm{elevation}}\cdot {M}_{\mathrm{boresight}}\cdot {M}_{\mathrm{tilt}}\\ & \cdot {\mathit{x}}^{\mathrm{local}\phantom{\rule{0.25em}{0ex}}\mathrm{horizon}}\end{array}\end{array}$$ 
${\mathcal{F}}_{\mathrm{grating}}\u27f6{\mathcal{F}}_{\mathrm{lens}}^{\mathit{\lambda}}$ can be determined with Eq. (2) and produces a vector for each wavelength λ.
From the vector $\overrightarrow{{v}^{\mathrm{lens}}\left(\mathit{\lambda}\right)}$ one could directly transform into the CCD frame and obtain the (x, y) coordinates of the image by normalizing the z component of $\overrightarrow{{v}^{\mathrm{lens}}\left(\mathit{\lambda}\right)}$ to 1. From the known field of view of the camera, one can construct vectors from the CCD center to the fieldofview corners assuming again that the z component equals 1. From these relations we can derive the x position and y position in the CCD frame of the meteor directly.
We follow a different approach and transform the obtained vectors back in an inverse transformation into ${\mathcal{F}}_{\mathrm{celestial}}^{\mathit{\lambda}}$ by applying the inverse matrices of Eqs. (2) and (3). As a result, we obtain the celestial coordinates for each wavelength. With knowledge of the plate constant of the camera, these coordinates can be directly transformed into the image (x, y) coordinates. We used the Python method all_world2pix() of the astropy.wcs class. To allow this last transformation, we use the web service offered by astrometry.net (Lang et al., 2010) (http://nova.astrometry.net, last access: 31 March 2023) to obtain a registered image. This registration was done once with an image under good viewing conditions. As the webservice registration also takes the geometric distortions into account, the algorithm does not need to apply any additional distortion correction routines.
The calibration step results for each frame in a list of wavelengths and its corresponding digital value at the computed image position and is stored in a single FITS file labeled as the level1 dataset (see Sect. 9 for a detailed description).
It is noted that the accuracy of our algorithm – determined through the correct identification of the (x, y) position of a maxima for a chemical element – was not always satisfactory. As the bright pixels representing the zeroorder meteor in ICC7 are spread over many pixels, a centroid algorithm is used in METREC that did not always correspond to the center pixel of the meteor as visually observed. The difference was within a few pixels, but the impact on the overall accuracy of the algorithm was too large to neglect. The problem was corrected in a later processing step (see Sect. 7.2) by adding the wavelength itself as one of the Bayesian parameters.
5.1 Spectral response of ICC8
The spectral response curve for ICC8 was derived from images containing the firstorder spectrum of Vega.
Vega is an A0type star with a wellknown emission spectrum in the visible range. A reference spectrum of Vega was taken from the STIS CALSPEC Archive at https://archive.stsci.edu/hlsps/referenceatlases/cdbs/current_calspec (last access: 31 March 2021). To obtain an ICC8 Vega spectrum, four frames acquired on 4 consecutive days were selected. Each frame is then radiometrically processed and the wavelengths are derived as provided in Sects. 3 and 4. The final observed Vega spectrum is the median of the four calibrated and atmosphericcorrected spectra. The ICC8 spectral response curve is the ratio between the final observed spectrum and the reference spectrum, fitted with a 12thorder polynomial and normalized to 1, as shown in Fig. 3. To validate the ICC8 spectral response, acquired spectra of Algol, Castor, Deneb, Elnath, and Vega underwent our calibration pipeline and were compared to the spectra provided in CALSPEC or HyperLeda (http://leda.univlyon1.fr/, last access: 31 March 2021).
5.2 Correction of atmospheric attenuation
As the light emitted from the ablation process undergoes atmospheric extinction, we apply a correction to the spectral response curves as obtained in Sect. 5. The scattering and absorption correction is applied to each individual frame and depends on the exact observation geometry, especially the zenith angle of the observed meteor θ_{met} and the optical thickness b_{met} of the atmosphere. The extinction law can be described as a correction vector (Appenzeller, 2013), as follows:
where $\overrightarrow{{F}_{\mathrm{met}}}$ represents the intensity of each feature of the observed meteor spectrum and $\overrightarrow{{F}_{\mathrm{0},\mathrm{met}}}$ the intensities of the true, atmosphericcorrected, spectrum. The atmospheric optical thickness is a coefficient specifying the attenuation of light through the atmosphere; $\overrightarrow{{b}_{\mathrm{met}}}$ is a vector varying with wavelength. In the case that we consider only scattering contributing to the atmospheric extinction and the aerosol absorption, $\overrightarrow{{b}_{\mathrm{met}}}$ =$\overrightarrow{{b}_{\mathrm{sca}}}$ is obtained as in Eq. (5), with N_{gas} the column number density of gas molecules (obtained using Eq. 6), σ_{sca} the scattering cross section of gas molecules (obtained with Eq. 7), and AOT the contribution of the aerosol optical thickness. The AOT is taken from NASA's AERONET (NASA's AERONET database is found at https://aeronet.gsfc.nasa.gov/new_web/aerosols.html, last access: 1 March 2023). The file which contained the data used was named “aerosol optical depth” (AOD) with precipitable water and Ångström parameter. database. The data were recorded by an AERONET sun photometer in Izana, Tenerife, located at around 1.17 km ground distance from our ICC8 camera. The data are provided at discrete, irregular intervals from 340 to 870 nm.
Here, m_{gas} is the air mass (obtained for dryair conditions, 28.96 g mol^{−1}) and g the gravitational acceleration (9.806 m s^{−2}), both considered constant through the atmosphere. p_{bottom} and p_{top} are the pressures at the altitudes of the atmosphere corresponding to the beginning and end of the meteor ablation. In Eq. (7), n is the gas refractive index (assumed constant), N_{L} is Loschmidt's number (2.54743 10^{25} m^{−3}), and δ is the gas depolarization constant (0.03 for terrestrial air, ignoring the λ dependence).
As several of the parameters depend on the line of sight between the station ICC8 and the meteor location (and thus the meteor altitude), the atmospheric correction is applied to each individual frame of a meteor event. The altitudevarying parameters are computed as follows: g and n are approximated as constants between ICC8 (the troposphere) and the meteor location (the thermosphere).
p_{top} is calculated using Eq. (8):
where p_{0} is the SI standard pressure (101325 Pa), L the temperature lapse rate (9.76 K km^{−1}, dry air), T_{0} the SI standard temperature (288.16 K), M the molar mass (0.029 kg mol^{−1}, dry air), and R_{0} the universal gas constant (8.314 J mol K^{−1}).
A detailed discussion of the assumptions, especially related to the atmospheric values, the improvements obtained, and potential further improvements of the correction using the atmospheric extinction correction, can be found in Vicinanza et al. (2021). Figure 4 shows a spectral profile before and after the atmospheric extinction correction was applied.
This calibration step results for each frame of a meteor event in a list of wavelengths and the corresponding spectralcalibrated and atmosphericcorrected value are stored in a single FITS file labeled as the level3 dataset (see Sect. 9 for a detailed description).
The aim of the intensity calibration is to convert the spectral intensity, measured in CCD pixel counts often expressed as arbitrary digital units (ADU), to SI physical units (W m^{−2} nm^{−1}). This allows inferring the true and absolute composition of elements in a meteoroid during the extraction of the line intensities of individual chemical elements (as will be described in Sect. 7). The intensity calibration is done in two steps following Jenniskens (2007): the calculation of the zeropoint bias and the estimation of the coefficient of intensity calibration.
The calculation of the zeropoint bias involves a photometric calibration. This is performed using the observations of A0V standard stars, simultaneously recorded in the background of meteor observations. We used Vega, which has a Vband magnitude m of approximately 0. Knowing the spectral intensity of Vega observed from ICC8, expressed as the sum of ICC8 pixel brightness ∑F, the zeropoint bias C is obtained using the relation
where ∑F is Vega's sum of the pixel brightness, and C is the zeropoint bias we want to obtain. By applying this step to our set of CILBO observations, a value of 10.81±0.1 was obtained for the zeropoint bias. Theoretically, this step should be performed separately for each individual meteor event using the stars observed in the meteor's background. However, since it was observed that the zeropoint bias remains approximately constant throughout different events, a constant term was applied during our research to the calibration of all events.
The zeropoint bias is then used for the estimation of the coefficient of intensity calibration, k, according to the relation below (Jenniskens, 2007):
where C is the zeropoint bias and EW is the equivalent width of the spectral response (measured as the fullwidth halfmaximum) of the camera (see Fig. 3) and allows performing a correction for the instrument aperture. Equation (10) is obtained knowing that the zeromagnitude A0V standard star Vega has a flux of $\mathrm{3.67}\times {\mathrm{10}}^{\mathrm{11}}$ W m ^{−2} nm^{−1} at 548.3 nm (V band). On average, the coefficient k, separately calculated for each event, was around $\mathrm{2}\times {\mathrm{10}}^{\mathrm{17}}$ W m^{−2} nm^{−1} ADU^{−1}. Once this coefficient is obtained, the calibrated meteor spectral intensity (at each wavelength λ) is obtained from the observed meteor flux F following the relation
F is the average of the flux over 10 binned pixel rows. Since the coefficient k is obtained using Vega star observations and we want to perform the intensity calibration on meteor observations, the correction factor f is needed. This accounts for the dependence of the intensity calibration coefficient k_{corr} on the exposure time of the star measurement compared to the exposure time of the meteor measurement.
The exposure time of a Vega observation is equal to the observation time for the entire meteor event. Since we perform the intensity calibration for each meteor frame (10 binned pixel rows), the exposure time of the meteor is a ratio of the time delay between consecutive frames. Our cameras have a frame rate of 25 FPS and thus a 0.04 s interval between consecutive frames (Koschny et al., 2013b). The exposure time of the meteor in each frame can thus be obtained as
where $\mathrm{\Delta}{y}_{i,i+\mathrm{1}}$ is the difference in y coordinate between consecutive MetRec recordings and depends on the meteor speed.
In conclusion, the coefficient of intensity calibration k – used to convert the spectral intensity of meteor spectra from some arbitrary digital units (ADU) to SI physical units (W m^{−2} nm^{−1}) – was considered constant for diverse observations of the same meteor event, while it changed between different meteor events. The change mostly depends on the different exposure times of different meteor events (related to the different meteor speed), the different conditions of atmospheric extinction encountered at different meteor entries, and the optical and camera characteristics (as obtained through the zeropoint bias). In our cases, the coefficient k was of the order of $\mathrm{2}\times {\mathrm{10}}^{\mathrm{17}}$ W m^{−2} nm^{−1} ADU^{−1}.
To determine the chemical elements and their contribution to a meteor spectral profile, a linebyline emission simulation has to be executed for each potential chemical atom, molecule, or oxide. The atomic behavior of each species needs to be modeled, taking into account a number of model input parameters, especially the ablation temperature and the chemical element number densities. As the selection of slightly different input parameters results in different spectral profiles that might or might not match a meteor spectral profile from the level3 datasets (see Sect. 9.3), we have to solve an optimization problem. The CILBO calibration pipeline uses the Plasma Radiation Database (PARADE) linebyline emission simulation tool (Smith, 2006). We use the Bayesian formalism to solve the optimization problem, and the PARADE tool is called from within the Bayesian method; both are discussed in the following sections.
7.1 The Plasma Radiation Database
The synthetic spectra are produced by ESA's Plasma Radiation Database (PARADE) tool, originally used to simulate variation of probe entries into planetary atmospheres (Smith, 2006; Pfeiffer et al., 2003). PARADE calculates the energy state transitions in atoms and molecules and provides the emission coefficient ϵ in radiance W m^{−3} sr^{−1} m^{−1}.
While the number of gaseous elements was originally restricted to the main atmospheric constituents of the atmospheres of Titan, Earth, and Mars, several new species have been added to support our spectral analysis activities for the CILBO acquired spectra. A detailed discussion is given in Loehle et al. (2021), extending the available species to include Na, K, Ti, V, Cr, Mn, Fe, Ca, Ni, Co, and Li. The positional information of an emission (or absorption) line is defined by the energy gap between the electron excitation levels as defined by, e.g., Herzberg (1950).
The energy levels of atoms are exclusively defined by their electron configuration and energy levels. These can be determined from the NIST database (Ralchenko and Kramida, 2020) directly. For molecules, the vibration and rotation of the nuclei are additional parameters to consider. A simple radiation model for diatomic molecules based on a Boltzmann distribution of the vibrational bands and under the assumption of thermal equilibrium is presented in Eq. (15) (Loehle et al., 2021).
Here, h denotes the Planck constant, k the Boltzmann constant, c the vacuum speed of light, A_{v,u} the Einstein coefficient for a spontaneous emission from a lower to a higher state N_{J}, J the rotational quantum number, and T the temperature. In addition, we have further extended the PARADE tool to the following molecules: CO, CO_{2}, AlO, CaO, FeO, MgO, NO, and TiO.
To avoid the repetitive execution of PARADE simulations in the upcoming processing step (see Sect. 7.2), a database is created containing for each species the emission profile over the spectral regime for temperatures in a temperature range from 2000 to 15 000 K in a 50 K granularity. When building up the synthetic emission lines, FeI is processed over the full spectrum, as well as over a limited spectral width that explicitly covers the FeI (15) multiplet response, ranging between 520 and 550 nm. The modeled spectra include the contribution of the Planck background at the individual temperatures; for details see Loehle et al. (2021). In summary, the following atoms, molecules, and ions are available in our database: NaI, MgI, FeI, FeI(15), CaI, SiI, AlI, MnI, CrI, Ar, C_{2}, CH, CN, CO_{2}, FeO, CaO, AlO, MgO, O_{2}, OH, TiO, Na^{+}, Mg^{+}, Fe^{+}, Ca^{+}, and Ar^{+}.
7.2 Optimization of the spectral fitting using the Bayesian methodology
The PARADE tool allows computing the radiative emission of gas species by calculating the line position as well as its intensity and distribution, i.e., the line profile (Herzberg, 1950; Loehle et al., 2021; Rudawska et al., 2020), for a given temperature and given number density of a selected gas species.
Taking the reduced line profiles as obtained after the calibration steps (the dataset level 3), a fit to the set of simulated results from the PARADE tool must be found. To solve this inverse problem we use the Bayesian treatment of the inference, Bayes' theorem, and apply it to the meteor ablation process.
Having priors is of key importance in the Bayesian framework as priors express our existing knowledge about the parameter space. Parameters for which a lot of knowledge is available are named informative priors in the Bayesian framework. They are updated by the algorithm based on new data provided. In our case, we have three informative priors, all being input parameters to the PARADE simulation. The first is the total number of atoms or molecules of a gaseous species directly relating to the maximum values of the simulated line profiles. We consider this number to be uniformly distributed. The temperature of a meteor ablation process is an informative prior and we assume a uniform distribution. Following the identification of two distinct thermal components in the meteor ablation by Jenniskens (2007), we take two distinct temperature components, T_{1} and T_{2}, as informative priors.
Uninformative or ignorance priors are very diffuse and do not constrain the parameters. These are used when we do not have a lot of knowledge about the parameter beforehand and are not used within the PARADE simulations. We use three uninformative priors in our framework.

The PARADE tool takes into account several linebroadening effects within its simulation like the Van der Waals, or natural or Stark broadening. As the errors associated with these models are not yet fully quantified, we have introduced the standard deviation of the instrumental broadening, σ, as an uninformative prior and consider its distribution function uniform.

We have reported in Sect. 5 an uncertainty in the position measurement of the meteor center that we are not able to characterize further and have therefore introduced as an uninformative prior the uncertainty along the wavelength axis, Δλ, of the measured line profiles. For this uninformative prior we also assume a Gaussian distribution.

Moreover, the line profile itself constitutes a level of signal noise that propagates through the algorithms (as described earlier) and is difficult to characterize. We thus introduce another uninformative prior, the line profile background b, assuming a Gaussian distribution.
Finally, we need to convert the modeled synthetic spectra in irradiance (W m^{−3} sr^{−1} nm^{−1}) into radiance (W m^{−2} nm^{−1}). This is done by performing a solid angle conversion assuming the meteor to have the shape of an isotropic sphere (Jenniskens, 2007). We introduce a final uninformative prior, the size of the ablating body, the meteor, in the direction of flight R, assuming a uniform distribution. From R we derive the luminous area of the radiating body as π*R^{2}.
Table 3 lists all priors and their initial settings.
The Bayesian treatment of the inference, Bayes' theorem, can be written as
P(θY) denotes the probability of a model applied over a parameter space θ, in our case the temperature, the instrumental line broadening, the uncertainty in the wavelength, and the background count rate, given a set of data samples $Y=({y}_{\mathrm{0}},{y}_{\mathrm{1}},\mathrm{\dots}{y}_{N})$. y_{i} denotes line intensities simulated for the individual species. In the literature, P(θY) is typically referred to as the posterior.
π(θ) represents the probabilities of the parameters, the prior. Finally, π(Y), typically referred to as the evidence, represents the probabilities of the simulated line profiles.
ℒ(Yθ), the likelihood, is the product of the probability of the individual data points,
given an independent set of data points $Y=({y}_{\mathrm{0}},{y}_{\mathrm{1}},\mathrm{\dots}{y}_{N})$ sampled over $X=({x}_{\mathrm{0}},{x}_{\mathrm{1}},\mathrm{\dots}{x}_{N})$.
Assume we have a parameterized model f(xθ) which models the data Y for given parameters θ and explanatory variable X, in our case the wavelength. With this model we can then calculate the probability of generating an observed data point P(yθ), also known as the likelihood ℒ(yθ). If we have an independent set of data points $Y=({y}_{\mathrm{0}},{y}_{\mathrm{1}},\mathrm{\dots}{y}_{N})$ sampled over all wavelength $X=({x}_{\mathrm{0}},{x}_{\mathrm{1}},\mathrm{\dots}{x}_{N})$, the total likelihood is the product of the probability of the individual data points.
As the measured line profiles from the gaseous species are derived from CCD images, and thus photon counts, the likelihood of the probability distribution line profiles over wavelength is given by a Poisson distribution.
Instead of using the likelihood as in Eq. (18), it is common to work with the loglikelihood, as that reduces the product to a sum over the data points. In the Poisson case, the negative loglikelihood is given by
Before we can start the Bayesian framework, an initial guess of the synthetic spectrum is computed based on the values in Table 3. As we study the optically thin case (optical depth τ≪1) and use discreet temperature values, the final synthetic spectrum is then a linear combination of the individual emission spectra of each species as described in Eq. (20):
where ${F}_{\mathit{\nu}}^{\text{source}}$ is the flux vector of the complete source spectrum, i.e., the sum of the spectra of all species “i” at the source, N_{i} the number density of each species, and j_{ν,i} the temperaturedependent emission coefficient of element i. The elements' emission coefficients are retrieved from PARADE (introduced in Sect. 7.1). The flux vector F_{ν} comes from the meteor observations. The only unknowns of the equation are the number densities N_{i}.
The estimation of the number densities of the meteor's radiating elements is articulated in three sequential steps: a linear regression retrieval, a nonlinear leastsquares initial solution, and the final Markov chain Monte Carlo (MCMC) parameter estimation.
For each species and each discreet temperature in our emission database (described in Sect. 7.1), the obtained spectral profile is further modeled into a set of synthetic spectra. These are obtained through different combinations of the parameters in Table 3, left free to vary within defined boundaries (“Prior” in Table 3). A linear solution of the parameters in Table 3 is then retrieved, estimating the Moore–Penrose pseudoinverse of the matrix having on its rows the species' synthetic spectra obtained at each combination. By applying a linear regression, considering the ICC8 observed meteor spectrum as the dependent variable, a first estimate of the number density of all the species considered is obtained.
Using these linear estimates as a first guess, a nonlinear leastsquares approach is used to obtain an even more accurate estimation of the number densities' initial solution. This provides the estimate of all the unknown parameters in Table 3, whose combination generates a complete synthetic meteor which best fits the ICC8 observed meteor spectrum. The nonlinear leastsquares solution for the number densities, and other parameters, is determined using the Levenberg–Marquardt algorithm. Our approach is very similar to Benneke and Seager (2012), where Bayesian inference is used to retrieve and constrain the composition of exoplanet atmospheres from their transit spectra. These bestfit leastsquares estimates are, however, only an initial solution and thus do not coincide with the final bestfit solution obtained at the end of the routine, after the Markov chain Monte Carlo (MCMC).
The leastsquares parameter estimates are used as initial guesses for the successive MCMC. Using the MCMC sampling we estimate the posterior distribution of the parameters of the spectral model. The MCMC sampler approximates the posterior distribution by generating samples with a probability that is proportional to the posterior:
where n(θ)dθ is the number of samples between θ and θ+dθ. In MCMC methods we generate socalled Markov chains, which follow this distribution. A Markov chain is a random process which obeys the Markov property that a successive data point is only based on the most recent data point.
We thus obtain not only a point estimate of the chemical composition of the meteor but also an estimate of the full probability distribution, giving a constraint on the range of abundances which could describe the meteor. During the MCMC, for each meteor event 1000 iterations are run. Each iteration models a different synthetic spectrum using parameters which are sampled around the nonlinear leastsquares solution via MCMC affineinvariance sampling. We use the emcee Python package (ForemanMackey et al., 2013, 2019) that implements the affineinvariant ensemble sampler by Goodman and Weare (2010). The algorithm performed well under all linear transformations and is insensitive to covariances among parameters. The complete synthetic spectrum which best fits the ICC8 observed spectrum is inferred using Bayesian inference: the bestfit parameters represent those for which the posterior has the highest value. Specifically, the inference returns the combination of parameters for which the modeled (synthetic) spectra best fit the observed spectra (data).
This section describes the individual processing steps and intermediate results based on the meteor event 20131214T041302. The meteor presented was classified by METREC to belong to the Geminid shower and was observed in December 2013. The meteor itself and its spectral signature were detected in 24 individual frames in cameras ICC7 and ICC8, respectively. The maximum apparent magnitude was −1.7. The level1 FITS file contains the individual frames and total image (shown in Fig. 5).
In the next processing step, the (x, y) position for each wavelength is determined frame by frame, resulting in an line plot representing the meteor spectrum intensity. The profile is further spectrally calibrated and corrected for atmospheric effects, resulting in a level2A FITS file. The level2B FITS file contains the final calibration product, after the intensity calibration is applied, in radiance (W m^{−2} nm^{−1}). The unprocessed and calibrated profile of the brightest frame, frame 20, is shown in Fig. 6.
The 1000 MCMC iterations of the ablation processes as part of the Bayesian interference algorithms are based on the priors and initial settings as given in Table 3 (see Sect. 7.2). In our setup – thus also applicable for this example – the initial settings used for the chemical elements are computed by a leastsquare fitting algorithm, and the resulting values are given in Table 4. The selection of chemical elements, the main temperature, and the secondary temperature (see Sect. 7.2) proved stable through the processing of all meteor events. In the future, we might apply a different larger selection of elements, i.e., when processing meteor events of a dedicated meteor shower.
After the Bayesian interference algorithm and MCMC run have finished simulations and the best fit is available, a check on the convergence of all parameters used in the simulation is needed. This is done by the visual analysis of the socalled corner plot that is produced by the data pipeline for each meteor event. The corner plot – see Fig. 7 in ForemanMackey (2016) – visualizes the posterior probability distribution of the parameters in diagonals and the covariances between the parameters of the Monte Carlo iterations on the nondiagonals. The distribution of most parameters indicates a Gaussian distribution around a best estimate. A clear peak is visible for MgI, but two side maxima are visible. This might indicate either a higher level of noise or interference with another chemical element. From Fig. 10, we could argue that CrI could interfere with MgI.
This validation step is currently the only manual step needed in the calibration pipeline. In the case that a parameter is not converging, the parameter needs to be taken out of the simulation and a rerun of the simulation is needed.

The two temperature regimes result in temperatures of 4000 and 9000 K. The variations in the chain plots are marginal and we consider the two temperatures to be an important result of our analysis.

The chemical elements Fe, Na, Mg, Si, Cr, Ti, K, CaO, FeO, N_{2}, Ca^{+}, Fe^{+}, MG^{+}, O, and N converge well.

The parameters background, wavelength shift, and radius_sphere are also contained in the corner plots (not shown in Fig. 7). The background is a proxy for the noise level in the spectral profiles, and the wavelength shift is an indicator for an additional correction of the wavelength to achieve a good result (in the presented case fewer than two pixels); radius_sphere is the radius of the meteor's luminous area. According to Jenniskens (2007), the meteor is considered an isotropic sphere; i.e., the radiating area perpendicular to the observer has the same value across all directions of the observer, and thus the luminous area equals π⋅radius_sphere^{2}.
At the end of our event analysis, the data allow the comparison of the observed spectrum to the simulated spectra obtained for the total images and the brightest frame. Figure 8 shows the obtained plots demonstrating the fitting algorithm. Most maxima in the observed spectrum are well fitted. Also, the broad bands are fitted well, but some overshooting is observed between 550 and 650 nm. The intensity of the individual frame is given on the right y axis and is less intense than in the total image. The reader shall note that the red line in Fig. 8 represents the calibrated spectra (W m^{−2} nm^{−1}) of the brightest frame and is thus much smaller than the corresponding black line that represents the energy of the overall meteor event (thus of all the frames).
Table 4 contains the resulting number densities in the right column that can be interpreted in the following way.

From the meteoric elements, the highest number density is obtained for Si.

The Na content is relatively small, especially in comparison to the Fe and Mg. As the sodium is a direct indicator of the fluffiness of the meteoroid and thus its parent body (Abe et al., 2020), this is an indication that the meteoroid might be composed of stronger material without small grains. Our observation confirms earlier studies from, e.g., Abe et al. (2020) and Vojáček et al. (2019).

Borovička (2010) and Abe et al. (2020) provided ratios of $\mathrm{Na}/\mathrm{Mg}$ of 0.023 and 0.067, respectively. The value we obtain for the event discussed is 0.0003 and deviates thus a lot.

Borovička (2010) and Abe et al. (2020) provided ratios of $\mathrm{Fe}/\mathrm{Mg}$ of 0.32 and 0.3, respectively. The value we obtain is 0.09, a factor of 3 smaller. Borovička (2010) and Abe et al. (2020) mention, however, that the values of individual meteor events differ largely, by a factor of 2 for $\mathrm{Fe}/\mathrm{Mg}$ and a factor of more than 10 for $\mathrm{Na}/\mathrm{Mg}$. Further investigations are needed to explain the large difference of the CILBOderived ratios from the ones provided in the cited literature.

$\mathrm{Mg}/\mathrm{Si}$ and $\mathrm{K}/\mathrm{Si}$ ratios of 0.046 and 0.005, respectively, are in the ranges of the COSIMA instrument on board Rosetta (Rubin et al., 2020).

The $\mathrm{Na}/\mathrm{Si}$ ratio of $\mathrm{1.8}\times {\mathrm{10}}^{\mathrm{5}}$ is, however, much smaller than the range given for COSIMA, which is 0.02 to 0.2 (Bardyn et al., 2017).

The $\mathrm{Ca}/\mathrm{Si}$ ratio of $\mathrm{8.6}\times {\mathrm{10}}^{\mathrm{6}}$ is also much smaller than the range given for COSIMA. It must be noted that the COSIMA instrument analyzed the bulk refractory material, which might explain some of the different ratios obtained. Madiedo et al. (2013) state that Ca might not fully vaporize during the ablation, which might also explain the low $\mathrm{Ca}/\mathrm{Si}$ ratio.

Cr is observed with high abundance, i.e., $\mathrm{Cr}/\mathrm{Na}$ = 3.1 and $\mathrm{Cr}/\mathrm{Ca}$ = 6.6. Chromium measurements were reported earlier by Russell et al. (1956) and recently by Matlovič et al. (2020). The measurements of Matlovič et al. (2020) were obtained from laboratory measurements on meteorite samples, and the high Cr value was attributed to the high values of daubreelite (FeCr2S4) and lower values of enstatite (MgSi3O).

No Ca and Fe ions are detected in the T=9000 K component. It is thus more surprising that the Mg ion number density is relatively large. Mg ions were reported back in Etheridge and Russell (1968). The observations discussed in Jenniskens et al. (2002) indicate that Ca, Fe, and Mg ions are a result of the hot T=10 000 K component, and one would thus expect to identify Ca, Fe, and Mg ions in the event. The detailed interpretation of these observations first needs to be confirmed by further detailed analysis of other events and supplemented by a better interpretation of the chemical ablation physics.
Table 5 summarizes the data processing levels used in the CILBO spectral calibration pipeline and in the archive.
As the processing is very intensive with respect to time and computing resources, the synthetic spectrum was computed for the total image of each event. This is reflected in the subdirectories and in the names of the FITS structures. The data processing levels 1 to 3 are made publicly available via the Guest Archive Facility of ESA's Planetary Science Archive (Besse et al., 2018). All data products are stored in the FITS format (Pence et al., 2010). Level 0 represents the raw, unmodified data; level1 data are the radiometrically calibrated images, and level2 data represent the spectral profiles after radiometric, camera spectral sensitivity, and altitudedependent atmospheric extinction corrections have been applied. Also, the response to the intensity calibration correction is part of the level2 data product. The level3 data contain the derived parameters obtained from the Bayesian framework and MCMC processing: the number densities of the chemical elements, several ablation parameters, and the fitted spectra. Level3 data are also accompanied by browse images in PNG format.
9.1 Level1 data
The level1 FITS files are radiometrically calibrated, and thus dark current subtraction and a flatfield multiplication correction have been applied. The file names are formatted as, e.g., CC8_20131214T041392_lev1.fits, combining all information of one particular meteor event into a single file. The primary header provides information on the observation, the camera, and the meteor. The primary data are the total image of the event, and thus each pixel represents the maximum pixel value of all frames at this position. The total image allows obtaining a quick indication of the brightness, completeness, and quality of the images available. Keyword–value pairs in the FITS header represent information on the observation (date, time, observer, location), camera details (name, focal length, manufacturer), or the meteor event itself (brightness, meteor shower identification). Each FITS data product has four extensions, each being a binary table. The individual tables are structured in a similar way starting with wavelength information, information on each frame, and finally the information on the total image.

ICC8FRAMESINFO contains for each frame in the meteor event information on the frame number, the Julian date, right ascension and declination in degrees (of the zeroorder location obtained from ICC7), the apparent magnitude, the measured signaltonoise ratio, a wavelength error, the number of CCD rows binned together, the distance between the camera and the meteor, the height of the meteor above the ground, and the intensity calibration coefficient. In the case that a parameter cannot be calculated, the “nan” value is used, which is, e.g., the case when there are no doublestation observations available and the height cannot be computed.

The algorithm has computed, based on the zeroorder information available from ICC7, the location (x and y position in frame) and the corresponding wavelength for each pixel coordinate. The pixel values of these (x, y) positions are read and combined into an array, or simply profile:

ICC8PROFILES contains for each frame the derived profile based on the zeroorder information.

ICC8WAVELENGTH contains for each frame and for each element in the profile the corresponding wavelength. The wavelengths are provided in meters.

ICC8SENSITIVITY contains the normalized instrument sensitivity curve. See Sect. 5 for further information.
9.2 Level2 data
The level2 FITS profiles are corrected for the spectral response of the camera (using the ICC8SENSITIVITY information) and corrected for the altitudedependent atmospheric correction (see Sect. 5). The result of these corrections is stored as level 2A. The result of the intensity calibration of the level2A profiles is stored as level2B profiles. Whereas level2A data are dimensionless, level2B data are provided in watts per square meter per nanometer (W m^{−2} nm^{−1}).

ICC8FRAMESINFO is identical to level 1.

ICC8LEVEL1 is identical to ICC8PROFILES from level 1.

ICC8LEVEL2A has an identical structure as in ICC8LEVEL1; the camera spectral calibration is applied to ICC8LEVEL1 profiles using the first array in the ICC8SPECTRALCAL extension as calibration curves.

ICC8LEVEL2B has an identical structure as in ICC8LEVEL1; the intensity (absolute) calibration is applied to ICC8LEVEL2B profiles.

ICC8SPECTRALCAL contains the ICC8 spectral calibration curve followed by the atmospheric correction factors per wavelength for each frame.
9.3 Level3 data
The level3 FITS products contain the simulation results of the Bayesian method applied to the level2 data products based on PARADE. The primary header repeats the general information on the meteor event (date, observation information, camera properties, and meteor shower information), a copy from the level1 and level2 data products. The two extension tables contain the gas species information for each species and considered frame.

TOT_IMAGE_BESTFIT contains the following for each chemical element analyzed with the Bayesian framework:

SPEC_TYPE that is either main or second,

NUM_DENSITY of floating type number in scientific notation (m^{−3}),

information on the temperature regime is given in the first two rows as two temperature regimes were analyzed, and

the Gaussian broadening, the background noise, the wavelength shift applied, and the spherical radius independent of the elements and provided in the first row only.


TOT_IMAGE_SPECTRA contains four arrays representing the wavelength, the observed spectral profile, the synthetic spectral profile, and the spectral response curve used. From these arrays the spectral profile plots can be regenerated.
For validation purposes the archive contains a subdirectory named “plots”, containing for each meteor event three different images in PNG format:

eventname_lev3_fit.png, a plot of the observed spectrum overlaid with the synthetic spectrum – see Fig. 9;

eventname_lev3_elementsfit.png, containing two plots of the observed spectrum, with the synthetic spectra of each chemical element overlaid (see Fig. 10); and

eventname_lev3_chain.png, containing the MCMC chain results from each chemical element used in the PARADE simulation – see Sect. 7.1 and Fig. 7.
We present an algorithm to process meteor events acquired using a camera system with a spectral grating attached to it from raw format into quantitative chemical information. The algorithm is based on the principle that the location of a specific diffracted wavelength in first order can be directly derived from the knowledge of the zeroorder position of the meteor. In our setup, we obtain the zeroorder information from an imageintensified camera that is mounted with an offset angle to another intensified camera with a transmission grating mounted on top. The presented algorithm includes the radiometric calibration of the images, the determination of the location of the firstorder information for each wavelength, the extraction of the meteor event spectral profile, and its spectral calibration. The spectrum derived is further corrected for atmospheric effects that take place between the observer and the location of the meteor ablation at the different heights during the event. In the next step, a Bayesian approach is used to map a simulated spectrum against the reduced spectrum obtained by our algorithm. To simulate the ablation we use the PARADE tool that was extended to allow the simulation for most meteoric and atmospheric atoms, diatomic molecules, ions, oxides, and a few complex molecules. As preknowledge for the Bayesian framework we use the assumption of two main temperature regimes during the ablation, as well as the typical chemical elements of meteors and their respective line intensities. After running typically 1000 Markov chain Monte Carlo iterations, we obtain the simulated spectrum as overall line intensities, but also obtain the line intensities of each chemical element. In a last processing step, the line intensities are calibrated themselves, resulting in the number densities of each chemical element.
Several hundred meteor events have been observed by CILBO in the past decade, and we provide full access to the data at their raw and intermediate calibrated data processing levels. The data are stored using the FITS standard, and an overview of the keyword–value pairs in the FITS header, as well as a description of the extended binary arrays used in the FITS extensions, is provided.
A discussion of one meteor event, its processing chain, and resulting spectral information is given.
The presented algorithm works autonomously – that means without user interaction. The processing time is, however, large, and even when using a multiprocessor, multicore system, we have not yet computed all individual frames of all meteor events. Instead, we have processed for each meteor event the information obtained from the total image as well as from the brightest frame. In the dataset published, we share the data representing the total image.
All events contained in the dataset have been checked for the convergence of the Bayesian parameters as a result of the MCMC iterations. The list of chemical elements used in the Bayesian approach was not extensive, and in the processing so far, we only checked for the typical chemical elements also reported by other observers in the past. It is our goal to reprocess all the data and use more chemical elements than PARADE is currently able to simulate. Such a step, however, needs a very careful study of the convergence of the chemical elements – which is still a manual step thus requiring dedicated efforts.
The meteor spectra data from CILBO are freely available at the Guest Archive Facility of ESA/ESAC: https://doi.org/10.57780/esatfnxgu7 (Zender and Koschny, 2023).
The activities reported in this paper were led and coordinated by JJZ, DVK, and RR. Major contributions to the software development were made by SV and RL. SL led the efforts in the extension of the PARADE simulator. LM gave advice on ablation physics and SL on the implementation of oxides and other chemical elements into PARADE. HS provided continuous software support and the installation of the CILBO computing system. DS gave advice on the correction of the atmospheric effects.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The work presented in this paper is the result of efforts over a 10year period executed on a continuous but timelimited basis. Many students supported parts of this activity, and each of them contributed to this final result presented. All of these efforts have either been funded by the Science Faculty of ESA's Science Directorate or through studentships provided for by ESA. Special thanks go to Cosette Molijn, Norika Bartholomae, Kevin Ravensberg, Imanol Uiarte, and Joe Lightfoot. The contributions from Jorge Diaz del Rio, Rico Landman, and Lionel Maraffa were highly valuable and appreciated. We would also like to thank Wouter van der Wals and Daphne Stams for their important advice on atmospheric physics and coleading students. We would like to thank the reviewers for the critical, constructive review, which led to many improvements of the paper.
This paper was edited by Jean Dumoulin and reviewed by Juraj Tóth.
Abe, S., Ogawa, T., Maeda, K., and Arai, T.: Sodium variation in Geminid meteoroids from (3200) Phaethon, Planet. Space Sci., 194, 105040, https://doi.org/10.1016/j.pss.2020.105040, 2020. a, b, c, d, e, f
Altwegg, K., Balsiger, H., Hänni, N., Rubin, M., Schuhmann, M., Schroeder, I., Sémon, T., Wampfler, S., Berthelier, J.J., Briois, C., Combi, M., Gombosi, T. I., Cottin, H., De Keyser, J., Dhooghe, F., Fiethe, B., and Fuselier, S. A.: Evidence of ammonium salts in comet 67P as explanation for the nitrogen depletion in cometary comae, Nat. Astro., 4, 533–540, https://doi.org/10.1038/s4155001909919, 2020. a
Appenzeller, I.: Introduction to Astronomical Spectroscopy, Cambridge University Press, 145–178, ISBN 1139779524, 9781139779524, 2013. a
Bardyn, A., Baklouti, D., Cottin, H., Fray, N., Briois, C., Paquette, J., Stenzel, O., Engrand, C., Fischer, H., Hornung, K., Isnard, R., Langevin, Y., Lehto, H., Le Roy, L., Ligier, N., Merouane, S., Modica, P., OrthousDaunay, F.R., Rynö, J., Schulz, R., Silén, J., Thirkell, L., Varmuza, K., Zaprudin, B., Kissel, J., and Hilchenbach, M.: Carbonrich dust in comet 67P/ChuryumovGerasimenko measured by COSIMA/Rosetta, Mon. Not. R. Astron. Soc., 469, S712–S722, https://doi.org/10.1093/mnras/stx2640, 2017. a
Benneke, B. and Seager, S.: Atmospheric retrieval for superearths: Uniquely constraining the atmospheric composition with transmission spectroscopy, Astrophys. J., 753, 100, https://doi.org/10.1088/0004637X/753/2/100, 2012. a
Besse, S., Vallat, C., Barthelemy, M., Coia, D., Costa, M., De Marchi, G., Fraga, D., Grotheer, E., Heather, D., Lim, T., Martinez, S., Arviset, C., Barbarisi, I., Docasal, R., Macfarlane, A., Rios, C., Saiz, J., and Vallejo, F.: ESA's Planetary Science Archive: Preserve and present reliable scientific data sets, Planet. Space Sci., 150, 131–140, https://doi.org/10.1016/j.pss.2017.07.013, 2018. a
Borovicka, J.: Two components in meteor spectra, Planet. Space Sci., 42, 145–150, https://doi.org/10.1016/00320633(94)900256, 1994a. a, b, c
Borovicka, J.: Line identifications in a fireball spectrum, A&AS, 103, 83–96, 1994b. a
Borovička, J.: Spectroscopic Analysis of Geminid Meteors, in: Proceedings of the International Meteor Conference, 26th IMC, Bareges, France, 2007, 42–51, 2010. a, b, c
Borovička, J., Koten, P., Spurný, P., Boček, J., and Štork, R.: A survey of meteor spectra and orbits: evidence for three populations of Nafree meteoroids, Icarus, 174, 15–30, https://doi.org/10.1016/j.icarus.2004.09.011, 2005. a
Brownlee, D.: The Stardust Mission: Analyzing Samples from the Edge of the Solar System, Annu. Rev. Earth Pl. Sci., 42, 179–205, https://doi.org/10.1146/annurevearth050212124203, 2014. a
Ceplecha, Z.: Determination of wavelengths in meteor spectra by using a diffraction grating, B. Astron. I. Czech., 12, 246, 1961. a, b
Ceplecha, Z.: Spectral data on terminal flare and wake of doublestation meteor No. 38421 (Ondřejov, April 21, 1963), B. Astron. I. Czech., 22, 219, 1971. a
Dubs, M. and Schlatter, P.: A practical method for the analysis of meteor spectra, WGN, Journal of the International Meteor Organization, 43, 94–101, https://doi.org/10.48550/arXiv.1509.07531, 2015. a, b
Etheridge, D. A. and Russell, J. A.: The Spectrum of a Meteor from the 1966 Leonid Shower, Publ. Astron. Soc. Pac., 80, 550, https://doi.org/10.1086/128676, 1968. a
ForemanMackey, D.: corner.py: Scatterplot matrices in Python, The J. Open Source Softw., 1, 24, https://doi.org/10.21105/joss.00024, 2016. a
ForemanMackey, D., Hogg, D. W., Lang, D., and Goodman, J.: emcee : The MCMC Hammer , Publ. Astron. Soc. Pac., 125, 306–312, https://doi.org/10.1086/670067, 2013. a
ForemanMackey, D., Farr, W., Sinha, M., Archibald, A., Hogg, D., Sanders, J., Zuntz, J., Williams, P., Nelson, A., de ValBorro, M., Erhardt, T., Pashchenko, I., and Pla, O.: emcee v3: A Python ensemble sampling toolkit for affineinvariant MCMC, J. Open Source Softw., 4, 1864, https://doi.org/10.21105/joss.01864, 2019. a
Goodman, J. and Weare, J.: Ensemble samplers with affine invariance, Comm. App. Math. Com. Sc., 5, 65–80, https://doi.org/10.2140/camcos.2010.5.65, 2010. a
Hadamcik, E., LevasseurRegourd, A. C., Hines, D. C., Sen, A. K., Lasue, J., and Renard, J. B.: Properties of dust particles in comets from photometric and polarimetric observations of 67P, Mon. Not. R. Astron. Soc., 462, S507–S515, https://doi.org/10.1093/mnras/stx030, 2016. a
Herzberg, G.: Molecular spectra and molecular structure, Vol.1: Spectra of diatomic molecules, Van Nostrand, Princeton, 212–240, ISBN 10 0442033850, ISBN 13 9780442033859, 1950. a, b
Jenniskens, P.: Quantitative meteor spectroscopy: Elemental abundances, Adv. Space Res., 39, 491–512, https://doi.org/10.1016/j.asr.2007.03.040, 2007. a, b, c, d, e, f, g
Jenniskens, P., Tedesco, E., Murthy, J., Laux, C., and Price, S.: Spaceborne ultraviolet 251384 nm spectroscopy of a meteor during the 1997 Leonid shower, Meteorit. Planet Sci., 37, 1071–1078, https://doi.org/10.1111/j.19455100.2002.tb00878.x, 2002. a
Jenniskens, P., Gural, P., and Berdeu, A.: CAMSS: A spectroscopic survey of meteoroid elemental abundances, in: Meteoroids 2013, edited by: Jopek, T. J., Rietmeijer, F. J. M., Watanabe, J., and Williams, I. P., A.M. University Press, 117–124, ISBN 9788323227267, 2014. a
Koschny, D. and Diaz del Rio, J.: Meteor Orbit and Trajectory Software (MOTS) – Determining the Position of a Meteor with Respect to the Earth Using Data Collected with the Software MetRec, WGN, J. Int. Meteor Organization, 30, 87–101, 2002. a
Koschny, D., Bettonvil, F., Licandro, J., Luijt, C. v. d., Mc Auliffe, J., Smit, H., Svedhem, H., de Wit, F., Witasse, O., and Zender, J.: A doublestation meteor camera setup in the Canary Islands – CILBO, Geosci. Instrum. Method. Data Syst., 2, 339–348, https://doi.org/10.5194/gi23392013, 2013a. a
Koschny, D., Bettonvil, F., Licandro, J., v. d. Luijt, C., Mc Auliffe, J., Smit, H., Svedhem, H., de Wit, F., Witasse, O., and Zender, J.: A doublestation meteor camera setup in the Canary Islands – CILBO, Geosci. Instrum. Method. Data Syst., 2, 339–348, https://doi.org/10.5194/gi23392013, 2013b. a, b
Koschny, D., Mc Auliffe, J., Drolshagen, E., Bettonvil, F., Licandro, J., van der Luijt, C., Ott, T., Smit, H., Svedhem, H., Witasse, O., and Zender, J.: CILBO – Lessons learned from a doublestation meteor camera setup in the Canary Islands, in: Proceedings of the International Meteor Conference, Giron, France, 18–21 September 2014, edited by: Rault, J. L. and Roggemans, P., 10–15, ISBN 9782873550288, 2014. a
Koschny, D., Albin, T., Drolshagen, E., Drolshagen, G., Drolshagen, S., Koschny, J., Kretschmer, J., van der Luijt, C., Molijn, C., Ott, T., Poppe, B., Smit, H., Svedhem, H., Toni, A., de Wit, F., and Zender, J.: Current activities at the ESA/ESTEC Meteor Research Group, in: International Meteor Conference Mistelbach, Austria, p. 204, ISBN 9782873550295, 2015a. a
Koschny, D., Albin, T., Drolshagen, E., Drolshagen, G., Drolshagen, S., Koschny, J., Kretschmer, J., van der Luijt, C., Molijn, C., Ott, T., Poppe, B., Smit, H., Svedhem, H., Toni, A., de Wit, F., and Zender, J.: Current activities at the ESA/ESTEC Meteor Research Group, in: International Meteor Conference Mistelbach, Austria, p. 204, ISBN 9782873550295, 2015b. a
Koschny, D., Drolshagen, E., Drolshagen, S., Kretschmer, J., Ott, T., Drolshagen, G., and Poppe, B.: Flux densities of meteoroids derived from optical doublestation observations, Planet. Space Sci., 143, 230–237, https://doi.org/10.1016/j.pss.2016.12.007, 2017. a
Koschny, D., Soja, R. H., Engrand, C., Flynn, G. J., Lasue, J., LevasseurRegourd, A.C., Malaspina, D., Nakamura, T., Poppe, A. R., Sterken, V. J., and TrigoRodríguez, J. M.: Interplanetary Dust, Meteoroids, Meteors and Meteorites, Space Sci. Rev., 215, 34, https://doi.org/10.1007/s1121401905977, 2019. a
Koschny, D. V., McAuliffe, J., Bettonvil, F. C. M., Gritsevich, M., van der Luijt, C., Ocaña, F., Smit, H., Svedhem, H., and Zender, J. J.: What happened at ESA's Meteor Research Group in 2010/11?, in: Proceedings of the International Meteor Conference, 30th IMC, Sibiu, Romania, 2011, 57–60, ISBN 29782873550233, 2012. a
Lang, D., Hogg, D. W., Mierle, K., Blanton, M., and Roweis, S.: Astrometry.net: Blind Astrometric Calibration of Arbitrary Astronomical Images, The Astro. J., 139, 1782–1800, https://doi.org/10.1088/00046256/139/5/1782, 2010. a
LevasseurRegourd, A.C., Agarwal, J., Cottin, H., Engrand, C., Flynn, G., Fulle, M., Gombosi, T., Langevin, Y., Lasue, J., Mannel, T., Merouane, S., Poch, O., Thomas, N., and Westphal, A.: Cometary Dust, Space Sci. Rev., 214, 64, https://doi.org/10.1007/s1121401804963, 2018. a
Loehle, S., Eberhart, M., Zander, F., Meindl, A., Rudawska, R., Koschny, D., Zender, J., Dantowitz, R., and Jenniskens, P.: Extension of the plasma radiation database PARADE for the analysis of meteor spectra, Meteorit. Planet Sci., 56, 352–361, https://doi.org/10.1111/maps.13622, 2021. a, b, c, d
Madiedo, J. M.: Spectroscopy of a κCygnid fireball afterglow, Planet. Space Sci., 118, 90–94, https://doi.org/10.1016/j.pss.2015.03.015, 2015. a
Madiedo, J. M., TrigoRodríguez, J. M., CastroTirado, A. J., Ortiz, J. L., and CabreraCaño, J.: The Geminid meteoroid stream as a potential meteorite dropper: a case study, Mon. Not. R. Astron. Soc., 436, 2818–2823, https://doi.org/10.1093/mnras/stt1777, 2013. a
Matlovič, P., Tóth, J., Rudawska, R., Kornoš, L., and Pisarčíková, A.: Spectral and orbital survey of mediumsized meteoroids, A&A, 629, A71, https://doi.org/10.1051/00046361/201936093, 2019. a
Matlovič, P., Tóth, J., Kornoš, L., and Loehle, S.: On the sodium enhancement in spectra of slow meteors and the origin of Narich meteoroids, Icarus, 347, 113817, https://doi.org/10.1016/j.icarus.2020.113817, 2020. a, b
Molau, S.: The Meteor Detection Software METREC, in: Proceedings of the International Meteor Conference, 17th IMC, Stara Lesna, Slovakia, 1998, 9–16, ISBN 2873550104, 1999. a
Pence, W. D., Chiappetti, L., Page, C. G., Shaw, R. A., and Stobie, E.: Definition of the Flexible Image Transport System (FITS), version 3.0, A&A, 524, A42, https://doi.org/10.1051/00046361/201015362, 2010. a
Pfeiffer, B., Fertig, M., Winter, M., and Auweter Kurtz, M.: PARADE – a program to calculate the radiation of atmospheric reentry in different atmospheres, edited by: Warmbein, B., ESA Special Publication, 533, 85–91, ISBN 9290928433, 2003. a
Ralchenko, Y. and Kramida, A.: Development of NIST Atomic Databases and Online Tools, Atoms, 8, 56, https://doi.org/10.3390/atoms8030056, 2020. a
Rivilla, V. M., Drozdovskaya, M. N., Altwegg, K., Caselli, P., Beltrán, M. T., Fontani, F., van der Tak, F. F. S., Cesaroni, R., Vasyunin, A., Rubin, M., Lique, F., Marinakis, S., Testi, L., Rosina Team, Balsiger, H., Berthelier, J. J., de Keyser, J., Fiethe, B., Fuselier, S. A., Gasc, S., Gombosi, T. I., Sémon, T., and Tzou, C. Y.: ALMA and ROSINA detections of phosphorusbearing molecules: the interstellar thread between starforming regions and comets, Mon. Not. R. Astron. Soc., 492, 1180–1198, https://doi.org/10.1093/mnras/stz3336, 2020. a
Rubin, M., Engrand, C., Snodgrass, C., Weissman, P., Altwegg, K., Busemann, H., Morbidelli, A., and Mumma, M.: On the Origin and Evolution of the Material in 67P/ChuryumovGerasimenko, Space Sci. Rev., 216, 102, https://doi.org/10.1007/s11214020007182, 2020. a, b, c
Rudawska, R., Zender, J., Jenniskens, P., Vaubaillon, J., Koten, P., Margonis, A., Tóth, J., McAuliffe, J., and Koschny, D.: Spectroscopic Observations of the 2011 Draconids Meteor Shower, Earth Moon Planet., 112, 45–57, https://doi.org/10.1007/s1103801494368, 2014. a
Rudawska, R., Zender, J., Koschny, D., Smit, H., Löhle, S., Zander, F., Eberhart, M., Meindl, A., and Latorre, I. U.: A spectroscopy pipeline for the Canary island long baseline observatory meteor detection system, Planet. Space Sci., 180, 104773, https://doi.org/10.1016/j.pss.2019.104773, 2020. a
Russell, J. A., Sadoski, Michael J., J., Sadoski, D. C., and Wetzel, G. F.: The Spectrum of a Geminid Meteor, Publ. Astron. Soc. Pac., 68, 64, https://doi.org/10.1086/126881, 1956. a
Smith, A. J.: Plasma radiation database parade v2.2, Technical Report TR 28/96 Issue 3, Tech. rep., Emsworth, UK: Fluid Gravity Engineering Ltd., 2006. a, b
Snodgrass, C., A'Hearn, M. F., Aceituno, F., Afanasiev, V., Bagnulo, S., Bauer, J., Bergond, G., Besse, S., Biver, N., Bodewits, D., Boehnhardt, H., Bonev, B. P., Borisov, G., Carry, B., Casanova, V., Cochran, A., Conn, B. C., Davidsson, B., Davies, J. K., de León, J., de Mooij, E., de ValBorro, M., Delacruz, M., DiSanti, M. A., Drew, J. E., Duffard, R., Edberg, N. J. T., Faggi, S., Feaga, L., Fitzsimmons, A., Fujiwara, H., Gibb, E. L., Gillon, M., Green, S. F., Guijarro, A., GuilbertLepoutre, A., Gutiérrez, P. J., Hadamcik, E., Hainaut, O., Haque, S., Hedrosa, R., Hines, D., Hopp, U., Hoyo, F., Hutsemékers, D., Hyland, M., Ivanova, O., Jehin, E., Jones, G. H., Keane, J. V., Kelley, M. S. P., Kiselev, N., Kleyna, J., Kluge, M., Knight, M. M., Kokotanekova, R., Koschny, D., Kramer, E. A., LópezMoreno, J. J., Lacerda, P., Lara, L. M., Lasue, J., Lehto, H. J., LevasseurRegourd, A. C., Licandro, J., Lin, Z. Y., Lister, T., Lowry, S. C., Mainzer, A., Manfroid, J., Marchant, J., McKay, A. J., McNeill, A., Meech, K. J., Micheli, M., Mohammed, I., Monguió, M., Moreno, F., Muñoz, O., Mumma, M. J., Nikolov, P., Opitom, C., Ortiz, J. L., Paganini, L., Pajuelo, M., Pozuelos, F. J., Protopapa, S., Pursimo, T., Rajkumar, B., Ramanjooloo, Y., Ramos, E., Ries, C., Riffeser, A., Rosenbush, V., Rousselot, P., Ryan, E. L., SantosSanz, P., Schleicher, D. G., Schmidt, M., Schulz, R., Sen, A. K., Somero, A., Sota, A., Stinson, A., Sunshine, J. M., Thompson, A., Tozzi, G. P., Tubiana, C., Villanueva, G. L., Wang, X., Wooden, D. H., Yagi, M., Yang, B., Zaprudin, B., and Zegmott, T. J.: The 67P/ChuryumovGerasimenko observation campaign in support of the Rosetta mission, Philos. T. Roy. Soc. A, 375, 20160249, https://doi.org/10.1098/rsta.2016.0249, 2017. a
Vaubaillon, J., Koten, P., Margonis, A., Toth, J., Rudawska, R., Gritsevich, M., Zender, J., McAuliffe, J., Pautet, P.D., Jenniskens, P., Koschny, D., Colas, F., Bouley, S., Maquet, L., Leroy, A., Lecacheux, J., Borovicka, J., Watanabe, J., and Oberst, J.: The 2011 Draconids: The First European Airborne Meteor Observation Campaign, Earth Moon Planet., 114, 137–157, https://doi.org/10.1007/s1103801494555, 2015. a
Vicinanza, S., Koschny, D., Rudawska, R., Stam, D., van der Wal, W., and Zender, J.: Spectral Calibration of Meteors: An ElevationDependent Atmospheric Correction, WGN, J. Int. Meteor. Org., 49, 201–210, 2021. a
Vojáček, V., Borovička, J., Koten, P., Spurný, P., and Štork, R.: Properties of small meteoroids studied by meteor video observations, A&A, 621, A68, https://doi.org/10.1051/00046361/201833289, 2019. a, b
Ward, W.: Video meteor spectroscopic and orbital observations, 2015 April to 2016 April, J. Br. Astro. Assoc., 127, 217–222, 2017. a
Zender, J.: Line identifications in a fireball spectrum, A&AS, 103, 83–96, 1994. a
Zender, J. and Koschny, D.: CILBO Meteor Spectra Database, ESA [data set] https://doi.org/10.57780/esatfnxgu7, 2023. a, b
Zender, J., Koschny, D., and Ravensberg, K.: Calibration of spectral video observations in the visual: theoretical overview of the ViDAS calibration pipeline, in: Proceedings of the International Meteor Conference, Poznan, Poland, 22–25 August 2013, edited by: Gyssens, M., Roggemans, P., and Zoladek, P., 126–129, ISBN 9782873550257, 2014. a
Zender, J. J., Witasse, O., Koschny, D. V., Trautner, R., Knöfel, A., Diaz del Rio, J., and Quicke, G.: First results of spectroscopic measurements during the ESA Leonid campaign 2001, in: Asteroids, Comets, and Meteors: ACM 2002, edited by: Warmbein, B., vol. 500 of ESA Special Publication, 121–125, ISBN 9290928107, 2002. a
 Abstract
 Introduction
 The Canary Island Longbaseline Observatory: CILBO
 Radiometric correction and background subtraction
 Wavelength determination
 Spectral calibration
 Intensity calibration
 Extraction of line intensities of individual chemical elements
 Spectral calibration discussed for one meteor event
 Available datasets
 Summary
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Review statement
 References
 Abstract
 Introduction
 The Canary Island Longbaseline Observatory: CILBO
 Radiometric correction and background subtraction
 Wavelength determination
 Spectral calibration
 Intensity calibration
 Extraction of line intensities of individual chemical elements
 Spectral calibration discussed for one meteor event
 Available datasets
 Summary
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Review statement
 References