A first assessment of the free energy in nongyrotropic plasmas

Introduction of gyrophase bunching in a particle population can enhance previously existing instabilities or destabilize a thoroughly passive magnetoplasma. The occurrence and intensity of these effects depend on the properties of the background medium and the nature of the gyrophase organization that brings about nongyrotropy. As a first step towards the characterization of the nongyrotropic free energy that drives the stimulated wave activity, this analysis studies the dependence of nongyrotropic growth rates on several gyrophase bunching and magnetoplasma parameters. In general, the unstable modes are reasonably rugged and endure sizeable modifications in their environment before quenching occurs.     

Stability of stationary and time-varying nongyrotropic particle distributions

The ubiquity of nongyrotropic particle populations in space plasmas warrants the study of their characteristics, in particular their stability. The unperturbed nongyrotropic distribution functions in homogeneous media without sources and sinks (closed phase space) must be rotating and time-varying (TNG), whereas consideration of open phase spaces allows for the occurrence of homogeneous and stationary distributions (SNG). The free energy brought about by the introduction of gyrophase organization in a particle population can destabilize otherwise thoroughly stable magnetoplasmas (or, a fortiori, enhance pre-existing gyrotropic instabilities) and feed intense wave growth both in TNG and SNG environments: The nongyrotropic (electron or ion) species can originate unstable coupling among the gyrotropic characteristic waves. The stability properties of these two types of homogeneous nongyrotropy shall be contrasted for parallel (with respect to the ambient magnetic field) and perpendicular propagation, and their potential role as wave activity sources shall be illustrated resorting to solutions of the appropriate dispersion equations and numerical simulations.     

Observation and interpretation of cometary low frequency waves

The first observations of cometary wave activity were carried out in 1985/1986 by several space missions (ICE, VEGAs 1 and 2, Suisei, Sakigake, Giotto) in the environments of comets Giacobini-Zinner and Halley. The interpretation of thesein situ field (and particle) measurements fostered investigations on (among other topics) wave generation that, leaving aside the inherently nonlinear (but related) problem of the eventual formation of a cometary bow shock wave, explored the free energy available in two specific features of the velocity distributions of the newborn particle populations: their parallel (with respect to the IMF direction) drift in the solar wind frame and perpendicular ring-like organization. Analytical and simulation works looked into the influence of the solar wind and cometary newborn parameters on the instabilities and the ensuing, or associated (as evidenced by wave observations), nonlinear phenomenology. Comprehensive reviews have described the experimental and theoretical results obtained in this cometary wave research until 1992 and identified outstanding problems warranting further attention. Here, only a cursory revisit to the Giacobini-Zinner/Halley era of low frequency wave observation and interpretation shall be made: rather, attention shall be predominantly focussed on the new implications to cometary wave research of the recent Giotto encounter with comet Grigg-Skjellerup on July 10 of 1992. The three visited comets, starting with their gas production rates, had different characteristics that showed up in thein situ observations. Yet, with the important exception of the Grigg-Skjellerup encounter, the interpretation of the wave activity measurements could be made in terms of common basic generation mechanisms adapted to the relevant properties of the appropriate plasma environment. New aspects emerged in the last Giotto cometary mission: the smaller gas production rates yield a scale length for the neutral gas density that is not (much) larger than the gyration distance of a heavy newborn ion (estimated by the product of the solar wind speed and the ion cyclotron period). As a consequence of this inhomogeneity, the velocity distribution of the heavy newborn ions exhibits gyrophase organization, i.e. nongyrotropy. This new source of free energy, albeit briefly mentioned in a few studies preceding the Grigg-Skjellerup mission, was not investigated in the context of the Giacobini-Zinner and Halley encounters. Since the last Giotto observations strongly suggest that nongyrotropy plays a prominent role in wave generation as the comet Gigg-Skjellerup nucleus is approached and its stability characteristics have only seldomly been analyzed, the review shall emphasize the wave generation capabilities of particle populations with gyrophase organization.     

等离子体片边界层内静电离子束流-密度漂移不稳定性

本文研究了地球磁尾等离子体片边界层内由离子束流和等离子体密度梯度联合作用产生的静电不稳定性.模型等离子体由向尾流动的冷离子束流、向地球流动的暖离子束流和背景暖电子组成,等离子体密度是非均匀的,等离子体β（热压强与磁压强之比值）很小,电子等离子体频率与电子退旋频率之比。ω_e/Ω_e》1.结果表明,斜传播的静电快、慢离子束流-密度漂移模能够被激发。     

Electrostatic radiation of plasma in the upper ionosphere in the inhomogeneous geomagnetic field

It has been indicated that the spectrum of electrostatic waves in the ionospheric plasma depends on the geophysical conditions and solar wind parameters. The wave field measurements in the frequency band 0.1–10 MHz in the topside ionosphere were used to analyze the electrostatic instabilities of the plasma electron content (the APEX satellite experiment). A change of the sign of one magnetic field component at the geomagnetic equator can result in the formation of the large-scale irregular plasma structure with a decay of the natural electrostatic oscillations and vortices in unstable plasma. The plasma particle polarization drift from the region of decay of electrostatic oscillations and vortices can cause large plasma density and temperature gradients across the geomagnetic field. New vortices can originate at these gradients. This mechanism of plasma vortex formation and decay can be important for mass and energy convection in the topside ionosphere.     

Nongyrotropic particle distributions in space plasmas

In nonstationary, strong inhomogeneous or open plasmas particle orbits are rather complicated. If the nonstationary time scale is smaller than the gyration period, if the inhomogeneity scale is smaller than the gyration radius, i.e. at magnetic plasma boundaries, or if the plasma has sources and sinks in phase space, then nongyrotropic distribution functions occur. The stability of such plasma configurations is studied in the framework of linear dispersion theory. In an open plasma nongyrotropy drives unstable waves parallel and perpendicular to the background magnetic field, whereas in the gyrotropic limit the plasma is stable. In nonstationary plasmas nongyrotropy drives perpendicular unstable waves only. Temporal modulation couples a seed mode with its side lobes and thus it renders unstable wave growth more difficult. As an example of an inhomogeneous plasma a magnetic halfspace is discussed. In a layer with thickness of the thermal proton gyroradius a nongyrotropic distribution is formed which may excite unstable parallel and perpendicular propagating waves.     

On the dispersion of two coexisting nongyrotropic ion species

Space observations in the solar wind and simulations of high Mach number bow-shocks have detected particle populations with two coexisting nongyrotropic ion species. We investigate the influence of these two sources of free energy on the stability of parallel (with respect to the ambient magnetic field) and perpendicular propagation. For parallel modes, we derive their dispersion equation in a magnetoplasma with protons and alpha particles that may exhibit stationary nongyrotropy (SNG) and discuss the characteristics of its solutions. Kinetic simulations study the behaviour of perpendicular electrostatic (Bernstein-like) waves in a plasma whose ion populations (positrons and fictitious singly-charged particles with twice the electron mass, for the sake of simulation feasability) can be time-varying nongyrotropic (TNG). The results show that the coexistence of two gyrophase bunched species does not significantly enhance the parallel SNG instability already found for media with only one nongyrotropic species, whereas it strongly intensifies the growth of Bernstein-like modes in TNG plasmas.     

Asymmetry effect on the development of instabilities in current sheets

The possibility of the development of tearing and kink instabilities (with symmetric and asymmetric modes) in anisotropic thin current sheets (CSs) is investigated. The profiles of perturbed vector potentials and the instability growth rates have been found for the Vlasov equations in the approximation of a linear perturbation theory. The effect of plasma anisotropy and CS asymmetry on the development of these instabilities in the sheet is investigated.     

Low-frequency anisotropy in fractured and layered media

Computing effective medium properties is very important when upscaling data measured at small scale. In the presence of stratigraphic layering, seismic velocities and anisotropy parameters are scale and frequency dependent. For a porous layer permeated by aligned fractures, wave-induced fluid flow between pores and fractures can also cause significant dispersion in velocities and anisotropy parameters. In this study, we compare the dispersion of anisotropy parameters due to fracturing and layering at low frequencies. We consider a two-layer model consisting of an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce anisotropy parameters dispersion due to fractures (meso-scale) in the sand layer. This intrinsic dispersion is added to anisotropy parameters dispersion induced by layering (macro-scale) at low frequencies. We derive the series coefficients that control the behaviour of anisotropy parameters at low frequencies. We investigate the influences of fracture length and fracture density on fracturing effect, layering effect and combined effect versus frequency and volume fraction of sand layer. Numerical modelling results indicate that the frequency dependence due to layering is not always the dominant effect of the effective properties of the medium. The intrinsic dispersion is not negligible compared with the layering effect while evaluating the frequency-dependent properties of the layered medium.     

Frequency‐dependent anisotropy of porous rocks with aligned fractures

Naturally fractured reservoirs are becoming increasingly important for oil and gas exploration in many areas of the world. Because fractures may control the permeability of a reservoir, it is important to be able to find and characterize fractured zones. In fractured reservoirs, the wave‐induced fluid flow between pores and fractures can cause significant dispersion and attenuation of seismic waves. For waves propagating normal to the fractures, this effect has been quantified in earlier studies. Here we extend normal incidence results to oblique incidence using known expressions for the stiffness tensors in the low‐ and high‐frequency limits. This allows us to quantify frequency‐dependent anisotropy due to the wave‐induced flow between pores and fractures and gives a simple recipe for computing phase velocities and attenuation factors of quasi‐P and SV waves as functions of frequency and angle. These frequency and angle dependencies are concisely expressed through dimensionless velocity anisotropy and attenuation anisotropy parameters. It is found that, although at low frequencies, the medium is close to elliptical (which is to be expected as a dry medium containing a distribution of penny‐shaped cracks is known to be close to elliptical); at high frequencies, the coupling between P‐wave and SV‐wave results in anisotropy due to the non‐vanishing excess tangential compliance.     

ULF/ELF monochromatic oscillations observed by Prognoz-8 and -10 spacecrafts during quasiperpendicular supercritical shock crossings

The present paper analyses electromagnetic turbulence in the frequency range 0.1-75 Hz, associated with the supercritical quasiperpendicular crossings of the Earth’s bow shock recorded by the Prognoz-8 and -10 satellites. The quasimonochromatic waves are identified in the shock transition region. Their frequencies, lying in the range 2–5 Hz upstream from the shock ramp, shift to a value less than 1 Hz in the downstream region. The amplitudes of these narrow emissions are great enough to provide the primary dissipation in the flow of the solar wind plasma. Electromagnetic oscillations with such properties are likely to be generated during non-linear evolution of the shock front, rather than by the instabilities driven by ion and electron drifts. Emissions with frequencies higher than 5 Hz have much smaller amplitudes and may be driven by lower hybrid-like instabilities.     

Plasma inhomogeneities in the topside ionosphere in the region of the geomagnetic equator and wave radiation according to the APEX satellite data

The measurements of the broadband wave radiation in the topside ionosphere in the region of the geomagnetic equator (the APEX satellite experiment) are presented. The region of unstable plasma with increased density was observed in the nightside topside ionosphere. This region could be formed by heating of the ionosphere from below. An asymmetric distribution of the frequency band width and electrostatic radiation intensity relative to the geomagnetic equator was registered. It has been indicated that a substantial effect of the geomagnetic equator on plasma diffusion from the heating region could be related to the generation, propagation, and damping of electrostatic oscillations and large-scale (as compared to the Larmor ion radius) plasma vortices. The anisotropy in the temperature of the plasma electron component can increase in the regions where the transverse electric field of disturbances damps. The intensity of the electromagnetic radiation, caused by the external sources, apparently, of an artificial origin at frequencies higher than the local plasma frequency, decreases to the radiation detection threshold level in the region of increased plasma density.     

非均匀磁等离子体中电磁波和电声波耦合传播理论

我们从磁流体力学基本方程出发,给出了计及外加常数磁场情况下规则非均匀等离子体中电磁波和电声波的耦合传播方程。在选用的电子密度剖面下,求解了该方程组的边界值问题,得到了|P_1|和|H_x|的数值解,进而证明了在此基础上我们得到的电磁波和电声波传播的简单近似耦合方程组是可用的。 最后,我们计算了Debye长度和等离子体波长,指出在等离子体鞘套中等离子体波不能辐射出去的主要原因是在鞘套外边界严重的Landau阻尼造成的,并非在激波线强间断处的多次反射结果。 同时,也对所得到的传播耦合方程的普遍性进行了讨论。     

中、低纬度地区电离层不同层结间的耦合对夜晚F区不规则结构的影响

电离层不规则结构的形成和演化与电离层等离子体不稳定性密切相关 .长期以来 ,在中低纬区扩展 -F现象的研究中 ,没有考虑电离层上下层结之间的相互作用 .本文从理论上全面探讨了在中纬度地区 ,E区可变的Pedersen电导率和Hall电导率 ,与F区可变的Pedersen电导率的共同作用对F区梯度漂移不稳定性的影响 ,导出了存在这种耦合时电离层梯度漂移不稳定性的统一表达式 .这一耦合理论不但解释了实际观测中发现的 ,在某些地方电离层F区顶部的不稳定性发生率要高于F区底部这一现象 ;同时还表明 ,电离层E区与F区的耦合对F区夜晚梯度漂移不稳定性的形成不仅会有阻碍作用 ,同时还使得中、低纬度地区的扰动增长具有了方向的选择性 .该理论的一个重要结论是 ,F区中某一地区能否发展出梯度漂移不稳定性 ,并不完全由当地电离层F区的状态决定 ,同一根磁力线连接的、位于不同纬度地区的E区层结对其发展和演化也会有相当大的影响 .     

Stimulated scattering of high-power radio waves in multi-component collisional plasmas

The nonlinear coupling of intense radio waves with the low-frequency electrostatic perturbations of multi-component collisional plasmas is considered. Assuming the presence of two distinct groups of electrons and singly charged ions, we obtain equations for the radio wave sidebands and plasma slow motions that are driven by the combined effects of the radiation pressure and the differential Joule heating of the electrons. The mode coupling equations are useful for studying various types of stimulated scattering instabilities and envelope soliton formation. The relevance of our investigation to the simultaneous generation of density and temperature fluctuations by high-power radio waves is pointed out.     

中高度极隙区电磁不稳定性

研究了中高度（离地心３－４个地球半径）极隙区上行电子束流和上行氧离子（ｏ^＋）锥引起的沿磁力线传播的电磁不稳定性．采用的物理模型假定：上行电子具有单能束流分布函数，而上行氧离子（ｏ^＋）锥可用单能环－束分布函数来描述．结果表明，左旋和右旋圆偏振的低频电磁模是不稳定的，激发不稳定性的自由能源主要由上行电子束流提供，而上行氧离子（ｏ^＋）锥因自由能太小只影响频率色散关系，上行粒子（电子和氧离子）与背景等离子体密度比的变化对电磁不稳定性有重要影响．这些结果对解释权隙区纬度地面站低频电磁波观测资料和理解极隙区动力学过程是很有益的．     

On possible fluid detection in equivalent transversely isotropic media

We consider a transversely isotropic medium that is long-wave equivalent to a stack of thin, parallel, isotropic layers and is obtained using the Backus average. In such media, we analyse the relations among anisotropy parameters; Thomsen parameters, ε and δ, and a new parameter ϕ. We discuss the last parameter and show its essential properties; it equals 0 in the case of isotropy of equivalent medium and/or constant Lamé coefficient λ in layers. The second property occurs to make ϕ sensitive to variations of λ in thin-bedded sequences. According to Gassmann, in isotropic media the variation of fluid content affects only the Lamé coefficient λ, not μ; thus, the sensitivity to changes of λ is an essential property in the context of possible detection of fluids. We show algebraically and numerically that ϕ is more sensitive to these variations than ε or δ. Nevertheless, each of these parameters is dependent on the changes of μ; to understand this influence, we exhibit comprehensive tables that illustrate the behaviour of anisotropy parameters with respect to specific variations of λ and μ. The changes of μ in layers can be presented by the Thomsen parameter γ that depends on them solely. Hence, knowing the values of elasticity coefficients of equivalent transversely isotropic medium, we may compute ϕ and γ, and based on the aforementioned tables, we predict the expected variation of λ; in this way, we propose a new method of possible fluid detection. Also, we show that the prior approach of possible detection of fluids, proposed by Berryman et al., may be unreliable in specific cases. To establish our results, we use the Monte Carlo method; for the range and chosen variations of Lamé coefficients λ and μ – relevant to sandstones – we generate these coefficients in thin layers and, after the averaging process, we obtain an equivalent transversely isotropic medium. We repeat that process numerous times to get many equivalent transversely isotropic media, and – for each of them  – we compute their anisotropy parameters. We illustrate ϕ, ε and δ in the form of cross-plots that are relevant to the chosen variations of λ and μ. Additionally, we present a table with the computed ranges of anisotropy parameters that correspond to different variations of Lamé coefficients.     

Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method

Seismic wave propagation in transversely isotropic (TI) media is commonly described by a set of coupled partial differential equations, derived from the acoustic approximation. These equations produce pure P‐wave responses in elliptically anisotropic media but generate undesired shear‐wave components for more general TI anisotropy. Furthermore, these equations suffer from instabilities when the anisotropy parameter ε is less than δ. One solution to both problems is to use pure acoustic anisotropic wave equations, which can produce pure P‐waves without any shear‐wave contaminations in both elliptical and anelliptical TI media. In this paper, we propose a new pure acoustic transversely isotropic wave equation, which can be conveniently solved using the pseudospectral method. Like most other pure acoustic anisotropic wave equations, our equation involves complicated pseudo‐differential operators in space which are difficult to handle using the finite difference method. The advantage of our equation is that all of its model parameters are separable from the spatial differential and pseudo‐differential operators; therefore, the pseudospectral method can be directly applied. We use phase velocity analysis to show that our equation, expressed in a summation form, can be properly truncated to achieve the desired accuracy according to anisotropy strength. This flexibility allows us to save computational time by choosing the right number of summation terms for a given model. We use numerical examples to demonstrate that this new pure acoustic wave equation can produce highly accurate results, completely free from shear‐wave artefacts. This equation can be straightforwardly generalized to tilted TI media.     

泥岩、页岩声速各向异性及其影响因素分析

在实验室超声波频率下(纵波主频为700kHz、横波为250kHz)对层理发育的页岩和泥岩的各向异性进行了研究,给出了在干燥和油饱和条件下,样品不同方向上纵、横波速度以及各向异性参数随压力的变化规律. 用X 射线衍射和扫描电镜分析了样品中引起各向异性的原因,指出平行于层理定向排列的粘土矿物和微裂隙是使样品显示出强弹性各向异性的内在原因. 随着压力的增高微裂隙逐渐闭合,样品的各向异性程度减弱. 孔隙流体的存在增强了孔隙（裂隙）的刚度,减弱了各向异性随压力增大而减小的趋势.     

Numerical simulation of the mutual conversion of ordinary and extraordinary waves in a magnetised plasma by artificial ionospheric turbulence

Some effects of normal mode coupling in weakly anisotropic inhomogeneous plasma are analysed on the basis of the numerical solution of transfer equations for the Stokes parameters. The numerical analysis of normal mode conversion on a set of isolated irregularities demonstrates the possibility of an effective polarisation transformation on such structures. It is shown that by the appropriate selection of discrete irregularities of the external magnetic field direction or plasma electron concentration resulting in radio-wave refraction, one can control electromagnetic radiation polarisation characteristics. Analytical expressions for mean Stokes parameters have been obtained for the rare isolated irregularities. The opportunity for the simulation of mutual wave conversion processes in the ionospheric heating experiments is discussed.