Different mechanisms of the ionospheric influence on GPS occultation signals

A local mechanism for strong ionospheric effects on radio occultation (RO) global positioning satellite system (GPS) signals is described. Peculiar zones centered at the critical points (the tangent points) in the ionosphere, where the gradient of the electron density is perpendicular to the RO ray trajectory, strongly influence the amplitude and phase of RO signals. It follows from the analytical model of local ionospheric effects that the positions of the critical points depend on the RO geometry and the structure of the ionospheric disturbances. Centers of strong ionospheric influence on RO signals can exist, for example, in the sporadic E-layers, which are inclined by 3–6° relative to the local horizontal direction. Also, intense F2 layer irregularities can contribute to the RO signal variations. A classification of the ionospheric influence on the GPS RO signals is introduced using the amplitude data, which indicates different mechanisms (local, diffraction, etc.) for radio waves propagation. The existence of regular mechanisms (e.g., local mechanism) indicates a potential for separating the regular and random parts in the ionospheric influence on the RO signals.     

GPS radio holography as a tool for remote sensing of the atmosphere, mesosphere, and terrestrial surface from space

GPS radio occultation (RO) signals are highly coherent and precise, and thus sufficient for holographic investigation of the atmosphere, ionosphere, and the Earth's surface from space. In principle, three-dimensional radio-holographic remote sensing is possible by using new radio holographic equations to retrieve the radio field within the atmosphere from a radio field known at some interface outside the atmosphere. A simplified two-dimensional form of the radio-holographic equations which are developed under an assumption of local spherical symmetry can be used to obtain two-dimensional radio images of the atmosphere and terrestrial surface. To achieve this, radio holograms recorded by a GPS receiver onboard a low earth orbit (LEO) satellite at two GPS frequencies can be used and focused synthetic aperture principle applied. Analysis of GPS/MET RO data is presented to show the effectiveness of a radio-holographic approach. It is shown that the amplitude of GPS radio signals (in addition to phase data) can be used to obtain detailed altitude profiles of the vertical gradient of refractivity in the atmosphere and electron density in the mesosphere. The results demonstrate the applicability of GPS radio holography for a detailed global study of the natural processes in the atmosphere and mesosphere. Electronic Publication     

Analysis of atmospheric and ionospheric structures using the GPS/MET and CHAMP radio occultation database: a methodological review

Since 1995, the global positioning system (GPS) has been exploited by the means of the radio occultation (RO) method to obtain the vertical profiles of refractivity, temperature, pressure, and water vapor in the neutral atmosphere and electron density in the ionosphere. Applying the RO method to the study of the Earths atmosphere was demonstrated for the first time with the GPS/MET experiment. Since then, several satellites with GPS receivers, suitable for RO experiments, have been launched including Oersted, SUNSAT, CHAMP, SAC-C, and GRACE. Future RO investigations that are planned now include FORMOSAT3/COSMIC and Terra-SAR missions. New elements in the RO technology are required to meet the goals of improving the accuracy and broadening the potential of the RO method. In this paper, a methodological review of RO investigations is presented to emphasize new directions in applying the RO method: measuring the vertical gradients of the refractivity in the atmosphere and electron density in the lower ionosphere, determination of the temperature regime in the upper stratosphere, investigation of the internal wave activity in the atmosphere, and study of the ionospheric disturbances on a global scale. These new directions may be relevant for investigating the relationships between processes in the atmosphere and mesosphere, the study of thermal regimes in the intermediate heights of the upper stratosphere—lower mesosphere, and the analysis of the influence of the space weather phenomena on the lower ionosphere.     

The M 2 ocean tide loading wave in Alaska: vertical and horizontal displacements, modelled and observed

Crustal deformations caused by surface load due to ocean tides are strongly dependent on the surface load closest to the observation site. In order to correctly model this ocean loading effect near irregular coastal areas, a high-resolution coastline is required. A test is carried out using two GPS sites located in Alaska, where the ocean tide loading effect is large and consequently observed easily by relative positioning with GPS. The selected sites are Fair (Fairbanks) and Chi3 (located on an island that separates Prince William Sound from the Gulf of Alaska). Processing of hourly baseline solutions between Fair and Chi3 over a period of 49 days yields a significant ocean tide loading effect. The data are processed using different strategies for the tropospheric delay correction. However, the best results are obtained when 1-h ZTD (Zenith Tropospheric Delay) parameters for hourly solutions are used. In this case ocean tide loading is not absorbed into the ZTD parameters. Hence, ocean tide loading can be well resolved in the GPS data analysis. In addition, the M ₂ ocean tide wave in the Gulf of Alaska has a very large amplitude. Although the horizontal M ₂ ocean tide loading amplitude in general is only about 1/4 of the vertical M ₂ ocean tide loading amplitude, the differential horizontal M ₂ ocean tide loading displacements are nevertheless measurable using differential GPS (DGPS). When using the GOT99.2 ocean tide model and taking the coastal structure into account, the predicted differential vertical M ₂ amplitude and Greenwich phase lag due to ocean tide loading are 19.3 mm and 110.2 degrees respectively, while GPS measurements yield 21.3 ± 1.0 mm and 99.7±2.8 degrees. Similarly, the predicted differential horizontal M ₂ amplitude and Greenwich phase lag (in the north–south direction) are 4.5 mm and –77.0 degrees, while GPS yields 5.4 ± 0.3 mm and –106.3±3.3 degrees. Only the north-south component of the differential horizontal M ₂ ocean tide loading wave is considered, because the east–west component is too small for the processed baseline and not detectable using DGPS.     

Quality assessment of COSMIC/FORMOSAT-3 GPS radio occultation data derived from single- and double-difference atmospheric excess phase processing

This study evaluates the quality of GPS radio occultation (RO) atmospheric excess phase data derived with single- and double-difference processing algorithms. A spectral analysis of 1 s GPS clock estimates indicates that a sampling interval of 1 s is necessary to adequately remove the GPS clock error with single-difference processing. One week (May 2–8, 2009) of COSMIC/FORMOSAT-3 data are analyzed in a post-processed mode with four different processing strategies: (1) double-differencing with 1 s GPS ground data, (2) single-differencing with 30 s GPS clock estimates (standard COSMIC Data Analysis and Archival Center product), (3) single-differencing with 5 s GPS clocks, and (4) single-differencing with 1 s GPS clocks. Analyses of a common set of 5,596 RO profiles show that the neutral atmospheric bending angles and refractivities derived from single-difference processing with 1 s GPS clocks are the highest quality. The random noise of neutral atmospheric bending angles between 60 and 80 km heights is about 1.50e−6 rad for the single-difference cases and 1.74e−6 rad for double-differencing. An analysis of pairs of collocated soundings also shows that bending angles derived from single-differencing with 1 s GPS clocks are more consistent than with the other processing strategies. Additionally, the standard deviation of the differences between RO and high-resolution European Center for Medium range Weather Forecasting (ECMWF) refractivity profiles at 30 km height is 0.60% for single-differencing with 1 and 5 s GPS clocks, 0.68% for single-differencing with 30 s clocks, and 0.66% for double-differencing. A GPS clock-sampling interval of 1 s or less is required for single- and zero-difference processing to achieve the highest quality excess atmospheric phase data for RO applications.     

FORMOSAT-3/COSMIC science mission update

The FORMOSAT-3/COSMIC mission is a microsatellite mission for weather forecast, climate monitoring, and atmospheric, ionospheric and geodesy research. This mission is a collaborative Taiwan-USA science experiment to deploy a constellation of six microsatellites in low Earth orbits. The mission life is 2 years with a goal of 5 years. The final mission orbit has an altitude of 750–800 km. Each satellite consists of three science payloads: global positioning system (GPS) occultation experiment (GOX) payload, tiny ionospheric photometer (TIP) and tri-band beacon (TBB). The GOX will collect the GPS signals for the study on atmosphere, ionosphere, and geodesy. The TIP and TBB can provide the electron distribution information for ionospheric research. The deployment of the FORMOSAT-3 constellation and the resulting influence on the occultation sounding distributions are reported. Details are also given on GOX, TIP, and TBB payload operations and the contributions of the Taiwan Science Team.     

Ionosphere probing with simultaneous GNSS radio occultations

GPS Solutions - Radio occultation (RO) is a powerful technique for providing vertical profiles of refractivity, temperature, pressure, and water vapor of the neutral atmosphere and electron density...     

Analysis on the global morphology of stratospheric gravity wave activity deduced from the COSMIC GPS occultation profiles

Monthly mean global morphologies of potential energy density E _p from stratospheric gravity waves are revealed by observations of COSMIC GPS radio occultation. The E _p is obtained from vertical wavelengths ranging from 2 to 10 km over cells of 1° × 2° in latitude and longitude. The computed values confirm previous results and obtain new ones. The large gravity wave E _p values found in the tropics between 25°N and 25°S could be mainly due to the strong tropical cumulus convection; July values are larger than those for January (2007). In mid and high latitudes, the most prominent features of the northern winter hemisphere are the enhanced densities above the Eurasian continent and the North Atlantic and the depressed E _p values above the North Pacific and North America for which topography, wind sources and wind filtering may be responsible. In southern winter hemisphere, large E _p values are found around 180° and 300° longitudes that are likely due to the topography of the Antarctic plateau, the Antarctic Peninsula and South America. Enhanced E _p values are found over Scandinavia. However, there is no clear evidence to show that gravity waves are localized over the Rocky Mountains, the Himalayas and the Andes. Topography and planetary wave modulations are proposed to interpret the large-scale longitudinal variations and inter-hemisphere asymmetry of the GW activity.     

Phase acceleration: a new important parameter in GPS occultation technology

Based on 40 years of radio-occultation (RO) experiments, it is now recognized that the phase acceleration of radio waves (equal to the time derivative of the Doppler shift), derived from analysis of high-stability Global Positioning System (GPS) RO signals, is as important as the Doppler frequency. The phase acceleration technique allows one to convert the phase and Doppler frequency changes into refractive attenuation variations. From such derived refractive attenuation and amplitude data, one can estimate the integral absorption of radio waves. This is important for future RO missions when measuring water vapor and minor atmospheric gas constituents, because the difficulty of removing the refractive attenuation effect from the amplitude data can be avoided. The phase acceleration technique can be applied also for determining the location and inclination of sharp layered plasma structures (including sporadic E _s layers) in the ionosphere. The advantages of the phase acceleration technique are validated by analyzing RO data from the Challenging Minisatellite Payload (CHAMP) and the FORMOSA Satellite Constellation Observing Systems for Meteorology, Ionosphere, and Climate missions (FORMOSAT-3/COSMIC).     

Use of the L2C signal for inversions of GPS radio occultation data in the neutral atmosphere

Results from processing FORMOSAT-3/COSMIC radio occultations (RO) with the new GPS L2C signal acquired both in phase locked loop (PLL) and open loop (OL) modes are presented. Analysis of L2P, L2C, and L1CA signals acquired in PLL mode shows that in the presence of strong ionospheric scintillation not only L2P tracking, but also L1CA tracking often fails, while L2C tracking is most stable. The use of L2C improves current RO processing in the neutral atmosphere mainly by increasing the number of processed occultations (due to significant reduction in the number of L2 tracking failures) and marginally by a reduction in noise in statistics. The latter is due to the combination of reduced L2C noise (compared to L2P) and increased L1CA noise in those occultations where L2P would have failed. This result suggests application of OL tracking for L1CA and L2C signals throughout an entire occultation to optimally acquire RO data. Two methods of concurrent processing of L1CA and L2C RO signals are considered. Based on testing of individual occultations, these methods allow: (1) reduction in uncertainty of bending angles retrieved by wave optics in the lower troposphere and (2) reduction in small-scale residual errors of the ionospheric correction in the stratosphere.     

GPS/INS navigation precision and its effect on airborne radio occultation retrieval accuracy

An airborne radio occultation (RO) system has been developed to retrieve atmospheric profiles of refractivity, moisture, and temperature. The long-term objective of such a system is deployment on commercial aircraft to increase the quantity of moisture observations in flight corridors in order to improve weather forecast accuracy. However, there are several factors important to operational feasibility that have an impact on the accuracy of the airborne RO results. We investigate the effects of different types of navigation system noise on the precision of the retrieved atmospheric profiles using recordings from the GNSS Instrument System for Multistatic and Occultation Sensing (GISMOS) test flights, which used an Applanix POS/AV 510 Global Positioning System (GPS)/Inertial Navigation System (INS). The data were processed using a carrier phase differential GPS technique, and then the GPS position and inertial measurement unit data were combined in a loosely coupled integrated inertial navigation solution. This study quantifies the velocity precision as a function of distance from GPS reference network sites, the velocity precision with or without an inertial measurement unit, the impact of the quality of the inertial measurement unit, and the compromise in precision resulting from the use of real-time autonomous GPS positioning. We find that using reference stations with baseline lengths of up to 760?km from the survey area has a negligible impact on the retrieved refractivity precision. We also find that only a small bias (less than 0.5% in refractivity) results from the use of an autonomous GPS solution rather than a post-processed differential solution when used in an integrated GPS/INS system. This greatly expands the potential range of an operational airborne radio occultation system, particularly over the oceans, where observations are sparse.     

GPS radio occultation measurements on ionospheric electron density from low Earth orbit

Since the proof-of-concept GPS/Meteorology (GPS/MET) experiment successfully demonstrated active limb sounding of the Earth’s neutral atmosphere and ionosphere via GPS radio occultation (RO) from low Earth orbit, the developments of electron density (n _e) retrieval techniques and powerful processing systems have made a significant progress in recent years. In this study, the researches of n _e profiling from space-based GPS RO observations are briefly reviewed. Applying to the Formosat-3/Constellation Observing System for Meteorology, Ionosphere and Climate (FS3/COSMIC) data, we also present a compensatory Abel inversion technique including the effects of large-scale horizontal gradients and/or inhomogeneous ionospheric n _e obtained from an improved near real-time phenomenological model of the TaiWan Ionospheric Model. The results were evaluated by the ionosonde foF2 and foE data and showed improvements of rms foF2 difference from 29.2 to 16.5% in relative percentage and rms foE difference from 54.2 to 32.7% over the standard Abel inversion.     

利用FY-3C折射率对大气边界层高度的反演与分析

2013年中国发射了首颗进行全球导航卫星系统(global navigation satellite system,GNSS)掩星观测的气象卫星风云3号C星(Fengyun-3C,FY-3C),且已发布自2014年6月以来的FY-3C掩星大气产品,但目前还未见将其应用于大气边界层的相关研究。首次尝试利用FY-3C折射率产品确定边界层高度并进一步进行空间分布分析。结果表明,在小波协方差变换法基础上,进行尖锐度约束,能够确定FY-3C掩星低层大气折射率廓线中可能存在的突变,反演边界层高度。所得到的2015―2018年各年边界层高度全球分布在不同纬度及海洋和陆地上的差异基本体现了边界层与地表气候及地形的关系,但FY-3C折射率产品在低层大气的精度和垂直分辨率相对较低。因此,反演成功率总体上较低,反演结果对边界层高度空间分布细节特征的呈现仍有待提升。     

GNSS掩星探测大气的研究现状及进展

与传统的无线电探空、雷达探测等手段相比,GNSS掩星技术为大气探测提供了一个强有力的工具,其具有无校准、全天候、精度高、垂直分辨率高、全球均匀覆盖等特点。介绍了利用GNSS掩星技术获取地球大气温、压、湿等相关参数大小的研究现状。同时,提出了GNSS掩星技术在气候研究领域的发展方向,将拓宽GNSS掩星技术在全球气候变化研究中的应用。     

Baltic Sea Level project with GPS

In order to study the Baltic Sea Level change and to unify national height systems a two week GPS campaign was performed in the region in Autumn 1990. Parties from Denmark, Finland, Germany, Poland and Sweden carried out GPS measurements at 26 tide gauges along the Baltic sea and 8 VLBI and SLR fiducial stations with baseline lengths ranging from 230 km to 1600 km. The observations were processed in the network mode with the Bernese version 3.3 software using orbit improvement techniques. To get rid of the scale error introduced by the ionospheric refraction from single-frequency data, the local models of the ionosphere were estimated using L4 observations. The tropospheric zenith corrections were also considered. The preliminary results show average root mean square (RMS) errors of about ±3 cm in the horizontal position and ±7 cm in the vertical position relative to the Potsdam SLR station in ITRF89 system. After transformation of the GPS results to geoid heights using the levelled heights, an absolute comparison with gravimetric geoid heights using the least squares modification of Stokes' formula (LSMS), the modified Molodensky and the NKG Scandinavian geoid 1989 (NGK-89) models gives a standard deviation of the difference of ±7cm to ±9cm for the NKG-89 model and of ±9cm to ±30cm for the LSMS and the modified Molodensky model. The Swedish height system is found to be about 8-37cm higher than those of the other Baltic countries for NKG-89 model.     

中国区域COSMIC数据掩星反射信号特征分析

掩星信号在传播过程中能够在地球表面发生反射,其反射信号中的低层大气信息在改进天气预报精准度和气候监测等方面存在很大的应用价值.基于COSMIC(constellation observing system for meteorology,ionosphere and climate)掩星廓线资料,提取了 2011-20...     

Gravity recovery using COSMIC GPS data: application of orbital perturbation theory

 COSMIC is a joint Taiwan–US mission to study the atmosphere using the Global Positioning System (GPS) occultation technique. Improved formulas are developed for the radial, along-track, and cross-track perturbations, which are more accurate than the commonly used order-zero formulas. The formulas are used to simulate gravity recovery using the geodetic GPS data of COSMIC in the operational phase. Results show that the EGM96 model can be improved up to degree 26 using 1 year of COSMIC data. TOPEX/POSEIDON altimeter data are used to derive a temporal gravity variation. COSMIC cannot reproduce this gravity variation perfectly because of data noise and orbital configuration, but the recovered field clearly shows the gravity signature due to mass movement in an El Ni?o. Received: 3 March 2000 / Accepted: 10 November 2000     

利用FY3C与FY3D GPS掩星资料研究中低纬地区偶发E层形态的比较

偶发E层(sporadic E, Es)是主要发生在90~120 km高度的电子密度显著增强的电离层薄层,Es层的存在会导致掩星观测中全球导航卫星系统信号强度和相位的强烈波动。利用2019-01—2021-12风云三号C（Fengyun-3C,FY3C）和风云三号D(Fengyun-3D,FY3D)卫星GPS（global positioning system）掩星观测的50 Hz信噪比数据提取Es层信息,进而对两颗卫星数据分别反演得到的60°S~60°N中低纬地区Es层发生率的时空分布及季节变化进行比较。结果发现,虽然两颗卫星掩星资料得到的Es层发生率分布形态基本一致,均反映了Es层的发生率与地磁场和中性大气背景风场的相关性,但在大部分季节和地区,由FY3D得到的Es层发生率低于由FY3C得到的结果,北半球夏季中纬地区尤为明显,而FY3C反演结果与基于电离层与气候星座观测系统掩星数据的反演结果更为接近。导致差异的可能原因包括两颗卫星信噪比廓线的上边界高度分布和地方时覆盖上的差异、两颗卫星掩星接收机噪声水平的差异等。上述结果表明,后续融合两颗卫星的掩星数据进行Es层相关研究时,可能需要顾及两颗卫星接收机的不同噪声水平,在Es层发生的判定策略上进行针对性调整。     

Single frequency processing of Ørsted GPS radio occultation measurements

The Global Positioning System (GPS) radio occultation measurements obtained using the TurboRogue GPS receiver on the Danish satellite Ørsted have been processed using the single frequency method. Atmospheric profiles of refractivity and temperature are derived and validated against numerical weather prediction data from the European Centre for Medium-Range Weather Forecast (ECMWF). Results from the Ørsted GPS measurement campaign in February 2000 indicate that the single frequency method can provide retrievals with accuracy comparable to that of using two frequencies. From comparisons between measured dry temperature profiles and corresponding dry temperature profiles derived from ECMWF analysis fields, we find a mean difference of less than 0.5 K and a standard deviation of 2–4 K between 500 and 30 hPa in height. Above 30 hPa the impact of the ionosphere becomes more dominant and more difficult to eliminate using the single frequency method, and the results show degraded accuracy when compared to previous analysis results of occultation data from other missions using the dual frequency method. At latitudes less than 40° (denoted low latitudes), the standard deviation is generally smaller than at latitudes higher than 40° (denoted high latitudes). A small temperature bias is observed centered at 200 hPa for low latitudes and at 300 hPa for high latitudes. This indicates that the ECMWF analyses do not adequately resolve the tropopause temperature minimum. In the lowest part of the troposphere an observed warm bias is thought to be due to erroneous tracking of the GPS signal in cases of atmospheric multipath propagation.     

Ionospheric electron density observed by FORMOSAT-3/COSMIC over the European region and validated by ionosonde data

This research is motivated by the recent IGS Ionosphere Working Group recommendation issued at the IGS 2010 Workshop held in Newcastle, UK. This recommendation encourages studies on the evaluation of the application of COSMIC radio occultation profiles for additional IGS global ionosphere map (GIM) validation. This is because the reliability of GIMs is crucial to many geodetic applications. On the other hand, radio occultation using GPS signals has been proven to be a promising technique to retrieve accurate profiles of the ionospheric electron density with high vertical resolution on a global scale. However, systematic validation work is still needed before using this powerful technique for sounding the ionosphere on a routine basis. In this paper, we analyze the properties of the ionospheric electron density profiling retrieved from COSMIC radio occultation measurements. A comparison of radio occultation data with ground-based measurements indicates that COSMIC profiles are usually in good agreement with ionosonde profiles, both in the F2 layer peak electron density and the bottom side of the profiles. For this comparison, ionograms recorded by European ionospheric stations (DIAS network) in 2008 were used.