Real Data Assimilation Using the Local Ensemble Transform Kalman Filter (LETKF) System for a Global Non-hydrostatic NWP model on the Cubed-sphere

An ensemble data assimilation system using the 4-dimensional Local Ensemble Transform Kalman Filter is implemented to a global non-hydrostatic Numerical Weather Prediction model on the cubed-sphere. The ensemble data assimilation system is coupled to the Korea Institute of Atmospheric Prediction Systems Package for Observation Processing, for real observation data from diverse resources, including satellites. For computational efficiency in a parallel computing environment, we employ some advanced software engineering techniques in the handling of a large number of files. The ensemble data assimilation system is tested in a semi-operational mode, and its performance is verified using the Integrated Forecast System analysis from the European Centre for Medium-Range Weather Forecasts. It is found that the system can be stabilized effectively by additive inflation to account for sampling errors, especially when radiance satellite data are additionally used.     

Development of a GPGPU-parallelized hybrid finite-discrete element method for modeling rock fracture

The hybrid finite-discrete element method (FDEM) is widely used for engineering applications, which, however, is computationally expensive and needs further development, especially when rock fracture process is modeled. This study aims to further develop a sequential hybrid FDEM code formerly proposed by the authors and parallelize it using compute unified device architecture (CUDA) C/C++ on the basis of a general-purpose graphics processing unit (GPGPU) for rock engineering applications. Because the contact detection algorithm in the sequential code is not suitable for GPGPU parallelization, a different contact detection algorithm is implemented in the GPGPU-parallelized hybrid FDEM. Moreover, a number of new features are implemented in the hybrid FDEM code, including the local damping technique for efficient geostatic stress analysis, contact damping, contact friction, and the absorbing boundary. Then, a number of simulations with both quasi-static and dynamic loading conditions are conducted using the GPGPU-parallelized hybrid FDEM, and the obtained results are compared both quantitatively and qualitatively with those from either theoretical analysis or the literature to calibrate the implementations. Finally, the speed-up performance of the hybrid FDEM is discussed in terms of its performance on various GPGPU accelerators and a comparison with the sequential code, which reveals that the GPGPU-parallelized hybrid FDEM can run more than 128 times faster than the sequential code if it is run on appropriate GPGPU accelerators, such as the Quadro GP100. It is concluded that the GPGPU-parallelized hybrid FDEM developed in this study is a valuable and powerful numerical tool for rock engineering applications.     

Impact of the spatial variability of daily precipitation on hydrological projections: A comparison of GCM‐ and RCM‐driven cases in the Han River basin,Korea

In this study, we investigate the impact of the spatial variability of daily precipitation on hydrological projections based on a comparative assessment of streamflow simulations driven by a global climate model (GCM) and two regional climate models (RCMs). A total of 12 different climate input datasets, that is, the raw and bias‐corrected GCM and raw and bias‐corrected two RCMs for the reference and future periods, are fed to a semidistributed hydrological model to assess whether the bias correction using quantile mapping and dynamical downscaling using RCMs can improve streamflow simulation in the Han River basin, Korea. A statistical analysis of the daily precipitation demonstrates that the precipitation simulated by the GCM fails to capture the large variability of the observed daily precipitation, in which the spatial autocorrelation decreases sharply within a relatively short distance. However, the spatial variability of precipitation simulated by the two RCMs shows better agreement with the observations. After applying bias correction to the raw GCM and raw RCMs outputs, only a slight change is observed in the spatial variability, whereas an improvement is observed in the precipitation intensity. Intensified precipitation but with the same spatial variability of the raw output from the bias‐corrected GCM does not improve the heterogeneous runoff distributions, which in turn regulate unrealistically high peak downstream streamflow. GCM‐simulated precipitation with a large bias correction that is necessary to compensate for the poor performance in present climate simulation appears to distort streamflow patterns in the future projection, which leads to misleading projections of climate change impacts on hydrological extremes.     

High sensitivity SiO maser survey for mira variables

Simultaneous observations of the six transitions of SiO for 106 late-type stars were made. The SiO maser emission was detected in 83 stars. Thev=3 maser emission was detected in eight stars, and the²⁹SiOv=0 emission in six stars. The²⁹SiOv=0 emission is stronger and narrower than that of²⁸SiO, suggesting that the²⁹SiO emission is masing.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.     

Impact of climate change on the persistent turbidity issue of a large dam reservoir in the temperate monsoon region

Long-term discharge of turbid water from reservoirs after flood events is a major socioenvironmental problem in many countries, including Korea. This study used a suite of mathematical models to simulate the fate of turbidity flows in the Soyanggang Reservoir in Korea, an important source of drinking water for the Seoul Capital Area, in response to extreme floods based on the Representative Concentration Pathway 4.5 climate scenario. It evaluated the effectiveness of the selective withdrawal facility (SWF), installed recently in the Soyanggang Reservoir to control persistent turbidity. Extreme floods with a maximum daily inflow rate greater than the historical maximum observed in 2006 were projected to occur four times in this century. The fate and transport of turbidity flows were highly influenced by both the thermal stability of the reservoir and the season in which the flood event occurred. Thus, SWF operations should consider the timing of extreme events (i.e., the imminence of the autumn turnover) to mitigate the impact of high turbidity on the water supply and downstream ecosystem. It was found to be ineffective under extreme events if these occurred in two consecutive years. Current reservoir operations, which rely heavily on the SWF, are likely to be inadequate to overcome the negative effects of extreme-turbidity events on reliably providing safe water supplies. Coping with the worst event expected to occur in the future would require additional countermeasures such as bypassing high-turbidity water.     

Effect of Data Assimilation Parameters on The Optimized Surface CO<Subscript>2</Subscript> Flux in Asia

In this study, CarbonTracker, an inverse modeling system based on the ensemble Kalman filter, was used to evaluate the effects of data assimilation parameters (assimilation window length and ensemble size) on the estimation of surface CO₂ fluxes in Asia. Several experiments with different parameters were conducted, and the results were verified using CO₂ concentration observations. The assimilation window lengths tested were 3, 5, 7, and 10 weeks, and the ensemble sizes were 100, 150, and 300. Therefore, a total of 12 experiments using combinations of these parameters were conducted. The experimental period was from January 2006 to December 2009. Differences between the optimized surface CO₂ fluxes of the experiments were largest in the Eurasian Boreal (EB) area, followed by Eurasian Temperate (ET) and Tropical Asia (TA), and were larger in boreal summer than in boreal winter. The effect of ensemble size on the optimized biosphere flux is larger than the effect of the assimilation window length in Asia, but the importance of them varies in specific regions in Asia. The optimized biosphere flux was more sensitive to the assimilation window length in EB, whereas it was sensitive to the ensemble size as well as the assimilation window length in ET. The larger the ensemble size and the shorter the assimilation window length, the larger the uncertainty (i.e., spread of ensemble) of optimized surface CO₂ fluxes. The 10-week assimilation window and 300 ensemble size were the optimal configuration for CarbonTracker in the Asian region based on several verifications using CO₂ concentration measurements.     

The influence of an excavation damaged zone on the thermal-mechanical and hydro-mechanical behaviors of an underground excavation

In Korea, a reference disposal system, KRS, was proposed in 2006 after 10 years of research and development. In the KRS, the high-level radioactive waste repository is considered to be located in a crystalline rock likes granite. For a validation of the feasibility, safety, and stability of the KRS, an underground research tunnel, KURT was constructed in Nov. 2006. During the construction of KURT by a controlled blasting, the size and characteristics of an excavation damaged zone(EDZ) were investigated by in situ as well as laboratory tests. The possible influences of an EDZ around a tunnel on the thermal, hydraulic and mechanical behaviors of the near field were investigated by using hydro-mechanical and thermo-mechanical coupling analyses. From this study, it was found that the existence of an EDZ can influence the thermal, hydraulic, and mechanical behaviors of the near field and it was recommended that an EDZ should be considered as an important parameter during the design of underground repositories.     

Effects of Topographic Slopes on Hydrological Proecsses and Climate

Based on previous research results on river re-distribution models, a modification on the effects of topographic slopes for a runoff parameterization was proposed and implemented to the NCAR＇s land sur face model (LSM). This modification has two aspects: firstly, the topographic slopes cause outflows from higher topography and inflows into the lower topography points; secondly, topographic slopes also cause decrease of infiltration at higher topography and increases of infiltration at lower topography. Then changes in infiltration result in changes in soil moisture, surface fluxes and then in surface temperature, and eventual ly in the upper atmosphere and the climate. This mechanism is very clearly demonstrated in the point bud gets analysis at the Andes Mountains vicinities. Analysis from a regional scale perspective in the Mackenzie GEWEX Study (MAGS) area, the focus of the ongoing Canadian GEWEX program, shows that the modi fied runoff parameterization does bring significant changes in the regional surface climate. More important ly, detailed analysis from a global perspective shows many encouraging improvements introduced by the modified LSM over the original model in simulating basic atmospheric climate properties such as thermodynamic features (temperature and humidity). All of these improvements in the atmospheric climate simulation illustrate that the inclusion of topographic effects in the LSM can force the AGCM to produce a more realistic model climate.     

INFLUENCES OF TOPOGRAPHY ON BAROCLINIC SOLITARY ROSSBY WAVES IN A MULTILEVEL MODEL

The KdV equation with topography included in an N-level model is derived. It is shown that if the topography ex-ists. the KdV equation may describe the solitary Rossby waves in the case of basic current without vertical shear, and itis no necessary to introduce the MKdV equation. The results of calculations show that the change of horizontal shearpattern of basic flow may cause an important change of the streamline pattern of the solitary waves with the oddmeridional wavenumber m, and has no effect for the even meridional wavenumber m. The vertical shear increases thesteepness of the barotropic solitary modes, and it has a complicated effect on the baroclinic modes. The influences oftopographic slope on the solitary waves are very great. The southern and northern slopes of topography may cause dif-ferent solitary wave patterns, with the effect of northern slope greater. The effect of Froude number on the solitarywaves is generally to steepen the solitary waves, however, the effect also depends on the meridional wavenumber m andthe modes of solitary wave.