Geosynchronous magnetic field response to the large and fast solar wind dynamic pressure change

We present a comparison of large and sharp solar wind dynamic pressure changes, observed by several spacecraft, with fast disturbances in the magnetospheric magnetic field measured by the GOES-8, 9 and 10 geosynchronous satellites. Almost 400 solar wind pressure changes in the period 1996–2003 were selected for this study. Using the large statistics we confirmed that increases (decreases) in the dynamic pressure always results in increases (decreases) in the magnitude of geosynchronous B_z component. The amplitude of the geosynchronous B_z response strongly depends on the location of observer relative to the noon meridian, from the value of solar wind pressure before the disturbance arriving and firstly – from the amplitude of the pressure change.     

地球磁场对太阳风扰动响应的研究

分别对行星际激波、太阳风动压增大事件和减小事件的地球磁场响应进行了比较. 分析结果表明, 同步轨道磁场对太阳风扰动在向阳面产生较强的正响应, 在背阳面 响应较弱且有时会出现负响应, 地磁指数SYM-H对太阳风扰动的响应为正响应. 同时还得出, 向阳侧同步轨道磁场响应幅度d Bz与地磁指数响应幅度d SYM-H、上下游动压均方差均具有较好的相关性. 地磁指数响应幅度与同步轨道磁场响应幅度相关关系在激波和动压增大事件中具有一致性, 动压减小事件出 现明显差异, 这说明激波和动压增大事件在影响地球磁场方面具有某种共性.      

Prediction of geosynchronous electron fluxes using an artificial neural network driven by solar wind parameters

There are hundreds of satellites operating at the geosynchronous (GEO) orbit where relativistic electrons can cause severe damage. Thus, predicting relativistic electron fluxes is significant for spacecraft safety. In this study, using GOES satellite data during 2011–2020, we propose two neural network models with two hidden layers to predict geosynchronous relativistic electron fluxes at two energy channels (>0.8 MeV and > 2 MeV). The number of input neurons of the two channels (>0.8 MeV and > 2 MeV) are determined to be 36 and 44, respectively. The > 0.8 MeV model has 22 and 9 neurons in the hidden layers, while the > 2 MeV model has 25 and 15 neurons in the hidden layers. The input parameters include the north–south component of the interplanetary magnetic field, solar wind speed, solar wind dynamic pressure and solar wind proton density. Through the analysis of different time delays, we determine that the optimal time delays of two energy channels (>0.8 MeV and > 2 MeV) are 8 days and 10 days, respectively. The training set and validation set (Jan 2011-Dec 2018) are divided by the 10-fold cross-validation method, and the remaining data (Jan 2019-Feb 2020) is used to analyze the model performance as a test set. The prediction results of both energy channels show good agreement with satellite observations indicated by low RMSE (～0.3 cm^-2sr^-1s^?1), high PE (～0.8) and CC (～0.9). These results suggest that only using solar wind parameters is capable of obtaining reasonable predictions of geosynchronous relativistic electron fluxes.     

Solar wind interaction with lunar crustal magnetic anomalies

Using Lunar Prospector data, we review the magnetic field and electron signatures of solar wind interaction with lunar crustal magnetic sources. Magnetic field amplifications, too large to represent direct measurements of crustal fields, appear in the solar wind over strong crustal sources, with the chance of observing these amplifications depending on upstream solar wind parameters. We often observe increases in low-energy (?100 eV) electron energy fluxes simultaneously with large magnetic field amplifications, consistent with an increase in plasma density across a shock surface. We also often observe low frequency wave activity in the magnetic field data (both broadband turbulence and monochromatic waves), often associated with electron energization, sometimes up to keV energies. Electron energization appears to be correlated more closely with wave activity than with magnetic amplifications. Detailed studies of the interaction region will be necessary in order to understand the physics of the Moon–solar wind interaction. At present, the Moon represents the only natural laboratory available to us to study solar wind interaction with small-scale crustal magnetic fields, though simulation results and theoretical work can also help us understand the physical processes at work.     

磁鞘等离子体密度统计建模

利用理想磁流体LFM模型的模拟数据,基于非参数统计方法对2004年11月14日03:00UT-07:00UT磁暴恢复相期间磁鞘等离子体平均密度进行建模.分析磁鞘等离子体平均密度与上游太阳风参数、行星际磁场参数及地磁扰动参数的统计关系,建立基于数据降维的经验模型.结果表明,电离层扰动强度因子、太阳风-磁层耦合强度因子和日地空间因果链耦合强度因子是影响磁鞘等离子体平均密度的三个主要方面.磁暴恢复相期间电离层上行离子是磁层环电流和磁尾等离子体的重要离子来源.建模分析过程表明,利用经验模型对空间物理过程开展建模,数据的严重多重共线性通常会导致模型的精度较差.而利用SIR和LPR建立的磁鞘等离子体平均密度随相关参数变化的经验模型可以有效解决该问题,具有较好的预测精度,统计特征显著.      

Interplanetary shocks and sudden impulses during solar maximum (2000) and solar minimum (1995–1996)

In this work a study is performed on the correlation between fast forward interplanetary shock parameters at 1 Astronomical Unit and sudden impulse (SI) amplitudes in the H-component of the geomagnetic field, for periods of solar activity maximum (year 2000) and minimum (year 1995–1996). Solar wind temperature, density and speed, and total magnetic field, were taken to calculate the static pressures (thermal and magnetic) both in the upstream and downstream sides of the shocks. The variations of the solar wind parameters and pressures were then correlated with SI amplitudes. The solar wind speed variations presented good correlations with sudden impulses, with correlation coefficients larger than 0.70 both in solar maximum and solar minimum, whereas the solar wind density presented very low correlation. The parameter better correlated with SI was the square root dynamic pressure variation, showing a larger correlation during solar maximum (r = 0.82) than during solar minimum (r = 0.77). The correlations of SI with square root thermal and magnetic pressure were smaller than with the dynamic pressure, but they also present a good correlation, with r > 0.70 during both solar maximum and minimum. Multiple linear correlation analysis of SI in terms of the three pressure terms have shown that 78% and 85% of the variance in SI during solar maximum and minimum, respectively, are explained by the three pressure variations. Average sudden impulse amplitude was 25 nT during solar maximum and 21 nT during solar minimum, while average square root dynamic pressure variation is 1.20 and 0.86 nPa^1/2 during solar maximum and minimum, respectively. Thus on average, fast forward interplanetary shocks are 33% stronger during solar maximum than during solar minimum, and the magnetospheric SI response has amplitude 20% higher during solar maximum than during solar minimum. A comparison with theoretical predictions (Tsyganenko’s model corrected by Earth’s induced currents) of the coefficient of sudden impulse change with solar wind dynamic pressure variation showed excellent agreement, with values around 17 nT/nPa^1/2.     

A polytropic model for the solar wind

The solar wind fills the heliosphere and is the background medium in which coronal mass ejections propagate. A realistic modelling of the solar wind is therefore essential for space weather research and for reliable predictions. Although the solar wind is highly anisotropic, magnetohydrodynamic (MHD) models are able to reproduce the global, average solar wind characteristics rather well. The modern computer power makes it possible to perform full three dimensional (3D) simulations in domains extending beyond the Earth’s orbit, to include observationally driven boundary conditions, and to implement even more realistic physics in the equations. In general, MHD models for the solar wind often make use of additional source and sink terms in order to mimic the observed solar wind parameters and/or they hide the not-explicitly modelled physical processes in a reduced or variable adiabatic index. Even the models that try to take as much as possible physics into account, still need additional source terms and fine tuning of the parameters in order to produce realistic results. In this paper we present a new and simple polytropic model for the solar wind, incorporating data from the ACE spacecraft to set the model parameters. This approach allows to reproduce the different types of solar wind, where the simulated plasma variables are in good correspondence with the observed solar wind plasma near 1 AU.     

Invariant Modulation of IMF Clock Angle on the Solar Wind Energy Input into the Magnetosphere

By use of the global PPMLR Magnetohydrodynamics (MHD) model, a serial of quasisteady- state numerical simulations were conducted to examine the modulation property of the interplanetary magnetic field clock angle θ on the solar wind energy input into the magnetosphere. All the simulations can be divided into seven groups according to different criteria of solar wind conditions. For each group, 37 numerical examples are analyzed, with the clock angle varying from 0° to 360° with an interval of 10°, keeping the other solar wind parameters (such as the solar wind number density, velocity, and the magnetic field magnitude) unchanged. As expected, the solar wind energy input into the magnetosphere is modulated by the IMF clock angle. The axisymmetrical bell-shaped curve peaks at the clock angle of 180°. However, the modulation effect remains invariant with varying other solar wind conditions. The function form of such an invariant modulation is found to be sin(θ/2)2.70 + 0.25.      

Hybrid simulations of the O ion escape from Venus: Influence of the solar wind density and the IMF x component

As an initial effort to study the evolution of the Venus atmosphere, the influence of the solar wind density and the interplanetary magnetic field (IMF) x component (the x-axis points from Venus towards the Sun) on the O⁺ ion escape rate from Venus is investigated using a three-dimensional quasi-neutral hybrid (HYB-Venus) model. The HYB-Venus model is first applied to a case of the high-density (100 cm⁻³) solar wind interaction with Venus selected from the Pioneer Venus Orbiter observations to demonstrate its capability for the study. Two sets of simulations with a wide range of solar wind densities and different IMF x components are then performed. It is found that the O⁺ ion escape rate increases with increasing solar wind density. The O⁺ ion escape rate saturates when the solar wind density becomes high (above 100 cm⁻³). The results also suggest that the IMF x component enhances the O⁺ ion escape rate, given a fixed IMF component perpendicular to the x-axis. Finally, the results imply a higher ion loss rate for early-Venus, when solar conditions were dramatically different.     

Dynamic processes in the outer heliosphere: voyager observations and models

We review recent Voyager 2 observations in the vicinity of 70 AU. The character of the solar wind plasma data between 2002 and 2003 changed to a regime in which the speed, density and magnetic field magnitude are positively correlated. The average speed of the solar wind at Voyager 2 increased between early 2003 and mid-2004, which we attribute to a return of fast coronal hole flow. We use solar wind data at Earth as input to numerical models which include the effect of pickup ions to model the radial evolution of the solar wind. The model reproduces the basic features of the observations. As a specific example, we investigate the propagation of the Halloween (Oct.–Nov.), 2003 storms in the outer heliosphere. The model predictions are in reasonable agreement with Voyager 2 observations.     

Solar wind structure in the outer heliosphere

A solar wind parcel evolves as it moves outward, interacting with the solar wind plasma ahead of and behind it and with the interstellar neutrals. This structure varies over a solar cycle as the latitudinal speed profile and current sheet tilt change. We model the evolution of the solar wind with distance, using inner heliosphere data to predict plasma parameters at Voyager. The shocks which pass Voyager 2 often have different structure than expected; changes in the plasma and/or magnetic field do not always occur simultaneously. We use the recent latitudinal alignment of Ulysses and Voyager 2 to determine the solar wind slowdown due to interstellar neutrals at 80 AU and estimate the interstellar neutral density. We use Voyager data to predict the termination shock motion and location as a function of time.     

Fluctuations of cosmic rays and IMF in the vicinity of interplanetary shocks

Fluctuations of cosmic rays and interplanetary magnetic field upstream of interplanetary shocks are studied using data of ground-based polar neutron monitors as well as measurements of energetic particles and solar wind plasma parameters aboard the ACE spacecraft. It is shown that coherent cosmic ray fluctuations in the energy range from 10 keV to 1 GeV are often observed at the Earth’s orbit before the arrival of interplanetary shocks. This corresponds to an increase of solar wind turbulence level by more than the order of magnitude upstream of the shock. We suggest a scenario where the cosmic ray fluctuation spectrum is modulated by fast magnetosonic waves generated by flux of low-energy cosmic rays which are reflected and/or accelerated by an interplanetary shock.     

Hybrid modelling the Pioneer Venus Orbiter magnetic field observations

This study presents comparisons between the Pioneer Venus Orbiter (PVO) magnetometer (OMAG) observations and the HYB-Venus hybrid simulation code. The comparisons are made near periapsides of four PVO orbits using the full resolution PVO/OMAG data. Also, the statistics of the solar wind and interplanetary magnetic field (IMF) conditions at Venus are studied using the PVO interplanetary dataset. The statistics include the histograms and the probability density maps of the selected upstream parameters. The confidence intervals derived from the upstream statistics demonstrate the nominal simulation input parameter space. Moreover, the probability density maps give the dependencies between the upstream parameters. The comparisons between the simulation code and the data along the spacecraft trajectory show that the basic, large scale, trends seen in the magnetic field can be understood by the current simulation runs. The discrepancies between the simulation and the data were found to arise at low altitudes close to the planetary ionosphere in the region which cannot be resolved in detail by the grid size of the runs.     

Case study of September 24–26, 1998 magnetic storm

We compute global magnetospheric parameters based upon solar wind data obtained from the WIND spacecraft upstream. Using the paraboloid magnetospheric model, calculations of the dynamic global magnetospheric current systems have been made. The solar wind dynamic pressure, the interplanetary magnetic field, the strength of the tail current, and the ring current control the polar cap and auroral oval size and location during the magnetic storm. The model calculations demonstrate that the polar cap and the auroral oval areas are mainly controlled by the tail current. The substorm onset at 0630 UT on September 25, 1998 happened near the minimum in the main phase field depression. The substorm expansion onset time is also marked by a sudden enhancement in the solar wind dynamic pressure and an enhancement in the tail current. The magnetic signatures of these two effects cancel each other, which explains why the D_st profile shows no strong time variation during the substorm. Evidence for the substorm expansion includes not only the signature in the AL index but also the strong asymmetry of the low latitude magnetic disturbances (substorm positive bay signature). Model calculations were checked by comparison with the GOES 8 and 10 magnetic field measurements.     

一个基于观测的太阳风背景模型

建立由太阳光球磁场和日冕偏振亮度等观测约束的单流体太阳风模型,包括日冕和太阳风的等离子体密度、速度和磁场,温度还有待于以后处理.这里采用高山观测台(HAO)MKⅢ的日冕偏振亮度(pB)在1.36Rs上的观测概图,根据Guhathakurta在1996年发展的日冕电子密度反演模型确定日冕的电子密度分布.同时采用Wilcox太阳观测台(WSO)的光球磁场视向分量的观测概图作为底部边界,根据Zhao等在1994年发展的水平电流-电流片(HCCS)模型得到全球磁场.Phillips在1995年及McComas在2003年分别用Ulysses第一次和第二次跨极飞行的观测发现,归一化到1 AU的太阳风动量流密度除了在10°～30°的纬度范围内略低以外几乎不变.根据这一结论,结合已经得到的密度数据,就可以得到日冕和太阳风的速度.将上面的模型应用于1918卡林顿自转周稳态太阳风的研究,结果与太阳活动极小期的观测基本相符,但是与观测相比较低速高密度区偏大,因此密度模型还有待改进.      

Variations in polar rain flux according to IMF geometries observed by STSAT-1

We examined polar rain flux observed by STSAT-1 in the northern polar cap and compared it with solar wind parameters. We found that the differential energy spectrum of polar rain was similar to that of the solar wind for the energy range 100 eV – 1 keV, although we cannot rule out the possibility of a small amount of acceleration. On the other hand, the low-energy component of the solar wind showed no correlation and, naturally, the solar wind density had only a weak correlation with the polar rain flux. Polar rain flux in the northern hemisphere is most significant for the condition of the interplanetary magnetic field components Bz < 0, Bx < 0, and By > 0, and in this case it correlated well with the magnitude of By and Bz. For other interplanetary magnetic field conditions, the correlation was insignificant. The results are consistent with those reported previously.     

On the long-term variability of the heliosphere–magnetosphere environment

Annual means of measured and reconstructed solar, heliospheric, and magnetospheric parameters are used to infer solar activity signatures at the Hale and Gleissberg cycles timescales. Available open solar flux, modulation strength, cosmic ray flux, total solar irradiance data, reconstructed back to 1700, solar wind parameters (speed and density) and the magnitude of the heliospheric magnetic field at 1 AU, reconstructed back to 1870, as well as the time series of geomagnetic activity indices (aa, IDV, IHV), going back to 1870, have been considered. Simple filtering procedures (successive 11-, 22-, and 88-year running averages and differences between them) and scaling by the standard deviation from the average value for the common interval covered by the data show that the long-discussed variation in the 20th century (a pronounced increase since ∼1900, followed by a depression in the ‘60s and a new, slower, increase) seen in the 11-year averages of parameters such as geomagnetic activity indices and reconstructed heliospheric magnetic field strength, solar wind speed, open solar flux, is a result of the superposition in data of solar activity signatures at Hale and Gleissberg cycles timescales. The Hale and Gleissberg signals were characterized and similarities and differences in the temporal behavior of the analyzed parameters at these timescales are discussed. The similarities in the studied parameters point to a common pacing source, the solar dynamo.     

多种观测数据驱动的三维行星际太阳风MHD模拟

三维磁流体力学（MHD）数值模拟是行星际太阳风研究的重要手段.本文发展了一种由多种观测数据驱动的三维行星际太阳风MHD数值模型.模型的计算区域为0.1AU到1AU附近，使用Lax-Friedrich差分格式在六片网格系统中进行数值求解.边界条件中磁场使用GONG台站观测的光球磁图外推获得，密度通过LASCO观测的白光偏振亮度反演得到，速度根据以上两种观测数据并利用一种基于人工神经网络技术（ANN）的方法得到，温度通过自洽方法根据磁场和密度导出.利用该模型模拟了第2062卡灵顿周（CR2062）时期的行星际太阳风，模拟结果显示出丰富的观测特征，并与OMNI以及Ulysses的实际观测值符合得较好.该模型可用于提供接近真实的行星际太阳风，有助于提高空间天气预报的精度.      

激波在行星际介质中的能量耗散

导出激波下游介质相对上游介质能流通量增量公式,并由HeliosA,B飞船太阳风观测资料得出不同流速太阳风流中各参量随日心距的幂律变化。以此作为背景值分别计算出磁能、内能、动能和总能在不同日心距离处的能量耗散率。结果指出激波后介质以动能增加为主,内能次之,磁能最少;总能耗率在近日处较大,但下降较快。从0.3-1.0AU,不同强度激波总能耗随初始Alfvén激波数A₁₀增大而增大,对A₁₀从2.0-10.0的计算结果与观测值一致。      

Relation between solar wind velocity and properties of its source region

Different kinds of coronal holes are sources of different kind of solar winds. A successful solar wind acceleration model should be able to explain all those solar winds. For the modeling it is important to find a universal relation between the solar wind physical parameters, such as velocity, and coronal physical parameters such as magnetic field energy. To clarify the physical parameters which control the solar wind velocity, we have studied the relation between solar wind velocity and properties of its source region such as photospheric/coronal magnetic field and the size of each coronal hole during the solar minimum. The solar wind velocity structures were derived by using interplanetary scintillation tomography obtained at Solar-Terrestrial Environment Laboratory, Japan. Potential magnetic fields were calculated to identify the source region of the solar wind. HeI 1083 nm absorption line maps obtained at Kitt Peak National Solar Observatory were used to identify coronal holes. As a result, we found a relation during solar minimum between the solar wind velocity and the coronal magnetic condition which is applicable to different kind of solar winds from different kind of coronal holes.