Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.     

径向圆跳动误差的最小二乘评定数学模型和微机数据处理

用传统的测量方法只能得到被测零件的径向圆跳动误差的近似值。为了得到该项误差的准确值，需要研究新的测量方法。建立了径向圆跳动误差最小二乘评定数学模型，并编制了高级语言数据处理程序。在万能工具显微镜上获得了采样数据，给出了计算结果。     

An empirical model for the variation of rain drop size distribution with rain rate at a few locations in southern India

A Joss–Waldvogel impact type disdrometer was installed at four different locations in the Indian peninsula during various periods from 2001 till date. The data are analysed to study the nature of rain drop size distribution (DSD) in this region. Out of the three well known distributions that describe DSD, namely, the Marshall–Palmer, Gamma and Lognormal, it has been found that Lognormal distribution fits the DSD in this region better than the other ones. Lognormal distributions for different rain rates were then derived by fitting the lognormal function to the data using a curve fitting software. Then the variation of fit parameters with rain rate was evaluated. Incorporating these variations, into the Lognormal distribution, an empirical equation that describes the DSD in this region for different rain rates was derived. Then this equation was tested with sample data from each of these stations. The data used for validation were not used for fitting lognormal equation to derive the fit parameters. The correlation between the DSD measured and derived using the empirical model was found to be quite good (0.9) except in some cases where the coefficient dropped to 0.75. The empirical model can be updated when more data are available.     

Analysis and improvement of ionospheric thin shell model used in SBAS for China region

The most frequently used mapping function for converting slant total electron content (STEC) to vertical total electron content (VTEC) uses a simplified ionospheric single layer model in satellite based augmentation system (SBAS). Aiming at ionospheric single layer model altitude variation, we analyzed the statistical characteristic of mapping function in systemic approach particular for the region of China, and then experimental mapping function was proposed. The experimental mapping function is used in SBAS ionospheric correction and the precision of vertical correction is compared with standard mapping function. The results show that conversion factor errors depend on latitude, time and elevation, experimental mapping function is helpful to the vertical TEC estimation and may improve the precision of SBAS ionospheric refraction correction especially for the low latitude region of China.     

A new method for improving the performance of an ionospheric model developed by multi-instrument measurements based on artificial neural network

There are remarkable ionospheric discrepancies between space-borne (COSMIC) measurements and ground-based (ionosonde) observations, the discrepancies could decrease the accuracies of the ionospheric model developed by multi-source data seriously. To reduce the discrepancies between two observational systems, the peak frequency (foF2) and peak height (hmF2) derived from the COSMIC and ionosonde data are used to develop the ionospheric models by an artificial neural network (ANN) method, respectively. The averaged root-mean-square errors (RMSEs) of COSPF (COSMIC peak frequency model), COSPH (COSMIC peak height model), IONOPF (Ionosonde peak frequency model) and IONOPH (Ionosonde peak height model) are 0.58 MHz, 19.59 km, 0.92 MHz and 23.40 km, respectively. The results indicate that the discrepancies between these models are dependent on universal time, geographic latitude and seasons. The peak frequencies measured by COSMIC are generally larger than ionosonde’s observations in the nighttime or middle-latitudes with the amplitude of lower than 25%, while the averaged peak height derived from COSMIC is smaller than ionosonde’s data in the polar regions. The differences between ANN-based maps and references show that the discrepancies between two ionospheric detecting techniques are proportional to the intensity of solar radiation. Besides, a new method based on the ANN technique is proposed to reduce the discrepancies for improving ionospheric models developed by multiple measurements, the results indicate that the RMSEs of ANN models optimized by the method are 14–25% lower than the models without the application of the method. Furthermore, the ionospheric model built by the multiple measurements with the application of the method is more powerful in capturing the ionospheric dynamic physics features, such as equatorial ionization, Weddell Sea, mid-latitude summer nighttime and winter anomalies. In conclusion, the new method is significant in improving the accuracy and physical characteristics of an ionospheric model based on multi-source observations.     

Ionospheric behavior of the F2 peak parameters foF2 and hmF2 at Hainan and comparisons with IRI model predictions

Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with hourly time interval resolution for the diurnal variation, obtained with DPS-4 digisonde observations at Hainan (19.4°N, 109.0°E) are used to study the low latitude ionospheric variation behavior. The observational results are compared with the International Reference Ionospheric Model (IRI) predictions. The time period coverage of the data used for the present study is from March 2002 to February 2005. Our present study showed that: (1) In general, IRI predictions using CCIR and URSI coefficients follow well the diurnal and seasonal variation patterns of the experimental values of foF2. However, CCIR foF2 and URSI foF2 IRI predictions systematically underestimate the observed results during most time period of the day, with the percentage difference ΔfoF2 (%) values changing between about −5% and −25%, whereas for a few hours around pre-sunrise, the IRI predictions generally overestimate the observational ones with ΔfoF2 (%) sometimes reaching as large as ∼30%. The agreement between the IRI results and the observational ones is better for the year 2002 than for the other years. The best agreement between the IRI results and the observational ones is obtained in summer when using URSI coefficients, with the seasonal average values of ΔfoF2 (%) being within the limits of ±10%. (2) In general, the IRI predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the observational results. However, when using the measured M(3000)F2 as input, the diurnal variation pattern of hmF2 given by IRI2001 has a much better agreement with the observational one with the detailed fine structures including the pre-sunrise and post-sunset peaks reproduced reasonably well. The agreement between the IRI predicted hmF2 values using CCIR M(30,000)F2 option and the observational ones is worst for the afternoon to post-midnight hours for the high solar activity year 2002. During daytime hours the agreement between the hmF2 values obtained with CCIR M(30,000)F2 option and the observational ones is best for summer season. The discrepancy between the observational hmF2 and that obtained with CCIR M(30,000)F2 option stem from the CCIR M(3000)F2 model, which does not produce the small scale structures observed in the measured M(3000)F2.     

A morphological study of GPS-TEC data at Agra and their comparison with the IRI model

The temporal and seasonal variations of Total Electron Content (TEC) are studied at Agra (Geographic Lat. 27.17°N, Long. 78.89°E, Dip: 41.4°), India, which is in the equatorial anomaly region, for a period of 12 months from 01 January to 31 December, 2007 using a Global Positioning System (GPS) receiver. The mean TEC values show a minimum at 0500 h LT (LT = UT + 5.5 h) and a peak value at about 1400 h LT. The lowest TEC values are observed in winter whereas largest values are observed in equinox and summer. Anomalous variations are found during the period of magnetic disturbances. These results are compared with the TEC derived from IRI-2007 using three different options of topside electron density, NeQuick, IRI01-corr, and IRI-2001. A good agreement is found between the TEC obtained at Agra and those derived from IRI models.     

Bottom side profiles for two close stations at the southern crest of the EIA: Differences and comparison with IRI-2012 and NeQuick2 for low and high solar activity

Bottom side electron density profiles for two stations at the southern crest of the Equatorial Ionization Anomaly (EIA), São José dos Campos (23.1°S, 314.5°E, dip latitude 19.8°S; Brazil) and Tucumán (26.9°S, 294.6°E, dip latitude 14.0°S; Argentina), located at similar latitude and separated by only 20° in longitude, have been compared during equinoctial, winter and summer months under low (year 2008, minimum of the solar cycle 23/24) and high solar activity (years 2013–2014, maximum of the solar cycle 24) conditions. An analysis of parameters describing the bottom side part of the electron density profile, namely the peak electron density NmF2, the height hmF2 at which it is reached, the thickness parameter B₀ and the shape parameter B₁, is carried out. Further, a comparison of bottom side profiles and F-layer parameters with the corresponding outputs of IRI-2012 and NeQuick2 models is also reported. The variations of NmF2 at both stations reveal the absence of semi-annual anomaly for low solar activity (LSA), evidencing the anomalous activity of the last solar minimum, while those related to hmF2 show an uplift of the ionosphere for high solar activity (HSA). As expected, the EIA is particularly visible at both stations during equinox for HSA, when its strength is at maximum in the South American sector. Despite the similar latitude of the two stations upon the southern crest of the EIA, the anomaly effect is more pronounced at Tucumán than at São José dos Campos. The differences encountered between these very close stations suggest that in this sector relevant longitudinal-dependent variations could occur, with the longitudinal gradient of the Equatorial Electrojet that plays a key role to explain such differences together with the 5.8° separation in dip latitude between the two ionosondes. Furthermore at Tucumán, the daily peak value of NmF2 around 21:00 LT during equinox for HSA is in temporal coincidence with an impulsive enhancement of hmF2, showing a kind of “elastic rebound” under the action of the EIA. IRI-2012 and NeQuick2 bottom side profiles show significant deviations from ionosonde observations. In particular, both models provide a clear underestimation of the EIA strength at both stations, with more pronounced differences for Tucumán. Large discrepancies are obtained for the parameter hmF2 for HSA during daytime at São José dos Campos, where clear underestimations made by both models are observed. The shape parameter B₀ is quite well described by the IRI-2012 model, with very good agreement in particular during equinox for both stations for both LSA and HSA. On the contrary, the two models show poor agreements with ionosonde data concerning the shape parameter B₁.     

Validation of B0 and B1 in the IRI 2001 model at low solar activity for Ilorin an equatorial station

A new set of data obtained at low solar activity from Ilorin, Nigeria (geog. latitude 8.5°N, geog longitude, 4.6°E, dip 4.1°S) is used to validate the IRI 2001 model at low solar activity. The results show in general a good agreement between model and observed B0 at night but an over estimation during daytime. The overestimation is greatest during the morning period (0600LT–1000LT). The model prediction for B1 is fairly good at night and during the day. A dependence of B0 on solar zenith angle χ is observed during the daytime. A formulation of the form B0 = A[cos(χ)ⁿ] is therefore proposed. Values of the constants n and A were determined for the period of low solar activity for this station.     

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.