Atmosphere sounding by GPS radio occultation: First results from TanDEM-X and comparison with TerraSAR-X

On 21 June 2010 the TerraSAR-X satellite was joined by the TanDEM-X satellite. A Global Positioning System (GPS) radio occultation (RO) experiment using the twin satellites has been carried out to estimate the precision of GPS atmospheric soundings. For the Day Of Year (DOY) 330–336, 2011, we analyze phase and amplitude data recorded by GPS receivers separated by a few hundred meters in a low earth orbit and derive collocated atmospheric refractivity profiles. In the altitude range 10–20 km the standard deviation between TerraSAR-X and TanDEM-X refractivity does not exceed 0.15%. The standard deviation is rapidly increasing for lower and higher altitudes; close to the surface and at an altitude of 30 km the standard deviation reaches 0.8% and 0.5%, respectively. Systematic deviations between TerraSAR-X and TanDEM-X refractivity in the considered altitude range (0–30 km) are negligible. The results confirm the anticipated high precision of the GPS RO technique. However, the difference in the retrieved refractivity in the lower troposphere for different Open Loop (OL) signal tracking parameters, altered onboard TanDEM-X for DOY 49–55, 2012, calls for an in depth analysis. At the moment we can not exclude that a potential bias in the OL Doppler model introduces a bias in our retrieved refractivity at altitudes <8$<8$ km.     

Ionospheric scintillations at Guilin detected by GPS ground-based and radio occultation observations

The occurrence of ionospheric scintillations with S4 ? 0.2 was studied using GPS measurements at Guilin, China (25.29°N, 110.33°E; geomagnetic: 15.04°N, 181.98°E), a station located near the northern crest of the equatorial anomaly. The results are presented for data collected from January 2009 to March 2010. The results show that nighttime amplitude scintillations only took place in February and March of the considered years, while daytime amplitude scintillations occurred in August and December of 2009. Nighttime amplitude scintillations, observed in the south of Guilin, always occurred with phase scintillations, TEC (Total Electron Content) depletions, and ROT (Rate Of change of TEC) fluctuations. However, TEC depletions and ROT fluctuations were weak during daytime amplitude scintillations, and daytime amplitude scintillations always took place simultaneously for most of the GPS satellites which appeared over Guilin in different azimuth directions. Ground-based GPS scintillation/TEC observations recorded at Guilin and signal-to-noise-ratio (SNR) measurements obtained from GPS-COSMIC radio occultation indicate that nighttime and daytime scintillations are very likely caused by ionospheric F region irregularities and sporadic E, respectively. Moreover, strong daytime amplitude scintillations may be associated with the plasma density enhancements in ionospheric E region caused by the Perseid and Geminid meteor shower activities.     

Change in refractivity of the atmosphere and large variation in TEC associated with some earthquakes,observed from GPS receiver

The present study reports the analysis of GPS based TEC for our station Surat (21.16°N, 72.78°E) located at the northern crest of equatorial anomaly region in India at times close to some earthquake events (M ? 5) during the year 2009 in India and its neighbouring regions. The TEC data used in the study are obtained from GPS Ionospheric Scintillation and TEC Monitoring (GISTM) system. The TEC data has been analysed corresponding to 11 earthquakes in low solar activity period and quiet geomagnetic condition. We found that, out of 11 cases of earthquakes (M > 5) there were seven cases in which enhancement in TEC occurred on earthquake day and in other four cases there was depletion in TEC on earthquake day. The variation in refractivity prior to earthquake was significant for the cases in which the epicentre lied within a distance of 600 km from the receiving station. By looking into the features on temporal enhancement and depletion of TEC a prediction was made 3–2 days prior to an earthquake (on 28 October 2009 in Bhuj – India). The paper includes a brief discussion on the method of potentially identifying an impending earthquake from ionospheric data.     

Ionospheric response of equatorial and low latitude F-region during the intense geomagnetic storm on 24–25 August 2005

In this investigation, we present and discuss the response of the ionospheric F-region in the South American and East Asian sectors during an intense geomagnetic storm in August 2005. The geomagnetic storm studied reached a minimum Dst of −216 nT at 12:00 UT on 24 August. In this work ionospheric sounding data obtained of 24, 25, and 26 August 2005 at Palmas (PAL; 10.2° S, 48.2° W; dip latitude 6.6° S), São José dos Campos (SJC, 23.2° S, 45.9° W; dip latitude 17.6° S), Brazil, Ho Chi Minh City, (HCM; 10.5° N, 106.3° E; dip latitude 2.9° N), Vietnam, Okinawa (OKI; 26.3° N, 127.8° E; dip latitude 21.2° N), Japan, are presented. Also, the GPS observations obtained at different stations in the equatorial and low-latitude regions in the Brazilian sector are presented. On the night of 24–25 August 2005, the h′F variations show traveling ionospheric disturbances associated with Joule heating in the auroral zone from SJC to PAL. The foF2 variations show a positive storm phase on the night of 24–25 August at PAL and SJC during the recovery phase. Also, the GPS-VTEC observations at several stations in the Brazilian sector show a fairly similar positive storm phase on 24 August. During the fast decrease of Dst (between 10:00 and 11:00 UT) on 24 August, there is a prompt penetration of electric field of magnetospheric origin that result in abrupt increase (∼12:00 UT) in foF2 at PAL, SJC (Brazil) and OKI (Japan) and in VTEC at IMPZ, BOMJ, PARA and SMAR (Brazil). OKI showed strong oscillations of the F-region on the night 24 August resulted to the propagation of traveling atmospheric disturbances (TADs) by Joule heating in the auroral region. These effects result a strong positive observed at OKI station. During the daytime on 25 August, in the recovery phase, the foF2 observations showed positive ionospheric storm at HCM station. Some differences in the latitudinal response of the F-region is also observed in the South American and East Asian sectors.     

基于STK的GPS/LEO无线电掩星技术仿真研究

GPS/LEO无线电掩星技术反演地球大气参数剖面已经具有较高的精度. 国外开展了多个GPS/LEO掩星项目, 但中国还尚未深入进行相关的实验, 这制约了中国掩星技术的发展. 本文提出基于STK进行GPS/LEO掩星技术研究的方法; 根据GPS/LEO掩星的原理, 推导出掩星事件发生的条件和掩星切点的计算公式; 利用STK对掩星过程进行模拟, 得到掩星数据. 在大气球对称假设和大气模型已知的条件下, 反演得到中性大气折射指数. 通过比较模型和反演数据, 表明反演数据精度较高, 验证了利用STK模拟GPS/LEO掩星实验方法的可行性.      

Difficulties in the study of cosmic radio noise absorption at 30 MHz using riometer at low latitude station,Kolhapur (Lat-16.8°N,Long-74.25°E)

A dual dipole antenna has been installed at low latitude station Kolhapur (Geographic 16.8°N, 74.25°E), Maharashtra, India for the study of cosmic radio noise absorption using Solid State Riometer (which operates at 30 MHz) during pre phase of 24th solar maxima. The aim for this type of study over Kolhapur was to know the response of lower (D region) ionosphere over low latitude by cosmic radio noise absorption using riometer technique during quite period as well as sudden ionospheric disturbances (SID). The observations are being taken for 3 years. Two different sites (∼40 km away from each other) were used for the installation of riometer equipment assuming minimum local noise. It is found that solar noise to cosmic radio noise hence resulting in signal saturation. The night time signal is relatively free of interference but sometimes local noise is responsible for spike-like signatures. Hence it is concluded that Kolhapur (a low latitude station) is not suitable for the study of cosmic radio noise absorption on 30 MHz with riometer and dual dipole antenna. Proper choice for operating frequency of riometer and antenna gain is suggested for low latitude use of this technique for ionospheric deviative and nondeviative absorption studies.     

The impact of microwave absorber and radome geometries on GNSS measurements of station coordinates and atmospheric water vapour

We have used microwave absorbing material in different geometries around ground-based Global Navigation Satellite System (GNSS) antennas in order to mitigate multipath effects on the estimates of station coordinates and atmospheric water vapour. The influence of a hemispheric radome – of the same type as in the Swedish GPS network SWEPOS – was also investigated. Two GNSS stations at the Onsala Space Observatory were used forming a 12 m baseline. GPS data from October 2008 to November 2009 were analyzed by the GIPSY/OASIS II software using the Precise Point Positioning (PPP) processing strategy for five different elevation cutoff angles from 5° to 25°. We found that the use of the absorbing material decreases the offset in the estimated vertical component of the baseline from ∼27 mm to ∼4 mm when the elevation cutoff angle varies from 5° to 20°. The horizontal components are much less affected. The corresponding offset in the estimates of the atmospheric Integrated Water Vapour (IWV) decreases from ∼1.6 kg/m² to ∼0.3 kg/m². Changes less than 5 mm in the offsets in the vertical component of the baseline are seen for all five elevation cutoff angle solutions when the antenna was covered by a hemispheric radome. Using the radome affects the IWV estimates less than 0.4 kg/m² for all different solutions. IWV comparisons between a Water Vapour Radiometer (WVR) and the GPS data give consistent results.     

Coordinated in situ measurements of plasma irregularities and ground based scintillation observations at the crest of equatorial anomaly

First comparison of in situ density fluctuations measured by the DEMETER satellite with ground based GPS receiver measurements at the equatorial anomaly station Bhopal (geographic coordinates (23.2°N, 77.6°E); geomagnetic coordinates (14.29°N, 151.12°E)) for the low solar activity year 2005, are presented in this paper. Calculation of the diurnal maximum of the strength of the equatorial electrojet, which can serve as precursor to ionospheric scintillations in the anomaly region is also done. The Langmuir Probe experiment and Plasma Analyzer onboard DEMETER measure the electron and ion densities respectively. Irregularities in electron density distribution cause scintillations on transionospheric links and there exists a close relationship between an irregularity and scintillation. In 40% of the cases, DEMETER detects the irregularity structures (dNe/Ne ? 5% and dNi/Ni (O⁺) ? 5%) and GPS L band scintillations (S4 ? 0.2) are also observed around the same time, for the low solar activity period. It is found that maximum irregularity intensity is obtained in the geomagnetic latitude range of 10–20° for both electron density and ion density. As the GPS signals pass through this irregularity structure, scintillations are recorded by the GPS receiver installed at the equatorial anomaly station, Bhopal it is interesting to note that in situ density fluctuations observed on magnetic flux tubes that pass over Bhopal can be used as indicator of ionospheric scintillations at that site. Many cases of density fluctuations and associated scintillations have been observed during the descending low solar activity period. The percentage occurrence of density irregularities and scintillations shows good correspondence with diurnal maximum of the strength of electrojet, however this varies with different seasons with maximum correspondence in summer (up to 66%) followed by equinox (up to 50%) and winter (up to 46%). Also, there is a threshold value of EEJ strength to produce density irregularities ((dNe/Ne)_max ? 5%) and for moderate to strong scintillations (S4 ? 0.3) to occur. For winter this value is found to be ∼40 nT whereas for equinox and summer it is around 50 nT.     

ROPP反演软件算法及其精度分析

介绍了ROPP反演软件中无线电掩星反演的算法与精度分析. 采用COSMIC卫星2008年1 月1日全天的附加相位数据, 反演得到折射率、温度、压强与湿度等参数, 并与CDAAC 相应结果进行对比. 实验结果表明, 在30km高度以下, 折射率、压强和湿度的相对 误差在2%以内, 温度误差不超过2K.      

Anomalous enhancement of electric field derived from ionosonde data before the great Wenchuan earthquake

The seismo-ionospheric coupling is a hot topic of modern scientific research. One of the most reasonable mechanisms is the hypothesis of quasistatic electric field. With a preliminary analytical method, electric field in F2-region before Wenchuan earthquake was derived from the five lower latitude ground-based ionosondes: Lhasa (29.63°N, 91.17°E), Chongqing (29.50°N, 106.40°E), Kunming (25.00°N, 102.70°E), Guangzhou (23.15°N, 113.35°E) and Haikou (23.15°N, 110.35°E). The results show that there was an anomalous enhancement of electric field in F2-region close to the epicenter with maximal amplitude ∼2 mV/m (about 10 times of the background) at 15:00 LT on 9 May, which is in line with the amplitudes of anomalous electric field related to some previous earthquakes observed by satellites.     

Occurrence rate of SAR arcs during the 23rd solar activity cycle

The occurrence rate of SAR arcs during 1997–2007 has been analyzed based on the photometric observations at the Yakutsk meridian (Maimaga station, corrected geomagnetic coordinates: 57°N, 200°E). SAR arcs appeared in 114 cases (∼500 h) during ∼370 nights of observations (∼3170 h). The occurrence frequency of SAR arcs increases to 27% during the growth phase of solar activity and has a clearly defined maximum at a decline of cycle 23. The SAR arc registration probability corresponds to the variations in geomagnetic activity in this solar cycle. The dates, intervals of UT, and geomagnetic latitudes of SAR arc observations at the Yakutsk meridian are presented.     

Temperature comparison between CHAMP radio occultation and TIMED/SABER measurements in the lower stratosphere

Global Positioning System (GPS) receiver on the CHAllenging Mini-satellite Payload (CHAMP) and the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument, one of four on board the TIMED satellite, provide middle atmosphere temperature profiles by Radio Occultation (RO) and limb viewing infrared emission measurements, respectively. These temperature profiles retrieved by two different techniques in the stratosphere are compared with each other using more than 1300 correlative profiles in March, September and December 2005. The over-all mean differences averaged over 15 and 35 km are approximately −2 K and standard deviation is less than 3 K. Below 20 km of altitude, relatively small mean temperature differences ∼1 K are observed in wide latitudinal range except for June (during the SABER nighttime observation). In the middle to low latitudes, between 30°S and 30°N, the temperature difference increases with height from ∼0–1 K at 15 km, to ∼−4 K at 35 km of altitude. Large temperature differences about −4 to −6 K are observed between 60°S and 30°N and 31–35 km of altitude for all months and between 0° and 30°N below 16 km during June (nighttime).     

An alternative Shell inversion technique - Analysis and validation based on COSMIC and ionosonde data

Multi-channel Global Positioning System (GPS) carrier phase signals, received by the six low Earth orbiting (LEO) satellites from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) program, were used to undertake active limb sounding of the Earth’s atmosphere and ionosphere via radio occultation. In the ionospheric radio occultation (IRO) data processing, the standard Shell inversion technique (SIT), transformed from the traditional Abel inversion technique (AIT), is widely used, and can retrieve good electron density profiles. In this paper, an alternative SIT method is proposed. The comparison between different inversion techniques will be discussed, taking advantage of the availability of COSMIC datasets. Moreover, the occultation results obtained from the SIT and alternative SIT at 500 km and 800 km, are compared with ionosonde measurements. The electron densities from the alternative SIT show excellent consistency to those from the SIT, with strong correlations over 0.996 and 0.999 at altitudes of 500 km and 800 km, respectively, and the peak electron densities (N_mF₂) from the alternative SIT are equivalent to the SIT, with 0.839 vs. 0.844, and 0.907 vs. 0.909 correlation coefficients when comparing to those by the ionosondes. These results show that: (1) the N_mF₂ and h_mF₂ retrieved from the SIT and alternative SIT are highly consistent, and in a good agreement with those measured by ionosondes, (2) no matter which inversion technique is used, the occultation results at the higher orbits (∼800 km) are better than those at the lower orbits (∼500 km).     

Assessment of precision in ionospheric electron density profiles retrieved by GPS radio occultations

The Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) is a six satellite radio occultation mission that was launched in April 2006. The close proximity of these satellites during some months after launch provides a unique opportunity to evaluate the precision of Global Positioning System (GPS) radio occultation (RO) retrievals of ionospheric electron density from nearly collocated and simultaneous observations. RO data from 30 consecutive days during July and August 2006 are divided into ten groups in terms of daytime or nighttime and latitude. In all cases, the best precision values (about 1%) are found at the F peak height and they slightly degrade upwards. For all daytime groups, it is seen that electron density profiles above about 120 km height exhibit a substantial improvement in precision. Nighttime groups are rather diverse: in particular, the precision becomes better than 10% above different levels between 120 and 200 km height. Our overall results show that up to 100–200 km (depending on each group), the uncertainty associated with the precision is in the order of the measured electron density values. Even worse, the retrieved values tend sometimes to be negative. Although we cannot rely directly on electron density values at these altitudes, the shape of the profiles could be indicative of some ionospheric features (e.g. waves and sporadic E layers). Above 200 km, the profiles of precision are qualitatively quite independent from daytime or latitude. From all the nearly collocated pairs studied, only 49 exhibited a difference between line of sight angles of both RO at the F peak height larger than 10°. After analyzing them we find no clear indications of a significant representativeness error in electron density profiles due to the spherical assumption above 120 km height. Differences in precision between setting and rising GPS RO may be attributed to the modification of the processing algorithms applied to rising cases during the initial period of the COSMIC mission.     

Towards a GPS-based TEC prediction model for Southern Africa with feed forward networks

In this paper, first results from a national Global Positioning System (GPS) based total electron content (TEC) prediction model over South Africa are presented. Data for 10 GPS receiver stations distributed through out the country were used to train a feed forward neural network (NN) over an interval of at most five years. In the NN training, validating and testing processes, five factors which are well known to influence TEC variability namely diurnal variation, seasonal variation, magnetic activity, solar activity and the geographic position of the GPS receivers were included in the NN model. The database consisted of 1-min data and therefore the NN model developed can be used to forecast TEC values 1 min in advance. Results from the NN national model (NM) were compared with hourly TEC values generated by the earlier developed NN single station models (SSMs) at Sutherland (32.38°S, 20.81°E) and Springbok (29.67°S, 17.88°E), to predict TEC variations over the Cape Town (33.95°S, 18.47°E) and Upington (28.41°S, 21.26°E) stations, respectively, during equinoxes and solstices. This revealed that, on average, the NM led to an improvement in TEC prediction accuracy compared to the SSMs for the considered testing periods.     

Oblique-incidence ionospheric soundings over Central Europe and their application for testing now casting and long term prediction models

After a first oblique-incidence ionospheric sounding campaign over Central Europe performed during the period 2003–2004 over the radio links between Inskip (UK, 53.5°N, 2.5°W) and Rome (Italy, 41.8°N, 12.5°E) and between Inskip and Chania (Crete, 35.7°N, 24.0°E), new and more extensive analysis of systematic MUF measurements from January 2005 to December 2006 have been performed. MUF measurements collected during moderately disturbed days (17 ? Ap ? 32), disturbed days (32 < Ap ? 50) and very disturbed days (Ap > 50), have been used to test the long term prediction models (ASAPS, ICEPAC and SIRM&LKW), and the now casting models (SIRMUP&LKW and ISWIRM&LKW). The performances of the different prediction methods in terms of r.m.s are shown for selected range of geomagnetic activity and for each season.     

Response of the EIA ionosphere to the 7-8 May 2005 geomagnetic storm

In this paper, response of low latitude ionosphere to a moderate geomagnetic storm of 7–8 May 2005 (SSC: 1920 UT on 7 May with Sym-H minimum, ∼−112 nT around 1600 UT on 8 May) has been investigated using the GPS measurements from a near EIA crest region, Rajkot (Geog. 22.29°N, 70.74°E, Geomag.14°), India. We found a decrease in total electron content (TEC) in 12 h after the onset of the storm, an increase during and after 6 h of Sym-H deep minimum with a decrease below its usual-day level on the second day during the recovery phase of the storm. On 8 May, an increase of TEC is observed after sunset and during post-midnight hours (maximum up to 170%) with the formation of ionospheric plasma bubbles followed by a nearly simultaneous onset of scintillations at L-band frequencies following the time of rapid decrease in Sym-H index (−30 nT/h around 1300 UT).     

Measurements of the upper troposphere and lower stratosphere during tropical cyclones using the GPS radio occultation technique

Water vapour transport to the upper troposphere and lower stratosphere by deep convective storms affects the radiation balance of the atmosphere and has been proposed as an important component of climate change. The aim of the work presented here is to understand if the GPS radio occultation technique is useful for characterization of this process. Our assessment addresses the question if severe storms leave a significant signature in radio occultation profiles in the upper troposphere/lower stratosphere. Radio occultation data from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) were analyzed, focusing on two particular tropical cyclones with completely different characteristics, the hurricane Bertha, which formed in the Atlantic Basin during July 2008 and reached a maximum intensity of Category 3, and the typhoon Hondo, which formed in the south Indian Ocean during 2008 reaching a maximum intensity of Category 4. The result is positive, suggesting that the bending angle of a GPS radio occultation signal contains interesting information on the atmosphere around the tropopause, but not any information regarding the water vapour. The maximum percentage anomaly of bending angle between 14 and 18 km of altitude during tropical cyclones is typically larger than the annual mean by 5–15% and it can reach 20% for extreme cases. The results are discussed in connection to the GPS radio occultation receiver which will be part of the Atomic Clock Ensemble in Space (ACES) payload on the International Space Station.     

Signatures of the Sudden Stratospheric Warming events of January–February 2008 in Seoul,S. Korea

The period January–February 2008 was characterized by four Sudden Stratospheric Warmings (SSWs) in the Northern Hemisphere, of which the last warming, at the end of February 2008, was a major warming. A significant decrease in mesospheric water vapour (H₂O) of more than 2 ppmv (∼40%) was observed by the ground-based microwave (GBMW) radiometer in Seoul, S. Korea [37.3°N, 126.3°E] during the major SSW. A comparison with ground-based mesospheric H₂O observations from the mid-latitude station in Bern [46.9°N, 7°E] revealed an anticorrelation in the mesospheric H₂O data during the major SSW. In addition, prior to the major warming, strong periodic fluctuations were recorded in the Aura MLS vertical temperature distribution between 15 and 0.05 hPa at Seoul. The mesospheric temperature oscillation was found to have a period of ∼10–14 days with a persistency of 3–4 cycles.     

Night-time total electron content enhancements at equatorial anomaly region in China

This paper presents small scale (duration ?1 h, ΔTEC ? 1TECU) night-time total electron content (TEC) enhancements observed at the equatorial anomaly region in China, for the first time. The data is from a GPS receiver chain established in 2005 by Institute of Center for Space Science and Applied Research, Chinese Academy of Sciences and a GPS receiver of International GPS Service (IGS), located between Fuzhou (26.1°N, 119.3°E) and Nanning (22.8°N, 108.3°E). Two other GPS observations of IGS taken at higher latitude were also used to investigate the localization of such phenomenon. The characteristics of the night-time TEC enhancement are examined with two case studies. The TEC increases about 1–3TECU, intermittently. While the night-time TEC enhancement mainly occurs at the equatorial anomaly region, it can be observed at middle latitude as well. The spatial size of the enhancement region is less than 5° in longitude. The primary statistical study shows that the TEC enhancement is more often observed in spring and autumn, but rarely in summer. It has no dependence on geomagnetic activity. The enhancement can occur both before and after midnight.