Impacts on Hubble Space Telescope solar arrays: Discrimination between natural and man-made particles

A Post-Flight Investigation was initiated by the European Space Agency to analyze impacts on solar arrays of the Hubble Space Telescope (HST), exposed to space for 8.25 years at approximately 600 km altitude. The solar cells deployed during the first Service Mission (SM-1 in December 1993) were retrieved in March 2002 as part of Service Mission 3B (SM-3B). A sub-panel of 2 m² was cut from the arrays for subsequent selection and removal of individual solar cells for analysis. Six cells (4.8 × 10⁻³ m²) were surveyed for flux of all craters of sizes greater than 5 microns. Analytical scanning electron microscopy was used to analyse residues in 111 features of 3–4000 micron conchoidal detachment diameter (D_co), examined on 23 solar cells. Eighty three show identifiable residue: 38 are Space Debris impacts and 45 Micrometeoroid impacts. Of the remaining 28, 2 contain residue of ambiguous origin, 1 is probably a minor manufacturing flaw, 1 is obscured by contamination, and 24 are unresolved, lacking recognizable residue.The majority of space debris impacts on the SM-3B cells are less than 80 microns D_co, dominated by Al-rich residue, probably of solid rocket motor origin, although three may be due to sodium metal droplet impacts. Three larger features include paint pigment and binder, ferrous alloy, and possible carbon-fibre composite material debris.Micrometeoroid residues are found across the entire crater size range and dominate features of between 100 and 1000 microns, their residues are similar to those found in earlier SM-1 surveys. Fe- and Mg-rich silicates dominate; Fe sulphides are common and there are occasional vesicular Ni- and S-bearing mafic silicates of hydrous phyllosilicate origin. A single sodium aluminosilicate residue and one Fe Ni metal residue were found; as well as enigmatic Mg- and S-bearing residues, all considered as probably of micrometeoroid origin. A few Fe-, O- and C-bearing residues were classified as of ambiguous origin.     

Observed flux density enhancement at submillimeter wavelengths during an X-class flare

We analyse the 30 October, 2004, X1.2/SF solar event that occurred in AR 10691 (N13 W18) at around 11:44 UT. Observations at 212 and 405 GHz of the Solar Submillimeter Telescope (SST), with high time resolution (5 ms), show an intense impulsive burst followed by a long-lasting thermal phase. EUV images from the Extreme Ultraviolet Imaging Telescope (SOHO/EIT) are used to identify the possible emitting sources. Data from the Radio Solar Telescope Network (RSTN) complement our spectral observations below 15 GHz. During the impulsive phase the turnover frequency is above 15.4 GHz. The long-lasting phase is analysed in terms of thermal emission and compared with GOES observations. From the ratio between the two GOES soft X-ray bands, we derive the temperature and emission measure, which is used to estimate the free-free submillimeter flux density. Good temporal agreement is found between the estimated and observed profiles, however the former is larger than the latter.     

Prospects for observing dark-matter remnants with GLAST

The satellite-based experiment, GLAST (Gamma-ray Large Area Space Telescope), is under construction and is planned to measure the cosmic γ-ray flux in the energy range 20 MeV to >300 GeV, with supporting measurements for γ-ray bursts from 10 keV to 25 MeV. With its launch in 2007, GLAST will open a new and important window on a wide variety of high-energy phenomena, including exotic relics from the Big Bang. Among these may be the decay/annihilation products of the hypothesized super symmetric image of the known particles. Single-photon energy thresholds for channels leading to such final states have been excluded in a model-dependent manner by accelerator searches to energies greater than 50 GeV. The ability of GLAST to set limits on this important component of cosmological evolution is presented along with an update on the present status of this mission.     

Impact effects from small size meteoroids and space debris

When the impact risk from meteoroids and orbital debris is assessed the main concern is usually structural damage. With their high impact velocities of typically 10–20 km/s millimeter or centimeter sized objects can puncture pressure vessels and other walls or lead to destruction of complete subsystems or even whole spacecraft. Fortunately chances of collisions with such larger objects are small (at least at present). However, particles in the size range 1–100 μm are far more abundant than larger objects and every orbiting spacecraft will encounter them with certainty. Every solar cell (8 cm² area) of the Hubble Space Telescope encountered on average 12 impacts during its 8.25 years of space exposure. Most were from micron sized particles.     

On long-term changes of electron density in the upper mesosphere

Recent review study done jointly by 19 experts of 17 institutes shows zero trend of temperature in the upper mesosphere. In the light of this latest development, we have examined the long-term changes in electron density, [e], in this region. The study has been concentrated at 80 km. At this altitude, electrons are mainly produced by the interaction of nitric oxide, NO, by solar Ly-α. Any long-term change in this flux will affect trend of [e]. Considering this flux proportional to 10.7 cm solar flux, analysis of available 10.7 cm solar flux data from 1948 to 2003 has been made. A decreasing trend up to about 1970 and then an increasing trend are found. The over-all increasing trend of Ly-α flux during the past five decades is ∼0.17% per year. This increase also gives a ∼0.17% increasing trend per year in [e]. This non-anthropogenic increase is much less compared to the observed increase in [e] which is reported to be >0.7% per year. The observed increase in [e] of this magnitude will then, predominantly, be due to the anthropogenic effect. In zero trend in temperature, significant change in electron loss coefficient, α_eff, and [NO] are unlikely to take place to cause a significant change in [e]. The increase in [e] > 0.7% per year then can be explained by considering a decreasing trend in [O₂].     

On the 27-day variations of the galactic cosmic ray anisotropy and intensity for different periods of solar magnetic cycle

The average amplitude of the 27-day variation of the galactic cosmic ray anisotropy calculated based on the neutron monitors experimental data is larger in the qA > 0 period than in the qA < 0 period of solar magnetic cycle. The amplitudes of the 27-day variation of the galactic cosmic rays anisotropy do not depend on the tilt angles of the heliospheric neutral sheet for different the qA > 0 and the qA < 0 periods of solar magnetic cycle. A good correlation has been revealed between the changes of the amplitudes of the 27-day variations of the galactic cosmic ray anisotropy and intensity versus the qA > 0 and the qA < 0 periods of solar magnetic cycle.     

Observations of “zebra” pattern in cm-range with spatial resolution

We present the results of the first observations of the solar microwave burst with fine spectral structure of zebra type at the frequency about 5.7 GHz. The burst has been detected simultaneously by the Siberian Solar Radio Telescope and by the spectropolarimeter of the National Astronomical Observatory of China. Zebra pattern consisted of three parallel stripes with complex frequency drift. The degree of circular polarization of emission reached 100%, the polarization sense corresponded to the extraordinary wave (X-mode). We have determined the plasma parameters in the emission source: plasma density about 10¹¹ cm⁻³, magnetic field strength 60–80 G. We argue that in the given event the most probable mechanism of the zebra pattern generation is non-linear coupling of harmonics of Bernstein modes.     

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.     

On solar 5.7 GHz subsecond burst source sizes

Due to the narrow bandwidth and the small size of the subsecond microwave pulses, we can use them as probe sources for study of propagation effects in the low corona. More than 160 microwave bursts with subsecond pulses (SSP) have been observed with the Siberian Solar Radio Telescope at 5.7 CHz for the period 2000–2004. Working with a large dataset of homogeneous observational material (spatial resolution from 15″ to 20″, temporal resolution 14 mc), we estimated sizes of SSP and studied relation between SSP sizes and sense of the polarization and their position on the Sun. Our results are in accordance with those obtained during the 22nd solar cycle. The apparent sizes of SSP increase toward the solar limb. The obtained dependence is in agreement with Bastian‘s (Bastian, T.S. Angular scattering of solar radio emission by coronal turbulence. ApJ 426, 774, 1994.) model calculations. The center-to-limb variation of the source size is explained by scattering on plasma turbulence along the ray path in the solar corona. The most events with high polarization occur near the central meridian (±30°). The polarization sense corresponds mainly to the ordinary mode.     

The temperature map of the loop and coronal sources of an occulted solar flare

An occulted solar flare occurred at about 06:07 UT on 2002, November 2. The RHESSI X-ray images show two separate parts. The lower part consists of a complete loop and the upper part a coronal source which well extends above the solar limb. The loop source shrank for about 3 min with a speed of ∼24 km s⁻¹ during the early impulsive phase and then expanded at ∼7 km s⁻¹, while the coronal source presented an upward motion at about 6 km s⁻¹. We obtained the temperature map of the loop source from RHESSI image spectrum. The temperature of the loop increases with altitude, indicating that the reconnection X-point of this flare is located above the loop source. However, the apparent coronal source is the top of another independent large-scale loop.     

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.     

Evaluation of solar proton spectra using balloon cosmic ray observations and Monte Carlo simulation results

We have developed a method to evaluate the spectrum of solar energetic protons at the top of the Earth’s atmosphere from the measurements of our balloon cosmic ray experiment. By using the Monte Carlo PLANETOCOSMICS code based on Geant4 we compute the interaction of solar protons [10 MeV–10 GeV] with the Earth’s atmosphere. We obtain the angular and energy distributions of secondary particles (p, e⁻, e⁺, photons, muons) at different atmospheric levels as a function of primary proton spectra. By comparing the calculated depth dependence of the particle flux with the data obtained by our balloon experiment we can deduce the parameters of the solar proton spectrum that best fit the observations. In this paper we discuss our solar proton spectrum estimation method, and present results of its application to selected solar proton events from 2001 to 2005.     

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.     

Modeling of effect of polarization on UV sky radiance during twilight

The spatial distribution of the vector of the Stokes parameters characterizing the radiance intensity and the radiance polarization describes the radiation field in the atmosphere. A simplified treatment of light as the scalar has only restricted application. A few studies compared previously results of the vector and scalar radiative transfer models and showed that scalar models are in error by up to 10% for many cases. Though several observational conditions were exploited, an effect of polarization on modeling of UV radiance has not been investigated yet for twilight. The paper presents a preliminary study of modeled UV radiance during twilight taking into account polarization. The intensity and the degree of linear polarization of the scattered UV radiance for two cases of the ground-based observations are discussed. In the first case, radiation incoming from the zenith for the solar zenith angles (SZA) from 90° to 98° is under investigation. Radiation in the solar principal plane for the beginning of twilight (SZA = 90.1°) was calculated in the second case. The study showed that the UV radiation field in the twilight atmosphere can be handled correctly only using the vector theory. The errors of scalar radiative transfer strongly depend on wavelength, line of an observation and solar position. The revealed distortion of the zenith radiance caused by using of the scalar approximation reaches maximum of 15% at 340 nm for the solar zenith angle (SZA) equal to 98°. The shorter wavelengths have the smaller errors, about 5% at 305 nm for SZA = 98°, due to the larger part of the single scattered radiance. The error of the scalar modeling may be as large as −17% for radiance incoming from the horizon for SZA = 90.1°. Scalar radiative transfer models underestimate the integral intensity in the principal plane up to 3–4% ± 0.5% at SZA = 90.1° for wavelengths from 320 to 340 nm. This should be taken into account in problems of radiative budget estimation and remote sensing of the atmosphere exploiting the twilight period.     

The total ozone and UV solar radiation over Stara Zagora,Bulgaria

The results from direct ground-based solar UV irradiance measurements and the total ozone content (TOC) over Stara Zagora (42° 25′N, 25° 37′E), Bulgaria are presented. During the period 1999–2003 the TOC data show seasonal variations, typical for the middle latitudes – maximum in the spring and minimum in the autumn. The comparison between TOC ground-based data and Global Ozone Monitoring Experiment (GOME) satellite-borne ones shows a seasonal dependence of the differences between them.A strong negative relationship between the total ozone and the 305 nm wavelength irradiance was found. The dependence between the two variables is significant (r = −0.62 ± 0.18) at 98% confidence level.The direct sun UV doses for some specific biological effects (erythema and eyes) are obtained. The estimation of the radiation amplification factor RAF shows that the ozone reduction by 1% increases the erythemal dose by 2.3%. The eye-damaging doses are more influenced by the TOC changes and in this case RAF = −2.7%.The amount of these biological doses depended on the solar altitude over the horizon. This dependence was not so strong when the total ozone content in the atmosphere was lower.     

Reconstruction of solar irradiance using the Group sunspot number

We present a reconstruction of total solar irradiance since 1610 to the present based on variations of the surface distribution of the solar magnetic field. The latter is calculated from the historical record of the Group sunspot number using a simple but consistent physical model. Our model successfully reproduces three independent data sets: total solar irradiance measurements available since 1978, total photospheric magnetic flux from 1974 and the open magnetic flux since 1868 (as empirically reconstructed from the geomagnetic aa-index). The model predicts an increase in the total solar irradiance since the Maunder Minimum of about 1.3 Wm⁻².     

A spectroscopic measurement of high velocity spray plasma from an M-class flare and coronal mass ejection

The M1.5-class flare and associated coronal mass ejection (CME) of 16 February 2011 was observed with the Extreme ultraviolet Imaging Spectrometer on board the Hinode spacecraft. Spray plasma associated with the CME is found to exhibit a Doppler blue-shift of 850 km s^?1 – one of the largest values reported from spectroscopy of the solar disk and inner corona. The observation is unusual in that the emission line (Fe xii 193.51 Å) is not observed directly, but the Doppler shift is so large that the blue-shifted component appears in a wavelength window at 192.82 Å, intended to observe lines of O v, Fe xi and Ca xvii. The Fe xii 195.12 Å emission line is used as a proxy for the rest component of 193.51 Å. The observation highlights the risks of using narrow wavelength windows for spectrometer observations when observing highly-dynamic solar phenomena. The consequences of large Doppler shifts for ultraviolet solar spectrometers, including the upcoming Multi-slit Solar Explorer (MUSE) mission, are discussed.     

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.     

Supernova remnant 1987A: The latest report from the Chandra X-ray Observatory

We continue monitoring supernova remnant (SNR) 1987A with the Chandra X-ray Observatory. As of 2004 January, bright X-ray spots in the northwest and the southwest are now evident in addition to the bright eastern ring. The overall X-ray spectrum, since 2002 December, can be described by a planar shock with an electron temperature of ∼2.1 keV. The soft X-ray flux is now 8 × 10⁻¹³ ergs cm⁻² s⁻¹, which is about five times higher than four years ago. This flux increase rate is consistent with our prediction based on an exponential density distribution along the radius of the SNR between the HII region and the inner ring. We still have no direct evidence of a central point source, and place an upper limit of L_X = 1.3 × 10³⁴ ergs s⁻¹ on the 3–10 keV band X-ray luminosity.     

Gamma-ray observations of explosive nucleosynthesis products

In this review I discuss the various γ-ray emission lines that can be expected and, in some cases have been observed, from radioactive explosive nucleosynthesis products. The most important γ-ray lines result from the decay chains of ⁵⁶Ni, ⁵⁷Ni, and ⁴⁴Ti. ⁵⁶Ni is the prime explosive nucleosynthesis product of Type Ia supernovae, and its decay determines to a large extent the Type Ia light curves. ⁵⁶Ni is also a product of core-collapse supernovae, and in fact, γ-ray line emission from its daughter product, ⁵⁶Co, has been detected from SN1987A by several instruments. The early occurrence of this emission was surprising and indicates that some fraction of ⁵⁶Ni, which is synthesized in the innermost supernova layers, must have mixed with the outermost supernova ejecta.Special attention is given to the γ-ray line emission of the decay chain of ⁴⁴Ti (⁴⁴Ti → ⁴⁴Sc → ⁴⁴Ca), which is accompanied by line emission at 68, 78, and 1157 keV. As the decay time of ⁴⁴Ti is ∼86 yr, one expects this line emission from young supernova remnants. Although the ⁴⁴Ti yield (typically 10⁻⁵–10⁻⁴M_⊙) is not very high, its production is very sensitive to the energetics and asymmetries of the supernova explosion, and to the mass cut, which defines the mass of the stellar remnant. This makes ⁴⁴Ti an ideal tool to study the inner layers of the supernova explosion. This is of particular interest in light of observational evidence for asymmetric supernova explosions.The γ-ray line emission from ⁴⁴Ti has so far only been detected from the supernova remnant Cas A. I discuss these detections, which were made by COMPTEL (the 1157 keV line) and BeppoSAX (the 68 and 78 keV lines), which, combined, give a flux of (2.6 ± 0.4 ± 0.5) × 10⁻⁵ ph cm⁻² s⁻¹ per line, suggesting a ⁴⁴Ti yield of (1.5 ± 1.0) × 10⁻⁴M_⊙. Moreover, I present some preliminary results of Cas A observations by INTEGRAL, which so far has yielded a 3σ detection of the 68 keV line with the ISGRI instrument with a flux that is consistent with the BeppoSAX detections. Future observations by INTEGRAL-ISGRI will be able to constrain the continuum flux above 90 keV, as the uncertainty about the continuum shape, is the main source of systematic error for the 68 and 78 keV line flux measurements. Moreover, with the INTEGRAL-SPI instrument it will be possible to measure or constrain the line broadening of the 1157 keV line. A preliminary analysis of the available data indicates that narrow line emission (i.e., Δv < 1000 km s⁻¹) can be almost excluded at the 2σ level, for an assumed line flux of 1.9 × 10⁻⁵ ph cm⁻² s⁻¹.