电离层中SEE现象的静电粒子模拟

利用网格粒子云（PIC）模型。在一定条件下，对电离层加热实验中静电波的激发和参量不稳定性的产生，作一维静电粒子模拟，模拟结果表明，SEE的特征与加热泵波频率ω0和电子回旋频率nωce之间的比值有关。入射泵波的电场与地磁场之间的夹角也是影响SEE的一个因素。     

An investigation of high frequency auroral plasma emissions with a tubular dipole antenna

A band of enhanced amplitudes which follows a local plasma frequency fn in raw high frequency (HF) noise spectra is usually related to plasma emissions in the upper hybrid band (fn, fu). The enhanced band in question occurs permanently in noise spectra recorded on the Intercosmos-19, APEX and CORONAS satellites in the altitude range of 500 km–3000 km. For moderately magnetized plasma with fn > 2fc (fc – electron gyro frequency), the band occurs below fn determined from the topside sounder and impedance data or from electron beam induced spectra. The simulations of an equivalent circuit composed of a dipole antenna in a cold plasma and its preamplifiers, determined the physical origin of the band as the passive circuit resonance, due to inductive character of the antenna in a frequency band (fc, fu). The resonance spectral content is highly structured due to an inflight variability of the circuit impedances. In this report we analyze the noise and impedance spectra which are the most typical in an auroral zone if fn > fc. We focus attention on determination of local electron plasma density, essential for provisional HF mode classification. We found that the natural plasma emission in the upper hybrid band does not manifest itself as the banded natural emission, which may be used for reliable determination of local plasma frequency in the altitude range of 500–3000 km. The fast magnetosonic mode predominates in the auroral emissions. The broadband and multi banded electromagnetic emissions extending from the fast magnetosonic band well above fn > fc are characteristic for the strong wave activity and are much less frequent.     

Equatorial F-region vertical ion drifts during quiet solar maximum

Median values of ionosonde h′F data acquired at Ibadan (Geographic:7.4°N, 3.9°E, Magnetic: dip 6°S, and magnetic declination, 3°W), Nigeria, West Africa, have been used to determine vertical ion drift (electric field) characteristics in the postsunset ionosphere in the African region during a time of high solar activity (average F_10.7 −208). The database spans from January and December 1958 during the era of International Geophysical Year (IGY) for geomagnetic quiet conditions. Bimonthly averaged diurnal variations patterns are very similar, but differ significantly in magnitude and in the evening reversal times. Also, monthly variations of F-region vertical ion drift reversal times inferred from the time of h′F maximum indicates early reversal during equinoxes and December solstice months except for the month of April. Late reversal is observed during the June solstice months. The equatorial evening prereversal enhancement in vertical ion drift (V_zp) occurs largely near 1900 LT with typical values 20–45 m/s. Comparison of Ibadan ionosonde V_zp with the values of prereversal peak velocity reported for Jicamarca (South America), Kodaikanal (India), and Scherliess and Fejer global model show considerable disparity. The changes of postsunset peak in virtual height of F-layer (h′F_P) with prereversal velocity peak V_zp are anti-correlated. Investigation of solar effects on monthly values of V_zp and h′F_P revealed that these parameters are independent of monthly averaged solar flux intensity during quiet-time sunspot maximum conditions.     

Physical and model interpretation of HF radio propagation on the St. Petersburg–Longyearbyen (Svalbard) path

HF radio wave observations have been carried out with an oblique ionospheric sounding (OIS) method on the radio path from St. Petersburg to Longyearbyen (Svalbard), and experimental ionograms were obtained for December 2001. These ionograms have been analysed to investigate the impact of the main ionospheric trough (MIT) and magnetic disturbances on the signals on this path. The observations during weakly disturbed (K_р = 2) magnetic conditions on 14–15 December 2001 were compared with predictions from ray-tracing through a numerical model of the ionosphere. The ray-tracing computer program synthesizes the OIS ionograms by means of the “shooting method”. This method calculates trajectories of HF radio waves for different values of elevation angle and transmission frequency. There was a variety of calculated trajectories, from which we choose those which reach the receiver, and the selected paths provide a synthesis of the oblique ionograms. To simulate HF radio wave propagation, we apply a three-dimensional distribution of the electron density calculated with the mathematical model of the high-latitude ionosphere developed in the Polar Geophysical Institute (PGI). These numerical simulations permit us to interpret specific peculiarities of the OIS data such as abnormal propagation modes, increased delays of signals, enhanced MOF (maximum observed frequency) values etc. New results of the study are summarised as follows. (1) An unusual feature of the propagation along the path is the change of propagation mechanism during substorms on entering a path midpoint (or 1-hop reflection point) to the MIT. (2) Even weak substorms, having the distinguished intensities, lead to the appearance of different types of irregularities observed by the CUTLASS radar and therefore to the different propagation modes and F2MOF values. (3) The PGI model of the ionosphere was first used for ray-tracing at high latitudes. The model results are basically in a good qualitative agreement with experimental observations. This model provides the satisfactory agreement between the calculated and experimental F2MOF values while not correctly representing the fine structure of the experimental OIS ionograms at night. An agreement between the calculated and experimental data is better for day and evening hours than at night.     

Transition radiation mechanism for broad-band component of stimulated electromagnetic emission

As is well known in the F-region of the ionosphere modified by high power HF radio waves broad-band electromagnetic stimulated emission (SEE) is observed. It was discovered both the beams of superthermal electrons and intensive small-scale irregularities in modified region. These magnetic field aligned irregularities have caviton shape — deep electron density holes. The presence of such irregularities and superthermal electrons create conditions for generation of transition emission. We consider this radiation mechanism for interpretation the SEE broad-band component.     

A linear stability analysis of large-Prandtl-number thermocapillary liquid bridges

A linear stability analysis is applied to determine the onset of oscillatory thermocapillary convection in cylindrical liquid bridges of large Prandtl numbers (4  Pr  50). We focus on the relationships between the critical Reynolds number Re_c, the azimuthal wave number m, the aspect ratio Γ and the Prandtl number Pr. A detailed Re_c–Pr stability diagram is given for liquid bridges with various Γ. In the region of Pr > 1, which has been less studied previously and where Re_c has been usually believed to decrease with the increase of Pr, we found Re_c exhibits an early increase for liquid bridges with Γ around one. From the computed surface temperature gradient, it is concluded that the boundary layers developed at both solid ends of liquid bridges strengthen the stability of basic axisymmetric thermocapillary convection at large Prandtl number, and that the stability property of the basic flow is determined by the “effective” part of liquid bridge.     

Study of correlations between waves and particle fluxes measured on board the DEMETER satellite

The topic of relativistic electron dynamics in the outer radiation belt has received considerable attention for many years. Nevertheless, the problem of understanding the physical phenomenon involved is far from being resolved. In this paper, we use DEMETER observations to examine the variations of the energetic electron fluxes and ELF/VLF wave intensities in the inner magnetosphere during the intense 8 November 2004 magnetic storm. Electron flux spectra and associated wave intensity spectra are analysed throughout the magnetic storm and common characteristics or differences to other storm events are retained. The overall objective of this study is to identify and derive parameters that are relevant for particle flux modelling; the time constant characterizing the persistent decay after particle enhancement was found to be one of these important model parameters.The analysis of the 8 November 2004 event reveals that for L-shell parameter higher than 4, an electron flux dropout is observed during the storm’s main phase for electrons in the energy range 0.1–1 MeV, as has been reported from other measurements. Characteristic wave spectra accompanying this phase are analysed. They show a typical enhancement in the frequency range 0.3–10 kHz at onset for all L-shell values under consideration (2 < L < 5). During the first stage of the recovery phase, the electron fluxes are increased to a level higher than the pre-storm level, whereas the level of wave intensity in the frequency range observed below 300 Hz is at its highest. In the second stage, the particle flux decrease goes hand in hand with a global wave activity decline, the relaxation time of the latter being smaller than the former’s one. In some other cases, long-lasting electron enhancement associated with constant wave activity has been observed during this latter stage. For the above mentioned storm, while at low L values the decay time constants are higher for low energy electrons than for high energy electrons, this order is reversed at high L values. At about L = 3.6 the time constant is independent of electron energy.     

Fast spectroscopy of the 4 November 1978 gamma-ray burst

The Franco-Soviet Signe experiments on Venera 11 and Venera 12 allow a spectral analysis of gamma-ray bursts with a time resolution of 250 ms. Evidence is presented for i) short annihilation flashes of up to 20 photons cm⁻²s⁻¹ and ii) rapid variations of the continuum, from a study of the intense 4 November 1978 event.     

Asymmetric precipitation in a coronal loop as explanation of a singular observed spectrum

Almost 10 years of solar submillimeter observations have shown new aspects of solar activity, such as the presence of rapid solar spikes associated with the launch of coronal mass ejections and an increasing submillimeter spectral component in flares. We analyse the singular microwave–submillimeter spectrum of an M class solar flare on 20 December, 2002. Flux density observations measured by Sun patrol telescopes and the Solar Submillimeter Telescope are used to build the radio spectrum, which is fitted using Ramaty’s code. At submillimeter frequencies the spectrum shows a component different from the microwave classical burst. The fitting is achieved proposing two homogeneous sources of emission. This theoretical fitting is in agreement with differential precipitation through a magnetically asymmetric loop or set of loops. From a coronal magnetic field model we infer an asymmetric magnetic structure at the flare location. The model proposed to quantify the differential precipitation rates due to the asymmetry results in a total precipitation ratio Q₂/Q₁≈10⁴–10⁵, where Q₁(Q₂) represents the total precipitation in the loop foot with the high (low) magnetic field intensity. This ratio agrees with the electron total number ratio of the two sources proposed to fit the radio spectrum.     

The iron Kα-line diagnostics of the rotating black hole metric in Seyfert galaxies

Observations of Seyfert galaxies in X-ray region reveal the wide emissive lines in their spectra, which can arise in inner parts of accretion disks, where the effects of General Relativity (GR) must be taken into account. A spectrum of a solitary emission line of a hot spot in Kerr accretion disk is simulated depending on the radial coordinate r and the angular momentum a = J/M of a black hole, under the assumption of equatorial circular motion of a hot spot. It is shown that the characteristic two-peak line profile with the sharp edges arises at a large distance (about r ≈ (3–10) r_g). The inner regions emit the line, which is observed with one maximum and extremely wide red wing. We present results of simulations for the isothermal and Shakura–Sunayev disks.     

Wave emissions at half electron gyroharmonics in the equatorial plasmasphere region: CLUSTER observations and statistics

Intense (n + 1/2) f_ce emissions are a common phenomenon observed in the terrestrial inner magnetosphere. One of their interests is their possible effect in the pitch angle scattering of plasmasheet keV-electron, leading to diffuse auroras. In this paper, we present CLUSTER’s point of view about this topic, in the equatorial region of the plasmasphere, via a statistical study using 3 years of data. Spectral characteristics of these waves, which represent an important clue concerning their generation mechanism, are obtained using WHISPER data near perigee. Details on the wave spectral signature are shown in an event study, in particular their splitting in fine frequency bands. The orbit configuration of the four spacecraft offers a complete sampling on all MLT sectors. A higher occurrence rate of the emissions in the dawn sector and their confinement to the geomagnetic equator, pointed out in previous studies, are confirmed and described with additional details. The proximity of emission sites, both to the plasmapause layer and to the geomagnetic equator surface, seems to be of great importance in the behaviour of the (n + 1/2) f_ce wave characteristics. Our study indicates for the first time, that both the intensity of (n + 1/2) f_ce emissions, and the number of harmonic bands they cover, are increasing as the observation point is located further away outside from the plasmapause layer. Moreover, a study of the wave intensity in the first harmonic band (near 3/2 f_ce) shows higher amplitude for these emissions than previous published values, these emissions can play a role in the scattering of hot electrons. Finally, geomagnetic activity influence, studied via time series of the D_st index preceding observations, indicates that (n + 1/2) f_ce emission events are observed at CLUSTER position under moderate geomagnetic activity conditions, no specific D_st time variation being required.     

Geomagnetic effects on mid-latitude railways: A statistical study of anomalies in the operation of signaling and train control equipment on the East-Siberian Railway

Intense geomagnetically induced currents (GIC) can hamper rail traffic by disturbing signaling and train control systems. GIC threats have been a concern for technological systems at high-latitude locations due to geomagnetic disturbances driven by substorm expansion electrojet or convection electrojet intensifications. However, other geomagnetic storm processes such as storm sudden commencement (SSC) and geomagnetic pulsations can also cause GIC concerns for technological systems. We present in this paper the first evidence based on statistical data for links between geomagnetic disturbances and faulty operations (anomalies) in the functioning of railway automatics and telemetry. We analyze anomalies of automatic signaling and train control equipment which occurred in 2004 on the East-Siberian Railway (corrected geomagnetic latitude _m = 46–51°N and longitude λ_m = 168–187°E). Our results reveal a seasonal effect in the number of anomalies per train similar to the one observed in geomagnetic activity (Kp, Ap, Dst indices). We also found an increase by a factor of 3 in the total duration of daily anomalies during intense geomagnetic storms (local geomagnetic index specific to Siberian Observatory A_max > 30), with a significant correlation between the daily sum of durations of anomalies with geomagnetic activity. Special attention was paid to failures not related to recognized technical malfunctions. We found that the probability of these failures occurring in geomagnetically disturbed periods was 5–7 times higher than the average anomaly occurrence.     

Multi-component synchrotron self-compton emission from NGC 1275

COS-B gamma-ray data (70–5000 MeV) in the latitude range 10°< |b| <90° are compared with the expected emission from cosmic-ray interactions with interstellar gas. An additional component is found to be necessary to explain the latitude dependence of the emission. Two possible origins for this component are discussed: a gamma-ray halo around the Galaxy and a local emission region.     

大功率高频无线电波对赤道E区双流不稳定性的影响--PIC静电粒子模拟研究

地面入射的大功率高频无线电波(泵波)和电离层等离子体之间的参数相互作用，能够引起静电波的激发，在一定条件下，产生不稳定性．本文用PIC静电粒子模拟方法，研究泵波与赤道电离层E区等离子体的相互作用．研究结果表明，泵波能够控制双流不稳定性的发生，在不同条件下，泵波对双流不稳定性起着稳定与不稳作用，模拟结果定性地与理论研究结果相符合，这为我们对不规则体产生的地面人工控制提供了依据．     

On the possibility to check the magnetosphere’s model by CR: The strong geomagnetic storm in November 2003

We report on a study of cosmic ray cutoff rigidity variations during the strong geomagnetic storm of 18–24 November 2003. We employed the most recent Tsyganenko magnetospheric model to represent the very strong disturbed Magnetosphere. We used this magnetic field for the cosmic ray trajectory calculations to determine the geomagnetic cutoff rigidity throughout this period of severe geomagnetic disturbance. We determine the cutoff rigidity changes during this period by two methods, by trajectory calculations and by the spectrographic global survey method. The values of the change of cutoff rigidities obtained by two different methods are correlated with the D_st and interplanetary magnetic field and plasma parameters and result in correlation coefficients in the range 0.63–0.84 for the various cosmic ray stations. The result of this study indicates that the most significant contributions to the cutoff rigidity changes are due to D_st variation although the influence of solar wind density and B_z and B_y components of IMF variations is significant.     

Propagation and origin of ultra high-energy cosmic rays

Propagation of UHE protons through CMB radiation leaves the imprint on energy spectrum in the form of Greisen–Zatsepin–Kuzmin (GZK) cutoff, bump (pile-up protons) and dip. The dip is a feature in energy range 1 × 10¹⁸–4 × 10¹⁹ eV, caused by electron–positron pair production on CMB photons. Calculated for power-law generation spectrum with index γ_g = 2.7, the shape of the dip is confirmed with high accuracy by data of Akeno-AGASA, HiRes, Yakutsk and Fly’s Eye detectors. The predicted shape of the dip is robust: it is valid for the rectilinear and diffusive propagation, for different discreteness in the source distribution, for local source overdensity, deficit, etc. This property of the dip allows us to use it for energy calibration of the detectors. The energy shift λ for each detector is determined by minimum χ² in comparison of observed and calculated dip. After this energy calibration the absolute fluxes, measured by AGASA, HiRes and Yakutsk detectors remarkably coincide in energy region 1 × 10¹⁸–1 × 10²⁰ eV. Below the characteristic energy E_c ≈ 1 × 10¹⁸ eV the spectrum of the dip flattens for both diffusive and rectilinear propagation, and more steep galactic spectrum becomes dominant at E < E_c. The energy of transition E_tr < E_c approximately coincides with the position of the second knee E_2kn, observed in the cosmic ray spectrum. The dip-induced transition from galactic to extragalactic cosmic rays at the second knee is compared with traditional model of transition at ankle, the feature observed at energy 1 × 10¹⁹ eV.     

A summary of the applications of a weighted average method determining times of solar cycle extrema

This paper is a summary of our recent researches on the applications of a weighted average method determining times of solar cycle extrema in the prediction of solar activity. Some correlation coefficients among the parameters in solar cycle according to this definition are higher than those according to the conventional definition. The descending time is found to be correlated (r = −0.77) with the ascending time 3 cycles earlier. The amplitude of solar cycle is found to be correlated (r = −0.77) with the max–max solar cycle length 2 cycles earlier. The ascending time is found to be correlated (r = −0.72) with the amplitude. A newly defined parameter called effective duration is found to be well correlated (r = 0.86) with the amplitude 5 cycles later. These correlations suggest that earlier cycles should influence later ones. The next (24th) solar cycle is estimated to start in March 2007 ± 7 months, reach its maximum in January 2011 ± 14 months, with a size of 150 ± 22, larger than those from some correlations according to the conventional definition.     

Gamma rays from accreting matter onto collapsed objects

In the spherical accretion onto massive objects, the matter may be heated up to temperatures as high as 10¹² °K. In such a hot plasma, the thermal bremsstrahlung (e-e and e-p) and π° decay from inelastic collisions of protons are the main γ-ray sources. We determined the γ -ray production spectra from the π° decay and from bremsstrahlung for different temperatures. The expected γ-ray spectra were evaluated too in order to fit experimental data. We have fitted COS B data from 3C 273 using a two temperatures plasma model. The best fit is for

Full-size image (1K)

(M₈ is the black hole mass in 10⁸ M) which gives . The hard X-ray measurements do not contradict the bremsstrahlung spectrum.     

Evaluation of gravity-dependent membrane potential shift in Paramecium

It is still debated whether or not gravity can stimulate unicellular organisms. This question may be settled by revealing changes in the membrane potential in a manner depending on the gravitational forces imposed on the cell. We estimated the gravity-dependent membrane potential shift to be about 1 mV G⁻¹ for Paramecium showing gravikinesis at 1–5 G, on the basis of measurements of gravity-induced changes in active propulsion and those of propulsive velocity in solutions, in which the membrane potential has been measured electrophysiologically. The shift in membrane potential to this extent may occur from mechanoreceptive changes in K⁺ or Ca²⁺ conductance by about 1% and might be at the limit of electrophysiological measurement using membrane potential-sensitive dyes. Our measurements of propulsive velocity vs membrane potential also suggested that the reported propulsive force of Paramecium measured in a solution of graded densities with the aid of a video centrifuge microscope at 350 G was 11 times as large as that for −29 mV, i.e., the resting membrane potential at [K⁺]_o = 1 mM and [Ca²⁺]_o = 1 mM, and, by extrapolation, that Paramecium was hyperpolarized to −60 mV by gravity stimulation of 100- G equivalent, the value corrected by considering the reduction of density difference between the interior and exterior of the cell in the graded density solution. The estimated shift of the membrane potential from −29 mV to −60 mV by 100- G equivalent stimulation, i.e., 0.3 mV G⁻¹, could reach the magnitude entirely feasible to be measured more directly.     

Registration of ELF waves in rocket-satellite experiment with plasma injection

Two rocket experiments KOMBI-SAMA with plasma injection at height 100–240 km were performed in August 1987 in the region of Brazilian magnetic anomaly (L = 1.25). The launching time of the rocket was determined so that plasma injection was at the time when satellite COSMOS 1809 passed as close as possible to magnetic tube of injection. Caesium plasma jet was produced during ≥ 300 s by electric plasma generator separated from the payload. By diagnostic instruments on board of the rocket and the satellite were registered energetic particle fluxes and plasma wave activities stimulated by plasma injection. When the satellite passed the geomagnetic tube intersecting the injection region an enhancement of ELF emission at 140 Hz, 450 Hz by 2 times was registered on board the satellite. An enhancement of energetic particles (E > 40 keV) flux by 4–5 times was registered on board the rocket. Observed ELV emission below 100 Hz is interpreted as generation of oblique electromagnetic ion-cyclotron waves due to drift plasma instability at the front of the plasma jet.