JUNO 지역 지오뉴트리노 예측을 위한 고해상도 중력 모델 GIGJ

초록

본 논문은 GOCE 위성 중력 데이터와 지역 지진·지구물리학 자료를 베이즈 방식으로 결합해 광동성(Guangdong) 지역을 50 m × 50 m × 100 m 해상도의 3차원 격자(≈46 000 voxel)로 모델링하였다. 모델은 평균 1 mGal 수준의 잔차를 보이며, 기존 전 지구 모델과 달리 상·중·하부 지각의 두께와 질량을 현장 맞춤형으로 재분배한다. 그 결과 중간 지각의 지오뉴트리노 기여는 21 % 감소, 하부 지각은 24 % 증가했으며, 각 지각층별 U·Th 함량에 대한 불확실성은 55 %~78 %까지 감소하였다.

상세 분석

이 연구는 지오뉴트리노 탐지를 위한 지역 지각 모델링의 새로운 패러다임을 제시한다. 먼저, GOCE 위성에서 제공하는 전역 중력장 데이터를 5 km × 5 km 격자로 전처리하고, 이를 50 km × 50 km × 0.1 km 크기의 정규 격자(voxel)로 변환하였다. 각 voxel은 물질 라벨(UC, MC, LC 등)과 고정 밀도를 갖는 두 파라미터로 정의되며, 전방 모델은 평면 근사(planar approximation)를 사용해 계산 효율성을 확보한다.

베이즈 역학을 기반으로 한 목표 함수는 관측 중력 신호와 모델링 중력 사이의 차이를 최소화하는 likelihood와, 사전 정보(prior)인 심층 지진 탐사(DSS), 수신 함수(RF), 텔레시즘 P파 속도 모델, Moho 깊이 지도 등을 결합한다. 사전 정보는 각 지각층의 평균 P파 속도와 두께에 대한 3σ 불확실성을 반영하며, 이를 통해 역문 문제의 비선형성·불안정성을 완화한다. 최적화는 마르코프 체인 몬테카를로(MCMC)와 유전 알고리즘을 혼합한 stochastic optimization으로 수행돼, 전역 최소값뿐 아니라 파라미터 공간의 불확실성 분포도 동시에 추정한다.

결과적으로 GIGJ 모델은 GOCE 관측값과의 잔차가 평균 1 mGal(σ≈1 mGal)로, 위성 데이터의 관측 오차와 일치한다. 기존 CRUST1.0·CRUST2.0 등 전 지구 모델은 평균적인 지각 두께와 밀도를 적용해 지역 변동성을 반영하지 못하지만, GIGJ는 중·하부 지각 두께를 현장 데이터에 맞게 재조정한다. 구체적으로, 중간 지각 두께가 평균보다 얇아져 지오뉴트리노 기여가 21 % 감소하고, 하부 지각 두께가 증가해 기여가 24 % 상승한다. 이러한 재배치는 지오뉴트리노 신호의 주요 불확실성 원인인 지각 내 U·Th 분포에 직접적인 영향을 미친다. 모델 적용 후, 상·중·하부 지각 각각에 대한 U·Th 함량 가정 시 신호 불확실성은 77 %, 55 %, 78 %로 크게 감소하였다. 이는 JUNO 실험이 지오뉴트리노를 이용해 맨틀 방사능 기여를 정밀히 분리할 수 있는 기반을 제공한다.

또한, 모델은 공개 데이터베이스(http://www.fe.infn.it/u/radioactivity/GIGJ)로 제공돼, 향후 다른 지오뉴트리노 실험이나 지역 지구물리학 연구에 재사용이 가능하다. 한계점으로는 퇴적층을 고정된 두께로 가정한 점과, 위성 중력 데이터의 공간 해상도가 30–50 km 수준이므로 0.1 km 세로 해상도는 지진 자료에 크게 의존한다는 점을 들 수 있다. 향후 지진 관측망 확대와 고해상도 중력 위성(예: GRACE‑FO) 데이터 통합을 통해 모델의 정밀도를 더욱 향상시킬 여지가 있다.