Home

도심녹지가 주변 지역의 기온 상승을 얼마나 완화할 수 있는지를 분석하기 위하여 최근 5년간(2011~2015년) 서울에서 가장 무더웠던 날들을 사례로 선정하였다(Table 1). 각 사례일들은 여름철 북태평양 고기압의 영향권에 본격적으로 들면서 덥고 습한 공기로 인하여 서울 지역에 폭염경보(2015년) 및 폭염주의보(2011~2014년)가 발령된 경우들이다. 또한, 빌딩이 밀집해 있고 번화한 지역을 대상으로 폭염 발생 시 도심녹지가 주변 지역 열 환경에 미치는 영향을 분석하기 위하여 서울시 강남구 선정릉 일대를 대상으로 연구를 수행하였다(Fig. 1). 선정릉은 전체 면적 약 20,000 m² 부지 내에 평균 약 15~20 m 높이의 수목과 초지로 이루어진 서울의 도심 속 대표적 녹지 중 한 곳이다.

Fig. 1. 
(Left) The location of the study area around the Seonjeongneung in Seoul, Korea, and (right) the areas and the resolutions of three domains used in the CAS simulation.

CAS는 선정릉과 같이 복잡한 도시 구조를 반영하여 기후영향을 진단 및 평가 할 수 있는 시스템으로써 지표 근처에서 나타나는 물리과정으로 인해 도시 캐노피 안에서 이루어지는 기상현상을 모의하고 그 결과를 수치지도(Digital Map)로 표출한다(Yi et al., 2011). 따라서 다양하고 복잡한 도시공간구조를 고려한 종합적이면서도 공간화된 상세한 도시기후영향 및 진단을 할 수 있는 정보를 제공할 수 있다(Yi et al., 2012a, b; Kwon et al., 2015).

선정릉 주변 열 환경을 모의하는데 있어 CAS는 Fig. 1과 같이 세 가지 모델 영역을 이용한다. 가장 외곽에 위치한 모델 지역(Model Region, MR)은 Yi et al. (2011)과 동일하게 수도권을 포함한 200 × 200 km 영역을 대상으로 공간 해상도 1 × 1 km로 설정하였다. 그러나 CAS 모의 결과의 정확도 향상을 위하여 본 연구에서는 서울 주변지역을 포함하여 50 × 40 km 영역에 대한 연구지역(Study Region, SR) 도메인의 해상도를 기존 25 × 25 m에서 5 × 5m로 개선하여 새로운 SR을 구축하였다. 마지막으로 상세지역(Detail Region, DR) 역시 선정릉을 포함하여 주변 5 × 5 km 영역을 대상으로 SR과 동일한 공간 해상도인 5 × 5 m로 설정하였다.

CAS에 입력되는 공간정보 자료는 수치표고모델(Digital Elevation Model, DEM)과 토지피복(Land Cover, LC)이다. 국립기상과학원에서는 2014년 환경부 토지피복도와 항공우주연구원의 KOMSAT-2 MSS 위성영상, 국토지리정보원의 항공라이다 자료를 이용하여 DEM과 토지피복 자료를 구축하였다(NIMR, 2014). CAS에서 DEM은 자연적인 지형에 대한 정보를 나타내는 수치지형모델(Digital Terrain Model, DTM)과 건물, 수목, 인공구조물을 나타내는 수치표면모델(Digital Surface Model, DSM)으로 세분화된다. 기존 CAS에서 DSM을 생성할 때 DTM과 건물 고도 정보만을 결합한 자료를 이용하였으나, 본 연구에서는 식생 고도의 정보까지 고려할 수 있도록 입력 자료를 개선하여 선정릉 주변에 대해 3가지 입력자료를 구축하였다(Fig. 2). 열 환경을 분석하기 위한 주요 입력자료 중 하나인 토지피복 자료는 KOMSAT-2 MSS 위성영상에 감독분류(Supervised Classification, SC) 기법을 적용하고 국토교통부의 도로명 주소 중 건물자료를 이용하여 Table 2와 같이 여섯 가지 유형으로 분류하였다.

Fig. 2. 
Three spatial information data (DTM, DSM, and LC) of the CAS simulation around the Seonjeongneung region.

이상의 공간 입력정보를 기반으로 CAS 내에 연동된 중규모 기상모형인 MetPhoMod (METeorology and atmospheric PHOtochemistry mesoscale MODel)에 입력할 알베도(albedo), 복사능(emissivity), 체적열용량(volumetric heat capacity), 열확산율(thermal diffusivity), 지표면온도(surface temperature), 지온(ground temperature), 거칠기 길이(roughness length) 등 7가지 parameter 들을 산출하였다. 모델 격자 내의 각 parameter들은 Table 3의 여섯 가지 주요 토지피복 종류들이 차지하는 비율을 퍼센트로 환산하여 fractional coverage (FC)로 계산하고, 각 격자 영역의 parameter들을 식 (1)과 같이 면적비율의 가중치 평균에 의해 구하였다.

여기서 f_LC는 각 격자 내에 여섯 가지 토지피복 형태들이 차지하는 비율로서 범위는 0~1이다. A_LC는 Table 2와 같이 각 토지피복 형태에 대한 속성값들로 Yi et al. (2011)이 독일의 선행연구를 기반으로 재구성한 것이다. 따라서 위 식을 이용하여 도시 공간정보로부터 계산된 열 확산, 열 용량 등 parameter들을 MetPhoMod에 입력하여 대상지역 주변의 기온과 바람장을 모의하고, 최종적으로 모의 결과를 수치지도로 표출 및 생산하였다(Perego, 1999).

Table 4는 선정릉 주변 기상 초기장 및 경계장 생성을 위한 MetPhoMod 모델의 입력값을 나타낸다. 모델 실행을 위해 필요한 초기값은 바람의 U와 V 성분, 가온위, 혼합비 등이 있다. 본 연구는 서울시 종로구 송월동에 위치한 서울기상대(ASOS No.108)의 관측자료를 기준으로 Table 4의 사례기간 중 평균값을 모델 초기값으로 입력하였다.

Figure 4a는 전체 분석사례들에 대해 CAS에서 모의된 선정릉 주변의 찬 공기 생성과 종합 기온편차이다. 찬 공기 생성량은 CAS에서 지면의 피복과 기복에 대한 분석자료와 중규모 기상모델인 MetPhoMod의 모의수행 분석결과를 이용하여 계산한다. 찬 공기생성량은 토지피복 분류 중 녹지를 나타내는 TV (Tall Vegetation)와 VS (Vegetated Surfaces)가 차지하는 비율에 따라 시간 및 면적당 차고 신선한 공기를 생성하는 양으로서 외부 유입에 의한 것이 아닌 녹지에 해당하는 토지피복에 의한 국지냉각(local cooling)을 의미한다. 그림에서 선정릉 내 녹지로 인해 생성되는 차고 신선한 공기가 평균 약 25.5 m³ m^-2 h^-1로 Mosimann et al. (1999)이 숲이나 수목이 많은 지역에 대해 제시한 30.0 m³ m^-2 h^-1와 유사한 값을 나타내었다.

찬 공기 생성량과 더불어 MetPhoMod 모델은 두 개의 공간 규모에 대해 독립적으로 각각 중규모 기온편차(Meso-scale air temperature Deviation, MD)와 국지규모 기온편차(Local-scale air temperature Deviation, LD)를 구한 후 이를 합하여 종합 기온편차(Total air temperature Deviation, TD)를 계산한다(Yi et al., 2011). 중규모 기온편차는 모델 영역 내 공기의 이동과 정체를 모두 계산한 후 야간(전날 20시부터 다음 날 6시)에 고도 10 m에서 산출된 전체 평균기온에서 각 격자별 평균기온을 뺀 차로 구해진다(Kim et al., 2014). Yi et al. (2016)은 현열 플럭스에 의한 지표와 대기 간 열 교환과 난류 혼합에 의한 열 확산, 지표층의 찬 공기 생성에 의한 냉각의 과정들을 아래의 식으로부터 국지규모 기온편차를 계산하였다:

위에서 dT_SHF는 야간에 건물 표면에서 대기로의 현열 방출로 인한 국지규모 기온 증가를 나타내는 것으로 건물이 없는 지역에서는 1의 값을 가지고 건물이 있는 지역에서는 1보다 큰 값을 가진다. 식 (2)에서 dT_SHF는 건물 전체 면적에 대한 종횡비(complete surface aspect ratio), f_CSAR과 현열 플럭스에 대한 empirical coefficient, c_SHF의 곱으로 계산된다(Christen and Vogt, 2004). 식 (3)로부터 f_CSAR은 모델 격자 당 평균 높이 h_B인 단일 건물이 갖는 표면적(지붕 1면과 4면의 builtup surface, f_BS)의 fractional coverage로 정의된다. 건물의 현열 플럭스에 의한 기온 상승과 대조적으로 local cooling, dT_ca는 식 (4)와 같이 찬공기 생성량, Q_ca에 empirical cooling factor, c_ca를 곱하여 구한다. 본 연구에서 c_SHF와 c_ca는 다른 선행연구와 동일하게 각각 2.0 Kh m^-1과 −0.05 Kh m^-1로 설정하였다(Kim et al., 2014; Yi et al., 2016).

공기의 이동과 정체와 관련한 중규모 기온편차와 각 격자별 피복과 건물의 열적 특성 등 국지규모 기온편차를 모두 고려한 종합 기온편차를 Fig. 3b와 같이 정량적으로 산출하였다. 찬 공기 생성량과 동일하게 종합 기온편차는 선정릉을 중심으로 전체 폭염사례에 대해 약 −2.54^oC의 값을 나타내었다. 각 사례일에 대해 산정된 중규모 및 국지, 종합 기온편차를 Table 5에 제시하였다. Table에서 중규모 기온편차는 종관풍에 의한 찬 공기가 주변지역으로의 수송을 고려하기 때문에 Table 4의 기상 입력자료에 따라 변화한다(Yi et al., 2011). 이에 반해 국지규모 기온편차는 토지피복률, 건물 밀도, 건축물의 높이 및 부피에 따라 분석된 상대적인 온도차이로서 모델 격자 별 지리환경정보에 따라 열 환경 정보를 나타내므로 선정릉에 대해서는 −1.2^oC의 동일한 기온편차를 나타내었다. 중규모와 국지규모 기온편차를 합한 종합 기온편차는 −2.0~−2.9^oC로 선정릉의 녹지로 인해 주변 상업지보다 상대적인 낮은 열 환경을 보여준다.

Fig. 3. 
Average distribution of (a) the cold air production (Q_ca) and (b) the total temperature deviation (TD) for the five heat-wave issued cases around the Seonjeongneung region.

CAS는 도시의 공간정보를 기반으로 상대적인 기온 등 열 환경 특성을 분석할 뿐만 아니라 모델 영역 내 공간 구조를 가상으로 변경함으로써 도시기후 변화특성을 정량적으로 모의할 수 있다. 본 연구에서는 선정릉을 고층건물로 대체하여 수목 등 녹지로 인해 주변보다 기온이 낮은 특성이 어떻게 변화하는지를 모의 실험하였다. Figure 4a와 같이 선정릉 영역을 삼성동 코엑스와 동일한 고층 빌딩으로 대체하고 다섯 번의 폭염 사례에 대해 Table 4와 같이 동일한 기상조건으로 CAS를 구동하였다. 그 결과 Fig. 4b와 같이 빌딩들로 인해 차고 신선한 공기 생성이 평균 약 1.1 m³ m^-2 h^-1 정도로 코엑스나 테헤란로 주변과 같이 매우 적은 량을 나타내었다. Table 6에서 각 사례에 대해 동일한 기상조건을 입력하였음에도 불구하고 빌딩의 영향으로 인한 공기의 흐름 및 열 확산 과정 등이 변함에 따라 녹지일 때 보다 중규모 기온편차가 약 50% 높게 나타났다. 국지규모 기온편차 역시 빌딩들의 면적에 따른 현열 플럭스 영향을 반영하여 녹지일 때 보다 3.1^oC 높은 편차를 보였다. 이러한 영향을 모두 고려한 종합 기온편차 역시 0.9~1.5^oC이고 평균 역시 1.14^oC 증가를 보였다(Fig. 4c). 또한 이러한 기온상승 구역이 선정릉뿐만 아니라 남동쪽의 테헤란로와 삼성동 코엑스 부근까지 확장되는 것을 Figs. 3b와 4c의 비교에서 확인할 수 있다. 따라서 선정릉이 녹지 대신 고층 건물로 대체될 경우 여름철 폭염 시 주변의 기온이 평균 3.68^oC 이상 상승하는 등 열 환경이 악화됨을 CAS 모의를 통해 알 수 있다.

Fig. 4. 
(a) LC, (b) Q_ca, and (c) TD of the Seonjeongneung when it is substituted for high-rise buildings.

Figure 5는 전체 사례일에 대해 선정릉과 이 지역이 빌딩으로 대체되었을 경우에 대한 종합 기온편차 차이를 1^oC 간격 등치선으로 나타낸 것이다. 이 그림에서 선정릉 내 격자점에서 산출된 편차값이 주변지역에 비해 4.0~5.0^oC 정도 낮고, 붉은색의 경계 지역도 약 3.0^oC의 낮은 값을 나타내었다. 또한 선정릉 녹지의 영향으로 상대적으로 1.0~2.0^oC 낮은 기온 분포가 선정릉 외곽에 까지 분포하였다. 따라서 선정릉 경계로부터 거리에 따른 종합 기온편차를 Fig. 6과 같이 구한 결과 약 200 m까지 뚜렷한 차이를 확인할 수 있다. 거리별 종합 기온편차는 50 m까지 약 1.6^oC를 나타내지만 선정릉에서 멀어짐에 따라 급속히 감소하여 200 m보다 먼 거리에서는 0.3^oC 이하의 매우 적은 값을 보였다. 따라서 선정릉 내 녹지가 주변 기온하강에 영향을 미치는 거리가 약 200 m 정도로써 Lee et al. (2013)이 선정릉에서 관측한 평균 거리 200~250 m와 잘 일치하였다. 또한 Myeong (2010)이 선정릉보다 10배 가량 규모가 큰 서울의 5개 도심 녹지에 대한 거리별 기온감소 효과 역시 평균 300 m로 제시한 결과와도 일관성을 보였다.

Fig. 5. 
The averaged TD difference distribution between the real land coverage (Fig. 3b) and its imaginary substitution of buildings (Fig. 4c) at the Seonjeongneung for the five heat-wave issued cases. The green type colors show the temperatures deviation and the red line presents the boundary of the urban green area.

Fig. 6. 
The calculated minimum, median, and maximum TD differences shown in Fig. 5 along with the distance from the boundary of the Seonjeongneung.