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.     

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

利用中国岢岚站（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₂振荡特性的研究结果一致.      

中国廊坊中间层和低热层大气平均风观测模拟

利用中国廊坊站（39.4°N,116.7°E）流星雷达在2012年4月1日至2013年3月31日的水平风场观测数据,分析廊坊上空80~100km的中间层与低热层（Mesosphere and Lower Thermosphere,MLT）大气平均纬向风和经向风的季节变化特征.结果表明平均纬向风和经向风都表现出明显的季节变化特征.平均纬向风在冬季MLT盛行西风,极大值位于中间层顶,随高度增加西风减弱;在夏季中间层为东风,低热层为强西风,风向转换高度约为82km.平均经向风在冬季以南风为主,在夏季盛行北风.纬向风和经向风在春秋两季主要表现为过渡阶段.流星雷达观测结果与WACCM4模式和HWM93模式模拟的气候变化特点基本一致,但WACCM4模式纬向风和经向风风速偏大,而HWM93模式纬向风和经向风风速偏小.      

Mid-latitude thermospheric wind changes during the St. Patrick’s Day storm of 2015 observed by two Fabry-Perot interferometers in China

In this work, we utilize thermospheric wind observations by the Fabry-Perot interferometers (FPI) from the Kelan (KL) station (38.7°N, 111.6°E, Magnetic Latitude: 28.9°N) and the Xinglong (XL) station (40.2°N, 117.4°E, Magnetic Latitude: 30.5°N) in central China during the St. Patrick’s Day storm (from Mar. 17 to Mar. 19) of 2015 to analyze thermospheric wind disturbances and compare observations with the Horizontal Wind Model 2007 (HWM07). The results reveal that the wind measurements at KL show very similar trends to those at XL. Large enhancements are seen in both the westward and equatorward winds after the severe geomagnetic storm occurred. The westward wind speed increased to a peak value of 75?m/s and the equatorward wind enhanced to a peak value of over 100?m/s. There also exist obvious poleward disturbances in the meridional winds during Mar. 17 to Mar. 19. According to the comparison with HWM07, there exist evident wind speed and temporal differences between FPI-winds and the model outputs in this severe geomagnetic storm. The discrepancies between the observations and HWM07 imply that the empirical model should be used carefully in wind disturbance forecast during large geomagnetic storms and more investigations between measurements and numerical models are necessary in future studies.     

Effects of CME and CIR induced geomagnetic storms on low-latitude ionization over Indian longitudes in terms of neutral dynamics

This paper presents the response of the ionosphere during the intense geomagnetic storms of October 12–20, 2016 and May 26–31, 2017 which occurred during the declining phase of the solar cycle 24. Total Electron Content (TEC) from GPS measured at Indore, Calcutta and Siliguri having geomagnetic dips varying from 32.23°N, 32°N and 39.49°N respectively and at the International GNSS Service (IGS) stations at Lucknow (beyond anomaly crest), Hyderabad (between geomagnetic equator and northern crest of EIA) and Bangalore (near magnetic equator) in the Indian longitude zone have been used for the storms. Prominent peaks in diurnal maximum in excess of 20–45 TECU over the quiet time values were observed during the October 2016 storm at Lucknow, Indore, Hyderabad, Bangalore and 10–20 TECU for the May 2017 storm at Siliguri, Indore, Calcutta and Hyderabad. The GUVI images onboard TIMED spacecraft that measures the thermospheric O/N₂ ratio, showed high values (O/N₂ ratio of about 0.7) on October 16 when positive storm effects were observed compared to the other days during the storm period. The observed features have been explained in terms of the O/N₂ ratio increase in the equatorial thermosphere, CIR-induced High Speed Solar Wind (HSSW) event for the October 2016 storm. The TEC enhancement has also been explained in terms of the Auroral Electrojet (AE), neutral wind values obtained from the Horizontal Wind Model (HWM14) and equatorial electrojet strength from magnetometer data for both October 2016 and May 2017 storms. These results are one of the first to be reported from the Indian longitude sector on influence of CME- and CIR-driven geomagnetic storms on TEC during the declining phase of solar cycle 24.     

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.     

Mid-latitude thermospheric zonal winds during the equinoxes

The diurnal variation of the mid-latitude upper thermosphere zonal winds during equinoxes has been studied using data recently generated from CHAMP measurements from 2002 to 2004 using an iterative algorithm. The wind data was separated into two geomagnetic activity levels, representing high geomagnetic activity level (Ap > 8) and low geomagnetic activity level (Ap ? 8). The data were further separated into two solar flux levels; with F10.7 > 140 for high and F10.7 ? 140 for low. Geomagnetic activity is a correlator just as significant as solar activity. The response of mid-latitude thermospheric zonal winds to increases in geomagnetic disturbances and solar flux is evident. With increase in geomagnetic activity, midday to midnight winds are generally less eastward and generally more westward after the about midnight transitions. The results show that east west transitions generally occurred about midnight hours for all the situations analyzed. The west to east transition occurs from 1400–1500 MLT. Enhanced westward averaged zonal wind speeds going above 150 ms⁻¹ are observed in the north hemisphere mid-latitude about sunrise hours (∼0700–1100 MLT). Nighttime winds in the north hemisphere are in good agreement with previous single station ground observations over Millstone Hill. Improved ground observations and multi satellite observations from space will greatly improve temporal coverage of the Earth’s thermosphere.     

A global study of meridional winds in the thermosphere

We present neutral meridional winds derived from existing ground-based measurements of the height of the F2-layer maximum electron density (h_max). The method of calculation uses an ionospheric model to determine the relationship between h_max and the neutral wind along a magnetic meridian. The meridional wind is derived from a comparison of the modeled layer height with the measured height. This method is used in the global mapping of thermospheric winds using radar and ionosonde measurements from the Global Thermospheric Mapping Study at the summer and winter solstices. The diurnal behavior of the meridional wind is calculated for the two seasons from Millstone Hill radar measurements and from several ionosonde stations. The major features of the winds calculated from Millstone Hill radar data are a 200 m/s southward wind between 0200 and 0400 hours local time at the summer solstice and a northward daytime wind in winter. A sampling of winds derived from ionosonde data shows the diurnal pattern to vary with geographic latitude and longitude. Nighttime equatorward winds are found to be larger in Europe than at other locations of similar latitude.     

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.     

Seasonal,annual, and interannual variability in MLT quasi-two-day waves over the low-latitude region Kolhapur (16.8°N; 74.2°E)

This study presents the quasi-two-day wave (Q2DW) characteristics of the mesosphere and lower thermosphere (MLT) region obtained by taking hourly mean values of horizontal wind velocities for 4? years (August 2013–July 2017) through continuous measurements using a medium-frequency (MF) radar (operating frequency – 1.98?MHz) located at the low-latitude Indian station Kolhapur (16.8°N; 74.2°E). The MF radar located at Kolhapur was upgraded in 2013, and these results of Q2DW have been reported for the first time after upgrading. The present study investigated variability in seasonal, annual, interannual, and solar indices of Q2DWs traveling in zonal (EW) and meridional (NS) components in the MLT region. The Q2DW activity is observed to be stronger during austral summer (January–February) (EW?=?～5?m/s and NS?=?～8–10?m/s) than during boreal summer (June–July) (EW = ～5 m/s and NS = ～6–8?m/s). The Q2DW amplitudes are larger in the meridional component than in the zonal one. A strong semiannual oscillation (SAO) has been observed in Q2DWs, with peak during January–February and June–July. In addition, small enhancement is seen in meridional Q2DW in October (～5–6?m/s). It is observed that the entire spectrum (40–60?h) measured between 86 and 94?km contributes to the SAO amplitudes during January–February and June–July, whereas the waves measured between 42?h and 52?h contribute to enhancement in October similar to that reported elsewhere. In general, the Q2DW amplitude shows large interannual variability. The easterlies developed in the global circulation model in Northern hemisphere during May intensify up to around summer solstice. Q2DW activity peaks during westerly shear zone and intensifies with time at a lower thermospheric altitude (above 90?km). Small positive correlations (r?=?0.2 for sunspot number and r?=?0.1 for 10.7?cm solar flux) have been observed between Q2DW amplitudes and solar activity.     

Development of wind observational models for the upper atmosphere of the earth based on the joint analysis of direct and indirect sounding data

Earlier latitudinal distribution models of zonal winds were developed mainly along the 80°W meridian [1, 2]. In this paper an attempt is made to take into account longitudinal differences in zonal and meridional wind distributions. These are considerable in winter prediods.     

Early winter meteor wind over Bologna (45°N, 11°E)

The Bologna meteor radar was operational during two winter campaigns, from 6 January 1982 to 1 February 1982 and from 10 December 1982 to 2 February 1983. As occurrence of minor stratospheric warmings has been reported for these intervals, possible effects on meteor wind over Bologna related to this kind of warming are pointed out. Zonal and meridional prevailing winds are found to exhibit the maximum peak to peak value in their oscillations when a minor stratospheric warming reaches such an intensity that ΔT_{(90°N–60°N)} at 10 mbar is reversed. Diurnal and semidiurnal tides exhibit usual amplitude variations, but the semidiurnal tide has a noticeable phase shift at the time of a minor warming occurrence, while a similar shift is less evident in the diurnal tide phase.     

The 16-day waves in the mesosphere and lower thermosphere over Wuhan (30.6°N, 114.5°E) and Adelaide (35°S, 138°E)

Winds from a meteor radar at Wuhan (30.6°N, 114.5°E) and a MF radar at Adelaide (35°S, 138°E) are used to study the 16-day waves in the mesosphere and lower thermosphere (MLT). The height range is 78–98 km at Wuhan and 70–98 km at Adelaide. By comparison, it is found that the zonal components at both sites are generally larger than the meridional ones, and eastward motion of the zonal background winds is favorable for the 16-day waves penetration to the MLT region. The zonal maximum amplitude appears in the autumn (September–October) around 86–98 km at Wuhan and in the winter months and early spring (July–October) around 72–82 km at Adelaide. Differences are found in wave amplitudes and time of appearance between the two years of 2002 and 2003. In 2003, the intensity of the wave amplitudes is relatively smaller than that for 2002 at both sites. The summer 16-day waves are comparatively weaker at Adelaide in both years, but stronger in 2002 at Wuhan near the mesopause and the lower thermosphere (86–98 km). The strong summer waves at Wuhan may come from the winter southern hemisphere.     

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.     

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.     

电离层水平电场与风场的耦合模拟研究

设计了一个将电离层水平电场与风场耦合的模拟方案,研究了电流函数和风场在耦合前后的变化与差异. 研究发现,水平电场与风场相互反馈后,风场的变化比电流函数小. 经向风在白天有较明显的差异,夜晚的差异比白天小,主要出现在中高纬地区,并随高度的增加而增大,300km左右达到最大值,其后几乎保持不变. 纬向风有与经向风相似的变化,但纬向风耦合前后的差异比经向风小. 电流函数在耦合后有较大改变,两个涡旋强度都有较强增加,并且北半球的增强大于南半球,而夜晚差异较小. 结果表明,在研究的高度范围内,风场对电场的控制作用大于电场对风场的影响.      

Some features of the day-to-day MLT wind variability in winter 2017–2018 as seen with a European/Siberian meteor radar network

We present results of wind measurements near the mesopause carried out with meteor radars (MRs) at Collm (51°N, 13°E), Obninsk (55°N, 37°E), Kazan (56°N, 49°E), Angarsk (52°N, 104°E) and Anadyr (65°N, 178°E) from October 1, 2017 till March 31, 2018. The Collm and Kazan MRs are SKiYMET radars with vertical transmission and radio echo height finding, while the other radars operate with horizontal transmission and without height finding. We paid particular attention to the meridional wind variability with periods of 4–6 days and 9–11 days. The waves with these periods are seen as spots of the wave activity in the wavelet spectra and include oscillations with different periods and different discrete zonal wavenumbers. These wave packets successively propagate as a group of waves from one site to another one in such a way that they are observed at one site and almost disappear at the previous one. The 4–6 wave group includes planetary-scale oscillations (individual spectral components) which have eastward phase velocities and mostly zonal wavenumbers 2 and 3, and the vertical wavelength exceeds 70 km at middle latitudes. The source of the oscillations is the polar jet instability. The wave group itself propagates westward, and the amplitudes of wind oscillations are approximately 5–6 m/s as obtained from the wind data averaged over the meteor zone. The 9–11 day wave set propagates westward as a group and mainly consists of spectral components which have westward phase velocity and zonal wavenumber 1. Amplitudes of these wind perturbations strongly vary from station to station and can reach, approximately, 8 m/s. The vertical wavenumber is 0.014 km⁻¹ as taken from the Kazan and 0.05 km⁻¹ according to the Collm data. We obtained a global view on the waves by using the AURA MLS geopotential data. We found a good correspondence between wave features obtained from the MR wind measurements and the MLS data. To our knowledge, such a wave propagation of planetary wave in the mesosphere/lower thermosphere (MLT) region has so far not obtained much attention.     

Climatology of the quasi two-day wave over Saskatoon (52°N, 107°W): 14 Years of MF radar observations

Oscillations in the wind field due to the “quasi two-day” wave (Q2DW) have been studied using 14 years of Medium Frequency (MF) radar data over Saskatoon (52°N). The results of analyses have shown that the Q2DW exhibits complex seasonal variability: in addition to the well-known strong summer activity the wave also has been found in winter months. In summer the amplitudes of the Q2DW in meridional (NS) and zonal (EW) winds have comparable values, while they tend to be larger in EW winds in winter. The heights of the maximum amplitude are lower in winter than in summer and there are indications of changes in the dominant period of oscillations throughout the year.     

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.