Noise analysis of continuous GPS coordinate time series for CMONOC

The Crustal Movement Observation Network of China (CMONOC) is one of the major scientific infrastructures, mainly using Global Positioning System (GPS) measurements, to monitor crustal deformation in the Chinese mainland. In this paper, decade-long coordinate time series of 26 continuous GPS sites of CMONOC are analyzed for their noise content using maximum likelihood estimation (MLE). We study the noise properties of continuous GPS time series of CMONOC for the unfiltered, filtered solutions and also the common mode signals in terms of power law plus white noise model. In the spatial filtering, we remove for every time series a common mode error that was estimated from a modified stacking of position residuals from other sites within ∼1000 km of the selected site. We find that the common mode signal in our network has a combination of spatially correlated flicker noise and a common white noise with large spatial extent. We demonstrate that for the unfiltered solutions of CMONOC continuous GPS sites the main colored noise is a flicker process, with a mean spectral index of ∼1. For the filtered solutions, the main colored noise is a general power law process, indicating that a major number of the filtered regional solutions have a combination of noise sources or local effects. The velocity uncertainties from CMONOC continuous GPS coordinate time series may be underestimated by factors of 8–16 if a pure white noise model is assumed. In addition, using a white plus flicker noise model, the median values of velocity errors for the unfiltered solutions are 0.16 (north), 0.17 (east) and 0.58 (vertical) mm/yr, and the median values for the filtered solutions are 0.09 (north), 0.10 (east) and 0.40 (vertical) mm/yr.     

Medium-scale traveling ionospheric disturbances (MSTID) modeling using a dense German GPS network

The MSTIDs are wave-like perturbations of the ionospheric plasma, which cause the most common ionospheric disturbances in mid-latitude regions. Generally the MSTIDs have velocities of several hundred meters per second and wavelengths of several hundred kilometers. The wave-like effect of the MSTID is one of the main obstacles for accurate interpolation of ionospheric corrections in a medium-scale reference GPS network. In this paper we show a new method of detecting and modeling MSTIDs using dense German GPS network. The between-epoch single difference ionospheric delays from a medium scale dense GPS network are used to estimate the parameter of the MSTID e.g. amplitude, wavelength and velocity. The efficiency of the approach is tested with data from about 320 GPS stations in and near Germany. A MSTID wave moving from east to west across Germany was observed at September 27 in 2009. Its wavelength is about 302 km, with a period of ∼7 min and velocity of about 700 m/s.     

Coronal mass ejections acceleration problem

A statistical study of acceleration and its error of coronal mass ejections (CMEs) observed by the Large Angle Spectrometric Coronagraph (LASCO) is performed. A total of 5594 CMEs events have been analyzed by using a least-square method and using the error in the height measures. We verify that slower CMEs (velocities in the interval from 200 to 500 km s⁻¹) tend to have a positive acceleration (about 1 m s⁻²) at heights above 5 solar radii, while less than 10% CMEs show an average negative acceleration (about −2.2 m s⁻²) as they propagate from 5 to 30 solar radii. For most individual CMEs one can not say if they are accelerated or decelerated, only for 8% of all observed CMEs events one can extract the sign of the acceleration in the 5–30 solar radii.     

Propagation of normal and faster CMEs in the interplanetary medium

We have analyzed 101 Coronal Mass Ejection (CME) events and their associated interplanetary CMEs (ICMEs) and interplanetary (IP) shocks observed during the period 1997–2005 from the list given by Mujiber Rahman et al. (2012). The aim of the present work is to correlate the interplanetary parameters such as, the speeds of IP shocks and ICMEs, CME transit time and their relation with CME parameters near the Sun. Mainly, a group of 10 faster CME events (V_INT > 2200 km/s) are compared with a list of 91 normal events of Manoharan et al. (2004). From the distribution diagrams of CME, ICME and IP shock speeds, we note that a large number of events tends to narrow towards the ambient (i.e., background) solar wind speed (∼500 km/s) in agreement with the literature. Also, we found that the IP shock speed and the average ICME speed measured at 1 AU are well correlated. In addition, the IP shock speed is found to be slightly higher than the ICME speed. While the normal events show CME travel time in the range of ∼40–80 h with a mean value of 65 h, the faster events have lower transit time with a mean value of 40 h. The effect of solar wind drag is studied using the correlation of CME acceleration with interplanetary (IP) acceleration and with other parameters of ICMEs. While the mean acceleration values of normal and faster CMEs in the LASCO FOV are 1 m/s², 18 m/s², they are −1.5 m/s² and −14 m/s² in the interplanetary medium, respectively. The relation between CME speed and IP acceleration for normal and faster events are found to agree with that of  and  except slight deviations for the faster events. It is also seen that the faster events with less travel time face higher negative acceleration (>−10 m/s²) in the interplanetary medium up to 1 AU.     

Ambiguity resolution in precise point positioning with hourly data for global single receiver

Integer ambiguity resolution (IAR) can improve precise point positioning (PPP) performance significantly. IAR for PPP became a highlight topic in global positioning system (GPS) community in recent years. More and more researchers focus on this issue. Progress has been made in the latest years. In this paper, we aim at investigating and demonstrating the performance of a global zero-differenced (ZD) PPP IAR service for GPS users by providing routine ZD uncalibrated fractional offsets (UFOs) for wide-lane and narrow-lane. Data sets from all IGS stations collected on DOY 1, 100, 200 and 300 of 2010 are used to validate and demonstrate this global service. Static experiment results show that an accuracy better than 1 cm in horizontal and 1–2 cm in vertical could be achieved in ambiguity-fixed PPP solution with only hourly data. Compared with PPP float solution, an average improvement reaches 58.2% in east, 28.3% in north and 23.8% in vertical for all tested stations. Results of kinematic experiments show that the RMS of kinematic PPP solutions can be improved from 21.6, 16.6 and 37.7 mm to 12.2, 13.3 and 34.3 mm for the fixed solutions in the east, north and vertical components, respectively. Both static and kinematic experiments show that wide-lane and narrow-lane UFO products of all satellites can be generated and provided in a routine way accompanying satellite orbit and clock products for the PPP user anywhere around the world, to obtain accurate and reliable ambiguity-fixed PPP solutions.     

The effect of geocenter motion on Jason-2 orbits and the mean sea level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards ( Lemoine et al., 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC’s SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a – 4.67 ± 3.40 mm error in the Z-component of the orbit frame which creates 1.06 ± 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.     

基于经验加速度的低轨卫星轨道预报新方法

研究将定轨过程中的经验加速度应用于地球低轨卫星轨道预报的新方法. 利用GPS伪距观测数据和简化动力学最小二乘批处理方法对地球低轨卫星定 轨, 其中卫星位置、速度及大气阻力系数和辐射光压系数可以直接用于轨道预报. 作为简化动力学最重要特征的经验加速度呈现准周期、余弦曲线特点, 可通过 傅里叶级数拟合建模. 确定性动力学模型与补偿大气阻力模型误差的切向经验 加速度级数拟合模型组成增强型动力学模型用于提高轨道预报精度. 应用 GRACE-A星载GPS伪距观测数据和IGS超快星历定轨并进行轨道预报, 结果表明 轨道预报初值位置精度达到0.2m, 速度精度达到1×10-4m·s-1, 预报3天位置精度优于60m, 比只利用确定性动力学模型进行预报精度平 均提高2.3倍. 先定轨后预报的模式可用在星上自主精确导航系统中.      

GNSS geodetic techniques for time and frequency transfer applications

We performed an initial analysis of the pseudorange data of the GIOVE-B satellite, one of the two experimental Galileo satellites currently in operation, for time transfer.¹ For this specific aim, software was developed to process the GIOVE-B raw pseudoranges and broadcast navigation messages collected by the Galileo Experimental Sensor Stations (GESS) tracking network, yielding station clock phase errors with respect to the Experimental Galileo System Time (EGST). The software also allows processing the Global Positioning System (GPS) P1 and P2 pseudorange data with broadcast navigation message collected at the same stations to obtain the station clock phase errors with respect to the GPS system time (GPST). Differencing these solutions between stations provides two independent means of GNSS time transfer. We compared these time transfer results with Precise Point Positioning (PPP) method applied to GPS data in combined carrier-phase and pseudorange mode as well as in pseudorange-only mode to show their relative merits. The PPP solutions in combined carrier-phase and pseudorange mode showed the least instability of the methods tested herein at all scales, at few parts in 10¹⁵ at 1 day for the stations processed, following a tau^−½ interval dependency. Conversely, the PPP solutions in pseudorange-only mode are an order of magnitude worst (few parts in 10¹⁴ at 1 day for the stations processed) following a tau⁻¹ power-law, but slightly better than the single-satellite raw GPS time transfer solutions obtained using the developed software, since the PPP least-squares solution effectively averages the pseudorange noise. The pseudorange noise levels estimated from PPP pseudorange residuals and from clock solution comparisons are largely consistent, providing a validation of our software operation. The raw GIOVE-B time transfer, as implemented in this work, proves to be slightly better than single-satellite raw GPS satellite time transfer, at least in the medium term. However, one of the processed stations shows a combined GPS P1 and P2 pseudorange noise level at 2 m, a factor 2 worst than usually seen for geodetic receivers, so the GPS time transfer results may not be at their best for the cases processed. Over the short term, the GPS single-satellite time transfer instability outperforms the GIOVE-B by an order of magnitude at 1 s interval, which would be due to the different characteristics of the tracking loop filters for GPS P1 and P2 on one hand and the GIOVE-B signals on the other. Even at this preliminary stage and using an experimental satellite system, results show that the GIOVE-B (and hence Galileo) signals offer interesting perspectives for high precision time transfer between metrological laboratories.     

GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

Ionospheric differential error determination using ray tracing for a short baseline

Since the United States government discontinued Selective Availability (SA) on 1 May 2000, ionospheric effects have been responsible for the largest errors in GPS systems. The standard Differential GPS (DGPS) method is incapable of completely eliminating the ionospheric error. This paper describes a new approach to determine the differential ionospheric error between geographically distributed receiver stations. The ray paths of GPS signals were simulated using a modified Jones 3D ray tracing programme that includes the effect of the geomagnetic field. A Nelder–Mead optimisation algorithm was embedded in the program to precisely determine the satellite-to-station path. A realistic ionospheric model is essential for accurate ray tracing results and for estimates of differential error that are accurate on sub-centimetre scales. Here, the ionospheric model used in the ray tracing programme was developed by fitting realistic ionosphere profiles with a number of exponential functions. Results were compared to the theoretical approach. Results show that the differential delay is about 1–5 cm at low elevation angles for a short baseline of 10 km, as reported in other literature. This delay is often neglected in DGPS application. The differential delay also shows a pattern similar to that predicted by the Klobuchar model. The method proposed here can be used to improve future GPS applications.     

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.     

Onboard orbit determination using GPS observations based on the unscented Kalman filter

Spaceborne GPS receivers are used for real-time navigation by most low Earth orbit (LEO) satellites. In general, the position and velocity accuracy of GPS navigation solutions without a dynamic filter are 25 m (1σ) and 0.5 m/s (1σ), respectively. However, GPS navigation solutions, which consist of position, velocity, and GPS receiver clock bias, have many abnormal excursions from the normal error range for space operation. These excursions lessen the accuracy of attitude control and onboard time synchronization. In this research, a new onboard orbit determination algorithm designed with the unscented Kalman filter (UKF) was developed to improve the performance. Because the UKF is able to obtain the posterior mean and covariance accurately by using the second-order Taylor series expansion through the sampled sigma points that are propagated by using the true nonlinear system, its performance can be better than that of the extended Kalman filter (EKF), which uses the linearized state transition matrix to predict the covariance. The dynamic models for orbit propagation applied perturbations due to the 40 × 40 geo-potential, the gravity of the Sun and Moon, solar radiation pressure, and atmospheric drag. The 7(8)th-order Runge–Kutta numerical integration was applied for orbit propagation. Two types of observations, navigation solutions and C/A code pseudorange, can be used at the user’s discretion. The performances of the onboard orbit determination were verified using real GPS data of the CHAMP and KOMPSAT-2 satellites. The results of the orbit determination were compared with the precision orbit ephemeris (POE) of the CHAMP and KOMPSAT-2 satellites.     

DORIS and GPS monitoring of the Gavdos calibration site in Crete

Due to its specific geographical location as well as its geodetic equipment (DORIS, GNSS, microwave transponder and tide gauges), the Gavdos station in Crete, Greece is one of the very few sites around the world used for satellite altimetry calibration. To investigate the quality of the Gavdos geodetic coordinates and velocities, we analyzed and compared here DORIS and GPS-derived results obtained during several years of observations. The DORIS solution is the latest ignwd11 solution at IGN, expressed in ITRF2008, while the GPS solution was obtained using the GAMIT software package. Current results show that 1–2 mm/yr agreement can be obtained for 3-D velocity, showing a good agreement with current geophysical models. In particular, the agreement obtained for the vertical velocity is around 0.3–0.4 mm/yr, depending on the terrestrial reference frame. As a by-product of these geodetic GPS and DORIS results, Zenith Tropospheric Delays (ZTDs) estimations were also compared in 2010 between these two techniques, and compared to ECMWF values, showing a 6.6 mm agreement in dispersion without any significant difference between GPS and DORIS (with a 97.6% correlation), but with a 13–14 mm agreement in dispersion when comparing to ECMWF model (with only about 90% correlation for both techniques). These tropospheric delay estimations could also provide an external calibration of the tropospheric correction used for the geophysical data of satellite altimetry missions.     

Synergistic use of multitemporal RAMP,ICESat and GPS to construct an accurate DEM of the Larsemann Hills region,Antarctica

An enhanced digital elevation model (DEM) of the Larsemann Hills region, east Antarctica, is constructed synergistically by using highly accurate ground-based GPS measurements, satellite-derived laser altimetry (GLAS/ICESat) and Radarsat Antarctic Mapping Project (RAMPv2) DEM-based point elevation dataset. Our DEM has a vertical accuracy of about 1.5 times better than RAMPv2 DEM and seven times better than GLAS/ICESAT-based DEM. The accuracy is improved by validating the RAMPv2 DEM elevation by supplementing with GLAS/ICESat and DGPS survey data, when compared to that of DEM constructed by using GLAS/ICESat or RAMPv2 alone. With the use of accurate GPS data as ground control points reference elevations, the DEM extracted is much more accurate with least mean RMSE of 34.5 m than that constructed by using a combination of GLAS/ICESat and RAMPv2 as true reference. The newly constructed DEM 7 achieves highest accuracy with the least average elevation difference of 0.27 m calculated using 46 ground reference points. Available DEMs of Antarctic region generated by using radar altimetry and the Antarctic Digital Database indicate elevation variations in the range of 50–100 m, which necessitates the generation of local DEM and its validation by using ground truth. This is our first attempt of fusing multi-temporal, multi-sensor and multi-source elevation data to generate a DEM of any part of Antarctica, in order to address the ice elevation change to infer the ice mass balance. Our approach focuses on the strengths of each elevation data source to produce an accurate DEM.     

The spatial distributions of the sources of UV solar Explosive Events at different velocities

We analyze five solar Explosive Events observed in the Si iv emission line at 139.37 nm during SUMER/SoHO raster scans near disk center. The Doppler velocities from −45 km s⁻¹ to +45 km s⁻¹ were sorted into eleven symmetrically organized velocity bins, five equally-sized bins on each side and one bin for the line at rest. The radiance values along the NS oriented slit around the maximum of each EE are used to build one-dimensional distributions. We study these distributions and their development in space and time for different velocity bins to unveil the 3D-structure and evolution of Explosive Events. The spatial radiance distributions for the EE maximum and around it were made. For some Explosive Events the dispersion direction is not oriented orthogonally to the slit direction; consequently, the blueshifted emission is in a linear way offset from the redshifted emission by several pixels. The largest offset (∼5″) is observed between distributions whose velocities differ by ∼36 km s⁻¹. The linear relation agrees with an scenario of flows inside magnetic arcs.     

Mid-latitude thermospheric zonal winds during the equinoxes

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

Generating precise and homogeneous orbits for Jason-1 and Jason-2

Driven by the GMES (Global Monitoring for Environment and Security) and GGOS (Global Geodetic Observing System) initiatives the user community has a strong demand for high-quality altimetry products. In order to derive such high-quality altimetry products, precise orbits for the altimetry satellites are a necessity. With the launch of the TOPEX/Poseidon mission in 1992 a still on-going time series of high-accuracy altimetry measurements of ocean topography started, continued by the altimetry missions Jason-1 in 2001 and Jason-2/OSTM in 2008. This paper contributes to the on-going orbit reprocessing carried out by several groups and presents the efforts of the Navigation Support Office at ESA/ESOC using its NAPEOS software for the generation of precise and homogeneous orbits referring to the same reference frame for the altimetry satellites Jason-1 and Jason-2. Data of all three tracking instruments on-board the satellites (beside the altimeter), i.e. GPS, DORIS, and SLR measurements, were used in a combined data analysis. About 7 years of Jason-1 data and more than 1 year of Jason-2 data were processed. Our processing strategy is close to the GDR-C standards. However, we estimated slightly different scaling factors for the solar radiation pressure model of 0.96 and 0.98 for Jason-1 and Jason-2, respectively. We used 30 s sampled GPS data and introduced 30 s satellite clocks stemming from ESOC’s reprocessing of the combined GPS/GLONASS IGS solution. We present the orbit determination results, focusing on the benefits of adding GPS data to the solution. The fully combined solution was found to give the best orbit results. We reach a post-fit RMS of the GPS phase observation residuals of 6 mm for Jason-1 and 7 mm for Jason-2. The DORIS post-fit residuals clearly benefit from using GPS data in addition, as the DORIS data editing improves. The DORIS observation RMS for the fully combined solution is with 3.5 mm and 3.4 mm, respectively, 0.3 mm better than for the DORIS-SLR solution. Our orbit solution agrees well with external solutions from other analysis centers, as CNES, LCA, and JPL. The orbit differences between our fully combined orbits and the CNES GDR-C orbits are of about 0.8 cm for Jason-1 and at 0.9 cm for Jason-2 in the radial direction. In the cross-track component we observe a clear improvement when adding GPS data to the POD process. The 3D-RMS of the orbit differences reveals a good orbit consistency at 2.7 cm and 2.9 cm for Jason-1 and Jason-2. Our resulting orbit series for both Jason satellites refer to the ITRF2005 reference frame and are provided in sp3 file format on our ftp server.     

Coronal radio-sounding detection of a CME during the 1997 Galileo solar conjunction

Frequency fluctuations of the Galileo S-band radio signal were recorded nearly continuously during the spacecraft’s solar conjunction from December 1996 to February 1997. A strong propagating disturbance, most probably associated with a coronal mass ejection (CME), was detected on 7 February when the radio ray path proximate point was on the west solar limb at about 54 solar radii from the Sun. The CME passage through the line of sight is characterized by a significant increase in the fluctuation intensity of the recorded frequency and by an increase in the plasma speed from about 234 km s⁻¹ up to about 755 km s⁻¹. These velocity estimates are obtained from a correlation analysis of frequency fluctuations recorded simultaneously at two widely-separated ground stations. The density turbulence power spectrum is found to steepen behind the CME front. The Galileo radio-sounding data are compared with SOHO/LASCO observations of the CME in the corona and with WIND spacecraft data near the Earth’s orbit.     

Advanced architectures for real-time Delay-Doppler Map GNSS-reflectometers: The GPS reflectometer instrument for PAU (griPAU)

In recent years Global Navigation Satellite System’s signals Reflectometry (GNSS-R) has stood as a potential powerful remote sensing technique to derive scientifically relevant geophysical parameters such as ocean altimetry, sea state or soil moisture. This has brought out the need of designing and implementing appropriate receivers in order to track and process this kind of signals in real-time to avoid the storage of huge volumes of raw data. This paper presents the architecture and performance of the Global Positioning System (GPS) Reflectometer Instrument for PAU (griPAU), a real-time high resolution Delay-Doppler Map reflectometer, operating at the GPS L1 frequency with the C/A codes. The griPAU instrument computes 24 × 32 complex points DDMs with configurable resolution (Δf_Dmin = 20 Hz, Δτ_min = 0.05 chips) and selectable coherent (minimum = 1 ms, maximum = 100 ms for correlation loss Δρ < 90%) and incoherent integration times (minimum of one coherent integration period and maximum not limited but typically <1 s). A high sensitivity (DDM peak relative error = 0.9% and DDM volume relative error = 0.03% @ T_i = 1 s) and stability (Δρ/Δt = −1 s⁻¹) have been achieved by means of advanced digital design techniques.     

Inertial aided BDS triple-frequency integer ambiguity rounding method

The integer ambiguity resolution (AR) of carrier phase is significant for Global Navigation Satellite System (GNSS) precise positioning. However, in kinematic case, single-epoch AR methods based on alone GNSS are usually not reliable due to the instable pseudorange accuracy. Moreover, the computation of classical AR method Least Squares Ambiguity Decorrelation Adjustment (LAMBDA) is large. Thus, the inertial measurement unit (IMU) is introduced, a new inertial-aided AR method that directly rounds the float ambiguity of BeiDou triple-frequency combined observations, which is characterized by long wavelength, low carrier-phase noise and ionospheric delay, is proposed. The mathematical model of the new method is derived first. Then the impacts of the carrier-phase noise, ionospheric delay and inertial navigation system (INS) position error on the AR success ratio of combined observation are analyzed through probabilistic approach. Based on above investigation, the combinations (0, ?1, 1), (1, 4, ?5) and (4, ?2, ?3) are selected to resolve the original ambiguity. A vehicular integrated navigation test is performed to demonstrate the proposed method. The results show that the average AR success ratios of the three selected combinations, whose float ambiguity errors are 0.041, 0.146, 0.279 cycle respectively, are above 97.25% without regard to low-elevation C05. With respect to positioning accuracy based on our AR method when compared with IE software, the east, north, up error RMS of position are 0.042, 0.024, 0.069 m, respectively. In terms of the AR recover after the BeiDou signals outage, as long as 62 s BeiDou signal complete outage, all the ambiguities of all satellites could be re-fixed immediately. Besides, during the 90 s signals partial outage, the AR is not influenced by the position error, since the float ambiguity errors are all below half-cycle. The research of this contribution demonstrates the effectiveness of the proposed new method, which indicates it is applicable to kinematic positioning, even in BDS degraded and denied environments.