This study attempts to improve estimation of ionospheric electron density profiles over Korea and adjacent areas by employing classical Kalman filtering technique to assimilate Total Electron Content (TEC) data from various sources into the NeQuick model. Successive corrections method was applied to spread the effect of TEC data assimilation at a given location to others that lacked TEC observations. In order to reveal that the assimilation results emulate the complex ionospheric changes during geomagnetic storms, the selected study days included both quiet (Kp ≤ 3) and disturbed geomagnetic conditions in the year 2015. The results showed that assimilation of TEC data derived from ground-based GPS receivers can improve the root mean squared error (RMSE) associated with the NeQuick model estimation of ionospheric parameters by ≥ 56 %. The improvement of RMSE achieved by assimilating TEC data that were measured using ionosondes was ~50 %. The assimilation of TEC observations made by the COSMIC radio occultation technique yielded results that depicted RMSE improvement of > 10 %. The assimilation of TEC data measured by GPS receiver onboard Low Earth Orbiting satellites yielded results that revealed deterioration of RMSE. This outcome might be due to either the fact that the receivers are on moving platforms and these dynamics might have not been accounted for during TEC computation or limitation of the assimilation process. Validation of our assimilation results with global ionosphere TEC data maps as processed at the center for orbit determination in Europe (CODE) revealed that both depicted similar TEC changes, showing response to a geomagnetic storm.