0 °C isotherm height for satellite communication in Malaysia

The 0 °C isotherm height is an important parameter for prediction of rain attenuation of microwave and millimeter wave for Earth-space communication. The variations of 0 °C isotherm heights for two monsoon seasons have been studied based on an analysis of radiosonde over three stations. The exceedence probability statistics of rain height are compared between the two seasons. The results on the 0 °C isotherm height can be utilized for the estimation of attenuation of microwave and millimeter wave due to rain over Earth-space paths. Attenuations of radio wave due to rain at frequencies above 10 GHz and above have also been estimated using the 0 °C isotherm height so derived.     

Analysis of rain attenuation prediction models at Ku-band in Thailand

The performance of the existing rain attenuation models in tropical zones is still a debated issue due to the lack of measurements reported from these areas of the world to develop and validate prediction models. A three-year (2003–2005) campaign of rainfall rate and rain attenuation measurements was conducted on a satellite beacon link located in a tropical region of Thailand. The cumulative distributions of rain attenuation derived from the measured data are presented and compared with those obtained with existing prediction models.     

Achievable debris orbit prediction accuracy using laser ranging data from a single station

Earlier studies have shown that an orbit prediction accuracy of 20 arc sec ground station pointing error for 1–2 day predictions was achievable for low Earth orbit (LEO) debris using two passes of debris laser ranging (DLR) data from a single station, separated by about 24 h. The accuracy was determined by comparing the predicted orbits with subsequent tracking data from the same station. This accuracy statement might be over-optimistic for other parts of orbit far away from the station. This paper presents the achievable orbit prediction accuracy using satellite laser ranging (SLR) data of Starlette and Larets under a similar data scenario as that of DLR. The SLR data is corrupted with random errors of 1 m standard deviation so that its accuracy is similar to that of DLR data. The accurate ILRS Consolidated Prediction Format orbits are used as reference to compute the orbit prediction errors. The study demonstrates that accuracy of 20 arc sec for 1–2 day predictions is achievable.     

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.     

Ku波段卫星通信雨衰预测

介绍３种计算Ｋｕ波段雨衰的数学模型,并且讨论其雨衰特性。根据３种数学模型预测各种情况下的雨衰值,并同实测数据进行比较。最后讨论我国雨区分布情况及Ｋｕ波段雨衰裕度。     

Modeling and initial assessment of the inter-frequency clock bias for COMPASS GEO satellites

The COMPASS system is a project established by China to develop an independent global satellite navigation system, which has five GEO (Geostationary Orbit) satellites and thirty Non-GEO satellites. An apparent inter-frequency clock bias (IFCB) for COMPASS GEO satellites is investigated using the real data. The bias also is modeled by the different models. Based on the 15 months (DOY 121, 2011–214, 2012) single-day-estimated results, the periodic variation of IFCBs of the COMPASS GEO satellite is studied using a harmonic analysis. The notable periods of 12 h and 8 h are noted. The harmonics-based models with different periods and different orders and quadratic function based model are used to describe the IFCB. The performances show that the 4-order harmonics-based model with the periods of 24, 12, 8 and 6 h is most optimal than others for describing the IFCB of COMPASS GEO satellite. Its amplitudes and phases estimated from a least square fit are used to study the features of the IFCB. The results show that the current amplitudes and phases do not present special features. Although the irregular amplitudes and phases of the model are disadvantageous for the long-term prediction of IFCB, it is obvious that the modeling IFCB can simple its service and a few of coefficients can replace the IFCB series. The performance of the model in short-term prediction IFCB is tested using the ten-day data (DOY 215-224, 2012).     

基于数值预报模式的Ka频段卫星系统降雨衰减短时预测

为精确地实现Ka频段星地链路降雨衰减的短时趋势和强度预测,提出了一种基于数值天气预报的降雨衰减短时预测方法,设计了降雨衰减预测流程,并应用上述方法对某地区发生的降雨过程进行研究,比较分析了降雨衰减的实测值和基于数值预报的降雨衰减预测值。试验结果表明,基于数值预报的降雨衰减和实测值变化趋势符合较好,误差较小,在20GHz频段平均误差约为1.31dB,在30GHz频段平均误差约为2.39dB,在20GHz频段平均误差约为1.31dB,在30GHz频段平均误差约为2.39dB。该方法可实现Ka频段卫星通信系统的短时降雨趋势预测。     

Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver

In order to establish a continuous GEO satellite orbit during repositioning maneuvers, a suitable maneuver force model has been established associated with an optimal orbit determination method and strategy. A continuous increasing acceleration is established by constructing a constant force that is equivalent to the pulse force, with the mass of the satellite decreasing throughout maneuver. This acceleration can be added to other accelerations, such as solar radiation, to obtain the continuous acceleration of the satellite. The orbit determination method and strategy are illuminated, with subsequent assessment of the orbit being determined and predicted accordingly. The orbit of the GEO satellite during repositioning maneuver can be determined and predicted by using C-Band pseudo-range observations of the BeiDou GEO satellite with COSPAR ID 2010-001A in 2011 and 2012. The results indicate that observations before maneuver do affect orbit determination and prediction, and should therefore be selected appropriately. A more precise orbit and prediction can be obtained compared to common short arc methods when observations starting 1 day prior the maneuver and 2 h after the maneuver are adopted in POD (Precise Orbit Determination). The achieved URE (User Range Error) under non-consideration of satellite clock errors is better than 2 m within the first 2 h after maneuver, and less than 3 m for further 2 h of orbit prediction.     

Comprehensive outage compensation of real-time orbit and clock corrections with broadcast ephemeris for ambiguity-fixed precise point positioning

The main challenge in real-time precise point positioning (PPP) is that the data outages or large time lags in receiving precise orbit and clock corrections greatly degrade the continuity and real-time performance of PPP positioning. To solve this problem, instead of directly predicting orbit and clock corrections in previous researches, this paper presents an alternative approach of generating combined corrections including orbit error, satellite clock and receiver-related error with broadcast ephemeris. Using ambiguities and satellite fractional-cycle biases (FCBs) of previous epoch and the short-term predicted tropospheric delay through linear extrapolation model (LEM), combined corrections at current epoch are retrieved and weighted with multiple reference stations, and further broadcast to user for continuous enhanced positioning during outages of orbit and clock corrections. To validate the proposed method, two reference station network with different inter-station distance from National Geodetic Survey (NGS) network are used for experiments with six different time lags (i.e., 5 s, 10 s, 15 s, 30 s, 45 s and 60 s), and one set of data collected by unmanned aerial vehicle (UAV) is also used. The performance of LEM is investigated, and the troposphere prediction accuracy of low elevation (e.g., 10–20degrees) satellites has been improved by 44.1% to 79.0%. The average accuracy of combined corrections before and after LEM is used is improved by 12.5% to 77.3%. Without LEM, an accuracy of 2–3 cm can be maintained only in case of small time lags, while the accuracies with LEM are all better than 2 cm in case of different time lags. The performance of simulated kinematic PPP at user end is assessed in terms of positioning accuracy and epoch fix rate. In case of different time lags, after LEM is used, the average accuracy in horizontal direction is better than 3 cm, and the accuracy in up direction is better than 5 cm. At the same time, the epoch fix rate has also increased to varying degrees. The results of the UAV data show that in real kinematic environment, the proposed method can still maintain a positioning accuracy of several centimeters in case of 20 s time lag.     

Integrating Non-Tidal Sea Level data from altimetry and tide gauges for coastal sea level prediction

The main objective of this paper is to integrate Non-Tidal Sea Level (NSL) from the joint TOPEX, Jason-1 and Jason-2 satellite altimetry with tide gauge data at the west and north coast of the United Kingdom for coastal sea level prediction. The temporal correlation coefficient between altimetric NSLs and tide gauge data reaches a maximum higher than 90% for each gauge. The results show that the multivariate regression approach can efficiently integrate the two types of data in the coastal waters of the area. The Multivariate Regression Model is established by integrating the along-track NSL from the joint TOPEX/Jason-1/Jason-2 altimeters with that from eleven tide gauges. The model results give a maximum hindcast skill of 0.95, which means maximum 95% of NSL variance can be explained by the model. The minimum Root Mean Square Error (RMSe) between altimetric observations and model predictions is 4.99 cm in the area. The validation of the model using Envisat satellite altimetric data gives a maximum temporal correlation coefficient of 0.96 and a minimum RMSe of 4.39 cm between altimetric observations and model predictions, respectively. The model is furthermore used to predict high frequency NSL variation (i.e., every 15 min) during a storm surge event at an independent tide gauge station at the Northeast of the UK (Aberdeen).     

A feasibility study for the detection of the diurnal variation of tropospheric NO2 over Tokyo from a geostationary orbit

We have conducted a feasibility study for the geostationary monitoring of the diurnal variation of tropospheric NO₂ over Tokyo. Using NO₂ fields from a chemical transport model, synthetic spectra were created by a radiative transfer model, SCIATRAN, for summer and winter cases. We then performed a Differential Optical Absorption Spectroscopy (DOAS) analysis to retrieve NO₂ slant column densities (SCDs), and after converting SCDs into vertical column densities (VCDs), we estimated the precision of the retrieved VCDs. The simulation showed that signal-to-noise ratio (SNR) ? 500 is needed to detect the diurnal variation and that SNR ? 1000 is needed to observe the local minimum occurring in the early afternoon (LT13–14) in summer. In winter, the detection of the diurnal variation during LT08–15 needs SNR ? 500, and SNR ? 1000 is needed if early morning (LT07) and early evening (LT16) are included. The currently discussed sensor specification for the Japanese geostationary satellite project, GMAP-Asia, which has a horizontal resolution of 10 km and a temporal resolution of 1hr, has demonstrated the performance of a precision of several percent, which is approximately corresponding to SNR = 1000–2000 during daytime and SNR ? 500 in the morning and evening. We also discuss possible biases caused by the temperature dependence of the absorption cross section utilized in the DOAS retrieval, and the effect of uncertainties of surface albedo and clouds on the estimation of precisions.     

Hardware-in-the-loop simulations of GPS-based navigation and control for satellite formation flying

A relative navigation and formation control algorithm for satellite formation flying was developed, and a hardware-in-the-loop (HIL) simulation testbed was established and configured to evaluate this algorithm. The algorithm presented is a relative navigation estimation algorithm using double-difference carrier-phase and single-difference code measurements based on the extended Kalman filter (EKF). In addition, a state-dependent Riccati equation (SDRE) technique is utilized as a nonlinear controller for the formation control problem. The state-dependent coefficient (SDC) form is formulated to include nonlinearities in the relative dynamics. To evaluate the relative navigation and control algorithms developed, a closed-loop HIL testbed is configured. To demonstrate the performance of the testbed, a test formation flying scenario comprising formation acquisition and keeping in a low earth orbit (LEO) has been established. The relative navigation results from the closed-loop simulations show that a 3D RMS of 0.07 m can be achieved for position accuracy. The targeted leader–follower formation flying in the along-track separation of 100 m was maintained with a mean position error of approximately 0.2 m and a standard deviation of 0.9 m. The simulation results show that the HIL testbed is capable of successful demonstration of the GPS-based satellite autonomous formation flying mission.     

Development and performance evaluation of tropospheric scintillation model on Ku-band satellite link over Akure,Nigeria

Existing amplitude scintillation prediction models often perform less satisfactorily when deployed outside the regions where they were formulated. This necessitates the need to evaluate the performance of scintillation models developed in one region using data data from other regions while documenting their relative errors. Due to its variation with elevation angle, frequency, other link parameters and meteorological factors, we employed three years (January 2016 to December 2018) of concurrently measured satellite radio beacons and tropospheric weather parameters to develop a location-based amplitude scintillation prediction model over the Earth-space path of Akure (7.17^oN, 5.18^oE), South-western Nigeria. The satellite beacon measurement used Tektronix Y400 NetTek Analyzer at 1 s integration time while meteorological parameters, namely; temperature, pressure and relative humidity were measured using Davis Vantage Vue weather station at 1 min integration time. Comparative study of the model’s performance with nine (9) existing scintillation prediction models indicates that the best and worst performing models, in terms of root mean square error (RMSE), are the Statistical Temperature and Refractivity (STN) and direct physical and statistical prediction (DPSP) models with values 11.48 and 51.03 respectively. Also, worst month analysis indicates that April, with respective enhancement and fade values of 0.88 and 0.90 dB for 0.01% exceedance, is the overall worst calendar month for amplitude scintillation.     

Cross testing of ionosphere models IRI-2001 and IRI-2007, data from satellite altimeters (Topex/Poseidon and Jason-1) and global ionosphere maps

In this study, we use a great body of statistical data covering the entire 23rd solar cycle to cross test data of satellite altimeters, Global Ionosphere Maps and the International Reference Ionosphere models, IRI-2001 and IRI-2007. It is revealed that experimental TEC values of the satellite altimeters regularly exceed the model ones by ∼3 TECU (1 TECU = 10¹⁶ m⁻²). The best possible value of difference between TECs obtained from altimeter and GIM-map data significantly differs for different laboratories: the maximum for CODG data falls on 2.5 TECU, ESAG – 3 TECU, JPLG – 0 TECU, UPCG – 2 TECU. The dependence of experimental and model data root-mean-square deviation on the F10.7 index is shown to be nearly linear. IRI-2001 and IRI-2007 relative errors are characterized by considerable 11-year and annual variations. Given the geomagnetic planetary index K_p under 7, IRI-2001 and IRI-2007 reproduce TEC in the ionosphere with an accuracy of ∼30% relative to measurement data from satellite altimeters. The amplitude of absolute error variations resulting from the difference in ionization enhancement between the model and the real ionosphere during the morning solar terminator transit is ∼5 TECU.     

A new model for atmospheric radiation under clear sky condition at various altitudes

Atmospheric radiation is one of the major factors that dominate the thermal behaviors of aerostats. A high-performance model is needed to evaluate the atmospheric radiation. Based on the atmospheric radiation database containing 24,862 data points compiled from 7 stations with the elevation from sea level to 2373 m and the reference code MODTRAN, a new atmospheric radiation model is proposed using regression and optimization software. It has excellent prediction accuracy with the coefficient of determination of 0.94, the root mean square error of 15.1 W/m², and the mean absolute percentage error of 4.13% for the database. Comparison with the well-known existing model shows that the new model has the highest prediction accuracy. The new model predictions agree with the MODTRAN calculations at various altitudes very well, and thus it can be used for estimating the thermal performances of a high altitude aerostat.     

A prediction model of short-term ionospheric foF2 based on AdaBoost

In this paper, the AdaBoost-BP algorithm is used to construct a new model to predict the critical frequency of the ionospheric F2-layer (foF2) one hour ahead. Different indices were used to characterize ionospheric diurnal and seasonal variations and their dependence on solar and geomagnetic activity. These indices, together with the current observed foF2 value, were input into the prediction model and the foF2 value at one hour ahead was output. We analyzed twenty-two years’ foF2 data from nine ionosonde stations in the East-Asian sector in this work. The first eleven years’ data were used as a training dataset and the second eleven years’ data were used as a testing dataset. The results show that the performance of AdaBoost-BP is better than those of BP Neural Network (BPNN), Support Vector Regression (SVR) and the IRI model. For example, the AdaBoost-BP prediction absolute error of foF2 at Irkutsk station (a middle latitude station) is 0.32 MHz, which is better than 0.34 MHz from BPNN, 0.35 MHz from SVR and also significantly outperforms the IRI model whose absolute error is 0.64 MHz. Meanwhile, AdaBoost-BP prediction absolute error at Taipei station from the low latitude is 0.78 MHz, which is better than 0.81 MHz from BPNN, 0.81 MHz from SVR and 1.37 MHz from the IRI model. Finally, the variety characteristics of the AdaBoost-BP prediction error along with seasonal variation, solar activity and latitude variation were also discussed in the paper.     

Real-time zenith tropospheric delays in support of numerical weather prediction applications

The Geodetic Observatory Pecný (GOP) routinely estimates near real-time zenith total delays (ZTD) from GPS permanent stations for assimilation in numerical weather prediction (NWP) models more than 12 years. Besides European regional, global and GPS and GLONASS solutions, we have recently developed real-time estimates aimed at supporting NWP nowcasting or severe weather event monitoring. While all previous solutions are based on data batch processing in a network mode, the real-time solution exploits real-time global orbits and clocks from the International GNSS Service (IGS) and Precise Point Positioning (PPP) processing strategy. New application G-Nut/Tefnut has been developed and real-time ZTDs have been continuously processed in the nine-month demonstration campaign (February–October, 2013) for selected 36 European and global stations. Resulting ZTDs can be characterized by mean standard deviations of 6–10 mm, but still remaining large biases up to 20 mm due to missing precise models in the software. These results fulfilled threshold requirements for the operational NWP nowcasting (i.e. 30 mm in ZTD). Since remaining ZTD biases can be effectively eliminated using the bias-reduction procedure prior to the assimilation, results are approaching the target requirements in terms of relative accuracy (i.e. 6 mm in ZTD). Real-time strategy and software are under the development and we foresee further improvements in reducing biases and in optimizing the accuracy within required timeliness. The real-time products from the International GNSS Service were found accurate and stable for supporting PPP-based tropospheric estimates for the NWP nowcasting.     

Development of one-minute rain-rate and rain-attenuation contour maps for satellite propagation system planning in a subtropical country: South Africa

Millimeter and microwave system design at higher frequencies require as input a 1-min rain-rate cumulative distribution function for estimating the level of degradation that can be encountered at such frequency bands. Owing to the lack of 1-min rain-rate data in South Africa and the availability of 5-min and hourly rainfall data, we have used rain-rate conversion models and the refined Moupfouma model to convert the available data into 1-min rain-rate statistics. The attenuation caused by these rain rates was predicted using the International Telecommunication Union (ITU) recommendations model. The Kriging interpolation method was used to draw contour maps over different percentages of time for spatial interpolation of rain-rate values into a regular grid in order to obtain a highly consistent and predictable inter-gauge rain-rate variation over South Africa. The present results will be useful for system designers of modern broadband wireless access (BWA) and high-density cell-based Ku/Ka, Q/V band satellite systems, over the desired area of coverage in order to determine the appropriate effective isotropically radiated power (EIRP) and receiver characteristics of this region.     

Geostationary orbit determination using SATRE

A new strategy of precise orbit determination (POD) for GEO (Geostationary Earth Orbit) satellite using SATRE (SAtellite Time and Ranging Equipment) is presented. Two observation modes are proposed and different channels of the same instruments are used to construct different observation modes, one mode receiving time signals from their own station and the other mode receiving time signals from each other for two stations called pairs of combined observations. Using data from such a tracking network in China, the results for both modes are compared. The precise orbit determination for the Sino-1 satellite using the data from 6 June 2005 to 13 June 2005 has been carried out in this work. The RMS (Root-Mean-Square) of observing residuals for 3-day solutions with the former mode is better than 9.1 cm. The RMS of observing residuals for 3-day solutions with the latter mode is better than 4.8 cm, much better than the former mode. Orbital overlapping (3-day orbit solution with 1-day orbit overlap) tests show that the RMS of the orbit difference for the former mode is 0.16 m in the radial direction, 0.53 m in the along-track direction, 0.97 m in the cross-track direction and 1.12 m in the 3-dimension position and the RMS of the orbit difference for the latter mode is 0.36 m in the radial direction, 0.89 m in the along-track direction, 1.18 m in the cross-track direction and 1.52 m in the 3-dimension position, almost the same as the former mode. All the experiments indicate that a meter-level accuracy of orbit determination for geostationary satellite is achievable.