Abstract

The - model has become widely used in statistical analyses of radio channels, due to the flexibility provided by its two degrees of freedom. Among several applications, it has been used in the characterization of low-latitude amplitude scintillation, which frequently occurs during the nighttime of particular seasons of high solar flux years, affecting radio signals that propagate through the ionosphere. Depending on temporal and spatial distributions, ionospheric scintillation may cause availability and precision problems to users of global navigation satellite systems. The present work initially stresses the importance of the flexibility provided by - model in comparison with the limitations of a single-parameter distribution for the representation of first-order statistics of amplitude scintillation. Next, it focuses on the statistical evaluation of the power spectral density of ionospheric amplitude scintillation. The formulation based on the - model is developed and validated using experimental data obtained in São José dos Campos (23.1°S; 45.8°W; dip latitude 17.3°S), Brazil, located near the southern crest of the ionospheric equatorial ionization anomaly. These data were collected between December 2001 and January 2002, a period of high solar flux conditions. The results show that the proposed model fits power spectral densities estimated from field data quite well.

1. Introduction

The presence of random irregularities in the ionospheric plasma density may cause amplitude and phase scintillation on signals of the global navigation satellite systems (GNSS). These effects may be strong enough to cause cycle slips, loss of signal tracking, and, consequently, degradation in the performance and positioning accuracy of ground-based receivers. Thus, all technological infrastructures that depend on GNSS signals for precise positioning and navigation can be directly affected.

Ionospheric scintillation at GNSS frequencies is a phenomenon that typically occurs after sunset and mainly during premidnight hours. Its occurrence is predominant at the polar and low-latitude regions. The present contribution will focus on scintillation which occurs within the latitude belt of ±20° centered at the geomagnetic equator. Climatological studies of equatorial and low-latitude ionospheric scintillation have shown combined seasonal and longitudinal dependence. For example, in the Brazilian longitudinal sector, the ionospheric irregularity activity is more intense near December solstice months [1]. This behavior departs from that in the western sector of South America, where the scintillation occurrence peaks during the equinoctial period. These studies also indicated that scintillation is more intense during solar maximum years (when the ionospheric layer is thicker and denser).

Adequate knowledge of the temporal properties and fading time scales of amplitude scintillation can be very useful in the development of appropriated mitigation tools of current positioning/navigation applications based on GNSS signals. Fremouw et al. [2] demonstrated that Nakagami- probability density function (PDF) describes well the variability of the amplitude scintillations. Later, Hegarty et al. [3] developed an amplitude scintillation signal model based on a Nakagami- distribution. In [4] it was presented that either the Nakagami- or the Rice PDF could describe well the distribution of amplitudes that were present in their measurements. More recently, Moraes et al. [5] showed that the ionospheric amplitude scintillation phenomenon can be modeled by using first-order statistics of the - distribution proposed by Yacoub [6]. Moraes et al. [7] also represented the level crossing rate (LCR) and average fading duration (AFD) due to equatorial ionospheric scintillation using the same distribution. This model explores the nonlinearity of the propagation medium, associating the physical fading phenomena with the and parameters. Moraes et al. [5, 7] empirically tested and parameterized this distribution as a function of scintillation severity, showing that this model provides a better fit to experimental data than those yielded by previous ones. Indeed, having two parameters that are physically described, instead of just one as in the previous fading models (such as those based on the Nakagami- distribution), makes the - distribution more flexible, assuring a better agreement with scintillation data. Additionally, they proposed a series of approximations to the estimation of the and parameters tailored to the severity of the ionospheric scintillation events.

In the communications field, the power spectral density (PSD) of the random refractive-index fluctuations in the propagation medium is an important feature for channel characterizations and performance evaluations. PSD characteristics are directly related to size, behavior, speed, and duration of ionospheric irregularities [810]. The PSD is widely used to study and model the ionospheric plasma irregularity processes. For example, by using these PSDs, it is possible to estimate the drift velocity and the average height of the irregularities in the ionosphere. The literature also lists a substantial number of works dealing with the PSDs of VHF and UHF transionospheric signals [1115]. The temporal variations in the received signal caused by the ionosphere have been analyzed using PSD models by many works in the literature [4, 16, 17]. The spectrum widths of received signals affected by scintillation provide important information that may be considered in the design of the bandwidths of tracking loops of GNSS receivers [16].

In continuation of previous work by Moraes et al. [5, 7, 18], the present study presents the experimental validation of the - model for the PSD of transionospheric signals received on the ground, by comparing results from the formulation of Dias and Yacoub [19] with those from real scintillation data.

The remainder of the present paper is organized in the following way. Section 2 describes the measurements and the initial processing of the data to estimate some fundamental parameters. Section 3 summarizes the benefits of the first-order statistical characterization of scintillation using the - distribution. Section 4 discusses the power spectral density formulation based on the - distribution and presents the obtained results from the model validation. Finally, Section 5 presents concluding remarks.

2. Measurements

The scintillation monitor (SCINTMON) used in the campaign was developed at Cornell University [20] with a 12-channel correlator and built using a GEC-Plessey GPS card with special firmware, which was designed to maintain lock even under extreme scintillation scenarios. The monitor tracks signals from up to 11 GPS satellites, simultaneously. One channel is dedicated to noise floor estimation. Its front end works on the L1 band (1575.42 MHz) and provides measurements of wide band power and wide band noise with a 50 Hz sampling rate. Unfortunately, SCINTMON does not provide carrier phase measurements. More information about SCINTMON can be found in [21].

The scintillation monitor was installed at the Instituto Nacional de Pesquisas Espaciais (INPE) headquarters in São José dos Campos (SJC), Brazil (geographic coordinates: 23.2°S, 45.9°W, dip latitude: 17.5°S), near the southern peaks of the equatorial anomaly, where the world’s strongest scintillation occurs during high solar flux years. The measurements were made between December 14, 2001, and January 14, 2002. During this period, the average sunspot number and the 10.7 cm solar flux density were 126 and 176 × 10−22 W/m2/Hz, respectively. This period is part of the equatorial spread F season in Brazil, which starts around September and lasts until April [22]. Scintillation observations started at 18:00 LT and continued until the next day at 06:00 LT. Scintillation at GHz frequencies is normally observed between sunset hours and local midnight, but cases of postmidnight scintillation have also been observed. In this study, only data from GPS satellites with elevation greater than 30° were considered, to avoid signal contamination by nongeophysical sources (such as tropospheric scattering and multipath). The geophysical conditions under which the measurements were made allow the study of broad range of scintillation levels. From the set of observations (over 370 hours during 32 days of measurements), approximately 102 hours of significant GPS L1 amplitude scintillation were obtained.

In Figure 1(a), the measurement site SJC is overlaid on the global map of the predicted maximum electron density (NmF2) of the ionosphere at 20:00 LT on December 31, 2001 (day of the year DOY = 365). This map, generated using the international reference ionosphere (IRI) model [23], clearly shows the enhancement of NmF2 in the equatorial anomaly region. Indeed, strong scintillation is known to occur more frequently during high solar flux conditions than during moderate or low solar flux conditions [24]. Figure 1(b) shows the azimuth-elevation track of satellite PRN 18 over SJC during the evening hours on December 15, 2001. The scintillation index will be defined in the next section. The time series of the ratio  (dBHz) between the received power and the noise power density corresponding to the link characterized by Figure 1(b) is shown in Figure 1(c), respectively. Additional details about the selected datasets and their preprocessing can be found in Moraes et al. [25].

2.1. Processing of the Scintillation Data

The strength of the amplitude scintillation is characterized by the index , defined as the normalized standard deviation of the received signal intensity (power). That is, it is given by [26] where is the intensity, is the amplitude of received signal, and the angular brackets denote an ensemble average (time average of one-minute records).

The index is computed by SCINTIMON using the wide band power and noise wide band power measurements, as well as their respective low-pass filtered versions and in a slightly less straightforward way. The filtered versions are obtained using a 6th order low-pass Butterworth filter with a cutoff frequency of 0.1 Hz [21]. The series of and are then used to estimate the real signal strength variance over a one-minute period following the approach described by [21]: where is the total number of samples during each 60-second period. The respective mean signal power is given by The computed index is finally given by Figure 1(d) shows the results from the calculations of the scintillation index for every one-minute record of the signal displayed in Figure 1(c).

Additionally, the normalized signal amplitude scintillation is obtained by [23] The autocorrelation coefficient of the normalized signal amplitude scintillation is given by where denotes the expected value operator (again implemented by a time average) and and are the mean and variance values of . Another important parameter used to characterize scintillation is the decorrelation time , defined as the time lag at which the autocorrelation function falls off by from its maximum (zero lag) value [4]: While the index is an indicator of the strength of amplitude fading, the decorrelation time is an indicator of the rapidity of the fades.

The two examples in Figures 2(a) and 2(b) illustrate the variability in the amplitude scintillation patterns and the decorrelation time estimated from the measurements made during the campaign. The two cases display approximately the same index (0.9), but very distinct values, 0.94 s and 0.18 s, respectively. Figure 2(c) shows the respective autocorrelation coefficients, , for both cases and as consequence their temporal differences.

3. The - Model for Scintillation

The - distribution is a general fading model proposed by Yacoub [6]. The use of this distribution for modeling ionospheric amplitude scintillation was proposed by Moraes et al. [5]. Assuming that the average signal intensity is equal to one, the - probability density function of the normalized amplitude envelope of the received signal is given by where is the Gamma function of the argument .

The - parameters can be estimated based on the equality that involves the moments of - envelope given by As suggested by Yacoub [6], the left hand side of (9) can be obtained from field data for arbitrarily selected values of . Assuming, for instance, and , a system with two nonlinear equations for the two desired unknowns and is established and can be solved. This is the ideal approach when real data are available. Alternatively, assuming in (9) and combining the result with (1), it is possible to establish the following relation between the scintillation index and the parameters of the - distribution:

The single-parameter Nakagami- distribution has been used for a long time as the model that best described the amplitude scintillation phenomenon [2, 16]. It can be obtained from the - distribution with and . Note that each value of the scintillation index provides a unique Nakagami- distribution. However, there are an infinite number of and values that satisfy (10) for each value of . Equation (10) is convenient, because it indicates that, differently from the Nakagami- distribution, the - model may describe different patterns of scintillation for the same scintillation index . This is especially interesting for strong scintillation, where the index alone is not considered a proper indication of the ionospheric perturbations. Indeed, it saturates to a value near unity under increasingly strong scattering conditions, in addition to being independent of the rate of signal fading [27].

To stress the importance of the flexibility provided by - model in comparison with the limitations of the Nakagami- distribution in this particular aspect, Figures 3(a), 3(b), and 3(c) show three examples of different scintillation patterns with approximately the same scintillation index (≈0.90). These cases, which have not been presented previously, were selected from the data base described in Section 2 [25]. Figure 3(d) displays the following: (i) the empirical PDF for each of these records, using color-coded symbols; (ii) the corresponding - PDFs, with parameters estimated by the application of the procedure proposed in the text immediately following (9) to each case; and (iii) the Nakagami- PDF (with ). As shown in Figure 3(d), the differences among the PDFs of these received signals, captured by the - model, cannot be displayed by the single-parameter distribution. Indeed, the curves in Figure 3(d) vary substantially for the same value, that is, for the same value on the left hand side of (10). This flexibility of - model, resulting from the existence of two parameters instead of just one, provides the capacity of a better fit to scintillation data. It is worth mentioning that the distributions tend to spread as increases, occupying the lower region of intensity values with higher probabilities. This means that the increased values, for a fixed value, represent a more severe scenario for propagation, with higher occurrence of deep fading of the received signal.

4. Power Spectral Density

The exact autocorrelation function for the - fading signal is given by [19] where for ; is the Gauss hypergeometric function with parameters , , and and argument ; and is the temporal autocorrelation coefficient. The literature provides several models of autocorrelation coefficients that attempt to describe different communications scenarios with their relative scattering conditions [28]. Mason [29] analyzed the model based on a second-order Butterworth filter, expressed by where and is the decorrelation time lag [4].

The power spectral density for the - envelope is the Fourier transform of (11). According to Dias and Yacoub [19], a closed-form expression for the PSD of the - envelope does not seem to exist. As a consequence, they proposed and validated the following approximation for (11) that only retains the first two lowest-order terms of the series expansion of the Gauss hypergeometric function [30]: By taking the Fourier transform of (13), the PSD can be written as where is the Dirac delta function and is the Fourier transform of the autocorrelation coefficients . It is important to note that (14) is transformed into the Nakagami- PSD of Dias et al. [31] for . Analyses of a limited set of representations for by Humphreys et al. [4] and Moraes et al. [32] performed in the frequency and time domains, respectively, verified that the Butterworth model provided the best spectral shape for ionospheric scintillation studies. The second-order Butterworth PSD is expressed by where is the 3 dB cutoff frequency and is the slope factor [4]. Figure 4 shows three examples of the Nakagami- (dashed lines) and - (solid lines) PSDs represented by (14), for different values of the scintillation index , assuming  s and thus  Hz. More information on and its typical values can be found in Moraes et al. [25].

5. Validation

In this section, the - power spectral density formulation represented by (14) is compared with corresponding estimations from one-minute records of real scintillation data using the well-known Welch method [33].

For calculations, the decorrelation time lag , which is the basic parameter of (15), is determined from the estimated autocorrelation coefficient according to (6). The values of , , and are estimated for each one-minute record as explained in Sections 2 and 3. The metric adopted by the validation of the power spectral density formulation is based on Carrano et al. [34], defined by where is the estimated PSD and is obtained from (14) and (15). Only the sections of the power spectral density located between  Hz and  Hz are used in the comparison. The value for was chosen with basis on the cutoff frequency of 0.1 Hz, according to the scintillation processing described in Section 2.1 [21]. The value for was selected to avoid any possible contamination of the results from the analysis by the noise floor of the estimated PSDs.

Figure 5 shows examples of empirical and - PSDs for different values of . Figure 6 displays two additional examples that also include the results from the Nakagami- formulation. A good agreement between empirical and calculated spectra is observed in all cases for , both visually and by the small values obtained for the parameter . This agreement is kept until the noise-dominated frequency band is reached. It is interesting to observe in Figure 6 that the differences between the two formulations are small. Note that only the positive halves of the empirical and calculated spectra have been plotted, since they are symmetrical around zero. However, all the PSD values corresponding to nonzero frequencies have been multiplied by two, to conserve the full power of the spectra. This procedure does not affect the calculations indicated by (16).

The results from the application of the validation procedure to the full data set are summarized in Table 1. The second line of the table displays the number of one-minute records in each of the classes of scintillation index   . In the next two lines, Table 1 displays the average error parameter for the - and Nakagami- formulations, which are larger than the ones displayed in Figures 5 and 6 for all classes of . It is observed that, on average, the - formulation provides slightly better fits to empirical PSDs than those provided by the Nakagami- model. This result was expected from the discussion in Section 3 and the studies of the second-order statistics of level crossing rate (LCR) and average fading duration (AFD) presented by Moraes et al. [5, 7]. Additionally, this table also shows the average values of the cutoff frequency and frequency for each class of scintillation index. The spectral broadening as increases, previously observed by other authors [16], is evident in the last two lines of Table 1.

Another important result from the above analysis is the performance of Nakagami- spectrum formulation in comparison with that of the - model. During the validation process, it was noted that the difference between both spectra is less than 1 dB. For spectral analysis purposes only, without taking first-order statistics considerations into account, the Nakagami- model proved to be fairly accurate.

6. Concluding Remarks

The ionosphere is known to affect Earth-space communications and GNSS applications operating in the L band and below. One of the effects on transionospheric radio waves is associated with amplitude scintillation caused by the presence of random plasma density irregularities. Thus, it is important to develop increasingly accurate statistical models to describe the temporal properties of amplitude scintillation.

The present work initially stressed the importance of the flexibility provided by - model in comparison with the limitations of the Nakagami- distribution for the representation of first-order statistics of scintillation.

Next, the - and its special case, Nakagami-, were validated for second-order statistics of amplitude scintillation by comparison of the associated power spectral density formulations with estimations from real data. The validation results show that the characterizations of the PSD of amplitude scintillation by both models are in good agreement with the experimental estimates. The - model is marginally better than that provided by the Nakagami- formulation. However, it should be observed that these models may be highly dependent on the assumed autocorrelation coefficient model.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful to Professor Fabiano Rodrigues (University of Texas at Dallas, TX, USA) for the discussions concerning the application of the - model for scintillation. Kelias Oliveira is grateful to Instituto Federal de Educação, Ciência e Tecnologia de Goiás, for supporting his doctoral studies at ITA through Award DINTER ITA/IFG 23038.044846/2009-76. Alison de Oliveira Moraes wishes to thank Instituto de Aeronáutica e Espaço (IAE), where he works as a Research Engineer, for assisting his cooperation research with ICEA, INPE, and ITA. Emanoel Costa was supported by AFOSR Award FA9550-12-1-0031 during the present work. Eurico Rodrigues de Paula is grateful for the partial supports from AFOSR FA9550-10-1-0564 and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) 305684/2010-8 Grants. Marcio Tadeu de Assis Honorato Muella would like to thank the support from CNPq under Process no. 308017/2011-0.