Using GPS-TEC data to calibrate VTEC computed with the IRI model over Nigeria

This paper presents the development of a Total Electron Content (TEC) map for the Nigerian ionosphere. In this work, TEC measurements obtained from the AFRL-SCINDA GPS (Air Force Research Laboratory-Scintillation Network Decision Aid, Global Positioning System) equipment installed at Nsukka (6.87°N, 7.38°E) are used to adapt the International Reference Ionosphere (IRI) model for the Nigerian Ionosphere. The map is being developed as a computer program (implemented in the MATLAB programming language) that shows spatial and temporal representations of TEC for the Nigerian ionosphere. The method is aimed at showing how the IRI model can be used to estimate VTEC over wide areas by incorporating GPS measurements. This method is validated by using GPS VTEC data collected from a station in Ilorin (8.50°N, 4.55°E).     

Estimation of differential code biases with Jason-2/3 onboard GPS observations

With the continuous deployment of Low Earth Orbit (LEO) satellites, the estimation of differential code biases (DCBs) based on GNSS observations from LEO has gained increasing attention. Previous studies on LEO-based DCB estimation are usually using the spherical symmetry ionosphere assumption (SSIA), in which a uniform electron density is assumed in a thick shell. In this study, we propose an approach (named the SHLEO method) to simultaneously estimate the satellite and LEO onboard receiver DCBs by modeling the distribution of the global plasmaspheric total electron content (PTEC) above the satellite orbit with a spherical harmonic (SH) function. Compared to the commonly used SSIA method, the SHLEO model improves the GPS satellite DCB estimation accuracy by 13.46% and the stability by 22.34%, respectively. Compared to the GPS satellite DCBs estimated based on the Jason-3-only observations, the accuracy and monthly stability of the satellite DCBs can be improved by 14.42% and 26.8% when both Jason-2 and Jason-3 onboard observations are jointly processed. Compared with the Jason-2 solutions, the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations have an improved consistency of better than 18.26% and 9.71% with the products provided by the Center for Orbit Determination in Europe (CODE) and Chinese Academy of Sciences (CAS). Taking the DCB products provided by the German Aerospace Center (DLR) as references, there is no improvement in accuracy of the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations than the Jason-2 solutions alone. A periodic variation is found in the time series of both the Jason-3 and Jason-2 onboard receiver DCB estimates. Preliminary analysis of the PTEC distribution based on the estimated SH coefficients are also presented.     

Predicting the probability of occurrence of spread-F over Brazil using neural networks

The probability of occurrence of spread-F can be modeled and predicted using neural networks (NNs). This paper presents a feasibility study into the development of a NN based model for the prediction of the probability of occurrence of spread-F over selected equatorial stations within the Brazilian sector. The input space included the day number (seasonal variation), hour (diurnal variation), sunspot number (measure of the solar activity), magnetic index (measure of the magnetic activity) and magnetic position. Twelve years of spread-F data from Brazil (covering the period 1978–1989) measured at the equatorial site Fortaleza (3.9°S, 38.45°W) and low latitude site Cachoeira Paulista (22.6°S, 45.0°W) are used in the development of an input space and NN architecture for the model. Spread-F data that is believed to be related to plasma bubble developments (range spread-F) was used in the development of the model. The model results show the probability of spread-F occurrence as a function of local time, season and latitude. Results from the Brazilian Sector NN (BSNN) based model are presented in this paper, as well as a comparative analysis with a Brazilian model developed for the same purpose.     

The response of the ionosphere over Ilorin to some geomagnetic storms

The effects of some geomagnetic storms on the F2 layer peak parameters over Ilorin, Nigeria (Lat. 8:53°N, Long. 4.5°E, dip angle, −2.96°) have been investigated. Our results showed that the highest intensity of the noon bite-out occurred during the March equinox and lowest during the June Solstice on quiet days. Quiet day NmF2 disturbances which appeared as a pre-storm enhancement, but not related to the magnetic storm event that followed were observed at this station. These enhancements were attributed to the modification of the equatorial electric field as a result of injection of the Auroral electric field to the low and equatorial ionosphere. For disturbed conditions, the morphology of the NmF2 on quiet days is altered. Daytime and nighttime NmF2 and hmF2 enhancements were recorded at this station. Decreases in NmF2 were also observed during the recovery periods, most of which appeared during the post-noon period, except the storm event of May 28–29. On the average, enhancements in NmF2 (i.e. Positive phases) are the prominent features of this station. Observations from this study also indicate that Dst, Ap and Kp which have been the most widely used indices in academic research in describing the behavior of geomagnetic storms, are not sufficient for storm time analysis in the equatorial and low latitude ionosphere.     

The use of ionospheric tomography and elevation masks to reduce the overall error in single-frequency GPS timing applications

Signals from Global Positioning System (GPS) satellites at the horizon or at low elevations are often excluded from a GPS solution because they experience considerable ionospheric delays and multipath effects. Their exclusion can degrade the overall satellite geometry for the calculations, resulting in greater errors; an effect known as the Dilution of Precision (DOP). In contrast, signals from high elevation satellites experience less ionospheric delays and multipath effects. The aim is to find a balance in the choice of elevation mask, to reduce the propagation delays and multipath whilst maintaining good satellite geometry, and to use tomography to correct for the ionosphere and thus improve single-frequency GPS timing accuracy. GPS data, collected from a global network of dual-frequency GPS receivers, have been used to produce four GPS timing solutions, each with a different ionospheric compensation technique. One solution uses a 4D tomographic algorithm, Multi-Instrument Data Analysis System (MIDAS), to compensate for the ionospheric delay. Maps of ionospheric electron density are produced and used to correct the single-frequency pseudorange observations. This method is compared to a dual-frequency solution and two other single-frequency solutions: one does not include any ionospheric compensation and the other uses the broadcast Klobuchar model. Data from the solar maximum year 2002 and October 2003 have been investigated to display results when the ionospheric delays are large and variable. The study focuses on Europe and results are produced for the chosen test site, VILL (Villafranca, Spain). The effects of excluding all of the GPS satellites below various elevation masks, ranging from 5° to 40°, on timing solutions for fixed (static) and mobile (moving) situations are presented. The greatest timing accuracies when using the fixed GPS receiver technique are obtained by using a 40° mask, rather than a 5° mask. The mobile GPS timing solutions are most accurate when satellites at lower elevations continue to be included: using a mask between 10° and 20°. MIDAS offers the most accurate and least variable single-frequency timing solution and accuracies to within 10 ns are achieved for fixed GPS receiver situations. Future improvements are anticipated by combining both GPS and Galileo data towards computing a timing solution.     

Near-real time monitoring of the TEC fluctuations over the northern hemisphere using GNSS permanent networks

Since the early 1990s, global positioning system measurements have been used to study of the state and rapid changes of the Total Electron Content in the ionosphere. Currently, the increasing number of permanent stations makes it possible to generate maps of the irregularities in the ionosphere for specified regions with sub-daily resolution. The main goal of this work was to apply global navigation satellite system observations to obtain information about ionospheric variability around the North Geomagnetic Pole. In order to detect the ionospheric disturbances, 30-s observation data was used. The Rate of Total Electron Content Index was applied as a measure of the variability in the ionosphere. The first analyses were executed using more than 100 permanent stations. The results show two kinds of products: 2-hour maps in spherical geomagnetic coordinates and daily maps presenting the occurrence of the strong Total Electron Content fluctuations as a magnetic local time function, for the most disturbed days of April 2010. Apart from the main product of the algorithm, the Rate of Total Electron Content time series for individual satellite tracks was presented. The results demonstrated very good sensitivity of the obtained maps, which can detect even quite weak disturbances. The presented algorithm developed at the Geodynamic Research Laboratory of the University of Warmia and Mazury, in cooperation with Institute of Terrestrial Magnetism, Ionosphere and Radiowave Propagation, will be applied in the near future to create near-real time service of the conditions in the ionosphere based on the Global Navigation Satellite Systems observations.     

Using the radial basis function neural network to predict ionospheric critical frequency of F2 layer over Wuhan

Neural networks (NNs) have been applied to ionospheric predictions recently. This paper uses radial basis function neural network (RBF-NN) to forecast hourly values of the ionospheric F₂ layer critical frequency(f_oF₂), over Wuhan (30.5N, 114.3E), China. The false nearest neighbor method is used to determine the embedding dimension, and the principal component analysis (PCA) is used to reduce noise and dimension. The whole study is based on a sample of about 26,000 observations of f_oF₂ with 1-h time resolution, derived during the period from January 1981 to December 1983. The performance of RBF-NN is estimated by calculating the normalized root-mean-squared (NRMSE) error, and its results show that short-term predictions of f_oF₂ are improved.     

Pre-storm behaviour of NmF2 and TEC (GPS) over equatorial and low latitude stations in the Indian sector

The pre-storm behavior of NmF₂ and TEC over an equatorial station, Trivandrum (8.47°N, 76.91°E, dip 0.6°S) and a low latitude station, Waltair (17.7°N,83.3°E, dip 20°N) has been studied for a total of 18 strong geomagnetic storms with DST ? −100 nT. The simultaneous measurements of GPS-TEC and NmF₂ over Trivandrum and Waltair during the period 2000–2005 have been considered for the present study. It is found that there is a substantial increase in NmF₂ and TEC before the onset of the storm over Waltair, while the increase is not present at Trivandrum. The origin of pre-storm enhancements in electron density still remains unresolved owing to several conditions in their potential sources and occurrence mechanisms. In the present study an attempt is made to identify the possible mechanisms responsible for such enhancements in electron density of the F-region.     

The variability and predictability of the IRI B0, B1 parameters over Grahamstown,South Africa

The International Reference Ionosphere (IRI) parameters B₀ and B₁ provide a representation of the thickness and shape, respectively, of the F2 layer of the bottomside ionosphere. These parameters can be derived from electron density profiles that are determined from vertical incidence ionograms. This paper aims to illustrate the variability of these parameters for a single mid latitude station and demonstrate the ability of the Neural Network (NN) modeling technique for developing a predictive model for these parameters. Grahamstown, South Africa (33.3°S, 26.5°E) was chosen as the mid latitude station used in this study and the B₀ and B₁ parameters for an 11 year period were determined from electron density profiles recorded at that station with a University of Massachusetts Lowell Center for Atmospheric Research (UMLCAR) Digisonde. A preliminary single station NN model was then developed using the Grahamstown data from 1996 to 2005 as a training database, and input parameters known to affect the behaviour of the F2 layer, such as day number, hour, solar and magnetic indices. An analysis of the diurnal, seasonal and solar variations of these parameters was undertaken for the years 2000, 2005 and 2006 using hourly monthly median values. Comparisons between the values derived from measured data and those predicted using the two available IRI-2001 methods (IRI tables and Gulyaeva, T. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7(6), 39–48, 1987.) and the newly developed NN model are also shown in this paper. The preliminary NN model showed that it is feasible to use the NN technique to develop a prediction tool for the IRI thickness and shape parameters and first results from this model reveal that for the mid latitude location used in this study the NN model provides a more accurate prediction than the current IRI model options.     

The effect of total solar eclipse of October 3, 2005, on the total electron content over Europe

GPS observations from EUREF permanent GPS network were used to observe the response of TEC (Total Electron Content) to the total solar eclipse on October 3, 2005, under quiet geomagnetic conditions of the daytime ionosphere. The effect of the eclipse was detected in diurnal variations and more distinctly in the variations of TEC along individual satellite passes. The trough-like variations with a gradual decrease and followed by an increase of TEC at the time of the eclipse were observed over a large region. The depression of TEC amounted to 3–4 TECU. The maximum depression was observed over all stations located at the maximum path of the solar eclipse. The delay of a minimum level of TEC with respect to the maximum phase of the eclipse was about 20–30 min.     

Application of satellite observations to study the changes of hypoxic conditions in near-bottom water in the western part of Peter the Great Bay (the Sea of Japan)

In this paper we explore the possibilities of applying satellite ocean colour (OC) observations and SST to study the changes in the conditions of hypoxia in the near-bottom water in the western part of Peter the Great Bay. Near-bottom water hypoxia occurs in water bodies with increased organic matter influx when the dissolved oxygen (DO) consumed at its oxidation is not restored. Consumption of most DO is usually attributed to the oxidation of organic matter formed as a result of increased algae growth during water eutrophication. Satellite data on indicators of phytoplankton (chlorophyll-a concentration (Chl) and fluorescence (FLH)) allow to analyze the spatial-temporal changes of this substation. Coloured dissolved organic matter (CDOM), non-algal particles (NAP) influence on satellite Chl estimates and also on near-bottom water hypoxia formation. This study analyzes daily, seasonal, and inter-annual changes in the distributions of indicators (Chl, FLH, the coefficients of light absorption by coloured detrital matter (a_CDM) and light backscattering by suspended particles (b_bp)), based on the instant satellite OC data from MODIS-Aqua. Data on the Chl, the sea surface temperature (SST) from the MODIS-Aqua, the precipitation from the TRMM satellite and the hydrometeorological stations (HMSs), the wind speed and direction from HMS “Vladivostok” are used to study the influence of hydrometeorological conditions on the Chl values. These distributions were compared with the literary information based on field observations of the hypoxia cases in the same area and with the changes in the vertical DO, Chl, temperature, salinity distribution obtained by coastal expeditions in October-November 2010 and February-March 2011. Significant interrelations within 95% confidence level between the satellite Chl, FLH values calculated at the MUMM atmospheric correction and in situ Chl values obtained in the autumn of 2010 were reached separately for the cases with winds of northern and southern directions with the correlation coefficients of 0.71, 0.48 and 0.49, 0.71, respectively. Significant dependences of Chl on SST and Chl on wind speed explained by the influence of continental runoff and water ventilation were obtained. Therefore, the changes of Chl reflect the changes of hypoxic conditions in the near-bottom water. In Amursky Bay the onset of hypoxia was at the Chl and SST values equal to 4 mg m^?3 and 13 °C (↑ – at increasing SST); near Furugelm Island it was at 1.6 mg m^?3 and 25 °C (↑), 1 mg m^?3 and 21 °C (↓). The difference in the Chl values was reflected in the hypoxia onset timings that were the beginning of June (2011), August (2013), and September (2014), respectively. The water flow from the eastern coast of Amursky Bay in early August of 2013 recorded from the OC and SST satellite imagers appeared in an additional hypoxic zone. Decreased OC characteristics in the runoff of the Razdolnaya River in August-September of 2014 were a sign of hypoxia at its mouth. Near Furugelm Island the hypoxia destruction (increase in the DO level from 1 to 4.5 ml L^?1) was observed at the Chl of 0.9 mg m^?3 and SST = 18 °C (↓). At the autumn maximum of Chl equal to 1.7 mg m^?3 and SST = 4 °C (↓) in mid-November the DO level here increased to 8 ml L^?1. In Amursky Bay, short-term destructions/weakening of hypoxia manifested themselves in sharp increases of Chl. At that, the ratio between the Chl value and the approximation level was equal to 2 and higher for SST equal to 22–25 °C (↑), to 0.9 and higher for SST equal to 5–13 °C (↓). With the water stratification destruction in temperature and the noticeable weakening of the stratification in salinity (mid-November), the hypoxia destructed (the DO level increased from 2 to 6 ml L^?1). In this case, Chl and SST were about 3 mg m^?3 and 5 °C (↓).     

Performance assessment of real-time orbit determination for the Haiyang-2D using onboard BDS-3/GPS observations

The Global Navigation Satellite System (GNSS) receivers equipped on the Haiyang-2D (HY-2D) satellite is capable of tracking the signals of both the third generation of BeiDou satellite navigation System (BDS-3) and the Global Positioning System (GPS), which make it feasible to assess the performance of real-time orbit determination (RTOD) for the HY-2D using onboard GNSS observations. In this study, the achievable accuracy and convergence time of RTOD for the HY-2D using onboard BDS-3 and GPS observations are analyzed. Benefiting from the binary-offset-carrier (BOC) modulation, the BDS-3 C1X signal includes less noise than the GPS C1C signal, which has the same signal frequency and chipping rate. The root mean squares (RMS) of the noises of C1X and C1C code measurements are 0.579 m and 1.636 m, respectively. Thanks to a ten-times higher chipping rate, the code measurements of BDS-3 C5P, GPS C1W and C2W are less noisy. The RMS of code noises of BDS-3 C5P, GPS C1W, and C2W are 0.044 m, 0.386 m, and 0.272 m, respectively. For the HY-2D orbit, the three-dimensional (3D) and radial accuracies can reach 31.8 cm and 7.5 cm with only BDS-3 observations, around 50 % better than the corresponding accuracies with GPS. Better performance of the BDS-3 in RTOD for the HY-2D is attributed to the high quality of its broadcast ephemeris. When random parameters are used to absorb ephemeris errors, substantial improvement is seen in the accuracy of HY-2D orbit with either BDS-3 or GPS. The 3D RMS of HY-2D orbit errors with BDS-3 and GPS are enhanced to 23.1 cm and 33.6 cm, and the RMS of the radial components are improved to 6.1 cm and 13.3 cm, respectively. The convergence time is 41.6 and 75.5 min for the RTOD with BDS-3 and GPS, while it is reduced to 39.2 and 27.4 min after the broadcast ephemeris errors are absorbed by random parameters. Overall, the achievable accuracy of RTOD with BDS-3 reaches decimeter level, which is even better than that with GPS, making real-time navigation using onboard BDS-3 observations a feasible choice for future remote sensing missions.     

Detection of Envisat RA2/ICE-1 retracked radar altimetry bias over the Amazon basin rivers using GPS

Altimetry is now routinely used to monitor stage variations over rivers, including in the Amazon basin. It is desirable for hydrologic studies to be able to combine altimetry from different satellite missions with other hydrogeodesy datasets such as leveled gauges and watershed topography. One requirement is to accurately determine altimetry bias, which could be different for river studies from the altimetry calibrated for deep ocean or lake applications. In this study, we estimate the bias in the Envisat ranges derived from the ICE-1 waveform retracking, which are nowadays widely used in hydrologic applications. As a reference, we use an extensive dataset of altitudes of gauge zeros measured by GPS collocated at the gauges. The thirty-nine gauges are spread along the major tributaries of the Amazon basin. The methodology consists in jointly modeling the vertical bias and spatial and temporal slope variations between altimetry series located upstream and downstream of each gauge. The resulting bias of the Envisat ICE-1 retracked altimetry over rivers is 1.044 ± 0.212 m, revealing a significant departure from other Envisat calibrations or from the Jason-2 ICE-1 calibration.     

Characteristics of recent dust storms over the Indian region using real time multi-satellite observations from the direct broadcast receiving system at IMD

In this study, observations from microwave satellites, visible and infrared instruments have been analyzed to detect dust storm over north and north–west part of India during 18–23 March 2012. This study investigated the approach to utilize the multi satellite data of Moderate Resolution Imaging Spectroradiometer (MODIS) on-board the Terra and Aqua satellite and the Advanced Microwave Sounding Unit (AMSU) on-board NOAA satellite to study the characteristics of dust storms from real time direct broadcast (DB) receiving system installed at three places of India Meteorological Department (IMD). The dust storm detection is based on the infrared brightness temperature (BT) difference between channels at 11 and 12 μm and polarized BT difference between two channels of 89 and 23.8 GHz. It is found that the significant differences between the BT of channel 89 and 23.8 can be used as a discriminator of identifying dust storm. The Total Ozone Mapping Spectroradiometer (TOMS) Aerosol Index (AI) and AMSU-A 23 GHz channel BT from NOAA satellite over the north and north-west part of India have also been analyzed. The result indicated the characteristic behavior between BT and AI during the different phases of the dust storm. Finally, the occurrence of dust outbreaks has also been validated with sky radiometer of IMD, which confirms the presence of a dust storm over the Indian region. Further, the findings of the study and its approaches apply to the other dust storm cases which occurred during the months of April and June 2012. The integrated approach suggested the potential to use high resolution data of microwave as well as thermal–infrared using multi-satellite observations from real time direct broadcast system for the detection of severe, moderate or weak dust storms very well. The approach is found to be promising for operational application.     

Short term variations of total electron content (TEC) fitting to a regional GPS network over the Kingdom of Saudi Arabia (KSA)

The ionosphere is a dispersive medium for radio waves with the refractive index which is a function of frequency and total electron content (TEC). TEC has a strong diurnal variation in addition to monthly, seasonal and solar cycle variations and small and large scale irregularities. Dual frequency GPS observations can be utilized to obtain TEC and investigate its spatial and temporal variations. We here studied short term TEC variations over the Kingdom of Saudi Arabia (KSA). A regional GPS network is formed consisting of 16 sites in and around KSA. GPS observations, acquired between 1st and 11th February 2009, were processed on a daily basis by using the Bernese v5.0 software and IGS final products. The geometry-free zero difference smoothed code observables were used to obtain two hour interval snapshots of TEC and their RMS errors at 0.5 × 0.5 degree grid nodes and regional ionosphere models in a spherical harmonics expansion to degree and order six. The equatorial ionized anomaly (EIA) is recovered in the south of 20°N from 08:00 to 12:00 UT. We found that day-by-day TEC variation is more stable than the night time variation.     

The impact of different background errors in the assimilation of satellite radiances and in-situ observational data using WRFDA for three rainfall events over Iran

Satellite radiances and in-situ observations are assimilated through Weather Research and Forecasting Data Assimilation (WRFDA) system into Advanced Research WRF (ARW) model over Iran and its neighboring area. Domain specific background error based on x and y components of wind speed (UV) control variables is calculated for WRFDA system and some sensitivity experiments are carried out to compare the impact of global background error and the domain specific background errors, both on the precipitation and 2-m temperature forecasts over Iran. Three precipitation events that occurred over the country during January, September and October 2014 are simulated in three different experiments and the results for precipitation and 2-m temperature are verified against the verifying surface observations. Results show that using domain specific background error improves 2-m temperature and 24-h accumulated precipitation forecasts consistently, while global background error may even degrade the forecasts compared to the experiments without data assimilation. The improvement in 2-m temperature is more evident during the first forecast hours and decreases significantly as the forecast length increases.     

The response of sporadic E-layer to the total solar eclipse of July 22, 2009 over the equatorial ionization anomaly region of the Indian zone

The digital ionosonde located in Bhopal (23.2°N, 77.2°E), India has been used to investigate the responses of the Es layer in the equatorial ionization anomaly (EIA) crest to the total solar eclipse (TSE) of July 22, 2009. Results show the presence of intense Es layer during and after the eclipse period. The gravity waves induced by the solar eclipse propagated upward in the Es layer and produced the periodic disturbance. The results of the wavelet analysis display the presence of dominant oscillation of about 24–32, 16–20 and 8 min. The appearance of intense sporadic-E concomitantly with the signatures of gravity wave suggests that the wind shear introduced by the solar eclipse induced gravity wave might be the plausible mechanism behind the intensification of Es-layer ionization.