一种全天空气辉成像仪数据预处理方法和初步结果

大气重力波是大气中的基本波动形式之一,在中高层大气动力和热力学过程中起着十分重要的作用.全天空气辉成像仪是一种以大气气辉辐射为示踪物,能够有效对大气重力波成像的仪器.本文针对中国科学院国家空间科学中心空间天气学国家重点实验室中高层大气组自主研制的全天空气辉成像仪所观测的数据,提出了一种气辉图像预处理方法,进行平场校正、方位校正、星光去除和坐标映射等数据订正.利用该方法处理2015年5月17日21:00BLT至次日05:00BLT西宁台站（36.6°N,101.7°E）的OH气辉数据,发现一次重力波事件,分析并获得了该重力波的水平波长、观测水平相速度和传播方向（分别为17.72km,47m·s-1,339°）.研究结果表明该方法是可行的.      

Optical studies of the ionospheric irregularities over the brazilian region by nocturnal images of the OI 630 nm emission

This study is an extension of previous statistical studies (Sobral et al., 1990, 1991, 1999) of both the local time and latitude variations of the zonal drift velocities of ionospheric plasma depletions, over the Brazilian low latitude station Cachoeira Paulista — CP (22.54°S, 45.00°W). The past studies were based on OI 630 nm scanning photometer data and the present one is based on digital OI 630 nm airglow images obtained by an all-sky imager system. These data were gathered between October 1998 and October 1999, at CP. The present results show that, in general, the velocities clearly tended to decrease with local time. Such a decrease should be associated with decreasing intensity of the vertical component of the ambient electric field which, in turn can be accounted for by recombination. All zonal drifts obtained for the 18 nights were eastwards. During equinox, the velocities clearly tended to decrease with local time at lower rates as compared with spring and summer. The highest and lowest zonal drift velocities, from all three seasons considered here, were observed to be in the summer ≈180 ms⁻¹ at 21:45 LT, and in the spring ≈25 ms⁻¹ at 03:15 LT, respectively. Ionospheric plasma bubbles were detected out to the maximum extra-tropical geographical latitude of ≈28° S, which was the highest latitude position analyzed in this study.     

Effect of vertical plasma transport on ionospheric F2-region parameters at equatorial latitude

The variability of the F2-layer even during magnetically quiet times are fairly complex owing to the effects of plasma transport. The vertical E × B drift velocities (estimated from simplified electron density continuity equation) were used to investigate the seasonal effects of the vertical ion drifts on the bottomside daytime ionospheric parameters over an equatorial latitude in West Africa, Ibadan, Nigeria (Geographic: 7.4°N, 3.9°E, dip angle: 6°S) using 1 year of ionsonde data during International Geophysical Year (IGY) of 1958, that correspond to a period of high solar activity for quiet conditions. The variation patterns between the changes of the vertical ion drifts and the ionospheric F2-layer parameters, especially; f_oF₂ and h_mF₂ are seen remarkable. On the other hand, we observed strong anti-correlation between vertical drift velocities and h′F in all the seasons. We found no clear trend between N_mF₂ and h_mF₂ variations. The yearly average value of upward daytime drift at 300 km altitude was a little less than the generally reported magnitude of 20 ms⁻¹ for equatorial F-region in published literature, and the largest upward velocity was roughly 32 ms⁻¹. Our results indicate that vertical plasma drifts; ionospheric F2-layer peak height, and the critical frequency of F2-layer appear to be somewhat interconnected.     

Ionospheric responses to the 21 August 2017 great American solar eclipse – A multi-instrument study

Herein, we report on the ionospheric responses to a total solar eclipse that occurred on 21 August 2017 over the US region. Ground-based GPS total electron content (TEC) data along with ground-based measurements (Millstone Hill Observatory (MHO) and digital ionosondes) and space-based measurements (COSMIC radio occultation (RO) technique) allowed us to identify eclipse-associated ionospheric responses. TEC data at ~20°, ~30°, and ~40°N latitudes from the west to east longitudes show not only considerable depression but also wave-like characteristics in TEC both in the path of totality and away from it, exclusively on the day of eclipse. Interestingly, the observed depressions are associated with lesser (higher) magnitudes at stations over which the solar obscuration percentage was meager (significant), a clear indication of bow-wave-like features. The MHO observes a 30% reduction in F₂-layer electron densities between 180 and 220 km on eclipse day. Ionosonde-scaled parameters over Boulder (40.4°N, 100°E) and Austin (30.4°N, 94.4°E) show a significant decrease in critical frequencies while an altitude elevation is seen in the virtual heights of the F-layer only during the eclipse day and that decreases are associated with wave-like signatures, which could be attributed to eclipse-generated waves. The estimated vertical electron density profile from the COSMIC RO-based technique shows a maximum depletion of 40%. Relatively intense and moderate depths of TEC depression, considerable reductions in the F₂-layer electron densities measured by the MHO and COSMIC RO-measured densities at the F₂-layer peak, and elevations in virtual heights and reduction in the critical frequencies measured by ionosondes during the eclipse day could be due to the eclipse-induced dynamical effects such as gravity waves (GWs) and their associated electro-dynamical effects (modification of ionospheric electric fields due to GWs).     

Registration of wave disturbances over Yakutia

The results of investigations of wave processes with periods 2 hours on their influence and on the night sky airglow intensity are given. The observations were carried out by multichannel spectrometer for three seasons of 1985–1988 at the optical testing ground Maimaga (γ = 63°N; λ = 129, 5°E). The synchronous detection of two and sometimes of three emissions of night sky airglow yielded the oppotunity to track a vertical travel of waves and to estimate their parameters. In most cases the waves propagate upward, i.e. the sources of waves were below mesosphere. The estimated vertical velocity change within 0,9-3,3 m/s and vertical wave length - within 18–85 km. A horizontal velocity varies from 83 to 330 m/s. The wave activity (the occurence frequency) and their amplitude in winter is higher than in spring. The estimated energies transfered by waves to the upper atmosphere are in winter 3.8·10⁻³ W/m² and in spring 2.7·10⁻³ W/m².     

Possible evidence for small-scale wave seeding of equatorial plasma bubbles

The small-scale wave-like structure (SSWS) of F region bottomside plasma density was proposed to be an important seeding for equatorial plasma bubble (EPB) generation, and employed in theoretical simulations of EPBs in recent years. The seeding role of SSWS, however, is waiting to be demonstrated by observation. Here we present two cases of SSWS and EPB observed by the Fuke all-sky airglow imager (19.3°N, 109.1°E; dip latitude 14.3°N). For each case, the results show that two large-scale wave-like structures (LSWSs) initially appeared around sunset in the longitude regions separated by 3–4°, but EPB irregularities were only generated in one of the LSWSs where SSWSs were seen riding on LSWS. For the other LSWS, no SSWS and EPB irregularities were seen. Considering that the two LSWSs were situated closely in longitude where the amplitude of pre-reversal enhancement of background eastward electric field should be similar, the observation that EPB was only generated in the longitude with simultaneous LSWS and SSWS could provide supporting evidence for SSWS seeding of EPB.     

Nighttime thermospheric meridional neutral winds inferred from ionospheric h′F and hpF2 data

Nighttime thermospheric meridional winds aligned to the magnetic meridian have been inferred using h′F and hpF2 ionosonde data taken from two equatorial stations, Manaus (2.9°S, 60.0°W, dip latitude 6.0°N) and Palmas (10.17°S, 48.2°W, dip latitude 6.2°S), and one low-latitude station, Sao Jose dos Campos (23.21°S, 45.86°W, dip latitude 17.26°S), during geomagnetic quiet days of August and September, 2002. Using an extension of the ionospheric servo model and a simple formulation of the diffusive vertical drift velocity, the magnetic meridional component of the thermospheric neutral winds is inferred, respectively, at the peak (hpF2) and at the base (h′F) heights of the F region over Sao Jose dos Campos. An approach has been included in the models to derive the effects of the electrodynamic drift over Sao Jose dos Campos from the time derivative of hpF2 and h′F observed at the equatorial stations. The magnetic meridional winds inferred from the two methods, for the months of August and September, are compared with winds calculated using the HWM-90 model and with measurements from Fabry–Perot technique. The results show varying agreements and disagreements. Meridional winds calculated from hpF2 ionospheric data (servo model) may produce errors of about 59 m/s, whereas the method calculated from the F-region base height (h′F) ionospheric data gives errors of about 69 m/s during the occurrence of equatorial spread-F.     

Electric field measurements at a southern mid-latitude station obtained using an HF digital ionosonde

It is important to understand the convection of the inner magnetosphere to fully describe the response of the low- to mid-latitude thermosphere-ionosphere system to geomagnetic storms. Realistic numerical simulations of mid-latitude electric fields suffer from limited knowledge of lower thermospheric winds and ionospheric conductivity on a global scale. Even empirical models of mid-latitude electric fields suffer from the paucity of measurements made by the handful of incoherent scatter radars concentrated in the American-European sector, and the intermittent satellite measurements made in other regions. Thus it would be very useful to show the extent to which Doppler velocity measurements made with the numerous digital ionosondes deployed around the globe can be used to infer F-region electric fields. The monthly average diurnal variation of Doppler velocity measured by a recently commissioned Digisonde at Bundoora (145.1°E, 37.7°S, geographic; 49°S magnetic) is seen to resemble the average diurnal variation of ion drift measured by the incoherent scatter radar at Millstone Hill (71.5°W 42.6°N; 57°N). Moreover, the Bundoora measurements exhibit the nighttime westward perturbation drifts found in Dynamics Explorer-2 ion drift measurements.     

Airglow OH emission height inferred from the OH temperature and meteor trail diffusion coefficient

Simultaneous observations of the airglow OH(6,2) band rotational temperature, TOH, and meteor trail ambipolar diffusion coefficient, D, were carried out at Shigaraki (35°N, 136°E), during PSMOS 2003 Campaign, January 28 to February 8, 2003. The OH emission height was estimated by cross correlation analysis of the TOH and D nocturnal variations. A good correlation between TOH and D was obtained at 85 km of altitude. From the nocturnal variations of TOH and D, it is found that the OH emission peak height varied from 88 km before the midnight to 84 km in the early morning. The height variation could be caused by an atmospheric tidal effect in the emission height.     

Supertyphoon Hagibis action in the ionosphere on 6–13 October 2019: Results from multi-frequency multiple path sounding at oblique incidence

The purpose of the present paper is to describe the observations of the variations in the parameters of HF radio waves propagating through the ionosphere when the action of the super typhoon Hagibis on 6–13 October 2019 occurred. The observations have been made with the Harbin Engineering University (the People's Republic of China) multi-frequency multiple path radio system involving the software-defined technology. The action of the super typhoon has been shown to be accompanied by enhanced atmospheric wave activity acting to generate wave processes with periods of 10 to 120 min. Coupling in the atmosphere–upper-atmosphere–ionosphere system has been confirmed to be carried out with atmospheric gravity waves. The ionosphere underwent the greatest impact on those days when the supertyphoon had maximum energy, on 8, 10, and especially 9 October 2019, and when it was found to be in an ～2,500–3,000-km distance range from the propagation path midpoints. Under the action of wave processes, the height of the reflection region was observed to oscillate within the ±(30–90 km) limits. The amplitude of the quasi-periodic variations in the ionospheric F-region electron density was estimated to be 10–12% for periods of ～20 min, and 30–60% for periods of ～60–120 min. The joint action of the dusk terminator and the supertyphoon has been confirmed to enhance wave activity in the ionosphere. Similar effects for the dawn terminator have not been detected.     

Hemisphere-coupled modeling of nighttime medium-scale traveling ionospheric disturbances

Nighttime medium-scale traveling ionospheric disturbances (MSTIDs), which have tilted frontal structures in the midlatitude ionosphere, are investigated by the midlatitude ionosphere electrodynamics coupling (MIECO) model in this study. It has been proposed that the electrodynamic coupling between the E and F regions plays an important role in generating MSTIDs within a few hours. An intriguing aspect of MSTIDs is that they were simultaneously observed at magnetic conjugate locations in the Northern and Southern Hemispheres. In order to study the hemisphere-coupled electrodynamics, the MIECO model has been upgraded to consist of two simulation domains for both hemispheres in which the electrostatic potential is solved by considering electrodynamics in both hemispheres. The simultaneous occurrence of MSTIDs at the magnetic conjugate stations has clearly been reproduced when the F-region neutral wind satisfies the unstable condition in both hemispheres and a sporadic-E layer is given only at the Northern (summer) Hemisphere. Even if the unstable condition is satisfied in the summer hemisphere, an unfavorable F-region neutral wind in the winter hemisphere largely suppresses the growth of MSTIDs in both hemispheres.     

Characteristics of nighttime E-region over Arecibo: Dependence on solar flux and geomagnetic variations

Electron concentration (Ne) inferred from Incoherent Scatter Radar (ISR) measurements has been used to determine the influence of solar flux and geomagnetic activity in the ionospheric E-region over Arecibo Observatory (AO). The approach is based on the determination of column integrated Ne, referred to as E-region total electron content (ErTEC) between 80 and 150 km altitude regions. The results discussed in this work are for the AO nighttime period. The study reveals higher ErTEC values during the low solar flux periods for all the seasons except for summer period. It is found that the E-region column abundance is higher in equinox periods than in the winter for low solar activity conditions. The column integrated Ne during the post-sunset/pre-sunrise periods always exceeds the midnight minima, independent of season or solar activity. This behavior has been attributed to the variations in the coupling processes from the F-region. The response of ErTEC to the geomagnetic variability is also examined for different solar flux conditions and seasons. During high solar flux periods, changes in Kp cause an ErTEC increase in summer and equinox, while producing a negative storm-like effect during the winter. Variations in ErTEC due to geomagnetic activity during low solar flux periods produce maximum variability in the E-region during equinox periods, while resulting in an increase/decrease in ErTEC before local midnight during the winter/summer periods, respectively.     

Response of OH,O2 and OI5577 airglow emissions to the mesospheric bore in the equatorial region of Brazil

An all-sky CCD imager capable of measuring wave structure in the airglow OH, O₂ and OI (557.7 nm) emissions was operated in the equatorial region at São João do Cariri (Cariri), Brazil (7°S, 36°W), in collaboration with the Instituto Nacional de Pesquisas Espaciais (INPE). Occurrence of mesospheric bore events was studied using the data from September 2000 to September 2002. Sixty-four bore events were detected during the observation period. Most of the bores showed the complementary effects suggested by Dewan and Picard [E.M. Dewan, R.H. Picard, Mesospheric bores. Journal of Geophysical Research 103, 6295–6305, 1998], except in a few cases where the relative variations were inconsistent with this model.     

Investigation of the characteristics of wavelike oscillations of post-sunset equatorial ionospheric irregularity by decomposing fluctuating TEC

Post-sunset ionospheric irregularities are common features of the equatorial ionosphere that affect radio communication and navigation systems; their triggering physical mechanism is not yet fully understood. Atmospheric gravity wave is considered as a seeding mechanism for the occurrence of ionospheric irregularities (Abdu et al., 2009). To understand the effects of atmospheric waves, characteristics of wavelike oscillation from ionospheric total electron content (TEC) fluctuation that can be obtained from superposition of different oscillation modes have been investigated. Decomposing fluctuating TEC into different oscillation modes and investigating oscillation characteristics of each component is also important to get insight about the characteristics of individual atmospheric waves that may cause TEC fluctuation. In this paper we have investigated characteristics of components of fluctuating TEC obtained from SCINDA GPS receiver installed at Bahir Dar, (geographic coordinate, 11.5°N, 37.6° E, and dip latitude of 2.5°N) Ethiopia during April 2012. First Empirical Mode Decomposition (EMD) has been applied to decompose TEC fluctuation into different oscillation modes that are known as Intrinsic Mode Function (IMF). Hilbert-Huang Transform (HHT) and Continuous Wavelet Transform (CWT) have been applied to investigate the characteristics of wave-like oscillations. Applying EMD on fluctuating vTEC corresponding to a GPS satellite, five components are found. Results from HHT and CWT have shown excellent agreement. In addition, it is found out that the median periods of oscillation of those five components are 9, 17, 47, 78, and 118 min. Of these periods, 17 and 47 min respectively are oscillation periods of components of TEC fluctuation with occurrence frequency of 92% and 91% that may be interpreted as the manifestation of two frequently occurring components of atmospheric gravity waves that are likely generated by the motion of solar terminator.     

基于多仪器观测的中国中纬区域MSTID长距离传播事件研究

电离层等离子体不规则结构通常会影响星地卫星的通信、导航及定位等,因此研究不规则体的结构特征和演化过程具有非常重要的科学意义和应用价值。中尺度电离层行进式扰动（MSTID）是一种常发于F层的电离层扰动,其演化过程十分复杂。本文利用伊春和兴隆台站全天空气辉成像仪、Swarm卫星、佳木斯高频雷达以及漠河和十三陵台站数字测高仪观测数据,对2018年10月17日夜间出现在中国东北区域上空的MSTID事件进行分析。该MSTID事件传播时间较长,在气辉观测中持续时间超过4 h（12:02－16:23 UT）,其波长范围为176.3～322.5 km,波速范围为67.0～154.1 m·s–1。研究结果显示,该MSTID可能产生于较高的纬度,自东北向西南往中纬传播,依次经过伊春和兴隆台站的气辉观测区域。      

Fine structure of the ionospheric plasma bubbles observed by the OI 6300 and 5577 airglow images

Simultaneous measurements of the ionospheric airglow OI 630.0 nm and OI 557.7 nm emissions have been carried out by means of an all-sky CCD imager system at Cachoeira Paulista, since October 1998. During a developed phase of plasma depletion (bubble) in the equatorial anomaly region, both emissions show intensity depletions along the geomagnetic North—South direction, and also bifurcation of the bubbles. It is frequently observed that the OI 557.7 image shows more fine structure of the bubble than the OI 630.0. The amplitude of the intensity depletion was also larger for OI557.7 than OI630.0. This might be due to the difference in life time between the O(¹D) and O(¹S) states, which are responsible for the OI 630.0 and OI 557.7 emissions, respectively. The O(¹D) might be affected by thermal relaxation and diffusion processes before the radiative transition.     

The Nonlinear Model of the Response of Airglow to Gravity Waves

In this paper, we develope a time-dependent, nonlinear, photochemical-dynamical 2-D model which is composed of 3 models: dynamical gravity wave model, middle atmospheric photochemical model, and airglow layer photochemical model. We use the model to study the effect of the gravity wave propagation on the airglow layer. The comparison between the effects of the different wavelength gravity wave on the airglow emission distributions is made. When the vertical wavelength of the gravity wave is close to or is shorter than the thickness of the airglow layer, the gravity wave can make complex structure of the airglow layer, such as the double and multi-peak structures of the airglow layers. However, the gravity wave that has long vertical wavelength can make large scale perturbation of the airglow emission distribution.     

Multi-station observation of ionospheric irregularities over South Africa during strong geomagnetic storms

This paper presents results pertaining to the response of the mid-latitude ionosphere to strong geomagnetic storms that occurred from 31 March to 02 April 2001 and 07–09 September 2002. The results are based on (i) Global Positioning Systems (GPSs) derived total electron content (TEC) variations accompanying the storm, (ii) ionosonde measurements of the ionospheric electrodynamic response towards the storms and (iii) effect of storm induced travelling ionospheric disturbances (TIDs) on GPS derived TEC. Ionospheric data comprising of ionospheric TEC obtained from GPS measurements, ionograms, solar wind data obtained from Advanced Composition Explorer (ACE) and magnetic data from ground based magnetometers were used in this study. Storm induced features in vertical TEC (VTEC) have been obtained and compared with the mean VTEC of quiet days. The response of the mid-latitude ionosphere during the two storm periods examined may be characterised in terms of increased or decreased level of VTEC, wave-like structures in VTEC perturbation and sudden enhancement in hmF2 and h′F. The study reveals both positive and negative ionospheric storm effects on the ionosphere over South Africa during the two strong storm conditions. These ionospheric features have been mainly attributed to the travelling ionospheric disturbances (TIDs) as the driving mechanism for the irregularities causing the perturbations observed. TEC perturbations due to the irregularities encountered by the satellites were observed on satellites with pseudo random numbers (PRNs) 15, 17, 18 and 23 between 17:00 and 23:00 UT on 07 September 2002.