Hypervelocity impacts on HST solar arrays and the debris and meteoroids population

Accurate debris and meteoroid flux models are crucial for the design of manned and unmanned space missions. For the most abundant particle sizes smaller than a few millimetres, knowledge of the populations can only be gained from in situ detectors or the analysis of retrieved space hardware. The measurement of impact flux from exposed surfaces improves with increased surface area and exposure time.A post-flight impact investigation was initiated by the European Space Agency to record and analyse the impact fluxes and any potential resulting damage on the two flexible solar arrays of the Hubble Space Telescope. The arrays were deployed during the first Hubble Space Telescope servicing mission in December 1993 and retrieved in March 2002. They have a total exposed surface area of roughly 120 m², including 42 m² covered with solar cells. This new Hubble post-flight impact study follows a similar activity undertaken after the retrieval of one of the first solar arrays, in 1993. The earlier study provided the first opportunity for a numerical survey of damage to exposed surfaces from more than 600 km altitude, and of impacts from particles larger than 1 mm. The results have proven very valuable in validation of important flux model regimes. The second set of Hubble solar arrays has again provided an unrivalled opportunity to measure the meteoroid and debris environment, now sampled during a long interval in low Earth orbit, and to identify changes in the space debris environment since the previous survey. The retrieved solar array wings exhibit thousands of craters, many of which are visible to the naked eye. A few hundred impacts have completely penetrated the 0.7 mm thick array. The largest impact features are about 7–8 mm in diameter. The cover glass of the solar cells is particularly well suited to the recognition of small impact features by optical and electron microscopy. In this paper, we present the first results of the impact survey. Data upon the abundance of craters of specific measured size ranges are plotted as cumulative flux curves, and compared to the results of model predictions. The most significant change to the particle flux since 1993 is a decrease in the small debris population.     

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.     

MULPEX: A compact multi-layered polymer foil collector for micrometeoroids and orbital debris

Detailed studies of preserved hypervelocity impact residues on spacecraft multi-layer insulation foils have yielded important information about the flux of small particles from different sources in low-Earth orbit (LEO). We have extended our earlier research on impacts occurring in LEO to design and testing of a compact capture device. MUlti-Layer Polymer EXperiment (MULPEX) is simple, cheap to build, lightweight, of no power demand, easy to deploy, and optimised for the efficient collection of impact residue for analysis on return to Earth. The capture medium is a stack of very thin (8 and 40 μm) polyimide foils, supported on poly-tetrafluoroethylene sheet frames, surrounded by a protective aluminium casing. The uppermost foil has a very thin metallic coating for thermal protection and resistance to atomic oxygen and ultra-violet exposure. The casing provides a simple detachable interface for deployment on the spacecraft, facing into the desired direction for particle collection. On return to the laboratory, the stacked foils are separated for examination in a variable pressure scanning electron microscope, without need for surface coating. Analysis of impact residue is performed using energy dispersive X-ray spectrometers. Our laboratory experiments, utilising buck-shot firings of analogues to micrometeoroids (35–38 μm olivine) and space debris (4 μm alumina and 1 mm stainless steel) in a light gas gun, have shown that impact residue is abundant within the foil layers, and preserves a record of the impacting particle, whether of micrometer or millimetre dimensions. Penetrations of the top foil are easily recognised, and act as a proxy for dimensions of the penetrating particle. Impact may cause disruption and melting, but some residue retains sufficient crystallographic structure to show clear Raman lines, diagnostic of the original mineral.     

A search for radio emission from the young 16-ms X-ray pulsar PSR J0537−6910

PSR J0537−6910 is a young, energetic, rotation-powered X-ray pulsar with a spin period of 16 ms located in the Large Magellanic Cloud. We have searched for previously undetected radio pulsations (both giant and standard) from this pulsar in a 12-h observation taken at 1400 MHz with the Parkes 64-m radio telescope. The very large value of the magnetic field at the light cylinder radius suggests that this pulsar might be emitting giant radio pulses like those seen in other pulsars with similar field strengths. No radio emission of either kind was detected from the pulsar, and we have established an upper limit of ∼25 mJy kpc² for the average 1400-MHz radio luminosity of PSR J0537−6910. The 5σ single-pulse detection threshold was ∼750 mJy for a single 80-μs sample. These limits are likely to be the best obtainable until searches with greatly improved sensitivity can be made with next-generation radio instruments.     

The hard X-ray emission from the Vela pulsar wind nebula

We present the results of a preliminary spectral analysis performed on the BeppoSAX and XMM observations of the Vela plerion. The broad energy range covered by the instruments on board the two observatories allows an evaluation of the spectral parameters of the high energy emission model and provides an indication on the morphology of the source emission above 10 keV. We confirm the softening of the PWN spectrum (3–10 keV band) at distances greater than 4′ from the pulsar and estimate the diameter of the high energy (>10 keV) emission region to be on the order of 25′–30′.     

High energy emission from pulsars: Outer gap scenario

I will give a brief review of the recent development in the emission models of isolated, rapidly rotating neutron stars, focusing on the γ-ray radiation mechanism in their outer magnetospheres. By examining the Poisson equation for the electrostatic potential, I show that an active particle accelerator must extend from the vicinity of the neutron star surface to the vicinity of light cylinder. Furthermore, combining the Poisson equation with the Boltzmann equations for electrons/positrons and γ-rays, and assuming that the gap trans-field thickness is large compared to the longitudinal width, I demonstrate that the energy distribution of ultra-relativistic particles cannot be described by a power-law but by a quasi-monoenergetic distribution at the terminal Lorentz factor. The particles are accelerated in the gap and escape from it with large Lorentz factors. Is is shown that such energetic particles migrating outside of the gap contribute significantly to the γ-ray luminosity and reproduce the observed soft γ-ray spectrum between 100 MeV and 3 GeV for the Vela pulsar.     

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).     

The role of drift and convection–diffusion mechanisms in small energy global cosmic ray modulation: Application to the hysteresis phenomenon of proton and alpha particle satellite data

The hysteresis effect for small energies of galactic cosmic rays is due to two effects. The first is the same as for neutron monitor energies – the delay of the interplanetary processes responsible for cosmic ray modulation with respect to the initiating solar processes, according to the effective velocity of solar wind and shock waves propagation. Then, the observed cosmic ray intensity is connected to the solar activity variations during many months before the time of cosmic ray measurement. The second is caused by the time delay of small energy cosmic ray diffusion from the boundary of modulation region to the Earth’s orbit. The model describing the connection between solar activity variation and cosmic ray convection–diffusion global modulation for neutron monitor energies is here developed by taking into account also the time-lag of the small energy particle diffusion in the Heliosphere. We use theoretical results on drifts and analytically approximate the dependences of drifts from tilt angle, and take into account the dependence from the sign of primary particles, and from the sign of polar magnetic field (A > 0 or A < 0). The obtained results are applied on proton and alpha-particle satellite data. We analyze satellite 5-min data of proton fluxes with energies >1 MeV, >2 MeV, >5 MeV, >10 MeV, >30 MeV, >50 MeV, >60 MeV, >100 MeV, and in intervals 10–30 MeV, 30–60 MeV, and 60–100 MeV during January 1986–December 1999. We exclude periods with great cosmic ray increases caused by particle acceleration in solar flare events. Then, we determine monthly averaged fluxes, as well as 5-month and 11-month smoothed data. We analyze also satellite 5-min data on alpha-particle fluxes in the energy intervals 60-160 MeV, 160–260 MeV and 330–500 MeV during January 1986–May 2000. We correct observation data for drifts and then compare with what is expected according to the convection–diffusion mechanism. We assume different dimensions of the modulation region (by the time propagation X₀ of solar wind from the Sun to the boundary of modulation region), for X₀ values from 1 to 60 average months, by one-month steps. For each value of X₀ we determine the correlation coefficient between variations of expected and observed cosmic ray intensities (the estimation of cosmic ray intensities values is given in Section 3 by Eq. (9), and the determination of correlation and regression coefficients in Section 3 by Eq. (8)). The dimension of modulation region is determined by the value of X_0 max, for which the correlation coefficient reaches the maximum value. Then the effective radial diffusion coefficient and residual modulation in small energy region can be estimated.     

Simulation study of motion of charged particles trapped in Earth’s magnetosphere

This article aims to understand the motion of the charged particles trapped in the Earth’s inner magnetosphere. The emphasis is on identifying the numerical scheme, which is appropriate to characterize the trajectories of the charged particles of different energies that enter the Earth’s magnetosphere and get trap along the magnetic field lines. These particles perform three different periodic motions, namely: gyration, bounce, and azimuthal drift. However, often, the gyration of the particle is ignored, and only the guiding center of the particle is traced to reduce the computational time. It is because the simulation of all three motions (gyro, bounce, and drift) together needed a robust numerical scheme, which has less numerical dissipation. We have developed a three-dimensional test particle simulation model in which the relativistic equation of motion is solved numerically using the fourth and sixth-order Runge-Kutta methods. The stability of the simulation model is verified by checking the conservation of total kinetic energy and adiabatic invariants linked with each type of motion. We found that the sixth-order Runge-Kutta method is suitable to trace the complete trajectories of both proton and electron of a wide energy range, 5 keV to 250 MeV for L = 2 – 6. We have estimated the bounce and drift periods from the simulations, and they are found to be in good agreement with the theory. The study implies that a simulation model with sixth-order Runge-Kutta method can be applied to the time-vary, non-analytical form of magnetic configuration in future studies to understand the dynamics of charged particles trapped in Earth’s magnetosphere.     

An analysis of solar radio burst events on December 1, 2004

We present the analysis of the radio observations of December 1, 2004 from 07:00 UT to 07:40 UT in the 1.100–1.340 GHz band by Solar Broadband Radio Dynamic Spectrometer (SBRS) in Huairou Station. There are three groups of radio fine structures during the impulsive phase of this flare denoted by N1, Z2, and Z3. N1 has several emission lines with mixed fast and slow frequency drift rate which may reflect the conditions of flare loop and fast flows out from reconnection site; Z2 and Z3 are zebra patterns. The radio observations combined with hard X-ray and other observations show that the fine structures are connected with energetic particles. The information about magnetic field and energetic particle during the burst are also estimated based on our model.     

Real-time space debris monitoring with EISCAT

Following a feasibility study in 2000–2001 on using the EISCAT ionospheric research radars to detect centimetre-sized space debris in the frame of an ESA contract, we are now finishing a continuation study, aimed at achieving debris detection and parameter estimation in real-time. A requirement is to “piggy-back” space debris measurements on top of EISCAT’s normal ionospheric work, without interfering with that work, and to be able to handle about 500 h of measurements per year. We use a special digital receiver back-end in parallel with EISCAT’s standard receiver. We sample fast enough to correctly band-pass sample the EISCAT analog frequency band. To increase detection sensitivity, we use coherent pulse-to-pulse integration. The coherent integration is built-in in our method of parameter estimation, which we call the match function (MF) method. The method is derived from Bayesian statistical inversion, but reduces, with standard assumptions about noise and prior, to minimizing the least squares norm ∥z(t) − bχ(R,v,a;t)∥, where z is the measured signal and {bχ} is a set of model signals. Because the model signals depend linearly on the amplitude b, it is sufficient to maximize the magnitude of the inner product (cross correlation) between z and χ, the amplitude estimate is then determined by direct computation. The magnitude of the inner product, when properly normalized, is the MF. To construct the set of model signals, we sample the EISCAT transmission, in the same way as we sample the received signal, and apply linearly changing Doppler-shifts to it. Our initial implementation of the MF-method in 2001 was about four orders of magnitude too slow for real-time applications, but we have now gained the required speed factors. A factor of ten comes from using faster computers, another factor of ten comes from coding our key algorithms in C instead of Matlab. The largest factor, typically 100–300, comes from using a special, approximative, but in practice quite sufficient, method of finding the MF maximum. Test measurements show that we get real-time speed already when using a single dual-processor 2 GHz G5 Macintosh to do the detection computations.     

Multi-year behaviour of the midnight OH1 temperature according to observations at Zvenigorod over 2000–2016

Using spectral measurements of the hydroxyl airglow at the Zvenigorod station (56° N, 37° E), Moscow region, over 2000–2016, we obtained the long-term set of data comprising 1822 midnight values of the OH¹ temperature in the mesopause region. These data revealed a 17-year series of its mean annual values, as well as amplitudes and phases of the first two harmonics of its annual variation. The obtained parameters were analyzed to determine statistically relevant characteristics of their long-term variations. As a result, we found that the long-term behaviour of the mean annual OH¹ temperature features a small negative linear trend (?0.07 ± 0.03 K/year) over the addressed period. Besides, its dependence on solar activity is shown to be 4.1 ± 0.5 K/100 SFU. Regarding the long-term behaviour of the mean annual OH¹ temperature, we revealed the existence of two oscillations with 3-year (the amplitude being 1.3 ± 0.2 K) and 4.1-year (the amplitude being 0.6 ± 0.2 K) periods. We obtained empirical relations describing year-to-year variations in the amplitudes and phases of the annual and semi-annual harmonics.     

Possible DC electric field in the ionosphere related to seismicity

This paper presents the method for calculation of DC electric field in the atmosphere and the ionosphere generated by model distribution of external electric current in the lower atmosphere. Appearance of such current is associated with enhancement of seismic activity that is accompanied by emanation of soil gases into the atmosphere. These gases transfer positive and negative charged aerosols. Atmospheric convection of charged aerosols forms external electric current, which works as a source of conductivity current in the atmosphere–ionosphere electric circuit. It is shown that DC electric field generated in the ionosphere by this current reaches up to 10 mV/m, while the long-term vertical electric field disturbances excited near the Earth surface do not exceed 100 V/m. Such limitation of the near-ground field is caused by the formation of potential barrier for charged particles at the Earth surface in a process of their transport from soil to atmosphere.     

Hard X-ray component in Sco X-1 spectrum: Synchrotron emission from a nano quasar

Sco X-1 is a low mass X-ray binary system and with the recent observations of a resolved radio jet, the source has been included in the list of galactic microquasars. The observed spectral data in the 2–20 keV energy band fits a thermal emission. Above 20 keV, a hard tail has been reported on occasions. During our continuing balloon borne X-ray survey in the 20–200 keV region using high sensitivity Large Area Scintillation counter Experiment, Sco X-1 was observed on two different occasions. Even though the total X-ray luminosity of the source was different, the spectral nature of the source did not show any variation. The presence of hard X-ray flux is unmistakable. We present the spectral and temporal data in the hard X-ray band and discuss the results in terms of geometrical characteristics of X-ray source and its observed temporal properties. We note that the jet activity is similar to the microquasars, however, the absence of the large magnitude abrupt changes in X-ray light curve compared to GRS1915 + 105 suggest that the quasar-like behaviour is at a nano scale.     

Comparison of genotoxic damage in monolayer cell cultures and three-dimensional tissue-like cell assemblies

Assessing the biological risks associated with exposure to the high-energy charged particles encountered in space is essential for the success of long-term space exploration. Although prokaryotic and eukaryotic cell models developed in our laboratory and others have advanced our understanding of many aspects of genotoxicity, in vitro models are needed to assess the risk to humans from space radiation insults. Such models must be representative of the cellular interactions present in tissues and capable of quantifying genotoxic damage. Toward this overall goal, the objectives of this study were to examine the effect of the localized microenvironment of cells, cultured as either 2-dimensional (2D) monolayers or 3-dimensional (3D) aggregates, on the rate and type of genotoxic damage resulting from exposure to Fe-charged particles, a significant portion of space radiation. We used rodent transgenic cell lines containing 50–70 copies of a LacI transgene to provide the enhanced sensitivity required to quantify mutational frequency and type in the 1100-bp LacI target as well as assessment of DNA damage to the entire 45-kbp construct. Cultured cells were exposed to high-energy Fe charged particles at Brookhaven National Laboratory’s Alternating Gradient Synchrotron facility for a total dose ranging from 0.1 to 2 Gy and allowed to recover for 0–7 days, after which mutational type and frequency were evaluated. The mutational frequency was found to be higher in 3D samples than in 2D samples at all radiation doses. Mutational frequency also was higher at 7 days after irradiation than immediately after exposure. DNA sequencing of the mutant targets revealed that deletional mutations contributed an increasingly high percentage (up to 27%) of all mutations in cells as the dose was increased from 0.5 to 2 Gy. Several mutants also showed large and complex deletions in multiple locations within the LacI target. However, no differences in mutational type were found between the 2D and the 3D samples. These 3D tissue-like model systems can reduce the uncertainty involved in extrapolating risk between in vitro cellular and in vivo models.     

Collaborative VLA,SOHO and RHESSI observations of evolving sources of energy release in the corona above active regions

Very Large Array (VLA) observations at 20 and 91 cm wavelength are compared with data from the SOHO (EIT and MDI) and RHESSI solar missions to investigate the evolution of decimetric Type I noise storms and Type III bursts and related magnetic activity in the photosphere and corona. The combined data sets provide clues about the mechanisms that initiate and sustain the decimetric bursts and about interactions between thermal and nonthermal plasmas at different locations in the solar atmosphere. On one day, frequent, low-level hard X-ray flaring observed by RHESSI appears to have had no clear affect on the evolution of two closely-spaced Type I noise storm sources lying above the target active region. EIT images however, indicate nearly continuous restructuring of the underlying EUV loops which, through accompanying low-level magnetic reconnection, might give rise to nonthermal particles and plasma turbulence that sustain the long-lasting Type I burst emission. On another day, the onset of an impulsive hard X-ray burst and subsequent decimetric burst emission followed the gradual displacement and coalescence of a small patch of magnetic magnetic polarity with a pre-existing area of mixed magnetic polarity. The time delay of the impulsive 20 and 91 cm bursts by up to 20 min suggests that these events were unlikely to represent the main sites of flare electron acceleration, but instead are related to the rearrangement of the coronal magnetic field after the main flare at lower altitude. Although the X-ray flare is associated with the decimetric burst, the brightness and structure of a long-lasting Type I noise storm from the same region was not affected by the flare. This suggests that the reconfiguration of the coronal magnetic fields and the subsequent energy release that gave rise to the impulsive burst emission did not significantly perturb that part of the corona where the noise storm emission was located.     

The accuracy of SST retrievals from AATSR: An initial assessment through geophysical validation against in situ radiometers,buoys and other SST data sets

The Advanced Along-Track Scanning Radiometer (AATSR) was launched on Envisat in March 2002. The AATSR instrument is designed to retrieve precise and accurate global sea surface temperature (SST) that, combined with the large data set collected from its predecessors, ATSR and ATSR-2, will provide a long term record of SST data that is greater than 15 years. This record can be used for independent monitoring and detection of climate change. The AATSR validation programme has successfully completed its initial phase. The programme involves validation of the AATSR derived SST values using in situ radiometers, in situ buoys and global SST fields from other data sets. The results of the initial programme presented here will demonstrate that the AATSR instrument is currently close to meeting its scientific objectives of determining global SST to an accuracy of 0.3 K (one sigma). For night time data, the analysis gives a warm bias of between +0.04 K (0.28 K) for buoys to +0.06 K (0.20 K) for radiometers, with slightly higher errors observed for day time data, showing warm biases of between +0.02 (0.39 K) for buoys to +0.11 K (0.33 K) for radiometers. They show that the ATSR series of instruments continues to be the world leader in delivering accurate space-based observations of SST, which is a key climate parameter.     

Modulations of broad-band radio continua and X-ray emissions in the large X-ray flare on 03 November 2003

The GOES X3.9 flare on 03 November 2003 at ∼09:45 UT was observed from metric to millimetric wavelengths by the Nançay Radioheliograph (NRH), the Radio Solar Telescope Network (RSTN) and by radio instruments operated by the Institute of Applied Physics (University of Bern). This flare was simultaneously observed and imaged up to several 100 keV by the RHESSI experiment. The time profile of the X-ray emission above 100 keV and of the radio emissions shows two main parts, impulsive emission lasting about 3 min and long duration emission (partially observed by RHESSI) separated in time by 4 min. We shall focus here on the modulations of the broad-band radio continua and of the X-ray emissions observed in the second part of the flare. The observations suggest that gyrosynchrotron emission is the prevailing emission mechanism even at decimetric wavelengths for the broad-band radio emission. Following this interpretation, we deduce the density and the magnetic field of the decimetric sources and briefly comment on possible interpretations of the modulations.     

SCIAMACHY solar irradiance observation in the spectral range from 240 to 2380 nm

The SCanning Imaging Absorption Spectrometer for Atmospheric CHartographY (SCIAMACHY) is part of the payload of ESA’s Environmental Satellite ENVISAT which was launched into a sun-synchronous polar orbit on 2002-03-01. It is the first spaceborne instrument covering a wavelength range from 240 to 2380 nm thus including ultraviolet, visible and near infrared spectral regions.The main purpose of SCIAMACHY is to determine the amount and distribution of a large number of atmospheric trace constituents by measuring the radiance backscattered from the Earth. In addition, several solar observations are performed with daily or orbital frequency.The presented results will cover the following topics: (a) comparison of the solar irradiance measured by SCIAMACHY with data from the instruments SOLSPEC/SOLSTICE/SUSIM and a solar spectrum derived by Kurucz; (b) comparison of the SCIAMACHY solar Mg II index with GOME and NOAA data; (c) correlation of the relative change of solar irradiance measured by SCIAMACHY with the sun spot index.The mean solar irradiance for each of the 8 SCIAMACHY channels agrees with the Kurucz data within ±2–3%. The presented analysis proves that SCIAMACHY is a valuable tool to monitor solar irradiance variations.     

Nonlinear processes in cosmic-ray precursor of strong supernova shock: Maximum energy and average energy spectrum of accelerated particles

The instability in the cosmic-ray precursor of a supernova shock is studied. The level of turbulence in this region determines the maximum energy of accelerated particles. The consideration is not limited by the case of weak turbulence. It is assumed that the Kolmogorov type nonlinear wave interactions together with the ion-neutral collisions restrict the amplitude of random magnetic field. As a result, the maximum energy of accelerated particles strongly depends on the age of a SNR. The average spectrum of cosmic rays injected in the interstellar medium in the course of adiabatic SNR evolution takes the approximate form E⁻² at energies larger than 10–30 GeV/nucleon with the maximum energy that is close to the position of the knee in cosmic-ray spectrum at 4 × 10¹⁵ eV. At an earlier stage of SNR evolution – the ejecta-dominated stage, the particles are accelerated to higher energies and have a rather steep power-law distribution. These results suggest that the knee may mark the transition from the ejecta-dominated to the adiabatic evolution of SNR shocks which accelerate cosmic rays.