Detection of meteors and sub-relativistic dust grains by the fluorescence detectors of ultra high energy cosmic rays

Fluorescence detectors of ultra high energy cosmic rays (UHECR) allow to record not only the extensive air showers, initiated by the UHECR particles, but also to detect light, produced by meteors and by the fast dust grains. It is shown that the fluorescence detector operated at the mountain site can register signals from meteors with kinetic energy threshold of about 25 J (meteor mass ∼ 5 × 10⁻⁶ g, velocity ∼ 3 × 10⁶ cm/s). The same detector might be used for recording of the dust grains of smaller mass (of about 10⁻¹⁰ g) but with velocity 10⁹ cm/s, close to the light velocity (sub-relativistic dust grains). The light signal from a sub-relativistic dust grain is expected in much shorter time scale (∼0.001 s), in comparison with the meteor signal (∼0.1–1 s), and much longer than duration of the UHECR signals (tens of μs). The fluorescence detector capable to register various phenomena: from meteors to UHECR – should have a variable pixel and selecting system integration time. A study of the new phenomenon of sub-relativistic grains will help to understand the mechanism of particle and dust grain acceleration in astrophysical objects (in SN explosions, for example).     

Velocity dependent response of a piezoelectric element to hypervelocity microparticles

Piezoelectric lead–zirconate–titanate (PZT) elements were studied by directly bombarding hypervelocity silver and iron particles in the velocity range of 2–63 km/s. This report focuses on the concept of “immediately after collision”. This concept was realized in the leading part of a pulsed signal that exhibited no effects due to reflection. The signal form was discussed from the viewpoint that information on impact was recorded by the waveform. The results indicated that the output amplitude was independent of the thickness of a PZT element, and that the waveform was explicitly related to the velocity at collision. Further, it appeared as an oscillating form at velocities less than 6 km/s. As the velocity increased, it gradually changed to a solitary pulse. The rise time of the single solitary waveform was linearly related to velocities above 6 km/s. The potential of a single PZT element is discussed as a real-time detector for hypervelocity microparticles.     

Redefinition of the meteor storm from the point of view of spaceflight security

ZHR¹ is defined as number of meteoroids passing 1000 km² zenith area per hour which can produce craters no less than 1 cm in diameter on an aluminum material. The relationship between ZHR and ZHR¹ is deduced. We evaluate the strong meteor showers since 1990s by ZHR¹ instead of ZHR and find the Giacobinid (Draconid) shower in 1998 is much stronger than the Leonid shower in 1999, 2001 and 2002.     

On the production of highest energy solar protons at 20 January 2005

On January 20, 2005, 7:02–7:05 UT the Aragats Multidirectional Muon Monitor (AMMM) located at 3200 m a.s.l. registered enhancement of the high energy secondary muon flux (threshold ∼5 GeV). The enhancement, lasting for 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES, and very fast (>2500 km/s) CME seen by SOHO, and the Ground Level Enhancements (GLE) #69 detected by the world-wide network of neutron monitors and muon detectors. The energetic and temporal characteristics of the muon signal from the AMMM are compared with the characteristics of other monitors located at the Aragats Space-Environmental Center (ASEC) and with other neutron and muon detectors. Since secondary muons with energies >5 GeV are corresponding to solar proton primaries with energies 20–30 GeV we conclude that in the episode of the particle acceleration at 7:02–7:05 UT 20 January 2005 solar protons were accelerated up to energies in excess of 20 GeV.     

Submicron lunar ejecta in the magnetosphere associated with meteor showers

Recent studies of the lunar ejecta from lunar impacts of interplanetary dust particles indicate that during favorable lunar phases, over 80% of the submicron ejecta enters the earth's magnetosphere. This “pulse” of lunar ejecta produced by the sporadic meteor background will follow the random variations of the sporadic flux. An additional enhancement of this flux can be related to major meteor showers. Since the annual periods of these showers occur during varying lunar phase angles, magnetosphere ejecta flux associated with major showers will vary depending on the coincidence of shower periods and favorable lunar phase angles. The results of an analysis of the “pulse” of ejecta flux in the magnetosphere during the Quadrantids, Geminids, Leonids, and Perseids meteor showers are presented. These results are compared to the satellite measurements of 1959 Eta and HEOS II.     

The 16-day waves in the mesosphere and lower thermosphere over Wuhan (30.6°N, 114.5°E) and Adelaide (35°S, 138°E)

Winds from a meteor radar at Wuhan (30.6°N, 114.5°E) and a MF radar at Adelaide (35°S, 138°E) are used to study the 16-day waves in the mesosphere and lower thermosphere (MLT). The height range is 78–98 km at Wuhan and 70–98 km at Adelaide. By comparison, it is found that the zonal components at both sites are generally larger than the meridional ones, and eastward motion of the zonal background winds is favorable for the 16-day waves penetration to the MLT region. The zonal maximum amplitude appears in the autumn (September–October) around 86–98 km at Wuhan and in the winter months and early spring (July–October) around 72–82 km at Adelaide. Differences are found in wave amplitudes and time of appearance between the two years of 2002 and 2003. In 2003, the intensity of the wave amplitudes is relatively smaller than that for 2002 at both sites. The summer 16-day waves are comparatively weaker at Adelaide in both years, but stronger in 2002 at Wuhan near the mesopause and the lower thermosphere (86–98 km). The strong summer waves at Wuhan may come from the winter southern hemisphere.     

Development of a 2.8 μm film for scientific balloons

Development of a balloon to fly at higher altitudes is one of the most attractive challenges in scientific balloon technologies. After reaching the highest record setting balloon altitude of 53.0 km using the 3.4 μm film in 2002, we tried to make a thinner balloon film. In 2003, we developed a forming die and an air-ring and succeeded in forming a film with a thickness of 3.0 μm and a width of 220 cm. Using this film, we manufactured a balloon with a volume of 5000 m³ and succeeded in flying the balloon up to an altitude of 46.0 km. We then searched for a good combination of resins to make a thinner and wider film and obtained films with widths of 280 cm, and a thickness of 3.0 μm at first, and then 2.8 μm. In 2004, we performed balloon experiments making a 30,000 m³ balloon with the 3.0 μm film and a 5000 m³ balloon with the 2.8 μm film. Both balloons were well manufactured and reached the highest altitudes of 50.7 and 42.6 km, respectively.     

Optical properties of dust and the opacity of the Martian atmosphere

Particulate component of the Mars atmosphere composed by micron-sized products of soil weathering and water ice clouds strongly affects the current climate of the planet. In the absence of a dust storm so-called permanent dust haze with τ ≈ 0.2 in the atmosphere of Mars determines its thermal structure. Dust loading varies substantially with the season and geographic location, and only the data of mapping instruments are adequate to characterize it, such as TES/MGS and IRTM/Viking. In spite of vast domain of collected data, no model is now capable to explain all observed spectral features of dust aerosol. Several mineralogical and microphysical models of the atmospheric dust have been proposed but they cannot explain the pronounced systematic differences between the IR data (τ = 0.05–0.2) and measurements from the surface (Viking landers, Pathfinder) which give the typical “clear” optical depth of τ ≈ 0.5 from one side, and ground-based observations in the UV–visible range showing much more transparent atmosphere, on the other side. Also the relationship between τ₉ and the visible optical depth is not well constrained experimentally so far. Future focused measurements are therefore necessary to study Martian aerosol.     

The superbubble origin for galactic cosmic rays

The source of energy for cosmic-ray acceleration appears to be shock waves driven by supernova (SNe) ejecta. The great majority (80–90%) of SNe (SNII and SNIb) are formed by the core collapse of young, massive O and B stars. However, it has been known for more than forty years that the births of such massive stars in stellar clumps, termed OB associations, are correlated in space and time. The combined ejecta of core-collapse SNe, occurring at the deaths of these massive stars, create low-density (∼5.0 × 10⁻³ cm⁻³) superbubbles that reach dimensions of several hundred pc. The occurrence of correlated SNe in superbubbles affects not just the source of cosmic-ray energy, SNe shock waves, it impinges as well on the elemental and isotopic source abundances of cosmic-ray nuclei. We argue that the well-known anomalous cosmic-ray ²²Ne/²⁰Ne ratio, a factor of five times the Solar System ratio, results from a mixing of freshly synthesized nucleosynthetic material in supernova active cores of superbubbles. Although diluted by mixing with older, lower metallicity interstellar gas, the mean metallicities in the superbubble, SNe-dominated cores are high ∼3 times the Solar System value.     

Analysis of electromagnetic and electrostatic effects of particle impacts on spacecraft

Particle impacts on spacecraft can cause considerable damage, even leading to complete failure. A theory for the resulting vehicle potential changes and the electromagnetic radiation from impact-induced plasma has been published by Close et al. (2010). Here we compare this theory to impacts registered by the Radio and Plasma Wave Science instrumentation on the Cassini spacecraft. We study both low-velocity (16 km/s) large particles (2.6 μm radius) detected in Saturn’s rings and high-velocity (450 km/s) small particles (1 nm radius) in the solar wind. The agreement with the theory is quite good. We also apply these results to earth orbit and conclude that both Electrostatic Discharge and Electromagnetic Pulse radiation are likely and could lead to spacecraft failure.     

On the quality of mutations in mammalian cells induced by high LET radiations

The deleterious effects of accelerated heavy ions as component of the space radiation environment on living cells are of increasing importance for long duration human space flight activities. The most important aspect of such densely ionizing particle radiation is attributed to the type and quality of biological damage induced by them. This issue is addressed by investigating cell inactivation and mutation induction at the Hprt locus (coding for hypoxanthine-guanine-phosphoribosyl-transferase) of cultured V79 Chinese hamster cells exposed to densely ionizing radiation (accelerated heavy ions with different LETs from oxygen to gold, specific energies ranging from 1.9 to 69.7 MeV/u, corresponding LET values range from 62 to 13,223 keV/μm) and to sparsely ionizing radiation (200 kV X-rays). 30 spontaneous, 40 X-ray induced and 196 heavy ion induced 6-thioguanine resistant Hprt mutant colonies were characterized by Southern technique using the restriction enzymes EcoRI, PstI and BglII and a full length Hprt cDNA probe isolated from the plasmid pHPT12. Restriction patterns of the spontaneous Hprt mutants were indistinguishable from the wild type pattern, as these mutants probably contain only small deletions or even point mutations in the Hprt locus. In contrast, the overall spectrum of heavy ion induced mutations revealed a majority of partial or total deletions of the Hprt gene. With constant particle fluence (3 × 10⁶ particles/cm²) the quality of heavy ion induced mutations in the Hprt locus depends on physical parameters of the beam (atomic number, specific energy, LET). This finding suggests a relationship between the type of DNA damage and track structure. The fraction of mutants with severe deletions in the Hprt locus after exposure to oxygen ions increases from 65% at 60 keV/μm up to a maximum (100%) at 300 keV/μm and declines with higher LET values to 75% at 750 keV/μm. With heavier ions (Ca- and Au-ions) and even higher LET-values this mutant fraction decreases to 58% at 13,200 keV/μm. Heavy ion induced DNA break points in the Hprt locus are not randomly distributed.     

Hard X-ray footpoint motions in solar flares: Comparing magnetic reconnection models with observations

Hard X-ray observations from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) of the October 29, 2003 GOES X10 two-ribbon flare are used together with magnetic field observations from the Michelson Doppler Imager (MDI) onboard SoHO to compare footpoint motions with predictions from magnetic reconnection models. The temporal variations of the velocity v of the hard X-ray footpoint motions and the photospheric magnetic field strength B in footpoints are investigated. The underlying photospheric magnetic field strength is generally higher (B ∼ 700–1200 G) in the slower moving (v ∼ 20–50 km s⁻¹) western footpoint than in the faster (v ∼ 20–100 km s⁻¹) moving eastern source (∼100–600 G). Furthermore, a rough temporal correlation between the HXR flux and the product vB² is observed.     

Detection of thermal X-ray emission in the halo of the plerionic supernova remnant G21.5-0.9

The detection of a soft thermal X-ray component in the spectrum of a bright knot in the halo of the plerion G21.5-0.9 is reported. Using a collisional ionization equilibrium model for an hot optically thin plasma, a temperature of kT 0.12–0.24 keV, a mass of 0.3–1.0 M_⊙ and a density of 1.6–6 cm⁻³ is derived. The spectral analysis suggests a possible overabundance of Silicon with respect to the solar value in the knot; if this will be confirmed this object may be a clump of shocked ejecta.     

Progress in laboratory performance verification of the R-process Isotope Observer

The source of galactic cosmic rays (GCR) remains unknown, despite many decades of research. It is generally accepted among high-energy astrophysicists that GCR are accelerated by supernova (SN) shocks in the interstellar medium, but the evidence for this remains indirect. However, the fact that the cosmic-ray spectrum continues far past the limit of SN shock acceleration is a serious problem for this model. Most SNae occur in superbubbles – regions of the interstellar medium containing massive stars and SN remnants. Thus, SNae should be accelerating the ejecta of other SNae. The local medium of the superbubble will be enhanced in freshly synthesized R-process material. We are currently developing the R-process Isotope Observer (RIO), which could make the first measurements of the isotopic abundances of the “ultraheavy” GCR (those in the range 32 ⩽ Z ⩽ 42) in order to determine the fractional contribution of freshly synthesized R-process material in GCR through the measurement of several key isotopic ratios. We have recently performed a series of tests of the RIO detector with 400 A MeV ⁸⁴Kr at the HIMAC accelerator in Japan. The results of these tests indicate that a mass resolution of ∼0.3 u can be achieved with RIO. We also discuss the prospects for accommodation of the RIO mission.     

Airglow OH emission height inferred from the OH temperature and meteor trail diffusion coefficient

Simultaneous observations of the airglow OH(6,2) band rotational temperature, TOH, and meteor trail ambipolar diffusion coefficient, D, were carried out at Shigaraki (35°N, 136°E), during PSMOS 2003 Campaign, January 28 to February 8, 2003. The OH emission height was estimated by cross correlation analysis of the TOH and D nocturnal variations. A good correlation between TOH and D was obtained at 85 km of altitude. From the nocturnal variations of TOH and D, it is found that the OH emission peak height varied from 88 km before the midnight to 84 km in the early morning. The height variation could be caused by an atmospheric tidal effect in the emission height.     

Comparison of vertical E × B drift velocities and ground-based magnetometer observations of DELTA H in the low latitude under geomagnetically disturbed conditions

In the present work, we analyzed the daytime vertical E × B drift velocities obtained from Jicamarca Unattended Long-term Ionosphere Atmosphere (JULIA) radar and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca, and Piura in Peru for 22 geomagnetically disturbed events in which either SC has occurred or Dst_max < ?50 nT during the period 2006–2011. The ΔH component of geomagnetic field is measured as the differences in the magnitudes of horizontal H component between magnetometer placed directly on the magnetic equator and one displaced 6–9° away. It will provide a direct measure of the daytime electrojet current, due to the eastward electric field. This will in turn gives the magnitude of vertical E × B drift velocity in the F region. A positive correlation exists between peak values of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of the events. It was observed that 45% of the events have daytime vertical E × B drift velocity peak in the magnitude range 10–20 m/s and 20–30 m/s and 20% have peak ΔH in the magnitude range 50–60 nT and 80–90 nT. It was observed that the time of occurrence of the peak value of both the vertical E × B drift velocity and the ΔH have a maximum (40%) probability in the same time range 11:00–13:00 LT. We also investigated the correlation between E × B drift velocity and D_st index and the correlation between delta H and D_st index. A strong positive correlation is found between E × B drift and D_st index as well as between delta H and D_st Index. Three different techniques of data analysis – linear, polynomial (order 2), and polynomial (order 3) regression analysis were considered. The regression parameters in all the three cases were calculated using the Least Square Method (LSM), using the daytime vertical E × B drift velocity and ΔH. A formula was developed which indicates the relationship between daytime vertical E × B drift velocity and ΔH, for the disturbed periods. The E × B drift velocity was then evaluated using the formulae thus found for the three regression analysis and validated for the ‘disturbed periods’ of 3 selected events. The E × B drift velocities estimated by the three regression analysis have a fairly good agreement with JULIA radar observed values under different seasons and solar activity conditions. Root Mean Square (RMS) errors calculated for each case suggest that polynomial (order 3) regression analysis provides a better agreement with the observations from among the three.     

Spectral filter for signal identification in the kHz SLR measurements of the fast spinning satellite Ajisai

The high repetition rate satellite laser ranging (SLR) measurements to the fast spinning satellites contain a frequency signal caused by the rotational motion of the corner cube reflector (CCR) array. The spectral filter, developed here, is based on the Lomb algorithm, and is tested with the simulated and the observed high repetition rate SLR data of the geodetic satellite Ajisai (spin period ∼2 s). The filter allows for the noise elimination from the SLR data, and for identification of the returns from the single CCRs of the array – even for the low return rates. Applying the spectral filter to the simulated SLR data increases the S/N ratio by a factor 40–45% for all return rates. Filtering out the noise from the observed data strengthens the frequency signal by factor of ∼25 for the low return rates, which significantly helps to determine the spin phase of the satellite. The spectral filter is applied to the Graz SLR data and the spin rates of Ajisai are determined by two different methods: the frequency analysis and the phase determination of the spinning retroreflector array.The analysis of more than 8 years of the Graz SLR measurements indicates an exponential spin rate trend: f = 0.67034 exp(−0.0148542 Y) [Hz], RMS = 0.085 mHz, where Y is the year since launch. The highly accurate spin rate information demonstrates periodic changes related to the precession of the orbital plane of Ajisai, as it determines the amount of energy received by the satellite from the Sun. The rate of deceleration of Ajisai is not constant: the half life period of the satellite’s spin oscillates around 46.7 years with an amplitude of about 5 years.     

The neutral gas in the environs of the Geminga gamma-ray pulsar

We present a high-resolution (24″) study of the HI interstellar gas distribution around the radio-quiet neutron star Geminga. Based on very large array and MPIfR Effelsberg telescope data, we analyzed a 40′ × 40′ field around Geminga. These observations have revealed the presence of a neutral gas shell, 0.4 pc in radius, with an associated HI mass of 0.8M_⊙, which surrounds Geminga at a radial velocity compatible with the kinematical distance of the neutron star. In addition, morphological agreement is observed between the internal face of the HI shell and the brightest structure of Geminga’s tail observed in X-rays. We explore the possibility that this morphological agreement is the result of a physical association.     

Heliospheric implications of structure in the interstellar medium

Equilibrium models of diffuse interstellar material (ISM) near the Sun show a range of cloud densities, ionization, and temperatures which are consistent with data, although the local ISM must be inhomogeneous over ∼2 pc scales. The ISM close to the Sun has properties that are consistent with the sheetlike warm neutral (and partially ionized) gas detected in the Arecibo Millennium Survey. Local interstellar magnetic fields are poorly understood, but data showing weak polarization for nearby stars indicate dust may be trapped in fields or currents in the heliosheath nose region. Implications of this dust capture are widespread, and may impact the interpretation of the cosmic microwave background data. Observations of interstellar H⁰ inside of the solar system between 1975 and 2000 do not suggest any variation in the properties or structure of local interstellar H⁰ over distance scales of ∼750 AU to within the uncertainties.     

USV test flight by stratospheric balloon: Preliminary mission analysis

The Unmanned Space Vehicle test flights will use a 7 m 1300 kg aircraft. The first three launches will take place at the Italian Space Agency ASI base in Trapani–Milo, Sicily, through a stratospheric balloon that will drop the aircraft at a predefined height. After free fall acceleration to transonic velocities, the parachute deployment will allow a safe splash down in the central Mediterranean Sea. The goal of this article is to show the preliminary analysis results for the first USV flight.We carried out a statistical study for the year 2000–2003, evaluating the typical summer and winter launch windows of the Trapani–Milo base.First, in the center Mediterranean, we define safe recovery areas. They cannot be reached during the balloon ascending phase so, after a sufficiently long floating part able to catch the open sea, the balloon will go down to the release height (24 km). The simulation foresees a 400,000 m³ balloon and 3 valves for the altitude transfer.A safe splash down must occur far enough from the nearest coast: the minimum distance is considered around 25 km. The vehicle should be released at a distance, from the nearest coast, greater than this minimum amount plus the USV model maximum horizontal translation, during its own trajectory from balloon separation to splash down. In this way we define safe release areas for some possible translations.Winter stratospheric winds are less stable. The winter average flight duration is 7 h and it is probably too long for the diurnal recovery requirement and its scheduled procedures.Comparing past stratospheric balloons flights and trajectories computed using measured meteorological data (analysis), with their predictions made using forecast models and soundings, we obtain the standard deviation of the trajectory forecast uncertainty at the balloon–aircraft separation. Two cases are taken into account: predictions made 24 and 6 h before the launch.Assuming a Gaussian latitudinal uncertainty distribution for the prediction 6 h before the launch, we are able to identify the forecast trajectories that have a probability greater than 97% to reach the safe release areas.Simulating the summer windows trajectories for the years from 2000 to 2003 and for the favorable ground wind days, we obtain the number of trajectories with the desired forecast probabilities.