Possible consequences of a disk around B0656+14 on e–e and near-1 Hz gravitational wave production

Nearby pulsars B0656+14 and Geminga were proposed in the literature as the main sources of cosmic-ray positrons observed near Earth above 10 GeV. B0656+14 has comparable distance from Earth, similar magnetic field and period of Geminga. However, observations in the R and I bands indicate the presence of a disk of approximately 10⁻⁴ M_⊙ around B0656+14. Radio and pulsed γ-ray flux observations from this pulsar are also consistent with supernova fallback material and disk entering the light cylinder and partially quenching the development of electromagnetic showers in the magnetosphere. If this is the case, B0656+14 has unlikely given any contribution to e⁺ and e⁻ observed near Earth. Absolute flux measurements and the level of anisotropy in the high energy electron and positron arrival directions above 50 GeV will help in revealing if none, one of both nearby pulsars are sources of these particles observed near Earth.     

Interplanetary spacecraft navigation using pulsars

We demonstrate how observations of pulsars can be used to help navigate a spacecraft travelling in the solar system. We make use of archival observations of millisecond pulsars from the Parkes radio telescope in order to demonstrate the effectiveness of the method and highlight issues, such as pulsar spin irregularities, which need to be accounted for. We show that observations of four millisecond pulsars every seven days using a realistic X-ray telescope on the spacecraft throughout a journey from Earth to Mars can lead to position determinations better than ∼20 km and velocity measurements with a precision of ∼0.1 ms⁻¹.     

Effects of solar modulation on the cosmic ray positron fraction

We implemented a 2D Monte Carlo model to simulate the solar modulation of galactic cosmic rays. The model is based on the Parker’s transport equation which contains diffusion, convection, particle drift and energy loss. Following the evolution in time of the solar activity, we are able to modulate a local interstellar spectrum (LIS), that we assumed isotropic beyond the termination shock, down to the Earth position inside the heliosphere. In this work we focused our attention to the cosmic ray positron fraction at energy below ∼10 GeV, showing how the particle drift processes could explain different results for AMS-01 and PAMELA. We compare our modulated spectra with observations at Earth, and then make a prediction of the cosmic ray positron fraction for the AMS-02 experiment.     

Where does the heliospheric modulation of galactic cosmic rays start?

The long outstanding question of where the heliospheric (solar) modulation of galactic cosmic rays actually begins, in terms of spatial position, as well as at what high kinetic energy, can now be answered. Both answers are possible by using the results of an advanced numerical model, together with appropriate observations. Voyager 1 has been exploring the outskirts of the heliosphere and is presently entering what can be called the very local interstellar medium. It has been generally expected, and accepted, that once the heliopause is crossed, the local interstellar spectrum (LIS) should be measured in situ by the Voyager spacecraft. However, we show that this may not be the case and that modulation effects on galactic cosmic rays can persist well beyond the heliopause. For example, proton observations at 100 MeV close to the heliopause can be lower by ∼$\sim$25% to 40% than the LIS, depending on solar modulation conditions. It is also illustrated quantitatively that significant solar modulation diminishes above ∼$\sim$50 GeV at Earth. It is found that cosmic ray observations above this energy contain less that 5%$5\%$ solar modulation effects and should therefore reflect the LIS for galactic cosmic rays. Input spectra, in other words the very LIS, for solar modulation models are now constrained by in situ observations and can therefore not any longer be treated arbitrarily. It is also possible for the first time to determine the lower limit of the very LIS from a few MeV/nuc to very high energies.     

Modelling and observing Jovian electron propagation times in the inner heliosphere

The propagation of Jovian electrons in interplanetary space was modelled by solving the relevant transport equation numerically through the use of stochastic differential equations. This approach allows us to calculate, for the first time, the propagation time of Jovian electrons from the Jovian magnetosphere to Earth. Using observed quiet-time increases of electron intensities at Earth, we also derive values for this quantity. Comparing the modelled and observed propagation times we can gauge the magnitude of the transport parameters sufficiently to place a limit on the 6 MeV Jovian electron flux reaching Earth. We also investigate how the modelled propagation time, and corresponding Jovian electron flux, varies with the well-known ∼13 month periodicity in the magnetic connectivity of Earth and Jupiter. The results show that the Jovian electron intensity varies by a factor of ∼10 during this cycle of magnetic connectivity.     

北斗卫星导航系统X射线脉冲星可见性分析

X射线脉冲星可以为深空乃至近地航天器提供全源、自主、高精度的导航服务.对近地航天器而言，由于中心天体的遮挡，脉冲星信号并非全程可见，所以X射线脉冲星的可见性成为近地导航应用中挑选脉冲星的关键因素，也是脉冲星导航的首要工作.本文对挑选出的30颗X射线脉冲星在北斗系统7个轨道面上的可见性进行了仿真分析，得出了一天内可见脉冲星数量的变化情况以及一年内太阳、地球和月球遮挡的仿真结果,总结了影响脉冲星可见性的天体的遮挡规律，为利用脉冲星对北斗卫星进行自主定轨和守时提供了数据支持.     

Intermediate-term periodicities in relativistic solar electron fluences during solar cycles 22 and 23

In this paper, we have investigated the intermediate-term periodicities of the relativistic (E > 10 MeV) solar electron flares measured by IMP-8 satellite of NASA for the time period of 1986–2001. This period of investigation includes the entire solar cycle 22; ascending, maximum and a part of descending phase of the current solar cycle 23. To determine accurately the occurrence rate of electron flux, we have employed three different spectral decomposition techniques, viz. fast Fourier transformation (FFT); maximum entropy method (MEM) and Lomb–Scargle periodogram analysis method. For solar cycle 22, in the low frequency range, power spectrum analysis exhibits statistically significant periodicities at ∼706, ∼504 and ∼392 days. In the intermediate frequency range, we have found a series of significant periodicities ∼294, ∼221, ∼153, ∼86, ∼73 and ∼66 days. For short term, periodicities of ∼21–23, ∼31 and ∼37 days were found in power spectrum. When solar cycle 23 is considered the significant periodicities are ∼20, ∼23, ∼29, ∼39, ∼54, ∼63, ∼118, ∼133 and ∼154 days. These results provide evidence that the best known Rieger period (∼153 days), appeared in the high energetic electron flux data for cycle 22 and also likely during maxima of cycle 23. The existence of these periodicities has been discussed in the light of earlier results.     

Response of the mid-latitude D-region ionosphere to the total solar eclipse of 22 July 2009 studied using VLF signals in South Korean peninsula

In this paper, we analyze VLF signals received at Busan to study the the D-region changes linked with the solar eclipse event of 22 July 2009 for very short (∼390 km) transmitter–receiver great circle path (TRGCP) during local noon time 00:36–03:13 UT (09:36–12:13 KST). The eclipse crossed south of Busan with a maximum obscuration of ∼84%. Observations clearly show a reduction of ∼6.2 dB in the VLF signal strength at the time of maximum solar obscuration (84% at 01:53 UT) as compared to those observed on the control days. Estimated values of change in Wait ionospheric parameters: reflection height (h′) in km and inverse scale height parameter (β) in km⁻¹ from Long Wave Propagation Capability (LWPC) model during the maximum eclipse phase as compared to unperturbed ionosphere are 7 km and 0.055 km⁻¹, respectively. Moreover, the D-region electron density estimated from model computation shows 95% depletion in electron density at the height of ∼71 km. The reflection height is found to increase by ∼7 km in the D-region during the eclipse as compared to those on the control days, implying a depletion in the Lyman-α flux by a factor of ∼7. The present observations are discussed in the light of current understanding on the solar eclipse induced D-region dynamics.     

Detection of the high energy component of Jovian electrons in Low Earth Orbit with the PAMELA experiment

The PAMELA experiment is devoted to the study of cosmic rays in Low Earth Orbit with an apparatus optimized to perform a precise determination of the galactic antimatter component of c.r. It is constituted by a number of detectors built around a permanent magnet spectrometer. PAMELA was launched in space on June 15th 2006 on board the Russian Resurs-DK1 satellite for a mission duration of 3 years. The characteristics of the detectors, the long lifetime and the orbit of the satellite, will allow to address several aspects of cosmic-ray physics. In this work we discuss the observational capabilities of PAMELA to detect the electron component above 50 MeV. The magnetic spectrometer allows a detailed measurement of the energy spectrum of electrons of galactic and Jovian origin. Long term measurements and correlations with Earth–Jupiter 13 months synodic period will allow to separate these two contributions and to measure the primary electron Jovian component, dominant in the 50–70 MeV energy range. With this technique it will also be possible to study the contribution to the electron spectrum of Jovian e⁻ reaccelerated up to 2 GeV at the Solar Wind Termination Shock.     

Prospects of Global Navigation Satellite System (GNSS) reflectometry for geodynamic studies

Innovative processing of satellite radar altimetry over solid Earth has been successfully applied for observing geodynamic process of glacial isostatic adjustment over the former Laurentide Ice Sheet in the present-day Hudson Bay land region. In this contribution, a simulation is conducted to study the prospects of the applications of space-/airborne and land-based Global Navigation Satellite System (GNSS) reflectometry to synoptically observe global-scale geodynamic processes with a vertical accuracy of ∼2 mm/yr.     

Stochastic particle acceleration by the forced interaction of relativistic current sheets

We address the problem of interacting relativistic current sheets in self-consistent kinetic plasma simulations within the framework of the Particle-In-Cell model. The interaction is enforced in head-on collisions of up to 10 current sheets at relativistic bulk speeds. The simulations are motivated by the general problem of Poynting flux dissipation in ‘striped wind’ configurations presumably governing the relativistic outflows pervasive in pulsar winds and gamma-ray bursts. We identify the generation of non-thermal particles and formation of a stable power-law shape in the particle energy distributions f(γ) dγ ∝ γ^−s dγ. In 1D, a spectral index s ∼ 2 is observed and attributed to a stochastic Fermi-type acceleration mechanism. In 2D, the generic index of s ∼ 3–4 is retained as in previous simulations of individual current sheets. Whereas in 2D the high energy cut-off is constrained by the limited dissipation of magnetic energy, in 1D the process converts the bulk motion of current sheets towards directed particle momentum of an exclusive class of non-thermal particles.     

Modulation of galactic cosmic rays in a north–south asymmetrical heliosphere

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock. Voyager 1 crossed this termination shock at ∼94 AU in 2004, while Voyager 2 crossed it in 2007 at a different heliolatitude, about 10 AU closer to the Sun. These different positions of the termination shock confirm the dynamic and cyclic nature of the shock’s position. Observations from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that apart from the dynamic nature caused by changing solar activity there also may exist a global asymmetry in the north–south (polar) dimensions of the heliosphere, in addition to the expected nose–tail asymmetry. This relates to the direction in which the heliosphere is moving in interstellar space and its orientation with respect to the interstellar magnetic field. In this paper we focus on illustrating the effects of this north–south asymmetry on modulation of galactic cosmic ray Carbon, between polar angles of 55° and 125°, using a numerical model which includes all four major modulation processes, the termination shock and the heliosheath. This asymmetry is incorporated in the model by assuming a significant dependence on heliolatitude of the thickness of the heliosheath. When comparing the computed spectra between the two polar angles, we find that at energies E < ∼1.0 GeV the effects of the assumed asymmetry on the modulated spectra are insignificant up to 60 AU from the Sun but become increasingly more significant with larger radial distances to reach a maximum inside the heliosheath. In contrast, with E > ∼1.0 GeV, these effects remain insignificant throughout the heliosphere even very close to the heliopause. Furthermore, we find that a higher local interstellar spectrum for Carbon enhances the effects of asymmetric modulation between the two polar angles at lower energies (E < ∼300 MeV). In conclusion, it is found that north–south asymmetrical effects on the modulation of cosmic ray Carbon depend strongly on the extent of the geometrical asymmetry of the heliosheath together with the assumed value of the local interstellar spectrum.     

Standing Alfvén waves with m ≫ 1 in a dipole magnetosphere with moving plasma and aurorae

The structure of standing Alfvén waves with large azimuthal wave numbers (m ? 1) is studied in a dipole model of the magnetosphere with rotating plasma. In the direction across magnetic shells the structure of such waves is determined by their dispersion associated with curvature of geomagnetic field lines and corresponds to the travelling wave localized between toroidal and poloidal resonant surfaces. In projection into the ionosphere (along geomagnetic field lines) this structure is similar to the structure of a discrete auroral arc. The azimuthal structure of an auroral arc is similar to azimuthal structure of Alfvén waves with m ∼ 100. Possible interaction mechanisms between the Alfvén waves and energetic electron fluxes forming auroral arcs are discussed.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

Studies of electron and ion temperatures at 500 km altitude during sunrise using Indian SROSS C2 satellite

Solar dependence of electron and ion temperatures (T_e and T_i) in the ionosphere is studied using RPA data onboard SROSS C2 at an altitude of ∼500 km and 77°E longitude during early morning hours (04:00–07:00 LT) for three solar activities: solar minimum, moderate and maximum during winter, summer and equinox months in 10°S–20°N geomagnetic latitude. In winter the morning overshoot phenomenon is observed around 06:00 LT (T_e enhances to ∼4000 K) during low-solar activity and to T_e ∼ 3800 K, during higher solar activity. In summer, it is observed around 05:30 LT, but the rate of T_e enhancement is higher during moderate solar activity (∼2700 K/hr) than the low-solar activity (∼1700 K/hr). During equinox, this phenomenon is delayed and is observed around 06:00 LT (∼4200 K) during all three activities.     

An alternative Shell inversion technique - Analysis and validation based on COSMIC and ionosonde data

Multi-channel Global Positioning System (GPS) carrier phase signals, received by the six low Earth orbiting (LEO) satellites from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) program, were used to undertake active limb sounding of the Earth’s atmosphere and ionosphere via radio occultation. In the ionospheric radio occultation (IRO) data processing, the standard Shell inversion technique (SIT), transformed from the traditional Abel inversion technique (AIT), is widely used, and can retrieve good electron density profiles. In this paper, an alternative SIT method is proposed. The comparison between different inversion techniques will be discussed, taking advantage of the availability of COSMIC datasets. Moreover, the occultation results obtained from the SIT and alternative SIT at 500 km and 800 km, are compared with ionosonde measurements. The electron densities from the alternative SIT show excellent consistency to those from the SIT, with strong correlations over 0.996 and 0.999 at altitudes of 500 km and 800 km, respectively, and the peak electron densities (N_mF₂) from the alternative SIT are equivalent to the SIT, with 0.839 vs. 0.844, and 0.907 vs. 0.909 correlation coefficients when comparing to those by the ionosondes. These results show that: (1) the N_mF₂ and h_mF₂ retrieved from the SIT and alternative SIT are highly consistent, and in a good agreement with those measured by ionosondes, (2) no matter which inversion technique is used, the occultation results at the higher orbits (∼800 km) are better than those at the lower orbits (∼500 km).     

Jovian periodicities (∼10 h, ∼40 min) on Ulysses’ Distant Jupiter Encounter observations around the Halloween CIR events

We analyzed data from four different instruments (HI-SCALE, URAP, SWOOPS, VHM/FGM) onboard Ulysses spacecraft (s/c) and we searched for possible evidence of Jovian emissions when the s/c approached Jupiter during the times of Halloween events (closest time approach/position to Jupiter: February 5, 2004/R = 1683 R_J,θ = ∼49°). In particular, we analyzed extensively the low energy ion measurements obtained by the HI-SCALE experiment in order to examine whether low energy ion/electron emissions show a symmetry, and whether they are observed at north high latitudes upstream from the jovian bow shock, as is known to occur in the region upstream from the south bow shock as well ( Marhavilas et al., 2001). We studied the period from October 2003 to March 2004, as Ulysses moved at distances 0.8–1.2 AU from the planet at north Jovicentric latitudes <75°, and we present here an example of characteristic Jovian periodicities in the measurements around a CIR observed by Ulysses on days ∼348–349/2003 (R = 1894 R_J,θ = 72°). We show that Ulysses observed low energy ion (∼0.055–4.7 MeV) and electron (>∼40 keV) flux and/or spectral modulation with the Jupiter rotation period (∼10 h) as well as variations with the same period in solar wind parameters, radio and magnetic field directional data. In addition, characteristic strong ∼40 min periodic variations were found superimposed on the ∼10 h ion spectral modulation. Both the ∼10 h and ∼40 min ion periodicities in HI-SCALE measurements were present in several cases during the whole period examined (October 2003 to March 2004) and were found to be more evident during some special conditions, for instance during enhanced fluxes around the start (forward shock) and the end (reverse shock) of CIRs. We infer that the Jovian magnetosphere was triggered by the impact of the CIRs, after the Halloween events, and it was (a) a principal source of forward and reverse shock-associated ion flux structures and (b) the cause of generation of ∼10 h quasi-periodic magnetic field and plasma modulation observed by Ulysses at those times.     

Photoelectron flux variations observed from the FAST satellite

This paper examines high resolution (ΔE/E = 0.15) photoelectron energy spectra from 10 eV to 1 keV, created by solar irradiances between 1.2 and 120 nm. The observations were made from the FAST satellite at ∼3000 km, equatorward of the auroral oval for the July–August, 2002 solar rotation. These data are compared with the solar irradiance observed by the Solar EUV Experiment (SEE) on the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite and fluxes calculated using the Field Line Interhemispheric Plasma (FLIP) code. The 41 eV photoelectron flux, which corresponds to solar EUV fluxes near 20 nm, shows a clear solar rotation variation in very good agreement with the EUV flux measurements. This offers the possibility that the 41 eV photoelectron flux could be used as a check on measured solar EUV fluxes near 20 nm. Because of unexpected noise, the solar rotation signal is not evident in the integral photoelectron flux between 156 and 1000 eV corresponding to EUV wavelengths between 0.1 and 7 nm measured by the SEE instrument. Examination of daily averaged photoelectron fluxes at energies between 25 and 500 eV show significant changes in the photoelectron spectra in response X and M class flares. The intensity of photoelectrons produced in this energy region is primarily due to two very narrow EUV wavelength regions at 2.3 and 3 nm driving Auger photoionization in O at 500 eV and N₂ at ∼360 eV. Comparison of calculated and daily averaged electron fluxes shows that the HEUVAC model solar spectrum used in the FLIP code does not reproduce the observed variations in photoelectron intensity. In principle, the 21 discrete photoelectron energy channels could be used to improve the reliability of the solar EUV fluxes at 2.3 and 3 nm inferred from broad band observations. In practice, orbital biases in the way the data were accumulated and/or noise signals arising from natural and anthropogenic longitudinally restricted sources of ionization complicate the application of this technique.     

XMM-Newton observations across the Cygnus Loop from northeastern rim to southwestern rim

We have observed the Cygnus Loop from the northeast (NE) to the southwest (SW) with XMM-Newton. We extracted spectra from 3′-spaced annular regions across the Loop and fitted them either with a one-kT_e-component non-equilibrium ionization (NEI) model or with two-kT_e-component NEI model. We found that the two-kT_e-component model yields significantly better fits in almost all the spectra than the one-kT_e-component model. Judging from the abundances, the high-kT_e-component in the two-temperature model must be fossil ejecta while the low-kT_e-component comes from the swept-up interstellar medium (ISM). The distributions of Ne, Mg, Si, and S for fossil ejecta appear to retain the onion-skin structure at the time of a supernova explosion, suggesting that significant overturning of the ejecta has not occurred yet. Comparing the relative abundances of fossil ejecta estimated for the entire Cygnus Loop with those from theoretical calculations for Type-II SN, the mass of the progenitor star is likely to be ∼13 M_⊙. The trends of the radial profiles of kT_e and emission integral for the swept-up ISM are adequately described by the Sedov model, suggesting that the swept-up ISM is concentrated in a shell-like structure. Comparing our data with the Sedov model, we found the ambient medium density to be ∼0.7 cm⁻³. Then, we estimated the total mass of the swept-up ISM and the age of the remnant to be ∼130 M_⊙ and 13,000 years, respectively. Note that if the explosion occurred within a stellar wind cavity, then the density in the cavity, the total swept-up mass in the cavity, and the age of the remnant are estimated to be ∼0.14 cm⁻³, ∼25 M_⊙, and ∼10,000 years, respectively.