Real-time global MHD simulation of the solar wind interaction with the earth’s magnetosphere

We have developed a real-time global MHD (magnetohydrodynamics) simulation of the solar wind interaction with the earth’s magnetosphere. By adopting the real-time solar wind parameters and interplanetary magnetic field (IMF) observed routinely by the ACE (Advanced Composition Explorer) spacecraft, responses of the magnetosphere are calculated with MHD code. The simulation is carried out routinely on the super computer system at National Institute of Information and Communications Technology (NICT), Japan. The visualized images of the magnetic field lines around the earth, pressure distribution on the meridian plane, and the conductivity of the polar ionosphere, can be referred to on the web site (http://www2.nict.go.jp/y/y223/simulation/realtime/).The results show that various magnetospheric activities are almost reproduced qualitatively. They also give us information how geomagnetic disturbances develop in the magnetosphere in relation with the ionosphere. From the viewpoint of space weather, the real-time simulation helps us to understand the whole image in the current condition of the magnetosphere. To evaluate the simulation results, we compare the AE indices derived from the simulation and observations. The simulation and observation agree well for quiet days and isolated substorm cases in general.     

A multiple type-II burst associated with a coronal transient and its MHD simulation

A large coronal transient took place on 8 May 1981. The transient was related to an M7.7/2B flare and was associated with at least two coronal type II bursts. The velocities of the type II bursts were in the range 1100–1800 kms^?1, in excess of the transient velocity of 500–1000 kms^?1. Two dimensional positions of the type II radio sources are available from both the Clark Lake and the Culgoora Radio Observatories. We carry out two dimensional MHD simulations of the event, taking into account the observed velocity, position, and size of the type II bursts. We simulate the multiple shocks observed during the event and their interaction, and discuss some results of the simulation.     

Application of an MHD simulation to the study of substorms

The substorm mechanism is studied by the numerical solutions obtained from a resistive magnetohydrodynamic (MHD) simulation. After a southward turning of the interplanetary magnetic field (IMF), the simulation results reproduce observed features of the growth phase. The numerical solutions show that the plasma sheet thinning during the growth phase is formed under the dynamic balance between the flux pileup from the midtail and the flux removal toward the dayside controlled by the convection in the magnetosphere-ionosphere (M-I) coupling system. After the growth phase, dipolarization is generated in the near-earth tail accompanied by a plasma injection into the inner magnetosphere, the formation of plasmoid in the midtail, and the enhancement of the nightside field-aligned currents (FACs). The direct cause of this onset is the state (phase space) transition of the convection system from a thinned state to a dipolarized state associated with a self-organization in the nonlinear system.     

MHD simulation of magnetic field configuration above the active region NOAA 10365

The current sheet (CS) creation before a flare in the vicinity of a singular line above the active region NOAA 10365 is shown in numerical experiments. Such a way the possibility of energy accumulation for a solar flare is demonstrated. These data and results of observation confirm the electrodynamical solar flare model that explains solar flares and CME appearance during CS disruption. The model explains also all phenomena observed in flares. For correct reproduction of the real boundary conditions the magnetic flux between spots should be taken into account. The full system of 3D MHD equations are solved using the PERESVET code. For setting the boundary conditions the method of photospheric magnetic maps is used. Such a method permits to take into account all evolution of photospherical magnetic field during several days before the flare.     

Relationships between energetic storm particle events and interplanetary shocks driven by full and partial halo coronal mass ejections

We have analysed energetic storm particle (ESP) events in 116 interplanetary (IP) shocks driven by front-side full and partial halo coronal mass ejections (CMEs) with speeds $>$400 km s^?1during the years 1996–2015. We investigated the occurrence and relationships of ESP events with several parameters describing the IP shocks, and the associated CMEs, type II radio bursts, and solar energetic particle (SEP) events. Most of the shocks (57 %) were associated with an ESP event at proton energies $>$1 MeV.The shock transit speeds from the Sun to 1 AU of the shocks associated with an ESP event were significantly greater than those of the shocks without an ESP event, and best distinguished these two groups of shocks from each other. The occurrence and maximum intensity of the ESP events also had the strongest dependence on the shock transit speed compared to the other parameters investigated. The correlation coefficient between ESP peak intensities and shock transit speeds was highest (0.73 $\pm$ 0.04) at 6.2 MeV. Weaker dependences were found on the shock speed at 1 AU, Alfvénic and magnetosonic Mach numbers, shock compression ratio, and CME speed. On average all these parameters were significantly different for shocks capable to accelerate ESPs compared to shocks not associated with ESPs, while the differences in the shock normal angle and in the width and longitude of the CMEs were insignificant.The CME-driven shocks producing energetic decametric–hectometric (DH) type II radio bursts and high-intensity SEP events proved to produce also more frequently ESP events with larger particle flux enhancements than other shocks. Together with the shock transit speed, the characteristics of solar DH type II radio bursts and SEP events play an important role in the occurrence and maximum intensity of ESP events at 1 AU.     

Comparison of global chlorophyll concentrations using MODIS data

Chlorophyll concentrations derived from satellite borne ocean color sensors data provide an idea of the distribution of phytoplanktons across the oceans which help us in understanding the spatial and temporal dynamics of phytoplanktons. The changes in the patterns of distribution and abundance of the planktons have significant impact on the entire ecosystem and play a key role in the global carbon cycle. In this paper, we have analyzed annual and seasonal chlorophyll concentrations retrieved from MODIS data for the periods March 2000–October 2003, which reveal the spatial and seasonal distribution of chlorophyll concentrations across the global oceans. Chlorophyll concentrations anomaly indicate that chlorophyll concentrations in almost all ocean regions responded similarly.     

MHD interaction of pulsar wind nebulae with SNRs and with the ISM

In the late 1960s the discovery of the Crab pulsar in its associated supernova remnant, launched a new field in supernova remnant research: the study of pulsar-driven or plerionic supernova remnants. In these type of remnants, the relativistic wind emitted by the pulsar, blows a pulsar wind nebula into the interior of its supernova remnant. Now, more then forty years after the discovery of the Crab pulsar, there are more then fifty plerionic supernova remnants known, due to the ever-increasing capacity of observational facilities. These observational studies reveal a Zoo of complex morphologies over a wide range of frequencies, indicating the significance of the interaction between a pulsar wind nebula with its surrounding supernova remnant. A pulsar which gained a kick velocity at birth, will ultimately break outside of its remnant, after which the pulsar wind nebula interacts directly with the interstellar medium. In general these pulsar wind nebulae are bounded by a bow shock, due to the supersonic motion of the pulsar. There are a few examples known of these pulsar-powered bow shocks, a number which is slowly increasing.I will review our current understanding of the different evolutionary stages of a pulsar wind nebula as it is interacting with its associated supernova remnant. Therefore, I will discuss both analytical and more recent numerical (M)HD models. The four main stages of a pulsar wind nebula are: the supersonic expansion stage, the reverse shock interaction stage, the subsonic expansion stage and ultimately the stage when the head of the bubble is bounded by a bow shock, due to the supersonic motion of the pulsar. Ultimately this pulsar wind nebula bow shock will break through its associated remnant, after which the pulsar-powered bow shock will interact directly with the interstellar medium. I will discuss recent numerical models from these type of pulsar wind nebulae and their morphology.     

Evaluation of substorm models with Cluster observations of plasma flow reversal in the magnetotail

We evaluate two prevailing substorm models with an event of plasma flow reversal from tailward to Earthward detected by Cluster at the downstream distance of ∼19 R_E in the magnetotail during a substorm on August 22, 2001. We use the unique capability of Cluster measurements in determining gradients to examine the associated current density, Lorentz force, and current dissipation/dynamo term. In association with plasma flow reversal, it is found that (1) there was no clear quadrupole magnetic perturbation signature, (2) the x-component of the Lorentz force did not change sign, (3) the y-component of the product of the current density and the electric field was occasionally negative indicative of a dynamo effect, and (4) the timing sequence of flow reversal from the Cluster configuration did not match tailward motion of a single plasma flow source. These observations are consistent with the near-Earth initiation model for substorms with multiple current disruption sites moving progressively tailward near the late stage of substorm expansion.     

Space storm measurements of the July 2005 solar extreme events from the low corona to the Earth

The Athens Neutron Monitor Data Processing (ANMODAP) Center recorded an unusual Forbush decrease with a sharp enhancement of cosmic ray intensity right after the main phase of the Forbush decrease on 16 July 2005, followed by a second decrease within less than 12 h. This exceptional event is neither a ground level enhancement nor a geomagnetic effect in cosmic rays. It rather appears as the effect of a special structure of interplanetary disturbances originating from a group of coronal mass ejections (CMEs) in the 13–14 July 2005 period. The initiation of the CMEs was accompanied by type IV radio bursts and intense solar flares (SFs) on the west solar limb (AR 786); this group of energetic phenomena appears under the label of Solar Extreme Events of July 2005. We study the characteristics of these events using combined data from Earth (the ARTEMIS IV radioheliograph, the Athens Neutron Monitor (ANMODAP)), space (WIND/WAVES) and data archives. We propose an interpretation of the unusual Forbush profile in terms of a magnetic structure and a succession of interplanetary shocks interacting with the magnetosphere.     

Ion dynamics associated with Alfven wave in the near-Earth magnetotail: Two-dimensional global hybrid simulation

Ion dynamics in the near-Earth magnetotail region is examined during periods of fast Earthward flow with a two-dimensional (2-D) global-scale hybrid simulation. The simulation shows that shear Alfven waves are generated at x ∼ −10R_E, where the strong earthward flow is arrested by the dipole field, and propagate along field lines from the equator to both southern and northern polar ionosphere. Non-gyrotropic ion velocity distributions occur where the large-amplitude Alfven waves are dominant. The simulation indicates that the Alfven waves are generated by interaction of the fast earthward flow with the stationary near-Earth plasma. Beam ions are found to be pitch-angle scattered and trapped in the wave field, leading to the non-gyrotropic ion distributions in the high-latitude plasma sheet boundary. In addition, significant particle heating and acceleration are found to occur behind the dipolarization front due to the effect of wave turbulence.     

Numerical simulation—assessment of FGGE data with regard to their assimilation in a global data set

The data-assimilation system at the European Centre for Medium Range Forecasts (ECMWF) is presented. The data assimilation system is used to process the FGGE level II-b data and to provide global 3-dimensional analyses for every 6th hour through the whole FGGE year. The quality and the excellent coverage of the data has meant a substantial improvement in our understanding of the global atmosphere with major implication for climate and extended weather forecasting.Results from an observing system experiment only using data observed and collected from space is presented. It is found that fully space-oriented systems are possible although forecast quality is higher in areas where radiosondes are available with a satisfactory density.     

Comparison of global morphologies of vertical ion convergence and sporadic E occurrence rate

This study employs the Horizontal Wind model 2014 to calculate the intensity of the Vertical ion convergence (VIC) to examine the role of neutral wind shear in the global distribution of the Es occurrence rate. The global S4 observations from the FORMOSAT-3/COSMIC GPS Occultation Experiment (GOX) are used to derive global information on the Es occurrence. Both the simulation and observation suggest that the Es highly occurs in the mid-latitude regions of the summer hemisphere.     

Multi-point observations of earthward fast flow in the plasma sheet by virtual satellites located in the MHD simulation domain

Time profiles of some physical values in earthward fast flows in the plasma sheet are observed at three dimensionally different positions by employing virtual satellites located in the three-dimensional magnetohydrodynamic simulation domain, and these simulations are done on the basis of the spontaneous fast reconnection model. In the spontaneous fast reconnection evolution, the width of the flow channel is narrow in the dawn-dusk direction, and it does not spread until the plasma collides with the magnetic loop. The enhancements in B_z and V_x are larger at the center of the fast flow channel than those at its dawn and dusk edges, reflecting the differences in the reconnection rate in the diffusion region. The enhancement in V_x is shorter near the plasma sheet boundary layer than that near the neutral sheet, reflecting the changes in the thickness of the flow channel.     

Effects of a “day-time” substorm on the ionosphere and radio propagation

Propagation mechanisms of lateral (non-great-circle) signals on a high-latitude HF radio path during magnetospheric substorms that occurred in the day-time have been considered. The path is equipped with oblique ionospheric sounding (OIS) from Murmansk to St. Petersburg. The OIS method gives the possibility to determine propagation modes, MOF (maximum observed frequency) values, signal delays, etc. Data of the CUTLASS radar, the IMAGE magnetometer system, the Finnish riometer chain, and the Tromso ionosonde were also used for the analysis. The main results are the following: (1) the lateral signal propagation takes place, as a rule, if the path midpoint is located near the irregularity region that moves sharply from high to low latitudes. The lateral signal propagation appearing during day-time is a new effect. (2) Formation of dense field-aligned irregularities during a substorm leads to decreasing F2MOF values on radio paths. These results can be useful for problems of radiolocation, HF communications and navigation.     

Significant proton events of solar cycle 22 and a comparison with events of previous solar cycles.

The sun has produced several high energy and large fluence solar proton events during solar cycle 22. This recent activity is similar to activity that occurred in the 19th solar cycle before the advent of routine space measurements. In a review of the recent events and a comparison with significant solar proton events of previous solar cycles, it appears that the 20th and 21st solar cycles were deficient in the total fluence of solar particles as detected at the earth. Therefore, when models of maximum solar proton fluxes to be encountered for deep space missions are developed, solar proton data acquired during the present solar cycle should be incorporated.     

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.     

MHD modeling of the outer heliosphere: Achievements and challenges

An overview is presented of magnetic-field-related effects in the solar wind (SW) interaction with the local interstellar medium (LISM) and the different theoretical approaches used in their investigation. We discuss the possibility that the interstellar magnetic field (ISMF) introduces north–south and east–west asymmetries of the heliosphere, which might explain observational data obtained by the Voyager 1 and Voyager 2 spacecraft. The SW–LISM interaction parameters that are responsible for the deflection of the interstellar neutral hydrogen flow from the direction of propagation of neutral helium in the inner heliosheath are outlined. The possibility of a strong ISMF, which increases the heliospheric asymmetry and the H–He flow deflection, is discussed. The effect of the combination of a slow-fast solar wind during solar minimum over the Sun’s 11-year activity cycle is illustrated. The consequences of a tilt between the Sun’s magnetic and rotational axes are analyzed. Band-like areas of an increased magnetic field distribution in the outer heliosheath are sought in order to discover regions of possible 2–3 kHz radio emission.     

Solar-interplanetary circumstances associated with the major events in March and June 1991 and comparison with similar events of previous solar cycles.

The solar activity and geomagnetic storm events of March and June 1991 were associated with the appearance of an enhanced particle flux in the trapped radiation belts as discovered by the CRRES satellite observations and later measured by shuttle radiation detectors. The solar-interplanetary conditions associated with these events appear to be a major sequence of activity near the sun's central meridian generating powerful fast interplanetary shocks resulting in major perturbations to the magnetosphere. The solar-interplanetary events in 1991 are discussed and compared to similar activity in the past such as the events in February 1986, August 1972, July 1961, November 1960, and July 1959.     

Cluster observations of earthward propagating plasmoid and flux ropes in the near-tail during the course of a substorm event

We analyze the magnetic structures in the near-tail at Xgsm = −17.5 Re on September 19, 2003 by Cluster. During the course of a substorm event, the earthward propagating plasmoid and flux ropes in the near-tail are observed. The earthward propagating plasmoid is associated with the bipolar Bz and By signatures. The two flux ropes are embedded within the earthward plasma flows, which might be referred to the population as ‘‘BBF-type’’ flux ropes. The first flux rope diameter is about 0.7 Re and duration based upon the Bz signature is ∼20 s, while the second one diameter is about 1.4 Re and duration is ∼30 s. The earthward propagating plasmoid and flux ropes could have influence upon the dipolarization and injection in inner magnetosphere. The Cluster observations of earthward propagating plasmoid and flux ropes can be interpreted as strong evidence for multiple X-lines. Our observations are consistent with that multiple plasmoids or flux ropes are formed repeatedly and ejected tailward in the course of geomagnetically active time.     

Scale analysis of pre- and post-midnight ESF bubbles at storm time and quiet time

This paper adopts a scale analysis technique to investigate the properties of intermediate-scale plasma structures observed by ROCSAT-1 in the equatorial F-region. A procedure of scale analysis that is developed via the empirical mode decomposition (EMD) method of Hilbert–Huang transform (HHT) technique allows the mutually correlated components in velocity, density and relative density gradient to be identified and extracted. Comparing the three parameters, good match in wave form is found for density and velocity in the scales between kilometers and hundred meters (few kilometers to 300 m). It implies that there are electric fields proportional to density fluctuation −δn/n in the form similar to what is expected for the generalized Rayleigh–Taylor instability. We find that such a one-to-one match holds for various pre- and post-midnight ESF bubbles during quiet and storm times. It, therefore, means that spatial structures of electric field in the intermediate-scale (300 m to few kilometers) correlates to the density structures in a manner of δE ∝  −δn/n that is not necessarily depending on the driving mechanism of ESF bubbles, although it is known that ESF bubbles can be driven by different mechanisms under different space weather conditions. In smaller scales (300–50 m), fluctuation patterns of density and velocity do not correlate to each other any more, the good match is then found in the density gradient ∇_xn/n and velocity. It is known as the manifestation of the Boltzmann relation. We note that the GRT instability related relationship δV_z ∝  −δn/n for irregularities in scale of kilometers holds only for ESF bubbles that occur within ±5 dip latitude, while the Boltzmann relation (δV_z proportional to ∇_xn/n) holds for small-scale irregularities without such a limitation.