[논문 리뷰] Modeling of path delay in the neutral atmosphere: a paradigm shift

Leonid Petrov|arXiv (Cornell University)|2015. 02. 24.

GNSS positioning and interference참고 문헌 3인용 수 3

한 줄 요약

이 논문은 항성 천문학을 위한 대기 경로 지연 모델링 방식의 전환을 제안한다. 고정밀 수치 기상 모델 출력물(GEOS 및 MERRA)을 활용하여 45 ps × cosec(기울기)의 정확도로 경로 지연과 흡수를 사전에 계산한다. 이 방법은 별 위치 갱신 수준과 동등한 빈도로 일상적이고 실시간으로 전파 효과를 계산할 수 있게 하며, 총 경로 지연의 경우 0.5% 정확도, 수증기 성분의 경우 10% 정확도를 달성하여 전통적인 보정 방법(예: 기울기 곡선)을 더 이상 필요로 하지 않는다.

ABSTRACT

Computation of propagation effects in the neutral atmosphere, namely path delay, extinction, and bending angle is a trivial task provided the 4D state of the atmosphere is known. Unfortunately, the mixing ratio of water vapor is highly variable and it cannot be deduced from surface measurements. That fact led to a paradigm that considers path delay and extinction in the atmosphere as a~priori unknown quantities that have to be evaluated from the radio astronomy data themselves. Development of our ability to model the atmosphere and to digest humongous outputs of these models that took place over the course of the 21st century changed the game. Using the publicly available output of operational numerical weather model GEOS run by NASA, we are in a position to compute path delay through the neutral atmosphere for any station and for any epoch from 1979 through now with accuracy of 45 ps * cosec elevation. We are in a position to compute extinction with accuracy better than 10 pro cents. We are in a position to do it routinely, in a similar way how we update apparent star positions for precession and nutation. Moreover, we are in a position to do it now. As a demonstration of current capabilities, I have computed time series of path delays for all radiotelecopes that I was aware of (220 sites) since 1979 with a step 3-6 hours. Results of the validation tests are presented. A new paradigm of data analysis assumes that we know the atmosphere propagation effects a priori with the accuracy higher that one could deduce them from radio astronomy observations.

연구 동기 및 목표

• 기존의 데이터 기반 대기 경로 지연 모델링을 수치 기상 모델 출력물을 활용한 물리 기반 사전 계산으로 대체하기 위해.
• 1979년 이후의 모든 관측소와 시점에 대해 경로 지연과 흡수를 일상적이고 정확하게 계산할 수 있음을 입증하기 위해.
• 기본선 길이 반복성과 천정 지연 잔차를 기반으로 사전 계산된 경로 지연 정확도를 검증하기 위해.
• 모델 기반 보정이 훨씬 뛰어난 정확도를 보이므로, 전통적인 보정 기법(예: 기울기 곡선 및 지질학적 블록)이 더 이상 필요로 하지 않음을 보여주기 위해.
• 예측 및 세차 운동 갱신 수준과 동등한 실시간 고정밀 대기 보정을 가능하게 하기 위해.

제안 방법

• 1979년 이래로 공개된 나사의 GEOS 및 MERRA 수치 기상 모델의 4차원 대기 상태 출력물(온도, 기압, 수증기 부분 기압)을 활용한다.
• 페르마 원리를 적용하고 굴절도 기울기에서 유도된 4차 비선형 미분방정식 시스템을 사용하여 광선 경로의 변분 문제를 해결한다.
• 높이, 경도, 위도, 시간에 대해 연속적이고 미분 가능한 굴절도 모델링을 가능하게 하기 위해 B-spline 함수의 4차원 텐서 곱으로 굴절도를 표현한다(차수 m).
• B-spline 전개를 사용하여 애자이머-기울기 격자에 걸쳐 경로 지연 적분을 계산하고, 효율적인 검색 및 적용을 위해 계수를 저장한다.
• 기본선 길이 반복성 테스트를 통해 결과를 검증하며, 가중 평균 제곱편차를 사용하고 지구의 지름으로 외삽하여 수직 위치 오차를 추정한다.
• 프리지 진폭을 캘리브레이션하기 위해 e^(-a) 요소를 사용하고, 모델 기반 대기 기여도를 이용해 시스템 온도에서 수신기 온도를 분리함으로써 대기 밝기 온도와 감쇠를 계산한다.

실험 결과

연구 질문

• RQ1 기존의 데이터 기반 보정을 대체할 수 있을 정도로 충분한 정확도로 사전에 대기 경로 지연을 계산할 수 있는가?
• RQ2 현대 수치 기상 모델을 활용할 경우 사전 경로 지연 및 흡수 모델링의 달성 가능한 정확도는 어느 정도인가?
• RQ3 사전 모델링의 정확도는 기울기 곡선 및 지질학적 블록과 같은 전통적 방법과 비교하여 어떻게 되는가?
• RQ4 모델 기반 밝기 온도를 활용해 수신기 및 대기 기여도를 분리함으로써 시스템 온도 이상치를 어느 정도 보정할 수 있는가?
• RQ5 새로운 패러다임은 예측 및 세차 운동 갱신과 동등한 방식으로 실시간, 일상적인 전파 효과 보정을 가능하게 하는가?

주요 결과

• 경로 지연는 45 ps × cosec(기울기)의 정확도로 계산 가능하며, 총 경로 지연의 경우 약 0.5%의 상대 정확도를 달성한다.
• 수증기 성분의 경로 지연은 10%의 정확도로 계산되며, 평균적으로 천정 수증기 지연에서 37 ps의 불확도를 가진다.
• 지구의 지름으로 외삽한 기본선 길이 반복성에서, 사전 경로 지연 모델의 오차로 인해 수직 관측소 위치에 추가로 1.45 cm(48 ps)의 분산이 발생한다.
• 천정 수증기 지연 잔차는 극지역 관측소에서 0.6 cm(37 ps), 열대 지역 관측소에서 2.2 cm(66 ps)이며, 평균은 1.25 cm(37 ps)이다.
• 대기 밝기 온도와 감쇠는 수증기 지배 주파수 대역에서 수증기 지연의 불확도와 유사한 불확도를 가지며, 이는 약 10%의 정확도를 의미한다.
• 기울기 곡선 및 지질학적 블록과 같은 전통적 보정 기법은 더 이상 필요로 하지 않으며, 모델 기반 보정이 이들의 정확도를 초월하기 때문이다.