武汉上空中层顶区域潮汐的MF雷达观测

利用武汉中频雷达观测数据进行分析,研究2至3月份武汉上空中层顶潮汐结构及其随高度和时间变化的特性.用Lomb-Scargle周期图方法计算的水平风场动态功率谱表明,武汉上空存在持续的周日潮汐,是中层顶区域风场结构的主要成分.周日潮的平均振幅随高度的增加呈先增后减的趋势.大多数情况下,潮汐谱峰对应的频率与定义值有一定的偏移.周日潮水平扰动速度矢量随时问和高度变化的轨迹表明,经向分量的相位比纬向分量的相位超前,潮汐能量向上传播,对应于向下的相位传播速度.计算得出的经向分量和纬向分量的垂直相速度分别为1.10和1.15 km/h.      

中纬度地磁暴期间热层垂直风响应机制的模拟

基于TIMEGCM模型,研究了2005年9月10日中纬度地磁暴期间热层(100～650 km)水平风场变化对垂直风的影响.通过连续性方程诊断分析了暴时引起垂直风场变化的机制,结果表明:250 km以上的垂直风场取决于水平风场的变化,而250 km以下的垂直风场由较高高度的垂直风拉动;在地磁暴初相开始时,经向风场相比纬向...     

Longitudinal MLT wind structure at higher mid-latitudes as seen by meteor radars at central and Eastern Europe (13°E/49°E)

The mid-latitude mesosphere and lower thermosphere (MLT) wind speeds measured by two SKiYMET meteor radars (MRs) at Collm (51°N, 13°E) and Kazan (56°N, 49°E) during 2016–2017 were analyzed to study longitudinal wind structures. The differences between monthly mean prevailing wind speeds and tidal amplitudes were compared with the corresponding differences obtained from TIMED/TIDI satellite winds and gradient wind speeds from the AURA/MLS instrument. It is shown that the MR wind difference between the two sites is statistically significant. The difference of the horizontal prevailing winds can be explained by a superposition of the background zonal flow, which is different at the two latitudes, with stationary planetary waves of different origin. Non-migrating tides contribute significantly to the difference of the semidiurnal tidal winds between the two sites.     

Mesosphere/lower thermosphere wind regime parameters using a newly installed SKiYMET meteor radar at Kazan (56°N, 49°E)

New meteor radar (MR) horizontal wind data obtained during 2015–2018 at Kazan (56°N, 49°E) are presented. The measurements were carried out with a state-of-the-art SKiYMET meteor radar. Monthly mean vertical profiles of zonal and meridional components of the prevailing wind speeds, also amplitudes and phases of the components of diurnal (DT) and semidiurnal tide (SDT) winds are displayed as contour plots for a mean calendar year over the four recent years and compared with distributions of these parameters provided by the previous multiyear (1986–2002) meteor radar (MR) measurements at Kazan and by the recent HWM07 empirical model. The analysis shows that the SKiYMET zonal and meridional prevailing wind speeds are generally in good agreement, sharing the same seasonal features, with the earlier MR seasonal winds. Comparisons with the HWM07 model are not favourable: eastward solstitial cells as modelled are significantly larger, >30?m/s compared to 15–20?m/s. Also, reversal lines are too variable with height, and the positions of modelled cells (positive and negative) are unlike those of either MRs at Kazan or other MLT radars. Both MR systems provide the large SDT amplitudes, approximately 30?m/s and vertical wavelengths, approximately 55?km, for both components at middle latitudes in winter. They also show the well known strong SDT September feature (heights 85–100?km, the vertical wavelength ～55–60?km), and the weak summer SDT for 80–91?km. HWM07 shows unrealistic amplitudes and phases above 90?km by height and month: minimal amplitudes in equinoxes and no September feature.The weak DT of middle to high latitudes provide similar amplitude and phase structures from both MRs, 1986–2002 and 2015–2017: largest amplitudes (10–12 or 8–10?m/s) for the evanescent meridional tide in summer, peaking in late July; weakest (0–2, 2–4?m/s) at 80 to 92–96?km, when the tide is vertically propagating (January, February, November, December) with a vertical wavelength near 40?km. Again, HWM07 differs in amplitude and phase structures: showing peak amplitudes in equinoxes: April, 15?m/s at 88?km; October, 21?m/s at 89?km.Coupling of the MR wind parameters with the ERA5 wind parameters is studied for a case in 2016. It is shown that the prevailing winds and DT amplitudes and phases of both datasets can be simply linked together, but that the ERA5 SDT amplitudes are significantly underestimated at the top model levels of the ERA5 reanalysis project.     

Equatorial M(3000)F2 estimation of F2-layer models peak heights during low solar activity

M(3000)F2 estimation of h_mF2 based on four different formulated models viz: (1) Shimazaki (1955) (2) Bradley and Dudeney (1973), (3) Dudeney (1974) and (4) Bilitza et al. (1979) at an equatorial station in West Africa during low solar activity period (1995) are used to validate its conformity with observed and International Reference Ionosphere (IRI) model. Local time analyses of data from fifteen (15) selected days during the January and July solstices and April and October equinoxes are used. The results obtained show that the M(3000)F2 estimation of h_mF2 from the ionosonde-measured values using the Ionospheric Prediction Service (IPS-42) sounder compared to the observed values which were deduced using an algorithm from scaled virtual heights of quiet day ionograms are highly correlated with Bilitza model. International Reference Ionosphere (IRI 2007) model for the equatorial region also agrees with the formulation developed by Bilitza et al. (1979) for the four different seasons of the year. h_mF2 is highest (425 km) in summer (June solstice) season and lowest (386 km) in autumn (September equinox) season with daytimes peaks occurring at 1100–1200 LT during the solstices and at 1000 LT during the equinoxes respectively. Also, the post-sunset peaks are highest (362 km) at the spring (March equinox) and lowest (308 km) at the summer (June solstice) both occurring between 1800 and 2000 LT.     

The effect of fluctuating ionospheric electric fields on Es-occurrence at cusp and polar cap latitudes

Theory predicts that in the high-latitude southern hemisphere, southwest (SW) electric fields will produce convergent ion flow and thereby create thin sporadic E (Es)-layers at node heights > 120 km, whilst northwest (NW) fields will produce downward ion flow and create thicker Es-layers at heights <110 km. To investigate this theory, Digisonde ionograms (giving the Es-occurrence) and drift measurements (giving electric field estimates) at two Antarctic stations were statistically analyzed. As previously found for the polar cap station Casey (81°S magnetic), more of the Es-traces were associated with SW fields than NW fields. However, new results for the cusp station Zhongshan (73°S) show that fewer Es-layers occur there, and NW fields play a slightly more important role than SW fields, similar to the results found at auroral latitudes in the northern hemisphere. To further our understanding of the occurrence distributions, we study the fluctuating properties of the electric fields at the two stations. It is found that the electric fields at Zhongshan fluctuate more than those at Casey. Thus we suggest that the field fluctuation is also an important consideration helping to explain the differences in the Es-occurrence at the two stations. This suggestion is confirmed by our numerical simulations which show that Es-layers are more effectively formed by steady SW fields than by steady NW fields, and less effectively by fluctuating SW fields than by fluctuating NW fields.     

Horizontal winds in the mesosphere at high latitudes

A total of 146 meteorological rocket flights applying the ‘falling sphere’ technique are used to obtain horizontal winds in the mesosphere at polar latitudes, namely at the Andøya Rocket Range (69°N, 125 flights), at Spitsbergen (78°N, 10 flights), and at Rothera (68°S, 11 January flights only). Nearly all flights took place around noon or midnight, i.e., in the same phase of the semidiurnal tide. Meridional winds at 69°N show a clear diurnal tidal variation which is not observed in the zonal winds. The zonal wind climatology shows a transition from summer to winter conditions with the zero wind line propagating upward from 40 km (end of August) to 80 km (end of September). Zonal winds are smaller at Spitsbergen compared to Andøya which is in line with a common angular velocity at both stations. Meridional winds at noon are of similar magnitude at all three stations and are directed towards the north and south pole, respectively. Horizontal and meridional winds generally agree with empirical models, except for the zonal winds at Antarctica which are similar to the NH, whereas there is a significant SH/NH difference in CIRA-1986.     

Estimation of the elastic Earth parameters (h2, l2) using SLR data

We present results for the global elastic parameters h₂ and l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two satellites LAGEOS 1 and LAGEOS 2 observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analysed using different approaches. The analysis was done separately for the two satellites and approaches to estimate the two elastic parameters independently and together were performed. We do a sequential analysis and study the stability of the estimates as a function of length of the data set used. The adjusted final values for h₂ equal to 0.6151 ± 0.0008 and 0.6152 ± 0.0008, and those for l₂ equal to 0.0886 ± 0.0003 and 0.0881 ± 0.0003 for LAGEOS 1 and LAGEOS 2 tracking data are compared to other independently derived estimates. These parameters and their errors achieve stability at about the 24 and 27 month time interval for h₂ and l₂, respectively.     

Zonal and meridional wind components derived from intercosmos-19 hmF2 measurements

The technique for the derivation of the meridional, V, and zonal, U, components of neutral wind from the longitudinal variations of vertical plasma drift, W, is developed. (Longitudinal variations of W were calculated by means of servo-model from Intercosmos-19 hmF2 data carefully selected for the ±40° invariant and geographical latitudes in the Northern and Southern hemispheres.) The technique is based on expansion of longitudinal variations of W, V, U and parameters of the geomagnetic field into finite Fourier series and on solution of the obtained equations set. The best solution of this problem is obtained by means of the Tikhonov regularization method. The most precise solution is derived for average value of meridional wind, the least precise one - for the longitudinal variations of zonal wind. The comparison with HWM and MWM wind models is carried out. The contributions of the different factors in the longitudinal variations of hmF2 are estimated.     

Structure of the migrating diurnal tide in the Whole Atmosphere Community Climate Model (WACCM)

As part of an ongoing effort to understand the migrating diurnal tide generated by the NCAR Whole Atmosphere Community Climate Model, version 3 (WACCM3), we compare the WACCM3 migrating diurnal tide in the horizontal wind and temperature fields to similar results from the Global Scale Wave Model (GSWM). The WACCM3 diurnal tidal wind fields are also compared to tropical radar measurements at Kauai (22°N, 200.2°E) and Rarotonga (21.3°S, 199.7°E). The large-scale features of the WACCM3 results, such as the global spatial structure and the semiannual amplitude variation are in general agreement with past tidal studies; however, several differences do exist. WACCM3 exhibits a much higher degree of hemispheric asymmetry, lower overall amplitudes around the equinoxes, and peaks which are more confined in latitude when compared with the GSWM. Factors which may contribute to such differences between WACCM3 and GSWM are the solar heating profiles from ozone and water vapor, dissipation, and the zonal mean zonal winds. We find that the internally generated heating in WACCM3 and eddy dissipation values are both smaller than the values specified in the GSWM; the eddy dissipation fields and zonal mean zonal winds of the two models also display measurable differences in spatial structure. Comparisons with radar data show several differences in spatial and seasonal structure. In particular, the diurnal tide zonal winds in WACCM3 above Kauai are considerably larger in amplitude than those observed in the radar data, due to contributions from nonmigrating tidal components including wave numbers eastward 1 through 3, westward 2, and stationary components, which interfere constructively with the migrating component around equinox in WACCM3.     

Four-peak longitudinal distribution of the equatorial plasma bubbles observed in the topside ionosphere: Possible troposphere tide influence

In this paper we consider an idea of the troposphere tide influence on the character of the longitudinal variations in the distribution of the equatorial plasma bubbles (EPBs) observed in the topside ionosphere. For this purpose, the obtained EPB longitudinal patterns were compared with the thermosphere and ionosphere characteristics having the prominent “wave-like” longitudinal structures with wave number 4, which are uniquely associated with the influence of the troposphere DE3 tides. The characteristics of the equatorial mass density anomaly (EMA), equatorial ionization anomaly (EIA), zonal wind and pre-reversal E?×?B drift enhancement (PRE) were used for comparison. The equinox seasons during high solar activity were under consideration. It was obtained that the longitudinal patterns of the EMA and zonal wind show the surprising similarity with the EPB distributions (R???0.8, R???0.72). On the other hand, the resemblance with the ionosphere characteristics (EIA, PRE) is rather faint (R???0.37, R???0.12). It was shown that the thermosphere zonal winds are the most possible transfer mediator of the troposphere DE3 tide influence. The most successful moment for the transfer of the troposphere DE3 tide energy takes place in the beginning of the EPB production, namely, during the seed perturbation development.     

Response of the lower atmosphere to intense geomagnetic storms

In this short paper we examine the possible connection between atmospheric parameters measured at low and middle altitudes and geomagnetic storms occurred in 2000 and 2003. For that, from a chain of stations located near the meridian 60°W we compare the storm time values of temperature and wind speed with their standard deviation 2σ obtained from quiet time values. We observed statistically significant variations at several altitudes during the storm recovery phase and after it, both in neutral wind speed and temperature. The results obtained suggest that atmospheric parameters could be affected by geomagnetic storms.     

基于中国岢岚地区法布里-珀罗干涉仪的风场及行星波观测特征

利用中国岢岚站（38.7°N,111.6°W）法布里-珀罗干涉仪2013年7月至2014年11月的水平风场数据,对87,97,250km风场长期变化和行星波特征进行了研究.通过分析年振荡（AO）和半年振荡（SAO）振幅相位,将午夜风场与HWM07数据对比发现:87km和97km处FPI纬向风变化趋势与HWM07相近,而经向风相位落后于HWM07,从振幅上看,HWM07振幅偏大;250km处风场月变化大,FPI与HWM07差异大,HWM07模式的准确性需进一步考虑太阳活动和行星际磁场的影响.利用Lomb-Scargle功率谱以及最小二乘谐波拟合提取了三个高度的行星波振幅,其特征表明87km和97km处纬向风16日波秋季及冬春季活动强,而6.5日波最强振幅出现在春季和秋季,在中间层顶附近两种行星波活动均较弱;250km处经向行星波活动略强于纬向,经向风不同周期带的行星波最强振幅主要出现在5-9月,与电离层f₀F₂振荡特性的研究结果一致.      

Thermospheric tides during thermosphere mapping study periods

Neutral exospheric temperatures at 53°, 43° and 33° latitude from Millstone Hill steerable-antenna Thomson scatter measurements, and at 19° latitude from the Arecibo Observatory, obtained during three Thermosphere Mapping Study (TMS) coordinated campaign intervals during 1984 and 1985, are analyzed for diurnal and semidiurnal tidal components. The resulting amplitude and phase latitudinal structures are compared with numerical simulations. The observed semidiurnal tidal components are thought to be significantly affected by tidal waves propagating upwards from below the thermosphere during these solar minimum periods. We speculate that current inadequacies in specifying F-region plasma densities and mean zonal winds at lower altitudes within the simulation model may account for certain discrepancies between observations and theory.     

Comparison of observed hmF2 and IRI 2007 model with M(3000)F2 estimation of hmF2 at low solar activity for an equatorial station

This is to investigate ways of improving the Equatorial F2-layer peak heights estimated from M(3000)F2 ionosonde data measured using the Ionospheric Prediction Service (IPS-42) sounder at Ouagadougou, Burkina Faso (Latitude +12.4°N, Longitude +1.5°W, Dip latitude +5.9°N) during a low solar activity year (1995). For this purpose, we have compared the observed h_mF2 (h_mF2_obs) deduced using an algorithm from scaled virtual heights of quiet day ionograms and the predicted h_mF2 values which is given by the IRI 2007 model (h_mF2_{IRI 2007}) with the ionosonde measured M(3000)F2 estimation of the h_mF2 values (h_mF2_est) respectively. The correlation coefficients R² for all the seasons were found to range from 0.259 to 0.692 for h_mF2_obs values, while it ranges from 0.551 to 0.875 for the h_mF2_{IRI 2007} values. During the nighttime, estimated h_mF2 (h_mF2_est) was found to be positively correlated with the h_mF2_obs values by the post-sunset peak representation which is also represented by the h_mF2_{IRI 2007} values. We also investigated the validity of the h_mF2_est values by finding the percentage deviations when compared with the h_mF2_obs and h_mF2_{IRI 2007}.     

Improvements in bottomside electron density definition in the Autoscala program

Some improvements introduced in the Autoscala program are presented. They include improvements in E valley modeling of the electron density profile N_e(h), and in the link between the E valley and bottom-side F regions. An abrupt variation in N_e(h) generated by the previous version of Autoscala under night conditions has been eliminated.A series of ionograms recorded by the Millstone Hill digisonde (42.6°, 288.5°) were automatically interpreted by the previous version of Autoscala and by the new one. Data from Incoherent Scatter Radar (ISR) were used to comparatively assess the performance of the two versions. For this purpose, the root mean square errors (RMSEs) of the N_e(h) provided by Autoscala were calculated relative to the corresponding values provided by ISR.A more accurate overall modeling of N_e(h) was achieved by the new Autoscala version (RMSE = 0.51 MHz for the new version against RMSE = 0.67 MHz for the previous one).     

Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.