太阳风湍流和磁层亚暴的一种机制

太阳风的动量涨落将通过磁层边界在磁尾激发磁流体力学波。快磁声波携带扰动能量传到等离子体片中,发展为激波,或者通过激波的相互作用而耗散能量,使等离子体加热。等离子体片中的随机费米加速机制,使麦克斯韦分布尾巴部分的高能量粒子被加速到更高能。在宁静态时,加热、加速与耗散过程平衡。当太阳风的动量或者其涨落较大时,整个加热和加速过程加剧,更多的高能粒子产生,并从等离子体片中逃逸,形成高速的等离子体流注入近地轨道和极区,表现为磁层亚暴过程。利用这种机制,可以解释地球磁层亚暴的定性特征。     

地球磁层顶湍动重联的数值模拟

用二维磁流体力学数值模拟研究了磁层顶的磁场湍动重联．提出了一个新的磁场湍动重联模型．数值模拟表明，如果磁层顶是一个开放系统并同时存在磁场剪切和流场剪切，当雷诺数和磁雷诺数超过某临界数值时，磁场重联具有很强的湍动特性，可产生许多不同尺度的磁岛和涡旋结构．随着雷诺数和磁雷诺数的增大和减小，大尺度结构能破碎成中小尺度结构，小尺度结构也能合并成大中尺度结构．湍动重联是涡旋诱发重联在一定条件下的过渡．依据本文的模拟结果，我们预言：磁层预可发生准定常重联、瞬时局地重联和湍动重联等多种重联过程；大中小不同尺度的结构都可以存在于磁层顶；湍动重联及其所产生的中小尺度结构在太阳风－磁层的能量、动量和质量耦合过程中可起重要作用．     

《地球物理学报》总目次

0001南极中山站电离层F2层临界频率变化特征…徐…中…华…,…刘…瑞…源,…刘…顺…林…,…张…北…辰,…胡…红…桥…,…綦……欣0009磁层相对论电子通量变化与磁暴/亚暴的关系………………………李柳元,曹晋滨,…周国成0016磁层顶边界层剪切流引起的MHD波模转化………………     

太阳风动量涨落激发磁层亚暴的机制

本文将太阳风涨落传输能量产生磁层亚暴的机制推广到无碰撞等离子体过程。太阳风的涨落在磁层顶激发压缩阿尔文波,并在磁尾的无碰撞等离子体中传播。尾瓣中满足条件β?1,而等离子体片中β≥1,其中β为等离子体压力与磁压之比。这样,快磁声波在尾瓣中几乎不衰减,而在等离子体片中很快衰减,将波动能量耗散在等离子体片中使等离子体加热或者粒子加速。这种机制还表明,磁尾等离子体片中的高能粒子可以由太阳风涨落动能耗散而被加速,不一定是直接源于太阳。     

地球磁层对磁云边界层的大尺度响应分析——个例研究

主要分析了WIND飞船2004年11月9日探测的磁云边界层引起的大尺度地球磁层活动.磁层响应主要包括以下3个方面:(1)磁云边界层内本身持续较强南向磁场驱动了一个强磁暴的主相.(2)由于磁云边界层内部较强南向磁场持续一段时间后发生向北偏转触发了一个典型磁层亚暴.文中详细分析了亚暴膨胀相发生时夜侧磁层各区域的观测现象,包括极光观测、高纬地磁湾扰、地球同步轨道无色散粒子注入现象、Pi2脉动突然增强以及等离子体片偶极化现象等.(3)磁云边界层和前面鞘区组成一个动压增强区,此动压增强区强烈压缩磁层,致使磁层顶进入地球同步轨道以内;当磁云边界层扫过磁层时,位于向阳侧地球同步轨道上的两颗GOES卫星大部分时间位于磁层磁鞘中,以致很长时间内直接暴露在太阳风中.利用Shue(1998)模型计算得到当磁云边界层扫过磁层时磁层顶日下点的位置被压缩至距地心最近距离为5.1RE,磁云边界层的强动压结构以及强间断面决定了磁云边界层对磁层的强压缩效应.强动压结构、多个强间断结构以及持续较长时间的强南向磁场是许多磁云边界层的共性,这里以此磁云边界层事件为例分析了磁云边界层的地球磁层响应.     

惯性区低-β等离子体中的动力学Alfvén孤波

运用双流体MHD方程描述电子和离子的行为, 采用非扰动理论, 把孤子看作经典粒子, 通过推导的赝势(也称Sagdeev势)方程, 结合数值计算来研究惯性区低-β等离子体中动力学Alfvén孤波的特性. 结果表明离子热效应对Alfvén孤波特性的影响不可忽略, 且在惯性区稀疏型孤波和压缩型孤波均存在. 此结论与Freja卫星在极区上空所观测结果吻合很好. Alfvén孤波携带有平行电场, 它对带电粒子的加速有重要作用, 这给极光粒子加速提供了一种可能的物理机制.     

阿尔文涨落与地球磁层亚暴

本文讨论了一种地球磁层的亚暴机制。当行星际磁场有大的南向分量时,磁层的位形可由基本闭式转变为开式。磁鞘中的阿尔文波可以携带超过10~(18)尔格/秒的能流传入磁层尾部,并将能量耗散于等离子体片中。等离子体片中的粒子被加热和加速后,注入近地空间,产生环电流和极区亚暴。计算了剪切流场中阿尔文波的传播过程,以及磁层中阿尔文波的耗散。将本文的结算与[4]中的结果合在一起,可以说明当行星际磁场转向南时,容易发生地球磁层亚暴,但这两者并非一一对应的关系,行星际磁场没有南向分量时也可以发生地球磁层亚暴。     

地球磁层对磁云边界层的大尺度响应分析——个例研究

主要分析了WIND飞船2004年11月9日探测的磁云边界层引起的大尺度地球磁层活动．磁层响应主要包括以下3个方面：（1）磁云边界层内本身持续较强南向磁场驱动了一个强磁暴的主相．（2）由于磁云边界层内部较强南向磁场持续一段时间后发生向北偏转触发了一个典型磁层亚暴．文中详细分析了亚暴膨胀相发生时夜侧磁层各区域的观测现象,包括极光观测、高纬地磁湾扰、地球同步轨道无色散粒子注入现荆、Pi2脉动突然增强以及等离子体片偶极化现象等．（3）磁云边界层和前面鞘区组成一个动压增强区,此动压增强区强烈压缩磁层,致使磁层顶进入地球同步轨道以内;当磁云边界层扫过磁层时,位于向阳侧地球同步轨道上的两颗GOES卫星大部分时间位于磁层磁鞘中,以致很长时间内直接暴露在太阳风中．利用Shue（1998）模型计算得到当磁云边界层扫过磁层时磁层顶日下点的位置被压缩至距地心最近距离为5．1RE,磁云边界层的强动压结构以及强间断面决定了磁云边界层对磁层的强压缩效应．强动压结构、多个强间断结构以及持续较长时间的强南向磁场是许多磁云边界层的共性,这里以此磁云边界层事件为例分析了磁云边界层的地球磁层响应．     

磁层顶的漂移动力学磁声-Alfvén模不稳定性

本文得出了高β等离子体中在速度梯度、密度梯度、磁场梯度和温度梯度联合作用下的一种新的不稳定模:漂移动力学磁声-Alfvèn模,导出了线性色散关系,研究了模的特性并讨论了高β和各种梯度项对该模的影响,估算了垂直扩散系数。计算表明,在地球磁层顶边界区中漂移动力学磁声-Alfvèn模总是不稳定的。随着β增加,增长率起初增大,而后减小,这是此模的一个显著特性;速度梯度是驱动此模的重要自由能源;当扰动磁场为1nT时,垂直扩散系数可达～10⁹m²/s。从而说明漂移动力学磁声-Alfvèn模不稳定性产生的反常输运,有可能维持准定常的磁层对流,并且在边界层的形成中起重要作用。     

K-H不稳定性在多电流片系统磁场重联中的效应

等离子体系统中存在两个或多个电流片时，电流片中发生的不稳定性可能会相互作用．行星际磁场北向时，背阳面碰层顶电流片与磁尾等离子体片之间可能发生相互作用，高纬边界层强烈的流场剪切可能促进磁场重联，产生磁层亚暴．本文运用二维可压缩磁流体模拟研究具有强流场剪切的多个电流片系统中磁场重联的演化．结果表明，Kelvin－Helmholtz不稳定性使多电流片系统的磁场重联过程明显加快；相邻电流片之间的距离越近，两者相互作用越强，重联增长率越大；在三电流片系统中，超Alfven速度强流场导致外侧两个电流片中出现强烈的磁场重联，并引发中心电流片的磁场重联．行星际磁场北向时，也可能发生磁层亚暴．     

非均匀动压冲击和停止后产生的瞬时重联过程

在非均匀动压冲击期间和冲击突然停止,可引起等离子体边界层的局部瞬时重联过程．本文用二维可压缩ＭＨＤ数值模拟方法研究了这两个过程．结果表明：当大尺度的均匀横向流从一侧边冲击边界层时,磁力线不弯曲,也不发生磁场重联,只是边界层被推着向下游运动;当局部的非均匀动压（特别是横向剪切流）冲击边界层时,被冲击的同向磁场区磁力线逐渐弯曲,在弯曲的反磁场区,出现磁岛,然后在电流片区发生磁场重联,且逐渐形成准稳态的“反Ｋ型”重联结构;当横向剪切流冲击停止后,边界层区变为非常不稳定的系统,产生多种流体涡旋和流型,并相应地产生多种类型的磁场重联结构,直到涡旋消失变为湍动状态时,磁场拓扑才逐渐恢复到未扰动状态。我们提出,外力作用的突然停止,可能是驱动重联的一种新机制,并对这种重联过程在磁层物理中可能的应用进行了讨论．     

Magnetosheath Interaction with the High Latitude Magnetopause

We present both statistical and case studies of magnetosheath interaction with the high-latitude magnetopause on the basis of Interball-1 and other ISTP spacecraft data. We discuss those data along with recently published results on the topology of cusp-magnetosheath transition and the roles of nonlinear disturbances in mass and energy transfer across the high-latitude magnetopause. For sunward dipole tilts, a cusp throat is magnetically open for direct interaction with the incident flow that results in the creation of a turbulent boundary layer (TBL) over an indented magnetopause and downstream of the cusp. For antisunward tilts, the cusp throat is closed by a smooth magnetopause; demagnetized ‘plasma balls’ (with scale ∼ few R_E, an occurrence rate of ∼25% and trapped energetic particles) present a major magnetosheath plasma channel just inside the cusp. The flow interacts with the ‘plasma balls’ via reflected waves, which trigger a chaotization of up to 40% of the upstream kinetic energy. These waves propagate upstream of the TBL and initiate amplification of the existing magnetosheath waves and their cascade-like decays during downstream passage throughout the TBL. The most striking feature of the nonlinear interaction is the appearance of magnetosonic jets, accelerated up to an Alfvenic Mach number of 3. The characteristic impulsive local momentum loss is followed by decelerated Alfvenic flows and modulated by the TBL waves; momentum balance is conserved only on time scales of the Alfvenic flows (1/f_A ∼12 min). Wave trains at f_A∼1.3 mHz are capable of synchronizing interactions throughout the outer and inner boundary layers. The sonic/Alfvenic flows, bounded by current sheets, control the TBL spectral shape and result in non-Gaussian statistical characteristics of the disturbances, indicating the fluctuation intermittency. We suggest that the multi-scale TBL processes play at least a comparable role to that of macro-reconnection (remote from or in the cusp) in solar wind energy transformation and population of the magnetosphere by the magnetosheath plasma. Secondary micro-reconnection constitutes a necessary chain at the small-scale (∼ion gyroradius) edge of the TBL cascades. The thick TBL transforms the flow energy, including deceleration and heating of the flow in the open throat, ‘plasma ball’ and the region downstream of the cusp.     

MHD model of magnetosheath flow: comparison with AMPTE/IRM observations on 24 October, 1985

We compare numerical results obtained from a steady-state MHD model of solar wind flow past the terrestrial magnetosphere with documented observations made by the AMPTE/IRM spacecraft on 24 October, 1985, during an inbound crossing of the magnetosheath. Observations indicate that steady conditions prevailed during this about 4 hour-long crossing. The magnetic shear at spacecraft entry into the magnetosphere was 15°. A steady density decrease and a concomitant magnetic field pile-up were observed during the 40 min interval just preceding the magnetopause crossing. In this plasma depletion layer (1) the plasma beta dropped to values below unity; (2) the flow speed tangential to the magnetopause was enhanced; and (3) the local magnetic field and velocity vectors became increasingly more orthogonal to each other as the magnetopause was approached (Phan et al., 1994). We model parameter variations along a spacecraft orbit approximating that of AMPTE/IRM, which was at slightly southern GSE latitudes and about 1.5 h postnoon Local Time. We model the magnetopause as a tangential discontinuity, as suggested by the observations, and take as input solar wind parameters those measured by AMPTE/IRM just prior to its bow shock crossing. We find that computed field and plasma profiles across the magnetosheath and plasma depletion layer match all observations closely. Theoretical predictions on stagnation line flow near this low-shear magnetopause are confirmed by the experimental findings. Our theory does not give, and the data on this pass do not show, any localized density enhancements in the inner magnetosheath region just outside the plasma depletion layer.     

Deformation patterns in a high-viscosity lava flow inferred from the crystal preferred orientation and imbrication structures: an example from Salina (Aeolian Islands,southern Tyrrhenian Sea,Italy)

Shape-preferred orientation and imbrication structures of crystals have been measured on samples representative of the base, centre and top of a highly viscous lava flow on Salina (Aeolian Islands, southern Tyrrhenian Sea). The data allow zones with different deformation patterns to be identified. In the base and top of the flow, deformation leads to the development of discrete preferred orientation and imbrication of the elongate crystals. The sense of shear is right-lateral at the base and left-lateral at the top of the flow. Shear strain can be estimated by the analysis of crystal preferred orientation. Deformation increases from the flow centre to the outer, more viscous boundary layers. Random orientation of crystals in the inner zone supports the presence of plug flow in a pseudoplastic lava. The textural features of the studied lava may be related to different mechanisms (i.e. lateral expansion). We conclude that the observed crystal alignments and imbrication structures may be related to a plug flow moving between two non-deforming walls. The walls are represented by the solidified, broken upper and basal crust of the flow. The low shear strain values calculated in the outer margins of the flow are indicative of the last deformation event. Crystal preferred orientation and imbrication structures may be related to the occurrence of velocity gradients existing between the inner zone of the flow and its solidus or near-solidus outer margins.     

A multi-layered kinematic dynamo model: implications of a stratified upper layer in the Earth's core

It has been suggested that there exists a stably stratified electrically conducting layer at the top of the Earth's outer fluid core and that lateral temperature gradients in the lower mantle is capable of a driving thermal-wind-type flow near the core–mantle boundary. We investigate how such a flow in a stable layer could influence the geomagnetic field and the geodynamo using a very simple two-dimensional kinematic dynamo model in Cartesian geometry. The dynamo has four layers representing the inner core, convecting lower outer core, stable upper core, and insulating mantle. An α² dynamo operates in the convecting outer core and a horizontal shear flow is imposed in the stable layer. Exact dynamo solutions are obtained for a range of parameters, including different conductivities for the stable layer and inner core. This allows us to connect our solutions with known, simpler solutions of a single-layer α² dynamo, and thereby assess the effects of the extra layers. We confirm earlier results that a stable, static layer can enhance dynamo action. We find that shear flows produce dynamo wave solutions with a different spatial structure from the steady α² dynamos solutions. The stable layer controls the behavior of the dynamo system through the interface conditions, providing a new means whereby lateral variations on the boundary can influence the geomagnetic field.     

地球磁层顶的漂移动力学Alfvn波不稳定性及其反常输运效应

本文用回旋动力学研究了磁层顶等离子体低频漂移动力学不稳定性.在β≥1附近发现两支不稳定的漂移动力学Alfvèn模（DKA）.它们可在▽T（∥-▽n）≠0时激发,速度剪切提供主要的自由能源.两支DKA均具有非零的（δE_y,δE_∥）和（δB_x,δB_∥）.在湍动的非线性饱和状态下,δB_x的起伏可导致很强的反常输运,当|δB_x|=1nT时,D_⊥可达到10⁹m²/s的数量级.因此,DKA可能在太阳风-磁层耦合过程中起重要作用.     

有限导电内核的驱动方式对发电机模型的影响

本文以两种绝缘内核的发电机为基准,设置内核电导率与外核相同,选择了以固定速度超速旋转和在外核驱动下发生旋转的两种内核旋转模式,比较分析不同模型间的能量差异、磁场强度、磁雷诺数、磁极翻转频率和四个类地发电机参数.结果表明:对于弱偶极子发电机模型,有限导电内核的引入会对其偶极子强度的相对变化量造成较大影响,最高达103.00%;由外核驱动旋转的有限导电内核模型对于本文其他目标研究量所带来的影响比较小,均小于5%;而固定旋转速度的有限导电内核的模型对磁极翻转频率、赤道对称性和纬向性均存在较明显影响,最大变化量达124.62%.综合本文所选用的发电机模型的特征和数值分析结果,可以发现虽然由外核驱动旋转的有限导电内核模型其转速不可控且存在较大波动,但各物理量变化量与实际内核与外核的能量比更为接近,因此可以推断其驱动机制较自驱动模式更为合理可靠.     

Geomagnetic cutoff rigidities of cosmic rays in a model of the bounded magnetosphere with the ring current

The cosmic ray geomagnetic cutoff rigidities are obtained by analytical calculations within an axisymmetric model of bounded magnetosphere, the magnetic field of which is created by the dipole field of the Earth and by two spheres located beyond the Earth with the currents that flow along the parallels and have a value proportional to the cosine of latitude. The inner sphere models the ring current flowing in the westerly direction; the outer sphere simulates the currents over the magnetopause, which flow in the easterly direction. The analytical results of calculations of variations in the geomagnetic cutoff rigidity for different levels of geomagnetic disturbances are given. The results are compared with the results of analytical calculations within the model of unbounded magnetosphere (when the outer sphere is absent).     

Determining the thickness of the low-latitude boundary layer in the Earth’s magnetosphere

We describe a method for determining the thickness of the low-latitude boundary layer (LLBL) of the Earth’s magnetosphere at the dayside near the equatorial plane based on the data gathered by a single satellite that traverses the layer and measures the plasma velocity. The method may be applied when the position of the magnetopause and the magnetosheath parameters fluctuate. The necessity of taking the presence of outer and inner LLBL regions into account is analyzed. The developed method is tested using the analysis results of two almost simultaneous close traverses of the magnetopause completed by the THEMIS mission satellites that provided relatively precise data on the LLBL thickness. It is shown that the developed method makes it possible to determine the LLBL thickness with an accuracy of ～10%.