skip to main content


Title: Investigation of geomagnetic reference models based on the Iridium$$^{\circledR }$$ constellation
Abstract

The World Magnetic Model (WMM) is a geomagnetic main field model that is widely used for navigation by governments, industry and the general public. In recent years, the model has been derived using high accuracy magnetometer data from the Swarm mission. This study explores the possibility of developing future WMMs in the post-Swarm era using data from the Iridium satellite constellation. Iridium magnetometers are primarily used for attitude control, so they are not designed to produce the same level of accuracy as magnetic data from scientific missions. Iridium magnetometer errors range from 30 nT quantization to hundreds of nT errors due to spacecraft contamination and calibration uncertainty, whereas Swarm measurements are accurate to about 1 nT. The calibration uncertainty in the Iridium measurements is identified as a major error source, and a method is developed to calibrate the spacecraft measurements using data from a subset of the INTERMAGNET observatory network producing quasi-definitive data on a regular basis. After calibration, the Iridium data produced main field models with approximately 20 nT average error and 40 nT maximum error as compared to the CHAOS-7.2 model. For many scientific and precision navigation applications, highly accurate Swarm-like measurements are still necessary, however, the Iridium-based models were shown to meet the WMM error tolerances, indicating that Iridium is a viable data source for future WMMs.

Graphical Abstract

 
more » « less
Award ID(s):
2002574
NSF-PAR ID:
10363195
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Earth, Planets and Space
Volume:
74
Issue:
1
ISSN:
1880-5981
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Characterization of Earth's magnetic field is key to understanding dynamics of the core. We assess whether Iridium Communications magnetometer data can be used for this purpose since. The 66 Iridium satellites are in 86° inclination, 780 km altitude, circular orbits, with 11 satellites in each of six orbit planes. In one day the constellation returns 300,000 measurements spanning the globe with <2° spacing. We used data from January 2010 through November 2015, and compared against International Geomagnetic Reference Field (IGRF‐11) to inter‐calibrate all data to the same model. Geomagnetically quiet 24‐h intervals were selected using the total Birkeland current, auroral electrojet, and ring current indices. Thez‐scores for these quantities were combined and the quietest 16 intervals from each quarter selected for analysis. Residuals between the data and IGRF‐11 yield consistent patterns that evolve gradually from 2010 to 2015. Residuals for each day were binned in 9° latitude by 9° longitude and the distributions about the mean in each bin are Gaussian with 1‐sigma standard errors of ∼3 nT. Spherical harmonic coefficients for each quiet day were computed and time series of the coefficients used to identify artifacts at the orbit precession (8 months) and seasonal (12 months) periods and their harmonics which were then removed by notch filtering. This analysis yields time series at 800 virtual geomagnetic observatories each providing a global field map using a single day of data. The results and CHAOS 7.4 generally agree, but systematic differences larger than the statistical uncertainties are present that warrant further exploration.

     
    more » « less
  2. Abstract

    This paper studies the ionosphere's response to the annular solar eclipse on 26 December 2019, utilizing the following ground‐based and space‐borne measurements: Global Navigation Satellite System (GNSS) total electron content (TEC) data, spectral radiance data from the Sentinel‐5P satellite, in situ electron density and/or temperature measurements from DMSP and Swarm satellites, and local magnetometer data. Analysis concentrated on ionospheric effects over low‐latitude regions with respect to obscuration, local time, latitude, and altitude. The main results are as follows: (1) a local TEC reduction of4–6 TECU (30–50%) was identified along the annular eclipse path, with larger depletion and longer recovery periods in the morning eclipse compared to midday. (2) The equatorial electrojet current was significantly weakened when the eclipse trajectory crossed the magnetic equator in the morning (India) sector, which contributed to large and prolonged TEC depletion therein. (3) At midday, equatorial ionization anomaly exhibited enhancements of 20–40% as well as poleward shifting of 3–4°, likely triggered by modified neutral wind and electrodynamics patterns. (4) The behavior of equatorial ionospheric electron density showed considerable altitudinal differences in the topside, exhibiting30% reduction around 500 km and30% enhancement with 300–500 KTereduction around 850 km, before the arrival of maximum eclipse. This may have been caused by the enhanced eastward electric field and equatorward neutral wind, and other possible factors are also discussed.

     
    more » « less
  3. Abstract Aim

    Ecological properties governed by threshold relationships can exhibit heightened sensitivity to climate, creating an inherent source of uncertainty when anticipating future change. We investigated the impact of threshold relationships on our ability to project ecological change outside the observational record (e.g., the 21st century), using the challenge of predicting late‐Holocene fire regimes in boreal forest and tundra ecosystems.

    <bold>Location</bold>

    Boreal forest and tundra ecosystems of Alaska.

    <bold>Time period</bold>

    850–2100 CE.

    <bold>Major taxa studied</bold>

    Not applicable.

    Methods

    We informed a set of published statistical models, designed to predict the 30‐year probability of fire occurrence based on climatological normals, with downscaled global climate model data for 850–1850 CE. To evaluate model performance outside the observational record and the implications of threshold relationships, we compared modelled estimates with mean fire return intervals estimated from 29 published lake‐sediment palaeofire reconstructions. To place our results in the context of future change, we evaluate changes in the location of threshold to burning under 21st‐century climate projections.

    Results

    Model–palaeodata comparisons highlight spatially varying accuracy across boreal forest and tundra regions, with variability strongly related to the summer temperature threshold to burning: sites closer to this threshold exhibited larger prediction errors than sites further away from this threshold. Modifying the modern (i.e., 1950–2009) fire–climate relationship also resulted in significant changes in modelled estimates. Under 21st‐century climate projections, increasing proportions of Alaskan tundra and boreal forest will approach and surpass the temperature threshold to burning, with > 50% exceeding this threshold by > 2 °C by 2070–2099.

    <bold>Main conclusions</bold>

    Our results highlight a high sensitivity of statistical projections to changing threshold relationships and data uncertainty, implying that projections of future ecosystem change in threshold‐governed ecosystems will be accompanied by notable uncertainty. This work also suggests that ecological responses to climate change will exhibit high spatio‐temporal variability as different regions approach and surpass climatic thresholds over the 21st century.

     
    more » « less
  4. Abstract

    Machine learning (ML) has been applied to space weather problems with increasing frequency in recent years, driven by an influx of in-situ measurements and a desire to improve modeling and forecasting capabilities throughout the field. Space weather originates from solar perturbations and is comprised of the resulting complex variations they cause within the numerous systems between the Sun and Earth. These systems are often tightly coupled and not well understood. This creates a need for skillful models with knowledge about the confidence of their predictions. One example of such a dynamical system highly impacted by space weather is the thermosphere, the neutral region of Earth’s upper atmosphere. Our inability to forecast it has severe repercussions in the context of satellite drag and computation of probability of collision between two space objects in low Earth orbit (LEO) for decision making in space operations. Even with (assumed) perfect forecast of model drivers, our incomplete knowledge of the system results in often inaccurate thermospheric neutral mass density predictions. Continuing efforts are being made to improve model accuracy, but density models rarely provide estimates of confidence in predictions. In this work, we propose two techniques to develop nonlinear ML regression models to predict thermospheric density while providing robust and reliable uncertainty estimates: Monte Carlo (MC) dropout and direct prediction of the probability distribution, both using the negative logarithm of predictive density (NLPD) loss function. We show the performance capabilities for models trained on both local and global datasets. We show that the NLPD loss provides similar results for both techniques but the direct probability distribution prediction method has a much lower computational cost. For the global model regressed on the Space Environment Technologies High Accuracy Satellite Drag Model (HASDM) density database, we achieve errors of approximately 11% on independent test data with well-calibrated uncertainty estimates. Using an in-situ CHAllenging Minisatellite Payload (CHAMP) density dataset, models developed using both techniques provide test error on the order of 13%. The CHAMP models—on validation and test data—are within 2% of perfect calibration for the twenty prediction intervals tested. We show that this model can also be used to obtain global density predictions with uncertainties at a given epoch.

     
    more » « less
  5. Abstract

    This paper examined the secular displacement of the dip equator and the geomagnetic poles as well as the variations of the global magnetic inclination and declination angles using magnetometer measurements onboard different low‐Earth orbit (LEO) satellite. The secular variation of the dip equator and geomagnetic poles has different impacts on different applications—from affecting the long‐term characterization of the low‐latitude ionosphere to degrading the precision of geomagnetic navigation. The strong displacement of the dip equator can result in a systematic error in the determination of the long‐term equatorial electric field variations and hence in the characterization of ionospheric density structure, especially in the region, where the displacement of dip equator is large enough (more than 20 km or 0.2°/year) within the time scale of a solar cycle or less. Similarly, the slowly moving locations of magnetic poles, estimated from magnetometer observations onboard LEO satellite, exemplify noticeable discrepancy with that of world magnetic model (WMM) and International Geomagnetic Reference Field (IGRF) values, indicating inevitable possible impact on the precise geomagnetic navigation for commercial and military applications. Thus, accurately estimated locations of the dip equator and magnetic poles, as well as declination angles, are critically important.

     
    more » « less