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.     

Simultaneous lidar observation of peculiar sporadic K and Na layers at São José dos Campos (23.1°S, 45.9°W), Brazil

We present the first simultaneous observation of mesopause sodium (Na) and potassium (K) layer by a lidar which has Na and K channel simultaneously at the South Hemisphere site, São José dos Campos (23.1°S, 45.9°W). Measurements reported here were conducted on two nights with 3.5 and 8?h of observations in November 2016. On 20 November 2016, sporadic layers in both Na and K layer occurred above 100?km, and the higher layers corresponded well with sporadic E (Es) layer. And the density of Na at 100–105?km is higher than that at main layer around 90?km, but K density around 100?km is at least 3 times smaller compared with the K density around 90?km for the similar period. On 21 November 2016, both sporadic layers occurred in main layer height with obvious descending variations with time, which seems like tidal induced. Notably, the peak K/Na ratio slowly increased with time. And Na layer and K layer showed different processes along with time with K density reaching its maximum 1?h later than that of Na. Correlations of Na/K density, Es, and winds were also discussed.     

武汉上空中层和低热层大气潮汐的流星雷达观测

武汉流星雷达是2002年元月建成的我国第一部全天空流星雷达，本文对2002年2月19日到7月31日流星雷达观测的潮汐的讨论表明，武汉中层顶以周日潮汐为潮汐运动的主要分量，它的强度远大于半日潮汐，周日潮汐和半日潮汐的波源都在80km以下．周日潮汐分量在3、4月份最强，并且经向分量略强于纬向分量．两个分量的峰值在约95km处出现，分别达到44m／s和60m／s．半日潮的最大值24m／s出现在4月初约93km处．周日潮汐和半日潮汐的振幅和相位随时间呈现出拟周期变化的特征，这可能是潮汐与行星波非线形相互作用的结果．观测结果与GSWM模型的比较表明，GSWM模型在相位随高度变化趋势上与观测结果一致，但模型的周日潮相位比观测约超前1—2h，半日潮相位约滞后1—4h．在周日潮汐较强的月份，模型与观测有较大的差异，观测的幅度通常在95km附近有极大值，而模型并没有极大值．GSWM模型对半日潮的幅度的估计通常过小，观测的半日潮汐幅度有时甚至超过模型值的一倍以上．     

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.     

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.     

The mesopause region wind field over Central Europe in 2003 and comparison with a long-term climatology

Using the D1 method in the LF range, monthly mean zonal and meridional winds in the mesosphere/lower thermosphere height range are measured continuously at Collm, Germany. The wind data are accompanied by reference height measurements since 1983, allowing the construction of an empirical long-term climatology of mean winds and tidal amplitudes and phases at 52N, 15E, which covers the height range of 80–110 km in winter, and about 83–107 km in summer. The climatology includes the time interval from 1983 to date, thus covering nearly two solar cycles. Vertical wind parameter profiles can also be constructed for individual years, so that from the time series interannual variability and long-term trends can be deduced. Here, we present the results for 2003 in comparison with the climatological means.     

A review on the upper atmospheric sodium observations from India: Insights

This review article briefly brings out the historical development of atmospheric sodium (Na) measurements over India and the importance of coordinated measurements with multiple techniques to address physical processes in the Earth's upper atmosphere. These measurements were initiated in the early 1970s by observing Na airglow emission intensities with broad band airglow photometer from Mt. Abu, a low-latitude hill station in India. Considerable amount of night-to-night variations in nocturnal emission intensities of the Na airglow were observed. Later, investigations regarding the dependence with the magnetic activities from the equatorial and low latitudes were carried out and double-humped structures in the nocturnal variation of intensities were reported. With the advent of Na lidar at Gadanki around 2005, the measurements of atmospheric neutral Na atom concentration became possible and more frequent occurrences of sporadic Na layers over the magnetic low latitude station compared to other latitudes were detected indicating the role played by electrodynamics. Later, a possible relationship between E-region field aligned plasma irregularities and the concentration of neutral Na atoms was investigated using coordinated measurements of VHF radar and Na lidar. Further, simultaneous measurements with Na lidar and a narrow band airglow photometer with narrow field of view brought out the importance of coordinated observation wherein the characterization of gravity waves could be carried out and also revealed the importance of collisional quenching due to ambient molecules in the Na airglow emission processes. In addition, combining the ground based measurements of Na lidar and meteor wind radar along with satellite measurements made possible to hypothesize the over-turning Kelvin–Helmholtz billow in the Na layer manifests “C-type” or inverted lambda shape structures in the height-time-concentration map of neutral Na atoms. This review paper presents a synoptic view mostly based on the previously reported observations of Na airglow emission, Na lidar and coordinated Na airglow and Na lidar observations from the Indian sector and highlights the importance of simultaneous measurements of mesospheric Na and its emissions along with satellite-borne measurements to address interesting geophysical processes in the Earth's upper atmosphere.     

大气中层顶区域波相互作用的一个观测个例

利用SOUSY VHF雷达的观测数据分析了极区中层顶83.4-91.2km范围内大气风场波动的非线性相互作用。大气风场的谱在不同高度上均有明显的潮汐分量峰值，纬向风分量中35h波、半日潮和8.9h惯必重力波构成共振相互作用对，经向风分量中33h波、半日潮和19h惯性重力波构成共振相互作用对。双谱分析表明，这些共振对在许多高度上都发生耦合，35h或33h波振幅的极小值与半日潮的极大值出现的高度几乎相同，呈现出明显的非线性相互作用在空间上不是局域的，而是存在于中层顶区域的几乎所有高度上，这种相互作用不仅导致半日潮振幅随时间的变化，也使半日潮的振幅随空间变化。35h和33h波动可能是在其他时段或其他位置通过行星波与周日潮相互作用产生的，然后传播到观测点并与半日潮发生相互作用。     

中纬度冬季低热层潮汐水平风分量相位关系的MF雷达观测

采用武汉(30°N,114°E)MF雷达在2001年冬季的风场观测数据研究中纬度低热层大气潮汐水平风分量之间的相位关系.统一用弧度定义的各潮汐经纬向分量的拟合初相位在三个连续的高度上分别显示出相同的时间变化倾向和相近的相位差,但是在绝大多数观测时间△ψ24和△ψ12准正交,而△ψ8出乎意料地准同相.周日、半日和8 h潮汐经纬向分量的二次相位耦合(QPC)方程被分别估计出来,利用它们相减还得到一个潮汐相位差相关方程.推测的8 h潮汐相位和相位差与相应的观测值很好地符合.在第14个时间窗内,三个潮汐一般表现为椭圆偏振而不是圆偏振或线偏振,但是△ψ24和△ψ12在三个连续的高度上准正交,而△ψ8在92.0和94.0 km上准同相.因此估计的潮汐QPC方程、推导的潮汐相位差相关方程、观测的8 h潮汐准同相相位差以及典型的潮汐偏振图都是观测的周日、半日和8 h潮汐之间真实QPC的反映.      

北京上空冬季潮汐波活动特征全天时激光雷达探测研究

利用北京延庆子午工程激光雷达对北京上空钠层进行长期连续观测,分析研究钠层及其相关参数的周日变化.提取钠层各个高度上的相位信息,与同时段经向风潮汐信息进行比较发现:在各个高度上,二者周日相位数值基本一致.半日相位对比结果表明,虽然存在差异,但整体仍保持较好的一致性.此外,从2014年至2016年每年10月到第二年1月共4个月的钠层连续观测数据中提取钠层周日和半日的振幅和相位,探究北京上空冬季的潮汐特征.结果显示:周日潮的相位自上向下传播,且无明显的季节变化特征,其垂直波长在40～50km的范围;周日潮较强,半日潮较弱.      

Simultaneous observation of dual-site airglow imagers and a sodium temperature-wind lidar,and effect of atmospheric stability on the airglow structure

The heights of horizontal structures in the OH airglow images were observed with two all-sky imagers at Platteville (40.2N, 104.7W) and Yucca Ridge Field Station (40.7N, 104.9W) in Colorado. The sodium temperature-wind lidar of CSU at Fort Collins (40.6N, 105.1W) observed temperature and wind velocities around 80–105-km altitude, providing time-height variations of atmospheric stability. Simultaneous observations with the imagers and the lidar were successful on two nights: November 18 and 19, 2003. The estimated height from OH airglow images were lower (80–86 km) than those observed in Japan, with significant time variation. Comparison with the temperature and wind observed by the lidar showed that the estimated height is likely to coincide with the height of small N² (close to 0) values. Ripples (<10-km horizontal scale) were simultaneously observed in the image.     

Comparisons between primrose lake (54°N, 110°W) ROCOB winds (20–60 km) and Saskatoon (52°N, 107°W) M. F. radar winds (60–110 km): 1978–1982

The development of the new CIRA will require the combination of winds from many sources, e.g. rockets (ROCOB) up to ～60 km, and radar winds ～60–110 km. Difficulties are that such rocket data have larger errors at 60–65 km, and tidal effects may become significant. Radar data for 60–80 km may also have tidal contamination, due to ? 16h of data per day: from 80–110 km tidal corrections are usually reliable.Comparisons are made between the unique Saskatoon MF radar set, which is continuous from mid 1978–1983, and the ROCOB data from Primrose Lake, which is only 340 km northwest. While the agreement is satisfactory, special care is required when matching the two regions: particular problems are the low rocket sampling rate, and the unexpectedly large amplitude of the diurnal tide. Important differences from the zonal winds of CIRA-72 emerge, especially in winter months. Meridional cross-sections differ from previous data models in the extent of the summer equatorward flow.     

Comparisons of winds at (44°S, 173°E), 65–110km,with CIRA 72

In the 95km height region of the atmosphere, ground-based techniques made an important contribution to the CIRA 72 [1] wind model. Recent wind measurements from a partial reflection experiment at 44S covering one and a half years are presented and compared with CIRA 72. The zonal wind component compares favourably, although the measured values are more easterly above 80km in autumn and winter; a feature of the autumn winds is a temporary easterly reversal above 90km. Winter mesospheric winds can be very disturbed. The summer mesosphere easterly maximum appears earlier in the season and at a higher altitude than the model. A much poorer comparison is shown between the measured meridional wind component and the 1969 model of Groves [2].     

用电离层特性参量提取等效风场信息

导出了利用中低纬电离层特性参量获取电离层F层峰区高度上等效风场（包含电场和风场信息在内）的基本方程，并尝试用该方法从电离层特性参量（峰高和临频）提取等效风场信息，利用武汉站DGS-256电离层数字测高仪数据及由美国Massachusetts Lowell大学最新版的剖面反演程序换算得到F层峰高，获得了武汉地区夏季至日点附近，冬季至日点附近，冬季地磁特别宁静的九天和冬季平均等效风场的初步特征，并利用Fejer经验电场模式计算冬季电场引起的垂直漂移，估计电场和风场对武汉地区的垂直等效风场的贡献大小，结果表明：等效风场呈现出白天与夜晚幅度和方向的差异。至日点附近冬季与夏季白天的幅度差异以及明显的凌晨凹陷现象；平均情况下，垂直等效风场幅度和方向的变化主要是由中性风引起，受电场的影响不大。     

Large-scale circulation of atomic oxygen in the upper mesosphere and lower thermosphere

The variability of the atomic oxygen green line airglow at 557.7 nm, originating from the O(¹S) level, has a long history of observation. Only recently, global observations of the oxygen airglow, interpreted with the help of global circulation models have provided some understanding. Satellite observations of winds and temperatures clearly demonstrate the dynamical influence of tides on the daily variations. Both annual and semi-annual components of the variation have been identified as occurring over the course of the year. From the large-scale circulation of the atmosphere one expects downwelling in the winter at high latitudes, causing enhanced atomic oxygen and enhanced airglow while the opposite is expected in summer. In the present investigation a search is made for this large-scale signature using data from the WIND Imaging Interferometer on the Upper Atmosphere Research Satellite, model results from the Thermosphere–Ionosphere–Mesosphere-Electrodynamics-Global Circulation Model, and earlier published results. A well-defined semi-annual variation of emission rate is found in the tropics, apparently the result of the semi-annual variation of the diurnal tide. Annual and semi-annual patterns are found at mid-latitudes, in satellite and model data, while ground-based observations detect only an annual variation, with a maximum in the autumn. At still higher latitudes the fall peak persists, but with a deep depletion of atomic oxygen in the springtime; this appears to be the signature of the large-scale circulation.     

Global structure,seasonal and interannual variability of the eastward propagating tides seen in the SABER/TIMED temperatures (2002–2007)

The present paper is focused on the global spatial (altitude and latitude) structure, seasonal and interannual variability of the most stable in amplitude and phase eastward propagating diurnal and semidiurnal tides with zonal wavenumbers 2 and 3 derived from the SABER/TIMED temperatures for full 6 years (January 2002–December 2007). The tidal results are obtained by an analysis method where the tides (migrating and nonmigrating) and the planetary waves (zonally travelling, zonally symmetric and stationary) are simultaneously extracted from the satellite data. It has been found that the structures of the eastward propagating diurnal tides with zonal wavenumbers 3 and 2 change from antisymmetric with respect to the equator below ∼85 km height, to more symmetric above ∼95 km. The seasonal behavior of the DE3 is dominated by annual variation with maximum in August–September reaching average (2002–2007) amplitude of ∼15 K, while that of the DE2 by semiannual variation with solstice maxima and with average amplitude of ∼8 K. These tides revealed some interannual variability with a period of quasi-2 years. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 2 in the southern hemisphere (SH) is dominated by annual variation with maximum in the austral summer (November–January) while that in the northern hemisphere (NH) by semiannual variation with equinoctial maxima. The SE2 maximizes near 115 km height and at latitude of ∼30° reaching an average amplitude of ∼6 K. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 3 in both hemispheres indicates a main maximum during June solstice and a secondary one during December solstice. The tide maximizes near 110–115 km height and at a latitude of ∼30° reaching an average amplitude of ∼4.8 K in the SH and ∼4 K in the NH. The tidal structures of the two eastward propagating semidiurnal tides are predominantly antisymmetric about the equator.     

Temporal variation of the volume emission rates between OI5577 and O2(0,1)

The intensities of the upper mesospheric airglow emissions, OI 557.7 nm and O₂(b) atmospheric (0,1) band at 864.5 nm, have been measured since October 1998 using a ground-based multichannel airglow photometer located near the equator at São João do Cariri (7S, 35W). The intensity ratio between the two emissions, OI5577/O₂b(0,1), was obtained as a function of time. The annual average of the ratio R in 1998 was 0.41 ± 0.14, and showed considerable day to day variations, ranging from 0.2 to 0.6. Nocturnal variations of the ratio also demonstrate a large amplitude of oscillation, varying between 0.3 to 0.6. These results were compared with a model calculation.     

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.     

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.     

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.