Impact of Soil Moisture Data Assimilation on Analysis and Medium-Range Forecasts in an Operational Global Data Assimilation and Prediction System

: Accurate initial soil moisture conditions are essential for numerical weather prediction models, because they play a major role in land–atmosphere interactions. This study constructed a soil moisture data assimilation system and evaluated its impacts on the Global Data Assimilation and Prediction System based on the Korea Integrated Model (GDAPS-KIM) to improve its weather forecast skill. Soil moisture data retrieved from the Advanced Scatterometer (ASCAT) onboard the Meteorological Operational Satellite was assimilated into GDAPS-KIM using the ensemble Kalman ﬁlter method, and its impacts were evaluated for the 2019 boreal summer period. Our results indicated that the soil moisture data assimilation improved the agreement of the observations with the initial conditions of GDAPS-KIM. This led to a statistically signiﬁcant improvement in the accuracy of the initial ﬁelds. A comparison of a ﬁve-day forecast against an ERA5 reanalysis and in situ observations revealed a reduction in the dry and warm biases of GDAPS-KIM over the surface and in the lower- and mid-level atmospheres. The temperature bias correction through the initialization of the soil moisture estimates from the data assimilation system was shown in the ﬁve-day weather forecast (root mean squared errors reduction of the temperature at 850 hPa by approximately 5% in East Asia).


Introduction
As the lower boundary condition of numerical weather prediction (NWP) models, soil moisture regulates the exchange process of water and energy between land and air on the local, regional, and global scales. Vertically, soil moisture anomalies have effects on the evolution of the surface flux, planetary boundary layer, and free atmosphere [1]. Additionally, it is well-known that there are feedbacks between the soil moisture and atmospheric variables (temperature and precipitation) [2]. Moreover, previous studies have shown that soil moisture affects both the predictability of seasonal forecasts and weather forecasts [1,[3][4][5][6][7][8]. For example, Dirmeyer et al. [6] demonstrated that accurate initial soil moisture conditions improved the forecast skill of land surface fluxes, nearsurface atmospheric states, and boundary layer properties. Soil moisture anomalies are more persistent over time than those of atmospheric variables [9,10], which are referred to as the soil moisture memory and can contribute to improving weather predictions [6,[10][11][12]. Thus, an accurate estimation of the initial soil moisture conditions is essential for the predictability of NWP models.
However, soil moisture is highly variable and depends on the precipitation, land cover, and soil texture [12]. For this reason, a satellite observation with a wide hori-application of the soil moisture data assimilation in GDAPS_KIM and, also, figure out how the corrected initial soil moisture conditions improve the skills of the forecasts in the NWP through these results. This manuscript is structured as follows: Section 2 describes the soil moisture data assimilation system for GDAPS-KIM. Section 3 demonstrates the impact of this system on the initial conditions and forecasts of GDAPS-KIM, and a summary and discussion of our findings are presented in Section 4.

LSM
In this study, Noah LSM (version 3.3) [29] was used as a stand-alone LSM for soil moisture data assimilation, because KIM's LSM is based on Noah LSM. In KIM physics packages, the land surface model was developed based on Noah LSM version 3.0, with several modifications [30]. Noah LSM includes the physical processes of soil hydrodynamics [31], penman potential evaporation [32], canopy resistance [33], surface infiltration [34], and subsurface heat flux [35] and consists of four soil layers (0.1, 0.3, 0.6, and 1.0 m). In a previous study [36], Advanced Very-High-Resolution Radiometer (AVHRR) University of Maryland (UMD) land use parameter data were replaced by the Moderate Resolution Imaging Spectroradiometer (MODIS) International Geosphere-Biosphere Program (IGBP), which is used for KIM's LSM, for the stand-alone LSM. Additionally, the soil moisture is fixed to 1.0 for the glacial region in the stand-alone LSM. Beck et al. [37] mentioned that Noah LSM was almost identical to KIM-LSM, except for calculating the snow amount. Additionally, there could be inherently some differences between the two models, because KIM's LSM as a coupled model could directly exchange the energy and mass with the atmospheric model. However, a stand-alone Joint UK Land Environment Simulator (JULES) for soil moisture data assimilation was a different version with UM-LSM [38]. The stand-alone Noah LSM will be replaced by stand-alone KIM-LSM that is being constructed.
The stand-alone LSM was spun up for 2007 with 10 consecutive loops, because a previous study demonstrated that a spin-up through 10-year cycles is sufficient for the soil moisture of the LSM to reach equilibrium conditions [36]. Subsequently, a retrospective spin-up was implemented with forcing data from 2008 to 2019. The atmospheric forcing data for the LSM spin-up is the three-hourly Global Land Data Assimilation System forcing data from NASA [39,40], with a horizontal resolution of 0.25 • .

ASCAT Soil Moisture Data
ASCAT is an active microwave remote sensing instrument with a C-band wavelength of 5.3 GHz. Surface soil moisture data of ASCAT Level-2 products were disseminated in near-real time to allow the NWP centers to use them for data assimilation by EUMET-SAT [41]. In this study, ASCAT, which was on-board the first and second meteorological operational (METOP-A and B) satellites, was used. The ASCAT at a 12.5-km swath grid was chosen instead of 25 km, because it has been used for the soil moisture data assimilation of UM [38], and the system we constructed in this study needed high-resolution data for operational KMA global NWP at a 12-km horizontal resolution.
A quality control (QC) process has been implemented to exclude low-quality data where the errors for estimating the ASCAT soil wetness are greater than or equal to 16%. Moreover, the observation data are excluded when the wetland rate exceeds 15%, the complexity of the topology exceeds 20%, or the estimated error is greater than or equal to 16%. In addition, ASCAT data are rejected for freezing temperatures, precipitation, high-density vegetation (i.e., monthly greenness > 70%), or snow cover [14,42]. ASCAT data are also rejected when the MODIS IGBP land cover type belongs to one of the four forest types, including evergreen needle leaf, evergreen broad leaf, deciduous needle leaf, and deciduous broad leaf. Figure 1a,b indicates the ASCAT soil wetness tracks before and after the thinning process, respectively. The observation error standard deviations for the ASCAT unscaled soil moisture are set to the spatial and temporal constants, 10% (soil wetness index) [20,43]. This value was scaled by the ratio of the soil moisture time series standard deviation of the model to that of the ASCAT soil moisture [25,44].
During the thinning and scaling procedure, the ASCAT soil wetness (Figure 1a) is converted to the volumetric soil moisture, which represents scaled ASCAT soil moisture data ( Figure 1b) equivalent to the LSM climate ( Figure 1c). In addition, the observation and model prediction should be unbiased towards each other so that the observation biases are corrected before the assimilation is performed using the cumulative density function (CDF) matching method [13,21]. CDF matching is to match the CDF of the ASCAT soil wetness to the CDF of the model's volumetric soil moisture. To account for the seasonal variability, monthly CDFs were calculated from the soil moisture data of the observations and model simulations during 2016 and 2018. Figure 2 shows the CDF estimates and soil moisture time series of the ASCAT and LSM at 44.625 • N, 120.875 • E. The CDFs indicate a higher ASCAT soil wetness than the LSM's volumetric soil moisture ( Figure 2a); however, the scaled ASCAT soil moisture is equivalent to the LSM simulation ( Figure 2b). In addition, the variability of the scaled values is similar to that of the LSM soil moisture series.

Data Assimilation Method
An EnKF [18,45] was applied to assimilate the ASCAT observations into GDAPS-KIM. The Kalman filter updates the model state with an observation based on the uncertainty information of the model and the observations. This uncertainty is referred to as the error covariance and is estimated from an ensemble model at the time when the EnKF is updated [46]. The number of ensemble members is 20. The overall soil moisture data assimilation system for GDAPS-KIM is based on LIS version 7.2 (LIS framework, Figure 3) [35,47].   (1), H is the observation operator, and K t is the Kalman gain matrix. The model analysis vector x a t is given by averaging the updated ensemble members (Figure 3).
To produce ensembles in the EnKF, perturbations were applied to the atmospheric forcing variables, indicating the uncertainty of forcing, and the LSM prognostic states, indicating the uncertainty of LSM's physics and parameters [18,19,45]. The parameters used for generation of the perturbations are described in Table 1, and these are based on previous studies [25,45,48]. Normally distributed additive perturbations with a zero mean were applied to the long-wave radiation and LSM soil moisture, and lognormally distributed multiplicative perturbations with a mean value of 1 were employed for shortwave radiation and precipitation. There was no consideration of the spatial correlation, because we applied one-dimensional EnKF. A first-order regression model (AR1) was used for temporal correlation with time scales of 24 h for the forcing variables and 12 h for the LSM prognostic states. The cross-correlations between the atmospheric forcing variables were employed to represent the physical balance in the atmospheric modeling system. Additionally, the LSM soil moisture variables for each layer were vertically correlated. The nonlinear feature of this perturbation process can make the ensemble background biased. Thus, this unintended perturbation bias was corrected using an unperturbed model simulation [49]. Table 1. Parameters for generating perturbations to the atmospheric forcing variables and soil moisture of LSM in EnKF. P is precipitation, SW is downward shortwave radiation, LW is downward longwave radiation, and SM1-SM4 are the LSM prognostic soil moisture at the first soil layer to the fourth soil layer.

Variables
Perturbation Types

Model Configuration
The GDAPS-KIM used in this study is a global atmospheric model system, in which a spectral-element nonhydrostatic dynamical core on a cubed sphere grid is used; its horizontal and vertical resolutions are 25 km (12 km in the operational run) and 91 layers (hybrid sigma-pressure coordinate), respectively [27].
To initialize the atmospheric variables, GDAPS-KIM employs H4DEV [28] with LETKF, together with the KIM Package of Observation Processing (KPOP) [50]. GDAPS-KIM assimilates the sonde, surface, aircraft, and satellite observations. Two experiments with (EXP) and without (CTL) soil moisture data assimilation were performed during the boreal summer season. These experiments were initiated with the atmospheric conditions at 00Z 15 June 2019 and then cyclically run using H4DEV.
In the CTL experiment, the land surface condition was initialized at the start using the GFS surface analysis data and then by the model background surface conditions every 6 h during the cycling run ( Figure 4a). Conversely, the EXP land surface conditions were initialized using the soil moisture data assimilation system implemented in this study, which was conducted for 6 h at 00, 06, 12, and 18 Z cycles (Figure 4b). GDAPS-KIM was performed for five-day forecasts at 00 and 12 Z cycles from 1 July 2019 to 10 August 2019.

Effects on Analysis Field
The soil moisture analysis fields from the EXP and CTL experiments were compared with each other and with those from the Soil Moisture Operational Product System (SMOPS) at NOAA, which is produced by using blended multi-satellite soil moisture data. The correlation coefficients were then determined to assess the skills required for analyzing the soil moisture in the experiments [25,51]. The anomaly values of the soil moisture analysis were computed by subtracting the period averages from the analysis values, according to the various climate regions [52]. In this climate region classification, tropical climates are the zone where the monthly mean temperature of the coldest month is higher than or equal to 18 • C, and temperate climates are the zone where the monthly mean temperature of the coldest month is higher than 3 • C and less than 18 • C. The correlation coefficients in the tropical (0.792) and temperate (0.615) regions between the EXP and SMOPS were higher than those (0.351 and 0.130) between the CTL and SMOPS (Table 2). In particular, the correlation coefficient of the EXP in the temperate region was statistically significant, but the correlation coefficient of the CTL was not.  Figure 5 shows the biases of CTL against the ECMWF fifth reanalysis (ERA5) [53] and the differences in the soil moisture analysis between the EXP and the CTL for the first soil layer and the third soil layer. The first soil layer's depths are 0.1 m and 0.07 m for Noah LSM and for ERA5, respectively, and the third soil layer's depth is 1.0 m (upper root zone). The second and fourth soil layers were similar to the first and third soil layers, respectively. For the first soil layer, the CTL was drier than ERA5, except for arid areas (Figure 5a). The EXP was wetter than the CTL in most regions of Eurasia, South Africa, and America ( Figure 5b). Conversely, the EXP was drier in several arid and high-latitude regions (from 50 • N to 60 • N) in North America (i.e., Western Canada). For the third soil layer, CTL's dry biases were larger than the first soil layer, particularly in the Brazilian highlands (Figure 5c).
The shrinking of the CTL's wet bias areas was seen, but the wet biases were larger at the third soil layer. EXP alleviated CTL's dry biases in East Asia, Eastern Europe, South Africa, the Central United States, and the Brazilian highlands, but the CTL's dry biases were worse in Western Canada in EXP (Figure 5d). These soil moisture differences affected the atmospheric analysis fields in the lower troposphere of the model. Figure 6 compares the 2-m specific humidity biases and root mean squared errors (RMSE) against ERA5 between the EXP and CTL experiments. A paired t-test was used to estimate whether the forecasts or errors from the two experiments were significantly different. This t-test is generally useful for two correlated samples [54]. Overall, a dry bias was dominant in the CTL experiment for 2-m-specific humidity, except for the region from 50 • N to 60 • N in Eurasia and Eastern and Northwest North America (Figure 6a). However, the EXP reduced the CTL's dry biases and RMSEs in the dry biased regions (Figure 6b,d). These regions were identical to those where the EXP was wetter than the CTL for the soil moisture field (Figure 5b). This result was attributed to more moisture being supplied in the EXP from the surface to the lower atmosphere and was particularly prominent in Asia from 30 • N to 50 • N, Africa, and South America from 0 • E to 30 • N. Moreover, this result was statistically significant at the 95% confidence level. In contrast, it was not apparent whether the wet biased regions in the CTL were improved in the EXP (Figure 6d).
The CTL shows warm biases for 2-m air temperatures for the Eurasian continent (predominant over latitudes higher than 30 • N), Africa (except for the Sahara Desert), and South America (0 • E to 30 • S), corresponding to dry biases (Figure 7a). Moreover, the mid-eastern part of North America and the Australian continent were warmer than ERA5. Overall, the warm biased regions were mostly consistent with the dry biased regions. In these areas, which comprised tropical and temperate regions, the EXP was cooler than the CTL, resulting in reduced RMSEs (Figure 7b,d). This is because increased soil moisture strengthens the latent heat flux and weakens the sensible heat flux from the ground to the lower atmosphere [3,12].

Effects on Forecast Field
The prediction of the NWP model at the mid-range period provides valuable information for weather forecasting. Changes in the analysis fields via soil moisture data assimilation affect the forecast field. For the 2-m-specific humidity and temperature, the effects of soil moisture assimilation on a five-day forecast were consistent with the analysis fields but were larger. In the EXP, dry biases at the 2-m-specific humidity and warm biases at the 2-m temperature were reduced by soil moisture assimilation, and the RMSEs also decreased (Figures 8b,d and 9b,d). These improvements were statistically significant. Owing to the memory characteristics of soil moisture, improvements in the dry and warm biases in a near surface could accumulate [10]. The EXP wetter regions in 2-m humidity than the CTL were similar to wetter regions in soil moisture at the third soil layer (Figure 5d) in East Asia, Eastern Europe, South Africa, the Central United States, and the Brazilian highlands. For example, it is noticeable that CTL's large, dry biases of 2-m humidity, indicating large RMSEs, grew from the analysis field (Figure 6a,c) in the Brazilian highlands (Figure 8a,c). These biases might be caused from drier soil moisture conditions from the third soil layer (Figure 5c). Additionally, it was likely that there was feedback between the CTL's dry biases (Figure 8a   In addition to the near-surface variables, we analyzed the vertical profiles of the atmospheric variables to examine how changes in the near surface would be translated into the lower troposphere and above. The CTL indicated dry biases for a specific humidity at the five-day forecast time in the lower and mid-atmospheric levels from 30 • S to 50 • N ( Figure 10a). Moreover, the EXP lessened the dry biases and RMSEs in the regions at the 95% confidence level because of the increased moisture transported from the surface (Figure 10b,d). More moisture from the surface would suppress the sensible heat flux and increase the latent heat flux, which correspond to a reduction in the warm biases (Figure 11b,d). The temperature improvement in the EXP also significantly improved the geopotential height in the lower and mid-levels from 30 • S to 50 • N (not shown). EXP increased the wet biases at 1000 hPa to 925 hPa in the tropical region but the EXP's RMSEs were reduced. This means that the CTL's errors have a larger variation than the EXP.
The overall skill scores of the atmospheric forecasts are displayed in a scorecard in Figure 12 [55], which indicates whether the performance of EXP is better (RMSE reduction; green color) or worse (RMSE increase; purple color) than that of the CTL against the ERA5 reanalysis at each forecast time. Table 3 defines the areas for the scorecards [56]. It is noticeable that the cells representing an improvement in the EXP (indicated by green colors) were dominant in the Northern Hemisphere, Asia, and East Asia at a specific humidity of 700 hPa and wind speed, geopotential height, and temperature of 850 hPa, 500 hPa, and 250 hPa, respectively. Moreover, the improvement rates tended to increase with the forecast time. Among these variables, the temperature at 850 hPa and the geopotential at 500 hPa showed significant improvements. The improvement rates at the five-day forecast were 2.92%, 5.49%, and 5.32% for the temperature at 850 hPa and 2.80%, 3.34%, and 1.18% for the geopotential height at 500 hPa in the Northern Hemisphere, Asia, and East Asia, respectively. However, several degradation cells (mostly > −1%) and neutrality were observed, except for the improvement metrics. ERA5 reanalysis data, though, are useful for verifying the NWP model forecasting skills; ERA5, like other reanalysis datasets, could have uncertainty, especially in the nearsurface variables [57][58][59][60][61]. Thus, we also analyzed the five-day forecasts of the experiments with the surface synoptic weather station network (SYNOP) observation data to ensure a robust verification. The SYNOP datasets are collected by the Global Telecommunication System, and the model forecast results were interpolated horizontally and vertically to the observation space, detailed in Lee et al. [62]. The evaluation results matched well with those against ERA5 in land. For a near-surface-specific humidity, the CTL was drier than the surface observations for most of the land surfaces, except for several parts of the Eurasian continent and North America from 50 • N to 60 • N and the Australian continent (Figure 13a). The EXP significantly reduced the dry biases and RMSEs in the dry bias region of the CTL (Figure 13b,d). Conversely, the EXP produced mixed effects on the wet bias region of the CTL. Moreover, the warm bias in the CTL was accompanied by a dry bias in the CTL for the land from 30 • N to 50 • N in Asia, and from 0 • E to 30 • S in Africa and South America (Figure 14a); the warm bias and RMSE in these regions were reduced in the EXP (Figure 14b,d). In Western Canada, a warm bias was accompanied by a wet bias in the CTL. In this region, the EXP was drier than the CTL, and it made the EXP warmer. Thus, there might be another factor for a warm bias for Western Canada, except for soil moisture in the CTL.  Scorecard of the RMSE improvement rates for the specific humidity at 700 hPa (s700 hPa); wind speed at 850 hPa, 500 hPa, and 250 hPa (w850 hPa, w500 hPa, and w250 hPa); geopotential height at 850 hPa, 500 hPa, and 250 hPa (z850 hPa, z500 hPa, and z250 hPa); and temperature at 850 hPa, 500 hPa, and 250 hPa (t850 hPa, t500 hPa, and t250 hPa) in the Northern Hemisphere (NH), Southern Hemisphere (SH), Tropics (TR), Asia (AS), and East Asia (EA) for a five-day forecast. Table 3. Area definitions for calculating the scorecards (Figure 12).

Areas Definitions
Northern The diurnal cycles of the near-surface specific humidity and air temperature were analyzed at 44.57 • N and 120.9 • E for an approximate five-day forecast ( Figure 15). The CTL was drier than the surface observations ( Figure 15a). However, the EXP did lessen (8.5% at 12Z cycle) the dry bias of the CTL. Moreover, the daytime improvement in the EXP was more significant than that at night. In addition to the bias reduction, the EXP indicated an improved performance in predicting the diurnal cycles of a specific humidity than that of the CTL. The near-surface temperature in the CTL was higher than the observation temperature during the daytime (Figure 15b), and this bias was attributed to excessive sensible heating from the surface in the CTL. Corresponding to a specific humidity bias reduction, the EXP reduced the CTL warm bias (5.0% at the 12Z cycle); however, the temperature bias of the EXP was similar to that of the CTL at night when the cooling effect was dominant.  If the number of sites was less than thirty, the mean value was not calculated. The variation of the near-surface-specific humidity clearly showed dry biases of the CTL and improvement of the EXP in Eurasia at almost the five-day forecast. The model forecasts at 120 h were averaged for the 00 and 12 UTC target times Figures 16a,c and 17a,c), and the ones at 114 h were averaged for 06 and 18 UTC (Figures 16b,d and 17b,d). EXP lessened the dry biases of the CTL and predicted a more similar near-surface-specific humidity to the observation than the CTL (Figure 16). Corresponding to the dry bias, there were the CTL's warm biases in variations of the near-surface temperature according to the longitude for Eurasia ( Figure 17). The EXP decreased the CTL's warm biases in Eurasia. However, the magnitude of the CTL's dry biases in North America were not as large as the ones in Eurasia, and there were wet biases at 00 UTC in Eastern North America ( Figure 16). Additionally, the differences between the EXP and CTL in North America were smaller. The CTL showed warm biases in Western North America and mixture bias patterns in Eastern North America (Figure 17). The EXP and CTL showed similar forecasts for the near-surface air temperatures in this region.

Summary and Discussion
In this study, we constructed a global soil moisture data assimilation system for GDASPS-KIM to improve its weather forecast skill and examined the impact of soil moisture data assimilation on GDAPS-KIM during the boreal summer season. In the soil moisture data assimilation system, the ASCAT soil moisture observations were assimilated every six hours using EnKF to improve the initial land surface conditions of GDAPS-KIM. Two experiments with and without soil moisture data assimilation were conducted for the period of 00Z 1 July 2019 to 12Z 10 August 2019.
For the initial fields, the soil moisture data assimilation (0.792 and 0.615) showed higher temporal correlation coefficients than those (0.351 and 0.130) without it in the tropical and temperate regions (compared to the SMOPS observations). In addition, the soil moisture content in Eurasia, South Africa, and America increased, which affected the low tropospheric initial fields. Compared to the ERA5 reanalysis fields, it statistically significantly reduced the dry biases and RMSEs in Asia from 30 • N to 50 • N and in Africa and South America from 0 • S to 30 • S. This result was attributed to a greater soil moisture supply from the ground, owing to soil moisture data assimilation. Corresponding to the dry biases, the experiment with the soil moisture data assimilation significantly reduced the warm biases and RMSEs of GDAPS-KIM. An augmented moisture content, which was associated with soil moisture data assimilation, would increase the latent heat flux and weaken the sensible heat flux from the surface to the lower atmosphere [12].
The influence of soil moisture data assimilation on a five-day weather forecast was also examined and compared to the ERA5 reanalysis fields. Soil moisture data assimilation decreased the temperature biases and RMSEs in similar regions, thereby positively impacting the initial field. Moreover, the magnitude of the improvement in the five-day forecast was stronger than that in the initial field. The effects of the realistic initial soil moisture conditions were accumulated because of the soil moisture memory characteristics. Specifically, the EXP's reduction for the five-day forecast was larger than the analysis field in the dry and warm errors of the CTL over the near surfaces in East Asia, Eastern Europe, South Africa, the Central United States, and the Brazilian highlands, where the soil moisture conditions were wetter than the CTL, especially in the third soil layer (root-zone). Ma et al. [63] and Jun et al. [64] also showed that the soil moisture was the cause of the 2-m air temperature over the Central United States and Asia for the five-day forecast in the NWP models in the boreal summer. However, employing ground observations would be needed for the robust evaluation. As shown in the zonally mean fields, GDAPS-KIM produced a dry bias and warm bias in the lower and mid-levels from 30 • S to 50 • N. However, the soil moisture data assimilation significantly decreased the biases and RMSEs of the specific humidity and temperature at the 95% confidence level.
It is evident from the scorecard, wherein the GDAPS-KIM forecast performance is indicated with a lead time, that the positive impact of soil moisture data assimilation is dominant in the Northern Hemisphere, Asia, and East Asia. It is remarkable that the positive impact tended to rise, along with the forecast period. This was likely because of soil moisture memory explained above. The soil moisture data assimilation led to higher forecast capabilities for temperatures at 850 hPa and geopotential heights at 500 hPa in the Northern Hemisphere (2.92% and 2.80%, respectively), Asia (5.49% and 3.34%, respectively), and East Asia (5.32% and 1.18%, respectively).
In conclusion, incorporating the soil moisture assimilation contributes to improving the initial surface condition and predictability of GDAPS-KIM during the boreal summertime. In particular, the soil moisture data assimilation alleviated the warm biases that are related to soil moisture and enhanced the weather forecasting. It demonstrated that soil moisture assimilation with satellite observation could improve not only the initial surface soil moisture conditions but the weather forecast skills through land-air interactions, as in previous studies [1,[3][4][5][6][7][8]24,26]. Based on these results, this assimilation system has been operational for the GDAPS-KIM at KMA since October 2020.
In a future work, another satellite product will be employed to further improve the quality of the initial soil moisture conditions. There are several satellite sensor products for soil moisture described above, but we will choose the one that can be constantly provided in near-real time.