Next Article in Journal
GWO-Based Multi-Stage Algorithm for PMDC Motor Parameter Estimation
Next Article in Special Issue
Analysis of 5G and LTE Signals for Opportunistic Navigation and Time Holdover
Previous Article in Journal
Adaptive Beam Divergence Control to Mitigate Scintillation Effect Caused by Pointing Error in Vertical FSO Transmissions
Previous Article in Special Issue
GNSS Data/Pilot Combining with Extended Integrations for Carrier Tracking
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Skywave Detection and Mitigation for the MF R-Mode Continuously Operating Reference Station

1
Korea Research Institute of Ships and Ocean Engineering, Daejeon 34103, Republic of Korea
2
Department of Ship and Ocean Engineering Major, University of Science and Technology, Daejeon 34113, Republic of Korea
3
School of Integrated Technology, Yonsei University, Incheon 21983, Republic of Korea
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(11), 5046; https://doi.org/10.3390/s23115046
Submission received: 21 March 2023 / Revised: 18 May 2023 / Accepted: 22 May 2023 / Published: 24 May 2023

Abstract

:
There is an increasing need for an independent terrestrial navigation system, owing to the increasing reliance on global navigation satellite systems (GNSS). The medium-frequency range (MF R-Mode) system is considered a promising alternative; however, the skywave effect caused by ionospheric changes at night can degrade its positioning accuracy. To address this problem, we developed an algorithm to detect and mitigate the skywave effect on MF R-Mode signals. The proposed algorithm was tested using data collected from Continuously Operating Reference Stations (CORS) monitoring the MF R-Mode signals. The skywave detection algorithm is based on the signal-to-noise ratio (SNR) induced by the groundwave and skywave composition, whereas the skywave mitigation algorithm was derived from the I and Q components of the signals obtained through IQ modulation. The results demonstrate a significant improvement in the precision and standard deviation of the range estimation using CW1 and CW2 signals. The standard deviations decreased from 39.01 and 39.28 m to 7.94 and 9.12 m, respectively, while the precision (2-sigma) increased from 92.12 and 79.82 m to 15.62 and 17.84 m, respectively. These findings confirm that the proposed algorithms can enhance the accuracy and reliability of MF R-Mode systems.

1. Introduction

The global navigation satellite system (GNSS) is a vital component of the modern world, enabling precise positioning and navigation for various applications such as transportation, mapping, and military operations [1,2,3]. GNSS is based on a constellation of satellites that orbit the Earth and transmit signals received by GNSS receivers on Earth. These signals are then used to calculate the position, velocity, and time of the receiver. GNSS can be used in navigation systems for aviation, shipping, and other modes of transportation, as well as in financial systems, communication networks, and energy grids [4]. The widespread reliance on GNSS has made it an important component of modern society, and its availability and accuracy are essential for ensuring the safety, efficiency, and security of critical systems and applications.
However, GNSSs are vulnerable to various types of interference, such as jamming and spoofing, which can compromise their accuracy and reliability, resulting in significant consequences for GNSS users [5,6]. For example, GNSS is critical for ensuring the safety and efficiency of vehicles in the transportation sector, particularly in autonomous systems such as self-driving cars. In the event of GNSS interference, these vehicles can become unguided, leading to potential accidents. Considering these security concerns, multiplexing the navigation system using a radio navigation system on the ground can ensure safe navigation for GNSS users [7].
Terrestrial radionavigation systems such as the enhanced long-range navigation (eLoran) system and the medium-frequency (MF) R-Mode system can provide a backup to the GNSS to mitigate its vulnerabilities [8,9,10]. These systems use ground-based radio signals to provide navigation information to users, offering a highly accurate and reliable alternative to GNSS [11]. Unlike GNSS signals, which are susceptible to jamming and interference, terrestrial navigation signals are much less vulnerable to such threats and can provide navigation information in areas where GNSS signals are weak or unavailable. Therefore, the deployment of terrestrial navigation signals can increase the resilience of navigation systems and ensure the safety and efficiency of navigation, particularly in critical applications such as aviation, shipping, and national defense [12].
Terrestrial radionavigation signals vary in propagation velocity depending on the characteristics of the propagation environment, causing changes in the time of flight of the navigation signals and becoming the main error source in position estimation [13,14]. This variation is influenced by regional and temporal factors such as terrain, weather, and season. To compensate for this, it is necessary to continuously monitor signals in the service area and generate correction information. The eLoran system can provide position estimation with an accuracy of approximately 10 to 20 m by applying real-time correction information [15]. This is similar to the concept of Continuously Operating Reference Stations (CORS) operated by the GNSS [16].
The skywave effect is a phenomenon that occurs when the D layer of the ionosphere disappears as the sun sets, causing eLoran or MF R-Mode signals to be refracted back to the surface of the Earth [17,18]. This effect can cause a significant degradation in the accuracy of terrestrial radionavigation systems, considering that the travel time of the navigation signals is affected by changes in the refractive index of the atmosphere [19]. The eLoran system is designed to be resistant to the skywave effect by using pulsed signals; however, the MF R-Mode system, which uses continuous wave signals, is susceptible to the skywave effect, thereby decreasing the performance of the system. Several MF R-Mode tests conducted in the Baltic Sea showed that the accuracy of a ship’s positioning is significantly reduced at night, owing to the impact of the skywave effect [20].
South Korea conducted a research project to develop MF R-Mode technology that could be used in case of GNSS failure by implementing four MF R-Mode transmission stations at Socheong-do (SC), Palmi-do (PM), Chungju (CJ), and Eocheong-do (EC), and one CORS to achieve a 10 m position accuracy in the nearby seas. The MF R-Mode CORS monitors changed across a broad region and generated correction information. However, the skywave effect can be difficult to predict as it occurs suddenly and irregularly in different areas of the ionosphere, making it difficult to consistently apply to users over a large area. Therefore, the impact of the skywave effect must be excluded to generate accurate correction information for the MF R-Mode CORS.
This paper presents an algorithm to eliminate the skywave effect that occurs at night, and verifies its performance. First, we describe the construction of the CORS for the MF R-Mode system in detail and analyze the phase change of the MF R-Mode signal according to the time variation at the CORS. Subsequently, a signal-to-noise ratio (SNR)-based skywave detection method and an IQ modulation-based skywave reduction algorithm are proposed. Finally, the verification of the effectiveness of the algorithm through actual experiments near the CORS is presented in detail.

2. Continuously Operating Reference Station of the Korean MF R-Mode System

The area near Daesan Port, which serves as the testbed in this study for validating the performance of the MF R-Mode, has various propagation paths with different distances and terrains to each of the four MF R-Mode transmission stations: SC, PM, CJ, and EC. In general, the range estimation accuracy for MF R-Mode signals tends to decrease with an increase in the propagation path length, owing to signal strength attenuation. However, this attenuation is relatively less pronounced when the signal travels over sea, compared to on land. Specifically, for an MF R-Mode signal with the same power amplitude, range estimation can be performed with higher precision when the distance to the user receiver is shorter and the propagation path has a lower percentage of land.
The propagation paths of signals from the SC and PM transmission stations were the longest and shortest, with approximate lengths of 171 and 39 km, respectively, mostly over the sea. The signal from the CJ transmission station, which had a propagation path length of approximately 119 km, was mostly over land and propagated through mountainous areas near the transmission station. The signal from the EC transmission station had the characteristic of passing through both land and sea, to some extent. The map on the left-hand side of Figure 1 shows the geometric arrangement of the four MF R-Mode transmitters and the Daesan CORS and the distances of the signals. The waveforms on the right-hand side show the elevation profiles of the signal propagation paths from each transmitter to the MF R-Mode CORS, which were produced from Google Earth.

2.1. Horizontal Dilution of Precision

The position dilution of precision (PDOP) is a term used in satellite navigation systems, such as the Global Positioning System (GPS), to describe the geometric quality of the satellite configuration as seen from a receiver [21]. The PDOP values provide information about the relative accuracy of the GPS position solution given the current satellite geometry. The smaller the PDOP value, the better the satellite geometry and the more accurate the positioning solution. Furthermore, the concept of PDOP can be applied to terrestrial-based navigation systems. A low PDOP between the transmitter and user implies that a higher accuracy in position estimation can be achieved when receiving the same signal quality. Consequently, analyzing the positioning accuracy performance based on the PDOP of CORS can help to predict the positioning accuracy performance that users can use in other marine regions.
Generally, the PDOP includes the distribution of satellites and receivers along the x-, y-, and z-axes based on the Earth-centered, Earth-fixed (ECEF) coordinate system. However, the propagation system based on terrestrial signals is different from the satellite navigation system, and only calculates the horizontal position based on the World Geodetic System 1984 (WGS84) coordinate system; therefore, the horizontal dilution of precision (HDOP) was used. The HDOP can be derived as follows:
A = [ sin θ 1 cos θ 1 1 sin θ 2 cos θ 2 1 sin θ n cos θ n 1 ] ,
Q = ( A T A ) 1 ,
H D O P = Q ( 1 , 1 ) 2 + Q ( 2 , 2 ) 2 ,
where θ i indicates the azimuth angle between the CORS position and the i-th MF R-Mode transmitter.

2.2. Equipment Installation for Data Collection

In the MF R-Mode CORS testbed, facilities were assembled to continuously receive and monitor MF R-Mode signals to generate correction information that compensates for the phase change of MF R-Mode signals changing over time. To achieve this, MF R-Mode and GNSS antennas were installed on the roofs of the CORS facilities, as shown in Figure 2, and Serco’s MF R-Mode receiver (Model: MFR-1a) was used to install the signal-collection and processing device. The receiver from Serco is the only commercially available MF R-Mode receiver whose performance has been verified through a Complementary Positioning, Navigation, and Timing (PNT) performance validation conducted by the U.S. Department of Transportation [22].
The MFR-1a specification features a marine beacon band front end, designed to enhance the signal by removing noise and interference from adjacent low-frequency (LF) and amplitude-modulated (AM) bands. This front end has a gain range of 0–20 dB and amplifies the signal in the frequency range of 285–325 kHz. The signal is then processed using an analog-to-digital converter (ADC) that records the entire marine beacon band as a single signal for synchronous processing. The ADC operates with a 10 MHz clock signal from an external rubidium clock and has a sample rate of 1 MHz, with 16-bit resolution for a wide dynamic range. Additionally, the MFR-1a specification includes a separate GNSS receiver, which provides a position for initial receiver calibration and ensures the long-term stability of the rubidium clock.

2.3. Diurnal Phase Variation at the MF R-Mode CORS

The MF R-Mode signals are produced by updating the differential global navigation satellite system (DGNSS) infrastructure, which transmits GNSS correction information. The signals had a bandwidth of 100 or 200 bps and were transmitted using minimum shift key (MSK) modulation centered at a frequency of approximately 300 kHz. To generate the MF R-Mode signal, a continuous wave (CW) was inserted into the null point in the spectrum of the MSK signal. In Europe, where the signal bandwidth is 100 bps, a null was created 225 Hz away from the center frequency. Meanwhile, in South Korea, where the signal bandwidth is 200 bps, a CW signal was added 250 Hz away from the center frequency.
The phenomenon of diurnal phase changes in LF or MF signals refers to the changes in the phases of these signals over the day owing to the effective conductivity of the surface and humidity of the atmosphere. In eLoran systems using LF signals, the time of arrival (TOA) changes over time; therefore, real-time correction information is generated to compensate for the changes in TOA based on the differential eLoran stations [9].
Figure 3 illustrates the results of converting the phases of the signals of the four MF R-Mode transmitters into units of meters over a 24 h period on 8 January 2023. The phases of the MF R-Mode signals ranged between 0 and 2 pi, and could be converted into a range within the wavelength range corresponding to the frequency of the signal, from 0 m. The blue dots represent the phase of the MF R-Mode signal output by the receiver, while the red dots correspond to the result of applying a 100-s moving average function to remove noise and observe changes over time.
The phases of the four MF R-Mode signals exhibit 30 m changes over the course of the day, and real-time correction information must be generated from CORS to compensate for the daily changes to achieve the project goal of a 10 m positioning accuracy with a 95% probability. Some of the transmitters located on the island face challenges in transmitting stable MF R-Mode signals owing to power issues. Consequently, the phases of the received signals had random values at certain times, as shown in the first and third figures of Figure 3. Therefore, CORS should alert users to avoid using these signals for position estimation.
Furthermore, one of the most significant observations from the 24 h MF R-Mode signal data was that the phase underwent large fluctuations between 09:00 and 23:00 UTC. Considering that sunrise and sunset in Korea occurred at 23:32 and 08:43 UTC, respectively, on that day, this was likely because of skywave interference. In addition, the intensity of skywave interference may differ depending on factors such as the distance to each transmission station and the environment of signal propagation; however, it seems necessary to reduce skywave influence to attain a positioning accuracy of 10 m.

2.4. Skywave Effect at the MF R-Mode CORS

The skywave effect of MF signals refers to the reflection of MF radio waves from the ionosphere back into the Earth. As the sun sets and the Earth moves into darkness, the D layer of the ionosphere disappears, increasing the height of the ionosphere. This indicates that MF signals can no longer be absorbed by the D layer and, instead, are reflected back to Earth by the higher E layer of the ionosphere [23]. This allows MF signals to travel over the horizon and reach greater distances, making the skywave effect more effective for nighttime MF signals in terms of long-distance communication.
In the MF R-Mode system, which calculates the distance through variations in the phases of signals, the presence of the skywave effect can create a new composite signal along with the ground wave that travels along the surface path. Consequently, the receiver receives erroneous measurements in terms of signal strength and phase, which may take the form of fading signals or multipath. Ultimately, the skywave effect has a detrimental impact on the positioning accuracy from a navigation system standpoint, rather than from a communication perspective.
Even if the D layer of the ionosphere vanishes at night, reflected signals cannot be received in all regions owing to the skywave effect. It is possible to receive reflected signals from the ionosphere only from distances greater than the skip distance. The skip distance refers to the maximum distance that a radio wave can travel along the Earth’s surface before being reflected back to Earth by the ionosphere. In contrast, the critical angle is the maximum angle at which a radio wave must be transmitted for it to be reflected back to Earth by the ionosphere. As shown in Figure 4 and Equation (4), the skip distance ( d s ) is determined by the critical angle ( θ c ) and the height of the E layer of the ionosphere ( h ).
d s = 2 h tan θ c
Previous studies on the impact of skywaves on MF signals show that the skip distance is 45 km from the transmitter [18], indicating that the positioning accuracy may be impacted at night in areas further than 45 km from the MF R-Mode signal transmitter. Recent experimental results indicate the absence of skywave effects at locations approximately 60 km from the transmitter [19]. However, there has also been evidence that skywave effects can be observed within an approximate range of 50 km, depending on the presence of abrupt mountainous regions along the transmission path [24].
A prior study [18] suggested that we can compare the signal attenuations experienced when MF signals travel 60 km and 500 km through paths consisting solely of land surface and solely of ocean surface, respectively, and verify the attenuations to be of the same degree. The composition of the propagation path, and not just the distance from the transmitter, can result in varying levels of amplitude attenuation. Consequently, it is more appropriate to assess the presence of skywave effects by considering the change in SNR and the variability of the groundwave reduction based on the features of the propagation path, rather than simply relying on the distance from the transmitter.

3. Skywave Detection and Mitigation Algorithm

3.1. SNR-Based Skywave Detection Algorithm

Assuming the received signal at the receiver ( S ) is composed solely of groundwaves and skywaves, without considering various other noises, it can be expressed as follows:
S ( t ) = A g ( t ) cos ( 2 π f c t + θ g ( t ) ) + A s ( t ) cos ( 2 π f c t + θ s ( t ) ) ,
where A g ( t ) denotes the amplitude of the groundwave, f c denotes the center frequency of the MF R-Mode signal, θ g ( t ) represents the phase according to the propagation path of the groundwave, A s ( t ) indicates the amplitude of the skywave, and θ s ( t ) represents the phase according to the propagation path of the one-hop skywave.
During the daytime, most skywaves are absorbed by the D layer of the ionosphere; therefore, A s ( t ) can be assumed to converge to 0. As a result, A g ( t ) is nearly the same as the amplitude of S ( t ) , and hence it would represent a stable value as long as the environment along the propagation path does not rapidly change. However, during nighttime, when the skywave reaches the surface, the ratio of A g ( t ) to A s ( t ) becomes significant and cannot be ignored.
Serco’s MF R-Mode receiver, which is installed at the MF R-Mode CORS, computes the SNR by comparing the magnitude of the center frequency of the MF R-Mode signal to that of the surrounding frequencies, with no signal in the fast Fourier transform (FFT) spectrum of the received data over a 1-s interval. This enables fluctuations in the signal strength to be reflected in the SNR of the MF R-Mode signal, providing a means to detect the presence of the skywave effect. Such multipath detection technology utilizing the detection of SNR changes is already used in fields such as GNSS-R [25].
Figure 5 shows the SNRs received from each transmitter, i.e., EC, CJ, PM, and SC, over 24 h. The data were collected over the same durations as those shown in Figure 3. The blue and red dots indicate the SNR output from the receiver and the results after passing through a moving average filter for 100 s, respectively.
The proposed skywave detection algorithm entails the following steps. The output SNR from the receiver is processed through a 100 s moving average filter (represented by the red points in Figure 5). If the standard deviation of these results, calculated over a 100 s time period, exceeds the standard deviation mean value of the 100 s moving average SNR over the course of the day, the skywave effect takes place as per Equation (6). The threshold for this detection and the length of the moving average filter were established experimentally using long-term data collected at the CORS.
s k y w a v e   f l a g = 1 i f   σ S N R ¯ ( t ) > σ S N R ¯ ¯
For example, the proposed algorithm was applied to the MF R-Mode signal from the CJ transmitter, which was expected to be the most affected by skywaves among the signals received at CORS, considering the distance and topographical features of the propagation path from the transmitter. In Figure 6, the statistical values for detecting the skywave effect are represented by blue dots, and the threshold is represented by a red line. The time when the statistical values exceeded the threshold was almost similar to the time when the phase changed abruptly, as shown in Figure 4.

3.2. Skywave Mitigation Algorithm Using IQ Modulation

In IQ modulation, the amplitude and phase of a carrier signal are modulated using two separate signals, referred to as in-phase (I) and quadrature (Q) signals. The I and Q signals are typically generated by mixing the carrier signal with two separate local oscillator signals, one of which is shifted in-phase by 90°. The I and Q signals can be expressed mathematically as
I ( t ) = S ( t ) cos ( 2 π f c t ) Q ( t ) = S ( t ) sin ( 2 π f c t ) .
Here, f c is the carrier frequency, and I ( t ) and Q ( t ) are the in-phase and quadrature signals, respectively. The in-phase signal controls the amplitude of the carrier signal, whereas the quadrature signal controls the phase of the carrier signal.
The received signal can be separated into its I and Q components by using IQ modulation, enabling us to analyze and manipulate each component independently. This is particularly useful in situations where there are multiple paths for the signal to reach the receiver, resulting in multipath interference and signal fading. By analyzing the I and Q components independently, we can detect and correct for the effects of multipath interference, leading to a more accurate and reliable signal.
The demodulation process using the phase angle of the complex baseband signal in IQ demodulation involves computing the arctangent of the ratio of quadrature and in-phase signals. This results in a phase-angle value that represents the phase of the original carrier signal. Typically, a complex baseband signal is represented as a complex number, with in-phase and quadrature signals forming the real and imaginary components of the complex number, respectively. The phase angle ( θ ) of the complex baseband signal can be computed as
θ = arctan ( Q ( t ) I ( t ) ) .
MFR-1a, a Serco MF R-Mode receiver, estimates the range of the signal by utilizing the phase angle calculated from the I and Q components of the complex baseband of the MF R-Mode signal, which is subsequently converted from radians to meters. MFR-1a can simultaneously receive up to five signals from MF R-Mode transmitters and output the calculated I and Q values each second. Figure 7 shows the phase and I/Q of the MF R-Mode signal from the CJ transmitter within an hour, with significant phase changes owing to the skywave effect within 24 h. The MF R-Mode signal comprises two CWs that are shifted by the same frequency difference from the center frequency in the positive and negative directions. In this study, signals with frequencies higher than the center frequency were defined as CW1, whereas those with frequencies lower than the center frequency were defined as CW2. The first row of Figure 7 shows the CW1 signal, and the second row shows the CW2 signal.
When observing the phase graphs in column 1 of Figure 7, it can be observed that both CW1 and CW2 experience significant phase changes of several hundred meters for 1 h owing to the skywave effect. Furthermore, the trends in these changes exhibit opposite directionality. However, this correlation was found to be more pronounced when the IQ graph in column 2 of Figure 7 was examined. The red and blue lines in the IQ graph represent the I and Q components, respectively. Furthermore, the I component of CW1 (the red line in the upper-right graph) exhibits a completely opposite trend to that of the Q component of CW2 (the blue line in the lower-right graph). In addition, the Q component of CW1 (blue line in the upper-right graph) exhibits a trend similar to that of the I component of CW2 (red line in the lower-right graph).
Based on long-term data collected at CORS, the aforementioned characteristics have been consistently observed without loss of generality. Based on these observations, we established two assumptions regarding the characteristics that occur when the skywave effect is present.
Assumption 1.
The variation in the I component of CW1 is negatively correlated with the variation in the Q component of CW2.
Assumption 2.
The variation in the Q component of CW1 is positively correlated with the variation in the I component of CW2.
An algorithm to implement these two assumptions can be expressed as follows:
θ ^ CW 1 = arctan 1 ( Q CW 1 Δ Q CW 1 + Δ I CW 2 2 I CW 1 Δ I CW 1 Δ Q CW 2 2 ) ,
θ ^ CW 2 = arctan 1 ( Q CW 2 Δ Q CW 2 Δ I CW 1 2 I CW 2 Δ I CW 2 + Δ Q CW 1 2 ) ,
where θ ^ CW 1 and θ ^ CW 2 denote the phases of the two MF R-Mode signals with mitigated skywave effects, respectively.

4. Results

The proposed algorithm was applied to the MF R-Mode signal from the CJ transmitter, which is most affected by skywaves at the CORS in Daesan Port, to verify its performance. Figure 8 shows the phases of the CW1 and CW2 signals from the MRF-1a receiver (black line), and the phases estimated again after applying the proposed algorithm (red line) provided by the receiver.
The measured phase values for the CW1 and CW2 signals exhibited standard deviations of 39.01 and 39.28 m, respectively, with a 95% level of precision, yielding values of 92.12 and 79.82 m, respectively. However, on applying the proposed algorithm, the standard deviations of the CW1 and CW2 signals were reduced to 7.94 and 9.12 m, respectively, whereas the 95% levels of precision were improved to 15.62 and 17.84 m, respectively. The results are presented in Table 1.
To confirm whether the proposed algorithm can be applied regardless of the receiver location, additional experiments were conducted at two points near the CORS. An antenna park was set up to install GNSS and MF R-Mode antennas on the vehicle loop, and the same signal-collection setup as in the CORS was implemented using an MFR-1a receiver. As shown in Figure 9, the two additional experimental locations were approximately 11 and 24 km away from the CORS and were selected at clear coastal locations to minimize the influence of the terrain.
The distances between the new measurement points (#1, #2) and the CJ transmitter were approximately 10 and 20 km longer, respectively, than the distance between the CORS and CJ transmitter; however, precision degradation did not increase proportionally with the distance, as shown in Table 2, owing to the complex terrain in Korea. Based on the additional validation at these two locations, we can confirm that our proposed algorithm can significantly reduce the impact of skywaves and improve precision. Notably, the effect of error reduction caused by the skywave in the second additional site is remarkable, with a reduction of 10 to 20 times.

5. Conclusions

We developed an algorithm to detect and mitigate the skywave effect of MF R-Mode signals caused by ionospheric changes at night, and evaluated its performance. The MF R-Mode system, a terrestrial navigation system that functions as an independent navigation system when GNSS is used in difficult circumstances in the ocean, uses signals in the medium-frequency band. However, the signal is negatively affected by the skywave reflected from the E layer when the ionospheric D layer disappears at night, which significantly decreases the system’s accuracy.
We installed CORS as a testbed to detect and minimize the skywave effect, and then collected signals from 2021 for our investigation. Because a skywave is a type of multipath signal, it combines with a groundwave that propagates through the ground, causing a rapid change in SNR. Therefore, we proposed a skywave detection algorithm using the SNR. Furthermore, because ionospheric changes are very irregular and unpredictable, it is difficult to secure the ground truth regarding the occurrence of skywaves. However, we confirmed that we could detect the moments when the phase changed significantly using the change in the SNR.
Additionally, by observing the data for a long period of time, we confirmed that we can effectively remove skywaves by using the I and Q signals generated in IQ modulation, and verified the performance of the algorithm using the data from the hour when the skywave effect was most severe. Furthermore, we confirmed that the proposed algorithm can be applied regardless of the receiver’s location, through experiments conducted at two additional locations along with the data from CORS. The algorithm proposed in this study can utilize signals that otherwise could not be used for position determination owing to the skywave effect, and this is expected to have a positive impact on the availability and accuracy performance of the MF R-Mode system.

Author Contributions

Conceptualization, P.-W.S. and T.H.F.; methodology, P.-W.S., J.P., J.Y. and S.J.; software, P.-W.S. and J.P.; validation, P.-W.S. and T.H.F.; formal analysis, P.-W.S. and T.H.F.; investigation, P.-W.S., J.P., J.Y. and S.J.; writing—original draft preparation, P.-W.S.; writing—review and editing, T.H.F.; visualization, P.-W.S.; supervision, T.H.F.; project administration, Y.H. and T.H.F.; funding acquisition, Y.H. and T.H.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was part of a project titled “Development of integrated R-Mode navigation system”, funded by the Ministry of Oceans and Fisheries, Korea (Grant No. 1525012261).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bakdi, A.; Glad, I.K.; Vanem, E. Testbed scenario design exploiting traffic big data for autonomous ship trials under multiple conflicts with collision/grounding risks and spatio-temporal dependencies. IEEE Trans. Intell. Transp. Syst. 2021, 22, 7914–7930. [Google Scholar] [CrossRef]
  2. Dardanelli, G.; Maltese, A. On the accuracy of cadastral marks: Statistical analyses to assess the congruence among GNSS-based positioning and official maps. Remote Sens. 2022, 14, 4086. [Google Scholar] [CrossRef]
  3. Hebeler, E.; McKneely, J.; Rigsbee, S. The application of human-systems integration: Designing the next generation of military global positioning system handheld devices. Johns Hopkins APL Tech. Dig. 2012, 31, 66–75. [Google Scholar]
  4. Manandhar, D. GNSS for Policy and Decision Makers GNSS Applications. Available online: https://www.unoosa.org/documents/pdf/icg/2021/Tokyo2021/ICG_CSISTokyo_2021_16.pdf (accessed on 24 February 2023).
  5. Grant, A.; Williams, P.; Ward, N.; Basker, S. GPS jamming and the impact on maritime navigation. J. Navig. 2009, 62, 173–187. [Google Scholar] [CrossRef]
  6. Seo, J.; Chen, Y.-H.; De Lorenzo, D.S.; Lo, S.; Enge, P.; Akos, D.; Lee, J. A real-time capable software-defined receiver using GPU for adaptive anti-jam GPS sensors. Sensors 2011, 11, 8966–8991. [Google Scholar] [CrossRef] [PubMed]
  7. Kim, E.; Seo, J. SFOL Pulse: A high accuracy DME pulse for alternative aircraft position and navigation. Sensors 2017, 17, 2183. [Google Scholar] [CrossRef] [PubMed]
  8. Son, P.W.; Rhee, J.H.; Seo, J. Novel multichain-based Loran positioning algorithm for resilient navigation. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 666–679. [Google Scholar] [CrossRef]
  9. Son, P.-W.; Park, S.G.; Han, Y.; Seo, K. eLoran: Resilient positioning, navigation, and timing infrastructure in maritime areas. IEEE Access 2020, 8, 193708–193716. [Google Scholar] [CrossRef]
  10. Gewies, S.; Dammann, A.; Ziebold, R.; Bäckstedt, J.; Bronk, K.; Wereszko, B.; Rieck, C.; Gustafson, P.; Eliassen, C.; Hoppe, M.; et al. R-Mode testbed in the Baltic Sea. In Proceedings of the 19th IALA Conference 2018, Incheon, Republic of Korea, 27 May–2 June 2018. [Google Scholar]
  11. Son, P.-W.; Rhee, J.H.; Hwang, J.; Seo, J. Universal kriging for loran ASF map generation. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 1828–1842. [Google Scholar] [CrossRef]
  12. Gewies, S.; Grundhöfer, L.; Hehenkamp, N. Availability of maritime radio beacon signals for R-Mode in the southern Baltic Sea. TransNav Int. J. Mar. Navig. Saf. Sea Transp. 2020, 14, 173–178. [Google Scholar] [CrossRef]
  13. Rhee, J.H.; Kim, S.; Son, P.-W.; Seo, J. Enhanced accuracy simulator for a future Korean nationwide eLoran system. IEEE Access 2021, 9, 115042–115052. [Google Scholar] [CrossRef]
  14. Kim, W.; Son, P.-W.; Park, S.G.; Park, S.H.; Seo, J. First demonstration of the Korean eLoran accuracy in a narrow waterway using improved ASF maps. IEEE Trans. Aerosp. Electron. Syst. 2020, 58, 1492–1496. [Google Scholar] [CrossRef]
  15. Offermans, G.; Johannessen, E.; Bartlett, S.; Schue, C.; Grebnev, A.; UrsaNav, Inc.; Bransby, M.; Williams, P.; Hargreaves, C. Loran initial operational capability in the United Kingdom-first results. In Proceedings of the International Technical Meeting of the Institute of Navigation (ION), Dana Point, CA, USA, 26–28 January 2015; pp. 27–39. [Google Scholar]
  16. Miya, M.; Fujita, S.; Sato, Y.; Ota, K.; Hirokawa, R.; Takiguchi, J. Centimeter Level Augmentation Service (CLAS) in Japanese Quasi-Zenith satellite system, its user interface, detailed design, and plan. In Proceedings of the 29th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2016), Portland, OR, USA, 12–16 September 2016; pp. 2864–2869. [Google Scholar]
  17. ITU. The Ionosphere and Its Effects on Radiowave Propagation. Available online: https://www.itu.int/dms_pub/itu-r/opb/hdb/R-HDB-32-1998-PDF-E.pdf (accessed on 24 February 2023).
  18. Poppe, D.C. Coverage and Performance Prediction of DGPS Systems Employing Radiobeacon Transmissions. Ph.D. Thesis, Bangor University, Bangor, UK, 1995. [Google Scholar]
  19. Rizzi, F.G.; Grundhöfer, L.; Gewies, S.; Ehlers, T. Performance assessment of the medium frequency R-Mode Baltic testbed at sea near Rostock. Appl. Sci. 2023, 13, 1872. [Google Scholar] [CrossRef]
  20. Jeong, S.; Son, P.-W. Preliminary analysis of skywave effects on MF DGNSS R-Mode signals during daytime and nighttime. In Proceedings of the 2022 IEEE International Conference on Consumer Electronics-Asia (ICCE-Asia), Yeosu, Republic of Korea, 26–28 October 2022; pp. 1–4. [Google Scholar] [CrossRef]
  21. Phillips, A.H. Geometrical determination of PDOP. Navigation 1984, 31, 329–337. [Google Scholar] [CrossRef]
  22. U.S. Department of Transportation. Complementary PNT and GPS Backup Technologies Demonstration Report. Available online: https://www.transportation.gov/administrations/assistant-secretary-research-and-technology/complementary-pnt-and-gps-backup (accessed on 24 February 2023).
  23. Higginson-Rollins, M.A.; Cohen, M.B. Measuring the electron density roughness of the D-region ionosphere. J. Geophys. Res. Space Phys. 2020, 125, 10. [Google Scholar] [CrossRef]
  24. Higginson-Rollins, M.A.; Cohen, M.B. Exploiting LF/MF signals of opportunity for lower ionospheric remote sensing. Geophys. Res. Lett. 2017, 44, 16. [Google Scholar] [CrossRef]
  25. Benton, C.J.; Mitchell, C.N. Isolating the multipath component in GNSS signal-to-noise data and locating reflecting objects. Radio Sci. 2011, 46, 1–11. [Google Scholar] [CrossRef]
Figure 1. Geometric distribution of the four MF R-Mode transmitters and the lengths of the propagation paths (left), and elevation profiles of each propagation path acquired from Google Earth (right). The yellow box in each profile presents the highest altitude along the propagation path.
Figure 1. Geometric distribution of the four MF R-Mode transmitters and the lengths of the propagation paths (left), and elevation profiles of each propagation path acquired from Google Earth (right). The yellow box in each profile presents the highest altitude along the propagation path.
Sensors 23 05046 g001
Figure 2. Installation of equipment at CORS to collect and monitor MF R-Mode signals. MF R-Mode E-field passive antenna and survey-grade GNSS antenna on the rooftop (left), Serco’s MF R-Mode receiver (MFR-1a), and display devices in the shelter (right).
Figure 2. Installation of equipment at CORS to collect and monitor MF R-Mode signals. MF R-Mode E-field passive antenna and survey-grade GNSS antenna on the rooftop (left), Serco’s MF R-Mode receiver (MFR-1a), and display devices in the shelter (right).
Sensors 23 05046 g002
Figure 3. Phase variation of MF R-Mode signals in meters during a 24 h campaign at CORS. The blue dots indicate the raw measurements from the receiver, and red dots indicate smoothed signals using a 100 s moving average filter.
Figure 3. Phase variation of MF R-Mode signals in meters during a 24 h campaign at CORS. The blue dots indicate the raw measurements from the receiver, and red dots indicate smoothed signals using a 100 s moving average filter.
Sensors 23 05046 g003
Figure 4. Illustration of the skywave propagation phenomenon and its key parameters, including critical angle, skip distance, and ionosphere.
Figure 4. Illustration of the skywave propagation phenomenon and its key parameters, including critical angle, skip distance, and ionosphere.
Sensors 23 05046 g004
Figure 5. SNR variation of MF R-Mode signals in meters received from EC, CJ, PM, and SC transmitters in sequence during the 24 h campaign at CORS. The blue dots indicate the raw measurements from the receiver, and red dots indicate smoothed signals using a 100 s moving average filter.
Figure 5. SNR variation of MF R-Mode signals in meters received from EC, CJ, PM, and SC transmitters in sequence during the 24 h campaign at CORS. The blue dots indicate the raw measurements from the receiver, and red dots indicate smoothed signals using a 100 s moving average filter.
Sensors 23 05046 g005
Figure 6. Variation in standard deviation of MF R-Mode signal’s SNR from CJ transmitter (denoted in blue) and the threshold of determining the presence of skywave (denoted in red).
Figure 6. Variation in standard deviation of MF R-Mode signal’s SNR from CJ transmitter (denoted in blue) and the threshold of determining the presence of skywave (denoted in red).
Sensors 23 05046 g006
Figure 7. Analysis of the phase and I/Q signals of the MF R-Mode signal from the CJ transmitter during a 1 h campaign at CORS: CW1 (upper) and CW2 (bottom) comparison.
Figure 7. Analysis of the phase and I/Q signals of the MF R-Mode signal from the CJ transmitter during a 1 h campaign at CORS: CW1 (upper) and CW2 (bottom) comparison.
Sensors 23 05046 g007
Figure 8. Improvement in the mitigating skywave effect at CORS by applying the proposed algorithm. The black and red dots show the raw measurements from the receiver and the result of the proposed algorithm, respectively.
Figure 8. Improvement in the mitigating skywave effect at CORS by applying the proposed algorithm. The black and red dots show the raw measurements from the receiver and the result of the proposed algorithm, respectively.
Sensors 23 05046 g008
Figure 9. Geometric distribution of CORS and two additional experiment sites (left), and potable experimental setting on the vehicle (right).
Figure 9. Geometric distribution of CORS and two additional experiment sites (left), and potable experimental setting on the vehicle (right).
Sensors 23 05046 g009
Table 1. Statistics about the improvement in the mitigation of the skywave effect at CORS.
Table 1. Statistics about the improvement in the mitigation of the skywave effect at CORS.
Receiver OutputProposed Algorithm Applied
CW1CW2CW1CW2
Standard Deviation39.01 m39.28 m7.94 m9.12 m
Precision (95%)92.12 m79.82 m15.62 m17.84 m
Table 2. Statistics about the improvement in the mitigation of the skywave effect at additional two sites.
Table 2. Statistics about the improvement in the mitigation of the skywave effect at additional two sites.
Standard DeviationReceiver OutputProposed Algorithm Applied
CW1CW2CW1CW2
#125.63 m20.78 m9.57 m13.78 m
#2130.49 m61.74 m4.97 m6.93 m
Precision (95%)Receiver OutputProposed Algorithm Applied
CW1CW2CW1CW2
#146.73 m41.36 m18.17 m25.43 m
#2253.99 m125.39 m9.82 m14.56 m
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Son, P.-W.; Park, J.; Yu, J.; Jeong, S.; Han, Y.; Fang, T.H. Skywave Detection and Mitigation for the MF R-Mode Continuously Operating Reference Station. Sensors 2023, 23, 5046. https://doi.org/10.3390/s23115046

AMA Style

Son P-W, Park J, Yu J, Jeong S, Han Y, Fang TH. Skywave Detection and Mitigation for the MF R-Mode Continuously Operating Reference Station. Sensors. 2023; 23(11):5046. https://doi.org/10.3390/s23115046

Chicago/Turabian Style

Son, Pyo-Woong, Jongmin Park, Jaewon Yu, Suhui Jeong, Younghoon Han, and Tae Hyun Fang. 2023. "Skywave Detection and Mitigation for the MF R-Mode Continuously Operating Reference Station" Sensors 23, no. 11: 5046. https://doi.org/10.3390/s23115046

APA Style

Son, P. -W., Park, J., Yu, J., Jeong, S., Han, Y., & Fang, T. H. (2023). Skywave Detection and Mitigation for the MF R-Mode Continuously Operating Reference Station. Sensors, 23(11), 5046. https://doi.org/10.3390/s23115046

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop