太阳风压力系数的研究

利用全球磁流体力学(MHD)的模拟结果,研究了太阳风压力系数与上游太阳风参数和日下点磁层顶张角的相关性.在识别出日下点附近磁层顶位置后,通过拟合得到日下点附近的磁层顶张角.在考虑上游太阳风中的磁压和热压以及磁层顶外侧的太阳风动压的情况下,计算了太阳风压力系数.通过分析行星际磁场不同方向时太阳风动压在日地连线上与磁压和热压的转化关系,详细研究了太阳风参数和日下点磁层顶张角对太阳风压力系数的影响,得到以下相关结论:(1) 在北向行星际磁场较大(B_z≥5 nT)时,磁层顶外侧磁压占主导,南向行星际磁场时磁层顶外侧热压占主导;(2) 太阳风压力系数随着行星际磁场的增大而增大,随着行星际磁场时钟角的增大而减小;并且在行星际磁场大小和其他太阳风条件相同时,北向行星际磁场时的太阳风压力系数要大于南向行星际磁场时的;北向行星际磁场时,太阳风压力系数随着太阳风动压的增大而减小,南向行星际磁场时,太阳风压力系数随着太阳风动压的增大而增大;以上结论是对观测结果的扩展;(3) 最后,我们还发现太阳风压力系数随着日下点磁层顶张角的增大而增大.     

Effects of the interplanetary magnetic field clock angle on the shape of bow shock

Using the global magnetohydrodynamics(MHD) simulation model, we investigated the effects of the interplanetary magnetic field(IMF) clock angle on the shape of bow shock, including its rotational asymmetry and subsolar point. For general northward IMF( z component Bz 0), the rotational symmetry of the bow shock is broken by the effects of fast magnetosonic Mach number(Mms), and the cross-sectional line of the bow shock is an ellipse with the semi-major axis along the direction perpendicular to the IMF. The ratio or D-value between semi-major and semi-minor axis can be used to illustrate the extent of asymmetry of the bow shock. On the basis of the multiple parameters fitting, we obtain the changing relationship of both semi-axes with the clock angle and the distance away from the Earth. For general southward IMF(Bz 0), the cross sectional line of the bow shock is highly asymmetrical under the multiple effects of magnetopause and Mms. The effects of IMF clock angle on subsolar point depend mainly on the changing subsolar point of magnetopause as an obstacle. The distance of subsolar point of bow shock from the Earth increases with the increasing IMF clock angle for Bz 0, and decreases with the increasing IMF clock angle for Bz 0.     

Effect of solar dynamics parameters on the formation of substorm activity

An algorithm for retrieving the AL index dynamics from the parameters of solar-wind plasma and the interplanetary magnetic field (IMF) has been developed. Along with other geoeffective parameters of the solar wind, an integral parameter in the form of the cumulative sum Σ[N*V ²] is used to determine the process of substorm formation. The algorithm is incorporated into a framework developed to retrieve the AL index of an Elman-type artificial neural network (ANN) containing an additional layer of neurons that provides an “internal memory” of the prehistory of the dynamical process to be retrieved. The ANN is trained on data of 70 eight-hour-long time intervals, including the periods of isolated magnetospheric substorms. The efficiency of this approach is demonstrated by numerical neural-network experiments on retrieving the dynamics of the AL index from the of solar wind and IMF parameters during substorms.     

Specific features of magnetic barrier formation near the magnetopause

The numerical three-dimensional MHD model is used to study the formation of the magnetic barrier in the inner part of the magnetosheath near the magnetopause. The set of the quasistationary solutions for several characteristic directions of the interplanetary magnetic field (IMF) has been obtained: for northward and southward IMF, for the direction along the Parker helix (at an angle of 45° with respect to the Sun-Earth line), and for the predominantly radial direction (at an angle of 20° with respect to the Sun-Earth line). The mechanism used to take into account the effect of magnetic reconnection at the magnetopause on a flow in the magnetosheath is introduced in the case of southward IMF. The results of the calculations indicate that the magnetic field absolute value in the magnetic barrier reaches its maximal value when IMF is northward. The introduction of magnetic reconnection at southward IMF can result in an insignificant decrease in the field value. However, the model predicts that a decrease in the magnetic field is much more substantial when the IMF direction is close to radial.     

The response of the high-latitude ionosphere to IMF variations

The convection of plasma in the high-latitude ionosphere is strongly affected by the interplanetary magnetic field (IMF) carried by the solar wind. From numerous statistical studies, it is known that the plasma circulation conforms to patterns that are characteristic of particular IMF states. Following a change in the IMF, the convection responds by reconfiguring into a pattern that is more consistent with the new IMF. Some early studies reported that the convection first begins to change near noon while on the dawn and dusk flanks and on the nightside it remains relatively unaffected for tens of minutes. Work by Ridley et al. (J. Geophys. Res. 103 (1998) 4023–4039) and Ruohoniemi and Greenwald (Geophys. Res. Lett. 25 (1998) 2913–2916) that was based on measurements with more global sets of instruments challenged this view. A debate ensued as to the true nature of the convection response. We follow the arguments of Lockwood and Cowley (J. Geophys. Res. 104 (1999) 4387–4391) and Ridley et al. (J. Geophys. Res. 104 (1999) 4393–4396) by reviewing recent results on the timing of the onset of the convection response to the changed IMF. We discuss the timing problem from the perspectives of observations and modeling. In our view, the onset of the ionospheric response to changed IMF is globally simultaneous on time scales of a few minutes. A physical basis for the rapid communication of effects in the dayside convection to the nightside has been demonstrated in magnetohydrodynamic simulations. We also offer some cautionary notes on the timing of convection changes and the use of global assimilative techniques to study local behavior.     

过渡内蕴模态函数对经验模态分解去噪结果的影响研究及改进算法

经验模态分解(Empirical Mode Decomposition, EMD)是一种具有较大应用潜力的去噪算法.目前,该算法存在的一个较大问题是过渡内蕴模态函数(Intrinsic Mode Function, IMF)中混叠噪声不能有效处理.过渡内蕴模态函数中混叠噪声不易剔除,限制了该算法的应用.本文针对此问题,通过研究过渡IMF的特点,首次提出一种有效去除过渡IMF中混叠噪声的方法.该方法首先对原信号进行一次EMD处理,得到包含过渡IMF的初步去噪结果,并将其与合适的余弦信号结合,改变其包络分布,然后对其结果再次进行EMD处理,仿真实验表明该方法在保留有效信号的同时,可以有效的去除过渡IMF中混叠的噪声,并将该方法用于实际地震资料随机噪声压制,处理效果令人满意.     

EMD方法在前兆资料分析处理中的应用

通过仿真实验,主要研究了EMD方法和算法,探讨EMD方法在前兆资料的应用;通过对乌什水管仪分钟值观测资料的分解,该方法能够突出震前高频异常信号,而抑制同震波效应的影响;对整点值数据分解,提取出的IMF1、IMF2分量为半日波和周日波,剩余分量能够精确刻画趋势变化,具有较强的物理意义;另外,该方法也可以应用在失真数据的分析上.研究结果表明该方法在前兆分析上具有较好的应用前景.     

基于经验模态分解的分数维地震随机噪声衰减方法

经验模态分解算法(EMD)是一种基于有效波和噪声尺度差异进行波场分离的随机噪声压制方法,但由于实际地震数据波场复杂,导致模态混叠较严重,仅凭该方法进行去噪很难达到理想效果.本文基于EMD算法对信号多尺度的分解特性,结合Hausdorff维数约束条件,提出一种用于地震随机噪声衰减的新方法.首先对地震数据进行EMD自适应分解,得到一系列具有不同尺度的、分形自相似性的固有模态分量(IMF);在此基础上,基于有效信号和随机噪声的Hausdorff维数差异,识别混有随机噪声的IMF分量,对该分量进行相关的阈值滤波处理,从而实现有效信号和随机噪声的有效分离.文中从仿真信号试验出发,到模型地震数据和实际地震数据的测试处理,同时与传统的EMD处理结果相对比.结果表明,本文方法对地震随机噪声的衰减有更佳的压制效果.     

Interplanetary magnetic field By control of dayside auroras

The effect of the interplanetary magnetic field (IMF) By component on the dayside auroral oval from Viking UV measurements for March–November 1986 is studied. Observations of dayside auroras from Viking UV images for large positive (15 cases) and negative (22 cases) IMF By (∣By∣>4 nT), suggest that: (1) the intensity of dayside auroras tends to increase for negative IMF By and to decrease for positive By, so that negative IMF By conditions seem preferable for observations of dayside auroras; (2) for negative IMF By, the auroral oval tends to be narrow and continuous throughout the noon meridian without any noon gap or any strong undulation in the auroral distribution. For positive IMF By, a sharp decrease and spreading of auroral activity is frequently observed in the post-noon sector, a strong undulation in the poleward boundary of the auroral oval around noon, and the formation of auroral forms poleward of the oval; and (3) the observed features of dayside auroras are in reasonable agreement with the expected distribution of upward field-aligned currents associated with the IMF By in the noon sector.     

Model of Transport of Large-Scale Magnetic Field Fluctuations in the Solar Wind

We consider a model that couples the magnetic field fluctuations in the heliosphere with random shifts of force line footpoints on the Sun. This model generalizes the Giacalone (2001) model by taking into account the large-scale inhomogeneity of the solar wind velocity. This generalization aims to explain a number of specific features of the distribution of IMF directions, such as the change in the asymmetry of the distribution of IMF directions as a function of heliographic latitude and the solar cycle phase and the correlation of azimuthal angles and inclinations of the IMF; the sign of this correlation changes during the solar magnetic cycle. The simulation results have shown that the gradients of the solar wind speed can actually explain these specific features of the distribution of IMF directions, at least qualitatively.     

亚暴期间高纬黄昏-子夜扇区极光弧增亮与衰减及其与IMF的关系

利用南极中山站极光全天空摄相、地磁、地磁脉动数据和Wind卫星的行星际磁场IMF观测数据，分析了7个亚暴期间高纬黄昏-子夜扇区极光弧的短暂增亮现象.极光弧特征是，短暂增亮随后很快衰减，历时10-20min，基本沿着日-地方向，有明显黄昏方向运动.这些事件大都发生在IMFB_z南转之后，亚暴增长相或膨胀相期间，极光浪涌到达之前10-73min消失.相应的IMFB_x＞0，IMFB_y＜0.这种极光弧和亚暴极光不同，它们与地磁活动及Pi2脉动不相关.这7个极光弧的形态和IMF特征表明，极光弧的增亮很可能由尾瓣重联产生，很快衰减归因于IMFB_z南向条件，而黄昏方向运动受IMFB_y控制.     

Poleward expansion of the westward electrojet depending on the solar wind and IMF parameters

The effect of the interplanetary parameters on the latitudinal position of the substorm westward electrojet is studied in the work. The data from the IMAGE chain of magnetic stations and POLAR and WIND satellites for the period close to the solar activity minimum (1995–1996) and for the period of the solar activity maximum (2000) have been used for this purpose. It has been indicated that the electrojet poleward edge reaches, on average, higher latitudes at a higher solar wind velocity and at a larger (B _s ) IMF southward component. It has been indicated that the average latitude of the westward electrojet center increases with increasing solar wind velocity and decreases with increasing IMF southward component, as a result of which the electrojet center is, specifically, not observed at high geomagnetic latitudes at large values of the IMF southward component.     

Nowcasting convection velocity in the high-latitude ionosphere using statistical models

The Weimer and IZMEM statistical convection models are driven with a time series of interplanetary magnetic field (IMF) measurements made onboard the Wind spacecraft. The model outputs are used to infer the ionospheric convection velocity at Casey, Antarctica (80.8°S geomagnetic latitude), and then compared with measurements of Doppler velocity made using a Digisonde, and measurements of F-region convection implied by a collocated magnetometer. During a single, representative campaign interval, 13–17 February 1996, the Weimer model explained 19% (42%) of the variation in Doppler speed (direction) observed by the Digisonde, and 21% (14%) of the equivalent convection components observed by the magnetometer. This compares with IZMEM which explained 16% (46%) of the variation in Doppler speed (direction) observed by the Digisonde, and 34% (32%) of the equivalent convection components observed by the magnetometer. In general, there was better agreement between convection direction than convection speed. Some of the disagreement was probably due to differences between the IMF measured by Wind located ∼170 R_E upstream in the solar wind and the IMF actually arriving at the magnetopause. However, the results of this study do show that measurements of ionospheric velocity using different experimental techniques need heavy averaging to identify a common component of velocity controlled by the IMF vector. The present time series approach was also used to estimate 16±5 min as the time required for the ionospheric convection to reconfigure in response to IMF changes occurring at the magnetopause.     

Some aspects of modelling the high-latitude ionospheric convection from Cluster/Edi data

Measurements onboard Cluster satellites are briefly described, which form the base for determining the intensity and direction of the electric field in the magnetosphere. The aim of this paper is to describe (1) the methodology of calculating the potential distribution at the ionospheric level and the results of constructing spatiotemporal convection patterns for different orientations of the IMF vector in the GSM YZ plane; (2) derivation of basic convection patterns (BCPs), which allow to deduce the statistical ionospheric convection pattern at high latitudes for any IMF Bz and By values (statistical convection model) using different sets of independent data; (3) the consequences of enlarging the amount of data used for analysis; (4) the results of potential calculations with various orders of the spherical harmonics describing them; (5) determination of the cross-polar cap potential with different IMF sector widths (α from 45° down to 10°); (6) the results of our trials to determine the contribution of the IMF Bx component to the convection pattern.     

Relationship between the IMF azimuthal angle and solar wind velocity

The relationship between the IMF azimuthal angle and plasma velocity has been studied independently for three types of solar wind streams (recurrent and transient high-speed streams and low-speed background wind) based on the interplanetary medium parameters measured in the near-Earth orbits in 1964–1996. The relationships between the IMF azimuthal angle cotangent and plasma velocity are close to linear but strongly differ from one another and from the theoretical relationship for all types of streams. These differences area caused by the magnetic field disturbance on the time scales smaller than a day, and the effect of this disturbance has been studied quantitatively. The effective periods of rotation of the IMF sources on the Sun, depending on the solar cycle phase, have been obtained from the relations between the IMF azimuthal angle cotangent and plasma velocity. During the most part of the solar cycle, the periods of rotation of the IMF sources are close to the period of rotation of the solar equator but abruptly increase to the values typical of the solar circumpolar zones in the years of solar minimums.     

模拟IMF北向且By分量占主导时磁层顶重联

本文基于自己开发的全球三维磁层模型, 模拟研究了IMF(Interplanetary Magnetic Field)北向并且B_y分量较大(时钟角为60°)时磁层顶三维结构及其重联图像. 结果发现, IMF B_y为正时, 在北极隙区附近尾-昏侧存在IMF与地磁场之间稳定持续的重联现象;参与重联的地球磁场既有闭合磁力线也有开放磁力线;IMF在北极隙区与地球闭合磁力重联后一端与南磁极相连的磁力线在尾向运动时还可能与北尾瓣的开放磁力线重联而重新闭合, 这种重联与磁力线循环过程不同于同一条IMF磁力线分别在南北半球与地磁场重联的模型. 南极隙区的重联发生在尾-晨侧, 其动力学过程与北极隙区情形类似. 我们的模拟结果表明, IMF B_y较大时不可能发生IMF同一条磁力线分别在南北极隙区重联的情形, 也不会因此而减少尾瓣的开放磁力线.     

地球弓激波的旋转非对称性

通过对太阳风-磁层-电离层系统的全球MHD模拟,研究地球弓激波相对日地连线的旋转非对称性．模拟限于太阳风速度沿日地连线、地球磁偶极矩和行星际磁场（IMF）与日地连线垂直的简单情况．模拟结果表明,即便对于IMF强度为零的情况,弓激波相对日地连线也不具备旋转对称性质：终端面（晨昏子午面）及其向阳侧的弓激波截线的东西宽度大于南北宽度（约9%~11％）,终端面尾侧的弓激波截线东西宽度小于南北宽度（约8％）．在存在IMF的情况下,弓激波的位形同时受到磁层顶的形状和快磁声波速度各向异性的影响．磁层顶向外扩张并沿IMF方向拉伸,且其扩张和拉伸程度随IMF由北转南而增强．在磁鞘中,垂直于磁场方向的快磁声波速度高于平行方向．因此,磁层顶拉伸方向与快磁声波速度最大方向垂直,它们对弓激波位置的效应恰好相反;弓激波的最终形状取决于何种效应占据主导地位．对于终端面尾侧,快磁声波速度的各向异性起主导作用,弓激波截线沿IMF垂直方向的宽度大于平行方向．对于终端面及其向阳侧,弓激波截线的形状与IMF取向有关：在准北向或晨昏向IMF情况下,弓激波截线沿IMF垂直方向的宽度仍大于平行方向;在准南向IMF情况下,弓激波截线沿IMF垂直方向的宽度小于平行方向的．鉴于弓激波形状同IMF取向之间的密切关系,我们提议以IMF为基准方向,提取弓激波截线的平行半宽度R_b∥和垂直半宽度R_b⊥作为尺度参数．这些尺度参数和通常引入的弓激波截线的东西半宽度y_b和南北半宽度z_b相比,更为合理地表征了弓激波的几何性质．模拟结果表明,在终端面上,y_b/z_b和R_b∥/R_b⊥在IMF各向同性取向下的统计平均值均低于1,与观测得到的结论一致.     

极隙区场向电流对高纬电离层电场和电流体系的影响

用Ｋａｍｉｄｅ－Ｍａｔｓｕｓｈｉｔａ方法，在行星际磁场具有较小的北向分量，且｜Ｂｙ｜＞＞｜Ｂｋ｜时，对磁语和磁扰状态以及Ｂｒ＞０和Ｂｙ＜０等不同情况，分别计算了场向电流引起的电离层电势、电场和电流体系．结果表明，极隙区场向电流的存在使高纬向日面区域的电势发生畸变，当Ｂｙ＞０时，无论是磁扰还是磁静日，极隙区电场具有显著的北向分量；等离子体对流有较大的西向分量；电离层电流为东向电流．当Ｂｙ＜０时，电场和等离子体对流的方向与Ｂｙ＞０时相反；电离层电流在磁抗日有西向分量，但在磁静日没有西向分量．电导率对电场和电流体系的影响十分明显，磁扰极光带电导率增强使电流涡从背阳面向向阳而漂移，与静日相比，磁扰时极隙区场向电流引起的电场畸变更为明显，但电场和电流强度的大小却基本保持不变．     

Does the interplanetary magnetic field and the solar activity contribute to the tropospheric circulation?

A significant tendency is shown for both Etesians and (+, –) sector boundaries of the interplanetary magnetic field (IMF) to occur on the same solar rotation days, during the main period of the Etesians effect (July–August). In addition, the solar activity seems to control the Etesians distribution within the IMF sector structure. In the epoch of maximum there is a significant tendency of Etesians to occur during toward IMF days. In contrast, in the epoch of minimum Etesians occur mainly during away IMF days. Finally, in the epoch of intermediate the Etesians are uniformly distributed in away and toward IMF days. Since these conclusions are statistically significant at high confidence levels, it is fair to assume that IMF and solar activity seem to contribute, to some extent, to the Etesians occurrences, as well as to their distribution within the solar rotation and the IMF sector structure; that is, some solar contribution to the tropospheric circulation is implied.     

Time–frequency characterization of sub-divisional scale seasonal rainfall in India using the Hilbert–Huang transform

Time–frequency characterization is useful in understanding the nonlinear and non-stationary signals of the hydro-climatic time series. The traditional Fourier transform, and wavelet transform approaches have certain limitations in analyzing non-linear and non-stationary hydro-climatic series. This paper presents an effective approach based on the Hilbert–Huang transform to investigate time–frequency characteristics, and the changing patterns of sub-divisional rainfall series in India, and explored the possible association of monsoon seasonal rainfall with different global climate oscillations. The proposed approach integrates the complete ensemble empirical mode decomposition with adaptive noise algorithm and normalized Hilbert transform method for analyzing the spectral characteristics of two principal seasonal rainfall series over four meteorological subdivisions namely Assam-Meghalaya, Kerala, Orissa and Telangana subdivisions in India. The Hilbert spectral analysis revealed the dynamic nature of dominant time scales for two principal seasonal rainfall time series. From the trend analysis of instantaneous amplitudes of multiscale components called intrinsic mode functions (IMFs), it is found that both intra and inter decadal modes are responsible for the changes in seasonal rainfall series of different subdivisions and significant changes are noticed in the amplitudes of inter decadal modes of two seasonal rainfalls in the four subdivisions since 1970s. Further, the study investigated the links between monsoon rainfall with the global climate oscillations such as Quasi Bienniel Oscillation (QBO), El Nino Southern Oscillation (ENSO), Sunspot Number (SN), Atlantic Multidecadal Oscillation (AMO) etc. The study noticed that the multiscale components of rainfall series IMF1, IMF2, IMF3, IMF4 and IMF5 have similar periodic structure of QBO, ENSO, SN, tidal forcing and AMO respectively. As per the seasonal rainfall patterns is concerned, the results of the study indicated that for Assam-Meghalaya subdivision, there is a likelihood of extreme rare events at ~0.2 cycles per year, and both monsoon and pre-monsoon rainfall series have decreasing trends; for Kerala subdivision, extreme events can be expected during monsoon season with shorter periodicity (~2.5 years), and monsoon rainfall has statistically significant decreasing trend and post-monsoon rainfall has a statistically significant increasing trend; and for Orissa subdivision, there are chances of extremes rainfall events in monsoon season and a relatively stable rainfall pattern during post-monsoon period, but both monsoon and post-monsoon rainfall series showed an overall decreasing trend; for Telangana subdivision, there is a likelihood of extreme events during monsoon season with a periodicity of ~4 years, but both monsoon and post-monsoon rainfall series showed increasing trends. The results of correlation analysis of IMF components of monsoon rainfall and five climate indices indicated that the association is expressed well only for low frequency modes with similar evolution of trend components.