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

기술 요약: 위도 측정을 위한 우주선 중성자 성분의 수신 계수 근사

문제 제기
은하 우주선(CR) 연구, 특히 해양 탐사 중의 위도 모니터링에 있어서는 관측 데이터로부터 1차 우주선의 변화를 분리해내는 것이 필수적이다. 전 지구적 조사 방법(GSM)을 통해 자기권 외부의 우주선 밀도와 이방성을 결정할 수 있지만, 이를 특정 임의의 위치(예: 이동 중인 선박)에 적용하기 위해서는 "직접 문제(direct problem)"를 해결해야 한다. 이는 알려진 이방성 벡터를 바탕으로 특정 지점에서의 예상 변화량을 계산하는 과정을 포함한다.

역사적으로 이러한 계산은 야스에 등(Yasue et al., 1982)이 기존 중성자 검출기 네트워크 내의 이산적인 지점들에 대해서만 계산한 수신 계수(자기권 외부의 변화와 지상 관측을 연결하는 계수)에 의존해 왔다. 1982년 이후 약 20개의 새로운 중성자 검출기가 설치되었고, 해양 탐사는 임의의 위도와 경도를 통과하므로 기존의 이산적인 데이터는 불충분하다. 매번 새로운 지점에 대해 노동 집약적인 궤적 계산을 수행하는 것은 비효율적이다. 따라서 위도 측정 및 신규 관측소에 대한 실시간 데이터 보정을 용이하게 하기 위해, 지자기 차단 강도($R_c$)의 연속 함수로서 수신 계수를 근사화할 필요가 있다.

방법론
저자들은 야스에 등(1982)이 전 세계 중성자 검출기 네트워크를 위해 계산한 0차, 1차 및 2차 조화 함수의 기존 수신 계수를 활용하였다. 이 계수들은 원래 이산적인 검출기 위치와 특정 태양 활동(SA) 조건(최소 및 최대)을 위해 도출된 것이다.

연속적인 모델을 만들기 위해 저자들은 우주선 위도 측정 처리에 효과적임이 이미 입증된 수학적 모델인 그라니츠키-도먼(Granitsky–Dorman) 근사 함수를 채택하였다. 사용된 구체적인 함수 형태는 다음과 같다:

• 0차 조화 함수 ($C_{00}$): $C_{00} = c [1 - \exp(-\alpha R^{-\kappa})]^{-1}$로 모델링됨.
• 1차 조화 함수 ($C_{10}$ 및 $A_{11}$): 적도 성분 진폭($A_{11}$)과 유효 입자 드리프트 각도($\phi_{11}$)에 대한 항을 포함하여 그라니츠키-도먼 함수의 변형을 사용하여 모델링됨.
• 남북 비대칭 ($C_{10}$): 점근적 방향에 의해 결정되는 곡선의 분리를 고려하여 북반구와 남반구를 별도로 취급함.

비선형 근사는 Python 기반 소프트웨어 패키지(lmfit)를 사용하여 수행되었다. 분석은 0.5, 1.0, 1.5의 멱법칙 스펙트럼 지수($\gamma$)를 고려하여 태양 활동 최소 및 최대 데이터를 모두 포함하였다. 근사의 품질은 결정 계수($R^2$)와 카이제곱($\chi^2$) 기준을 사용하여 평가되었다.

주요 결과
본 연구는 지자기 차단 강도의 함수로서 수신 계수에 대한 연속 근사치를 성공적으로 생성하였다:

1. 0차 조화 함수 ($C_{00}$): 근사는 태양 활동 최소 시 $R^2 \approx 0.97$의 높은 정확도를 보였다. 추정 오차는 약 2%이다. 이 모델은 해수면 검출기의 민감도를 효과적으로 설명하지만, 깊이 변화로 인해 고산 지대 관측소는 주요 적합(fit)에서 제외되었다.
2. 남북 비대 asymmetry ($C_{10}$): 근사 결과 $R^2 \approx 0.83$을 얻었다. 이 낮은 수치는 이 특정 변화 유형에 민감한 고위도 관측소의 수가 제한적이기 때문이다. 추정 오차는 약 8%이다.
3. 적도 성분 진폭 ($A_{11}$): 이것이 근사하기 가장 어려운 부분이었으며, $R^2 \approx 0.71$을 기록했다. 저자들은 $R^2$가 낮더라도 파라미터들이 평균값을 나타내며 잔차가 체계적인 편향을 보이지 않으므로 모델이 여전히 유의미하다고 언급하였다. 추정 오차는 약 10%이다.
4. 드리프트 각도 ($\phi_{11}$): 유효 드리프트 각도에 대한 근사는 약 5%의 추정 오차와 함께 R^2 \approx 0.79}를 달성하였다.

저자들은 연구선 아카데믹 쿠르차토프(Akademik Kurchatov) 호의 34차 탐사(1981–1982년 블라디보스토크에서 칼리닝라드까지) 데이터를 사용하여 모델을 검증하였다. 선박 경로를 따른 계산된 계수들은 인근 관측소 데이터와 오차 범위 내에서 일치하였으며, 일본의 특정 관측소에서 나타난 미세한 편차는 기존 이산 데이터의 계산 오류에 기인하였다.

의의 및 주장
본 논문은 1차 조화 함수의 수신 계수를 지자기 차단 강도의 함수로 근사한 최초의 사례라고 주장한다. 이 개발의 주요 의의는 두 가지 측면에서 나타난다:

1. 신규 관측소: 1982년 이후 설치된 약 20개의 새로운 중성자 검출기 중, 이전에 계산된 계수가 없었던 관측소들에 대한 수신 계수 추정을 가능하게 한다.
2. 해양 탐사: 선박의 경로를 따라 존재하는 임의의 지점에서 1차 우주선 변화에 대한 관측 데이터를 보정할 수 있는 방법을 제공함으로써, 해양 탐사의 위도 측정 데이터 처리를 가능하게 한다.

저자들은 통계적 적합성에 대해 겸허한 어조를 유지하며, 일부 결정 계수가 낮더라도(예: $A_{11}$의 0.71) 잔차의 무작위성과 탐사 데이터에 대한 성공적인 검증을 통해 모델이 의도된 응용 분야에 적절함을 인정하였다. 본 연구는 새로운 실험 설계를 제안하는 것이 아니라, 기존 관측 데이터를 보정하기 위한 필수적인 수학적 도구를 제공하는 것이다.