A calibration neutron monitor: Statistical accuracy and environmental sensitivity

This paper is a follow-on of that of Krüger et al. [Krüger, H., Moraal, H., Bieber, J.W., Clem, J.M., Evenson, P.A., Pyle, K.R., Duldig, M.L., Humble, J.E. A calibration neutron monitor: energy response and instrumental temperature sensitivity. J. Geophys. Res. 113, A08101, doi:10.1029/2008JA013229, 2008], that describes the characteristics of a pair of calibration neutron monitors that were developed to intercalibrate the count rates of the world’s neutron monitors against each other. Such an intercalibration will allow the calculation of energy (rigidity) spectra, which will enhance the quality of the neutron monitor data. Krüger et al. (2008) investigated the energy and temperature response of the calibrators. This paper studies the statistical accuracy of the calibration procedure, its repeatability, and the sensitivity to its environment. The paper concludes with a calibration procedure that can minimise the uncertainties caused by these five effects, or at least correct for them.     

CME-geometry and cosmic-ray anisotropy observed by a prototype muon detector network

We analyze the cosmic-ray anisotropy observed by a prototype network of muon detectors during geomagnetic storms associated with coronal mass ejections (CMEs). The network currently consists of multidirectional surface muon detectors at Nagoya (Japan) and Hobart (Australia), together with a prototype detector at São Martinho (Brazil) which has been in operation since March, 2001. In this report, we analyze the anisotropy recorded in both the muon detector and neutron monitor (the Spaceship Earth) networks and find significant enhancements of cosmic-ray anisotropy during geomagnetic storms. Following the analysis by Bieber and Evenson [Bieber, J.W., Evenson, P. CME geometry in relation to cosmic ray anisotropy. Geophys. Res. Lett. 25 (1998) 2955–2958] for the neutron monitor data at 10 GeV, we also derive cosmic-ray density gradients from muon data at higher-energy (50 GeV), possibly reflecting the larger-scale geometry of CMEs causing geomagnetic storms. We particularly find in some events the anisotropy enhancement clearly starting prior to the storm onset in both the muon and neutron data. This is the first result of the CME-geometry derived from simultaneous observations of the anisotropy with networks of multidirectional muon detectors and neutron monitors.     

A heuristic derivation of an improved analytical theory for perpendicular diffusion of charged particles

A heuristic derivation of an improved nonlinear guiding center theory for perpendicular diffusion of charged particles is presented. This new derivation complements previous work which is based on an indirect solution of the Fokker–Planck equation. The new derivation confirms the improved theory for diffusion of charged particles and makes the validity of the theory plausible.     

Effective non-vertical and apparent cutoff rigidities for a cosmic ray latitude survey from Antarctica to Italy in minimum of solar activity

In this paper we will report the results of the computation of cutoff rigidities of vertical and non-vertical incident cosmic ray particles. Non-vertical effective cutoff rigidities have been computed by tracing particle trajectories through the “real” geomagnetic magnetic field comprising the International Geomagnetic Reference Field model (IGRF95, IAGA Division 5 Working Group 8, 1996: Sabaka, T.J., Langel, R.A., Baldwin, R.T., Conrad, J.A. The geomagnetic field, 1900–1995, including the large scale fields from magnetospheric sources and NASA candidate models for the 1995 IGRF revision. J. Geomag. Geoelect. 49, 157–206, 1997.) and the Tsyganenko [Tsyganenko, N.A. A magnetospheric magnetic field model with a warped tail current sheet. Planet. Space Sci. 37, 5–20, 1989.] magnetosphere model. The computation have been done for the backward route (from Antarctica to Italy) of the Italian Antarctic ship survey 1996–1997, for geographic points corresponding to the daily average coordinates of the ship; for zenith angles 15°, 30°, 45° and 60°, and azimuth angles from 0° to 360° in steps of 45°. By means of the obtained non-vertical cutoffs the apparent cutoff rigidities have been calculated. The information on integral multiplicities of secondary neutrons detected by the neutron monitor in dependence of the zenith angle of incoming primary cosmic ray particles have also been used. This information is based on the theoretical calculations of meson-nuclear cascades of primary protons with different rigidities arriving to the Earth’s atmosphere at the zenith angles of 0°, 15°, 30°, 45°, 60° and 75°. The difference between the computed apparent and vertical cutoff rigidities reaches ∼1 GV at rigidities >7–8 GV. At rigidities of 10–16 GV, the difference between the apparent and vertical cutoff rigidities is larger than that obtained earlier by Clem et al. [Clem, J.M., Bieber, J.W., Duldig, M., Evenson, P., Hall, D., Humble, J.E. Contribution of obliquely incident particles to neutron monitor counting rate. J. Geophys. Res. 102, 26919–26926, 1997.] and Dorman et al. [Dorman, L.I., Villoresi, G., Iucci, N., Parisi, M., Tyasto, M.I., Danilova, O.A., Ptitsyna, N.G. Cosmic ray survey to Antarctica and coupling functions for neutron component near solar minimum (1996–1997), 3. Geomagnetic effects and coupling functions. J. Geophys. Res. 105, 21047–21056, 2000.].     

Perpendicular transport of charged particles: Results for the unified nonlinear transport theory derived from the Newton–Lorentz equation

The perpendicular diffusion coefficient is calculated by combining a recently developed Unified Nonlinear Transport (UNLT) theory with the Newton–Lorentz equation. The total perpendicular mean free path can be described as a combination of a guiding center contribution and a microscopic contribution. It is shown that the total mean free path depends strongly on the energy range of the turbulence power spectrum and on particle energy. Further, a slab/2D composite model is used to investigate the influence of each contribution to the total mean free path for a quasi-3D turbulence model. For pure 2D turbulence the UNLT reduces to the NLGC-theory. For pure slab turbulence the guiding center contribution is subdiffusive in accordance with simulations and the theorem on reduced dimensionality. Conversely, the microscopic contribution is non-zero, which has to be interpreted as normal diffusion.     

Particle transport in a vortex medium

We study numerically particle transport in a two-dimensional coherent vortex field. Reasonable agreement exists between previously derived radial transport coefficients for energetic particles (Verkhoglyadova, O.P., le Roux, J.A. Particle diffusion on vortices in nearly incompressible magnetohydrodynamics. Astrophys. J. 602, 1002–1005, 2004a; Verkhoglyadova, O.P., le Roux, J.A. Cosmic ray transport in a vortex flow. IGPP-UCR Conf. Physics of the Outer Heliosphere (Riverside), AIP Conf. Proc., pp. 243–248, 2004b) and results of numerical simulations. Different physical factors controlling particle momentum change and drifts are analysed. It is shown that the vortex electric field is the main cause of trapped particle motion. Drift due to magnetic field inhomogeneity predominantly disturbs free particle gyroorbit along background magnetic field in the vicinity of a vortex. Our simulations show the development of a subdiffusion regime.     

Flux upper limits of diffuse TeV gamma rays from the Galactic plane using the effective area of the Tibet-II and -III arrays

We obtained new upper limits on the diffuse gamma rays from the inner Galactic (IG) and outer Galactic (OG) planes in 3–10 TeV region, using the Tibet air shower data and new Monte Carlo simulation results. A difference of the effective area of the air-shower array for observing gamma rays and cosmic rays was carefully taken into account in this analysis, resulting in that the flux upper limits of the diffuse TeV gamma rays were reduced by factors of 4.0–3.7 for 3–10 TeV than those in our previous results (Amenomori, M., Ayabe, S., Cui, S.W., et al. Observation of multi-TeV diffuse gamma rays from the Galactic plane with the Tibet air shower array. Astrophys. J. 580, 887–895, 2002.). This new result suggests that the inverse power index of the energy spectrum of source electrons responsible for generating diffuse TeV gamma rays through inverse Compton effect should be steeper than 2.2 and 2.1 for IG and OG planes, respectively, with 99%C.L.     

GRB 980425, SN1998BW and the EMBH model

The EMBH model, previously developed using GRB 991216 as a prototype, is here applied to GRB 980425. We fit the luminosity observed in the 40–700 keV, 2–26 keV and 2–10 keV bands by the BeppoSAX satellite. In addition we present a novel scenario in which the supernova SN1998bw is the outcome of an “induced gravitational collapse” triggered by GRB 980425, in agreement with the GRB-Supernova Time Sequence (GSTS) paradigm [Ruffini, R., Bianco, C.L., Chardonnet, P., Fraschetti, F., Xue, S.-S. On a possible GRB-supernova time sequence. Astrophys. J. 555, L117–L120, 2001c]. A further outcome of this astrophysically exceptional sequence of events is the formation of a young neutron star generated by the SN1998bw event. A coordinated observational activity is recommended to further enlighten the underlying scenario of this most unique astrophysical system.     

On the effects of dynamical turbulence on the perpendicular diffusion of low-energy cosmic ray electrons

Diffusion perpendicular to the heliospheric magnetic field plays an integral role in the transport of charged particles in the heliosphere. In this study the perpendicular diffusion coefficient of low-energy cosmic ray electrons is calculated, using an equation derived from the random ballistic decorrelation interpretation of nonlinear guiding centre theory. An observationally motivated 2D turbulence power spectrum is assumed and the effects of various turbulence inputs on the resulting perpendicular diffusion coefficient are investigated. The perpendicular diffusion coefficients are first determined at 1 AU, for both magnetostatic and dynamical turbulence conditions. These solutions are also evaluated for radial distances of 0.1 AU to 10 AU to further investigate the values of the perpendicular diffusion coefficients in the very inner heliosphere. The results of this study show that the dissipation range of the turbulence power spectrum provides a negligible contribution towards the perpendicular diffusion coefficient, and that solutions derived using only the energy containing range serve as good approximations for solutions derived assuming the full 2D turbulence power spectrum. Finally, it is shown that the effects of dynamical turbulence, as considered in the present study, do not affect the perpendicular diffusion coefficients derived from the scattering theory considered here.     

Galactic cosmic ray gradients in the ecliptic plane and at high latitudes during two solar cycles (meteorite data)

Results of investigations of cosmogenic isotope radioactivity in chondrites fallen to earth during two solar cycles are presented. The data obtained on radial and latitudinal gradients cover the period 1955–1976, heliocentric distances from 1.03 AU to 3.33 AU, and heliographic latitudes from 23°S to 16°N. The dependence of radial and latitudinal gradients on the phase of solar activity is established, as well as a north-south asymmetry during a certain period after the inversion of the general solar magnetic field in 1969.     

The heliospheric modulation of 3–10 MeV electrons: Modeling of changes in the solar wind speed in relation to perpendicular polar diffusion

The discrepancy between cosmic ray model predictions representing solar minimum conditions in the heliosphere and the 3–10 MeV post-1998 electrons observations by the Kiel Electron Telescope (KET) onboard Ulysses suggests the need for consistent changes in model parameters with increasing solar activity. In order to reduce this discrepancy, an effort is made to model the KET observations realistically during periods of increased solar activity by applying an advanced three-dimensional, steady-state electron modulation model based on Parker’s transport equation including the Jovian electron source. Some elements of the diffusion tensor which were not previously emphasized are revisited. A new relation is also established between the latitudinal dependence of the solar wind speed and the perpendicular polar diffusion. Based on this relation, a transition of an average solar wind speed from solar minimum to solar maximum conditions, as observed on board the Ulysses spacecraft, is modeled on the concept of the time-evolution of large polar coronal holes. These changes are correlated to different scenarios of the enhancement of perpendicular polar diffusion. Effects of these scenarios are illustrated, as a series of steady-state solutions, on the computed 7 MeV Jovian and galactic electrons in comparison with 3–10 MeV electrons observed from the period 1998 to the end of 2003. It is shown that this approach improves compatibility with the KET observations but it also points to the need for a time-dependent electron modulation model to fully describe modulation during moderate to extreme solar maximum conditions.     

Modeling of the outer heliosphere with the realistic solar cycle

Time-dependent kinetic-continuum model of the solar wind interaction with the two-component local interstellar cloud (LIC) has been developed recently [Izmodenov, V., Malama, Y.G., Ruderman, M.S. Solar cycle influence on the interaction of the solar wind with local interstellar cloud. Astron. Astrophys. 429, 1069–1080, 2005a.]. Here, we adopted this model to the realistic solar cycle, when the solar wind parameters at the Earth’s orbit are taken from space data. This paper focuses on the results related to the termination shock (TS) excursion with the solar cycle that may help to understand Voyager 1 data obtained at and after the crossing of the termination shock and to predict the time of the TS crossing by Voyager 2.     

Evaluation of the IRI-2007 model options for the topside electron density

The International Reference Ionosphere (IRI) 2007 provides two new options for the topside electron density profile: (a) a correction of the IRI-2001 model, and (b) the NeQuick topside formula. We use the large volume of Alouette 1, 2 and ISIS 1, 2 topside sounder data to evaluate these two new options with special emphasis on the uppermost topside where IRI-2001 showed the largest discrepancies. We will also study the accurate representation of profiles in the equatorial anomaly region where the profile function has to accommodate two latitudinal maxima (crests) at lower altitudes but only a single maximum (at the equator) higher up. In addition to IRI-2001 and the two new IRI-2007 options we also include the Intercosmos-based topside model of Triskova, Truhlik, and Smilauer [Triskova, L., Truhlik, V., Smilauer, J. An empirical topside electron density model for calculation of absolute ion densities in IRI. Adv. Space Res. 37 (5), 928–934, 2006] (TTS model) in our analysis. We find that overall IRI-2007-NeQ gives the best results but IRI-2007-corrected provides a more realistic representation of the altitudinal–latitudinal structure in the equatorial anomaly region. The applicability of the TTS model is limited by the fact that it is not normalized to the F2 peak density and height.     

Modelling of galactic Carbon in an asymmetrical heliosphere: Effects of asymmetrical modulation conditions

Observations of galactic cosmic rays (GCRs) from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that in addition to a possible global asymmetry in the north–south dimensions (meridional plane) of the heliosphere, it is also possible that different modulation (turbulence) conditions could exist between the two hemispheres of the heliosphere. We focus on illustrating the effects on GCR Carbon of asymmetrical modulation conditions combined with a heliosheath thickness that has a significant dependence on heliolatitude. To reflect different modulation conditions between the two heliospheric hemispheres in our numerical model, the enhancement of both polar and radial perpendicular diffusion off the ecliptic plane is assumed to differ from heliographic pole to pole. The computed radial GCR intensities at polar angles of 55° (approximating the Voyager 1 direction) and 125° (approximating the Voyager 2 direction) are compared at different energies and for both particle drift cycles. This is done in the context of illustrating how different values of the enhancement of both polar and radial perpendicular diffusion between the two hemispheres contribute to causing differences in radial intensities during solar minimum and moderate maximum conditions. We find that in the A > 0 cycle these differences between 55° and 125° change both quantitatively and qualitatively for the assumed asymmetrical modulation condition as reflected by polar diffusion, while in the A < 0 cycle, minute quantitative differences are obtained. However, when both polar and radial perpendicular diffusion have significant latitude dependences, major differences in radial intensities between the two polar angles are obtained in both polarity cycles. Furthermore, significant differences in radial intensity gradients obtained in the heliosheath at lower energies may suggest that the solar wind turbulence at and beyond the solar wind termination shock must have a larger latitudinal dependence.     

Modeling gamma radiation from millisecond pulsars

We present spectral and directional properties of gamma radiation obtained by numerical modeling of two millisecond pulsars, PSR J0218+4232 and PSR J0437-4715. The general relativistic space charge limited flow model of Muslimov and Harding [Toward the quasi-steady state electrodynamics of a neutron star. Astrophys. J. 485, 735–746, 1997] in its version of an unscreened accelerating electric field was used. The results give both general ideas and constraints on the specific character of magnetospheric activity in each of these two pulsars.     

Converging motion of conjugate flaring kernels during two large solar flares

In this paper, we analyze the footpoint motion of two large solar flares using observations made by the Transition Region and Coronal Explorer (TRACE) and Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The two flares are the M5.7 flare of March 14, 2002 and the X10 flare of October 29, 2003. They are both classical two-ribbon flares as observed in TRACE 1600 or 171 Å images and have long-duration conjugate hard X-ray (HXR) footpoint emission. We use the ‘center-of-mass’ method to locate the centroids of the UV/EUV flare ribbons. The results are: (1) The conjugate UV/EUV ribbons and HXR footpoints of the two flares show a converging (inward) motion during the impulsive phase. For the two flares, the converging motion lasts about 3 and 10 min, respectively. The usual separation (outward) motion for the flare ribbons and footpoints take place only after the converging motion. (2) During the inward and the outward motion, the conjugate ribbons and footpoints of the two events exhibit a strong unshear motion. In obtaining above results, TRACE UV/EUV and RHESSI HXR data show an overall agreement. The two events demonstrate that the magnetic reconnection for the flares occurs in highly sheared magnetic field. Furthermore, the results support the magnetic model constructed by Ji et al. [Ji, H., Huang, G., Wang, H. Astrophys. J. 660, 893–900, 2007], who proposed that the contracting motion of flaring loops is the signature of the relaxation of sheared magnetic fields.     

Midday bottomside electron density profiles during moderate solar activity and comparison with IRI-2001

Bottomside electron density (Ne-h) profiles during midday (10–14 h) are analyzed using modern digital ionosonde observations at a low-middle latitude station, New Delhi (28.6N, 77.2E, dip 42.4N), for the period from January 2003 to December 2003, pertaining to moderate solar activity (MSA). Each individual profile is normalized with respect to the peak height and density (hmF2, NmF2) of the F2-region. These profiles are compared with those obtained from the International Reference Ionosphere (IRI-2001) model. Bilitza [Bilitza, D. International Reference Ionosphere 2000. Radio Sci. 36 (2), 261–275, 2001] using both the options namely: Gulyaeva’s model [Gulyaeva, T.L. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7 (6) 39–48, 1987] and B0 Tab. option [Bilitza, D., Radicella, S.M., Reinisch, B.W., Adeniyi, J.O., Mosert Gonzalez, M.E., Zhang, S.R., Obrou, O. New B0 and B1 models for IRI. Adv. Space Res. 25 (1), 89–95, 2000]. The study reveals that during summer and equinox, the IRI model with B0 Tab. option in general, produces better agreement with the observed median profiles, while the IRI predictions using Gulyaeva’s option, overestimate the electron density distribution at all the heights below the F2-peak. However, during winter, in general, the IRI model, using both the options, reveals shows fairly good agreement with the observations.     

Radial and latitudinal gradients of galactic cosmic rays in the heliosphere at solar maximum

Our understanding of galactic cosmic ray (GCR) modulation has advanced greatly in the last three decades. However, we still need an appropriate knowledge of the GCR intensity gradient. Numerical simulations of the transport particle equation allow interpretation of cosmic ray intensities in the heliosphere. We use the numerical solution of the GCR transport equation during solar maximum epoch to compute the radial and latitudinal gradients. Our analysis indicates that adiabatic energy loss plays an important role in the radial distribution of GCR in the inner heliosphere, while in the outer region the diffusion and convection are the relevant processes. The latitudinal gradient is small.     

Numeric simulations of in situ observations of interacting ejecta near 1 AU

We use a simple numerical model (González-Esparza, J.A., Santillán, A., Ferrer, J. A numerical study of the interaction between two ejecta in the interplanetary medium: one and two dimensional hydrodynamic simulations, Ann. Geophys. 22, 3741–3749, 2004) to study the evolution of three events in the solar wind reported by Wang et al. (Wang, Y.M., Ye, P.Z., Wang, S. Multiple magnetic clouds: several examples during March–April 2001. J. Geophys. Res. 108, 1370, 2003, doi:10.1029/2003JA009850), where two interacting ejecta detected in situ by ACE near 1 AU were related to CMEs observed previously by SOHO-LASCO. The study is based on a 1-D hydrodynamic model using the ZEUS code (Stone, J.M., Norman, M. ZEUS 2-D: A radiation magnetohydrodynamics code for astrophysical flows in two dimensions, I, the hydrodynamics algorithms and tests, Astrophys. J. 80, 753, 1992). Although this model cannot address either the magnetic field dynamics or the complex geometrical effects intrinsic in the three-dimensional nature of the phenomena, it illuminates the transferring of momentum and evolution of interacting large-scale solar wind disturbances in those cases where there is no merging (magnetic reconnection) between the two ejecta. This model can reproduce, in some cases, characteristics of the events such as transit times and flow signatures as inferred from the two-point measurements from spacecraft.