Sensing snow height and surface temperature variations in Greenland from GPS reflected signals

The in situ measurements of snow surface temperature (SST) and snow height (SH) are very difficult with high costs, particularly in Greenland Ice Sheet (GrIS). Since the snow depth variations coupling with surface temperature are related to GPS multipath, it is possible to estimate the snow depth and surface air temperature variations by incorporating GPS-Reflectometry (GPS-R). In this paper, the reflected signals from ground GPS receivers are used to sense the SST and SH variations based on the thermophysical behavior and variations of snow layer from April to June 2010 at SMM1 site and from March to December 2010 at MARG site in Greenland. The results show that the mean daily changes in the ionospheric geometrical-free linear combination (GPS-L4) of dual-frequency GPS signals are related to daily SST and SH variations. The nonparametric bootstrapping model in direct (forward) and inverse models are developed and applied to estimate the SST and SH variations. The mean biases of SST and SH estimates are 0.18 °C and 0.23 m at SMM1 site, respectively, and 3.8 °C and 0.13 m at MARG site, respectively.     

Direct and indirect estimates of turbulence around the turbopause

Direct and indirect estimates of turbulence parameters at altitudes 80–120 km are reviewed. It is shown that there are contradictions in data on absolute value, shape, seasonal variations etc. of eddy diffusion coefficient K_t revealed or accepted by different authors. The strongest contradiction is in the views on the seasonal variations of turbulence. Possible explanation may be connected with the correct account for mean motions. Data on turbopause height obtained from neutral composition data (Ar, N₂) are discussed. The revealed reverse connection of this height with the temperature at 120 km is discussed in terms of K_t dependence on temperature gradient.     

Height and time variability of planetary wave activity

The height–season and year-to-year regularities of parameters of first and second spatial harmonics determine the structure of the stratosphere and mesosphere circulation and its variability. In the period 1992–2002 at heights 0–55 km, the amplitudes and phases of the first and second spatial harmonics in the field of temperature, geopotential height, zonal and meridional wind were calculated by the method of harmonic decomposition. Dispersion (standard or mean square deviation) of their day-to-day and year-to-year variations was calculated by their wavelength constants. Height and season patterns of variability have been estimated. The difference in height–longitude variability for wave numbers m = 1 and 2 has been discovered. At the same time, the intensity of wave disturbances for m = 1 is less than for m = 2 excluding the polar areas, where a significant variability appears at the heights 0–55 km. There is also a tendency for the intensity of year-to-year variations to decrease in comparison with day-to-day variations. In cold and warm periods the amplitude of perturbation waves with m = 2 both for day-to-day and year-to-year variations is greater than for waves with m = 1. Transient height areas in the interval of 20–30 km are more distinct for day-to-day variations of polar area.     

Comparison of topside ionosphere scale height modeled by the Topside Sounder Model and incoherent scatter radar ionospheric model

The topside ionosphere scale height extracted from two empirical models are compared in the paper. The Topside Sounder Model (TSM) provides directly the scale height (H_T), while the incoherent scatter radar ionospheric model (ISRIM) provides electron density profiles and its scale height (H_R) is determined by the lowest gradient in the topside part of the profile. H_T and H_R are presented for 7 ISR locations along with their dependences on season, local time, solar flux F10.7, and geomagnetic index ap. Comparison reveals that H_T values are systematically lower than respective H_R values as the average offset for all 7 stations is 55 km. For the midlatitude stations Arecibo, Shigaraki, and Millstone Hill this difference is reduced to 43 km. The range of variations of H_R is much larger than that of H_T, as the H_T range overlaps the lower part of the H_R range. Dependences on ap, DoY and LT are much stronger in the ISRIM than in TSM. This results in much larger values of H_R at higher ap. Diurnal amplitude of H_R is much larger than that of H_T, with large maximum of H_R at night. The present comparison yields the conclusion that the ISR measurements provide steeper topside Ne profiles than that provided by the topside sounders.     

Comparison of mesopause region meteor radar winds,medium frequency radar winds and low frequency drifts over Germany

During 2004 and 2005 measurements of mesospheric/lower thermospheric (80–100 km) winds have been carried out in Germany using three different ground-based systems, namely a meteor radar (36.2 MHz) at the Collm Observatory (51.3°N, 13°E), a MF radar (3.18 MHz) at Juliusruh (54.6°N, 13.4°E) and the LF D1 measurements using a transmitter (177 kHz) at Zehlendorf near Berlin and receivers at Collm with the reflection point at 52.1°N, 13.2°E. This provides the possibility of comparing the results of different radar systems in nearly the same measuring volume. Meteor radar winds are generally stronger than the winds observed by MF and especially by LF radars. This difference is small near 80 km but increases with height. The difference between meteor radar and medium frequency radar winds is larger during winter than during summer, which might indicate an indirect influence of gravity waves on spaced antenna measurements.     

The propagation features of LF radio waves at topside ionosphere and their variations possibly related to Wenchuan earthquake in 2008

The artificial low frequency (LF) signals recorded at topside ionosphere from DEMETER satellite were analyzed in this paper, and the typical diurnal and seasonal variations were illustrated around 162?kHz in electric field spectra. The larger power spectrum density (PSD) values in electric field at local nighttime and in winter season all demonstrate the correlation feature of lower plasma content to higher penetration of LF waves into ionosphere. Around Wenchuan earthquake, the comparison of signal-noise-ratio (SNR) values in electric field with each half month during January to May in 2008 and the same half month in May from 2005 to 2007 revealed their lowest values and small covering area around the preparation region of Wenchuan earthquake in 2008. Combined with other researches in VLF radio waves and geochemical observations from satellite, the interaction of ion accumulation and upward movement from gas-water release at surface might be a key factor to disturb the ionospheric plasma density, and then possibly leading to the decrease of low energy penetration of LF radio waves from the artificial transmitted source at ground.     

Variations in the ionospheric scale height parameter at the F2 peak over Grahamstown,South Africa

The Grahamstown, South Africa (33.3°S, 26.5°E) ionospheric field station operates a Lowell digital pulse ionospheric sounder (Digisonde) whose output includes scaled parameters derived from the measured ionogram. One of these output parameters is the ionospheric scale height parameter (H), and this paper presents an analysis of the seasonal, diurnal, and solar activity variations of this parameter over the Grahamstown station. Ionosonde data from three years 2002, 2003, and 2004 were used in this study. The data was subjected to a general trend analysis to remove any outliers and then the monthly median data were used to explore the different variations. The results of this analysis were found to be similar to what has already been presented in the literature for low latitude stations, and are presented as well as the correlation at this mid-latitude station between the H parameter, the IRI shape parameter (B0), and the peak electron density (NmF2).     

Reconstruction of topside density profile by using the topside sounder model profiler and digisonde data

To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (H_s), O⁺–H⁺ transition height (H_T), and their ratio Rt = H_s/H_T, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.     

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.     

A generalization of the Hines” dispersion relation

A simple dispersion relation which forms a generalization to the Hines' relation has been developed for the real atmosphere which fully allows for reflection and transmission at every height level. The relation is an analytic function of the Brunt frequency, the speed of sound profile and its derivatives and the variations in mean molecular mass (using CIRA 1972). The correction terms to the Hines' dispersion relation become important for gravity waves with frequencies close to the Brunt frequency.     

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.     

Response of low latitude D-region ionosphere to the total solar eclipse of 22 July 2009 deduced from ELF/VLF analysis

Response of the D-region of the ionosphere to the total solar eclipse of 22 July 2009 at low latitude, Varanasi (Geog. lat., 25.27° N; Geog. long., 82.98° E; Geomag. lat. = 14° 55’ N) was investigated using ELF/VLF radio signal. Tweeks, a naturally occurring VLF signal and radio signals from various VLF navigational transmitters are first time used simultaneously to study the effect of total solar eclipse (TSE). Tweeks occurrence is a nighttime phenomena but the obscuration of solar disc during TSE in early morning leads to tweek occurrence. The changes in D-region ionospheric VLF reflection heights (h) and electron density (n_e: 22.6–24.6 cm⁻³) during eclipse have been estimated from tweek analysis. The reflection height increased from ∼89 km from the first occurrence of tweek to about ∼93 km at the totality and then decreased to ∼88 km at the end of the eclipse, suggesting significant increase in tweek reflection height of about 5.5 km during the eclipse. The reflection heights at the time of totality during TSE are found to be less by 2–3 km as compared to the usual nighttime tweek reflection heights. This is due to partial nighttime condition created by TSE. A significant increase of 3 dB in the strength of the amplitude of VLF signal of 22.2 kHz transmitted from JJI-Japan is observed around the time of the total solar eclipse (TSE) as compared to a normal day. The modeled electron density height profile of the lower ionosphere depicts linear variation in the electron density with respect to solar radiation as observed by tweek analysis also. These low latitude ionospheric perturbations on the eclipse day are discussed and compared with other normal days.     

Plasma transport process in the equatorial/low-latitude ionosphere

The behaviour of the equatorial/low-latitude ionosphere and the transport processes during magnetic disturbed and quiet periods of a high solar activity year, 2014, in the American sector are investigated. Parameters used include vertical drift (Vz), transport term (W), NmF2, hmF2 and scale-height (H). The F2 plasma variations followed the diurnal local solar pattern, being higher at daytime. The sunset maximum and sunrise minimum peaks of hmF2 were directly opposite to the scale height (H) pattern. The plasma distribution was basically controlled by combined actions of the electrodynamic convection/thermospheric composition, which is geomagnetic activity dependent. The annual, semi-annual and winter-anomalies of the F2 parameters were higher at the dip equator in comparison with the low-latitude. The Vz pre-reversal peak magnitude coincided with hmF2 peak and the effects are more pronounced during geomagnetic disturbed conditions. The transport term pattern was similar to that of the scale height and it is suggested as a proxy parameter for quantifying low-latitude plasma irregularities and distribution of thermospheric composition.     

武昌低电离层LF观测的季节变化

通过分析武昌低电离层LF观测数据,得到了中低纬度低电离层电子密度剖面的季节变化特征及其与太阳天顶角周年变化的关系,给出LF相位周年振荡幅度的年际变化曲线．     

Comparing topside and bottomside-measured characteristics of the F2 layer peak

The ionospheric characteristics of the F2 layer peak have been measured with ionosondes from the ground or with satellites from space. The most common characteristics are the F2-peak density NmF2 and peak height hmF2. In addition to these two parameters this paper studies the F2-peak scale height. Comparing the median values of hmF2 and NmF2 obtained from topside and bottomside sounding shows good agreement in general. The Chapman scale height values for the F2 layer peak derived from topside profiles, H_m,top, are generally several times larger than H_m,bot derived from bottomside profiles.     

Mean solar quiet daily variations in the earth’s magnetic field along East African longitudes

Solar quiet daily (Sq) variation in the earth’s magnetic field along the East African meridian was studied using data of the H, D and Z components recorded with Magnetic Data Acquisition System of SERC. One year data recorded at ten African geomagnetic observatories was used in the analysis of worldwide solar quiet daily variation (Wsq). The study revealed that the focus of Sq (H) in the southern hemisphere lies at the boundary of low and middle latitude region. Noon-time enhancement of Sq (H) was generally noticed at all stations along the meridian, though it is latitudinal dependent in terms of magnitude as it reduces with distance from dip equator. In addition, night-time variations also occur in small magnitude along African meridian in Sq (H) and Sq (Z) which could be attributed to non-ionospheric sources. Semi-diurnal variation was noticed in Sq (D) at all stations except in AAB that is under the influence of electrojet current. Dusk sector calm condition of Sq (D) current was notice in some stations and the same condition was also noticed at dawn sector in some other stations. The usual sunrise maximum and sunset minimum for D component at stations north of dip equator as well as sunrise minimum and sunset maximum was found to increase with distance away from dip equator. Day-time perturbation of Sq current was noticed to be more pronounced in all the three field elements. Mass plots of annual mean hourly value show contrasting phase pattern about the focus in H element and the results of the variations at each region with the associated standard error. It was concluded from the result of correlation coefficients computed that different currents system flowing in opposite directions could be responsible for contrasting patterns.     

Temperature variability in the tropical mesosphere during the northern hemisphere winter

Temperature observations at 20–90 km height and 5°N–15°N during the period of December 1992–March 1993 from the WINDII and MLS experiments on the UARS satellite are analysed together with MF radar winds and UKMO assimilated fields of temperature and zonal and meridional winds. The correlation between the different datasets at the tropics and zonal mean wind data at mid latitudes is examined for period February–March 1993, when series of stratospheric warming events were observed at middle and high latitudes. Wavelet analysis is applied to investigate coupling between stationary and travelling planetary waves in the stratosphere and the upper mesosphere. Planetary waves m = 1 with periods of 4–7 days, 8–12 days and 13–18 days are found to dominate the period. Westward 7- and 16–18 day waves at the tropics appear enhanced by stationary planetary waves during sudden stratospheric warming events.     

北京地区大气温度及重力波活动的季节变化

利用瑞利激光雷达观测数据,分析了北京地区35~70km高度范围内大气温度和重力波活动的季节变化.发现北京地区30~70km高度范围内的大气温度有明显的年周期变化:平流层顶最高温度出现在6,7月份,大约为270K;中间层70km高度最低温度也出现在6,7月份,大约为200K.以2014年10月14日晚数据为例,分析重力波势能密度,发现50km以下重力波势能存在耗散,而在50km以上重力波近乎无耗散地向上传播.通过对比35~50km高度范围内的平均势能密度,对北京地区重力波活动强弱的季节变化进行了研究.研究结果表明,北京上空重力波活动强度具有明显的年周期变化,冬季平均势能密度为18J·kg-1,夏季为8J·kg-1,且冬季重力波活动强度约为夏季的两倍.此外,还分析了春夏秋冬四个季节重力波势能密度随高度的变化.结果表明,不同季节和不同高度的重力波势能密度不同.      

Detection of long-term trends in the mesosphere/lower thermosphere from ground-based radio propagation measurements

Mesospheric temperature trends can be derived from LF phase-height observations in mid-latitudes supported by ionospheric absorption and ionosonde observations. Analysing the full observation period from 1959 until 2003, a mean yearly temperature trend has been derived with −0.25 K/yr for the height interval from 48 to 82 km. Subdividing the whole observation interval in two parts before and after 1979, the trend is markedly stronger in the second period with −0.38 K/yr compared with −0.20 K/yr in the first part before 1979. These differences can at least partly be explained by a steeper CO₂ increase and ozone decrease in the second interval. The differences in the mesospheric temperature trends are most evidently expressed during winter months and are markedly smaller during summer season. The reason of this seasonal difference is not quite clear; it may be related with detected ozone trends which are clearly stronger during winter months on both hemispheres.