Aims & Scope

본 연구에서 약 78 m 코어에 대해 전기 전도를 측정하였다(그림 3). 전기 전도도는 최소값 0.4 μA, 최대값 7 μA, 평균값 1.3 μA을 나타냈으며, 눈에 띄게 높은 5개의 피크를 확인 할 수 있었다. 지난 203년간의 역사적으로 기록된 대표적인 화산활동으로는 Pinatubo (1991년), Agung (1963년), Tarawera (1886년), Krakatau (1883년), Coseguina (1835년), Galungung (1822년), Tambora (1815년) 등이 있다. 전기 전도도로 나타나는 피크는 염기성 성분들의 유입량이 증가할 경우 낮아질 수 있으므로(Moore et al., 1994), 남극에서 나타나는 화산폭발지수(Volcanic Explosivity Index, VEI) 5이상의 화산들로 선정하였다.

Fig. 3. 
Electrical conductivity versus depth. Uncertainty (standard deviation) is shown in shaded gray line. More than 6 volcanic signals were identified.

1815년 인도네시아 Lesser Sunda Island에 위치한 Tambora 화산폭발은 화산폭발지수 7에 해당하며 역사적으로 기록된 화산 중 2번째로 큰 화산폭발로써 “그 해 여름이 없었다”라는 역사적인 기록이 남아 있을 정도로 전 지구적으로 영향을 미친 화산폭발이다. 지금까지 시추 된 대부분의 빙하코어에서 그 기록이 나타나고 있어 빙하코어의 연대측정에 널리 활용되어 왔다. Tambora와 같은 대폭발은 많은 양의 H₂SO₄를 방출하며 대기에서 구름 응결핵(cloud condensation nuclei, CCN)으로 작용하여 태양에너지를 반사하여 전지구적인 온도하강을 가져오는 등 기후변화와 밀접하게 연관된다(Zielinski et al., 1994). Tambora 화산폭발은 적어도 4년 동안 강설의 산성도를 증가시켰다고 보고된 바 있으며(Simkin and Siebert, 1994), 남극에서는 1816년부터 침적되어 나타나는데(Dai et al., 1991), 알려지지 않은 1809년의 큰 화산활동(VEI=6)에 뒤이어 관찰된다고 보고 되었다(Yalcin et al., 2006). 본 연구에서는 펀 압축 모델로 1823년(74 m)에 해당하는 지점에 약 5.4 μA의 전기 전도도 피크가 나타난다. 그러나 전지구적으로 영향을 미친 화산활동처럼 피크 이 후의 산성도 증가가 나타나지 않으며, 전기 전도도 값이 최대값을 갖지 않아 Tambora 화산에 해당하는 피크라고 보기 어렵다(그림 4e). 반면, 펀 압축 모델로 1839년(69 m)으로 추정된 구간의 피크는 가장 높은 7 μA의 전기 전도도 값과 피크 이후에 약 1.5년 동안의 뚜렷한 산성도 증가를 나타내므로 Tambora 화산 피크로 보기에 적합하다고 판단되며(그림 4d), 펀 압축 모델로 계산된 연도와 약 23년의 오차를 보인다. 1823년(74 m)에 해당하는 지점에 나타난 피크는 Tambora 화산 피크에서 약 5 m정도 떨어져 있으며, 연간적설량 기준으로 약 20년 전인 1796년에 해당하여, 이는 Tambora 화산이 동반하는 1809년에 알려지지 않은 화산활동으로 보기 어렵다.

Fig. 4. 
Volcanic peaks identified in Fig. 3.

1822년에 인도네시아에 위치한 Galunggung화산은 VEI 5에 해당하는 폭발을 하였다. 이는 동남극에 위치한 East Dronning Maud Land 빙하코어의 전기 전도도에서도 나타난다(Moore et al., 1991). 이 외에 13년뒤인 1835년에 니카라과에 위치한 Coseguina 화산이 VEI 5에 해당하는 폭발을 하여 전기 전도도로 나타난다. 본 연구에서 Galungggung 화산 피크와 Coseguina 화산 피크는 펀 압축 모델로 각각 1845년(67 m), 1855년(64 m)으로 서로 약 10년 정도의 차이를 두고 피크가 나타나며, 두 화산활동 모두 전기 전도도 값은 4.58 μA이다(그림 4d).

본 연구에서 57 m 구간에 4번째로 높은 전기 전도도 값(5.4 μA)이 나타난다(그림 4c). 펀 압축 모델로 계산된 연대는 1874년에 해당한다. Tambora와 Krakatau 사이에 큰 화산폭발로는 러시아 캄차카반도에 위치한 Shiveluch 화산이 VEI 5에 해당하는 폭발을 하였다. 그러나 북반구 고위도에 위치한 Shiveluch화산은 연구지역에서부터 약 14,000 km 떨어져 있어 남극에서 높은 산성도를 초래하기에는 다소 무리가 있다고 판단된다.

1883년 8월 인도네시아에 위치한 Krakatau화산의 화산폭발은 VEI 6에 해당하는 큰 폭발이다. 이는 몇몇의 다른 남극의 빙하코어에서도 나타난다(Langway et al., 1988). Krakatau 화산 폭발 이 후 약 3년 뒤인 1886년 6월에 뉴질랜드의 Tarawera 화산이 VEI 5에 해당하는 폭발을 하였다. 위에서 확인된 Tambora 화산의 폭발 연대와 펀 압축 모델의 차이를 바탕으로 3년 간격으로 나타난 전기 전도도 피크를 확인하였다. Krakatau와 Tawawera는 각각 49 m와 48 m지점에서 5.4 μA와 4.7 μA의 전기 전도도 값을 나타냈다(그림 4b). 펀 압축 모델 연대로 각각 1899년과 1902년으로 펀 압축 모델에서 계산된 연대와 약 15년의 차이가 난다.

10-12 m 구간에 약 4 μA 피크 3개가 나타난다(그림 4a). Pinatubo화산은 1991년에 폭발한 대표적인 화산으로서, VEI 6에 해당한다. 여러 개의 피크가 나타나는 원인으로는 같은 기간에 칠레에 위치한 Cerro Hudson화산이 VEI 5에 해당하는 폭발을 하여 영향을 주었을 것으로 보여진다. 같은 기간에 폭발한 화산이 서로 다른 피크로 나타나는 원인으로는 연구지역과 화산폭발 지역 간의 거리와 화산폭발 지수가 서로 다르기 때문이다. 이 밖에 다른 원인으로는 1991- 1993년의 엘리뇨의 확장으로 MSA의 유입량 증가에 따른 원인이 있으나, 이는 이온 성분 분석을 통해 명확하게 규명할 필요가 있다(Karlöf et al., 2000). 펀 압축 모델로 계산된 연대는 1992-1995년으로 큰 오차를 보이지 않았다.

화산활동에 의한 영향이 폭발 1년뒤에 남극에 나타난다고 가정하고 깊이-연대 모델과 비교해보았다(그림 5). 인도네시아에 위치한 Agung 화산의 경우 1963년에 VEI 5에 해당하는 폭발을 하였다. 약 27.5 m지점에서 약 4.1 μA의 피크가 나타난다. 그러나 펀 압축 모델에서 계산된 연대에서 약 -7년의 차이를 보이며, 그림 5와 같이 Agung은 깊이-연대 모델에서 벗어나 있으므로 Agung 화산 피크로 보기 어렵다. 추후 Agung 화산 피크와 알려지지 않은 피크들에 대해서는 물안정동위원소 조성비 분석과 이온성분 분석을 통해 보다 정밀하게 연대측정을 할 필요성이 있다. 또한 펀 압축 모델로 계산된 깊이-연대와 화산활동 피크가 나타난 깊이-연대는 상부 구간에서 일치하나 하부 구간으로 갈수록 그 오차가 커졌다. 이러한 오차를 초래할 수 있는 이유 중 하나는 시간에 따른 적설량 변화이다. Herron-Langway 펀 밀도 모델에서는 전기간 동안 일정한 연간적설량을 가정하지만 실제로는 적설량은 변화를 보인다.

Fig. 5. 
Depth-age relationships built for different accumulation rates. The accumulation rates were estimated from the Herron-Langway firn densification model, previous study and volcanic peaks.

본 연구에서 측정한 전기 전도도는 대부분 일정한 주기성을 나타내었다. 빙하코어의 배경 전기 전도도 특성을 조사하기 위하여 화산활동의 영향이 없는 깊이 구간(20~25 m)을 선택하여 물안정동위원소 조성비와 비교하여 계절별 특성을 알아보았다. 전기 전도도 값의 변화는 물안정동위원소 조성비의 변화와 비슷한 거동을 보인다(그림 6). 봄-여름철에 2~4 μA 수준까지 증가 하며, 겨울철은 0.5~1 μA 수준으로 감소하여, 겨울철 값들에 비해 봄-여름철이 약 2~8배 높은 양상을 보인다. 그러나 모든 구간에서 물안정동위원소 조성비와 일치하지 않았다. 대부분의 피크들은 봄-여름 또는 여름-가을 시기에 전기 전도도가 증가하였으며 24~25 m 구간에서는 여름철 피크가 나타나지 않기도 하였다.

Fig. 6. 
Seasonal variation in conductivity observed by comparison with δ¹⁸O.

빙하코어 내의 산성도는 강설로 유입된 산성 성분과 염기성 성분의 양에 의해 결정되므로 전기 전도도의 계절 변화 특성을 이해하기 위해서는 빙하 시료에 유입된 산성 및 염기성 성분들의 발생 기작과 수송과정에 대한 이해가 필요하다.

H₂SO₄은 남극대륙 주변의 해빙이 소멸하는 봄철과 여름철 시기에 남극해에 번성하는 식물성 플랑크톤에 의해 대기 중으로 방출되는 DMS의 대기산화반응을 거쳐 생산된다(Kettle and Andreae, 2000). HCl은 H₂SO₄와 해염기원 미세입자에 존재하는 NaCl 성분과의 화학반응에 의해 생성되므로 H₂SO₄의 농도가 증가하는 봄철과 여름철에 농도가 증가한다. HNO₃의 생성원과 수송과정에 대한 이해는 H₂SO₄와 HCl에 비해 부족한 편이다. Legrand and Delmas (1986)은 남극 대기환경의 HNO₃은 성층권 응결핵의 강하, 성층권과 대류권간의 대기 교환, 저위도와 중위도 지역에서 발생한 NOX의 장거리 수송 과정 등에 의해 영향을 받는다고 하였다. 또한, NO₃의 계절별 유입량의 변화에 따라 HNO₃의 농도도 변화한다고 보고된 바 있다(Taylor et al., 1992).

H⁺의 농도는 산성가스 성분들의 유입에 의해 증가하지만 대기화학 반응에 참여하는 염기성 성분들의 농도가 증가할 경우 산성도는 감소한다. 남극 강설의 NH₄⁺의 농도는 대기 화학 반응을 통해 H⁺ 농도를 감소시키는 대표적인 염기 성분인 NH₃ 성분으로부터 유래되기 때문에 강설의 NH₄⁺ 농도가 증가하는 경우 산성도는 감소한다고 보고된 바 있다(Wolff et al., 1997). 또한, 겨울철 해염기원 성분들(Na⁺, Cl^-)의 농도가 증가하는 경우에는 오히려 산성도에 기인한 고체 얼음 표면의 전기 전도도는 감소하게 된다(Wolff et al., 1997).

빙하코어는 무수히 많은 얼음 결정으로 이루어져 있는 다결정 물질이며, 고체상태(밀도 0.917 Mg m^-3)인 빙하코어는 하나의 얼음결정 주변에 14개에서 16개의 다른 얼음결정으로 꽉 들어차있다(Shabtaie and Bentley, 1994). 반면, 밀도가 감소할수록 주변의 얼음 결정들의 수는 줄어들고 얼음 결정 사이에 공기방울로 채워지게 된다.

본 연구에서 전기 전도도 값은 약 11 m까지는 현저하게 낮은 값을 나타내었다. 전기 전도도와 밀도와의 관계를 확인하기 위해 펀 압축 모델로 계산된 밀도와 전기 전도도 평균값을 계산하여 상관성을 확인해 보았다(그림 7). 계산된 밀도와 측정된 밀도의 오차는 ±0.002이었다. 밀도가 0.55 Mg m^-3에 이르기 까지 밀도와 전기 전도도의 상관계수는 r=0.78으로 나타났다. 밀도와 전기 전도도가 상관관계를 나타내는 원인은 입자경계에 맥 형태로 연결되어 있던 이온성분들이 얼음 결정 사이의 공기방울 또는 가스에 의해 끊기게 되므로 전도가 정상적으로 일어나지 않는 것으로 보인다(Wolff et al., 1988). 그러나 밀도가 0.60 Mg m^-3에 이르기까지의 밀도와 전기 전도도의 상관계수는 r=0.87으로 더 높게 나타났다. 0.55 Mg m^-3보다 0.60 Mg m^-3에서 더 높은 상관관계를 보이는 이유는 임계밀도에서 눈의 결정들이 입자들의 재분배, 입자경계 미끄럼, 재결정화, 이동, 공기방울의 축소와 같은 과정을 거치면서 시추지역마다 임계밀도 구간이 변화하기 때문이다(Hörhold et al., 2011). 밀도가 0.60 Mg m^-3이상 구간에서는 밀도와 전기 전도도의 상관관계는 명확하게 나타나지 않았다. 이를 바탕으로 Herron-Langway 펀 밀도 모델의 임계밀도를 수정하여 0.55 Mg m^-3로 설정한 모델의 결과와 0.60 Mg m^-3로 설정한 모델의 결과를 비교하였다(그림 5). 0.60 Mg m^-3로 설정한 모델의 연간적설량은 0.23 Mg m^-2y^-1으로 0.55 Mg m^-3로 설정한 모델의 연간적설량(0.23 Mg m^-2y^-1)보다 적게 나타났다. 그림 5와 같이 0.55 Mg m^-3로 설정한 Herron-Langway 모델에서 Pinatubo, Tawawera, Krakatau, Coseguina, Galungung, Tambora 화산활동은 각각 1992(0)년, 1887(+15)년, 1884(+15)년, 1836(+19), 1822(+22) 1816(+23)년의 차이를 보인다. 반면, 0.60 Mg m^-3로 설정한 Herron-Langway 모델에서는 각각 1992(-3)년, 1887(-3)년, 1884(-2)년, 1836(-5)년, 1823(-3)년, 1816(-2)년으로 하부로 갈수록 심해지는 오차가 줄어 들었으며, 화산활동의 폭발 시기가 0.60 Mg m^-3로 설정한 Herron-Langway 모델의 연대 결과와 유사하게 나타났다.

Fig. 7. 
The relationship between conductivity and density above the critical density. The solid and dotted line indicate the linear regression lines for ρ < 0.55 Mg m^-3 and ρ < 0.60 Mg m^-3, respectively.