대기층에서의 난류 확산과 난류 열 확산에 의한 에어로졸 이동

초록

본 논문은 난류 열 확산(TTD) 현상을 이론적으로 재조명하고, 실험 및 위성 관측을 통해 그 존재를 확인한다. TTD는 온도 구배 방향으로 비확산성 에어로졸 플럭스를 발생시켜 온도 역전층 근처에 장기간 지속되는 에어로졸 층을 형성한다. GOMOS 위성 자료와의 비교를 통해 열역학적 프로파일과 일치하는 에어로졸 수직 분포를 설명하고, TTD와 난류 확산 계수의 비율을 이용해 대류권 상부의 난류 확산 계수를 독립적으로 추정한다. 또한 수치 확산 모델에 TTD를 포함시켜 하층 대류권에서의 에어로졸 상승 흐름을 재현했으며, 입자 크기가 클수록 TTD 효과가 크게 나타남을 확인하였다.

상세 분석

본 연구는 난류 열 확산(Turbulent Thermal Diffusion, TTD)이 대기 중 에어로졸 분포에 미치는 영향을 다각도로 검증한다. 먼저 기존 이론에서 제시된 에어로졸 플럭스 식 J = ‑K∇C ‑ K_TTD (C/T)∇T 를 재정립하고, 여기서 K는 전통적인 난류 확산 계수, K_TTD는 TTD에 대응하는 계수이다. 두 계수의 비율 α = K_TTD/K는 난류 강도와 온도 구배에 따라 달라지는 함수로, 실험실 풍동 실험에서 입자 크기와 레이놀즈 수에 따른 α 값을 측정함으로써 이론적 예측과 일치함을 확인하였다.

위성 관측 측면에서는 GOMOS(글로벌 오존 모니터링 by 스타 옥클루전) 자료를 활용하였다. GOMOS는 별의 광흡수를 통해 대기 중 에어로졸 소실률과 온도 프로파일을 동시에 제공한다. 논문에서는 열역학적 역전층 근처에서 에어로졸 농도가 온도 최소점 양쪽에 대칭적으로 높게 나타나는 현상을 발견했으며, 이는 TTD가 온도 구배 방향으로 에어로졸을 끌어당겨 역전층에 축적시키는 메커니즘과 일치한다.

특히 중요한 점은 α와 K 사이의 관계식을 이용해 관측된 에어로졸 프로파일로부터 K 값을 역산할 수 있다는 것이다. 기존에는 라디오소닉 측정이나 풍동 실험을 통해서만 추정되던 대류권 상부의 난류 확산 계수가, 이번 연구에서는 위성 기반 독립적인 방법으로 산출되었다.

수치 모델링에서는 3차원 라그랑지안 입자 확산 모델(FLEXPART 등)에 TTD 항을 추가하였다. 모델 실험에서는 입자 직경 0.1 µm에서 2 µm까지의 크기 구간을 고려했으며, 결과는 하층 대류권(0–3 km)에서 평균적으로 10–30 % 정도의 상승 플럭스를 유발한다. 특히 1 µm 이상 입자는 중력 침강이 강하지만, TTD에 의해 상승 흐름이 보강되어 실제 수직 분포가 평탄해지는 효과가 관측되었다.

이러한 결과는 기후 모델에서 에어로졸의 방사능 효과와 구름 응결핵(CCN) 공급을 보다 정확히 평가하기 위해 TTD를 포함시켜야 함을 시사한다. 또한, 온도 역전이 빈번히 발생하는 대기 경계층이나 고산 지역에서도 유사한 에어로졸 축적 현상이 기대되며, 이는 지역 기후와 대기 화학에 중요한 영향을 미칠 수 있다.