1. 서론
기후변화로 유발되는 평균 해수면 변화 양상은 광역에서 장기적으로 점진적으로 발생하는 현상이지만, 그 정도나 특성은 국지적이고 일시적인 영향을 받을 수 있다. 어떤 특정 지역의 평균 해수면 변화 양상은 가용한 장기 자료를 이용하는 것이 적절할 수 있으나, 우리나라 모든 연안 해역의 해수면 변화 양상 분석은 가용한 자료 각각의 분석과 더불어 상호간의 상관관계 분석도 요구된다. 특히 평균 해수면 자료는 관측기간(자료의 개수), 변화 양상, 결측(missing) 양상 및 정도 등이 지점에 따라 다르기 때문에 가용 자료의 특성, 국지적인 특성 및 기본 통계정보를 고려하는 분석이 요구된다. Jung and Jeong(2013)과 Harker et al.(2019)는 서해 연안해역과 호주 주위의 해수면 변화를 추정하기 위해 조화분석을 통한 해수면 변화 분석을 수행했다. 한편, 우리나라의 평균해수면 분석에 관한 연구는 다양한 방법 및 관점에서 수행됐다. Hwang et al.(2016)은 인공위성 자료를 활용한 해수면 변동에 관한 연구를 수행했으며, Heo et al.(2018)은 IPCC AR5에서 사용된 CMIP5 모델 자료를 사용해 우리나라 주변 해역을 6개 구역으로 분리하여 해수면 상승의 특징을 분석했다. 한편, Yoon and Kim(2012)은 방재 관점에서 검조소 17개 지점의 평균해수면 변화 정도를 분석하였으며, Jung(2014), Lim et al.(2017)은 기후변화에 따른 영향분석의 관점에서 분석을 수행했다. 특히 Jung(2014)은 서해 연안의 평균해수면 변화 양상을 다양한 측면에서 분석하였다. 또한, Kim et al. (2016)은 특정 모델을 적용하여 극치해수면 변화양상(비정상성) 분석을 수행하였다. 대부분의 평균해수면 변화 분석은 분석 내용이 추세 분석으로 제한되어 있으며, 사용하는 자료의 기간, 지점, 이상자료 처리과정을 포함한 분석에 사용한 자료의 범위 등은 분석시점에 따라 서로 상이하고, 명확한 처리과정이 생략되어 있다. Wahl et al.(2010)은 평균 해수면 시계열을 생성하고 분석함에 있어 SSA(Singular System Analysis)와 같은 비모수적 데이터 분석 및 모수적 Fitting 접근방법을 통해 관측된 해수면 상승을 분석했으며, Cho et al.(2011)은 우리나라 조위자료의 결측 양상을 정리-분석하여 제시한 바 있다. MSL 변화는 특정 지역에 국한되는 국지적인(local) 변화가 아니라 광역적인(global) 변화이기 때문에 우리나라 전체의 MSL 자료를 이용한 변화 양상 분석이 필요하며, 광역적인 변화양상이 반영되어야 한다는 가정에 대한 통계적인 확인이 필요하다.
본 연구에서는 국립해양조사원에서 제공하는 모든 평균해수면 자료를 이용하여 자료의 통계적인 특성 및 상관관계를 분석하였다. 또한 추정 상관계수의 신뢰구간 분석, 계수 추정에 영향을 미치는 이상자료(outliers) 진단 및 평균해수면 자료의 품질에 대한 확률적인 진단도 수행하였다. 본 연구결과로 제시되는 공간적인 상관관계 분석 결과는 평균해수면 정보가 단기간으로 제한되어 있는 지점의 최적 평균해수면 추정에 활용할 수 있으며, 신뢰수준 향상을 위한 자료의 품질 관리 및 결측구간의 자료보충(filling-in) 영향 분석에도 적용 가능하다. 어떤 특정 지점, 특정 연도의 자료가 전체적인 변화 양상에서 크게 벗어나는 경우는 자료의 QC 관점에서 이상자료로 판단할 수 있으며, 그 이상자료에 대한 적절한 원인분석을 수행하여 자료의 사용여부를 판단할 수 있다. 한편 결측구간의 자료보충은 자료의 변동양상을 파악하여 최적의 방법을 찾는 과정이 필요하지만, 간단한 방법으로는 지점 또는 연도에 따른 상관관계 정보 등을 이용하여 신뢰할만한 인접 지점, 인접 기간의 자료로 보충할 수 있다.
2. 자료 및 방법
본 연구에서 사용한 평균해수면 자료는 국립해양조사원에서 공개되는 자료로, 누구나 접근하여 받을 수 있는 공용자료로 분석 결과의 상호검정을 위한 공유 자료로서 큰 의미를 갖고 있다. 평균해수면 자료는 국립해양조사원(www.khoa.go.kr) -> 해양정보포털(실시간해양관측시스템, http://www.khoa.go.kr/koofs/kor/tide/tide.do -> 조석예보 및 성과 -> 조석성과 -> 평균해면 성과표 자료이며, 이 자료는 월평균 해수면 자료와 연평균 해수면 자료, 조위관측소 위치정보를 포함하고 있다. 본 연구에서 사용한 평균 해수면 자료가 제공되는 조위관측소 위치는 Fig. 1에 도기하였으며, 명칭 및 위치, 자료 개수 등은 Table 1과 같다. 우리나라에서는 총 50개 지점에서 평균 해수면 자료를 제공하고 있으나, 순천만, 포항 지점은 관측기간이 매우 짧아 MSL 자료가 없기 때문에 분석에서 제외하였다. 이 외에도 경인항, 삼천포 등과 같이 자료 개수도 적고, 결측자료가 많은 관측지점도 존재한다. 각각의 지점에 관측한 MSL 변동특성을 고려하여 분석하는 과정이 매우 중요하지만, 특정 지점에 국한되거나 명확한 변동 특성 파악이 필요한 경우(지반침하 및 지진으로 인한 일시적인 변화 등)가 발생할 수 있다. 그러나 본 연구에서는 이러한 국지적인 분석보다는 MSL 자료의 시간· 공간적인 전체적인 변화 양상 분석에 중점을 두고 수행하였다.
본 연구에서 사용한 해수면 자료는 연 평균 해수면 자료이며, 연간 결측 기간이 상당한 경우에는 부분적으로 월 평균 해수면 자료가 가용함에도 불구하고, Table 1에 정리된 바와 같이 연 평균 해수면 추산자료를 제공하지 않는 경우가 빈번하다. 평균 해수면 자료의 가용여부를 결측인자 행렬(missing indicator matrix)을 도시하였으며(Fig. 2 참고), 분석 시간구간은 목포 검조소의 1959년을 시점, 2018년을 종점으로 하고 있다. 조위관측지점마다 관측 시점이 다르고, 가용한 기간이 달라 자료 개수가 차이가 나기 때문에 신뢰수준에 차이가 발생할 수 있다. 상관계수 분석은 두 지점의 MSL 자료가 모두 존재하는 시기의 자료가 존재하는 경우에만 계산할 수 있으며, 자료개수에 의한 영향은 추정 상관계수의 신뢰구간에 반영된다. 또한 어떤 지점의 MSL 추세 분석도 자료의 개수 영향을 받기 때문에 다른 지점의 추세와 동등한 수준으로 비교하는 것은 곤란할 수 있으나, 실질적으로 불가피한 제약사항이다. 그림에서 볼 수 있는 바와 같이, 조위 관측소는 지속적인 증가양상을 보이고 있으며, 2010년 이후 결측이 급격하게 감소하고 있다. 평균 해수면 추세 분석 등을 위한 가용 자료의 개수가 30개 이상인 지점은 전체 48개 지점에서 총 18개 지점으로 38% 정도이다(Fig. 3 참조).
한편 조위 관측 지점간의 MSL 자료 상관관계 분석을 위한 조합의 개수는 nC2 = 49C2 = 1,176 가지의 경우가 발생한다. 그 경우에 지점 상관분석을 위한 동일 연도의 자료가 12개 이상 존재하는 경우로 한정하여 계산을 수행하였으며, 이 경우의 상관계수 추정이 가능한 조합은 304개, 비율로는 26% 정도이다. 상관분석을 위한 자료개수를 30개 이상으로 조정하는 경우 신뢰구간은 감소하지만, 상관관계 분석을 위한 조합은 175개, 비율로는 9% 수준으로 감소하기 때문에 분석에 상당한 제약이 따른다. 따라서 본 연구에서는 상관분석에 사용한 두 지점의 공동 기간 자료개수를 12개로 선정하여 상관분석을 수행하였다. 최소 자료의 개수를 12개로 제한한 이유는 추정 상관계수의 신뢰구간(오차범위)을 적절한 범위로 제한하고, 작은 표본(small sample)의 자료 개수 추천기준을 반영한 것이다. 신뢰구간 또는 오차구간은 이론적으로 자료의 표준편차가 일정한 경우, 추정에 사용되는 자료의 개수(n)가 증가하면 1 / n
3. 결과 및 토론
3.1 상관계수 분석
두 지점 사이의 상관 계수 추정이 가능한 조합에서 자료 개수 12개 이상 조건을 만족하는 304개 경우에 대한 상관계수는 다음 그림과 같다(Fig. 4 참조). 그림에서 볼 수 있는 바와 같이 약 10개(3% 정도)는 음의 상관을 보이고 있으나, 전체적으로는 0.5 이상의 수치를 보이고 있다. 상관계수 자료의 빈도분포를 보면 상관계수 0.5 이상이 75% 정도를 차지하고 있다(Fig. 5 참조). 추정 상관계수의 신뢰구간은 변동성분과 자료개수의 영향을 받아 조합에 따라 큰 차이를 보이고 있음을 알 수 있다(Fig. 6 참고).
상관계수의 신뢰구간 추정에 영향을 미치는 또 하나의 중요한 인자는 자료에 따라 차이를 보이지만, 자료에 포함되어 있는 소수의 이상자료이다. 이 이상자료는 MSL 자료의 품질 검정 과정에서 면밀한 검토와 적절한 처리가 요구된다. 본 연구에서는 Rosner 방법을 이용하여 이상자료를 진단하여 분석하였으며, 그 이상자료가 특정 지점의 자료에 국한되거나 특정시기로 제한(대산, 2003; 진도, 2006; 여수, 1980; 도농탄, 2014, 2015; 흑산도, 1968 등)되어 있는 것으로 파악되었다. 따라서 전체적인 MSL 분석 및 본 논문의 핵심주제에 해당하는 공간 상관계수 분석에서는 이상자료가 차지하는 비율이 매우 낮아 그 영향은 미미한 수준이다.
아래 그림(Figs. 7~8)은 자료 개수에 따른 신뢰구간과 상관 계수 변화를 도시한 그림이다. 그림(Fig. 7)에서 볼 수 있는 바와 같이 자료의 개수 증가에 따라 신뢰구간의 크기가 감소(높은 신뢰수준)하는 일반적인 경향을 보이고 있으나, 자료의 개수가 20개 이하인 경우에서는 크고 작은 신뢰구간의 크기가 발생하고 있는 것으로 보아, 소수 자료의 변화 양상에 크게 영향을 받는 것으로 판단된다. 한편 다음 그림(Fig. 8)에서 보이는 바와 같이, 12개 이하의 작은 개수의 자료에서 낮은 상관계수를 보이거나 30개 이상의 다수 자료에서 낮은 상관계수를 보이는 결과는 특정 자료(이상자료)의 영향으로 파악되었다. 이러한 영향은 상관분석에 사용한 자료의 원시자료의 검토 및 적절한 처리, 품질 관리가 요구된다.
3.2 공간 상관관계 분석
지점간의 상관계수 및 신뢰분석에 더불어 지점간의 거리를 고려한 공간 상관관계를 분석하였다. 해수면 변화가 광역적인 현상이라면 거리에 따른 공간 상관계수의 뚜렷한 변화를 기대할 수 없다. 이러한 가정에서 보면, 우리나라에서의 거리에 따른 MSL 자료의 상관계수는 거리에 무관한 광역적인 형상으로 파악된다. 이는 MSL 장기 변화가 국지적으로는 일시적인 해안개발 영향 및 지역적인 특성을 반영할 수 도 있으나, 전반적으로는 거리와 무관하게 그리고 해역 특성과도 무관하게 0.8 정도 수준의 높은 상관정도를 보이고 있는 것으로 파악되었으며, 이는 한반도 공간 규모의 광역적인 현상임을 자료로부터 확인할 수 있다(Fig. 9). 그러나 우리나라 연안의 정량적인 MSL 변화는 지역적 특성이 매우 강하게 나타날 수 있다. 역시 소수의 특이한 자료는 적절한 해석보다는 원시자료의 품질 검정이 수행되어야 할 것으로 판단된다. 특정 지점의 정량적인 MSL 변화 양상은 공간적으로 정도나 양상의 차이를 보일 수 가 있으나, 본 연구에서는 정성적인 분석으로 제한하였으며, 선형 관계로 분석이 제한되는 한계를 가지고 있다.
3.3 MSL 자료의 변화양상 분석
48개 검조소 자료의 평균해면 변화 양상을 분석한 결과는 Fig. 10에 도시하였다. 일반적으로 평균해수면의 변화 양상은 선형 추세변화를 가정하여 수행하고 있으며, 해수면 변화 정도는 연간 mm 단위로 제시하는 경우가 일반적이다. 그러나 우리나라 모든 검조소의 평균해수면 변화양상을 도시한 결과를 살펴보면, Jung(2014)도 이미 보고한 바와 같이 선형 추세 가정을 벗어나는 경우(21, 48 - 서귀포, 포함) 또는 변동이 크게 나타나는 자료가 빈번하게 발생하고 있음을 알 수 있다. 따라서 조위관측소간의 평균해수면 변화 상관관계 분석을 비선형의 경우로 확장하는 연구가 필요하며, 비모수적인 방법을 이용한 방법도 수행할 필요가 있을 것으로 판단된다.
4. 결론 및 제언
우리나라 조위 관측소의 연평균 해수면 변화 자료를 이용하여 상관분석을 수행한 결과 다음과 같은 결론이 도출되었다.
(1) 지점별로 상관계수는 0.75 이상이 대부분을 차지하는 높은 수준으로 파악되었다.
(2) 추정 상관계수의 신뢰구간 크기는 전반적으로 자료의 개수 증가에 따라 감소하는 경향을 보이고 있는 것으로 파악되었다.
(3) 두 검조소간의 상관계수는 거리와 무관하게 0.8 정도의 높은 상관정도를 보이고 있는 것으로 파악되었으며, 이는 해수면 변화가 한반도 공간 규모의 광역적인 현상임을 자료로부터 파악이 가능함을 입증하였다.
(4) 조위 관측소의 MSL 변화 양상은 뚜렷한 선형 증가/감소 추세를 보이지 않는 자료가 적지 않은 부분을 차지하고 있기 때문에, 선형 추세 가정에 한계가 있는 것으로 판단된다.
또한 MSL 변화 자료를 이용한 상관계수 추정 및 추세분석 결과는 인접한 조위 관측소와 높은 상관을 보이고 있기 때문에, 다른 인접한 관심 지점의 MSL 변화 자료가 없거나 단기 관측 자료만이 가능한 상황에서도 상관분석 결과를 이용한 해수면 변화 양상의 적절한 추정이 가능할 것으로 판단되나, 가용 자료의 한계로 인한 불확실한 변동성분을 고려한 정량적인 오차분석이 병행되어야 한다.
