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.     

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.     

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.     

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⁻¹.     

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

Giant radio pulses from the Crab pulsar

Individual giant radio pulses (GRPs) from the Crab pulsar last only a few microseconds. However, during that time they rank among the brightest objects in the radio sky reaching peak flux densities of up to 1500 Jy even at high radio frequencies. Our observations show that GRPs can be found in all phases of ordinary radio emission including the two high frequency components (HFCs) visible only between 5 and 9 GHz [Moffett, D.A., Hankins, T.H. Multifrequency radio observations of the Crab pulsar. Astrophys. J. 468, 779–783, 1996]. This leads us to believe that there is no difference in the emission mechanism of the main pulse (MP), inter pulse (IP) and HFCs. High resolution dynamic spectra from our recent observations of giant pulses with the Effelsberg telescope at a center frequency of 8.35 GHz show distinct spectral maxima within our observational bandwidth of 500 MHz for individual pulses. Their narrow band components appear to be brighter at higher frequencies (8.6 GHz) than at lower ones (8.1 GHz). Moreover, there is an evidence for spectral evolution within and between those structures. High frequency features occur earlier than low frequency ones. Strong plasma turbulence might be a feasible mechanism for the creation of the high energy densities of ∼6.7 × 10⁴ erg cm⁻³ and brightness temperatures of ∼10³¹ K.     

Features of the 11-year variation of galactic cosmic rays in different periods of solar magnetic cycles

Variations of galactic cosmic ray intensity have been studied based on the neutron monitors and interplanetary magnetic field experimental data for different ascending and descending epochs of solar activity. The dependence of the diffusion coefficient on the cosmic ray particles rigidity R is stronger in the maxima epoch than in the minima epoch of solar activity. For the period of 1977–1981 (qA > 0) the diffusion coefficient for the minimum epoch is, χ_min ∝ R^{0.7 ± 0.04} and for the maximum χ_max ∝ R^{1.3 ± 0.05}; for the period of 1987–1990 (qA < 0), χ_min ∝ R^{0.8 ± 0.05} and χ_max ∝ R^{1.1 ± 0.04}. The exponents ν_y and ν_z of the power spectral density of the B_y and B_z components of the IMF in the region of the frequencies (10⁻⁶– 4 × 10⁻⁶) Hz are larger for the minimum epoch of 1987 (ν_y ≈ 2.0 and ν_z ≈ 1.93) than for the maximum epoch of 1990 (ν_y ≈ 1.43 and ν_z ≈ 1.27).     

Energetic particle transport in anisotropic magnetic turbulence

We study energetic particle transport in a magnetic field configuration which models the solar wind magnetic turbulence plus the background field. A power-law Fourier amplitude is used for the fully 3D turbulence model, and in order to model anisotropic turbulence, the constant amplitude surfaces in k space are ellipsoids. The turbulence correlation lengths parallel (perpendicular) to the background magnetic field l_∥ (l_⊥) are varied in a wide range, and proton energies from 1 MeV to 10 GeV are assumed. Considering propagation on a distance corresponding to 1 AU, it is found that transport parallel and perpendicular to the background field heavily depends on the turbulence anisotropy, that is on the ratio l_∥/l_⊥. The spatial distribution of energetic particle follows the shape of magnetic flux tube up to about 10 MeV, while for larger energies the structure of the magnetic flux tube is progressively washed out. The scatterplots of particle distribution show intermittent, non Gaussian structures for l_∥ ≪ l_⊥ (quasi slab turbulence), while a more diffusive, Gaussian structure is obtained for l_∥ ≫ l_⊥ (quasi 2D turbulence). The long time behavior of transport shows that anomalous (subdiffusive perpendicular and superdiffusive parallel) transport regimes are obtained for l_∥ ≪ l_⊥, while Gaussian diffusive transport is obtained for both l_∥ ≫ l_⊥ and the isotropic turbulence case.     

Regional tropospheric responses to long-term solar activity variations

The influence of ∼200-year solar activity variations (de Vries cyclicity) on climatic parameters has been analyzed. Analysis of palaeoclimatic data from different regions of the Earth for the last millennium has shown that ∼200-year variations in solar activity give rise to a pronounced climatic response. Owing to a nonlinear character of the processes in the atmosphere–ocean system and the inertia of this system, the climatic response to the global influence of solar activity variations has been found to have a regional character. The regions where the climatic response to long-term solar activity variations is stable and the regions where the climatic response is unstable, both in time and space, have been revealed. It has also been found that a considerable lag of the climatic response and reversal of its sign with respect to the solar signal can occur. Comparison of the obtained results with the simulation predictions of the atmosphere–ocean system response to long-term solar irradiance variations (T > 40 years) has shown that there is a good agreement between experimental and simulation results.     

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.     

Characterization of the semi-annual-oscillation in mesospheric temperatures at low-latitudes

Novel measurements of the seasonal variability in mesospheric temperature at low-latitudes have been obtained from Maui, Hawaii (20.8°N, 156.2°W) during a 25-month period from October 2001 to January 2004. Independent observations of the OH (6, 2) Meinel band (peak height ∼87 km) and the O₂ (0–1) atmospheric band emission (∼94 km) were made using the CEDAR Mesospheric Temperature Mapper. The data revealed a coherent oscillation in emission intensity and rotational temperature with a well-defined periodicity of 181 ± 7 days. The amplitude of this oscillation was determined to be ∼5–6 K in temperature and ∼8–9% in intensity for both the OH and O₂ data sets. In addition, a strong asymmetry in the shape of the oscillation was also observed with the spring maximum significantly larger than the fall peak. These data provide new evidence in support of a semi-annual-oscillation in mesospheric temperature (and airglow emission intensities) and help quantify its seasonal characteristics.     

XMM-Newton observations of the Mouse,SLX 1744–299 and SLX 1744–300

We observed the radio and X-ray source G359.23–0.82, also known as “the Mouse”, with XMM-Newton. The X-ray image of this object shows a point-like source at the Mouse’s “head”, accompanied by a “tail” that extends for about 40″ westward. The morphology is consistent with that observed recently with Chandra [Gaensler, B.M., van der Swaluw, E., Camilo, F., et al. The Mouse that soared: high resolution X-ray imaging of the pulsar-powered bow shock G359.23–0.82, ApJ 616, 383–402, 2004]. The spectrum of the head can be described by a power-law model with a photon index Γ ≃ 1.9. These results confirm that the Mouse is a bow-shock pulsar wind nebula (PWN) powered by PSR J1747–2958. We found that the hydrogen column density toward the Mouse, N_H = (2.60 ± 0.09) × 10²² cm⁻², is 20%–40% lower than those toward two serendipitously detected X-ray bursters, SLX 1744–299 and SLX 1744–300. At a plausible distance of 5 kpc, the X-ray luminosity of the Mouse, L(0.5–10 keV) = 3.7 × 10³⁴ erg s⁻¹, is 1.5% of the pulsar’s spin-down luminosity. We detected a Type I X-ray burst from SLX 1744–300 and found a possible decrease of N_H and persistent luminosity for this source, in comparison with those observed with ROSAT in 1992.     

Propagation of secondary antiprotons and cosmic rays in the Galaxy

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., B/C) produce too few antiprotons. In the present paper, we discuss one possibility to overcome these difficulties. Using the measured antiproton flux and B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local “unprocessed” component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z ⩽ 28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three “radioactive clock” isotopes in CR, ¹⁰Be, ²⁶Al, ³⁶Cl, are all consistent and indicate a halo size z_h ∼ 4 kpc based on the most accurate data taken by the ACE spacecraft.     

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.     

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.     

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.     

Scaling and universality of the thermal conductivity of liquid 4He near the superfluid transition and in restricted geometries

We present measurements of the thermal conductivity λ(t, P, L) = l/ρ(t, P, L) near the superfluid transition of ⁴He at saturated vapor pressure and confined in cylindrical geometries with radii L = 0.5 and 1.0 μm (t ≡ T/T_λ(P) − 1). For L = 1.0 μm measurements at six pressures P are presented. At and above T_λ the data are consistent with a universal scaling function F(X) = (L/ξ_o)^x/ν(ρ/ρ₀), X = (L/ξ_o)^1/νt valid for all P (ρ₀ and x are the pressure-dependent amplitude and effective exponent of the bulk resistivity ρ(t, P, ∞) = ρ₀t^x and ξ = ξ₀t^−ν is the correlation length). Indications of breakdown of scaling and universality are observed below T_λ.     

A Gemini observation of the anomalous X-ray pulsar 1RXS J170849-400910

The anomalous X-ray pulsars (AXPs) represent a growing class of neutron stars discovered at X-ray energies. While the nature of their multi-wavelength emission mechanism is still under debate, evidence has been recently accumulating in favor of their magnetar nature. Their study in the optical and infrared (IR) wavelengths has recently opened a new window to constrain the proposed models. We here present a brief overview of AXPs and our Gemini-South observation of 1RXS J170849-400910, which is a relatively bright AXP discovered with ROSAT and later found to be an 11 s X-ray pulsar by ASCA. The observation was taken with the near-IR imager Flamingos in J (1.25 μm), H (1.65 μm), and K_s (2.15 μm). We confirm the recent detection by (ApJ, 589, L93–L96) of a source coincident with the CHANDRA source (candidate ‘A’). Our derived magnitudes of J = 20.6 (0.2), H = 18.6 (0.2), and K_s = 17.1 (0.2) are consistent with those derived by (ApJ, 589, L93–L96), and indicate that if this source is indeed the IR counterpart to 1RXS J170849-400910, then there is no evidence of variability from this AXP. However, given the lack of IR variability and the relatively high IR to X-ray flux of this source when compared to the other AXPs, we conclude that this source is unlikely the counterpart of the AXP, and that the other source (candidate ‘B’) within the CHANDRA error circle should not be ruled out as the counterpart. Further monitoring of these sources and a deep observation of this complex field are needed to confirm the nature of these sources and their association with the AXP.     

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.