Issue |
A&A
Volume 630, October 2019
|
|
---|---|---|
Article Number | A127 | |
Number of page(s) | 20 | |
Section | Atomic, molecular, and nuclear data | |
DOI | https://doi.org/10.1051/0004-6361/201935816 | |
Published online | 07 October 2019 |
Main spectral features of meteors studied using a terawatt-class high-power laser
1
J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague 8, Czech Republic
e-mail: martin.ferus@jh-inst.cas.cz
2
Charles University in Prague, Faculty of Science, Department of Physical and Macromolecular Chemistry, Albertov 2030, 12840 Prague 2, Czech Republic
3
Valašské Meziříčí Observatory, Vsetínská 78, 75701 Valašské, Meziříčí, Czech Republic
4
Department of Radiation and Chemical Physics, Institute of Physics, Czech Academy of Sciences, Na Slovance 1999/2, 18221 Prague 8, Czech Republic
5
Institute of Plasma Physics, Czech Academy of Sciences, Za Slovankou 1782/3, 182 00 Prague 8, Czech Republic
e-mail: krus@ipp.cas.cz
6
National Technical University of Athens, School of Mining and Metallurgical Engineering, 9 Heroon Polytechneiou str., 15780, Zografou Athens, Greece
7
Czech Technical University in Prague, Jugoslávských partyzánů 1580/3, 160 00 Prague 6, Czech Republic
Received:
1
May
2019
Accepted:
17
August
2019
Context. Meteor spectra are commonly interpreted using data from databases and tables. Several studies have demonstrated very sophisticated calculations of elemental compositions of meteoroid bodies based on the computation of synthetic meteor spectra or on the spectral analysis of airglow plasma containing evaporated, atomized, and ionized meteoroid matter. However, considering accuracy, reliability of computations, lack of laboratory experimental data in this field, as well as the complicated physical structure of meteor plasma, such qualitative assignment or quantitative calculations are still extensively discussed in the scientific community. Even on the laboratory level, many studies have shown the high complexity of the acquisition and interpretation of the data that are recorded with techniques of emission spectroscopy that are in fashion and philosophy similar to the spectral analysis of meteor plasma, that is, detection and quantification of the elements that are ablated from complicated multicomponent matrices.
Aims. The current study is focused on the application of terawatt-class laser-induced breakdown spectroscopy (TC-LIBS) of real samples of chondritic meteorites. We recorded emission spectra with high resolution and high precision that contain spectral lines that are typical for real meteoric spectra. Experimental data were compiled in a form that is convenient for the meteoric spectra interpretation and calibration.
Methods. TC-LIBS was carried out by a high-power terawatt-class laser facility, the Prague Asterix Laser System (PALS). The spectra were simultaneously recorded by an echelle high-resolution spectrograph in the UV/VIS spectral ranges and by a low-resolution spectrograph that was used for real observation of meteor spectra. We also present calculated synthetic spectra based on data from the NIST atomic spectra database.
Results. We assembled etalon qualitative tables of major meteoric spectral features that can be used both for the spectral wavelength calibration of low-resolution observational instruments and for the exact interpretation of meteor spectra. The data are compared with real meteor spectra.
Key words: atomic data / meteorites / meteors / meteoroids
© ESO 2019
1. Introduction
Meteor observation and spectroscopy is a highly important astronomical discipline. To date, most observational data are interpreted theoretically, and experimental work is lacking in this field. An in-depth analysis of meteor spectra can be used for example for basic qualitative or quantitative elemental analyses and characterization of meteoroids and their parent bodies, i.e., asteroids and comets. Asteroids and comets are remnants of the protoplanetary disk that formed the planetesimals and planets. Therefore, meteorites allow us to discover details about the properties of this disk as well as its history, that is, the physical and chemical evolution of the solar system. However, even in the rare cases when a piece of meteorite is found, the fundamental problem of linking specific meteorites to their parent bodies remains (interplanetary matter, asteroids, or comet nuclei). Most bodies are evaporated and completely disintegrated during their descent, and their emission spectra, measured using spectrographs, are the only record of their chemical composition. When the trajectory is recorded alongside with the emission spectrum and if chemical composition is correctly interpreted, these data can provide detailed information about the chemical properties of its source in interplanetary space. Determination of the elemental composition provided by emission spectra analysis of the plasma formed above the surface of the evaporated object (meteor) is practically similar to laboratory laser-induced breakdown spectroscopy (LIBS) elemental spectral analysis. In our recent pilot studies, we focused on the application of the calibration-free LIBS method in the interpretation of meteor spectra and the calculation of the elemental composition of meteoroids (Ferus et al. 2017a, 2018a). The method has been identified as feasible, but relatively time-demanding computation is still necessary. Ciucci et al. (1999), Tognoni et al. (2010), Giacomo (2011), Dell’Aglio et al. (2014), Horzňáčková et al. (2014), Ozdín et al. (2015) and Takahashi et al. (2015) provided a LIBS chemical analysis of meteorites and identified many advantages of this method, such as the applicability for a real-time in situ analysis without any preceding treatment, preparation, or isolation. On the other hand, as concluded in Tognoni et al. (2010) and references therein, even if the LIBS method is operated under strictly laboratory controlled conditions, the emission intensity of a particular spectral line, which is the main factor that influences the quantitave analysis, depends not only on the physical parameters of the spectral transition and the quantity of the emitting element, but also on the matrix in which the element is embedded. This influence of the matrix means that calibration curves or matrix-matched standards are required, which are not available in some practical situations, including meteorite samples. A coincidence of spectral lines in multicomponent analysis also complicates or even renders impossible a qualitative analysis of the spectra when low-resolution data are used. Although several studies have been published by only a few experts, such as Borovička et al. (1999), Jenniskens (2007), and Madiedo et al. (2014, 2013a,b), it was Ferus et al. (2018b) who demonstrated that a qualitative and quantitative analysis of mainly high-quality emission spectra that are assigned to the most abundant class of ordinary chondrites is possible. Many spectra of meteors are still classified as “unusual” and difficult to interpret.
The simulation of the chemical and physical consequences analogous to plasma formed by hypervelocity meteoroids or asteroids entering the Earth’s atmosphere by terawatt-class laser introduces a novel experimental approach also in the research of meteor spectra. In the past two decades, our team demonstrated several experiments that mostly focused on the chemical transformations of the atmosphere or on the interaction with solid or liquid surfaces (Babánková et al. 2006a). Most studies have been focused on impact-induced synthesis of biomolecules (Šponer et al. 2016; Ferus et al. 2017b) such as canonical nucleobases in Ferus et al. (2012, 2014a,b, 2015, 2017b), sugars in Civiš et al. (2016a), and aminoacids in Civiš et al. (2004), or the transformation of atmospheric molecules on early terrestrial planets (Civiš et al. 2008) such as the formation or decay of prebiotic substances, for instance, formamide (Ferus et al. 2011), isocyanic acid (Ferus et al. 2018a), and the transformations of hydrogen cyanide (Ferus et al. 2017c), acetylene (Civiš et al. 2016b), methane (Civiš et al. 2017), or carbon monoxide (Civiš et al. 2008; Ferus et al. 2009). Pyrometric measurement of dielectric breakdown induced by high-power terawatt-class lasers in the gas phase has shown that airglow temperatures of 4500 K (Babánková et al. 2006b) are very close to the low-temperature component in meteor spectra that exhibits 3800–6000 K independent of the mass of the impacting body or its velocity (specifically in the ranges of 35 and 72 km s−1 and masses between 1025 g and 1 g, according to Jenniskens et al. 2004). A spectroscopic survey of high-power terawatt-class laser plasma interaction with samples of meteorites has not been explored spfar. Several studies that focused on the simulation of high-density energy events connected with asteroid or meteoroid bodies have been published. Zakharov (2003) concluded that among other phenomena such as barium releases from the Earth’s magnetosphere, collisionless deceleration of supernova remnants or related shock-wave generation in the interstellar medium, near-Earth anti-asteroid explosions could be also reproduced in the laboratory using laser experiments. Milley et al. (2007) simulated meteor luminosity through laser ablation of meteorites. The first comprehensive work proposing that laboratory-based laser ablation techniques can be used not only to study the size of the luminous region, but also to predict spectral features, estimate the luminous efficiency factor, and assess the role of chemically differentiated thermal ablation of meteoroids was published by Hawkes et al. (2007). Ebert et al. (2017) simulated the virtually instantaneous melting of target rocks during meteorite impacts. They discovered that the entropy changes for laser-melting of a sandstone and an iron meteorite correspond to a minimum impact velocity of 6 km s−1, inducing peak shock pressures of 100 GPa in the target.
Alongside with these studies, the simulation of space weathering together with the deflection of asteroids by lasers remains the most frequently studied research topic connected with application of lasers in this field (Park & Mazanek 2003). Recently, Aristova et al. (2018) studied impact physics of nuclear explosion on hazardous asteroids, Moroz et al. (1996) simulated using laser ablation optical effects of impact melting and repeated crystallization on asteroidal surfaces, and Kurahashi et al. (2002) conducted laser ablation laboratory simulation of space weathering focused on finding the source of the difference between reflectance spectra of ordinary chondrites and their parent bodies, i.e., S-type asteroids. Loeffler et al. (2008) studied the effect of the redeposition of impact-ejecta on mineral surfaces using laser ablation.
To facilitate the understanding of meteor spectra and in order to increase the precision of the qualitatively assigned emission, we employed a high-power terawatt-class laser that can generate in one shot large volumes of plasma under strictly defined conditions. Alongside with the more than 1000× higher power in comparison with classical laboratory lasers, our terawatt-class laser-induced breakdown spectroscopy (TC-LIBS) also provides a great advantage in evaporation of a significantly larger area of the meteorite specimen. This occurs in a manner entirely to that induced by the plasma that is generated during meteoroid descent. Spectra have been recorded in situ by modern high-resolution echelle spectrograph alongside with a low-resolution astronomical meteor spectrograph that was used for the observation of real meteors. The current study provides a catalog of emission lines that we recorded in our experiments and highlights how low-resolution data can be interpreted.
2. Experimental details
Samples of chondrite metorites (see Table 1) were ablated by sub-nanosecond laser pulses (pulse duration 300 ps, wavelength 1.315 μm) with an energy of 600 J generated by the high-power TW-class facility Prague Asterix Laser System (PALS; Jungwirth et al. 2001). The intensity of the laser radiation focused on the surface of a meteorite specimen reached 2.2 TWcm−2. The experimental setup is depicted in Fig. 1. TC-LIBS spectra were simultaneously recorded by a high-resolution echelle spectrograph and the astronomical meteor spectrograph that was employed in spectroscopic surveys of meteors in direct observation. The instrumental equipment and methods are described in the following two subsections. The TC-LIBS spectra of these chondritic meteorites were then averaged in order to obtain the appropriate signal-to-noise ratio and the qualitative average analysis that is representative for a spectrum of a chondrite meteorite. These data qualitatively provide spectral features that are expected in the emission from the meteor plasma. These spectra and their recorded dominant features were compared with the average synthetic spectrum calculated for a wide range of temperatures and averaged for the chemical composition of a chondrite in Table 1. Elemental composition of the particular chondrite meteorite studied in this work by TC-LIBS was examined using an energy dispersive analysis (EDS). The physical constants for this calculation were taken from the NIST (National Institute of Standards and Technology) database. We would like to highlight that the intention of this work is in the first place to provide a qualitative catalog of spectral features, their experimental measurement, and their comparison with synthetic spectra. A quantitative comparison of spectra recorded for a particular specimen that is ablated by high-power TW-class laser, but also by a standard laboratory Nd:YAG laser, an excimer laser and a femtosecond Ti:Sapphire TW-laser will be provided in a follow-up study.
Elemental composition of the meteorites used to calculate the synthetic spectra.
Fig. 1. Panel A: vacuum interaction chamber with ablation set-up. A photography of a spectrum is shown in the embedded picture (left). Panel B: detailed experimental set-up with a collimator of Echelle spectrograph. Panel C: low resolution astronomical camera for comparative measurement. Panel D: sample of a meteorite prepared for the TC-LIBS experiment (left), high-power laser induced dielectric breakdown inside the interaction chamber (right). |
2.1. Simulation of meteor plasma using the PALS laser facility
The plasma creation and radiation during the meteoroid descent in the atmosphere was simulated using the high-power laser PALS. The laser beam was focused by a plano-convex lens (f/2, f = 60 cm) on the target, which was positioned out of focus; the spot size was 1 cm in diameter. All the experiments were performed under pressure of 2 mbar, corresponding to the altitudes that are typical for the very early stages of ablation of the meteoroid body ablation (above 110 km).
2.2. Optical emission spectroscopy of laser-induced breakdown plasma
Complex spectra within the entire UV/VIS range have been recorded by the High Resolution Echelle Spectra Analyzer ESA 4000 (LLA Instruments GmbH, Germany). The optical analyzer unit enables a spatially and temporally resolved image at the lowest spectral intensities. The resolution ranges from a few pm in the range 200–780 nm with resolution of 0.005 nm (200 nm) to 0.019 nm (780 nm). Simultaneously, the spectra were recorded for direct calibration by our astronomical meteor spectrograph that is directly employed in the observation of real meteors for comparative measurements. The spectrograph is based on the high-resolution camera PointGrey Grasshoper3 GS3-U3-32S4M-C with high quantum efficiency (QE = 76%, 525 nm), dynamical range (71.34 dB), with CMOS chip Sony Pregius 2048 × 1536 px, and a grating of 1000 lines mm−1, which allows a resolution of 0.48 nm px−1. The spectral intensity recorded by the instrument is calibrated using standard sources, such as a deuterium lamp, a tungsten source and calibration using standard spectra of Venus. Wavelength calibration is achieved using high-resolution data and standard wavelengths of calibration sources (deuterium lamp). In our measurements, for every meteorite sample, the echelle spectrograph was set to trigger 1 μs after the laser pulse with the gate open for 4 μs for a total accumulation of one signal after one large laser shot. The PALS high-power laser provides a single 600 J-laser shot every 30 min. The low-resolution astronomical meteor spectrograph opened the gate for 1 s after it was triggered by the laser system.
2.3. Bulk elemental analysis of chondrite samples
Chemical analyses of meteorites were provided by secondary and backscatter electron imaging (SEI/BSE) conducted with the JEOL JSM-6380 LV SEM system at the National Technical University of Athens, Greece. This system is equipped with the Energy Dispersive X-Ray Analysis (EDX) from Oxford Instruments, and it is controlled by the INCA software. Prior to any quantitative chemical analyses the instrument was calibrated using a set of mineral standards. Analyses were performed in high vacuum, with an electron beam size of about 1 μm and accelerating voltage of 20 kV. The beam current was generally set between 1–2 nA, and the analytical integration time was 120 s. The bulk analysis was performed by acquiring an X-ray spectrum while scanning and averaging an area of 100 mm2.
3. Calculating the synthetic spectra
The measured high-resolution emission spectra were compared with the synthetic spectra calculated from high-precision data obtained from the NIST atomic spectral database (Kramida et al. 2018). This approach greatly facilitates the orientation in the experimental data because the user can visually quickly compare the measured spectra with the standard spectra from the database; in addition to this advantage, the positions of the simulated peaks can be also used as standard values for the wavelength calibration of the observed spectra. This subsection describes how the synthetic spectra were obtained.
The synthetic spectra were generated as a sum of the individual peaks corresponding to the energy transitions that are included in the simulation. The intensity of each peak was calculated according to the equation
where Iul, h, c, λul, Aul, gu, Ni, Eu, k, T, and Q(T)i are the emission intensity, Planck’s constant, the speed of light in vacuum, the transition wavelength, the Einstein A-coefficient of the transition, the g-factor of the upper energy level, the number of the particles in ionization state i, the energy of the upper level, the Boltzmann constant, temperature, and the partition function of the ionization state, respectively. The simulation assumed a local thermodynamic equilibrium and optically thin conditions where no self-absorption effects are present. The Ni value was calculated as the product of the abundance of the specific element and the relative abundance of the ionization states produced by ionization of the neutral atoms of the element. The relative concentrations of the ionization states were calculated according to the Saha equation. Synthetic spectra were calculated based on elemental abundances measured in our laboratory for several samples of chondritic meteorites by means of energy dispersive analysis (EDS). These specimens were also used for ablation by the PALS laser. The average elemental composition is shown in the Table 1 and compared with values from the current literature.
The values of the electron density were set to 1 × 1014 cm−3, which is a typical value for the meteoric plasma (Borovička 1994). The synthetic spectra are provided in three different temperatures: 4000 K, which is typical for the lower temperature component of the meteoric spectra, 7000 K, which is often observed in our ablation experiments, and 10 000 K, which is near the value estimated for the so-called higher temperature component of the meteoric spectra (Borovička 1994). The maximum ionization degree of the species considered in this simulations was 1 (e.g., Ca II). Because the spectral resolution of most of the observation cameras used today is still lower than 1 nm px−1, the observed spectral peak profiles are still strongly influenced by the instrumental function of these spectrometers. For this reason, the simplest possible (Gaussian) peak profile was chosen as the profile function. The FWHM of the peaks was set to the same value as the spectral resolution of our observation camera used in this experiment (0.48 nm px−1).
4. Results and discussion
4.1. Compilation of the spectral features
Although the meteoric spectra show considerable variability, certain groups of commonly observed spectra can be identified. The spectral features in this work are presented in the form of groups of spectral lines that occur in narrow spectral regions between approximately 370 and 660 nm.
The purpose of this data compilation is to provide a detailed description of the selected spectral features in terms of assignment to the energy transitions that contribute significantly to the appearance of the most frequently observed spectral lines in the meteor spectra. As candidates for the most prominent meteoric spectral features, the lines listed in the literature (Vojáček et al. 2015; Ceplecha 1971, 1966; IMC 20151) were selected. For an overview of all the selected spectral features, see Fig. 2, where our laboratory low-resolution spectrum acquired with the observation camera is compared to the real spectra of several meteors (Ferus et al. 2018a).
Fig. 2. Panel A: high-resolution (gray) together with the low-resolution (red) spectrum. The most significant spectral lines are labeled. Panel B: examples of spectra recorded by observational spectographs. Starting with the upper spectrum of Perseid 20150811_014658 (brown), Sporadic 20161227_020734 (red), Perseid 20150812_231001 (black), Leonide 20161117_042009 (blue), Ursae Majorid 20161213_014212 (green), Alpha Camelopardalid 20161002_013415 (purple), and Perseid 20150812_232102 (wine). Overview of the selected spectral features (marked with gray lines). The detailed description of the individual spectral lines is provided in Table A.1. |
Because this study uses extensively the experimental data measured by the TC-LIBS of real meteoritic samples, we mainly focussed on the spectral lines assigned to the metalic species. The spectral lines of atmospheric gases that are also observed in meteoric spectra are included, but they are not paid much attention.
All the spectral features are listed in the Table A.1. Each feature is composed of one or more spectral lines that can include several energy transitions. The features are separated by horizontal lines that reach over all the table columns. Each feature is labeled by a tag (F1, F2, etc.) that serves as its reference in the text and figures. Each feature in the table is accompanied by a figure (Figs. B.1–B.4). All the marked peaks (their positions and heights) are listed in the corresponding part the Table A.1. In addition to the experimental and simulated spectra, the individual energy transitions are plotted (as a stick diagram). The columns of Table A.1 can be divided into two main groups. The first group contains the data of the energy transitions, and the second group collects the information about the selected spectral lines (peaks). Each spectral line can be composed of several energy transitions. The table lines (transitions) contributing to a specific peak are separated by one column spanning horizontal lines in the last four columns of the table (peak wavelengths and intensities). Each peak begins at the row of the table where the values of its wavelength and height are written and ends at the nearest following horizontal line. The beginning and end of each peak were determined by the nearest left and nearest right local minima, respectively. Lines that can be observed too far from the listed transitions are not included in the table, nor are they marked in the figures.
Only transitions with known Einstein A-coefficients were added into the table. Moreover to simplify the present table we included only the energy transitions gathered from the NIST Atomic Spectral Database Kramida et al. (2018) whose intensities (calculated using the Eq. (1)) are higher than 10% of the maximum calculated intensity within the particular spectral feature. Similarly, only the peaks exceeding the height of 10% of the highest peak in the feature are selected.
4.2. Description of spectral features
This section contains a brief description of the selected spectral features. Only interesting or doubtful cases are listed.
-
F1
(371.58–375.35 nm): almost only lines of Fe I.
-
F2
(381.64–389.03 nm): the low-temperature lines of Mg I and Fe I, which partially overlap. The lines of Si II can also become important at higher temperatures.
-
F3
(392.95–397.19 nm): typical for fast meteors. The most significant lines belong to the Ca II doublet. These lines are overlapped by two Al I lines of comparable intensity, however. This may cause difficulties when the lines are used for the quantitative analysis.
-
F4
(402.69–406.74 nm): the low-temperature lines of Fe I and Mn I.
-
F5
(412.68–413.50 nm): the fast meteor feature, the Si II line (higher temperatures), is very close to a line of Fe I, which is much more intensive at lower temperatures, however.
-
F6
(419.67–432.97 nm): low-temperature lines of Fe I, Ca I, and Cr I.
-
F7
(437.18–440.89 nm): low-temperature and weak lines of Fe I; at higher temperatures, lines of Mg II and Fe II can apear.
-
F8
(441.65–448.44 nm): fast meteors, contains the low-temperature lines of Ca I and Fe I, at higher temperatures, lines of Mg II can appear. The measured high-resolution spectrum also contains one significant peak that can be assigned to Mn I energy transitions.
-
F9
(451.61–454.40 nm): low-temperature lines of Ca I and Fe I. At high temperature, the lines of Fe II become more significant. The measured high-resolution spectrum also contains one significant peak that was not selected and included in the table because the Einstein A-coefficient of the nearby Fe I transition is missing.
-
F10
(469.90–470.70 nm): low-temperature lines of Mg I.
-
F11
(486.80–496.16 nm): low-temperature lines of Fe I. At high temperature, the line of Fe II becomes more significant.
-
F12
(500.74–502.11 nm): fast meteors, Fe I and Fe II lines, together with transitions of Ti I, can also contribute to the lines observed in this spectral region.
-
F13
(510.59–511.49 nm) contains Fe I and Fe II transitions.
-
F14
(516.33–522.41 nm): low-temperature and weak lines of Mg I, Fe I, and Cr I.
-
F15
(526.49–527.34 nm): low-temperature Fe I lines, lines of Ca I are also of comparable intensity in this spectral region.
-
F16
(532.60–533.60 nm): lines of atmospheric O I (at higher temperatures) and Fe I (lower temperatures).
-
F17
(536.75–545.96 nm): in this spectral region and at low temperature Fe I, Cr I lines are most significant. At high temperature, the lines of Fe II become more significant.
-
F18
(552.10–553.40 nm): the dominant Mg I line.
-
F19
(557.29–561.94 nm): the dominant lines in this spectral region belong to Ca I and Fe I. At high temperature, the lines of Fe II can appear.
-
F20
(588.60–589.98 nm): low-temperature Na I lines (usually a very intensive feature).
-
F21
(623.10–624.42 nm): low-temperature lines of Fe I. At high temperature, the lines of Fe II and Si II can appear.
-
F22
(634.32–637.52 nm): fast meteors, high-temperature lines of Si II.
-
F23
(645.24–649.88 nm): lines of atmospheric O I and N I, but the lines of Ca I and Fe I are also present in this region.
-
F24
(655.88–656.69 nm): fast meteors, H I lines. At the experimental spectrum only the left wing of the spectral line can be seen.
-
F25
(674,55–679,41 nm): low-temperature lines of Fe I and Ni I. At higher temperature, the lines of Mg II and N I can appear.
-
F26
(693.72–698.59 nm): in this spectral region and at low temperature Fe I, Cr I lines are most significant. At high temperature, the lines of Fe II and N I become more significant.
4.3. Meteor spectra in the laboratory: advantages, challenges, and limitations
Evaluating meteor emission spectra is a very complex scientific problem. Our experimental spectra depicted in panel A of Fig. 2, as well as several examples of spectra recorded by observational spectrographs in panel B, exhibit a series of overlapping bands that are typical for complex multicomponent matrices. Our experimental data as well as their comparison with synthetic spectra calculated using the NIST database show that position, range, and intensity of many spectral features can depend not only on the chemical composition of the meteoroid body, but also on their overlapping and on the temperature of the plasma. Meteor spectra recorded by the low-resolution observational camera were obtained by the analysis of a spectrum photography. For this analysis, theoretical tables cannot serve as calibration spectroscopic standard. Experimental LIBS of chondrites by terawatt-class laser and their assignment introduced in this study can therefore help with the precise assignment of spectral lines. In Table A.1 we show that wavelenghts assigned to features observed in meteor spectra can differ. Experimental data can provide a standard for the spectral assignment for such complicated and complex matrices. Experimental data also show that considering intensity and number of spectral features in UV region between 115–400 nm, accuracy and sensitivity of meteor spectroscopy could be improved by measuring in this range. However, this is beyond the capabilities of the current instrumentation; detection outside Earth’s atmosphere using very high-resolution spectrographs are required. Standard observational cameras operate in the 400 nm to 780 nm region, and emission lines above 300 nm are significantly attenuated by absorption in atmospheres as well as by the decreasing sensitivity of detectors. On the other hand in regions of longer wavelengths we performed in our laboratory systematic studies of the Rydberg states of atoms in the near- and mid-infrared ranges (NIR and MIR) (Civiš et al. 2011, 2012a,b,c,d, 2013a,b) (e.g., Fourier transform type) (Kawaguchi et al. 2008; Civiš et al. 2012e,c,f). Molecules and atomic species exhibit strong emission spectra under these conditions. In this manner, the elemental and molecular compositions of the objects can be simultaneously studied in the NIR and MIR regions. However, there are several limitations: the emission spectra of atoms are either not known or they are not assigned to Rydberg states that emit infrared radiation.
We supply a comparison of averaged spectra recorded by TC-LIBS of chondrite specimens with a selection of several spectra of meteors recorded at the Observatory Valašské Meziříčí (Czech Republic, Moravia - Zlín Region, GPS 49.46373, 17.97366). Our TC-LIBS measurements exhibit a wide range of emission lines. The most intense include Fe I, II, Mg I, II, Ca I, II, Mn I, Si II, Cr I, II, Ti I, Al I, Na I, H I, and N I, II (as described in Sect. 4.2, Figs. 2 and B.1–B.5, and Table A.1). However, we should note that several studies reported a number of rather extraordinary spectra: For instance, our experiments always exhibit intense sodium lines. Meteors with expected chondrite composition exhibit either a lack of sodium (Borovička et al. 2005) or a high concentration of sodium is described (Trigo-Rodriguez & Llorca 2007). The effect is explained by the depletion in Na through the exposure to solar radiation. Other meteors are classified based on the spectra by iron (in agreement with common meteorite classification). They are classified into mainstream and normal (with expected composition close to chondritic), and Fe poor (Vojáček et al. 2015). The TC-LIBS experiments will also facilitate a more precise classification, when the results are also connected with an abundance determination of particular elements. In these investigations, the simulation of relatively optically thin plasma around the path of the meteor is another challenging experimental problem.
Laboratory studies focused on meteor plasma are complicated by the fundamental problem connected with every experimental work: achieving experimental conditions as close to the real system as possible. The large volume of plasma produced by a terawatt-class laser exhibits a temperature of 9600 K and expands to 6 km s−1 (based on the Doppler shift of the spectral lines). Furthermore, regarding the high-energy output of 600 J of a terawatt-laser, this device is powerful enough to ablate significantly large areas of about 1 cm2 of meteorite. This is similarly to the ablation of small meteoroid bodies in the atmosphere. This is not possible with common laboratory lasers. These conditions can serve as a kind of simulation of the frontal shock wave region of the meteor.
This frontal shockwave contributes to the high-temperature component in meteor emission spectra. Two meteor spectral components were described by Borovička et al. for the first time in 1993 (Borovička 1994; Borovička & Betlem 1997; Borovička et al. 1999; Berezhnoy & Borovička 2010). The intense low-temperature component exhibits an excitation temperature ranging from 2900 K (lower extreme) to 6200 K (higher extreme) with a typical value of about 4500 K. Jenniskens et al. (2004) discovered a slight increase in the temperature with decreasing meteor altitude, but otherwise nearly constant values for meteoroids with speeds between 35 and 72 km s−1 and masses between 1025 g and 1 g. However, Ceplecha (1965) observed a decrease in the excitation temperature as a function of meteor luminosity and also reported a temperature of atomic oxygen lines in airglow plasma of about 14 000 K. Typical lines of the second component make up only 0.02% of the meteor vapor envelope in slow meteors, but account for more than 5% in fast meteors (Borovička 1994). These lines consist of high excitation of H I, Mg II, Si II, or Ca II that is also observed in the spectra of TC-LIBS recorded in this work. However, they are also mixed with a series of neutral lines in our case. We should note that global spectroscopic observations still oversimplify the much more complicated structure as well as temperature distribution in the plasma fireball that surrounds the meteoroid. The problem is described in detail in a recent comprehensive review by Silber et al. (2018). Our pilot study of TC-LIBS again opens this topic. TC-LIBS will certainly offer a unique opportunity to study the complicated phenomenon of meteoroid surface ablation under strictly controlled laboratory conditions on the macroscopic scale of large laser sparks.
5. Conclusions
We presented a compilation of the dominant spectral features of the plasma produced by ablation of chondrite meteoroids. Laboratory simulations of meteor plasma were performed by ablating of several specimens of chondrite meteorites by the high-power terawatt-class laser facility PALS. They reached the energy of 600 J and an output intensity of 2.2 TWcm−2 when focused on an area of 1 cm2 of the meteorite specimen. The produced ablation plasma attained the temperature of 9600 K and plasma expanded with a velocity of 6 km s−1. Spectra were recorded by the high-resolution echelle spectrograph, and for comparison and calibration, spectra were also recorded by the astronomical meteor spectrograph. All the spectral features were interpreted by synthetic spectra that were calculated based on the data from the NIST atomic spectra database. In this manner, we provide a data set for the evaluation of observational spectra that are recorded by the meteor spectrographs. Experimental spectra were also compared with selected meteor spectra.
Acknowledgments
This paper has been published as part of research series supported by the Czech Science Foundation within the project reg. no. 18-27653S, Programme of Regional Cooperation between the Regions and the Institutes of the Czech Academy of Sciences in 2018 (projects no.: R200401721 and R200401801), ERDF/ESF “Centre of Advanced Applied Sciences” (No. CZ.02.1.01/0.0/0.0/16_019/0000778) and Ministry of Education, Youth and Sports – PALS RI (LM2015083). Moreover, we would like thank to PALS crew: J. Skála, J. Mare, P. Prchal, M. erveák, J. Golasowski, J. Hebíek for their help and support. Contribution Statement. MF invetended the research, conducted the experiments, evaluated the data, and wrote the paper, PK conducted the experiments, conducted the spectra simulation, and wrote the paper, LP, LL, JK, AK, AP, JH, TK, SC participated in the experiments, AK evaluated the data and wrote the paper, HS wrote the paper, LJ, RD contributed from the side of the PALS facility, EC conducted EDS and wrote the paper, MK supervised and conducted experiments in the PALS facility and wrote the paper.
References
- Aristova, E. Y., Aushev, A. A., Baranov, V. K., et al. 2018, J. Exp. Theor. Phys., 126, 132 [NASA ADS] [CrossRef] [Google Scholar]
- Babánková, D., Civiš, S., & Juha, L. 2006a, Prog. Quantum Electron., 30, 75 [NASA ADS] [CrossRef] [Google Scholar]
- Babánková, D., Civiš, S., Juha, L., et al. 2006b, J. Phys. Chem. A, 110, 12113 [CrossRef] [Google Scholar]
- Berezhnoy, A. A., & Borovička, J. 2010, Icarus, 210, 150 [NASA ADS] [CrossRef] [Google Scholar]
- Borovička, J. 1994, Planet. Space Sci., 42, 145 [Google Scholar]
- Borovička, J., & Betlem, H. 1997, Planet. Space Sci., 45, 563 [NASA ADS] [CrossRef] [Google Scholar]
- Borovička, J., Stork, R., & Bocek, J. 1999, Meteorit. Planet. Sci., 34, 987 [NASA ADS] [CrossRef] [Google Scholar]
- Borovička, J., Koten, P., Spurný, P., Boček, J., & Štork, R. 2005, Icarus, 174, 15 [NASA ADS] [CrossRef] [Google Scholar]
- Ceplecha, Z. 1965, Bull. Astron. Inst. Czechosl., 16, 88 [NASA ADS] [Google Scholar]
- Ceplecha, Z. 1966, Bull. Astron. Inst. Czechosl., 17, 195 [NASA ADS] [Google Scholar]
- Ceplecha, Z. 1971, Bull. Astron. Inst. Czechosl., 22, 219 [NASA ADS] [Google Scholar]
- Ciucci, A., Corsi, M., Palleschi, V., et al. 1999, Appl. Spectros., 53, 960 [NASA ADS] [CrossRef] [Google Scholar]
- Civiš, S., Juha, L., Babánková, D., et al. 2004, Chem. Phys. Lett., 386, 169 [NASA ADS] [CrossRef] [Google Scholar]
- Civiš, S., Babánková, D., Cihelka, J., Sazama, P., & Juha, L. 2008, J. Phys. Chem. A, 112, 7162 [Google Scholar]
- Civiš, S., Matulková, I., Cihelka, J., et al. 2011, J. Phys. B: At. Mol. Opt. Phys., 44, 025002 [NASA ADS] [CrossRef] [Google Scholar]
- Civiš, S., Ferus, M., Kubelík, P., Chernov, V. E., & Zanozina, E. M. 2012a, J. Phys. B: At. Mol. Opt. Phys., 45, 175002 [NASA ADS] [CrossRef] [Google Scholar]
- Civiš, S., Ferus, M., Kubelík, P., Chernov, V. E., & Zanozina, E. M. 2012b, A&A, 545, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Civiš, S., Ferus, M., Kubelík, P., Jelinek, P., & Chernov, V. E. 2012c, A&A, 541, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Civiš, S., Ferus, M., Kubelík, P., et al. 2012d, A&A, 542, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Civiš, S., Ferus, M., Kubelík, P., et al. 2012e, J. Opt. Soc. Amer. B, 29, 1112 [NASA ADS] [CrossRef] [Google Scholar]
- Civiš, S., Kubelík, P., & Ferus, M. 2012f, J. Phys. Chem. A, 116, 3137 [CrossRef] [Google Scholar]
- Civiš, S., Ferus, M., Chernov, V., Zanozina, E., & Juha, L. 2013a, J. Quant. Spectrosc. Radiat. Transfer, 129, 324 [NASA ADS] [CrossRef] [Google Scholar]
- Civiš, S., Ferus, M., Chernov, V. E., & Zanozina, E. M. 2013b, A&A, 554, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Civiš, S., Szabla, R., Szyja, B. M., et al. 2016a, Sci. Rep., 6 [Google Scholar]
- Civiš, M., Ferus, M., Knížek, A., et al. 2016b, Phys. Chem. Chem. Phys., 18, 27317 [CrossRef] [Google Scholar]
- Civiš, S., Knížek, A., Ivanek, O., et al. 2017, Nat. Astron., 1, 721 [NASA ADS] [CrossRef] [Google Scholar]
- Dell’Aglio, M., Giacomo, A. D., Gaudiuso, R., Pascale, O. D., & Longo, S. 2014, At. Spectrosc., 101, 68 [NASA ADS] [CrossRef] [Google Scholar]
- Drouard, A., Vernazza, P., Loehle, S., et al. 2018, A&A, 613, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ebert, M., Hecht, L., Hamann, C., & Luther, R. 2017, Meteorit. Planet. Sci., 52, 1475 [NASA ADS] [CrossRef] [Google Scholar]
- Ferus, M., Matulková, I., Juha, L., & Civiš, S. 2009, Chem. Phys. Lett., 472, 14 [NASA ADS] [CrossRef] [Google Scholar]
- Ferus, M., Kubelík, P., & Civiš, S. 2011, J. Phys. Chem. A, 115, 12132 [CrossRef] [Google Scholar]
- Ferus, M., Civiš, S., Mládek, A., et al. 2012, J. Am. Chem. Soc., 134, 20788 [CrossRef] [Google Scholar]
- Ferus, M., Michalčíková, R., Shestivská, V., et al. 2014a, J. Phys. Chem. A, 118, 719 [CrossRef] [Google Scholar]
- Ferus, M., Nesvorný, D., Šponer, J., et al. 2014b, Proc. Nat. Acad. Sci., 112, 657 [Google Scholar]
- Ferus, M., Knížek, A., & Civiš, S. 2015, Proc. Nat. Acad. Sci., 112, 7109 [NASA ADS] [CrossRef] [Google Scholar]
- Ferus, M., Koukal, J., Lenza, L., et al. 2017a, 2017 19th International Conference on Transparent Optical Networks (ICTON) (IEEE) [Google Scholar]
- Ferus, M., Pietrucci, F., Saitta, A. M., et al. 2017b, Proc. Nat. Acad. Sci., 114, 4306 [CrossRef] [Google Scholar]
- Ferus, M., Kubelík, P., Knížek, A., et al. 2017c, Sci. Rep., 7 [Google Scholar]
- Ferus, M., Laitl, V., Knizek, A., et al. 2018a, A&A, 616, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ferus, M., Koukal, J., Lenza, L., et al. 2018b, A&A, 610, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Giacomo, A. D. 2011, Spectrochim. Acta Part B, 66, 661 [NASA ADS] [CrossRef] [Google Scholar]
- Hawkes, R. L., Milley, E. P., Ehrman, J. M., et al. 2007, Earth Moon Planets, 102, 331 [NASA ADS] [CrossRef] [Google Scholar]
- Horzňáčková, M., Plavčan, J., Rakovský, J., et al. 2014, Eur. Phys. J. Appl. Phys., 66, 10702 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jenniskens, P. 2007, Adv. Space Res., 39, 491 [NASA ADS] [CrossRef] [Google Scholar]
- Jenniskens, P., Laux, C. O., Wilson, M. A., & Schaller, E. L. 2004, Astrobiol., 4, 81 [NASA ADS] [CrossRef] [Google Scholar]
- Jungwirth, K., Cejnarova, A., Juha, L., et al. 2001, Phys. Plasmas, 8, 2495 [NASA ADS] [CrossRef] [Google Scholar]
- Kawaguchi, K., Sanechika, N., Nishimura, Y., et al. 2008, Chem. Phys. Lett., 463, 38 [NASA ADS] [CrossRef] [Google Scholar]
- Kramida, A., Ralchenko, Y., Reader, J., & NIST ASD Team (2018) 2018, NIST Atomic Spectra Database (version 5.5.6) [Google Scholar]
- Kurahashi, E., Yamanaka, C., Nakamura, K., & Sasaki, S. 2002, Earth Planets Space, 54, e5 [CrossRef] [Google Scholar]
- Loeffler, M., Baragiola, R., & Murayama, M. 2008, Icarus, 196, 285 [NASA ADS] [CrossRef] [Google Scholar]
- Madiedo, J. M., Ortiz, J. L., Trigo-Rodríguez, J. M., et al. 2014, Icarus, 231, 356 [NASA ADS] [CrossRef] [Google Scholar]
- Madiedo, J. M., Trigo-Rodríguez, J. M., Castro-Tirado, A. J., Ortiz, J. L., & Cabrera-Caño, J. 2013a, MNRAS, 436, 2818 [NASA ADS] [CrossRef] [Google Scholar]
- Madiedo, J. M., Trigo-Rodríguez, J. M., Lyytinen, E., et al. 2013b, MNRAS, 431, 1678 [NASA ADS] [CrossRef] [Google Scholar]
- Milley, E. P., Hawkes, R. L., & Ehrman, J. M. 2007, MNRAS, 382, L67 [NASA ADS] [CrossRef] [Google Scholar]
- Moroz, L., Fisenko, A., Semjonova, L., Pieters, C., & Korotaeva, N. 1996, Icarus, 122, 366 [NASA ADS] [CrossRef] [Google Scholar]
- Ozdín, D., Plavčan, J., Horzňáčková, M., et al. 2015, Meteorit. Planet. Sci., 50, 864 [NASA ADS] [CrossRef] [Google Scholar]
- Park, S.-Y., & Mazanek, D. D. 2003, J. Guidance Control Dyn., 26, 734 [NASA ADS] [CrossRef] [Google Scholar]
- Silber, E. A., Boslough, M., Hocking, W. K., Gritsevich, M., & Whitaker, R. W. 2018, Adv. Space Res., 62, 489 [NASA ADS] [CrossRef] [Google Scholar]
- Šponer, J. E., Szabla, R., Góra, R. W., et al. 2016, Phys. Chem. Chem. Phys., 18, 20047 [CrossRef] [Google Scholar]
- Takahashi, T., Thornton, B., Ohki, K., & Sakka, T. 2015, At. Spectrosc., 111, 8 [NASA ADS] [CrossRef] [Google Scholar]
- Tognoni, E., Cristoforetti, G., Legnaioli, S., & Palleschi, V. 2010, Spectrochim. Acta Part B, 65, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Trigo-Rodriguez, J., & Llorca, J. 2007, Adv. Space Res., 39, 517 [NASA ADS] [CrossRef] [Google Scholar]
- Trigo-Rodriguez, J., Llorca, J., & Fabregat, J. 2004a, MNRAS, 348, 802 [NASA ADS] [CrossRef] [Google Scholar]
- Trigo-Rodriguez, J. M., Llorca, J., Borovička, J., & Fabregat, J. 2004b, Earth Moon Planets, 95, 375 [Google Scholar]
- Vojáček, V., Borovička, J., Koten, P., Spurný, P., & Štork, R. 2015, A&A, 580, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zakharov, Y. 2003, IEEE Trans. Plasma Sci., 31, 1243 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Energy transitions and lines of the spectral features
Spectral features of the meteors.
Appendix B: Figures of the spectral features
Fig. B.1. Spectral features F1–F6. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
Fig. B.2. Spectral features F7–F12. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
Fig. B.3. Spectral features F13–F18. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
Fig. B.4. Spectral features F19–F24. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
Fig. B.5. Spectral features F25–F26. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
Appendix C: Simulation spectra of chondrite meteorites
Fig. C.1. Low-resolution simulation spectra of chondrite meteorites for gradually increasing temperatures. |
All Tables
Elemental composition of the meteorites used to calculate the synthetic spectra.
All Figures
Fig. 1. Panel A: vacuum interaction chamber with ablation set-up. A photography of a spectrum is shown in the embedded picture (left). Panel B: detailed experimental set-up with a collimator of Echelle spectrograph. Panel C: low resolution astronomical camera for comparative measurement. Panel D: sample of a meteorite prepared for the TC-LIBS experiment (left), high-power laser induced dielectric breakdown inside the interaction chamber (right). |
|
In the text |
Fig. 2. Panel A: high-resolution (gray) together with the low-resolution (red) spectrum. The most significant spectral lines are labeled. Panel B: examples of spectra recorded by observational spectographs. Starting with the upper spectrum of Perseid 20150811_014658 (brown), Sporadic 20161227_020734 (red), Perseid 20150812_231001 (black), Leonide 20161117_042009 (blue), Ursae Majorid 20161213_014212 (green), Alpha Camelopardalid 20161002_013415 (purple), and Perseid 20150812_232102 (wine). Overview of the selected spectral features (marked with gray lines). The detailed description of the individual spectral lines is provided in Table A.1. |
|
In the text |
Fig. B.1. Spectral features F1–F6. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
|
In the text |
Fig. B.2. Spectral features F7–F12. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
|
In the text |
Fig. B.3. Spectral features F13–F18. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
|
In the text |
Fig. B.4. Spectral features F19–F24. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
|
In the text |
Fig. B.5. Spectral features F25–F26. Upper part: synthetic spectra calculated under different temperatures: blue shows 4000 K, green shows 7000 K, and red shows 10 000 K. The intensity of each of these spectra is scaled to 1. Lower part: high-resolution emission spectrum of the meteorite laser ablation plasma. Each marked peak is listed in Table A.1 |
|
In the text |
Fig. C.1. Low-resolution simulation spectra of chondrite meteorites for gradually increasing temperatures. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.