• Atmosphere
• /
• v.23 no.1
• /
• pp.33-37
• /
• 2013

In meteorology, exploitation of Kalman filter as a data assimilation system is virtually impossible due to simultaneous requirements of adjoint model and large computer resource. The other substitute of utilizing ensemble Kalman filter is only affordable by compensating an enormous usage of computing resource. Furthermore, the latter employs ensemble integration sets for evolving the background error covariance matrix by compensating the dynamical feature of the temporal evolution of weather conditions. We propose a new implementation method that works without the adjoint model by utilizing the explicit expression of the background error covariance matrix in backward evolution. It will also break a barrier in the evolution of the covariance matrix. The method may be applied with a slight modification to the real time assimilation or the retrospective analysis.

• Ocean and Polar Research
• /
• v.30 no.2
• /
• pp.129-140
• /
• 2008

As a preliminary effort to establish a data assimilative ocean forecasting system, we reviewed the theory of the Ensemble Kamlan Filter (EnKF) and developed practical techniques to apply the EnKF algorithm in a real ocean circulation modeling system. To verify the performance of the developed EnKF algorithm, a wind-driven double gyre was established in a rectangular ocean using the Regional Ocean Modeling System (ROMS) and the EnKF algorithm was implemented. In the ideal ocean, sea surface temperature and sea surface height were assimilated. The results showed that the multivariate background error covariance is useful in the EnKF system. We also tested the sensitivity of the EnKF algorithm to the localization and inflation of the background error covariance and the number of ensemble members. In the sensitivity tests, the ensemble spread as well as the root-mean square (RMS) error of the ensemble mean was assessed. The EnKF produces the optimal solution as the ensemble spread approaches the RMS error of the ensemble mean because the ensembles are well distributed so that they may include the true state. The localization and inflation of the background error covariance increased the ensemble spread while building up well-distributed ensembles. Without the localization of the background error covariance, the ensemble spread tended to decrease continuously over time. In addition, the ensemble spread is proportional to the number of ensemble members. However, it is difficult to increase the ensemble members because of the computational cost.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.25-25
• /
• 2023

The ensemble optimal interpolation (EnOI) scheme is a sub-optimal alternative to the ensemble Kalman filter (EnKF) with a reduced computational demand making it potentially more suitable for operational applications. Since only one model is integrated forward instead of an ensemble of model realizations, online estimation of the background error covariance matrix is not possible in the EnOI scheme. In this study, we investigate two Gaussian noise based ensemble generation strategies to produce dynamic covariance matrices for assimilation of water level observations into a distributed hydrological model. In the first approach, spatially correlated noise, sampled from a normal distribution with a fixed fractional error parameter (which controls its standard deviation), is added to the model forecast state vector to prepare the ensembles. In the second method, we use an adaptive error estimation technique based on the innovation diagnostics to estimate this error parameter within the assimilation framework. The results from a real and a set of synthetic experiments indicate that the EnOI scheme can provide better results when an optimal EnKF is not identified, but performs worse than the ensemble filter when the true error characteristics are known. Furthermore, while the adaptive approach is able to reduce the sensitivity to the fractional error parameter affecting the first (non-adaptive) approach, results are usually worse at ungauged locations with the former.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.21 no.3
• /
• pp.403-412
• /
• 2018

In this study, a three-dimensional variational(3DVAR) data assimilation system was developed for the operational numerical weather prediction(NWP) system at the Republic of Korea Air Force Weather Group. The Air Force NWP system utilizes the Weather Research and Forecasting(WRF) meso-scale regional model to provide weather information for the military service. Thus, the data assimilation system was developed based on the WRF model. Experiments were conducted to identify the nested model domain to assimilate observations and the period appropriate in estimating the background error covariance(BEC) in 3DVAR. The assimilation of observations in domain 2 is beneficial to improve 24-h forecasts in domain 3. The 24-h forecast performance does not change much depending on the estimation period of the BEC in 3DVAR. The results of this study provide a basis to establish the operational data assimilation system for the Republic of Korea Air Force Weather Group.

• Journal of Korea Multimedia Society
• /
• v.18 no.7
• /
• pp.886-894
• /
• 2015

In this paper, we conduct sensitivity analysis of parameters used for inverse modeling in order to estimate the PM₁₀ emissions from the 16 areas in East Asia accurately. Parameters used in sensitivity analysis are R, the observational error covariance matrix, and B, a priori (background) error covariance matrix. In previous studies, it was used with the predetermined parameter empirically. Such a method, however, has difficulties in estimating an accurate emissions. Therefore, an automatically determining method for the most suitable value of R and B with an error measurement criteria and posteriori emissions accuracy is required. We determined the parameters through a sensitivity analysis, and improved the accuracy of posteriori emissions estimation. Inverse modeling methods used in the emissions estimation are pseudo inverse, NNLS (Nonnegative Least Square), and BA(Bayesian Approach). Pseudo inverse has a small error, but has negative values of emissions. In order to resolve the problem, NNLS is used. It has a unrealistic emissions, too. The problems are resolved with BA(Bayesian Approach). We showed the effectiveness and the accuracy of three methods through case studies.

• Atmosphere
• /
• v.32 no.4
• /
• pp.263-276
• /
• 2022

The observation error of satellite radiation data that assimilated into the Korean Integrated Model (KIM) was diagnosed by applying the Hollingsworth and Lönnberg and Desrozier techniques commonly used. The magnitude and correlation of the observation error, and the degree of contribution for the satellite radiance data were calculated. The observation errors of the similar device, such as Advanced Technology Microwave Sounder (ATMS) and Advanced Microwave Sounding Unit-A shows different characteristics. The model resolution accounts for only 1% of the observation error, and seasonal variation is not significant factor, either. The observation error used in the KIM is amplified by 3-8 times compared to the diagnosed value or standard deviation of first-guess departures. The new inflation value was calculated based on the correlation between channels and the ratio of background error and observation error. As a result of performing the model sensitivity evaluation by applying the newly inflated observation error of ATMS, the error of temperature and water vapor analysis field were decreased. And temperature and water vapor forecast field have been significantly improved, so the accuracy of precipitation prediction has also been increased by 1.7% on average in Asia especially.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.6
• /
• pp.458-467
• /
• 2019

The purpose of this study is the suggestion of optimized parameters in OI (Optimal Interpolation) by experimental study. The observation of applying optimal interpolation is ADCP (Acoustic Doppler Current Profiler) data at the southwestern sea of Korea. FVCOM (Finite Volume Coastal Ocean Model) is used for the barotropic model. OI is to the estimation of the gain matrix by a minimum value between the background error covariance and the observation error covariance using the least square method. The scaling factor and correlation radius are very important parameters for OI. It is used to calculate the weight between observation data and model data in the model domain. The optimized parameters from the experiments were found by the Taylor diagram. Constantly each observation point requires optimizing each parameter for the best assimilation. Also, a high accuracy of numerical model means background error covariance is low and then it can decrease all of the parameters in OI. In conclusion, it is expected to have prepared the foundation for research for the selection of ocean observation points and the construction of ocean prediction systems in the future.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.34 no.5C
• /
• pp.527-539
• /
• 2009

In this paper, we propose an asynchonous ICI (Inter-Cell Interference) mitigation techniques for 3GPP LTE MIMO-OFDMA down-link receiver. An increasing in symbol timing misalignments may occur relative to sychronous network as the result of BS (Base Station) timing differences. Such symbol synchronization errors that exceed the guard interval or the cyclic prefix duration may result in MAI (Multiple Access Interference) for other carriers. In particular, at the cell boundary, this MAI becomes a critical factor, leading to degraded channel throughput and severe asynchronous ICI. Hence, many researchers have investigated the interference mitigation method in the presence of asynchronous ICI and it appears that the knowledge of the SCM (Spatial Covariance Matrix) of the asynchronous ICI plus background noise is an important issue. Generally, it is assumed that the SCM estimated by using training symbols. However, it is difficult to measure the interference statistics for a long time and training symbol is also not appropriate for MIMO-OFDMA system such as LTE. Therefore, a noise reduction method is required to improve the estimation accuracy. Although the conventional time-domain low-pass type weighting method can be effective for noise reduction, it causes significant estimation error due to the spectral leakage in practical OFDM system. Therefore, we propose a time-domain sinc type weighing method which can not only reduce the noise effectively minimizing estimation error caused by the spectral leakage but also implement frequency-domain moving average filter easily. By using computer simulation, we show that the proposed method can provide up to 3dB SIR gain compared with the conventional method.

• Journal of information and communication convergence engineering
• /
• v.8 no.1
• /
• pp.23-28
• /
• 2010

According to increasing projection order, the AP algorithm bas noise amplification problem in large background noise. This phenomenon degrades the performances of the AP algorithm. In this paper, we analyze convergence characteristic of the AP algorithm and then suggest a noise robust modified AP algorithm for reducing this problem. The proposed algorithm normalizes the update equation to reduce noise amplification of AP algorithm, by adding the multiplication of error power and projection order to auto-covariance matrix of input signal. By computer simulation, we show the improved performance than conventional AP algorithm.

• Atmosphere
• /
• v.17 no.3
• /
• pp.217-230
• /
• 2007

Using the Lorenz-95 simple model, which can simulate many atmospheric characteristics, we compare the performance of ensemble strategies such as bred vectors, the bred vectors rotated (to be orthogonal to each bred member), and the Ensemble Transform Kalman Filter (ETKF). The performance metrics used are the RMSE of ensemble means, the ratio of RMS error of ensemble mean to the spread of ensemble, rank histograms to see if the ensemble member can well represent the true probability density function (pdf), and the distribution of eigen-values of the forecast ensemble, which can provide useful information on the independence of each member. In the meantime, the orthogonal bred vectors can achieve the considerable progress comparing the bred vectors in all aspects of RMSE, spread, and independence of members. When we rotate the bred vectors for orthogonalization, the improvement rate for the spread of ensemble is almost as double as that for RMS error of ensemble mean compared to the non-rotated bred vectors on a simple model. It appears that the result is consistent with the tentative test on the operational model in KMA. In conclusion, ETKF is superior to the other two methods in all terms of the assesment ways we used when it comes to ensemble prediction. But we cannot decide which perturbation strategy is better in aspect of the structure of the background error covariance. It appears that further studies on the best perturbation way for hybrid variational data assimilation to consider an error-of-the-day(EOTD) should be needed.