본 연구는 2020년 3월 18일부터 20일까지 영동지역에 강풍이 발생했던 사례(남고북저형, 대류권계면 접힘에 의해 급격하게 발달하는 저기압)의 종관 및 열역학적/운동학적특성을 조사하기 위해 AWS 관측 자료, 종관 일기도, ECMWF 재분석 자료, 레윈존데, 윈드프로파일러 자료를 이용하였다. 분석결과, 사례 기간 영동지역 5개소에서 관측된 최대 순간 풍속은 20 m s−1 이상으로 나타났으며 대관령(27.7 m s−1)에서 가장 강하게 나타났다. 종관분석에서는 남고북저형의 기압배치와 함께 영동지역으로 등압선의 모양이 사인(sin)파 형태를 보이며 강한 기압경도력에 의해 강풍이 발달 하다가 3월 19일부터는 한반도 북부지역에서 하루 내에 19 hPa 이상의 기압 하강과 함께 발달하는 저기압에 의해 지속적인 강풍이 발달했다. 북강릉 단열선도에서 역전층의 고도는 하층 강풍대와 함께 산 정상의 약 1-3 km 고도에 위치 하였고, 레윈존데 및 수직 측풍 장비(윈드프로파일러)의 연직 바람장 분석 결과와 일치함을 확인할 수 있었다. 특히 열역학적 및 운동학적 연직 분석에서, 하층에서 온위의 연직 경도에 의한 강한 바람과 대류권계면 접힘에 의한 위치 소용돌이도의 발달이 영동지역 강풍 발생에 큰 역할을 한 것으로 사료된다.

기상현상관측은 기상청에서 다양한 방법(지상, 고층, 해양, 항공, 등)으로 관측되고 있다. 하지만, 인간생활에 많은 영향을 미치는 대기경계층 관측에는 한계가 있다. 특히, 존데 또는 항공기를 이용한 기상관측은 경제적인 측면에서 상당 한 비용이 필요하다. 따라서 본 연구의 목적은 기상드론을 이용하여 국지기상현상 중 해륙풍 연직분포에 대한 기상 인자들을 측정하고 분석하는 것이다. 해륙풍의 공간적 분포를 연구하기 위해 보성지역 표준기상관측소의 보성종합기상탑을 포함한 다른 세 지점(해안가, 산기슭, 산중턱)에 동일한 통합기상센서를 각 드론에 탑재하였다. 2018년 8월 4일 1100 LST부터 1800 LST까지 30분 간격으로 최대 400 m 고도까지 기온, 상대 습도, 풍향, 풍속, 기압의 연직 프로파일 관측이 수행되었다. 기온, 상대 습도, 기압에 대한 기상현상의 공간적 특성은 네 지점에서 보이지 않았다. 강한 일사량 시간 대에 중간지점(~100 m)에서 강한 바람(~8 m s−1 )이 관측되었고, 오후에는 풍향이 내륙지역의 상층부터 서풍으로 바뀌었다. 기상드론을 이용하여 관측한 하부 대기층의 분석결과는 보다 정확한 기상예보 향상에 도움이 될 것으로 기대된다.

본 연구는 서울에서 강한 강수와 관련된 대기연직구조를 객관적으로 분류하고 대표 종관장과 강수 특성을 제시하고자 하였다. 이를 위해 2009년부터 2018년까지 여름철 (6~8월) 서울에서 강한 강수(>15 mm hr−1 ) 시 오산에서 비양 된 레윈존데 자료에 객관적 방법을 적용하여 대기연직구조를 분류하였다. 그 결과 대기 전체가 습윤한 형태인 Thin Tube (TT) 형이 34.7% (17회), 건조한 하층 위로 습윤한 층이 존재하는 Inverted V (IV) 형이 20.4% (10회), 습윤한 하층 위로 건조한 공기가 침투하는 Loaded Gun (LG)이 20.4% (10회)로 분석되었다. TT형의 경우 SRH 값이 357.6 J kg−1으로 역학적 불안정이 큰 반면, LG와 IV형의 경우 1000 hPa부터 600 hPa까지 열적 불안정이 큰 특징을 보였다. NCEP/FNL 자료를 사용한 합성장 분석에서 TT형의 경우 기압골 전면(500 hPa)인 서해상에 저기압이 위치하여 (850 hPa) 저기압이 강화될 수 있는 종관 패턴이 형성되었다. IV와 LG형의 경우 북만주와 중국의 북동에 강한 저기압이 위치하는 종관 패턴을 보이며, 기압경도에 의한 남서기류의 유입이 상대적으로 약하였다. 강수 전반부에 강수가 집중되는 형태가 모든 유형에서 나타났으며, 특히 IV형의 경우 강수 전반 높은 강도로 강수가 집중되어 내린다. TT형의 경우 가장 많은 강수량(123.9 mm)을 보였지만 다른 유형과 비교하였을 때, 강수가 전·후반 고르게 오랜 시간 지속되는 특징을 보였다. 본 연구 결과는 서울에서 강수와 관련된 고층관측자료의 이해도를 높이는 동시에 강수 예보기술 발전에 기여할 수 있을 것이다.

이 연구의 목적은 강남 선정릉지역에서 전산유체역학모델(CFD)을 사용하여 도시지역의 흐름 및 열 환경 모의를 검증하는 것이고, CFD 모델의 모의결과와 선정릉 지역의 관측 자료와 비교하는 것이다. CFD 모델은 국립기상과학원과 서울대가 공동으로 연구 개발된 모델이다. CFD_NIMR_SNU 모델은 기상청 현업 모델인 국지예보모델(LDAPS)의 바람성분과 온도성분을 초기 및 경계조건으로 적용되었고 수목효과와 지표 온도를 고려하여 2015년 8월 4일에서 6일까지 강남 선정릉 지역을 대상으로 수치실험을 진행하였다. 선정릉지역에서 수목효과 적용 전후의 풍속을 비교하였을 때 평균 제곱근 오차(RMSE)는 각각 1.06, 0.62 m s−1로 나타났고 수목효과 적용으로 풍속 모의정확도가 향상되었다. 기온은 LDAPS 과소 모의하는 경향을 나타내고 CFD_NIMR_SNU 모델에 의해 향상된 것을 확인하였다. CFD_NIMR_SNU 모델을 이용하여 복잡한 도시지역의 흐름과 열 환경을 자세하고 정밀한 분석이 가능하며, 도시 환경 및 계획에 대한 정보를 제공 할 수 있을 것이다.

본 연구에서는 구미에 위치한 낙동강에서 발생한 안개를 분류하였고 안개 발생 전후의 기상특성을 분석하였다. 안개는 2013년부터 2015년까지 시정계를 이용하여 관측되었다. 안개는 총 74회 발생하였고 대부분 증기안개로 분류되 었다. 관측된 증기 안개는 내륙에서 발생한 안개보다 지속시간이 길게 나타났는데 이는 지형적 특징으로 인해 야간 증 발이 다른 지역보다 강하게 나타났기 때문이다. 안개 지속시간에 대한 하천 증발 효과를 분석하기 위해 Penman- Monteith (Penman법)와 Bulk aerodynamic (Bulk법) 방법을 사용하여 증발량을 추정하였다. 이 중 Bulk법은 실제 수면 에서 측정한 증발량과 유사하게 나타났다. 따라서 Bulk법이 실제 수면 증발량 추정에 적합한 방법임을 확인할 수 있었 다. Bulk법으로 추정한 증발량은 안개 비발생일 보다 안개 발생일에 06 LST와 07 LST에 더 높게 나타났다. 안개 발 생일에 하천의 증발은 대기에 증발잠열 에너지를 공급하고 안개 내부의 난류를 유지하는 에너지원으로 작용한다. 이와 같은 결과는 안개내부의 풍속, 기온, 그리고 난류운동에너지의 증가를 통해 확인하였다.

중규모 기상 모델을 이용하여 안개와 같은 미세규모 국지현상을 정확히 재현하는 것은 매우 어려운 실정이다. 특히, 수치모델의 초기 입력 자료의 불확도는 수치모델의 예측 정확도에 결정적인 영향을 미치기 때문에 이를 보완하기 위한 자료동화 과정이 요구되어진다. 본 연구에서는 WRF (Weather Research and Forecasting) 모델을 이용하여 낙동강 지역에서 발생한 여름철 안개사례 재현실험을 대상으로 중규모 기상 모델의 한계를 검증하였다. 중규모 기상 모델에서 초기 및 경계장으로 사용되는 KLAPS (Korea Local Analysis and Prediction System)와 LDAPS (Local Data Assimilation and Prediction System) 분석장 자료를 이용하여 수치모델 모의 정확도 민감도 분석을 수행하였다. 또한 AWS (Automatic Weather System) 자료를 이용한 자료동화(Four-Dimensional Data Assimilation)에 의한 수치모델의 정확도 개선 정도를 평가하였다. 초기 및 경계장 민감도 분석 결과에서 LDAPS 자료를 입력 자료로 사용한 경우가 KLAPS 자 료 보다 기온과 이슬점온도, 상대습도에서 높은 정확도를 보였고, 풍속은 더 낮은 수준을 나타내었다. 특히, 상대습도에 서 LDAPS의 경우는 RMSE (Root Mean Square Error)가 15.9%, KLAPS는 35.6%의 수준을 보여 그 차이가 매우 크 게 나타났다. 또한 자료동화를 통하여 기온, 풍속, 상대습도의 RMSE가 각각 0.3 o C, 0.2ms−1 , 2.2% 수준으로 개선되었다.

본 연구에서는 겨울철 대관령지역의 지형성 구름에 대해 인공증설을 위한 구름씨뿌리기(이하 시딩) 영향을 알아보기 위해서 2013년 3월 13일 실험사례를 분석하였다. 지상연소기를 이용하여 기온 −4oC 이하, 풍향 45-130o, 풍속 5 m s−1 이하일 때 AgI 입자를 시딩 하였으며 대관령지역에서 적절한 시딩량을 알아보기 위해 38 g h−1 (SR1)과 113 g h−1 (SR2)에 대해 실험을 수행하였다. AgI point-source 모듈을 추가한 WRF (Weather Research and Forecast) 수치모의실험을 통해 시딩 물질의 확산장을 알아보았다. 수치모의 결과 과냉각수적이 충분히 존재한 상태에서 실험이 실시되었으며 시딩 물질은 주풍에 따라 이동하는 경향을 보였다. 시딩 효과를 알아보기 위해 안개입자측정기, 강수입자측정기와 광학우적계에서 관측된 자료를 분석하였다. 본 연구사례에서는 빙정핵 시딩에 의해 1 mm 이하 크기의 강수입자 수농도의 증가가 나타났으며 대관령지역에는 SR1 시딩이 더 적절하다고 판단된다.

태풍 발생빈도에서 1999-2013년 동안의 태풍 발생빈도는 1977-1998년 동안의 태풍보다 열대 및 아열대 서태평양의 북서해역에서 더 많이 발생하는 경향이 확인되었다. 또 1977-1998년 동안의 태풍은 주로 필리핀 동쪽 먼 해상에서 필리핀 및 남중국해를 지나 인도차이나 반도를 향해 서쪽으로 이동하거나 필리핀 동쪽 먼 해상에서 일본 동쪽 먼 해상으로 북상하는 경향을 보였다. 반면에 1999-2013년 동안에 태풍들은 주로 동아시아 중위도 지역으로 북상하는 패턴을 나타내었다. 따라서 1999-2013년 동안의 태풍들이 1977-1998년 동안의 태풍들보다 훨씬 고위도로 이동하는 경향이 있으며, 결국 후자의 기간보다 전자의 기간에 태풍 최대강도가 고위도에서 나타날 가능성이 높음을 알 수 있었다. 두 기간 사이에 500 hPa 유선에 대한 차에서 30-50˚N에서는 고기압성 순환 아노말리가 강화되어 있는 반면 남중국해의 북쪽에는 몬순 기압골 아노말리가 위치해 있으며, 이 몬순 기압골 아노말리는 145˚E까지 동쪽으로 확장되어 있었다. 이고기압성 순환 아노말리와 몬순 기압골 아노말리에 의해 동아시아 중위도 지역은 남동풍 아노말리의 영향을 받고 있으며, 이 남동풍 아노말리는 태풍들을 동아시아 중위도 지역으로 향하게 하는 지향류 아노말리의 역할을 하게 되어 1999-2013년 동안의 태풍들이 1977-1998년 동안의 태풍들보다 최대 강도의 위도가 증가할 수 있었다.

이 연구는 최근 37년 동안 여름철 한국 부근 지역에 영향을 준 태풍빈도와 북서태평양 몬순(western North Pacific monsoon index, WNPMI)과의 상관을 분석하였다. 두 변수 사이에는 뚜렷한 양의 상관관계가 존재하였으며, 엘니뇨-남방진동(El Niño-Southern Oscillation, ENSO) 해를 제외하여도 높은 양의 상관관계는 변하지 않았다. 이러한 두 변수 사이에 양의 상관관계의 원인을 알아보기 위해 ENSO해를 제외하고 가장 높은 북서태평양 몬순지수를 갖는 8해(양의 북서태평양 몬순지수 해)와 가장 낮은 북서태평양 몬순지수를 갖는 8해(음의 북서태평양 몬순지수 해)를 선정하여 두 그룹 사이에 평균 차를 분석하였다. 양의 북서태평양 몬순지수 해에는 태풍들이 열대 및 아열대 북서태평양의 동쪽해역에 주로 발생하여 동중국해를 지나 한국 및 일본을 향해 북상하는 경향을 나타내었다. 음의 북서태평양 몬순지수 해에는 태풍들이 열대 및 아열대 북서태평양의 서쪽해상에 주로 발생하여 남중국해를 지나 중국 남동부 해안 및 인도 차이나 반도지역을 향해 서진하는 패턴을 보였다. 따라서 한국 부근 지역까지 먼 거리를 이동하면서 바다로부터 충분한 에너지를 얻을 수 있는 양의 북서태평양 몬순지수 해에의 태풍강도가 더 강하였다. 또한 양의 북서태평양 몬순지수 해에 태풍들이 더 많이 발생하는 특성을 보였다. 850 hPa과 500 hPa 유선에 대한 두 그룹 사이에 차에서는 열대 및 아열대 북서태평양에서 저기압성 아노말리가, 동아시아 중위도 지역에는 고기압성 아노말리가 강화되었다. 이 두 기압계 아노말리로 인해 한국 부근 지역에서는 남동풍 아노말리가 발달하였으며, 이 남동풍 아노말리가 태풍들을 한국 부근 지역으로 향하게 하는 지향류 아노말리의 역할을 하였다. 또한 열대 및 아열대 북서태평양에서 발달한 저기압성 아노말리로 인해 양의 북서태평양 몬순 지수해에 태풍들이 좀 더 많이 발생할 수 있었다.

이 연구는 6월 남극진동이 한국 6월 강우량에 영향을 주는지 알아보기 위해 두 변수 사이에 상관분석이 이루어졌고 높은 양의 상관관계가 있음이 밝혀졌다. 이는 한국 6월 강우량이 같은 시기에 남반구에서 강화되는 머스커렌 고기압과 호주 고기압의 영향을 받음을 의미하는 것이다. 이 두 고기압이 발달할 때 호주 주위지역으로부터 적도방향으로 적도 횡단류가 강화되며, 이 적도 횡단류의 강화는 북태평양고기압의 북쪽으로의 강화로 이어진다. 이는 결국 장마전선을 한국으로 북상시키는 역할을 한다. 더욱이 북태평양 고기압의 북쪽으로의 강화는 한국에 상륙하거나 영향을 주는 태풍의 빈도를 증가시켜 6월 강우량의 증가에 중요한 역할을 한다.

저온 습윤한 1993년 여름철과 고온 건조한 1994년 여름철과 관련된 동아시아에서 대조적인 여름철 몬순순환의 특성을 상, 하층 대류권의 대기순환 특성과 함께 전구 해수온도 및 적도 대류성 강수장을 분석함으로써 조사하였다. 1993년의 경우, 동아시아, 중앙 북태평양 및 미국 서부지역에서 500hpa 면과 200hpa 면의 음의 지위고도 편차가 나타났지만, 1994년의 경우, 이들 지역들에서 양의 편차를 보였다. 1993년의 아열대 제트류는 평년보다 다소 남쪽에 치우쳐져 한반도 북쪽에 위치하였다. 서태평양 아열대 고기압이 남쪽으로 이동하여 동아시아지역에는 평년보다 많은 여름철 강수와 낮은 여름철 기온이 나타났다. 이는 오오츠크해로부터 동해로 저온 습윤한 기단의 확장에 때문으로 판단된다. 대조적으로 1994년의 아열대 제트류는 평년보다 다소 북쪽에 위치하였고, 서태평양 아열대 고기압의 갑작스런 북상은 동아시아 여름철 강수대의 북상을 동반하였다. 따라서, 아시아 여름철 강수 및 기온 편차는 1993년과 반대 양상을 보였다. 적도 태평양상의 해수온 편차에서는, 1993년은 엘니뇨가, 1994년은 라니냐가 각각 나타났다. 오스트레일리아 고기압과 마스카렝 고기압의 북서 연변을 따른 하층 적도 횡단류와 관련된 서태평양과 인도양에서 이상적인 대류성 강수는 이들 대조적인 동아시아 여름철의 대규모 대기순환에 영향을 준 것으로 판단된다. 또한, 동아시아 여름몬순지역에서 평균된 200 hPa 면의 동서바람편차는 한반도 여름 기온편차와 음의 상관을 보였다.

상층과 하층의 대기구조 특성을 알아내기 위하여, 해양의 특성을 잘 표현한다고 사료되는 섬(충남 보령시 오천면 외연도리)을 관측지점으로 선택하여 라디오존데를 이용하여 대기의 연직구조를 관측하였다. 관측결과를 이용하여, 저 고기압이 통과할 때의 상대습도, 기온과 바람의 변화를 분석하였으며 각각의 경우에 대해 결과치들을 비교하였다. 해양과 육상에서의 대기구조를 비교하기 위하여 오산의 관측 결과를 분석하였다. 그 결과, 상대습도가 고 저기압의 통과에 따라 변화함을 알 수 있었다. 즉, 저기압이 다가올 때는 상대습도는 하층보다 상층에서 먼저 증가하였으나 고기압이 다가올 때는 상층에서 10%까지 급격하게 감소하였다. 또, 상대습도 변화폭이 육상에서보다 큰 것으로 분석되었다. 앞으로 해양과 육상에서의 수증기간의 관계를 계절별, 그리고 일별로 동시간 관측을 통해 찾고자 한다.

Better understanding the mechanism of black ice occurrence on the road in winter is necessary to reduce the socio-economic damage it causes. In this study, intensive observations of road weather elements and surface information under the influence of synoptic high-pressure patterns (22nd December, 2020 and 29th January, and 25th February, 2021) were carried out using a mobile observation vehicle. We found that temperature and road surface temperature change is significantly influenced by observation time, altitude and structure of the road, surrounding terrain, and traffic volume, especially in tunnels and bridges. In addition, even if the spatial distribution of temperature and road surface temperature for the entire observation route is similar, there is a difference between air and road surface temperatures due to the influence of current weather conditions. The observed road temperature, air temperature and air pressure in Nongong Bridge were significantly different to other fixed road weather observation points.

The synoptic structural characteristics associated with heavy snowfall (Bukgangneung: 31.3 cm) that occurred in the Yeongdong area on 20 January 2017 was investigated using surface and upper-level weather charts, European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis data, radiosonde data, and Moderate Resolution Imaging Spectroradiometer (MODIS) cloud product. The cold dome and warm trough of approximately 500 hPa appeared with tropopause folding. As a result, cold and dry air penetrated into the middle and upper levels. At this time, the enhanced cyclonic potential vorticity caused strong baroclinicity, resulting in the sudden development of low pressure at the surface. Under the synoptic structure, localized heavy snowfall occurred in the Yeongdong area within a short time. These results can be confirmed from the vertical analysis of radiosonde data and the characteristics of the MODIS cloud product.

We developed a small sensor observation system (SSOS) at a relatively low cost to observe the atmospheric boundary layer. The accuracy of the SSOS sensor was compared with that of the automatic weather system (AWS) and meteorological tower at the Korea Meteorological Administration (KMA). Comparisons between SSOS sensors and KMA sensors were carried out by dividing into ground and lower atmosphere. As a result of comparing the raw data of the SSOS sensor with the raw data of AWS and the observation tower by applying the root-mean-square-error to the error, the corresponding values were within the error tolerance range (KMA meteorological reference point: humidity ± 5%, atmospheric pressure ± 0.5 hPa, temperature ± 0.5℃. In the case of humidity, even if the altitude changed, it tends to be underestimated. In the case of temperature, when the altitude rose to 40 m above the ground, the value changed from underestimation to overestimation. However, it can be confirmed that the errors are within the KMA’s permissible range after correction.

To analyze the cooling effect of urban green areas, we conducted micrometeorological measurements in these areas and their surroundings in Seoul, Korea. From the average hourly temperature measurements through each month for the last two years (March 2013 to February 2015), we found that the maximum temperature difference between urban and green areas was about 2.9℃ at 16:00 LST in summer, and the minimum was about 1.7℃ at 22:00 LST in winter. In summer, the temperature difference was the largest during the day, rather than at night, due mainly to shading by the tree canopy. The specific humidity difference between the two areas was about 1.5 g kg-1 in summer, and this decreased in the winter. The specific humidity difference between urban and green areas in summer is relatively large during the day, due to the higher evapotranspiration level of biologically active plants.

We developed the Parameter-elevation Regressions on Independent Slopes Model(PRISM)-based Dynamical downscaling Error correction(PRIDE)-Wind speed(WS) model version 3.0 to produce highresolution( 1km) grid data at a monthly time scale by using observation and Regional Climate Model(RCM) wind speed data. We consequently produced monthly wind speed grid data during the observation period(2000~2014) and the future period(2021~2100) for Representative Concentration Pathway(RCP) 2 type scenarios by using the PRIDE-WS model. The PRIDE-WS model is constructed by combining the MK(Modified Korean)-PRISM-Wind, the RCM anomaly and Cumulative Density Function(CDF) fitting, basically based on Kim et al.(2016)’s algorithm applied for daily precipitation. The upper level wind(80m altitude) was estimated by Deacon equation using surface wind speed that was produced by the PRIDE-WS model. The results show that the wind speed at the upper level generally increased during the summer season while it decreased during the spring, autumn and winter seasons.

The occurrence of heat waves estimated on historical runs of climate change was compared to that on reanalysis data from 1981 to 2005. Heat waves in the future then were predicted on the basis of climate change scenarios from 2006 to 2100. For the past period, the heat wave days predicted from the climate change scenarios data overestimated and than those by the reanalysis data. For the future period, the heat wave days increased until the mid-21st century and then stay stagnant by the RCP 2.6 scenario. However, the yearly heat wave days steadily increased until 2100 by the RCP 8.5 scenario. The synoptic cause of the most severe year of the heat wave days was analyzed as a strong high pressure developed around the Korean peninsula. The high pressure under the RCP 2.6 scenario was caused by the high level jet stream in the border area between China and Russia, whereas the high pressure under the RCP 8.5 scenario was caused by the strong high level jet stream and pressure ridge in the East Sea.

We analyzed diurnal variations in the surface air temperature using the high density urban climate observation network of Daegu in summer, 2013. We compared the time elements, which are characterized by the diurnal variation of surface air temperature. The warming and cooling rates in rural areas are faster than in urban areas. It is mainly due to the difference of surface heat capacity. In addition, local wind circulation also affects the discrepancy of thermal spatiotemporal distribution in Daegu. Namely, the valley and mountain breezes affect diurnal variation of horizontal distribution of air temperature. During daytimes, the air(valley breeze) flows up from urban located at lowlands to higher altitudes of rural areas. The temperature of valley breeze rises gradually as it flows from lowland to upland. Hence the difference of air temperature decreases between urban and rural areas. At nighttime, the mountains cool more rapidly than do low-lying areas, so the air(mountain breeze) becomes denser and sinks toward the valleys(lowlands). As the result, the air temperature becomes lower in rural areas than in urban areas.

Observation data (1981-2014) and climate change scenario data (historical: 1981-2005; RCP 2.6 and 8.5: 2006-2100) were used to analyze occurrence and future outlook of the extreme heat days and tropical nights in Daegu and Jeju. Then we compared the mortality and observations data (1993-2013). During 1981-2014, the average of extreme heat days (tropical nights) was 24.41 days (12.47 days) in Daegu, and 6.5 days (22.14 days) in Jeju. Extreme heat days and tropical nights have been similarly increased in Daegu, but tropical nights increased more than extreme heat days in Jeju. Extreme heat days and tropical nights in both, Daegu and Jeju showed high correlation with daily mortality, specifically Daegu’s correlation was higher than that of jeju. The yearly increasing rate of extreme heat of the future (2076-2100) was 1.7-3.6 times and 7.8-37.7 times higher than the past (1981-2005) in Daegu and Jeju, respectively. The yearly increase rate of tropical nights of future was 2.6-5.0 times and 2.9-5.6 times higher in Daegu and Jeju, respectively. During 2006-2100 periods, the trend of extreme heat days was observed both in Daegu and Jeju. On the average, extreme heat days and tropical nights in Jeju increased more than that of Daegu. However, the trend of extreme heat days increase in Daegu was higher than that in Jeju, whereas, the trend of tropical nights in Jeju was higher than that in Daegu.