Trapped Proton Environment in Medium-Earth Orbit (2000-2010)

This report presents a new empirical model derived from Polar mission data and scaled by GPS in-situ measurements to characterize the dynamic trapped proton environment in Medium-Earth Orbit (2000–2010), demonstrating that it provides more accurate statistical flux estimates and captures temporal variations better than the standard AP8 model.

The Gist

The Big Picture: The GPS "Highway" in a Stormy Sea

Imagine the Earth is surrounded by a giant, invisible, and very dangerous storm belt. This isn't a storm of rain and wind, but a storm of high-speed, tiny particles called protons. These particles are trapped by Earth's magnetic field, swirling around like bees in a hive.

This report is about a specific stretch of space called Medium-Earth Orbit (MEO). This is where our GPS satellites live. They fly at about 20,000 km up, right through the middle of this proton storm.

The Problem:
GPS satellites are like precious cargo ships. If they get hit by too many of these "proton bees," their electronics can get fried, leading to navigation failures. To protect them, scientists need to know exactly how dangerous the storm is.

However, for a long time, the scientists' "weather maps" were outdated. They were using a map from the 1970s (called the AP8 model) that was based on old data. It was like trying to navigate a modern city using a map from the 1950s; it missed new roads and didn't account for traffic jams.

The Solution: A New Map and a Real-Time Sensor

The authors of this report (Yue Chen and colleagues from Los Alamos National Laboratory) decided to build a brand-new, better map and combine it with real-time data. They used a 3-step recipe to do this:

Step 1: Building the "PolarP" Map (The Historical Archive)

First, they looked at a massive library of data collected by the Polar satellite between 1996 and 2007.

• The Analogy: Imagine the Polar satellite was a photographer taking millions of pictures of the proton storm over a decade.
• The Result: They turned these photos into a new statistical model called PolarP. Unlike the old map, which just showed the "average" storm, this new map shows the range of possibilities. It tells you not just what the weather is usually like, but what it looks like on the worst possible days (the "90th percentile"). It's like knowing that while it usually rains lightly, sometimes it pours.

Step 2: The "GPS ns41" Sensor (The Real-Time Check)

Next, they looked at a specific GPS satellite, ns41, which has been flying since 2000.

• The Analogy: This satellite is like a weather station floating right in the middle of the storm. It has a sensor that counts the protons hitting it.
• The Catch: This sensor is a bit limited. It can only "see" protons in a narrow energy range (like a camera that only takes photos in black and white, but the storm is full of colors). It can't see the whole spectrum of particle energies.

Step 3: Mixing the Map with the Sensor (The Magic Sauce)

This is the clever part. The scientists took the PolarP map (which has all the colors and details but is just a statistical average) and scaled it using the ns41 sensor (which is real-time but limited).

• The Analogy: Imagine you have a high-definition, 3D movie of a storm (PolarP), but you only have a black-and-white thermometer reading from a specific spot (ns41).
• The Fix: You adjust the brightness and contrast of the 3D movie so that the temperature in the movie matches the thermometer reading exactly. Now, you have a 3D movie that is both detailed (from the Polar data) and accurate for that specific time (from the GPS data).

What Did They Find?

1. The Old Map Was Wrong: The old AP8 model was too simple. It assumed the storm was steady and predictable. The new data shows the storm is wild and dynamic. Sometimes the proton levels are much higher than the old map predicted, and sometimes lower.
2. The "Error Factor" Myth: Scientists used to think the old map was only off by a factor of 2 (e.g., if it said 10 protons, the real number was between 5 and 20). This report shows that the reality is much more complex. The variations are huge, and the old "factor of 2" rule is an oversimplification that could lead to underestimating the danger.
3. Solar Cycles Matter: They noticed that the storm gets worse or better depending on the Sun's activity. During quiet times, the storm builds up; during solar storms, the particles get knocked around. The new model captures these changes day-by-day.

Why Does This Matter?

This report provides a daily log of the proton storm from 2000 to 2010.

• For Engineers: It helps them design better shields for future satellites so they don't get fried.
• For Scientists: It proves that space weather isn't static; it changes constantly, and we need models that reflect that chaos.

The Bottom Line

The authors took old, detailed data (Polar) and mixed it with fresh, real-time data (GPS) to create a super-accurate, dynamic weather forecast for the proton belt. They showed that the old maps were too simple and that the space environment around our GPS satellites is much more unpredictable and powerful than we previously thought.

Uncertainty Note: Even with this new method, there is still some guesswork. The "weather forecast" might be off by a factor of 3 due to the model, or a factor of 5 due to the sensor limitations. But it is a much better forecast than what we had before.

Technical Summary

1. Problem Statement

The Medium-Earth Orbit (MEO) region, particularly the Global Positioning System (GPS) orbit at ~20,000 km altitude, is a critical zone for navigation and communication satellites. However, it remains underexplored due to the intense radiation environment of the trapped proton and electron belts.

• Limitations of Existing Data: In-situ measurements from GPS satellites (specifically the BDD-IIR instrument) have limited energy coverage (typically >1 MeV), failing to provide a full-spectrum description of the radiation belt.
• Limitations of Existing Models: The de-facto standard empirical model, AP8 (based on 1960s–70s data), suffers from significant deficiencies:
  • It relies on outdated data with limited L-shell coverage.
  • It provides only long-term mean values, failing to capture the dynamic temporal variations of proton fluxes modulated by geomagnetic storms.
  • The commonly accepted "error factor of 2" for AP8 is an oversimplification that underestimates the true variability of the proton belt, potentially leading to incorrect radiation dose accumulation estimates.

2. Methodology

The authors developed a 3-step hybrid method to derive dynamic, wide-energy-range proton fluxes along the GPS ns41 orbit for the period 2000–2010 (covering nearly a full solar cycle).

• Step 1: Development of the "PolarP" Empirical Model

  • Data Source: Long-term (1996–2007) in-situ measurements from the CEPPAD experiment on the Polar satellite.
  • Scope: The model covers a wide energy range (20 keV to >10 MeV), L-shells (up to ~10), and equatorial pitch angles.
  • Statistical Approach: Unlike AP8, PolarP generates occurrence distributions (percentiles) rather than just mean values. It calculates flux values at specific percentiles (e.g., 50th, 90th) to represent worst-case scenarios and statistical variability.
  • Validation: The model was validated against quiet and disturbed days, confirming that observed fluxes fall within the predicted percentile ranges.

• Step 2: Derivation of Scaling Ratios

  • The authors utilized GPS ns41 in-situ measurements (specifically the lowest energy channel: 1.27–5.3 MeV) to calculate a time-dependent scaling ratio ($R$).
  • Formula: $R(t) = \frac{J_{measured}(t)}{J_{PolarP}(t)}$
  • This ratio accounts for the dynamic nature of the outer radiation belt during geomagnetic storms, which the static statistical model (PolarP) alone cannot capture.

• Step 3: Generation of Extended Fluxes

  • Assuming the spectral shape derived from the PolarP model is valid (due to the limited energy coverage of GPS), the authors scaled the PolarP model output across the entire energy spectrum (50 keV – 6 MeV) using the time-dependent ratio $R(t)$.
  • Equation: $J_{derived}(E, t) = J_{PolarP}(E, t) \times R(t)$
  • This process preserves the temporal features from the GPS measurements while expanding the energy coverage using the PolarP spectral shape.

3. Key Contributions

• New Empirical Model (PolarP): Introduced a modern, statistically robust proton radiation belt model based on Polar/CEPPAD data, offering percentile-based flux predictions (e.g., 90th percentile for worst-case analysis) rather than simple averages.
• Dynamic Flux Derivation: Successfully generated a 10-year time series of daily proton fluxes for the GPS orbit that captures temporal variations driven by geomagnetic activity, overcoming the static nature of AP8.
• Wide Energy Coverage: Bridged the gap between low-energy (50 keV) and high-energy (6 MeV) protons, providing a complete spectral description previously unavailable for MEO.
• Inter-Instrument Calibration: Identified and addressed systematic differences between the Polar CEPPAD instrument and the GPS BDD-IIR instrument, noting that CEPPAD measurements were systematically 2–10 times higher than GPS measurements. The study chose to trust and scale the GPS data due to its long-term continuity in MEO.

4. Key Results

• Discrepancy with AP8: Comparisons between PolarP and AP8 revealed significant differences.
  • For protons < 1 MeV, AP8 fluxes were significantly higher than PolarP.
  • For protons > 1.5 MeV, AP8 fluxes were significantly lower than PolarP.
  • The differences were often orders of magnitude, suggesting AP8 is unreliable for accurate dose accumulation in MEO.
• Temporal Variability: The derived daily fluxes showed that proton levels fluctuate significantly. For example, the average proton flux in the second half of the study period (2005–2010) was ~5 times higher than the first half, attributed to the solar minimum reducing storm-induced de-trapping of protons.
• Fluence Accumulation: Calculations of accumulated fluence (2000–2010) showed that AP8 (a static model) predicts a constant accumulation rate. In contrast, the new model showed variable accumulation rates, with high-energy fluence eventually exceeding AP8 predictions in later years.
• Error Factors:
  • Uncertainty from the PolarP model: Error factor < ~3.
  • Uncertainty from in-situ measurements (instrument calibration differences): Error factor ~5.

5. Significance

• Satellite Protection: The derived daily flux data files provide a more accurate basis for designing radiation-hardened electronics and predicting single-event upsets (SEUs) for GPS and other MEO satellites.
• Operational Safety: By providing worst-case scenarios (e.g., 90th percentile) and dynamic time-series data, operators can better predict radiation hazards during geomagnetic storms.
• Model Evolution: The report demonstrates that the "error factor of 2" convention for AP8 is insufficient. It establishes a new methodology for combining statistical models with real-time in-situ scaling to create dynamic radiation environment specifications.
• Data Availability: The report delivers the actual daily proton flux data files for the GPS ns41 orbit (2000–2010), serving as a valuable resource for the space weather and satellite engineering communities.

In conclusion, this report moves beyond the static limitations of the AP8 model by integrating modern statistical modeling with long-term in-situ observations, offering a significantly more accurate and dynamic characterization of the trapped proton environment in Medium-Earth Orbit.