2505.01100

# Meteoroid rotation and quasi-periodic brightness variation of meteor light curves **Authors**: Salvatore Mancuso, Dario Barghiniand Daniele Gardiol > Istituto Nazionale di Astrofisica, Osservatorio Astrofisico di Torino, via Osservatorio 20, Pino Torinese 10025, Italy ## Abstract Meteor light curves are sometimes known to display flickering: rapid, quasi-periodic variations in brightness. This effect is generally attributed to the rotational modulation of the ablation rate, which is caused by the time-varying cross section area presented by a nonspherical rotating meteoroid to the oncoming airflow. In this work we investigate the effects that the rotation of a meteoroid of given shape (spherical, cubic, or cylindrical) has on the meteor’s light curve, given state-of-the-art experimental laboratory estimates of the drag and lift coefficients of hypersonic flow (Mach number $>5$ ) around various shaped objects. The meteoroid’s shape is important in determining these two forces, due to the different response of the drag and lift coefficients according to the angle of attack. As a case study, the model was applied to a fireball observed on 2018 April 17 by the PRISMA network, a system of all-sky cameras that achieves a systematic monitoring of meteors and fireballs in the skies over the Italian territory. The results show that this methodology is potentially able to yield a powerful diagnostic of the rotation rate of meteoroids prior to their encounter with the atmosphere, while also providing essential information on their pre-fall actual shapes. meteorites, meteors, meteoroids / astrometry / methods: data analysis / techniques: image processing ## 1 Introduction A meteoroid passing through the Earth’s atmosphere, at a wide range of hypersonic velocities from about 11 km s -1 to above 73 km s -1 (Peña-Asensio et al. 2024), generates a visual fireball or bolide as a result of frictional heating induced by inelastic collisions with air molecules (Opik 1958; McKinley 1961; Ceplecha et al. 1998). Part of the meteoroid kinetic energy is then transformed into visible radiation, and some of its momentum is also transferred to the air; the resultant atmospheric drag decelerates the meteoroid (Hill et al. 2004; Rogers et al. 2005; Vinković 2007). When the ablated material interacts through collisions with the atmospheric or evaporated meteoric atoms and molecules, it excites them, such that we observe the luminosity emitted in spectral bands due to the decay of these excited states (Vondrak et al. 2008). Thermal ablation occurs intensively only after the meteoroid reaches the melting point, at the height of intensive vaporization (Plane 2012) occurring in a range of altitudes between about 80 and 125 km. While most of the falling meteoroids are too small to survive thermal ablation, the surviving ones are also subject to fragmentation during their descent trough the atmosphere before reaching the ground (Ceplecha et al. 1998; Fisher et al. 2000; Rogers et al. 2005). For the surviving meteorites, when ablation ceases, the body experiences a rapid exponential cooling with time, and the crust thus solidifies, impacting the ground at speeds of tens to hundreds of m s -1, according to their leftover mass. It is thus difficult to determine the actual physical characteristics of meteoroids entering the atmosphere since only a few of them are able to reach the ground, while small meteoroids ablate completely due to the removal of material from their surfaces. The visual brightness of the falling meteoroid can be recorded by ground-based cameras. Since the intensity of the light produced by a fireball is indicative of the almost instantaneous mass-loss rate, these light curves represent ideal indicators of the mode of meteor ablation and fragmentation, and even of the implied meteoroid structure. Classically, the vast majority of light curve shapes can be represented by a steady rise to a maximum brightness followed by a steep drop. However, faint meteors show a wide variety of light curve shapes (Fleming et al. 1993; Murray et al. 1999; Koten & Borovicka 2001; Beech & Murray 2003). A few of them also display flickering, a distinct quasi-periodic brightness modulation, whose amplitude typically changes with time and rarely exceeds 1 mag above or below the classical curve (Beech 2001). Flickering is present in about 5% of observed fireballs and is a distinct phenomenon from the effect of flaring, which is due to stochastic meteoroid fragmentation. Beech & Brown (2000) identified the rotation of nonspherical meteoroids into the Earth’s atmosphere as the most likely underlying physical mechanism leading to flickering. In this work we examine the effects that the rotation of a nonfragmenting meteoroid of spherical and nonspherical shapes has on the meteor’s light curve, given state-of-the-art laboratory estimates of the drag and lift coefficients of hypersonic flow (Mach number $>5$ ) around objects of various shapes (spherical, cubic, and cylindrical). We show that the meteoroid’s shape has an important impact on these two aerodynamic forces due to the different response of the drag and lift coefficients according to the body’s angle of attack. In particular, this analysis will yield an apt method to measure the rotation rate of meteoroids prior to their encounter with the atmosphere by examining the modulation of their brightness and speed dependence on height, while also providing essential information on the actual shape of meteoroids. This paper is organized as follows. In Section 2 we present the PRISMA network, a system of all-sky cameras that achieves a systematic monitoring of meteors in the skies over the Italian territory. Section 3 develops the time-dependent equations describing the meteoroid kinematics and ablation and discusses state-of-the-art laboratory estimates of the drag and lift coefficients of hypersonic flow around various shaped objects. In Section 4, as a case study, we apply the model to a specific fireball observed on 2018 April 15 by the PRISMA network. A discussion and our conclusions are presented in Section 5. ## 2 The PRISMA network The PRISMA network (Gardiol et al. 2016, 2019) is a system of tens of all-sky cameras that has been active since 2016; its aim is to achieve a systematic monitoring of meteors and fireballs in the skies over the Italian territory. PRISMA is a part of the international collaboration initiated by the Fireball Recovery and InterPlanetary Observation Network (FRIPON) project (Colas et al. 2014; Jeanne et al. 2019; Colas et al. 2020). As a partner of the FRIPON collaboration, PRISMA shares the same technology of the network. Each station is equipped with a CCD camera (6 mm diagonal, $1296× 966$ pixel) coupled with a short focal lens objective (1.25 mm) in order to obtain an all-sky field of view. The camera, connected to a Linux operating mini-PC via LAN and controlled by the open-source FREETURE software (Audureau et al. 2014), is operated at 30 frames per second (1/30 s exposure time) in order to sample the meteor trail at a suitable rate. The meteor detection is triggered locally in each node by a frame difference method, and cross-correlated with respect to the data of other nodes to check for multiple detections of the same event. The video stream of the detected meteor is saved locally as FITS files and, in the case of a multiple detection, the stream is collected by the FRIPON central server located at the Laboratoire d’Astrophysique de Marseille (LAM). Station monitoring, network security, software maintenance, real-time data processing, and sharing tasks are in charge of FRIPON and LAM teams. PRISMA data are synchronized and stored at the Italian Center for Astronomical Archives (IA2) INAF archiving facilities in Trieste (Knapic et al. 2014). The PRISMA reduction pipeline is implemented in IDL8 v8.7 and MATLAB9 Release 2015b and it was developed at the Osservatorio Astrofisico di Torino and Osservatorio di Astrofisica e Scienze dello Spazio in Bologna, which are part of the Italian National Institute for Astrophysics (INAF). The astrometric calibration of the all-sky cameras of the PRISMA fireball network has been thoroughly described in Barghini et al. (2019). To date, PRISMA has observed hundreds of meteors, and even helped to locate the fragments of two meteoroids in Italy near Cavezzo, on 2020 January 4 (Gardiol et al. 2021; Boaca et al. 2022; Bizzarri et al. 2023), and near Matera, on 2023 February 14. ## 3 Model The dynamics and ablation of a self-similar nonfragmenting meteoric body of mass $m$ traveling through the Earth’s atmosphere at a speed $v$ can be modeled by means of two kinematic differential equations that describe the change in velocity (deceleration or drag) and mass loss (ablation) of a meteoroid traveling at hypersonic speed (e.g., Bronshten 1983): $$ \displaystyle\frac{dv}{dt}=-\frac{1}{2m}C_DSρ_av^2+g\sin{θ}, \tag{1} $$ $$ \displaystyle\frac{dm}{dt}=-\frac{1}{2}\frac{C_H}{Q}Sρ_av^3. \tag{2} $$ In the above equations, $ρ_a$ is the density of the atmosphere through which the particle moves, $v$ is the magnitude of the velocity vector $v$ , $θ$ is the inclination angle between $v$ and the perpendicular to the gravitational acceleration $g$ , and $S$ is the middle cross sectional area of the frontal bow shock presented to the oncoming airflow and projected in the direction perpendicular to the flow. In eq. (1) the coefficient of drag $C_D$ describes the contribution on falling meteoroids of drag forces $F_D=\frac{1}{2}C_Dρ_av^2S$ opposite to the motion. This aerodynamic coefficients critically depends on the shape of the meteorite, but also on its speed and size and the properties of the atmospheric environment (and thus the flow regime) through the Reynolds ( $Re$ ) and Mach ( $Ma$ ) numbers (e.g., Charters & Thomas 1945; Bailey & Hiatt 1972; Henderson 1976; Spearman 1991). A thorough discussion on flow regimes appropriate for meteor studies can be found in Moreno-Ibáñez et al. (2018). In eq. (2) $C_H$ is a dimensionless heat transfer coefficient, which measures the efficiency of the collision process in converting kinetic energy into heat, while $Q$ is the heat of ablation or the energy required to ablate a unit mass of the meteoroid, which is a function of the material type and the specific process of ablation. According to Bronshten (1983), based on the best fit to the observational data, this value is about $C_H≈ 0.1$ or even lower by about an order of magnitude (e.g., Gritsevich & Koschny 2011), although this coefficient actually depends on speed $v$ , altitude $h$ , atmospheric density $ρ_a$ , and body size. On the other hand, a significant uncertainty is associated with this value (e.g., Johnston et al. 2018), and in many works $C_H$ is not even assumed, being identified as an integral part of a mass-loss parameter $β$ (see Gritsevich 2007, 2009). The value of $Q$ is usually taken to be constant and equal to the heat of vaporization; a typical value for an ordinary chondrite meteoroid of density $ρ_m=3.5· 10^3$ kg/m 3 is $Q≈ 8· 10^6$ J/kg (Chyba et al. 1993). The ratio of the two parameters $C_H$ and $Q$ corresponds to the ablation parameter, defined as $C_ab≡ C_H/Q$ (e.g., McMullan & Collins 2019; Brykina & Bragin 2020; Bragin & Brykina 2021). We note that while $C_ab$ is clearly not fixed with velocity and radius, since this parameter is only poorly constrained in the literature, we follow previous works and assume that it remains constant during the meteoroid’s fall, at least in the range of heights where the fireball was observed. In our model, $C_ab$ is thus left as a free parameter. The assumption of negligible lift for meteoroids is usually reasonable for first-order evaluations of aerodynamic effects on quasi-spherical falling bodies. However, if the shape of the meteoroid is nonspherical or elongate, a lift force $F_L=\frac{1}{2}C_Lρ_av^2S$ is expected when the object moves at an angle of attack $α$ (calculated with respect to the semimajor axis of the elongated object) with respect to the resultant velocity vector $v$ of the body. The coefficients of lift $C_L$ also depends on the shape of the meteorite and on the Reynolds ( $Re$ ) and Mach ( $Ma$ ) numbers. A maximum value for $C_L$ of order $∼ 0.001$ has been estimated by Passey & Melosh (1980) on the basis of the observed final distribution of impact craters, but this upper limit strictly applies only to massive ( $>100$ kg) meteoroids, while it is expected to be more pronounced for smaller meteors with shallow entry angles. For nonspherical bodies, the above system of differential equations is thus necessarily completed by adding two other differential equations that describe the variations in the inclination angle $θ$ with time along a meteoroid’s trajectory and the variation in a meteoroid’s height above the Earth’s surface $h$ of radius $R_E$ as a function of time: $$ \displaystyle\frac{dθ}{dt}=≤ft(\frac{g}{v}-\frac{v}{h+R_E}\right)\cos {θ}-\frac{1}{2m}C_LSρ_av, \tag{3} $$ $$ \displaystyle\frac{dh}{dt}=-v\sin(θ). \tag{4} $$ While the gravitational force $F_G=mg$ is always directed towards the Earth’s center, the aerodynamic drag $F_D$ and lift $F_L$ are respectively anti-parallel and perpendicular to the velocity vector. The above system of four differential equations can be solved numerically by means of a fourth-order Runge–Kutta scheme if we provide an apt model for the atmospheric density $ρ_a$ as a function of altitude $h$ and appropriate boundary conditions. Here we assume an isothermal atmosphere $ρ_a(t)=ρ_o\exp(-h(t)/H)$ , where $ρ_o≈ 1.29$ kg m -3 is the atmospheric density at sea level and $H≈ 7.16$ km is the constant scale height. We note however that more reliable results can be obtained by using real atmospheric measurements provided by national weather services for the times, heights, and locations of the specific fireballs under study (e.g., Lyytinen & Gritsevich 2016). A different approach to the solution of the meteor’s dynamics was proposed by Gritsevich & Koschny (2011), based on the work by Stulov et al. (1995). This method has the advantage that is entirely based on the interpretation of variations in speed and height as consequences of braking and mass loss, with no assumption on meteoroid density or shape. In this approach the basic equations of meteor physics are rewritten using dimensionless parameters, and it is assumed that the mass $m$ and cross section $S$ are connected through $S/S_e=(m/m_e)^μ$ , where $μ$ is a coefficient ranging from 0 to 2/3, representing the effect of the body’s change in shape, and the subscripts $e$ indicate the parameters at the entry to the atmosphere. This is related to its rotation, which prevents shape change and may distribute heat all around the surface: a value of $μ=0$ corresponds to the case where the maximum heating and evaporation occurs in the front of the meteoroid, while $μ=2/3$ represents a uniform mass loss of the meteoroid over the entire surface. In this different approach, the shape change coefficient, $μ$ , may be estimated through light curve analysis by fitting the meteor’s luminosity. In synthesis, this method allows us to infer the initial mass, luminous efficiency, and ablation coefficient from a least-squares adjustment of the luminosity and velocity observations to the solution of the equations of atmospheric entry. Further details of the method can be found in Gritsevich & Koschny (2011). The amount of light produced by the meteor per second in the instrument passband, or meteor luminosity $I$ , is assumed to be proportional to the rate of change in the meteoroid’s kinetic energy ${E_kin}$ and is described by the luminosity equation: $$ \displaystyle I=-τ\frac{E_kin}{dt}=-τ≤ft(\frac{v^2}{2}\frac{dm}{dt }+mv\frac{dv}{dt}\right)\hskip 14.22636pt[\rm J\hskip 5.69046pts^-1]. \tag{5} $$ This expression takes into account both the the mass-loss rate of ablation $\frac{dm}{dt}$ and the contribution from the deceleration $\frac{dv}{dt}$ . The coefficient $τ$ in eq. (5) is the (dimensionless) instantaneous luminous efficiency, a parameter that yields the fraction of the meteoroid’s kinetic energy that is converted into visible radiation. The visual meteor luminosity, the portion of the total radiation in the optical bandpass (400–700 nm), can be finally expressed in terms of absolute magnitude $M_v$ given by $$ M_v=-2.5(\log_10I-C), \tag{6} $$ where $C$ is a calibration constant, here taken as $C=3.185$ (e.g., Ceplecha & ReVelle 2005; Gritsevich & Koschny 2011; Drolshagen et al. 2021). Since the typical speed in the atmosphere along the meteor’s trajectory is much higher than 11 km s -1, the flow regime is always expected to be hypersonic ( $Ma=v/c_s>5$ ) with a Reynolds number $Re$ above $∼$ $10^4$ . Depending on the speed of the object, its size, and mass-loss rate, a bow shock wave may occur ahead of the meteoroid’s trajectory. The front bow shock separates a uniform flow from the field of perturbed flow in the shock layer around the point of braking or the stagnation point where the flow becomes subsonic. In the hypersonic case, the values of the drag and lift coefficients are relatively simple to estimate and do not critically depend on the atmospheric conditions. In the literature the drag coefficient $C_D$ is usually set as a constant (between about 0.5 and 2), although it actually depends on many parameters, such as the shape of the meteoroid, its velocity, and the air density. Since we model the meteoroid entering the atmosphere as a rotating body subject to both drag and lift effects, its shape is a major factor in determining these two forces, so that its cross sectional area $S$ has to be time dependent, thus somehow influencing the motion, also due to the different response of the drag and lift coefficients according to the angle of attack. In recent decades, the drag coefficient has been calculated for several restricted ranges of $Re$ and $Ma$ , and its estimation has often relied on ad hoc interpolations between different regimes. Spherical bodies are adequately studied aerodynamic test objects (Charters & Thomas 1945; Hodges 1957; Sighard 1965; Bailey & Hiatt 1972; Henderson 1976; Spearman 1991). These studies experimentally ascertained that the drag coefficient $C_D$ of a spherical body does not depend on the Mach number $Ma$ when the flow is hypersonic and has an approximately constant value. In general, for spheres traveling at hypersonic speeds in the continuum regime, the value of the drag coefficient is thus independent from $Re$ and $Ma$ , and assumes an approximate value of $C_D≈ 0.92$ (Sighard 1965). <details> <summary>extracted/6406308/Fig1.png Details</summary>  ### Visual Description ## Aerodynamic Coefficient Charts: Cube vs. Cylinder ### Overview The image contains two vertically stacked line charts comparing the aerodynamic coefficients of a cube and a cylinder as a function of the angle of attack. Both charts share the same x-axis and y-axis scales and legend structure, allowing for direct comparison. The data appears to be from a computational or experimental study of bluff body aerodynamics. ### Components/Axes * **Chart Type:** Two 2D line plots. * **X-Axis (Both Charts):** Label: `angle of attack [deg]`. Scale: Linear, from -30 to 390 degrees. Major tick marks are every 30 degrees (0, 30, 60, ..., 360). * **Y-Axis (Both Charts):** Unlabeled, but represents the magnitude of the coefficients. Scale: Linear, from 0.0 to 2.5. Major tick marks are every 0.5 units (0.0, 0.5, 1.0, 1.5, 2.0, 2.5). * **Legend (Both Charts):** Located in the top-left corner of each plot area. * **Top Chart (Cube):** * `C_D (cube)`: Solid blue line. Represents the drag coefficient for the cube. * `S`: Dashed red line. Likely represents a side force coefficient or a related parameter. * `C_L`: Dash-dot green line. Represents the lift coefficient. * **Bottom Chart (Cylinder):** * `C_D (cylinder)`: Solid blue line. Represents the drag coefficient for the cylinder. * `S`: Dashed red line. Same parameter as in the top chart. * `C_L`: Dash-dot green line. Same parameter as in the top chart. ### Detailed Analysis **Top Chart: Cube Aerodynamics** * **C_D (cube) - Blue Solid Line:** Exhibits a strong, periodic pattern with sharp peaks and broad minima. Peaks occur at approximately 0°, 90°, 180°, 270°, and 360°, with a value of ~1.7. Minima occur at approximately 45°, 135°, 225°, and 315°, with a value of ~1.2. The trend is perfectly symmetric about 180°. * **S - Red Dashed Line:** Shows a smoother, sinusoidal oscillation. Peaks at ~0.8 occur near 45°, 135°, 225°, and 315°. Minima at ~0.5 occur near 0°, 90°, 180°, 270°, and 360°. This pattern is 45° out of phase with the C_D peaks. * **C_L - Green Dash-Dot Line:** Also shows a sinusoidal oscillation, centered around 0. Peaks at ~0.2 occur near 60°, 150°, 240°, and 330°. Minima at ~-0.2 occur near 30°, 120°, 210°, and 300°. The pattern is symmetric about the x-axis. **Bottom Chart: Cylinder Aerodynamics** * **C_D (cylinder) - Blue Solid Line:** Shows a more complex periodic pattern. It features very sharp, high peaks at 0°, 180°, and 360°, reaching ~1.7. There are secondary, broader peaks at 90° and 270°, reaching ~1.2. Deep minima occur at approximately 45°, 135°, 225°, and 315°, with a value of ~0.8. The pattern is symmetric about 180°. * **S - Red Dashed Line:** Exhibits a smooth oscillation. Peaks at ~0.55 occur near 90° and 270°. Minima at ~0.2 occur near 0°, 180°, and 360°. The amplitude is smaller than for the cube. * **C_L - Green Dash-Dot Line:** Shows a complex oscillation. It has peaks at ~0.25 near 60°, 120°, 240°, and 300°. It has minima at ~-0.1 near 30°, 150°, 210°, and 330°. The pattern is symmetric about the x-axis but has a different shape compared to the cube's C_L. ### Key Observations 1. **Fundamental Periodicity:** Both shapes show aerodynamic coefficients that are periodic with the angle of attack, repeating every 90° for the cube and every 180° for the cylinder, reflecting their geometric symmetries. 2. **Drag Coefficient (C_D) Behavior:** The cube's drag is highest when a flat face is perpendicular to the flow (0°, 90°, etc.) and lowest when a corner is presented (45°, 135°, etc.). The cylinder's drag is highest at 0° and 180° (broadside) and has a secondary maximum at 90° and 270°. 3. **Phase Relationships:** For the cube, the peaks in `S` and the minima in `C_D` are aligned (e.g., at 45°). For the cylinder, the peaks in `S` align with the secondary peaks in `C_D` (at 90° and 270°). 4. **Lift Coefficient (C_L) Symmetry:** The C_L for both shapes oscillates symmetrically around zero, indicating no net lift over a full rotation, as expected for symmetric bodies. ### Interpretation These charts illustrate the fundamental differences in flow separation and pressure distribution between a sharp-edged bluff body (cube) and a smooth, curved bluff body (cylinder). * **Cube:** The sharp corners fix the flow separation points. The drag is maximized when the projected frontal area is largest (flat face to flow) and minimized when the flow can more smoothly navigate around a leading corner. The `S` coefficient likely represents a side force that peaks when the flow asymmetry is greatest (at 45° angles). * **Cylinder:** The smooth surface allows the separation point to move with the angle of attack, leading to a more complex drag curve. The very sharp drag peaks at 0° and 180° suggest a sudden change in the wake structure, possibly related to the onset of vortex shedding or a transition in the separation pattern. The secondary peaks at 90°/270° correspond to the orientation where the cylinder's curvature presents the most abrupt change to the oncoming flow. **Notable Anomaly:** The most striking feature is the extremely sharp, narrow peak in the cylinder's C_D at exactly 180°. This suggests a highly sensitive, possibly unstable, flow condition at that precise orientation, which could be an important design consideration for structures subject to wind loads. The data implies that small changes in angle near 180° for a cylinder result in large changes in drag force. </details> Figure 1: Drag and lift coefficients as a function of the angle of attack for the cases of a rotating cube (upper panel) and a rotating cylinder (lower panel), as derived by Seltner et al. (2019, 2021). The red curve in the middle of each panel shows, in arbitrary units, the varying cross sectional area. The drag coefficient for nonspherical objects has also been a field of study for several decades. Theoretical models for the drag force on objects have however been derived only for a limited range of conditions, often subject to empirical limitations. Fewer studies have been carried out on cube-shaped bodies (Sighard 1965; Hansche & Rinehart 1952; Gerasimov et al. 2015). In particular, Hansche & Rinehart (1952) studied the dependence of the Mach number on the aerodynamic drag coefficient for rotating cubes free-flying in sub-, trans-, and supersonic flows, and a rough convergent value of $C_D≈ 1.66$ for Mach numbers $>4$ has been given by Sighard (1965). These earlier studies, however, dealt only with the drag, whereas the lift was neglected. Furthermore, the aerodynamic coefficients were not related to the orientation angles of the cubes. In addition, these works provide little information about the aerodynamic behavior of cubes in hypersonic flows. In the first half of the last century studies on the aerodynamic behavior over a wide range of Mach numbers of cylindrical objects with both blunt and streamlined ends were conducted (e.g., Vennard 1940; Rouse 1946; Gowen & Perkins 1952; Sighard 1965). These studies evidenced that the drag coefficient significantly depends on Mach number and is affected by the forebody shape in axial flows at transonic and supersonic speeds. However, these earlier studies tended to focus on the drag of cylindrical bodies in axial flows and crossflows at supersonic speed, while the issue of inclined right circular cylinders in hypersonic flows was neglected. Seltner et al. (2021) addressed this gap by publishing the results of experimental and numerical investigations of pitched cylindrical bodies in hypersonic flow, ascertaining both the aerodynamic drag and lift force. It has been recently proven by Seltner et al. (2019, 2021) that the aerodynamic drag and lift coefficients of rotating cubes and cylindrical bodies in hypersonic flow have a significant dependence on the body orientation. Seltner et al. (2019) calculate the aerodynamic coefficients of a rotating cube using the free-flight technique in a hypersonic flowfield. A rotating cube of side $l$ yields a pitch-angle-dependent projected frontal area given by $S(θ)=l^2(|\sin(θ)|+|\cos(θ)|).$ The drag coefficient $C_D$ is thus not constant, but a periodic function of the pitch angle $θ$ (see Figure 1). The drag coefficient for a pitch-angle-dependent projected frontal area ranges between a minimum of about $C_D=1.16$ and a maximum of $C_D=1.71$ at a pitch angle equal to zero. Moreover, a periodic lift force also arises for nonzero angles, with a lift coefficient $C_L$ continuously ranging between values of about +0.2 and -0.2 with the characteristics of negative sine functions. The peak of the lift coefficient for a rotating cube is thus about 12% of the maximum drag coefficient. Seltner et al. (2021) also calculated the drag and lift coefficients for a rotating cylinder in hypersonic flow, obtaining a detailed dependence of $C_D$ and $C_L$ on the angle of attack, also shown in Figure 1. In particular, the drag coefficient of cylinders of length $L$ base area $S=π D^2/4$ , where $D$ is the cylinder’s diameter and $L=2D$ shows a more important influence of the inclination angle than in the case of cubes, with $C_D$ almost doubling between a pitch angle of 0 ∘ and 90 ∘. In the following section we apply the above-described model, complemented with the theoretical drag and lift coefficients calculated by Seltner et al. (2019) and Seltner et al. (2021), to a specific fireball observed by the PRISMA network on 2018 April 17. ## 4 A case study: The meteor observed on 2018 April 17 On 2018 April 17, at around 19:45 UT, three stations of the PRISMA network (ITER04-Bedonia, ITER05-Piacenza, ITVA01-Lignan) detected a brilliant fireball in the skies of northern Italy. The meteor was tracked from an initial height of 71.56 km to a final height of 44.45 km, with an initial velocity of about 17.720 km s -1 and a time of flight of about 1.97 s. The light curve (see middle panel of Figure 2), which reached a minimum absolute magnitude of -3.87, did not show sudden brightenings that might be attributed to important fragmentation episodes. We solved the system of differential equations by using a standard fourth-order Runge-Kutta algorithm (RK4), as supplied in the Interactive Data Language (IDL), using a time step of 2 ms. Runge-Kutta methods offer higher accuracy than the Euler integration method and are known for their accuracy, especially at higher orders. RK4, in particular, is widely used and allows a good balance between accuracy and computational efficiency. We applied the model described in the previous section to this particular event as a case study to explore the implications of the changing aerodynamic coefficients of a spherical or nonspherical rotating body as a function of the angle of attack during its fall through the Earth’s atmosphere. For the ablating meteoroid, we considered rotating spherical, cubic, and cylindrical shapes, as discussed in the previous section, by using the state-of-the-art theoretical aerodynamic coefficients for hypersonic flows given by Seltner et al. (2019, 2021); we neglected the possible changes in shape due to ablation. For modeling the absolute magnitude as a function of height, we had to choose a functional form for the luminous efficiency. The coefficient $τ$ is not a well-known quantity, and it critically depends on other parameters, such as $m$ , $v$ , and even $ρ_m$ . The uncertainty on $τ$ derives mostly from the fact that the majority of light emitted by fireballs is restricted to spectral lines that depend on the specific meteoroid’s chemical composition (Ceplecha et al. 1998; Campbell-Brown et al. 2012). On the other hand, Gritsevich & Koschny (2011) were able to compute $τ$ with only very few assumptions, starting from the drag and mass-loss equations. In their approach, considering the change in the meteoroid’s speed and mass during its trajectory and taking the geometrical relation along the path of the meteor into account, they solved the formulas for the meteoroid’s dynamical behavior. The results were then compared to the observed meteor’s drag rate and light curve, and $τ$ was finally computed on the basis of this comparison. Many values for $τ$ are quoted in the literature according to the different assumptions the authors made and the methods applied (e.g., Verniani 1965; Ceplecha & McCrosky 1976; Halliday et al. 1996; Hill et al. 2005; Weryk & Brown 2013). All these studies found a dependency of $τ$ on the speed $v$ of the meteoroid. However, according to a recent study of luminous efficiency values determined with modern high-resolution instruments and obtained by directly comparing dynamic and photometric meteoroid masses, there is no clear dependence of $τ$ on initial velocity (Subasinghe & Campbell-Brown 2018). Moreover, most nonfragmenting meteoroids appear to have luminous efficiencies of less than 1% (Subasinghe & Campbell-Brown 2018). Given the large uncertainty on the luminous efficiency $τ$ , this parameter is left as a free parameter in our model in order to match the observed absolute magnitude as a function of height. We thus solved the system of differential equations presented in the previous section by using the theoretical drag and lift coefficients given by Seltner et al. (2019, 2021) that varied as a function of time in dependence on the supposed rotation of the meteoroid along the trajectory and the angle of attack. The criteria for the best fit were a set of initial parameters (size of the body $h_0$ , $C_ab$ , and $τ$ ) that best fit in a least-squares sense both the fireball’s speed and the absolute magnitude dependence on height. The initial size of the body $h_0$ was identified with the radius in the case of the sphere, the side length in the case of the cube, and the height in the case of the cylinder. For the cubic and cylindrical bodies two more parameters, the period of rotation $P_0$ and a phase factor $φ_0$ , were also introduced in order to model the modulation of the light curve of the fireball and of, at least in principle, the speed during the meteoroid’s fall. The best-fit models for the three cases are displayed in Figure 2, together with an evaluation of the respective reduced chi-squared $χ_ν^2$ . Our best-fit model for the ideal spherical case, obtained by assuming a typical value for an ordinary chondrite meteoroid density ( $ρ_m=3.5· 10^3$ kg/m 3), implies an initial instantaneous luminous efficiency $τ=0.134$ %, an initial radius $h_0=0.0055$ m, and a constant ablation parameter $C_ab=1.74· 10^-8$ kg J -1. Given the assumed meteoroid density, the initial mass was $m=0.00247$ kg. As for the two nonspherical cases, we found $τ=0.041$ %, $h_0=0.0132$ m, $C_ab=2.45· 10^-8$ kg J -1 for the cubic case, and $τ=0.077$ %, $h_0=0.0184$ m, $C_ab=1.89· 10^-8$ kg J -1 for the cylindrical case. We also inferred a possible periodicity of $P_0=1.778$ s for the rotating cubic body and a possible periodicity of $P_0=0.964$ s for the rotating cylindrical body. Interestingly, the best-fitting model for the cylindrical case clearly displays a peculiar slight modulation of the speed curve, which is, at least in principle, in view of more refined future fireball observations, potentially discernible from the observed speed dependence with height. For the absolute magnitudes, by visual inspection, the curve corresponding to the rotating cylindrical shape appears to nicely follow the observed modulation although, perhaps counterintuitively, the minimum reduced chi-squared $χ_ν^2$ is actually achieved in the case of the spinning cubic shape, thus hinting to the possibility that an intermediate brick-shaped body rotating with a period of $P_0≈ 1.8$ s could have been responsible for the observed modulation. In the bottom panel of Figure 2, for a better visual evaluation of the effects of the varying aerodynamic coefficients on the modulation of the speed and absolute magnitude, we also show the height dependence of the theoretical drag and lift coefficients for the three selected cases. From the above case study it is thus clear that our model is able to satisfactorily explain the observed modulations in the light curves (and, at least potentially, in the speed curves) as being attributable to the nonsphericity of rotating falling meteoroids, even pointing to the possibility of determining, with some degree of confidence, their probable shapes and rotation periods. <details> <summary>extracted/6406308/Fig3.png Details</summary>  ### Visual Description \n ## Scientific Plot Analysis: Atmospheric Entry Simulation Data ### Overview The image contains three vertically stacked scientific plots comparing the behavior of three geometric shapes (sphere, cube, cylinder) during what appears to be atmospheric entry or descent. The plots show speed, absolute magnitude, and model coefficients as a function of decreasing altitude (h). All text is in English. ### Components/Axes **Common Elements:** - **X-axis (all plots):** Labeled "h [km]" (altitude in kilometers). The axis is inverted, with values decreasing from left (70 km) to right (45 km), indicating descent. - **Data Series:** Three shapes are consistently compared: - **Sphere:** Green line - **Cube:** Blue line - **Cylinder:** Red line - **Legends:** Positioned within each plot's frame (top-left for plots 1 & 2, bottom-left for plot 3). Each entry includes the shape name and a chi-squared (χ²) goodness-of-fit value in parentheses. **Plot 1 (Top): Speed vs. Altitude** - **Y-axis:** Labeled "speed [km s⁻¹]" (kilometers per second). Scale ranges from 12 to 18. - **Legend:** Located in the top-left quadrant. - sphere (χ² = 11.47) - cube (χ² = 11.65) - cylinder (χ² = 11.44) - **Data:** Black points with vertical error bars represent measured data. Colored lines are model fits. **Plot 2 (Middle): Absolute Magnitude vs. Altitude** - **Y-axis:** Labeled "Absolute magnitude". Scale ranges from -1.0 to -4.0 (note: more negative values indicate greater brightness). - **Legend:** Located in the bottom-left quadrant. - sphere (χ² = 7.19) - cube (χ² = 5.64) - cylinder (χ² = 7.42) - **Data:** Black points with vertical error bars represent measured data. Colored lines are model fits. **Plot 3 (Bottom): Model Coefficients vs. Altitude** - **Y-axis:** Unlabeled, but numerical scale ranges from -1.5 to 2.0. - **Legend:** Located in the bottom-left quadrant. It defines two line styles and the three shapes. - **Line Styles:** - Solid line: C₀ - Dash-dotted line: C₁ - **Shapes:** sphere (green), cube (blue), cylinder (red) - **Data:** Six lines are plotted (C₀ and C₁ for each shape). No error bars or data points are shown. ### Detailed Analysis **Plot 1: Speed vs. Altitude** - **Trend Verification:** All three model lines show a consistent downward slope from left to right, indicating speed decreases as altitude decreases (descent). - **Data Points:** The measured data (black points) cluster tightly around the model lines, especially between 65-50 km. Scatter increases slightly at the highest (~70 km) and lowest (~45 km) altitudes. - **Model Comparison:** The sphere (green) and cylinder (red) fits are nearly identical and track the central tendency of the data closely. The cube (blue) fit is slightly higher than the others between 60-50 km. The χ² values are very close (11.44 to 11.65), indicating similarly good fits for all three shapes. **Plot 2: Absolute Magnitude vs. Altitude** - **Trend Verification:** The data and models show a non-monotonic trend. Magnitude becomes more negative (brighter) from 70 km to a peak brightness around 55 km, then becomes less negative (dimmer) from 55 km to 45 km. - **Data Points:** The measured data shows significant scatter and larger error bars compared to the speed plot. There is a notable dip in brightness (less negative magnitude) around 60 km. - **Model Comparison:** The cube model (blue, χ²=5.64) provides the best fit to the data, particularly capturing the peak brightness near 55 km. The sphere (green) and cylinder (red) models underestimate the peak brightness. The cylinder model shows a sharp, unexplained dip near 52 km. **Plot 3: Model Coefficients (C₀ and C₁)** - **Component Isolation:** This plot is segmented into two distinct behaviors. - **C₀ (Solid Lines):** The sphere's C₀ is constant at 1.0. The cube's C₀ oscillates between ~1.0 and ~1.7 with a regular period. The cylinder's C₀ shows sharper, narrower peaks reaching ~1.8, also periodic. - **C₁ (Dash-dotted Lines):** All C₁ values are small (between -0.3 and 0.3). The sphere's C₁ is constant at 0.0. The cube's C₁ oscillates in anti-phase with its C₀. The cylinder's C₁ also oscillates with a smaller amplitude. - **Spatial Grounding:** The oscillations in C₀ and C₁ for the cube and cylinder are perfectly synchronized in altitude, suggesting they are coupled parameters in the model. ### Key Observations 1. **Shape-Dependent Fit Quality:** No single shape is the best fit for all measured parameters. The cube best fits the magnitude data, while sphere and cylinder are marginally better for speed. 2. **Brightness Peak:** The absolute magnitude data shows a clear peak in brightness (most negative value) at an altitude of approximately 55 km. 3. **Oscillatory Coefficients:** For non-spherical shapes (cube, cylinder), the model coefficients C₀ and C₁ are not constant but oscillate with altitude, indicating a periodic physical process (e.g., tumbling or changing cross-section). 4. **Data Scatter:** The magnitude measurements have significantly more noise/scatter than the speed measurements. ### Interpretation This data likely comes from a study simulating or observing the atmospheric entry of objects with different geometries. The plots collectively demonstrate that an object's shape critically influences its observable characteristics during descent. - **Physical Meaning:** The speed plot shows expected deceleration due to atmospheric drag. The magnitude plot suggests the object heats up and glows brightest around 55 km altitude, likely due to peak aerodynamic heating. The coefficient plot reveals that for non-spherical objects, the model parameters governing the interaction (perhaps related to drag or luminous efficiency) are not static but vary periodically, strongly implying the objects are tumbling or rotating, changing their presented cross-sectional area and orientation relative to the flow. - **Why It Matters:** The analysis shows that using a simple spherical model to interpret data from an irregularly shaped object (like an asteroid or spacecraft debris) could lead to errors, particularly in estimating brightness or interpreting light curves. The oscillatory coefficients provide a potential signature to infer shape or rotational state from remote observations. - **Anomaly:** The sharp dip in the cylinder's magnitude model near 52 km is an artifact not supported by the data points and may indicate a numerical instability or an oversimplification in the model for that specific geometry at that altitude. </details> Figure 2: Speed (top panel) and absolute magnitude (middle panel) vs height for the 2018 April 17 fireball, as seen from the PRISMA cameras. Best-fit model curves for the three chosen meteoroid shapes (sphere, cube, and cylinder) are also shown superposed to the data. The prominent oscillations in the light curve are interpreted as being due to the asymmetric shape of the spinning meteoroid. The bottom panel shows the height dependence of the theoretical drag and lift coefficients for the three selected shapes. ## 5 Discussion and conclusions The drag and lift coefficients of nonspherical bodies are of crucial importance for the correct computations of meteor-body parameters from observations. In this work the effects on the modulation of an observed meteor’s light curve due to the rotation of a rotating meteoroid of given shape that enters the Earth’s atmosphere was investigated by applying state-of-the-art experimental laboratory estimates of the drag and lift coefficients of hypersonic flow (Mach number $>5$ ) around various shaped objects (spherical, cubic, and cylindrical), as given by Seltner et al. (2019, 2021). The underlying major assumption is that the shape of the meteoroid does not change, and that during its descent through the atmosphere there are no episodes of fragmentation. We showed that the meteoroid’s shape has an important impact on the drag and lift forces due to the different response of the respective time-dependent aerodynamic coefficients according to the given angle of attack. These time-dependent coefficients are quite different from those generally used in the literature related to the studies of meteor events to calculate the motion of meteor bodies, which usually use an assumed guess value and/or a spherical shape. Our model was applied to specific light curve data associated with a fireball observed on 2018 April 17 by the PRISMA network, which showed a clear modulation that, excluding a priori the possibility of fragmenting episodes, could be attributed to the rotation of a meteoroid. Among the three cases, the cylindrical shape showed a peculiar speed modulation with height and appeared to match, at least visually, the observed absolute magnitude dependence with height of the meteor. The cubic case, however, was the best-fitting of the three shapes, and displayed the best result for the reduced chi-square, hinting to the possibility that an intermediate brick-shaped body rotating with a period of $P≈ 1.8$ s might have been responsible for the observed modulation with height of the absolute magnitude. As a final remark, we would like to point out that our enquiry, although it used accurate aerodynamic coefficients, was limited to specific perfectly shaped bodies, which may not represent the actual shapes of the vast majority of nonspherical meteoric bodies falling into Earth’s atmosphere. For example, the rough surface of a quasi-spherical meteoroid would make the value of the drag coefficient depart from its theoretical value of 0.92. Moreover, even when not accounting for the body’s surface roughness, if the meteor body is of a similar shape to those modeled in this work (cubic or cylindrical), the sharp edges of this nonspherical meteoroid moving through the atmosphere are expected to become rounded because of the effect of both melting and evaporation. This implies that the drag coefficient could considerably decrease with time due to the increasing rounding of the edges (e.g., Zhdan et al. 2007). Even with the above caveats, we demonstrated that our model yields a powerful means to estimate the rotation rate of meteoroids prior to their encounter with the atmosphere by examining the observed light and speed curves, while also providing essential information on their actual shapes. ## Acknowledgments We thank the referee for carefully reading our paper and providing very helpful and thoughtful comments on our manuscript. PRISMA is the Italian Network for Systematic surveillance of Meteors and Atmosphere. It is a collaboration initiated and coordinated by the Italian National Institute for Astrophysics (INAF) that counts members among research institutes, associations and schools. The complete list of PRISMA members is available at http://www.prisma.inaf.it. PRISMA was partially funded by Research and Education grants from Fondazione CRT (years 2016, 2019, 2020, 2022, 2023). PRISMA data are hosted by the INAF research e-infrastructure project IA2 (Italian Center for Astronomical Archives). The initial FRIPON hardware and software has been developed by the FRIPON-France core team under a French ANR grant (2014-2018). ## References - Audureau et al. (2014) Audureau, Y., Marmo, C., Bouley, S., et al. 2014, in Proceedings of the International Meteor Conference, Giron, France, 18-21 September 2014, ed. J. L. Rault & P. Roggemans, 39–41 - Bailey & Hiatt (1972) Bailey, A. & Hiatt, J. 1972, Aiaa Journal, 10, 1436 - Barghini et al. (2019) Barghini, D., Gardiol, D., Carbognani, A., & Mancuso, S. 2019, Astronomy & Astrophysics, 626, A105 - Beech (2001) Beech, M. 2001, MNRAS, 326, 937 - Beech & Brown (2000) Beech, M. & Brown, P. 2000, Planetary and Space Science, 48, 925 - Beech & Murray (2003) Beech, M. & Murray, I. S. 2003, Monthly Notices of the Royal Astronomical Society, 345, 696 - Bizzarri et al. (2023) Bizzarri, I., Barghini, D., Colombetti, P., et al. 2023, Applied Radiation and Isotopes, 194, 110651 - Boaca et al. (2022) Boaca, I., Gritsevich, M., Birlan, M., et al. 2022, ApJ, 936, 150 - Bragin & Brykina (2021) Bragin, M. & Brykina, I. 2021, Fluid Dynamics, 56, 566 - Bronshten (1983) Bronshten, V. A. 1983, Physics of Meteoric Phenomena - Brykina & Bragin (2020) Brykina, I. & Bragin, M. 2020, Planetary and Space Science, 187, 104942 - Campbell-Brown et al. (2012) Campbell-Brown, M. D., Kero, J., Szasz, C., Pellinen-Wannberg, A., & Weryk, R. J. 2012, Journal of Geophysical Research (Space Physics), 117, A09323 - Ceplecha et al. (1998) Ceplecha, Z., Borovička, J., Elford, W. G., et al. 1998, Space Sci. Rev., 84, 327 - Ceplecha & McCrosky (1976) Ceplecha, Z. & McCrosky, R. 1976, Journal of Geophysical Research, 81, 6257 - Ceplecha & ReVelle (2005) Ceplecha, Z. & ReVelle, D. O. 2005, Meteoritics & Planetary Science, 40, 35 - Charters & Thomas (1945) Charters, A. C. & Thomas, R. N. 1945, Journal of the Aeronautical sciences, 12, 468 - Chyba et al. (1993) Chyba, C. F., Thomas, P. J., & Zahnle, K. J. 1993, Nature, 361, 40 - Colas et al. (2020) Colas, F., Zanda, B., Bouley, S., et al. 2020, A&A, 644, A53 - Colas et al. (2014) Colas, F., Zanda, B., Bouley, S., et al. 2014, in Asteroids, Comets, Meteors 2014, ed. K. Muinonen, A. Penttilä, M. Granvik, A. Virkki, G. Fedorets, O. Wilkman, & T. Kohout, 110 - Drolshagen et al. (2021) Drolshagen, E., Ott, T., Koschny, D., et al. 2021, A&A, 650, A159 - Fisher et al. (2000) Fisher, A., Hawkes, R., Murray, I., Campbell, M., & LeBlanc, A. 2000, Planetary and Space Science, 48, 911 - Fleming et al. (1993) Fleming, D. E., Hawkes, R. L., & Jones, J. 1993, in Meteoroids and their parent bodies, 261 - Gardiol et al. (2019) Gardiol, D., Affaticati, F., Albani, M., et al. 2019, in Proceedings of the 37th International Meteor Conference Pezinok-Modra, Slovakia, 30 August - 2 September 2018, ed. R. Rudawska, J. Rendtel, C. Powell, R. Lunsford, C. Verbeeck, & A. Knöfel, 81–86 - Gardiol et al. (2021) Gardiol, D., Barghini, D., Buzzoni, A., et al. 2021, Monthly Notices of the Royal Astronomical Society, 501, 1215 - Gardiol et al. (2016) Gardiol, D., Cellino, A., & Di Martino, M. 2016, in International Meteor Conference Egmond, the Netherlands, 2-5 June 2016, ed. A. Roggemans & P. Roggemans, 76 - Gerasimov et al. (2015) Gerasimov, S., Erofeev, V., Kikeev, V., et al. 2015, Nauchnaya vizualizaciya [Scientific Visualization], 7, 44 - Gowen & Perkins (1952) Gowen, F. E. & Perkins, E. W. 1952, Drag of Circular Cylinders for a Wide Range of Reynolds Numbers and Mach Numbers. Technical report, Ames Aeronautical Laboratory, Moffett Field, California., Tech. rep. - Gritsevich & Koschny (2011) Gritsevich, M. & Koschny, D. 2011, Icarus, 212, 877 - Gritsevich (2007) Gritsevich, M. I. 2007, Solar System Research, 41, 509 - Gritsevich (2009) Gritsevich, M. I. 2009, Advances in Space Research, 44, 323 - Halliday et al. (1996) Halliday, I., Griffin, A. A., & Blackwell, A. T. 1996, Meteoritics & Planetary Science, 31, 185 - Hansche & Rinehart (1952) Hansche, G. E. & Rinehart, J. S. 1952, Journal of the Aeronautical Sciences, 19, 83 - Henderson (1976) Henderson, C. B. 1976, AIAA journal, 14, 707 - Hill et al. (2004) Hill, K., Rogers, L., & Hawkes, R. 2004, Earth, Moon, and Planets, 95, 403 - Hill et al. (2005) Hill, K., Rogers, L., & Hawkes, R. 2005, Astronomy & Astrophysics, 444, 615 - Hodges (1957) Hodges, A. J. 1957, Journal of the Aeronautical Sciences, 24, 755 - Jeanne et al. (2019) Jeanne, S., Colas, F., Zanda, B., et al. 2019, A&A, 627, A78 - Johnston et al. (2018) Johnston, C. O., Stern, E. C., & Wheeler, L. F. 2018, Icarus, 309, 25 - Knapic et al. (2014) Knapic, C., De Marco, M., Smareglia, R., & Molinaro, M. 2014, in Astronomical Society of the Pacific Conference Series, Vol. 485, Astronomical Data Analysis Software and Systems XXIII, ed. N. Manset & P. Forshay, 131 - Koten & Borovicka (2001) Koten, P. & Borovicka, J. 2001, in Meteoroids 2001 Conference, Vol. 495, 259–264 - Lyytinen & Gritsevich (2016) Lyytinen, E. & Gritsevich, M. 2016, Planet. Space Sci., 120, 35 - McKinley (1961) McKinley, D. W. R. 1961, New York - McMullan & Collins (2019) McMullan, S. & Collins, G. S. 2019, Icarus, 327, 19 - Moreno-Ibáñez et al. (2018) Moreno-Ibáñez, M., Silber, E. A., Gritsevich, M., & Trigo-Rodríguez, J. M. 2018, The Astrophysical Journal, 863, 174 - Murray et al. (1999) Murray, I. S., Hawkes, R. L., & Jenniskens, P. 1999, Meteoritics & Planetary Science, 34, 949 - Opik (1958) Opik, E. J. 1958, New York - Passey & Melosh (1980) Passey, Q. & Melosh, H. 1980, Icarus, 42, 211 - Peña-Asensio et al. (2024) Peña-Asensio, E., Visuri, J., Trigo-Rodríguez, J. M., et al. 2024, Icarus, 408, 115844 - Plane (2012) Plane, J. M. 2012, Chemical Society Reviews, 41, 6507 - Rogers et al. (2005) Rogers, L. A., Hill, K. A., & Hawkes, R. L. 2005, Planet. Space Sci., 53, 1341 - Rouse (1946) Rouse, H. 1946, Elementary mechanics of fluids, 1st edn. (Dover Publications, New York) - Seltner et al. (2019) Seltner, P. M., Willems, S., & Gülhan, A. 2019, Journal of Spacecraft and Rockets, 56, 1725 - Seltner et al. (2021) Seltner, P. M., Willems, S., Gülhan, A., et al. 2021, Experiments in Fluids, 62, 182 - Sighard (1965) Sighard, F. H. 1965, Hoerner Fluid Dynamics - Spearman (1991) Spearman, M. 1991, in 11th Aerodynamic Decelerator Systems Technology Conference, 894 - Stulov et al. (1995) Stulov, V., Mirskii, V., & Vislyi, A. 1995, Moscow: Science. Fizmatlit - Subasinghe & Campbell-Brown (2018) Subasinghe, D. & Campbell-Brown, M. 2018, AJ, 155, 88 - Vennard (1940) Vennard, J. K. 1940, Elementary fluid mechanics, 1st edn. (Wiley, New York) - Verniani (1965) Verniani, F. 1965, Smithsonian Contributions to Astrophysics, 8, 141 - Vinković (2007) Vinković, D. 2007, Advances in Space Research, 39, 574 - Vondrak et al. (2008) Vondrak, T., Plane, J., Broadley, S., & Janches, D. 2008, Atmospheric Chemistry and Physics, 8, 7015 - Weryk & Brown (2013) Weryk, R. J. & Brown, P. G. 2013, Planetary and Space Science, 81, 32 - Zhdan et al. (2007) Zhdan, I. A., Stulov, V. P., Stulov, P. V., & Turchak, L. I. 2007, Solar System Research, 41, 505