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.     

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.     

Estimates of mesospheric gravity wave activity over convection from a global model

A new convective gravity wave source spectrum parameterization has been implemented in the Whole Atmosphere Community Climate Model version 2 (WACCM2). This parameterization specifies the momentum flux phase speed spectrum of gravity waves in the Tropics based on the properties of underlying convection; Hence, this parameterization provides realistic global estimates of gravity wave activity. In this paper, we show the estimated gravity wave phase speed spectra in the Tropics from a WACCM2 simulation, at the source level and at 85 km. Spatial distribution of gravity wave activity at 85 km is also presented. Subsequently, we discuss the factors that are primarily responsible for the estimated differences in gravity wave distribution across phase speeds with latitude and asymmetries in direction of gravity wave propagation in the mesosphere. We also examine which of the model assumptions can lead to uncertainties in our estimates of mesospheric gravity wave activity and we discuss how these assumptions provide challenges for comparison with observations of gravity waves in the mesosphere.     

Centre-to-limb variation of photospheric facular radiance and image resolution

We study the effect of the angular resolution on the determination of the angular properties of the facular radiance. We analyze photospheric intensity in the continuum, around the Ni 676.8 nm line, and longitudinal magnetic field along the line of sight, measured by the MDI instrument aboard SOHO with two spatial resolutions, 4″ and 1.2″ (2″ and 0.6″ pixels, respectively). The effect of the limited photometric sensitivity of the instrument and the limited information on the angular structure of the magnetic field tubes are considered. Our study of the high-resolution data shows that intensity contrast of magnetic features between 80 and 600 Gauss increases from centre to limb up to a maximum that occurs at higher heliocentric angles (θ) when obtained with higher resolution data than for lower resolution data. There is a suggestion that at heliocentric angles below about 75° there is only a monotonic increase in the contrast as one goes from cos (θ) = 1 to cos (θ) = 0.2.     

A radio study of the mouse,G359.23 − 0.82

The recent detection of a young pulsar powering “the Mouse”, G359.23 − 0.82, as well as detailed imaging of surrounding nebular X-ray emission, have motivated us to investigate the structural details and polarization characteristics of the radio emission from this axisymmetric source with a supersonic bow shock. Using polarization data at 3.6 and 6 cm, we find that the magnetic field wraps around the bow-shock structure near the apex of the system, but downnstream runs parallel to the inferred direction of the pulsar’s motion. The rotation measure (RM) distribution of the Mouse also suggests that the low degree of polarization combined with a high RM ahead of the pulsar result from internal plasma within the bow-shock region. In addition, using sub-arcsecond radio image of the Mouse, we identify modulations in the brightness distribution of the Mouse that may be associated with the unshocked pulsar wind behind the pulsar. Lastly, we discuss the relationship between the Mouse and its neighboring shell-type supernova remnant G359.1 − 0.5 and argue that these two sources could potentially have the same origin.     

Reconstruction of solar UV irradiance

Understanding solar influence on the Earth’s climate requires a reconstruction of solar irradiance for the pre-satellite period. Considerable advances have been made in modelling the irradiance variations at wavelengths longer than 200 nm. At shorter wavelengths, however, the LTE approximation usually taken in such models fails, which makes a reconstruction of the solar UV irradiance a rather intricate problem. We choose an alternative approach and use the observed SUSIM UV spectra to extrapolate available models to shorter wavelengths.     

The contribution of radar interferometry to the assessment of landslide hazards

Interferometric analysis of repeat pass radar images (InSAR) enables the mapping of slope movements with high accuracy, being of great relevance for the assessment of landslide risk. The InSAR method is briefly explained, with emphasis on data processing for landslide motion mapping. Constraints of the method are discussed, mainly resulting from temporal decorrelation of the radar signal in vegetated areas. Application examples for InSAR mapping of very slow slope movements, using SAR data of the European Earth Observation Satellite ERS, are shown for a mass waste above an alpine traffic route and for slope deformation in a mountain village. The comparison with GPS measurements of surface motion shows RMSE = 3.7 mm a⁻¹ between the InSAR and GPS retrieved motion, where part of this difference is attributed to the different observation periods. The investigations demonstrate the good capability of InSAR for detecting and monitoring surface deformation on mountain slopes, thus providing valuable support for landslide inventories and risk assessment.     

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.     

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.     

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.     

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 temporal variation of mean gravity of atmosphere

One of the precise widely used global Zenith Hydrostatic Delay (ZHD) model is based on the gravity value at the centroid of the atmospheric column at the station of observation and gravity value at the centroid is constant in this model for a specific location throughout the year. However, as the content and extent of atmosphere varies temporally, its centroid and consequently gravity value at the centroid also varies. Apart from this, the actual atmospheric condition of different region is not alike. Therefore, there is a need to develop a regional mean gravity model and development of such model has been discussed in this paper. To obtain the mean gravity model, first a regional model of centroid height of atmospheric column was developed as a function of the surface pressure and temperature. It was developed by multiple regressions between estimated centroid of the atmosphere and surface pressure, surface temperature using radiosonde data of five radiosonde stations spread over the Indian subcontinent. Three years radiosonde data from 2006 to 2008 was used for each station. The root mean square error in estimating centroid of the atmospheric column is about ±326 m, which is negligible considering the variability of the atmosphere and its content. The centroid height model has been used to formulate the mean gravity model, considering uniform lapse rate in gravity with height. It is found that proposed mean gravity model provides temporal variation of mean gravity values at the centroid and thus matches with the reality. The interesting advantage of the developed model is that the model shows diurnal variation of mean gravity. The accuracy of ZHD has shown of the order of about 0.3 mm using the developed regional mean gravity model. However, already developed ZHD model has shown a slight inferior result compared to the developed model. These models have shown accuracy of about 0.8 mm and 0.6 mm.     

Active galactic nuclei imaging programs of the RadioAstron mission

Imaging relativistic jets in active galactic nuclei (AGN) at angular resolution significantly surpassing that of the ground-based VLBI at centimetre wavelengths is one of the key science objectives of the RadioAstron space-VLBI mission. There are three RadioAstron imaging key science programs that target both nearby radio galaxies and blazars, with one of the programs specifically focusing on polarimetry of the jets. The first images from these programs reach angular resolution of a few tens of microarcseconds and reveal unprecedented details about the jet collimation profile, magnetic field configuration, and Kelvin-Helmholtz instabilities along the flow in some of the most studied AGN (3C 84, BL Lac, 3C 273, S5 0836 + 710). Here we give an overview of the goals and strategy of these three ongoing programs, highlight their early results, and discuss the challenges of space-VLBI imaging.     

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.     

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.     

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.     

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.     

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

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.