🔭 astrophysics

Approximation of the reception coefficients of cosmic rays neutron component for latitude measurement

This study develops and validates a high-accuracy method for approximating cosmic ray reception coefficients as a function of latitude, enabling the correction of monitoring data to facilitate precise latitude measurements during marine expeditions.

Original authors: Kobelev P. G., Yanke V. G

Published 2026-02-02

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Kobelev P. G., Yanke V. G

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Tuning the Radio to Hear the Stars

Imagine the Earth is surrounded by a constant, invisible rain of tiny particles from deep space called Cosmic Rays. Most of the time, this rain falls evenly from all directions, like a gentle, uniform drizzle. However, sometimes the "wind" in space (the solar wind) pushes these particles, creating a slight tilt or "anisotropy" in the rain—like the rain seeming to come more from the east than the west.

Scientists want to measure this tilt to understand what's happening in our solar system. But here's the problem: Earth has a giant magnetic shield (the magnetosphere) and a thick blanket of air (the atmosphere). These act like a complex filter. Depending on where you are standing on Earth (your latitude), your "window" to space looks different. A detector in the Arctic sees a very different slice of the cosmic ray sky than a detector near the equator.

The Problem: The "Translation" Gap

For decades, scientists have had a network of detectors (Neutron Monitors) all over the globe. They know exactly how each specific detector reacts to the cosmic rays. However, these calculations were only done for specific, fixed locations (like a city or a mountain peak).

The authors of this paper faced a practical problem: What if you are on a moving ship?
Imagine a research vessel sailing from Vladivostok to Kaliningrad. As the ship moves, its "view" of the cosmic rays changes constantly. To analyze the data correctly, the scientists needed to know the "reception coefficient" (a fancy way of saying "sensitivity factor") for any point along that route, not just the fixed cities where detectors happen to sit.

Previously, if a ship moved to a spot without a nearby detector, scientists had to do incredibly difficult, manual calculations to figure out how that specific spot would react. It was like trying to translate a book by hand, word by word, every time you moved to a new room.

The Solution: Creating a "Universal Translator"

The goal of this paper was to create a mathematical "shortcut" or a universal formula.

The authors took all the existing, hard-won data from the global network of detectors and asked: "Can we find a smooth curve that connects all these dots?"

They used a mathematical tool called the Granitsky–Dorman function. Think of this function as a master key. Instead of needing a unique key for every single lock (every specific location on Earth), they created one master key that fits almost all locks, provided you know the "lock size" (which, in this case, is the geomagnetic cutoff rigidity—a measure of how strong the Earth's magnetic shield is at that specific spot).

How They Did It (The Recipe)

Gathered the Data: They looked at the known sensitivity of about 100 detectors around the world.
Found the Pattern: They noticed that a detector's sensitivity changes predictably based on how strong the local magnetic shield is.
Built the Model: They fitted their "master key" formula to this data.
- Zero Harmonic (The General Rain): They modeled how sensitive detectors are to the overall amount of cosmic rays. This was a very good fit (97% accuracy).
- First Harmonic (The Directional Wind): They modeled how sensitive detectors are to the direction the particles are coming from. This was trickier because the Earth's magnetic field splits the view into "North" and "South" branches, but they still found a good formula for it.
- Drift Angle: They also figured out how the Earth's magnetic field twists the path of the particles, creating a "drift" angle, and modeled that too.

The Results: A Map for Moving Ships

The paper presents a set of equations (found in Table 1 of the original text) that allow scientists to plug in a location's magnetic "strength" (rigidity) and instantly get the correct sensitivity factors.

They tested this new method on a real-world scenario: a 4-month voyage of the research vessel Akademik Kurchatov.

The Test: They compared the new formula's predictions against actual data from detectors located near the ship's path.
The Outcome: The formula worked. The predictions matched the real-world data almost perfectly, staying well within the margin of error.

Why This Matters (According to the Paper)

New Detectors: Since 1982, about 20 new detectors have been built. This new formula allows scientists to instantly calculate how these new, untested detectors will behave without waiting years to gather data.
Marine Expeditions: It solves the "moving ship" problem. Now, when a ship sails across the ocean, scientists can accurately strip away the Earth's magnetic influence to see the true cosmic ray variations happening in space.

In Summary

The authors took a complex, location-specific puzzle and solved it by finding a simple, smooth mathematical pattern. They turned a process that required laborious, point-by-point calculations into a simple "plug-and-play" formula. This allows scientists to accurately measure cosmic rays from anywhere on Earth, whether standing still in a lab or sailing across the ocean.

Technical Summary: Approximation of the Reception Coefficients of Cosmic Rays Neutron Component for Latitude Measurement

Problem Statement
In the study of Galactic Cosmic Rays (CR), particularly for latitude monitoring during marine expeditions, it is essential to isolate primary CR variations from observation data. While the Global Survey Method (GSM) allows for the determination of CR density and anisotropy outside the magnetosphere, applying these results to specific, arbitrary locations (such as a moving vessel) requires solving a "direct problem." This involves calculating expected variations at a specific point based on known anisotropy vectors.

Historically, this calculation relied on reception coefficients (linking variations outside the magnetosphere to ground observations) calculated by Yasue et al. (1982) only for discrete points within the existing global network of neutron monitors. Since approximately 20 new neutron detectors have been established since 1982, and marine expeditions traverse arbitrary latitudes and longitudes, the existing discrete data is insufficient. Performing laborious trajectory calculations for every new point is inefficient. Therefore, there is a need to approximate these reception coefficients as a continuous function of geomagnetic cutoff rigidity ( $R_c$ ) to facilitate real-time data correction for latitude measurements and new stations.

Methodology
The authors utilized the existing reception coefficients for zero, first, and second harmonics calculated by Yasue et al. (1982) for the global network of neutron monitors. These coefficients were originally derived for discrete detector locations and specific solar activity (SA) conditions (minimum and maximum).

To create a continuous model, the authors employed the Granitsky–Dorman approximation function, a mathematical model previously proven effective for processing CR latitude measurements. The specific functional forms used were:

Zero Harmonic ( $C_{00}$ ): Modeled as $C_{00} = c [1 - \exp(-\alpha R^{-\kappa})]^{-1}$ .
First Harmonic ( $C_{10}$ and $A_{11}$ ): Modeled using variations of the Granitsky–Dorman function, including terms for the equatorial component amplitude ( $A_{11}$ ) and the effective particle drift angle ( $\phi_{11}$ ).
North-South Asymmetry ( $C_{10}$ ): Treated separately for Northern and Southern hemispheres due to the splitting of curves determined by asymptotic directions.

The non-linear approximation was performed using a Python-based software package (lmfit). The analysis covered data for both SA minimum and maximum, considering power-law spectrum indices ( $\gamma$ ) of 0.5, 1.0, and 1.5. The quality of the approximation was assessed using the coefficient of determination ( $R^2$ ) and the chi-squared ( $\chi^2$ ) criterion.

Key Results
The study successfully generated continuous approximations for the reception coefficients as functions of geomagnetic cutoff rigidity:

Zero Harmonic ( $C_{00}$ ): The approximation showed high accuracy with a coefficient of determination $R^2 \approx 0.97$ for SA minimum. The estimation error is approximately 2%. The model effectively describes the sensitivity of sea-level detectors, though high-altitude mountain stations were excluded from the primary fit due to depth variations.
North-South Asymmetry ( $C_{10}$ ): The approximation resulted in $R^2 \approx 0.83$ . The lower value is attributed to the limited number of high-latitude stations sensitive to this specific variation type. The estimation error is approximately 8%.
Equatorial Component Amplitude ( $A_{11}$ ): This proved the most difficult to approximate, yielding $R^2 \approx 0.71$ . The authors note that while the $R^2$ is lower, the model remains significant as the parameters represent average values, and the residuals do not show systematic bias. The estimation error is approximately 10%.
Drift Angle ( $\phi_{11}$ ): The approximation for the effective drift angle achieved $R^2 \approx 0.79$ with an estimation error of about 5%.

The authors validated these models using data from the 34th expedition of the research vessel Akademik Kurchatov (Vladivostok to Kaliningrad, 1981–1982). The calculated coefficients along the vessel's route agreed with nearby station data within the error margins, with minor deviations at specific Japanese stations attributed to calculation errors in the original discrete data.

Significance and Claims
The paper claims that this work represents the first time reception coefficients for the first harmonic have been approximated as a function of geomagnetic cutoff rigidity. The primary significance of this development is twofold:

New Stations: It allows for the estimation of reception coefficients for the approximately 20 new neutron detectors established since 1982, which previously lacked calculated coefficients.
Marine Expeditions: It enables the processing of latitude measurement data from marine expeditions by providing a method to correct observation data for primary CR variations at arbitrary points along a vessel's route, rather than being restricted to the locations of existing ground stations.

The authors maintain a modest tone regarding the statistical fit, acknowledging that while some coefficients of determination are lower (e.g., 0.71 for $A_{11}$ ), the models are adequate for their intended application, as evidenced by the random nature of the residuals and successful validation against expedition data. The work does not propose new experimental setups but rather provides a necessary mathematical tool for correcting existing observational data.