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.     

The energy release of substorm and magnetic storm: Its present state and future improvements by KuaFu mission

Considering the KuaFu mission, state of the energy release of substorm and storm is simply presented and it’s improvements by KuaFu mission are investigated. The KuaFu mission will provide us an opportunity to improve our understanding of the energy release during the storm and the substorms. The two KuaFu-B satellites flying in 180° phase-lagged formation in a polar orbit will allow synoptic observations of the auroral oval, central plasma sheet, ring current and other regions. It can monitor the polar region 24/7 continuously. The advantage of the KuaFu mission is to provide the data during all phases of storm and substorm time that can be used to study the global energy release during all phases continuously. The data from auroral imager and other in-situ instruments on board KuaFu-B can be used to study the auroral dynamics and Joule heating during a storm and substorm. The data from the neutral atom imager instrument can be used to study the dynamics and the energy release in the ring current region from sudden commencement to complete storm recovery. Furthermore the data from KuaFu-A, which is around L1 point, can be used to study the interplanetary conditions along with the data from the plasma sheet to study the triggering process and energy release during a substorm. So, KuaFu mission with its continuous time monitoring facilities would enable us to make much progress towards solving the underlying problems.     

超声速磁流体管道流动的数值模拟

应用可同时适用于描述磁流体MHD (MagnetoHydroDynamic)强化超燃冲压发动机磁流体发电器与磁流体加速器的,包含磁流体N-S(Navier-Stokes)方程与磁感应方程的可压缩磁流体方程组对超声速磁流体管道流动进行了数值模拟研究,以验证本项工作所提出的对2组方程进行交替迭代的分裂式算法.给出了算例的计算结果,并与无磁场及恒磁场下相同管道的流动计算结果进行了对比分析.     

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.     

再入飞行器表面磁流体发电装置数值模拟

提出了钝锥型再入飞行器上的表面磁流体发电装置方案,采用低磁雷诺数近似下的磁流体力学模型对其进行了流动与静电场耦合的数值模拟,由此建立了表面磁流体发电装置中的物理图像.数值模拟结果表明,在典型再入条件(飞行高度46 km,速度7 km/s)下,所提出表面磁流体发电装置方案能够实现兆瓦级(电功率1.28MW)能量输出,电能提取导致再入飞行器阻力增大13.7%,对飞行器壁面总热流值的影响并不显著,但发电区域及其下游壁面的热流密度分布发生明显变化,该区域内热流密度峰值发生于电极的前、后缘.     

Mini-magnetosphere: Laboratory experiment,physical model and Hall MHD simulation

Magnetosphere with a size comparable to the ion kinetic scales is investigated by means of laboratory experiment, analytical analysis and Hall MHD simulation. In experiment a specific magnetic field was observed which is non-coplanar to dipole field, does not change sign at dipole moment inversion and could be generated only via the quadratic Hall term. Magnetopause position and plasma stand off distance were found to be profoundly different between the experimental regimes with small and large ion inertia length. In the previous studies of a mini-magnetosphere by kinetic codes such novel features were observed as absence of the bow shock and plasma stopping at the Stoermer particle limit instead of the pressure balance distance. Proposed analytical model explains these features by Hall currents which tend to cancel magnetic field convection by ions. Performed numerical simulation shows a good agreement with experiment and analytical model. It gives detailed spatial structure of the Hall field and reveals that while ions penetrate deep inside mini-magnetosphere electrons overflow around it along magnetopause boundary.     

MHD oscillations in solar and stellar coronae: Current results and perspectives

Wave and oscillatory activity is observed with modern imaging and spectral instruments in the visible light, EUV, X-ray and radio bands in all parts of the solar corona. Magnetohydrodynamic (MHD) wave theory gives satisfactory interpretation of these phenomena in terms of MHD modes of coronal structures. The paper reviews the current trends in the observational study of coronal oscillations, recent development of theoretical modelling of MHD wave interaction with plasma structures, and implementation of the theoretical results for the mode identification. Also the use of MHD waves for remote diagnostics of coronal plasmas is discussed. In particular, the applicability of this method to the estimation of the coronal magnetic field is demonstrated.     

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.     

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.     

Refinement of global ionospheric coefficients for GNSS applications: Methodology and results

We developed the methodology for the optimal estimation of global ionospheric coefficients of the current Global Navigation Satellite Systems (GNSSs), including the eight- and ten-parameter Klobuchar-like as well as NeQuick models. The ionospheric coefficients of those correction models are calculated from two sets of globally distributed tracking stations of the International GNSS Services (IGS). Performance of the re-estimated Klobuchar-like and NeQuick coefficients are validated during 2002–2014 over the continental and oceanic areas, respectively. Over the continental areas, GPS TECs derived from 40 ground GPS receivers are selected as reference. The eight-, ten-parameter Klobuchar-like and NeQuick models can mitigate the ionospheric delay by 65.8, 67.3 and 75.0%, respectively. Over the global oceans, the independent TECs derived from Jason-1&2 altimeters are used as reference. The re-estimated ionospheric correction models can mitigate 56.1–66.7% of the delay errors. Compared to the original GPS Ionospheric Correction Algorithm (ICA), performance of those eight-, ten-parameter Klobuchar-like and NeQuick models has improved 3.4, 5.9 and 13.4% during the whole test period, respectively. The methodology developed here takes the advantage of high-quality ionospheric TECs derived from the global network of GNSS receivers. The re-estimated ionospheric coefficients can be used as precise ionospheric products to monitor and assess GNSS broadcast ionospheric parameters and to improve the performance of various single-frequency GNSS applications.     

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.     

强磁暴期间TIEGCM模式与CHAMP卫星热层大气密度的比较分析

利用NCAR-TIEGCM模式计算了2003年11月20—21日强磁暴期间410km高度上的大气密度,并与CHAMP/STAR加速度计反演数据进行对比和分析. 结果表明,模式结果能够准确反映磁暴期间大气密度的分布和变化情况,与实测结果在变化趋势和量级上具有较好的一致性,但在精细结构和数值大小上仍存在一定差异. 模式低估了磁暴期间大气密度的增幅,实测大气密度增幅高达250%～400%,而模式结果为100%～125%. 模式结果与实测数据的偏差在高纬地区高于低纬地区,日侧高于夜侧. 通过模式和实测数据的分析发现,磁暴期间大气密度扰动具有日夜侧和南北半球不对称性. 此外,模式能够准确反映磁暴期间大气密度扰动从高纬向低纬的传播以及大气密度对SYM-H指数响应的延迟特性.     

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.     

Three-dimensional global simulation of multiple ICMEs’ interaction and propagation from the Sun to the heliosphere following the 25–28 October 2003 solar events

This study performs simulations of interplanetary coronal mass ejection (ICME) propagation in a realistic three-dimensional (3D) solar wind structure from the Sun to the Earth by using the newly developed hybrid code, HAFv.2+3DMHD. This model combines two simulation codes, Hakamada–Akasofu–Fry code version 2 (HAFv.2) and a fully 3D, time-dependent MHD simulation code. The solar wind structure is simulated out to 0.08 AU (18 Rs) from source surface maps using the HAFv.2 code. The outputs at 0.08 AU are then used to provide inputs for the lower boundary, at that location, of the 3D MHD code to calculate solar wind and its evolution to 1 AU and beyond. A dynamic disturbance, mimicking a particular flare’s energy output, is delivered to this non-uniform structure to model the evolution and interplanetary propagation of ICMEs (including their shocks). We then show the interaction between two ICMEs and the dynamic process during the overtaking of one shock by the other. The results show that both CMEs and heliosphere current sheet/plasma sheet were deformed by interacting with each other.     

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.