본 연구의 목적은 위성의 밝기온도를 기반으로 한 정량적 강우량 추정기법의 개선을 위함이다. 우리나라 여름철 강우사례를 이용하여 강우추정을 위한 비선형 관계식을 개선하였다. 분석을 위해 레이더 자료로 기상청 기상레이더 관측망의 고도 1.5 km와 CMAX 반사도 합성자료를 사용하였으며, 위성자료는 천리안 위성의 가시, 적외, 수증기 채널의 자료를 이용하였다. 새롭게 도출된 알고리즘은 A-E method, CRR v4.0 analytic function의 결과와 비교를 하였다. 검증을 위해 우리나라 ASOS에서 관측한 지상강우량 자료를 사용하였다. 공간검증을 위해 검증지수로 POD, FAR, CSI를 계산하였으며 각각 0.67, 0.76, 0.21로 나타났다. 정량적 강우검증을 위해 MAE와 RMSE를 계산하였으며 각각 2.49, 6.18 mm/h였다. A-E에 비하여 정량적인 오차가 줄어들었으며 CRR에 비하여 공간적인 정확도가 증가하였다. 개선한 관계식을 적용한 방법이 두 알고리즘의 부족한 부분을 보완할 수 있는 것으로 판단된다. 개선한 관계식을 통해 강우를 추정하는 방법은 복잡한 알고리즘을 거치지 않고 짧은 시간에 강우추정이 가능함으로써 현업용 실시간 초단기 예보에 활용될 수 있다.

레이더 반사도를 이용한 강수추정의 개선을 위해 새로운 접근 방식인 경북대학교에서 개발한 하이브리드 고도면을 이용한 강수량 추정기법(Hybrid Surface Rainfall, KNU-HSR)을 사용하였다. KNU-HSR기법은 지형에코와 레이더 빔차폐의 영향을 받지 않는 2차원 하이브리드 고도면에서의 반사도를 이용하여 강수량을 추정한다. 본 연구에서는 정적 HSR 및 동적 HSR기법이 사용되었으며 비교 검증되었다. 정적 HSR은 빔차폐지도와 지형에코지도를 사용하며, 동적 HSR은 정적 HSR에 추가적으로 실시간 퍼지로직 품질관리를 통한 품질지수지도를 사용한다. 검증을 위해 상관계수(correlation coefficient), 총비율(total ratio), 평균편의(mean bias), 정규화된 표준편차(normalized standard deviation), 평균 상대오차(mean relative error)를 사용하였으며, 10개 강우사례의 지상우량계 강우자료를 이용하여 두 HSR의 강우추정 성능을 평가하였다. 모든 검증지수에서 동적 HSR은 반사도 보정을 하지 않은 정적 HSR에 비해 더 우수한 성능을 보였다. 동적 HSR은 레이더로부터 근거리에서는 과대추정하였으며 원거리에서는 빔 폭 확장 및 빔 고도증가로 인해 과소추정하였다. 동적 HSR의 정규화된 표준편차와 평균상대오차는 레이더로부터의 거리에 관계없이 가장 좋은 결과를 보였다. 정적 HSR은 약한 강우강도에서 상당히 과대추정하였으나 동적 HSR은 모든 강우강도에서 1.0에 총비율을 보였다. 반사도의 시스템오차 보정 후, 동적 HSR의 정규화된 표준편차와 평균상대오차는 각각 약 20%와 15%로 개선되었다.

R-Z 관계식은 레이더 강우추정의 정확도를 결정하는 중요한 요소이다. 본 연구에서는 캐나다 궤벡주의 셍레미에서 홍수를 야기한 폭우사례에서 관측된 우적계 및 레이더 자료를 이용하여 레이더 강우추정 시 우적계 자료에서 도출된 R-Z 관계식의 효과를 분석하였다. 이를 위하여 맥길 S-밴드 레이더에서 시간 분해능 2.5분과 공간 분해능 1˚×250m로 관측된 레이더 반사도를 사용하였다. 레이더 반사도 자료에서는 폭우를 동반한 강우세포가 셍레미를 통과한 것으로 관측되었지만 우량계 관측망에서는 낮은 공간 분해능으로 인하여 이 세포가 관측되지 않았다. 셍레미에서 30분과 1시간 최대 누적 강우량은 각각 39 mm와 42 mm였다. 강우사례 동안 두 개의 우적계(POSS; Precipitation Occurrence Sensor System)가 사용되었다. 하나의 우적계는 레이더 반사도와 우적계 반사도를 비교하여 레이더 반사도를 보정하고 다른 우적계는 R-Z 관계식을 유도하는데 사용되었다. 기후학적 R-Z 관계식을 사용하였을 때 보다 반사도에 의존적인 우적계에서 유도된 관계식을 사용하였을 때 강우 추정 오차가 크게 줄었다. 일 누적 강우량에 대하여 편차는 +12%에서 -2%, 평균제곱근오차가 16%에서 10%로 줄었다. 우적계에서 도출된 R-Z 관계식으로 추정된 레이더 강우장을 이용하였을 때 홍수사례에 대하여 강우 발생 시간 및 강우량이 잘 일치하였다.

수문 ․ 기상레이더는 강우량을 바로 추정하지 못하고 여러 단계의 정량적 강우량 추정과정을 거치게 되므로 많은 불확실성 발생요소가 존재한다. 불확실성 관련한 기존 연구들은 정량적 레이더기반 강우량 추정과정에서 보정방법을 이용하여 각 단계별 불확실성을 줄이는 연구들을 수행하였다. 하지만 기존 연구들은 전체 과정에 대한 포괄적인 불확실성을 나타내지 못하고 각 단계별 불확실성의 상대적인 비율도 제시하지 못하는 단점이 있다. 본 연구에서는 정량적 레이더강우량 추정과정의 각 단계별 불확실성을 정량화하고 불확실성 전파를 나타낼 수 있는 적합한 방법을 제시하였다. 첫 번째로 초기와 최종 불확실성, 각 단계별 불확실성의 변동과 상대적인 비율을 나타낼 수 있는 새로운 개념을 제안하였다. 두 번째로 레이더기반 추정과정의 불확실성 정량화와 전파과정을 분석하기 위해 Maximum Entropy Method (MEM)와 Uncertainty Delta Method (UMD)를 적용하였다. 세 번째로 레이더기반 강우량 추정과정의 불확실성 정량화를 위해 2개 품질관리 알고리즘, 2개 강우량 추정방법, 2개 후처리 강우량 보정방법을 2012년 여름철 18개 사례에 대하여 사용하였다. 적용결과, MEM에서 최종 불확실성(후처리 강우량 보정 불확실성: ME = 3.81)이 초기 불확실성(품질관리 불확실성: ME = 4.28)보다 작게 나타났으며, UMD에서도 최종 불확실성(UMD = 4.75)이 초기 불확실성(UMD = 5.33)보다 작게 나타나 불확실성이 감소하는 것으로 나타났다. 하지만 레이더강우량 추정단계의 불확실성은 증가하는 것으로 나타났다. 또한 레이더강우량 추정과정에서 각 단계별로 적합한 방법을 선정하는 것이 각 단계별로 불확실성이 감소시킬 수 있음을 확인하였다. 따라서 본 연구는 새로운 방법이 명확히 불확실성을 정량화할 수 있으며 정확한 정량적 레이더 강우추정에 기여할 것으로 판단한다.

기상청, 국토교통부 및 국방부는 각 기관의 특성에 맞게 기상, 수문 및 항공 관측용 레이더를 운영하고 있다. 2015년까지 국내에는 총 8대의 이 중편파레이더(백령도, 용인테스트베드, 진도, 면봉산, 비슬산, 소백산, 모후산, 서대산 레이더)가 도입되어 운영되고 있다. 레이더를 운영 중인 세 개 부처는 관측목적이 달라 레이더에 대한 운영방식, 자료 처리 방식 및 활용 방안 등에서 서로 많은 차이를 보여 왔으며 이러한 차이는 각 부처에서 운영하는 레이더 자료의 정확도가 달라지게 하는 원인이 되고 있다. 이에, 세 부처는 2010년 “기상-강우 레이더 자료 공동 활용에 관한 합의서”를 체결하고, 서로 다른 부처에서 생산된 모든 레이더 자료를 공동 활용하는 체계를 마련하였다. 범부처 이중편파레이더 자료의 합성을 위해서는 부처 별 레이더 간의 정확도가 서로 유사해야한다. 서로간의 정확도가 다른 상태에서 레이더 강우 자료를 합성하게 된다면 공간적으로 품질이 다른 합성 강우장이 생산될 수밖에 없다. 이에 본 연구에서는 범부처 이중편파레이더 강우량의 정확도 향상과 통합을 위해 2015년도에 용인테스트베드 레이 더에서 시행되었던 경험적 방법을 이용하였다. 그 결과 2015년 5월부터 10월까지의 백령도, 비슬산 및 소백산 레이더 자료에 대한 강우량 정확도 (1-NE)가 70% 수준으로 향상되었다. 또한 레이더 편파변수 조절 전에는 정확도가 30∼60%까지 넓은 범위를 보이고 있으나 조절 후 65~70% 로 그 범위가 줄어들었다.

본 연구에서는 확률 대응법의 적용성을 평가하고, 이를 적용하여 실제 호우사상에 적합한 레이더 강우 추정관계식을 결정하였다. 먼저, 확률 대응법의 과정을 새로 정리하여 zT나 rT 등의 임계값을 결정하였다. 둘 째, 확률 대응법에 적용할 자료의 수가 매개변수 결정에 미치는 영향에 대하여 모의 자료를 이용하여 민감도 분석을 수행하였다. 셋 째, 확률을 나타낼 반사도와 강우강도 히스토그램의 구간 간격이 매개변수 결정에 미치는 영향에 대해 민감도 분석을 수행하였다. 그 결과, 확률대응법으로 매개변수를 결정할 경우 확률대응법에 2시간 정도 누적된 반사도와 강우 쌍을 적용하는 것이 적합함을 확인하였다. 아울러 확률대응에는 반사도와 강우강도의 누가확률 차이보다 1차 모멘트 차이를 이용하는 것이 양호한 결과를 얻을 수 있고, 30~100% 구간을 대응하는 것이 적합하다는 것을 확인하였다. 본 연구에서는 이러한 결과를 바탕으로 비슬산 레이더로 관측한 실제 호우사상에 대해 확률대응법으로 매개변수를 결정하였다. 또한, 매개변수를 이용하여 레이더 강우를 추정하였고, 이를 지상강우와 비교하였다. 전체적으로 레이더 강우가 지상강우보다 크게 산정되었지만, 레이더 강우와 지상 강우를 쉽게 대응하기 어려운 기존의 연구 결과를 고려하여 본 연구에서 평가한 확률대응법이 비교적 호우사상에 적합한 레이더 강우를 산정할 수 있다고 판단하였다.

One of main benefits of a dual polarization radar is improvement of quantitative rainfall estimation. In this paper, performance of two representative rainfall estimation methods for a dual polarization radar, JPOLE and CSU algorithms, have been compared by using data from a MOLIT S-band dual polarization radar. In addition, this paper presents evaluation of specific differential phase (Kdp) retrieval algorithm proposed by Lim et al. (2013). Current Kdp retrieval methods are based on range filtering technique or regression analysis. However, these methods can result in underestimating peak Kdp or negative values in convective regions, and fluctuated Kdp in low rain rate regions. To resolve these problems, this study applied the Kdp distribution method suggested by Lim et al. (2013) and evaluated by adopting new Kdp to JPOLE and CSU algorithms. Data were obtained from the Mt. Biseul radar of MOLIT for two rainfall events in 2012. Results of evaluation showed improvement of the peak Kdp and did not show fluctuation and negative Kdp values. Also, in heavy rain (daily rainfall > 80 mm), accumulated daily rainfall using new Kdp was closer to AWS observation data than that using legacy Kdp, but in light rain(daily rainfall < 80 mm), improvement was insignificant, because Kdp is used mostly in case of heavy rain rate of quantitative rainfall estimation algorithm.

본 연구에서는 범용 매개변수 최적화 모형인 PEST를 이용하여 분포형 수문모형인 GRM(grid based rainfall-runoff model) 모형의매개변수및불확실성범위를추정하였다. 특히, 레이더강우및지상관측강우를각각적용하여, 입력자료 차이가매개변수추정에미치는영향을분석하였다. 자동보정모형은GUI (graphic user interface)에대한접근없이모형구동이가능하도록개선된GRM-MP (multiple projects) 버전과병렬PEST버전을결합하여매개변수추정에소요되는시간을단축시켰다. 이를낙동강수계금호강유역과감천유역에대해적용하여, 초기포화도, 지표면조도계수및토양투수계수의보정계수에 대해 매개변수 최적화 및 불확실성 추정을 수행하였다. 강우자료 분석 결과, 레이더와 지상 강우의 유역평균누적시계열은비슷하거나지상강우가조금큰경향을보였으나, 공간분포에있어서는지상강우에비해레이더강우에서큰 변동성이 확인되었다. 보정된 수문모의 결과는 레이더 강우 적용 시, 지상 강우에 비해 비슷하거나 더 나은 정확도를보였다. 추정된매개변수는레이더강우적용시, 토양투수계수의보정계수가일관되게1보다작은경향을보였으며, 이는강우강도가 강한 격자가 상당수 존재하기 때문으로 판단되었다. 초기 포화도 및 지표면 조도계수의 보정계수는 레이더및 지상 강우에서 일정한 경향성을 보이지 않았다. 본 연구의 대상 유역 및 호우사상에 대한 PEST의 최적화 모의 결과,동일유역및호우사상에대해서도강우추정방법에따라서로다른최적매개변수값을갖는것을알수있었으며, 이는향후 레이더 강우 자료의 수문 모의 활용 시 유의해야할 점으로 판단된다.

본 연구에서는 관악산 레이더 자료로부터 소양강 유역에서 최적 추정된 레이더 강우로 준분포형 TOPMODEL의 적용성을 평가하고 Monte Carlo난수발생 기법과 연계한 앙상블 유량모의 기법을 구축하여 레이더 추정강우가 가질 수 있는 불확실성이 수문모형에 미치는 영향을 모형의 매개변수와 연계하여 분석하였다. 레이더 추정강우를 활용하여 소양강 유역의 유량계산을 2003년 7월의 사례기간 동안 수행한 결과 실시간 편차보정 기법으로 수정된 레이더 입력강우의

본 연구의 목적은 기상 레이더의 수문학적 활용성을 높이고자 최적의 레이더 강우를 추정하기 위해 관악산 레이더 자료를 대상으로 POD 분석기법을 활용하여 지형클러터 및 빔 차폐영역 등을 제거하고 Marshall-Palmer의 Z-R 관계식으로 레이더 강우를 추정한 결과 강우계 관측강우와 비교하여 시공간적으로 과소 모의하는 것으로 나타났다. 이러한 결과를 바탕으로 강우계 관측강우와의 실시간 보정기법을 소양강 유역을 대상으로 적용한 결과 평균 G/R 값은 로