Articles | Volume 9, issue 3
https://doi.org/10.5194/amt-9-1303-2016
© Author(s) 2016. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
https://doi.org/10.5194/amt-9-1303-2016
© Author(s) 2016. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
Statistical framework for estimating GNSS bias
Haystack Observatory, Massachusetts Institute of Technology, Route 40 Westford, 01469 MA, USA
Anthea J. Coster
Haystack Observatory, Massachusetts Institute of Technology, Route 40 Westford, 01469 MA, USA
William C. Rideout
Haystack Observatory, Massachusetts Institute of Technology, Route 40 Westford, 01469 MA, USA
Philip J. Erickson
Haystack Observatory, Massachusetts Institute of Technology, Route 40 Westford, 01469 MA, USA
Johannes Norberg
Finnish Meteorological Institute, P.O. Box 503, 00101 Helsinki, Finland
Related authors
Kian Sartipzadeh, Andreas Kvammen, Björn Gustavsson, Njål Gulbrandsen, Magnar Gullikstad Johnsen, Devin Huyghebaert, and Juha Vierinen
EGUsphere, https://doi.org/10.5194/egusphere-2025-3070, https://doi.org/10.5194/egusphere-2025-3070, 2025
Short summary
Short summary
Knowing charged particle densities high above Earth is key for forecasting space weather effects on satellites and communications, but they are difficult to estimate at high latitudes because of auroras. We built an artificial intelligence model for northern Norway using radar observations, magnetic field measurements, geophysical indices and solar activity. It produces more accurate estimates than existing methods, even during auroral events, and can be adapted to other regions.
Devin Huyghebaert, Juha Vierinen, Björn Gustavsson, Ralph Latteck, Toralf Renkwitz, Marius Zecha, Claudia C. Stephan, J. Federico Conte, Daniel Kastinen, Johan Kero, and Jorge L. Chau
EGUsphere, https://doi.org/10.5194/egusphere-2025-2323, https://doi.org/10.5194/egusphere-2025-2323, 2025
Short summary
Short summary
The phenomena of meteors occurs at altitudes of 60–120 km and can be used to measure the neutral atmosphere. We use a large high power radar system in Norway (MAARSY) to determine changes to the atmospheric density between the years of 2016–2023 at altitudes of 85–115 km. The same day-of-year is compared, minimizing changes to the measurements due to factors other than the atmosphere. This presents a novel method by which to obtain atmospheric neutral density variations.
Spencer Mark Hatch, Ilkka Virtanen, Karl Magnus Laundal, Habtamu Wubie Tesfaw, Juha Vierinen, Devin Ray Huyghebaert, Andres Spicher, and Jens Christian Hessen
EGUsphere, https://doi.org/10.5194/egusphere-2025-1768, https://doi.org/10.5194/egusphere-2025-1768, 2025
Short summary
Short summary
This study addresses the design of next-generation incoherent scatter radar experiments used to study the ionosphere, particularly with systems that have multiple sites. We have developed a method to estimate uncertainties of measurements of plasma density, temperature, and ion drift. Our method is open-source, and helps to optimize radar configurations and assess the effectiveness of an experiment. This method ultimately serves to enhance our understanding of Earth's space environment.
Devin Huyghebaert, Björn Gustavsson, Juha Vierinen, Andreas Kvammen, Matthew Zettergren, John Swoboda, Ilkka Virtanen, Spencer M. Hatch, and Karl M. Laundal
Ann. Geophys., 43, 99–113, https://doi.org/10.5194/angeo-43-99-2025, https://doi.org/10.5194/angeo-43-99-2025, 2025
Short summary
Short summary
The EISCAT_3D radar is a new ionospheric radar under construction in the Fennoscandia region. The radar will make measurements of plasma characteristics at altitudes above approximately 60 km. The capability of the system to make these measurements at spatial scales of less than 100 m using multiple digitised signals from each of the radar antenna panels is highlighted. There are many ionospheric small-scale processes that will be further resolved using the techniques discussed here.
Theresa Rexer, Björn Gustavsson, Juha Vierinen, Andres Spicher, Devin Ray Huyghebaert, Andreas Kvammen, Robert Gillies, and Asti Bhatt
Geosci. Instrum. Method. Data Syst. Discuss., https://doi.org/10.5194/gi-2023-18, https://doi.org/10.5194/gi-2023-18, 2024
Preprint under review for GI
Short summary
Short summary
We present a second-level calibration method for electron density measurements from multi-beam incoherent scatter radars. It is based on the well-known Flat field correction method used in imaging and photography. The methods improve data quality and useability as they account for unaccounted, and unpredictable variations in the radar system. This is valuable for studies where inter-beam calibration is important such as studies of polar cap patches, plasma irregularities and turbulence.
Johann Stamm, Juha Vierinen, Björn Gustavsson, and Andres Spicher
Ann. Geophys., 41, 55–67, https://doi.org/10.5194/angeo-41-55-2023, https://doi.org/10.5194/angeo-41-55-2023, 2023
Short summary
Short summary
The study of some ionospheric events benefit from the knowledge of how the physics varies over a volume and over time. Examples are studies of aurora or energy deposition. With EISCAT3D, measurements of ion velocity vectors in a volume will be possible for the first time. We present a technique that uses a set of such measurements to estimate electric field and neutral wind. The technique relies on adding restrictions to the estimates. We successfully consider restrictions based on physics.
Knut Ola Dølven, Juha Vierinen, Roberto Grilli, Jack Triest, and Bénédicte Ferré
Geosci. Instrum. Method. Data Syst., 11, 293–306, https://doi.org/10.5194/gi-11-293-2022, https://doi.org/10.5194/gi-11-293-2022, 2022
Short summary
Short summary
Sensors capable of measuring rapid fluctuations are important to improve our understanding of environmental processes. Many sensors are unable to do this, due to their reliance on the transfer of the measured property, for instance a gas, across a semi-permeable barrier. We have developed a mathematical tool which enables the retrieval of fast-response signals from sensors with this type of sensor design.
Carsten Baumann, Antti Kero, Shikha Raizada, Markus Rapp, Michael P. Sulzer, Pekka T. Verronen, and Juha Vierinen
Ann. Geophys., 40, 519–530, https://doi.org/10.5194/angeo-40-519-2022, https://doi.org/10.5194/angeo-40-519-2022, 2022
Short summary
Short summary
The Arecibo radar was used to probe free electrons of the ionized atmosphere between 70 and 100 km altitude. This is also the altitude region were meteors evaporate and form secondary particulate matter, the so-called meteor smoke particles (MSPs). Free electrons attach to these MSPs when the sun is below the horizon and cause a drop in the number of free electrons, which are the subject of these measurements. We also identified a different number of free electrons during sunset and sunrise.
Derek McKay, Juha Vierinen, Antti Kero, and Noora Partamies
Geosci. Instrum. Method. Data Syst., 11, 25–35, https://doi.org/10.5194/gi-11-25-2022, https://doi.org/10.5194/gi-11-25-2022, 2022
Short summary
Short summary
When radio waves from our galaxy enter the Earth's atmosphere, they are absorbed by electrons in the upper atmosphere. It was thought that by measuring the amount of absorption, it would allow the height of these electrons in the atmosphere to be determined. If so, this would have significance for future instrument design. However, this paper demonstrates that it is not possible to do this, but it does explain how multiple-frequency measurements can nevertheless be useful.
Ryan Volz, Jorge L. Chau, Philip J. Erickson, Juha P. Vierinen, J. Miguel Urco, and Matthias Clahsen
Atmos. Meas. Tech., 14, 7199–7219, https://doi.org/10.5194/amt-14-7199-2021, https://doi.org/10.5194/amt-14-7199-2021, 2021
Short summary
Short summary
We introduce a new way of estimating winds in the upper atmosphere (about 80 to 100 km in altitude) from the observed Doppler shift of meteor trails using a statistical method called Gaussian process regression. Wind estimates and, critically, the uncertainty of those estimates can be evaluated smoothly (i.e., not gridded) in space and time. The effective resolution is set by provided parameters, which are limited in practice by the number density of the observed meteors.
Johann Stamm, Juha Vierinen, and Björn Gustavsson
Ann. Geophys., 39, 961–974, https://doi.org/10.5194/angeo-39-961-2021, https://doi.org/10.5194/angeo-39-961-2021, 2021
Short summary
Short summary
Measurements of the electric field and neutral wind in the ionosphere are important for understanding energy flows or electric currents. With incoherent scatter radars (ISRs), we can measure the velocity of the ions, which depends on both the electrical field and the neutral wind. In this paper, we investigate methods to use ISR data to find reasonable values for both parameters. We find that electric field can be well measured down to 125 km height and neutral wind below this height.
Torbjørn Tveito, Juha Vierinen, Björn Gustavsson, and Viswanathan Lakshmi Narayanan
Ann. Geophys., 39, 427–438, https://doi.org/10.5194/angeo-39-427-2021, https://doi.org/10.5194/angeo-39-427-2021, 2021
Short summary
Short summary
This work explores the role of EISCAT 3D as a tool for planetary mapping. Due to the challenges inherent in detecting the signals reflected from faraway bodies, we have concluded that only the Moon is a viable mapping target. We estimate the impact of the ionosphere on lunar mapping, concluding that its distorting effects should be easily manageable. EISCAT 3D will be useful for mapping the lunar nearside due to its previously unused frequency (233 MHz) and its interferometric capabilities.
Johann Stamm, Juha Vierinen, Juan M. Urco, Björn Gustavsson, and Jorge L. Chau
Ann. Geophys., 39, 119–134, https://doi.org/10.5194/angeo-39-119-2021, https://doi.org/10.5194/angeo-39-119-2021, 2021
Harikrishnan Charuvil Asokan, Jorge L. Chau, Raffaele Marino, Juha Vierinen, Fabio Vargas, Juan Miguel Urco, Matthias Clahsen, and Christoph Jacobi
Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2020-974, https://doi.org/10.5194/acp-2020-974, 2020
Preprint withdrawn
Short summary
Short summary
This paper explores the dynamics of gravity waves and turbulence present in the mesosphere and lower thermosphere (MLT) region. We utilized two different techniques on meteor radar observations and simulations to obtain power spectra at different horizontal scales. The techniques are applied to a special campaign conducted in northern Germany in November 2018. The study revealed the dominance of large-scale structures with horizontal scales larger than 500 km during the campaign period.
Kian Sartipzadeh, Andreas Kvammen, Björn Gustavsson, Njål Gulbrandsen, Magnar Gullikstad Johnsen, Devin Huyghebaert, and Juha Vierinen
EGUsphere, https://doi.org/10.5194/egusphere-2025-3070, https://doi.org/10.5194/egusphere-2025-3070, 2025
Short summary
Short summary
Knowing charged particle densities high above Earth is key for forecasting space weather effects on satellites and communications, but they are difficult to estimate at high latitudes because of auroras. We built an artificial intelligence model for northern Norway using radar observations, magnetic field measurements, geophysical indices and solar activity. It produces more accurate estimates than existing methods, even during auroral events, and can be adapted to other regions.
Devin Huyghebaert, Juha Vierinen, Björn Gustavsson, Ralph Latteck, Toralf Renkwitz, Marius Zecha, Claudia C. Stephan, J. Federico Conte, Daniel Kastinen, Johan Kero, and Jorge L. Chau
EGUsphere, https://doi.org/10.5194/egusphere-2025-2323, https://doi.org/10.5194/egusphere-2025-2323, 2025
Short summary
Short summary
The phenomena of meteors occurs at altitudes of 60–120 km and can be used to measure the neutral atmosphere. We use a large high power radar system in Norway (MAARSY) to determine changes to the atmospheric density between the years of 2016–2023 at altitudes of 85–115 km. The same day-of-year is compared, minimizing changes to the measurements due to factors other than the atmosphere. This presents a novel method by which to obtain atmospheric neutral density variations.
Spencer Mark Hatch, Ilkka Virtanen, Karl Magnus Laundal, Habtamu Wubie Tesfaw, Juha Vierinen, Devin Ray Huyghebaert, Andres Spicher, and Jens Christian Hessen
EGUsphere, https://doi.org/10.5194/egusphere-2025-1768, https://doi.org/10.5194/egusphere-2025-1768, 2025
Short summary
Short summary
This study addresses the design of next-generation incoherent scatter radar experiments used to study the ionosphere, particularly with systems that have multiple sites. We have developed a method to estimate uncertainties of measurements of plasma density, temperature, and ion drift. Our method is open-source, and helps to optimize radar configurations and assess the effectiveness of an experiment. This method ultimately serves to enhance our understanding of Earth's space environment.
Devin Huyghebaert, Björn Gustavsson, Juha Vierinen, Andreas Kvammen, Matthew Zettergren, John Swoboda, Ilkka Virtanen, Spencer M. Hatch, and Karl M. Laundal
Ann. Geophys., 43, 99–113, https://doi.org/10.5194/angeo-43-99-2025, https://doi.org/10.5194/angeo-43-99-2025, 2025
Short summary
Short summary
The EISCAT_3D radar is a new ionospheric radar under construction in the Fennoscandia region. The radar will make measurements of plasma characteristics at altitudes above approximately 60 km. The capability of the system to make these measurements at spatial scales of less than 100 m using multiple digitised signals from each of the radar antenna panels is highlighted. There are many ionospheric small-scale processes that will be further resolved using the techniques discussed here.
Theresa Rexer, Björn Gustavsson, Juha Vierinen, Andres Spicher, Devin Ray Huyghebaert, Andreas Kvammen, Robert Gillies, and Asti Bhatt
Geosci. Instrum. Method. Data Syst. Discuss., https://doi.org/10.5194/gi-2023-18, https://doi.org/10.5194/gi-2023-18, 2024
Preprint under review for GI
Short summary
Short summary
We present a second-level calibration method for electron density measurements from multi-beam incoherent scatter radars. It is based on the well-known Flat field correction method used in imaging and photography. The methods improve data quality and useability as they account for unaccounted, and unpredictable variations in the radar system. This is valuable for studies where inter-beam calibration is important such as studies of polar cap patches, plasma irregularities and turbulence.
Johann Stamm, Juha Vierinen, Björn Gustavsson, and Andres Spicher
Ann. Geophys., 41, 55–67, https://doi.org/10.5194/angeo-41-55-2023, https://doi.org/10.5194/angeo-41-55-2023, 2023
Short summary
Short summary
The study of some ionospheric events benefit from the knowledge of how the physics varies over a volume and over time. Examples are studies of aurora or energy deposition. With EISCAT3D, measurements of ion velocity vectors in a volume will be possible for the first time. We present a technique that uses a set of such measurements to estimate electric field and neutral wind. The technique relies on adding restrictions to the estimates. We successfully consider restrictions based on physics.
Knut Ola Dølven, Juha Vierinen, Roberto Grilli, Jack Triest, and Bénédicte Ferré
Geosci. Instrum. Method. Data Syst., 11, 293–306, https://doi.org/10.5194/gi-11-293-2022, https://doi.org/10.5194/gi-11-293-2022, 2022
Short summary
Short summary
Sensors capable of measuring rapid fluctuations are important to improve our understanding of environmental processes. Many sensors are unable to do this, due to their reliance on the transfer of the measured property, for instance a gas, across a semi-permeable barrier. We have developed a mathematical tool which enables the retrieval of fast-response signals from sensors with this type of sensor design.
Carsten Baumann, Antti Kero, Shikha Raizada, Markus Rapp, Michael P. Sulzer, Pekka T. Verronen, and Juha Vierinen
Ann. Geophys., 40, 519–530, https://doi.org/10.5194/angeo-40-519-2022, https://doi.org/10.5194/angeo-40-519-2022, 2022
Short summary
Short summary
The Arecibo radar was used to probe free electrons of the ionized atmosphere between 70 and 100 km altitude. This is also the altitude region were meteors evaporate and form secondary particulate matter, the so-called meteor smoke particles (MSPs). Free electrons attach to these MSPs when the sun is below the horizon and cause a drop in the number of free electrons, which are the subject of these measurements. We also identified a different number of free electrons during sunset and sunrise.
Derek McKay, Juha Vierinen, Antti Kero, and Noora Partamies
Geosci. Instrum. Method. Data Syst., 11, 25–35, https://doi.org/10.5194/gi-11-25-2022, https://doi.org/10.5194/gi-11-25-2022, 2022
Short summary
Short summary
When radio waves from our galaxy enter the Earth's atmosphere, they are absorbed by electrons in the upper atmosphere. It was thought that by measuring the amount of absorption, it would allow the height of these electrons in the atmosphere to be determined. If so, this would have significance for future instrument design. However, this paper demonstrates that it is not possible to do this, but it does explain how multiple-frequency measurements can nevertheless be useful.
Ryan Volz, Jorge L. Chau, Philip J. Erickson, Juha P. Vierinen, J. Miguel Urco, and Matthias Clahsen
Atmos. Meas. Tech., 14, 7199–7219, https://doi.org/10.5194/amt-14-7199-2021, https://doi.org/10.5194/amt-14-7199-2021, 2021
Short summary
Short summary
We introduce a new way of estimating winds in the upper atmosphere (about 80 to 100 km in altitude) from the observed Doppler shift of meteor trails using a statistical method called Gaussian process regression. Wind estimates and, critically, the uncertainty of those estimates can be evaluated smoothly (i.e., not gridded) in space and time. The effective resolution is set by provided parameters, which are limited in practice by the number density of the observed meteors.
Johann Stamm, Juha Vierinen, and Björn Gustavsson
Ann. Geophys., 39, 961–974, https://doi.org/10.5194/angeo-39-961-2021, https://doi.org/10.5194/angeo-39-961-2021, 2021
Short summary
Short summary
Measurements of the electric field and neutral wind in the ionosphere are important for understanding energy flows or electric currents. With incoherent scatter radars (ISRs), we can measure the velocity of the ions, which depends on both the electrical field and the neutral wind. In this paper, we investigate methods to use ISR data to find reasonable values for both parameters. We find that electric field can be well measured down to 125 km height and neutral wind below this height.
Torbjørn Tveito, Juha Vierinen, Björn Gustavsson, and Viswanathan Lakshmi Narayanan
Ann. Geophys., 39, 427–438, https://doi.org/10.5194/angeo-39-427-2021, https://doi.org/10.5194/angeo-39-427-2021, 2021
Short summary
Short summary
This work explores the role of EISCAT 3D as a tool for planetary mapping. Due to the challenges inherent in detecting the signals reflected from faraway bodies, we have concluded that only the Moon is a viable mapping target. We estimate the impact of the ionosphere on lunar mapping, concluding that its distorting effects should be easily manageable. EISCAT 3D will be useful for mapping the lunar nearside due to its previously unused frequency (233 MHz) and its interferometric capabilities.
Johann Stamm, Juha Vierinen, Juan M. Urco, Björn Gustavsson, and Jorge L. Chau
Ann. Geophys., 39, 119–134, https://doi.org/10.5194/angeo-39-119-2021, https://doi.org/10.5194/angeo-39-119-2021, 2021
Harikrishnan Charuvil Asokan, Jorge L. Chau, Raffaele Marino, Juha Vierinen, Fabio Vargas, Juan Miguel Urco, Matthias Clahsen, and Christoph Jacobi
Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2020-974, https://doi.org/10.5194/acp-2020-974, 2020
Preprint withdrawn
Short summary
Short summary
This paper explores the dynamics of gravity waves and turbulence present in the mesosphere and lower thermosphere (MLT) region. We utilized two different techniques on meteor radar observations and simulations to obtain power spectra at different horizontal scales. The techniques are applied to a special campaign conducted in northern Germany in November 2018. The study revealed the dominance of large-scale structures with horizontal scales larger than 500 km during the campaign period.
Cited articles
Bust, G. S. and Mitchell, C. N.:
History, current state, and future directions of ionospheric imaging,
Rev. Geophys., 46, 1–23, 2008.
Carrano, C. S. and Groves, K.:
The GPS Segment of the AFRL-SCINDA Global Network and the Challenges of Real-Time TEC Estimation in the Equatorial Ionosphere,
Proceedings of the 2006 National Technical Meeting of The Institute of Navigation, Monterey,
CA, 2006.
Coster, A., Williams, J., Weatherwax, A., Rideout, W., and Herne, D.:
Accuracy of GPS total electron content: GPS receiver bias temperature
dependence, Radio Sci., 48, 190–196, https://doi.org/10.1002/rds.20011, 2013.
Coster, A. J., Gaposchkin, E. M., and Thornton, L. E.:
Real-time ionospheric monitoring system using the GPS,
MIT Lincoln Laboratory, Technical Report, 954, 1992.
Davies, K.:
Ionospheric Radio Propagation,
National Bureau of Standards, 278–279, 1965.
Dyrud, L., Jovancevic, A., Brown, A., Wilson, D., and Ganguly, S.:
Ionospheric measurement with GPS: Receiver techniques and methods, Radio Sci., 43, RS6002, https://doi.org/10.1029/2007RS003770, 2008.
Feltens, J.:
Chapman profile approach for 3-D global TEC representation, IGS Presentation, 1998.
Gaposchkin, E. M. and Coster, A. J.:
GPS L1-1,2 Bias Determination,
Lincoln Laboratory Technical Report, 971, 1993.
Kaipio, J. and Somersalo, E.:
Statistical and Computational Inverse Problems,
Springer, New York, USA, 2005.
Komjathy, A., Sparks, L., Wilson, B. D., and Mannucci, A. J.: Automated
daily processing of more than 1000 ground-based GPS receivers for studying
intense ionospheric storms, Radio Sci., 40, RS6006, https://doi.org/10.1029/2005RS003279, 2005.
Mannucci, A. J., Wilson, B. D., Yuan, D. N., Ho, C. H., Lindqwister, U. J., and
Runge, T. F.: A global mapping technique for GPS-derived ionospheric
total electron content measurements, Radio Sci., 33, 565–582, 1998.
Paige, C. C. and Michael, A.: Saunders,
LSQR: An algorithm for sparse linear equations and sparse least squares,
ACM Transactions on Mathematical Software (TOMS), 8, 43–71, 1982.
Rideout, W. and Coster, A.: Automated GPS processing for global
total electron content data, GPS Solutions, 10, 219–228, 2006.
Vierinen, J., Norberg, J., Lehtinen, M. S., Amm, O., Roininen, L., Väänänen, A.,
and McKay-Bukowski, D. J.: Beacon satellite receiver for ionospheric tomography, Radio Sci.,
49,
1141–1152, 2014.
Short summary
We present a statistical framework for estimating GNSS receiver bias by using a weighted linear least squares of independent differences (WLLSID) model to examine differences of a large number of TEC measurements. This allows a consistent way for treating elevation-dependent model errors and spatiotemporal distance-dependent geophysical differences arising in ionospheric GNSS measurements. The method is also applicable to other GNSS system than GPS, supporting, e.g., GLONASS.
We present a statistical framework for estimating GNSS receiver bias by using a weighted linear...