GOCI-II 대기상한 반사도와 기계학습을 이용한 남한 지역 시간별 에어로졸 광학 두께 산출 (Retrieval of Hourly Aerosol Optical Depth Using Top-of-Atmosphere Reflectance from GOCI-II and Machine Learning over South Korea)
-
- 대한원격탐사학회지
- /
- 제39권5_3호
- /
- pp.933-948
- /
- 2023
대기 중 에어로졸은 인체에 악영향을 끼칠 뿐 아니라 기후 시스템에도 직간접적인 영향을 미치므로 에어로졸의 특성과 시공간적인 분포에 대한 이해는 매우 중요하다. 이를 위해 위성기반 관측을 통해 에어로졸 광학 두께(Aerosol Optical Depth, AOD)를 산출하여 에어로졸을 모니터링하는 다양한 연구가 수행되어 왔다. 하지만 이는 주로 조견표를 활용한 역 산출 알고리즘에 기반하여 이루어지기 때문에 많은 계산량을 요구하며 불확실성이 존재한다. 따라서, 본 연구에서는 Geostationary Ocean Color Imager-II (GOCI-II)의 대기상한반사도와 30일 동안의 대기상한반사도 중 최솟값과 관측 시점 값의 차이 값, 수치 모델 기반 기상학적 변수 등을 활용하여 기계학습 기반 고해상도 AOD 직접 산출 알고리즘을 개발하였다. Light Gradient Boosting Machine (LGBM) 기법이 사용되었으며, 추정된 결과는 지상 관측 자료인 Aerosol Robotic Network (AERONET) AOD를 활용하여 랜덤, 시간 및 공간별 N-fold 교차검증을 통해 검증되었다. 세 가지 교차검증 결과 R2=0.70-0.80, RMSE=0.08-0.09, 기대오차(Expected Error, EE) 안에 있는 비율은 75.2-85.1% 수준으로 안정적인 성능을 보였다. Shapley Additive exPlanations (SHAP) 분석에서는 반사도 관련 변수들이 기여도의 상위권 대부분을 차지하고 있는 것을 통해 반사도 자료가 AOD 추정에 많은 기여를 하는 것을 확인하였다. 서울과 울산 지역에 대한 시간 별 AOD의 공간 분포를 분석한 결과, 개발된 LGBM 모델은 시간의 흐름에 따라 AERONET AOD 값과 유사한 수준으로 AOD를 추정하고 있었다. 이를 통해 높은 시공간 해상도(i.e., 시간별, 250 m)에서의 AOD 산출이 가능함을 확인하였다. 또한, 산출 커버리지 비교에서 LGBM 모델의 평균 산출 빈도가 GOCI-II L2 AOD 산출물 대비 8.8%가량 증가한 것을 통해 기존 물리모델기반 AOD 산출 과정에서 발생하던 밝은 지표면에 대한 과도한 마스킹의 문제점을 개선시킨 것을 확인하였다.
A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.
공동충전 효과를 검증하기 위하여 실시한 시차 공대공 탄성파 탐사자료로부터 지하공동 부존 지역에서 충전 전과 후에 매질의 탄성파 전파속도의 변화를 확인하였다. 시차 공대공 탄성파 탐사자료에 나타난 반응과 시추조사 결과에 의하면 본 지역의 공동은 규모가 극히 소규모이거나 또는 폐석 등으로 충전된 것으로 보인다. 공동충진 효과는 토모그래피로부터 도출된 속도단면상의 탄성파 속도의 증가량을 분석함으로써 평가하였다. 시추공용 에어건을 진원으로 24-채널 하이드로폰을 수진기로 하여 자료를 취득하였다. 취득한 자료에는 무시할 수 없을 정도의 source statics를 확인할 수 있었다. 본 논문에서 제시한 보정방법은 2단계로; 1) 불규칙한 발파시점에 의한 영향 보정과 2) 잔여 정보정으로 이는 진원의 부정확한 위치에 대한 정보정이다. 본 논문에서는 고주파수 성분의 수치잡음이 억제되고 관심대상 부분에서 비교적 고분해능 영상을 도출할 수 있는 다단계 역산 방안을 제시하였다. 일반적으로 최소자승 주시토모그래피로는 평활화된 속도 영상을 얻을 수 있다. 따라서 이러한 역산으로는 비교적 소규모의 구간에서 발생한 적은 속도변화를 영상화하기에는 어려운 면이 있다. 본 논문에서는 속도모델의 파라메터를 변화시킨 2단계 제어 역산법으로 도출한 시차 토모그램으로부터 채굴 영향대에서 발생한 매질의 속도변화를 시각화 할 수 있었다. 2단계 역산법은 1-단계에서는 적정한 크기의 균일 격자로 구성된 모델을 사용하여 토모그램을 작성하고 이 토모그램에 2차원 중위수 필터를 적용하여 대략적인 속도구조 모델을 작성한다. 2-단계 역산시는 1-단계에서 작성한 속도모델을 수정하여 초기 모델로 한다. 모델 수정은 관심대상 부분만을 작은 크기의 균일격자로 재구성하는 것이다. 기준조사 토모그램을 2차 조사자료 역산의 초기 속도모델로 사용하였다. 속도변화는 공동대 부근에서만 예상되므로 그 이외 부분의 속도는 기준 토모그램과 동일하게 고정시키고 역산을 수행하였다.
명시적 액체-얼음상 미시물리 과정을 포함하는 구름 분해 모형(ARPS: Advanced Regional Prediction System)을 이용하여 2차원 그리고 주변 바람이 없는 경우에 뇌우 유출의 구조와 진화를 조사하였다. 고해상도 격자 간격(50m)을 이용하여 유출의 난류 구조를 명시적으로 분해하였다. 모사된 단세포 스톰과 스톰과 연관된 Kelvin-Helmholtz(KH) 빌로우(billow)는 발달, 성숙, 소멸의 생애 단계를 가졌다. 구름 역학과 미시물리 사이의 상호작용으로 야기된 이차 맥동과 대류 세포의 분활이 관측되었다. 상대적으로 건조한 하층 대기를 낙하하는 빗방울과 우박의 증발에 기인한 찬 하강류는 뇌우 찬 공기 유출을 야기시켰다. 유출 머리는 거의 일정한 속도로 이동하였다. KH 불안정에 의해 생성된 KH 빌로우는 유출 상부에서 난류 혼합을 야기하였으며 유출의 구조를 지배하였다. KH 빌로우는 유출 머리에서 생성되었고 돌풍 전선에 상대적으로 뒤쪽으로 이동함에 따라 성장하고 소멸하였다. 가장 빨리 성장하는 섭동의 수평 파장과 임계 시어층 깊이의 비 그리고 KH 빌로우의 수평 파장과 최대 진폭의 비에 대한 수치 모사 결과는 다른 연구 결과와 잘 일치하였다.
본 연구에서는 GDAPS(T213) 중기 기상 수치예보 자료를 활용한 ESP (Ensemble Streamflow Prediction) 기법을 개발하여 미래에 발생할 수 있는 댐 유입량의 중장기적 확률예측을 위해 초과 확률구간별 댐 유입량을 예측하고 RPSS 검증기법으로 예측결과의 정확도를 분석하였다. 개발된 ESP시스템을 적용한 결과 일단위 개념의 확률예보는 높은 불확실성을 내포할 수 있고, 중장기 확률예보에 초점을 맞추어 1, 3, 7일 등의 예측시간 해상도에 대한 ESP정확도의 민감도를 분석한 결과 예측시간 해상도 간격이 증가할수록 예측결과의 불확실성이 감소하면서 그 정확도가 전반적으로 증가함을 살펴볼 수 있었다. 이러한 결과를 바탕으로 GDAPS 자료를 활용한 1주 단위의 한달(28일)예보를 수행한 ESP 결과는 각 초과 확률구간 분포의 적절한 증가 및 감소로 인하여 그 시간적 변동성이 안정적으로 예측되고 예측결과의 불확실성을 감소시킬 수 있어 그 활용가치가 높은 것으로 나타났다. 이러한 관점에서 본 연구의 ESP 시스템은 중장기적 측면에서 GDAPS 자료의 활용가치를 높일 수 있고, 기존 ESP 결과보다 향상된 정확도로 댐 유입량을 예측할 수 있으므로 실시간 댐 유입량 예측에 적용한다면 수자원 관리 차원에서 유용한 수단이 될 수 있을 것이다.
최근 환경, 과학, 사회 등 다양한 분야에서 원격 탐사를 활용한 연구가 활발해지고 있다. 원격탐사를 활용한 연구의 결과가 수치적인 결과뿐만 아니라 사회 및 과학적 문제 해결과 예방하는데 있어 중요한 역할을 하고 있다. 이 논문은 갯벌 환경을 대상으로 원격탐사 기술을 활용하여 제공되는 데이터와 각 데이터의 연관성을 말하고자 하는데 목적을 두고 있다. 원격 탐사 중 인공위성을 활용한 고해상도 데이터를 얻어 갯벌 조류로에 대한 정보를 알아보고자 한다. 침식과 퇴적, 간조, 만조 등의 원인으로 생성되는 갯벌의 조류로는 갯벌 경사에 대한 정보와 생태계에 대한 정보를 내재하고 있다. 따라서 갯벌의 전체적인 조류로를 분석하는 것은 매우 중요한 연구 중에 하나라고 볼 수 있다. 본 논문은 TanDEM-X를 통해 DEM(Digital Elevation Model, 수치표고모델)을 생성하며, DEM은 조류로를 생성하는데 가장 기본적인 데이터로 활용이 된다. 이번 연구 지역은 강화도 갯벌로 서해안 중부에 위치한 갯벌이다. 연구 결과로 생성된 조류로를 분석해보면 강화도 갯벌의 경사 방향과 조류로의 형태 등 다양한 정보를 파악 할 수 있다. 또한 갯벌의 입도분포, 경사도, 밀도도 등의 데이터와 본 연구에서 추출한 조류로와의 연관성이 높은 점을 고려하여 해당 연구의 결과가 앞으로 갯벌 연구에 DEM 데이터의 활용에 대한 중요성과 필요성이 높아지고 양질의 결과를 얻을 수 있다라는 기대를 가져본다.
미국 해양대기청의 CCM3 기후모델을 이용하여 마지막최대빙하기 (Last Glacial Maximum (LGM))의 동아시아기후변화를 시뮬레이션 하였다. 이용된 수치모델의 수평해상도는 약 75km로 동아시아의 비교적 상세한 지형 및 기후 기작들이 표현된다. 현재 기후는 미국 해양대기청 에서 제공한 월평균 표층해수가 경계조건으로 이용되었으며, LGM 실험은 CLIMAP 프로젝트에서 복원된 표층해수가 경계조건으로 이용되었다. 또한 LGM실험에서는 현재 의 이산화탄소농도(330 ppm) 보다 낮았던 대기 이산화탄소농도 (200 ppm)가 적용되었고, 대륙빙하를 포함한 LGM 지표지형이 표현되었다. LGM 경계조건하에서 한국 서해와 동해의 대륙사면에서는 표층 온도가 겨울철에 약 18도 정도 현재보다 낮게 나타나며, 여름철에는 약 2도 정도 높게 나타나는데, 이와 같은 온도 변화는 한국 서해가 빙하기때 해수면 하강에 따라 지표로 드러남으로서 표층 열용량의 차이가 나타나기 때문이다. 전체적으로 마지막최대빙하기에 동아시아 내륙에서는 표층온도가 약 4-6도 감소하며, 주변해역을 포함하여 약 7.1도 하강하는 것으로 시뮬레이션 되는데, 이와 같은 온도감소는 이산화탄소농도 감소에 따른 장파복사량의 증가에 기인한다. 표층의 온도감소는 물수지의 약화를 초래하는데, 겨울철 아시아의 남동부에서 강수량이 약 1-4 mm/day 감소하며, 여름철에는 중국에서의 감수감소가 나타난다. 대체로 마지막최대 빙하기동안 동북아시아에서 현재보다 약 50%정도 감수가 감소하는 것으로 시뮬레이션 된다. 증발량도 감소하는 것으로 나타나지만, 강수의 감소가 더 커서 지난빙하기동안 동아시아는 대체로 건조한 기후가 나타나며, 이와 같은 결과는 다른 지역에서 보이고 있는 결과와 잘 일치한다. 결론적으로 고해상도 수치모델은 전 지구 영역 하에서 동아시아의 기후변화를 잘 시뮬레이션 하고 있다.mu}M$에서
해저 진흙화산은 유동화/기화된 퇴적물이 표층으로 분출하여 만들어진 화산과 유사한 지형이다. 진흙화산은 지하의 열, 퇴적물이나 탄화수소를 지상으로 공급하는 공급원으로 관심을 받고 있다. 북극 캐나다 보퍼트 해의 대륙사면에는 다양한 수심에서 진흙화산이 존재한다고 알려져 있다. 수심 420 m에 위치한 MV420 진흙화산은 현재 분출하고 있는 활동성 진흙화산으로 많은 연구 대상이 되고 있다. 극지연구소에서는 쇄빙연구선 아라온호를 이용하여 MV420을 통과하는 고해상도 다중채널 탄성파 자료를 획득하였고, 진흙 분출구 주변에서 지열 관측을 수행하였다. 탄성파 자료에서는 가스하이드레이트에 의한 해저면모사반사파(bottom simulating reflector, BSR)로 추정되는 역위상 반사파를 확인하였다. 탄성파 자료의 BSR이 가스하이드레이트에 의한 반사파인지 확인하기 위하여, 정상 상태의 열방정식을 바탕으로 MV420 내부의 열구조를 수치적으로 모사하였다. 그리고 모사한 지열온도 모델을 이용하여 가스하이드레이트 안정영역의 하부 경계를 추정하였다. BSR의 깊이와 가스하이드레이트 안정영역의 하부 경계를 비교한 결과, 두 자료가 일치하며 이는 가스하이드레이트의 부존을 암시하는 지구물리학적 증거 중의 하나이다. 선행 연구 결과는 MV420의 분출구에서 표층에 가스하이드레이트가 부존한다는 것을 보여주었으며, 이번 연구의 결과는 분출구와의 거리에 따라 최대 50 m 깊이까지 가스하이드레이트가 존재할 수 있음을 시사한다.
The wall shear stress in the vicinity of end-to end anastomoses under steady flow conditions was measured using a flush-mounted hot-film anemometer(FMHFA) probe. The experimental measurements were in good agreement with numerical results except in flow with low Reynolds numbers. The wall shear stress increased proximal to the anastomosis in flow from the Penrose tubing (simulating an artery) to the PTFE: graft. In flow from the PTFE graft to the Penrose tubing, low wall shear stress was observed distal to the anastomosis. Abnormal distributions of wall shear stress in the vicinity of the anastomosis, resulting from the compliance mismatch between the graft and the host artery, might be an important factor of ANFH formation and the graft failure. The present study suggests a correlation between regions of the low wall shear stress and the development of anastomotic neointimal fibrous hyperplasia(ANPH) in end-to-end anastomoses. 30523 T00401030523 ^x Air pressure decay(APD) rate and ultrafiltration rate(UFR) tests were performed on new and saline rinsed dialyzers as well as those roused in patients several times. C-DAK 4000 (Cordis Dow) and CF IS-11 (Baxter Travenol) reused dialyzers obtained from the dialysis clinic were used in the present study. The new dialyzers exhibited a relatively flat APD, whereas saline rinsed and reused dialyzers showed considerable amount of decay. C-DAH dialyzers had a larger APD(11.70
The wall shear stress in the vicinity of end-to end anastomoses under steady flow conditions was measured using a flush-mounted hot-film anemometer(FMHFA) probe. The experimental measurements were in good agreement with numerical results except in flow with low Reynolds numbers. The wall shear stress increased proximal to the anastomosis in flow from the Penrose tubing (simulating an artery) to the PTFE: graft. In flow from the PTFE graft to the Penrose tubing, low wall shear stress was observed distal to the anastomosis. Abnormal distributions of wall shear stress in the vicinity of the anastomosis, resulting from the compliance mismatch between the graft and the host artery, might be an important factor of ANFH formation and the graft failure. The present study suggests a correlation between regions of the low wall shear stress and the development of anastomotic neointimal fibrous hyperplasia(ANPH) in end-to-end anastomoses. 30523 T00401030523 ^x Air pressure decay(APD) rate and ultrafiltration rate(UFR) tests were performed on new and saline rinsed dialyzers as well as those roused in patients several times. C-DAK 4000 (Cordis Dow) and CF IS-11 (Baxter Travenol) reused dialyzers obtained from the dialysis clinic were used in the present study. The new dialyzers exhibited a relatively flat APD, whereas saline rinsed and reused dialyzers showed considerable amount of decay. C-DAH dialyzers had a larger APD(11.70