Research Article | Open Access
M. V. Tinin, "Influence of Ionospheric Irregularities on GNSS Remote Sensing", Advances in Meteorology, vol. 2015, Article ID 532015, 10 pages, 2015. https://doi.org/10.1155/2015/532015
Influence of Ionospheric Irregularities on GNSS Remote Sensing
We have used numerical simulation to study the effects of ionospheric irregularities on accuracy of global navigation satellite system (GNSS) measurements, using ionosphere-free (in atmospheric research) and geometry-free (in ionospheric research) dual-frequency phase combinations. It is known that elimination of these effects from multifrequency GNSS measurements is handi-capped by diffraction effects during signal propagation through turbulent ionospheric plasma with the inner scale being smaller than the Fresnel radius. We demonstrated the possibility of reducing the residual ionospheric error in dual-frequency GNSS remote sensing in ionosphere-free combination by Fresnel inversion. The inversion parameter, the distance to the virtual screen, may be selected from the minimum of amplitude fluctuations. This suggests the possibility of improving the accuracy of GNSS remote sensing in meteorology. In the study of ionospheric disturbances with the aid of geometry-free combination, the Fresnel inversion eliminates only the third-order error. To eliminate the random TEC component which, like the measured average TEC, is the first-order correction, we should use temporal filtering (averaging).
At present, remote sensing methods relying on ground-based and low-orbit observations of signals from global navigation satellite systems (GNSS) such as GLONASS and GPS are finding ever-widening application in environmental research [1–8]. These methods allow remote diagnostics of wide sectors of overground space, which makes them attractive for use especially in the Arctic region. However, the spatially and temporally variable ionosphere significantly affects the accuracy and reliability of GNSS remote sensing. Therefore, it is necessary to study effects of ionospheric irregularities on GNSS measurements and to find possibilities of their elimination. Unfortunately, the ionospheric variability impedes elimination of ionospheric errors through modeling. Particularly, the ionospheric error reduction with the aid of the well-known Klobuchar model [9, 10] does not exceed 50% . A more precise approach makes use of the frequency dependence of ionospheric error. Here we will restrict our consideration to phase measurements as being more precise versus measurements of time delay of signal (code measurements). Let us have a phase path measured taking into account the whole number of wavelengths determined by methods for phase ambiguity resolution (see, e.g., ). Then the first geometrical optics approximation yields [9, 10]Here the pseudo distance comprises, besides the distance between a satellite and a ground-based observer, all other frequency-independent corrections (errors opposite in sign) to phase path (e.g., delay in the troposphere); consideris the electron density integral along a ground-based-receiver-satellite line, which is called total electron content (TEC); is the radiation frequency. The first-order correction, the second term in (1) proportional to , describes ionospheric plasma effects on GNSS signal propagation, ignoring diffraction effects, geomagnetic field effect (second-order correction), and ray bending (the third-order correction). Here we use a generally accepted (see, e.g., [13–15]) notation of the -order correction for the term in observation equations (1) and below, which is proportional to .
GNSS remote sensing methods (see, e.g., [1–8]) rest on the relation of terms on the right-hand side of (1) with environmental parameters. The distance included into depends on the observation point. Therefore, if the observer is on the Earth’s surface, motions of this surface can be seen in variations. On the other hand, the atmospheric phase delay of signal included into is determined by the refractive index which depends on atmospheric temperature and density and vapor pressure. Thus, measurements can provide parameters important for meteorology. Estimating the second term in (1) allows us to acquire data on ionospheric plasma disturbances from the spatial and temporal TEC behavior.
In order to find the pseudo distance and TEC from (1), we need to perform dual-frequency measurements. These measurements give the following system:By solving system (3) for two unknowns and , we obtain the result of the dual-frequency measurements:Equations (4) and (5) are so-called ionosphere-free and geometry-free observation combinations, respectively .
The dual-frequency measurements from formula (4) reduce the ionospheric distance-measurement error by 99% . In practice, this is usually quite sufficient. However, in atmospheric investigations, in studies of crust movements during earthquakes and volcanic activity and other precisional GNSS measurements, this decimetric accuracy is not sufficient. Therefore, many researchers have been exploring the possibility of improving the accuracy of dual-frequency measurements.
The second-order correction which takes into account geomagnetic effects on the ionospheric refractive index is ignored in dual-frequency measurements (3). This correction has been examined in many studies (see, e.g., [1, 13–15, 18–21]). It has been revealed that in the thin layer approximation whose error is 1-2 mm this correction is expressed through the first-order one [13, 15, 22]. Therefore, it can be taken into account by calculating with the familiar geomagnetic field model [1, 13–15] or by considering directly in dual-frequency measurements through modification of the ionosphere-free combination [22, 23].
After accounting for the second-order residual error, further accuracy improvement is associated with the consideration of the third-order error. In the geometrical optics approximation, this error is largely attributed to the path deviation from the straight line and can be eliminated in triple-frequency measurements [23, 24]. Yet, like dual-frequency measurements, triple-frequency measurements yield a system of three observation equations from which we can find , first-order correction (TEC), and third-order correction (ray bending). This can only be done, however, when scales of ionospheric irregularities exceed the Fresnel radius, and the geometrical optics approximation is applicable. Unfortunately, scales of ionospheric irregularities range widely from tens of meters to tens of kilometers [25–27]. In this case, the lower bound (so-called inner scale) of the ionospheric turbulent spectrum may be less than the Fresnel radius which varies with the elevation angle in the 200–700 m range (see, e.g., ). As a result, in triple-frequency measurements, as compared to dual-frequency ones, only the average error decreases, whereas the random error variance increases [23, 28].
The above results of studies of GNSS ionospheric errors largely refer to the results obtained with ionosphere-free combination (4), that is, during GNSS measurements of frequency-independent characteristics such as crust movement and tropospheric parameters. At present, however, GNSS measurements are widely used for studying traveling ionospheric disturbances as well as the ionosphere’s behavior during earthquakes and magnetic storms [7, 29]. Yet TEC measurements (2) obtained from geometry-free combination (5) are utilized. It is obvious that in the study of this large-scale structure of the ionosphere and its dynamics, a significant role in formation of the measurement error will be played by ionospheric irregularities of less scale. In this paper, therefore, we will present results of investigations into residual ionospheric errors of dual-frequency GNSS measurements for both combinations (4) and (5).
2. Ionospheric Effects in Single-Frequency GNSS Measurements
Consider the behavior of GNSS signal in the form of harmonic wave with the time dependence . Since the GNSS satellite is far from the ionosphere, our task reduces to the study of incidence of the plane-wave with ( is the wavelength) on an inhomogeneous anisotropic ionosphere. Let the -axis of our Cartesian coordinate system pass through the source point and the point located in the vicinity of the receiving point . At GNSS frequencies, sizes of ionospheric irregularities usually exceed With condition (6) fulfilled, when solving the wave equation, we use the second Rytov approximation and obtain the following solution [30, 31]: whereis the complex phase whose first and second approximations have the forms:Because the geomagnetic field effect on the velocity of GNSS signal propagation in the ionosphere can be taken into account, say, in the modified dual-frequency reception [22, 23], we do not discuss it here. In (9)is the 2D spectrum of electron density of ionospheric plasma andThe imaginary part of the complex phase determines the GNSS signal phase which, with the point of observation on the axis (), isThe real part of the complex phase is associated with the level Equations (12), (8), and (9) for the phase path yield the observation equation :where The second term on the right-hand side in (14) is the first-order correction which, unlike (2), accounts for diffraction effects; the third term is the diffraction variant of the third-order correction. When the minimum size of irregularities exceeds the Fresnel radius (14) and (15) give the geometrical optics approximation [23, 31]whereIn this case, unlike first approximation formula (1), besides the first-order correction proportional to the total electron content , the third-order correction (the final term on the right-hand side in (17)) associated with the path deviation from the straight line in the inhomogeneous ionosphere is taken into consideration. When ignoring the third-order correction, (17) yields (1) from which follow formulas (4) and (5) of dual-frequency measurements.
Thus, the errors to be considered are associated in the geometrical optics approximation with the GNSS ray bending. Moreover, it is necessary to correctly describe the interaction of the wave field with the irregularities for which conditions of the geometrical optics approximation (16) are not fulfilled. Expression (14), derived in the second Rytov approximation, accounts both for the ray bending and for the diffraction effects.
The electron density of ionospheric plasma is normally represented as sum of two componentsHere is the average electron density, that is, density of the background ionosphere and its medium-scale disturbances being of the utmost interest for ionosphere researchers [7, 29]; is the random component of electron density associated with small-scale irregularities.
The electron density of the large-scale background ionosphere and its medium-scale disturbances largely fulfills the condition of validity of the geometrical optics approximation (16) at GNSS frequencies. Besides, ray bending effects in the background ionosphere are sufficiently small [13, 24, 31]. Therefore, diffraction effects in GNSS measurements are associated first and foremost with random small-scale irregularities caused by ionospheric plasma turbulence.
Random ionospheric irregularities are usually represented by a random, quasi-homogeneous, normally distributed field with a given spectrum related to the correlation function by the relation The spectrum of turbulence irregularities is often as follows :where is the dispersion of electron density fluctuations, , .
Spectrum (21) is bounded by the outer scale in the order of several tens of kilometers and by the inner scale m.
By substituting (19) in (15), we can express characteristics of the random () through corresponding characteristics of the random field . For spectrum (21) the average and the correlation functions are determined by the following relations :where ,
Note that contrary to the second term, the first term on the right-hand side in (25) accounts for the regular-refraction influence on the random component of the third-order error. According to the numerical simulation, the contribution of the first term is much less than that of the second term .
3. Ionospheric Effects in Dual-Frequency GNSS Measurements
By substituting (14) in geometry-free combination (5), we get where is, according to (22), the mean first-order TEC value,is the random component of the first-order TEC, and is the third-order TEC correction.
Now it is easy to find statistical characteristics of dual-frequency measurement errors. For measurements with ionosphere-free combination (4), the average correction is equal to the average third-order correction :where the average is determined by (23). The correction variance of this combination is where the first-order correction variance is and the third-order correction variance isThe correlation functions and are defined by (24) and (25)–(28), respectively.
For measurements with geometry-free combination (5), the average correction isVariance of this correction iswhere is the first-order correction variance andis the third-order correction variance.
To assess probable errors in dual-frequency measurements, we have performed numerical simulation with (35)–(42). For the model of background ionospheric layer we took the Chapman layerwhere and Hz is the layer critical frequency; km is the Earth’s radius and km is the characteristic scale of the layer; km is the height of its maximum. Such parameters corresponded to a vertical electron content of 81 TECU ( m2). Random ionospheric irregularities were specified by spectrum (21). Yet we took values for the outer scale km, for the index , and for the coordinate dependence of rms deviation of electron density fluctuation we assumed the proportionality to the background electron density : . In the calculations, we set , and for radiation frequencies we took GPS frequencies Hz, Hz.
Figures 1 and 2 illustrate elevation angle dependence of average (a) and standard deviations (b) of corrections of first (dashed line) and third (solid line) orders at inner sizes of 50 m (blue line), 200 m (green line), and 1000 m (red line), using the ionosphere-free (Figure 1) and geometry-free (Figure 2) combinations.
A characteristic property of absolute average corrections and standard deviations is their increase with decreasing elevation angle (see Figures 1 and 2) associated with the wave-path extension in the inhomogeneous layer. What stands out here is (see Figure 1(a)) the weak dependence of the average correction in ionosphere-free combination (4) on the inner scale . This is due to the fact that the greatest contribution to this correction, as is clear from (23), is made by the irregularities located at a distance where . In this area, the dependence of the integral in (23) on is offset by the factor before the integral. Figure 1(b) shows that the standard deviation of the third-order correction (solid line) in measurements can exceed the standard deviation of the first-order correction (dashed line). The weak dependence of the standard deviation of the first-order (dotted lines in Figure 2(b)) on the inner scale is not associated with diffraction effects, as in the case of the average correction in Figure 1(a). Here we deal with phase variance that, as is known , is generally determined by large-scale irregularities and slightly depends on the high-frequency part of irregularity spectrum which is associated with the inner scale . The results of Figure 2 suggest that a probable error in TEC measurements may reach 1 TECU. Third-order corrections increase with decreasing inner sizes, and, as their averages (Figure 2(a)) are other than zero, they cannot be eliminated after time averaging. This, of course, must be taken into account in studying features of ionospheric disturbances of the same order via GNSS. The results of Figure 1 also suggest probable errors in measurements to be of order of several centimeters, which is comparable with effects measured with GNSS in atmosphere .
In , numerical simulation was used to explore the possibility of improving the accuracy of GNSS measurements with increasing number of frequencies. It emerged that when the inner scale exceeds the Fresnel radius, the transition from dual-frequency to triple-frequency measurements reduces the ionospheric correction. However, if the Fresnel radius exceeds the inner scale, such a transition decreases only the average correction. The correction variance in this case only rises. Let us, therefore, consider possibilities of improving the accuracy of GNSS measurements with the aid of spatial processing.
4. Eliminating GNSS Measurement Errors via Fresnel Inversion
As is well known, with the Green function at a given wave field on the plane where the wave leaves the inhomogeneous medium, it is possible to find the field on another, distant plane. When the medium between these planes is homogeneous, we can take the field of point source in free space as the Green function. This Kirchhoff approximation is widely used for analyzing wave propagation behind the slightly inhomogeneous layer. Yet the wave field where the wave leaves this layer can be described using the geometrical optics approximation. The Kirchhoff formula allows us to calculate this field in the distant region exhibiting various diffraction effects, caustics, and strong fluctuations for which the geometrical optics approximation does not work. Since the small-angle approximation is applicable to large-scale irregularities, the Kirchhoff formula takes the form of Fresnel transform :
When the inhomogeneous layer is sufficiently thin, it modulates only the transmitted wave phase and from (44) we can obtain the phase screen model which has been used for analyzing various problems of wave propagation in a randomly inhomogeneous ionosphere.
On the other hand, inversion of Fresnel transform (44) (Fresnel inversion) allows us, by measuring the field on the plane , to find the field on the plane nearby the irregularity where diffraction effects are absent and thus increases the accuracy of dual-frequency combinations (4) and (5) in determining and . Besides reduction of diffraction effects, ray path variations are low and consequently the third-order error is small nearby the irregularity. Let us analyze the action of the inverse operator in the form of inverse Fresnel transform :
Substituting second Rytov approximation (7)–(11) in (45) yields (for details see ) formulas analogous to the above ones. When the virtual screen plane is in a homogeneous medium between an irregularity and an observer, it is the natural result of using the Kirchhoff formula in determining the field on the plane from values of this field on the plane . However, formally we may apply transform (45) at any position of the virtual screen. This is also suggested by the results of using the double weighted Fourier transform to study wave propagation through a distant inhomogeneous region .
Figure 3 presents the results of calculation of the average correction (а) and the standard deviation (b) of first- and third-order corrections in determining the pseudo distance from ionosphere-free combination (4) after Fresnel processing (45) of GNSS measurements. These results are represented as virtual screen position functions at m for elevation angles of (blue line), (green line), and (red line). After spatial processing (45), the condition of validity of geometrical optics approximation (16) takes the formAs is seen from (46), bringing the virtual screen nearer to the irregularity () leads to the elimination of diffraction effects. Note that corresponds to the above case of absence of spatial processing. Figure 3 shows that the absolute average corrections (see Figure 3(a)) and standard deviations (see Figure 3(b)) decrease with increasing distance between the observer and the screen and have minimums at the defined . These minimums for a narrow region with irregularities are close to zero. However, in the case of large extent of the region with irregularities, condition (46) works only for some irregularities and error values in increase. According to the analysis, these minimums are located at distances corresponding roughly to the points of intersection of the observer-satellite line and the layer maximum which for our model of irregularities corresponds to the maximum. With this position of the virtual screen, the error is minimal because condition (46) is fulfilled for the most part of the inhomogeneous layer. From the well-known relation between wave amplitude fluctuations and diffraction effects  follows an analogous dependence of the amplitude scintillation index on the screen position after spatial processing (45).
Figure 4 for the L1 frequency ( Hz) with the same parameters presents calculations of the amplitude scintillation index of the wave processed by the algorithm of Fresnel inversion (45). This figure implies that processing (45) reduces the scintillation index and the minimum of amplitude fluctuations is very close to the minimum of GNSS measurement errors.
In Figure 5 are the results of the simulation for dual-frequency TEC measurements with geometry-free combination (5). Referring to Figure 5(a), in this case the Fresnel inversion can also reduce the average correction of TEC measurements through the proper choice of distance to the virtual screen . However, the standard deviation of the first-order correction (dashed lines in Figure 5(b)) varies slightly during spatial processing. This is attributed to the said (Figure 2(b)) relationship between the variance of the first-order correction and the variance of the first-approximation phase which weakly depends on diffraction effects . However, the standard deviation of the third-order correction (solid lines in Figure 5(b)) in optimum Fresnel processing (at ) reduces due to elimination of the diffraction effects and decrease in the ray bending.
In this study, ionospheric errors of GNSS remote sensing using both ionosphere-free and geometry-free combinations were analyzed. In the former case, we demonstrated the possibility of reducing residual ionospheric errors in ionosphere-free combination by Fresnel inversion. The inversion parameter, the distance to the virtual screen, may be selected from the minimum of amplitude fluctuations. This suggests the possibility of improving the accuracy of GNSS remote sensing in meteorology.
In investigations into ionospheric disturbances through geometry-free combinations, the Fresnel inversion eliminates only the third-order correction. To eliminate the random TEC component which, like the measured average TEC, is the first-order correction, we should use temporal filtering (averaging). However, measurements of this random TEC component allow us to examine the fine structure of ionospheric plasma via GNSS.
The Fresnel spatial processing requires constructing a phase GNSS array of sufficiently large sizes. Despite the difficulty in constructing this array, it seems necessary to develop the array, taking into account feasibility of using it in many fields of research on the troposphere, ionosphere, earthquakes, and volcanic activity.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
The author is very grateful to O. A. Kulish for her assistance in preparing the English version of the paper. This study was supported by the grant from the Russian Scientific Foundation (Project no. 14-37-00027).
- S. Kedar, G. A. Hajj, B. D. Wilson, and M. B. Heflin, “The effect of the second order GPS ionospheric correction on receiver positions,” Geophysical Research Letters, vol. 30, no. 16, 2003.
- E. R. Kursinski, G. A. Hajj, J. T. Schofield, R. P. Linfield, and K. R. Hardy, “Observing Earth's atmosphere with radio occultation measurements using the global positioning system,” Journal of Geophysical Research D: Atmospheres, vol. 102, no. 19, pp. 23429–23465, 1997.
- M. D. Moïtensen and P. Høeg, “Inversion of GPS occultation measurements using Fresnel diffraction theory,” Geophysical Research Letters, vol. 25, no. 13, pp. 2441–2444, 1998.
- R. E. Schaal and A. P. C. Larocca, “A methodology for monitoring vertical dynamic sub-centimeter displacements with GPS,” GPS Solutions, vol. 5, no. 3, pp. 15–18, 2002.
- H. C. Baker, A. H. Dodson, N. T. Penna, M. Higgins, and D. Offiler, “Ground-based GPS water vapour estimation: potential for meteorological forecasting,” Journal of Atmospheric and Solar-Terrestrial Physics, vol. 63, no. 12, pp. 1305–1314, 2001.
- V. Janssen, “GPS volcano deformation monitoring (GPS on the Web),” GPS Solutions, vol. 6, pp. 128–130, 2002.
- S. Jin, E. Cardellach, and F. Xie, GNSS Remote Sensing. Theory, Methods and Applications, Springer, Dordrecht, The Netherlands, 2014.
- V. Hoyle, S. Skone, and S. Shrestha, “The use of the global positioning system to derive atmospheric water vapor distribution for environmental applications,” Environmental Informatics Archives, vol. 1, pp. 512–523, 2003.
- M. S. Grewal, L. R. Weill, and A. P. Andrews, Global Positioning Systems, Inertial Navigation, and Integration, Wiley, New York, NY, USA, 2007.
- G. Xu, GPS Theory, Algorithms and Applications, Springer, New York, NY, USA, 2nd edition, 2007.
- J. B. Tsui, Fundamentals of Global Positioning System Receivers: A Software Approach, John Wiley & Sons, Hoboken, NJ, USA, 2nd edition, 2005.
- P. J. G. Teunissen, “Theory of integer equivariant estimation with application to GNSS,” Journal of Geodesy, vol. 77, no. 7-8, pp. 402–410, 2003.
- S. Bassiri and G. A. Hajj, “Higher-order ionospheric effects on the global positioning system observables and means of modeling them,” Manuscripta Geodaetica, vol. 18, pp. 280–289, 1993.
- M. Fritsche, R. Dietrich, C. Knöfel et al., “Impact of higher-order ionospheric terms on GPS estimates,” Geophysical Research Letters, vol. 32, no. 23, 2005.
- M. M. Hoque and N. Jakowski, “Estimate of higher order ionospheric errors in GNSS positioning,” Radio Science, vol. 43, no. 5, Article ID RS5008, 2008.
- G. Beutler, E. Brockmann, S. Fankhauser et al., Bernese GPS Software Version 4.0, University of Berne, Bern, Switzerland, 1996.
- J. A. Klobuchar and J. M. Kunches, “Comparative range delay and variability of the earth's troposphere and the ionosphere,” GPS Solutions, vol. 7, no. 1, pp. 55–58, 2003.
- M. Elsobeiey and A. El-Rabbany, “On modelling of second-order ionospheric delay for GPS precise point positioning,” The Journal of Navigation, vol. 65, no. 1, pp. 59–72, 2012.
- M. Hernández-Pajares, J. M. Juan, J. Sanz, and R. Orús, “Second-order ionospheric term in GPS: implementation and impact on geodetic estimates,” Journal of Geophysical Research: Solid Earth, vol. 112, Article ID B08417, 2007.
- H. Munekane, “A semi-analytical estimation of the effect of second-order ionospheric correction on the GPS positioning,” Geophysical Journal International, vol. 163, no. 1, pp. 10–17, 2005.
- K. Palamartchouk, “Apparent geocenter oscillations in Global Navigation Satellite Systems solutions caused by the ionospheric effect of second order,” Journal of Geophysical Research B: Solid Earth, vol. 115, no. 3, Article ID B03415, 2010.
- M. V. Tinin and E. V. Konetskaya, “Influence of geomagnetic field on the ionospheric error of satellite navigation systems,” Geomagnetism and Aeronomy, vol. 53, pp. 741–749, 2013.
- B. C. Kim and M. V. Tinin, “Potentialities of multifrequency ionospheric correction in Global Navigation Satellite Systems,” Journal of Geodesy, vol. 85, no. 3, pp. 159–169, 2011.
- B.-C. Kim and M. V. Tinin, “Effect of ionospheric irregularities on accuracy of dual-frequency GPS systems,” Geomagnetism and Aeronomy, vol. 47, no. 2, pp. 238–243, 2007.
- K. C. Yeh and C. H. Liu, “Radio-wave scintillations in the ionosphere,” Proceedings of the IEEE, vol. 70, no. 4, pp. 324–360, 1982.
- B. N. Gershman, L. M. Erukhimov, and Y. Y. Yashin, Wave Phenomena in the Ionosphere and Cosmic Plasma, Nauka, Moscow, Russia, 1984, (Russian).
- A. D. Wheelon, Electromagnetic Scintillation II. Weak Scattering, Cambridge University Press, Cambridge, UK, 2001.
- M. V. Tinin, “Eliminating diffraction effects during multi-frequency correction in global navigation satellite systems,” Journal of Geodesy, vol. 89, 2015.
- E. L. Afraimovich, N. P. Perevalova, and S. V. Voyeikov, “Traveling wave packets of total electron content disturbances as deduced from global GPS network data,” Journal of Atmospheric and Solar-Terrestrial Physics, vol. 65, no. 11–13, pp. 1245–1262, 2003.
- V. I. Tatarskii, The Effect of a Turbulent Atmosphere on Wave Propagation, National Technical Information Service, Springfield, Va USA, 1971.
- B. C. Kim and M. V. Tinin, “The second-order Rytov approximation and residual error in dual-frequency satellite navigation systems,” Waves in Random and Complex Media, vol. 19, no. 2, pp. 284–304, 2009.
- S. M. Rytov, Y. A. Kravtsov, and V. I. Tatarskii, Introduction to Statistical Radiophysics, Vol 4. Wave Propagation through Random Media, Springer, New York, NY, USA, 1989.
- M. V. Tinin, “Integral representation of the field of the wave propagating in a medium with large-scale irregularities,” Radiophysics and Quantum Electronics, vol. 55, no. 6, pp. 391–398, 2012.
Copyright © 2015 M. V. Tinin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.