Enhanced Simulation of an Asian Dust Storm by Assimilating GCOM-C Observations

: Dust aerosols have great effects on global and regional climate systems. The Global Change Observation Mission-Climate (GCOM-C), also known as SHIKISAI, which was launched on 23 December 2017 by the Japan Aerospace Exploration Agency (JAXA), is a next-generation Earth observation satellite that is used for climate studies. The Second-Generation Global Imager (SGLI) aboard GCOM-C enables the retrieval of more precious global aerosols. Here, the ﬁrst assimilation study of the aerosol optical thicknesses (AOTs) at 500 nm observed by this new satellite is performed to investigate a severe dust storm in spring over East Asia during 28–31 March 2018. The aerosol observation assimilation system is an integration of the four-dimensional local ensemble transform Kalman ﬁlter (4D-LETKF) and the Spectral Radiation Transport Model for Aerosol Species (SPRINTARS) coupled with the Non-Hydrostatic Icosahedral Atmospheric Model (NICAM). Through veriﬁcation with the independent observations from the Aerosol Robotic Network (AERONET) and the Asian Dust and Aerosol Lidar Observation Network (AD-Net), the results demonstrate that the assimilation of the GCOM-C aerosol observations can signiﬁcantly enhance Asian dust storm simulations. The dust characteristics over the regions without GCOM-C observations are better revealed from assimilating the adjacent observations within the localization length, suggesting the importance of the technical advances in observation and assimilation, which are helpful in clarifying the temporal–spatial structure of Asian dust and which could also improve the forecasting of dust storms, climate prediction models, and aerosol reanalysis. the observed AOTs was symmetric around 1.2, while the probability distribution of the simulated AOTs in the DA experiment was generally consistent with that of the GCOM-C-observed ones; however, the simulated AOTs in the FR experiment peaked at around 0.5. The probability distribution of AOTs biases without GCOM-C data assimilation showed signiﬁcantly negative biases, while the one in the DA experiment with a peak of 25.39% was generally symmetrical and more compressed at the value of 0. Merely 3.71% and 21.40% of the AOTs biases were within ± 0.1 and ± 0.5 in the FR experiment, whereas 20.68% and 59.34% of the biases were achieved within ± 0.1 and ± 0.5 in the DA experiment. These diagnostics of the posterior AOTs in the DA experiment demonstrate that assimilating GCOM-C AOTs can signiﬁcantly improve modelled AOTs, even if the model tends to underestimate the AOTs in that region. AOTs simulated by the two experiments in the four days both show the dust transportation from the Gobi and Taklimakan Deserts to downwind areas, including the North China Plain, Korean Peninsula, and Western Japan. The spatial distributions of dust AOTs on 30 March revealed that dusts emitted from the Gobi Desert moved eastwards towards Eastern China and the North China Plain, and there was evidence of a dust front in the Yellow Sea. On 31 March, the dust storms from the Sahara Desert and Central Asia were transported across Eurasia and probably mixed with the Asian dust storms. During the four days, the daily mean dust AOTs after GCOM-C data assimilation were signiﬁcantly higher than the ones in the FR experiment. Dust AOTs on 29 March around the Gobi Desert were signiﬁcantly strengthened with aerosol data assimilation. Figure shows that the total AOTs simulated by DA experiment were more consistent with the AERONET-observed ones, indicating that GCOM-C aerosol data assimilation contributed to the reproduction of the aerosol processes during this dust event. On 28 March, the total AOTs around the Korean Peninsula in the FR experiment were overestimated compared with AERONET AOTs, which is probably due to the anthropogenic aerosol emission inventory used in this study was developed in 2010 rather than in 2018. Data assimilation corrected these overestimations and drove the simulations closer to observations. From 29 to 31 March, simulated total AOTs in FR experiments over the Beijing site are clearly lower than those observed by the AERONET. AOT underestimations were signiﬁcantly reduced in the DA experiment, indicating that the simulated dust processes in the DA experiment were more authentic than those in the FR experiment. Moreover, the simulated AOTs over South Japan on 30 and 31 March were also signiﬁcantly improved by the GCOM-C AOT data assimilation.


Introduction
Through absorbing and scattering solar and Earth radiation, dust considerably influences the radiative budget of the atmosphere [1][2][3][4][5]. Absorbing aerosols, such as dust, can heat the atmosphere by absorbing the solar radiation [6]. Dust can also serve as the ice nuclei (IN) and the cloud condensation nuclei (CCN) involved in the development of clouds, which can change a cloud's lifetime and albedo [7,8]. In addition, dust is an important factor affecting local air quality [9][10][11], tropospheric chemistry, and phytoplankton growth in the Pacific Ocean [12,13].
The Global Change Observation Mission-Climate (GCOM-C), also known as SHIKI-SAI, is a next-generation Earth-observing satellite for climate studies, which was successfully launched by Japan Aerospace Exploration Agency (JAXA) on 23 December 2017. The

Assimilated GCOM-C Observations
Global aerosol optical thicknesses (AOTs) at 500 nm over land with a spatial resolution of 4 km in SGLI Version 2 Level 2 Aerosol by Non-Polarization (ARNP) [26] are assimilated in this study. All of the products can be downloaded through the JAXA G-Portal (https://gportal.jaxa.jp/, last access: 28 June 2021). Aerosol by SGLI polarimetry is not used here because its retrieval algorithm is currently under development. The accuracy of the monthly mean GCOM-C AOTs over land was estimated to be 0.15, while that over the Remote Sens. 2021, 13, 3020 3 of 19 ocean was 0.1. The retrieval uncertainty σ i of the original GCOM-C AOTs was estimated according to the empirical formulas as follows [27]: 0.166, AOT > 0.6 0.0995, AOT ≤ 0.6 0.0895, AOT ≤ 0.4 0.0770, AOT ≤ 0.3 0.0650, AOT ≤ 0.2 0.0355, AOT ≤ 0.1 (1) In order to ensure the quality of these originally available data, the quality assurance flags (QA_Flags), which include several bit fields for the aerosol discrimination algorithm, were applied to eliminate observations with high uncertainties (Table 1). Quality control procedures for the available data include the following: (1) the coastal flag is an indicator that excludes the coastal pixels with complicated structure; (2) the cloud flag separates "clear" aerosol pixels from "cloudy" ones; (3) we only used pixels with very good AOT confidence flag, indicating high confidence in AOT values; (4) sun glint and stray light flags had to be 0; (5) to decrease cloud shadow influences and large surface reflectance estimation error, the cloud shadow possibility flag and uncertain surface reflectance flag had to be also 0. Then, AOTs were aggregated hourly into a 0.5 × 0.5 • model grid with additional buddy checks to better ensure the spatiotemporal representation of the observations and serve as the assimilated observations [28]. The aggregation procedures that we used were as follows: first, only AOT observations within µ ± σ were used to aggregate in one model grid, where µ represents the mean value of the subgrid observations and σ represents their standard deviation in the aggregation cell; second, the number of subgrid observations within the 0.5 • grid cell had to be more than 20%; third, the observation over the grid cell with the ratio of the AOT standard deviation and mean value higher than 0.5 should be excluded for assimilation [27,29]. This methodology can work on finer domains according to Liu et al. [30]. Gridded total observation uncertainties are calculated as the square root of representative error variances σ 2 r coupled with instrumental error variances σ 2 i . After quality assurance and aggregation, the total number of available hourly GCOM-C AOTs from 28 to 31 March 2018 over land are shown in Figure 1. Over East Asia, most of the GCOM-C AOTs with high confidence were concentrated in Northeast China and the middle part of China. As shown in Figure S1, although the original GCOM-C-observed-AOTs were concentrated in the Gobi Desert, Mongolia, Northeast China, and the middle part of China, the observed-AOTs in the Gobi Desert and Mongolia were mostly eliminated by the strictly quality control procedures. The strictest quality control procedures used in this study eliminated a large number of AOTs observations over bright-reflecting surfaces to ensure the best observations for assimilation. The independent retrievals of optical aerosol properties from the Aerosol Robotic Network (AERONET) and Asian Dust and Aerosol Lidar Observation Network (AD-Net) were used to validate the simulated aerosol optical properties. The locations of the AERONET and AD-Net stations are marked in Figure 1.

Independent
Observations Used for Evaluation 2.2.1. AERONET AERONET (https://aeronet.gsfc.nasa.gov/, last access: 28 June 2021) is a global distributed network of autonomously operated CIMEL Sun photometers measure the sun and sky radiances, which provide a continuously long-term database of retrieved aerosol properties [31,32]. Over the past 25 years, the number of the AERONET sites has increased to more than 700 sites, which cover most parts of the world. AERONET Version 3 Level 2 database includes spectral aerosol optical depth, absorption optical depth, angstrom exponent, single scattering albedo, complex refractive index, and particle size distribution data. Near-real-time retrieved AOTs have a high temporal frequency of approximately 3-15 min for multiple wavelengths (340, 380, 440, 500, 675, 870, 1020, and 1640 nm) [32]. As most of the AERONET sites can provide the AOTs at 440 nm, the AERONET AOTs Remote Sens. 2021, 13, 3020 4 of 19 at 440 nm with the highest data quality levels were directly used to validate the model simulated AOTs at 440 nm. Instantaneous AOTs at 440 nm were averaged within each hour to obtain the hourly mean AOTs for independent validation.

AD-Net
AD-Net (http://www-lidar.nies.go.jp/AD-Net, last access: 28 June 2021) is continuously operated by the National Institute for Environmental Studies (NIES). It is distributed over 20 locations, including in China, Japan, Mongolia, Thailand, and Korea. The twowavelength (532, 1064 nm), polarization-sensitive (532 nm) Mie scattering lidar is the standard lidar system used for AD-Net, which can measure vertical distributions of dust and other aerosols over East Asia for the quantitative understanding of influences of aerosols on the environment, climate, and health [33]. The vertical profiles of the volume depolarization ratio at 532 nm and the attenuated backscatter coefficients at 532 and Remote Sens. 2021, 13, 3020 5 of 19 1064 nm are provided by AD-net. The depolarization ratio and attenuated backscatter coefficients are then used to estimate the extinction coefficients for spherical and non-spherical aerosols [9,10]. The vertical resolution of the lidar observations is 30 m and the temporal resolution is 15 min. AD-net also provides East Asian aerosol components for the Global Atmosphere Watch (GAW) Aerosol Lidar Observation Network [34].

Forward Model
The forward model of the aerosol assimilation system is the global three-dimensional aerosol transport model named the Spectral Radiation Transport Model for Aerosol Species (SPRINTARS) [35][36][37][38], which is coupled online with a new dynamical atmospheric model called the Non-Hydrostatic Icosahedral Atmospheric Model (NICAM) [39][40][41][42], which predicts the spatiotemporal evolution of aerosols. Utilizing the icosahedral grid configuration and the non-hydrostatic dynamic core, NICAM can be run with a subkilometer horizontal resolution [43] and can employ finer meshes in the target regions as the stretched grid system [44]. In this study, the model horizontal resolution was 56 km due to limitations of our computational resources, while the vertical resolution was 38 layers, including 10 layers below 2 km. The aerosol-coupled NICAM considers the main aerosol processes, including emission, advection, diffusion, and deposition. It predicts the mass mixing ratios of the major tropospheric aerosol components, including sea salt, soil dust, sulfate, and carbonaceous aerosols [45,46]. The sulfate and carbonaceous aerosols are represented using the bulk scheme, while the soil dust and sea salt are divided into 10 bins and 4 bins from 0.1 to 10 µm by radius. In this study, the Global Fire Emissions Database (GFEDv3.1) [47] and the Hemispheric Transport of Air Pollution (HTAP) [48] are used as the biomass burning and anthropogenic emission inventories, respectively. Using the referred scheme from [49], dust emission fluxes are calculated online when the three following conditions are satisfied: (1) 10 m wind speeds higher than 6.5 m/s; (2) soil moisture lower than the threshold value; (3) snow cover lower than 1 kg/m 2 [38,50]. Aerosol sink processes are parameterized, including dry deposition, wet deposition, and gravitational settling processes [36]. To reduce the influence of meteorology field errors on the model simulation, the analytical National Centers for Environmental Prediction (NCEP) Final (FNL) data are used to nudge the predicted wind and temperature fields during the simulation.

Data assimilation Methodology
We implemented the local ensemble transform Kalman filter (LETKF) with a fourdimensional extension [25,51,52] to assimilate GCOM-C AOTs observations, which permits the optimization of model-predicted aerosol trajectories over the assimilation time window with the asynchronous observation. The model state was optimized with the LETKF assimilation scheme according to the four formulas below: where x b and x a represent the background-state estimate and analysis-state estimate, respectively; matrix X b represents the background ensemble perturbation; the weight vector w a specifies a linear combination of the background ensemble perturbation to obtain the analysis state estimate in model space by adding to the background state estimate; the vectors y b and y o represent the vector of the background observation mean and the vector of the assimilated observation, respectively; R, H, and I represent the observation error co-variance matrix, observation operator, and identity matrix, respectively; and the ensemble background observation y b(i) is calculated by applying H to the background ensemble as y b(i) = H(x b(i) ). The matrix Y b represents the ensemble background observation perturbation. The analysis ensemble perturbation X a is obtained as the multiple of matrices X b and W a , which continuously depends on analysis error covariance P a . To avoid frequently switching between background ensemble predictions and assimilation and asynchronous assimilating GCOM-C observations, the 4D-LETKF scheme considers approximate model trajectories using linear combinations of the background ensemble trajectories with similar initial conditions over a time interval [51]. Cheng et al. [29] and Dai et al. [25] give the details of this data assimilation system. Three-dimensional aerosol mass mixing ratios in fine and coarse modes were the optimized state vectors used in the data assimilation system [30,50]. The fine mode included sulfate and carbonaceous aerosols, while the coarse mode included the soil dust and sea salt aerosols. As the observables were spectral optical thicknesses that were different from optimized state vectors, the observation operator H was used to convert aerosol mass mixing ratios into AOT space. The simulated AOT at wavelength λ was calculated where M d j is the dry mass concentration of the aerosol in each model layer and β is the prescribed aerosol mass extinction coefficient calculated using Mie scattering theory. Then, the optimal weight, which was determined by the model error covariance matrix and observational error covariance, combined the model-predicted state vector and the observed state vector to produce the updated aerosol mass mixing ratios. Aerosol mass mixing ratios in fine and coarse modes were determined after data assimilation according to their proportions in the background ensemble [46]. Lastly, the updated aerosol mass mixing ratios were used as initial conditions for the next assimilation cycle.
To generate an ensemble simulation with 20 members, the aerosol emissions at every integration step for each model grid were randomly perturbed in a spatially dependent manner to create the model ensemble member. Perturbations for each ensemble member were generated from lognormal distribution, with a mean value of 1 and variance of assumed uncertainty. Sea salt emission was randomly perturbed with 500% uncertainty [25,[53][54][55] and 100% uncertainty was applied to natural dust as in our previous study [29,53]. The emission uncertainties were 130% for sulfate precursors SO 2 and 350% for carbonaceous aerosols [25,55]. Perturbations were assumed to be constant in space and time to generate a sufficient model spread [45]. Using a localization technique, 4D-LETKF was used to explore a much higher-dimensional space and reduce the effects of spurious covariances with distance [56]. The horizontal localization factor f (r) is defined as f (r) = exp −r 2 /2σ 2 with Gaussian features, where r is the distance between the model grid and the position of its assimilated observations in the grid space and σ is the horizontal localization scale [52]. Following Gaspari and Cohn [57], a fifth-order piecewise rational function was used to truncate the infinitely long tails of the Gaussian function. The influences of the observations decreased to zero at 3.65 times the localization scale and observations beyond this distance were not used for data assimilation.
The aerosol-coupled NICAM with the 4D-LETKF was adopted to investigate the aerosol processes and features over East Asia during March 2018, when a strong dust storm event was recorded [13,23]. Two experiments were conducted to reveal the influences of assimilating GCOM-C AOT retrievals on aerosol simulations. One was a simulation without data assimilation, which was used as the reference (FR simulation), while the other was a simulation with assimilating hourly mean AOTs derived from GCOM-C retrievals (DA simulation). The experiments were conducted for 4 days from 00:00 UTC on 28 March to 00:00 UTC on 1 April. The experiments were first initialized from 00:00 UTC on 1 March 2018 using the unmodified emissions without assimilation as a cold start for 22 days, then the experiments were spun up using perturbed emissions with 20-member ensemble assimilation from 00:00 UTC on 23 March 2018 as a warm start for 5 days. According to the 6 h assimilation window, the system performed the ensemble forecast for 6 h with Remote Sens. 2021, 13, 3020 7 of 19 an output every hour. According to Cheng et al. [29], 200 km was set as the horizontal localization scale in this study.

Sanity Check with GCOM-C AOTs
In this section, we first present the simulated AOTs in the FR and DA experiments and their differences during the whole assimilation period (28-31 March 2018) in Figure 2. The simulated AOTs with data assimilation were higher than those without data assimilation over the dust transport path, such as for Tianshan Mountain, the Mongolian Plateau, and Northeast China. The spatial distribution of the AOT increments with data assimilation is generally consistent with that of dust AOTs, indicating the dust dominants the aerosols in these areas during this period; therefore, we focus on these key areas for the studies of the Asian dust storm (shown as a green rectangle in Figure 2). forecast for 6 h with an output every hour. According to Cheng et al. [29], 200 km as the horizontal localization scale in this study.

Sanity Check with GCOM-C AOTs
In this section, we first present the simulated AOTs in the FR and DA expe and their differences during the whole assimilation period (28-31 March 2018) i 2. The simulated AOTs with data assimilation were higher than those witho assimilation over the dust transport path, such as for Tianshan Mountain, the Mo Plateau, and Northeast China. The spatial distribution of the AOT increments w assimilation is generally consistent with that of dust AOTs, indicating the dust do the aerosols in these areas during this period; therefore, we focus on these key a the studies of the Asian dust storm (shown as a green rectangle in Figure 2).  To assess the influences of assimilating GCOM-C AOTs on the modelled fields, Figure 3 shows the biases and root mean square errors (RMSEs) between the simulated and GCOM-C-observed AOTs. Such comparisons are considered to be "sanity checks" to directly determine the effects of aerosol assimilation on the model simulations [29]. Bias and RMSE between the modelled and observed AOTs were reduced from −0.610 and 1.342 to 0.139 and 0.703 with DA. The correlation between simulated and observed field was significantly increased from 0.137 to 0.687 by leveraging new GCOM-C satellite information, indicating the more reasonable horizontal distribution of aerosols. The probability distribution of the observed AOTs was symmetric around 1.2, while the probability distribution of the simulated AOTs in the DA experiment was generally consistent with that of the GCOM-Cobserved ones; however, the simulated AOTs in the FR experiment peaked at around 0.5. The probability distribution of AOTs biases without GCOM-C data assimilation showed significantly negative biases, while the one in the DA experiment with a peak of 25.39% was generally symmetrical and more compressed at the value of 0. Merely 3.71% and 21.40% of the AOTs biases were within ±0.1 and ±0.5 in the FR experiment, whereas 20.68% and 59.34% of the biases were achieved within ±0.1 and ±0.5 in the DA experiment. These diagnostics of the posterior AOTs in the DA experiment demonstrate that assimilating GCOM-C AOTs can significantly improve modelled AOTs, even if the model tends to underestimate the AOTs in that region.
significantly increased from 0.137 to 0.687 by leveraging new GCOM-C satell information, indicating the more reasonable horizontal distribution of aerosols. T probability distribution of the observed AOTs was symmetric around 1.2, while t probability distribution of the simulated AOTs in the DA experiment was genera consistent with that of the GCOM-C-observed ones; however, the simulated AOTs in t FR experiment peaked at around 0.5. The probability distribution of AOTs biases witho GCOM-C data assimilation showed significantly negative biases, while the one in the D experiment with a peak of 25.39% was generally symmetrical and more compressed at t value of 0. Merely 3.71% and 21.40% of the AOTs biases were within ±0.1 and ±0.5 in t FR experiment, whereas 20.68% and 59.34% of the biases were achieved within ±0.1 a ±0.5 in the DA experiment. These diagnostics of the posterior AOTs in the DA experime demonstrate that assimilating GCOM-C AOTs can significantly improve modelled AO even if the model tends to underestimate the AOTs in that region.   Figures 4 and 5 give the spatial distributions of daily mean dust AOTs and total AOTs simulated in the FR and DA experiments, respectively. The daily mean AOTs in "dusty" stations from AERONET are also shown in Figure 5. Here, daily mean dust AOTs represent the daily mean AOTs of dust aerosols, while total AOTs represent the sum of the daily mean AOTs of all dust, sea salt, sulfate, and carbonaceous aerosols. The spatial distributions of AOTs simulated by the two experiments in the four days both show the dust transportation from the Gobi and Taklimakan Deserts to downwind areas, including the North China Plain, Korean Peninsula, and Western Japan. The spatial distributions of dust AOTs on 30 March revealed that dusts emitted from the Gobi Desert moved eastwards towards Eastern China and the North China Plain, and there was evidence of a dust front in the Yellow Sea. On 31 March, the dust storms from the Sahara Desert and Central Asia were transported across Eurasia and probably mixed with the Asian dust storms. During the four days, the daily mean dust AOTs after GCOM-C data assimilation were significantly higher than the ones in the FR experiment. Dust AOTs on 29 March around the Gobi Desert were significantly strengthened with aerosol data assimilation. Figure 5 shows that the total AOTs simulated by DA experiment were more consistent with the AERONET-observed ones, indicating that GCOM-C aerosol data assimilation contributed to the reproduction of the aerosol processes during this dust event. On 28 March, the total AOTs around the Korean Peninsula in the FR experiment were overestimated compared with AERONET AOTs, which is probably due to the anthropogenic aerosol emission inventory used in this study was developed in 2010 rather than in 2018. Data assimilation corrected these overestimations and drove the simulations closer to observations. From 29 to 31 March, simulated total AOTs in FR experiments over the Beijing site are clearly lower than those observed by the AERONET. AOT underestimations were significantly reduced in the DA experiment, indicating that the simulated dust processes in the DA experiment were more authentic than those in the FR experiment. Moreover, the simulated AOTs over South Japan on 30 and 31 March were also significantly improved by the GCOM-C AOT data assimilation.

Verification with Independent AERONET Observations
dust transportation from the Gobi and Taklimakan Deserts to downwind areas, including the North China Plain, Korean Peninsula, and Western Japan. The spatial distributions of dust AOTs on 30 March revealed that dusts emitted from the Gobi Desert moved eastwards towards Eastern China and the North China Plain, and there was evidence of a dust front in the Yellow Sea. On 31 March, the dust storms from the Sahara Desert and Central Asia were transported across Eurasia and probably mixed with the Asian dust storms. During the four days, the daily mean dust AOTs after GCOM-C data assimilation were significantly higher than the ones in the FR experiment. Dust AOTs on 29 March around the Gobi Desert were significantly strengthened with aerosol data assimilation. Figure 5 shows that the total AOTs simulated by DA experiment were more consistent with the AERONET-observed ones, indicating that GCOM-C aerosol data assimilation contributed to the reproduction of the aerosol processes during this dust event. On 28 March, the total AOTs around the Korean Peninsula in the FR experiment were overestimated compared with AERONET AOTs, which is probably due to the anthropogenic aerosol emission inventory used in this study was developed in 2010 rather than in 2018. Data assimilation corrected these overestimations and drove the simulations closer to observations. From 29 to 31 March, simulated total AOTs in FR experiments over the Beijing site are clearly lower than those observed by the AERONET. AOT underestimations were significantly reduced in the DA experiment, indicating that the simulated dust processes in the DA experiment were more authentic than those in the FR experiment. Moreover, the simulated AOTs over South Japan on 30 and 31 March were also significantly improved by the GCOM-C AOT data assimilation.       65° N, 135.59° E) sites are respectively located in the southern and western parts of Japan. The time series of the aerosol extinctions vertical profiles observed for non-spherical particles and the simulated dust are also shown in Figure 8. Comparing the extinction coefficients of total and dust aerosols, it is clear that dust aerosols are the dominant aerosol type in the Ulaanbaatar site. Without the aerosol data assimilation, aerosol extinctions were obviously underestimated during the whole period. Both total aerosol extinctions and the extinctions of dust aerosols in the DA experiment were more consistent with the observed ones, especially over the dust event. From 29 March 12:00 UTC to 30 March 12:00 UTC, the extinctions of dust aerosols observed by lidar and simulated by the DA experiment above 0.5 km were generally up to 0.1 km −1 . During this period, aerosol extinctions simulated by the DA experiment and the observed ones both increased with height and peaked at around 3 km. At the Fukue site, the FR and DA experiments both showed a phenomenon whereby aerosols below 0.5 km were controlled by anthropogenic aerosols (spherical particles) over the whole time period, which was captured by lidar measurements; however, both total and non-spherical particle extinction coefficients in the 2-4 km altitude in the FR experiment on 30-31 March were obviously smaller than the observed ones, indicating that dust aerosols from source regions generally transport above the planetary boundary layer, while the FR experiment underestimated the transported dust aerosols at the Fukue site. The vertical distribution of aerosols, even in FR, should be improved by the model development, as pointed out by Goto et al. [58]. The GCOM-C DA successfully reduced the underestimations, showing comparable vertical patterns of both total and non-spherical particle aerosol extinctions derived from AD-Net.   Figure 8. Comparing the extinction coefficients of total and dust aerosols, it is clear that dust aerosols are the dominant aerosol type in the Ulaanbaatar site. Without the aerosol data assimilation, aerosol extinctions were obviously underestimated during the whole period. Both total aerosol extinctions and the extinctions of dust aerosols in the DA experiment were more consistent with the observed ones, especially over the dust event. From 29 March 12:00 UTC to 30 March 12:00 UTC, the extinctions of dust aerosols observed by lidar and simulated by the DA experiment above 0.5 km were generally up to 0.1 km −1 . During this period, aerosol extinctions simulated by the DA experiment and the observed ones both increased with height and peaked at around 3 km. At the Fukue site, the FR and DA experiments both showed a phenomenon whereby aerosols below 0.5 km were controlled by anthropogenic aerosols (spherical particles) over the whole time period, which was captured by lidar measurements; however, both total and non-spherical particle extinction coefficients in the 2-4 km altitude in the FR experiment on 30-31 March were obviously smaller than the observed ones, indicating that dust aerosols from source regions generally transport above the planetary boundary layer, while the FR experiment underestimated the transported dust aerosols at the Fukue site. The vertical distribution of aerosols, even in FR, should be improved by the model development, as pointed out by Goto et al. [58]. The GCOM-C DA successfully reduced the underestimations, showing comparable vertical patterns of both total and non-spherical particle aerosol extinctions derived from AD-Net. Over the Osaka site, temporal variations of the vertical distributions of the aerosol extinctions in the DA experiment were more comparable to the lidar observations, especially from 16:00 (UTC) on 30 March to 12:00 (UTC) on 31 March, whereas the aerosol extinctions in the FR experiment were evidently lower than the observed ones due to the underestimations of the dust aerosols.

Discussion
To generalize this study, two additional experiments with the same configurations

Discussion
To generalize this study, two additional experiments with the same configurations as before were further conducted from 1 to 3 April 2018 to investigate the effects of GCOM-C AOT assimilation on simulated aerosol optical properties. Figure S2 shows the simulated AOTs in the FR and DA experiments and their differences from 1 to 3 April 2018. Both simulated total AOTs and dust AOTs in the DA experiment were higher than the ones in the FR experiment over northern Xinjiang province, the Mongolian Plateau, and Northeast China. The significantly negative differences of the simulated total AOTs between the DA and FR experiments in the Tarim Basin were dominated by reductions in dust AOTs. Although the simulated AOTs in the DA experiment were significantly lower than the ones in the FR experiment around the middle part of China, Korea, and the Sea of Japan, the simulated dust AOTs were slightly increased with GCOM-C AOT assimilation.
To quantitatively evaluate the influences of GCOM-C AOT assimilations on aerosol simulations, we first compared AOTs simulated by the FR and DA experiments with the AERONET-retrieved ones over the Fukuoka and Osaka sites in Figure 9. RMSEs over the Fukuoka and Osaka sites were reduced from 0.255 and 0.266 to 0.178 and 0.251, while CORRs over these two sites were increased from −0.120 and 0.055 to 0.567 and 0.605. On 2 April, simulated AOTs before assimilation were obviously lower than the observed AOTs over the Osaka site, whereas the ones after assimilation generally agreed with the AERONET observations. FOR PEER REVIEW 13 of 18 the ones in the FR experiment over northern Xinjiang province, the Mongolian Plateau, and Northeast China. The significantly negative differences of the simulated total AOTs between the DA and FR experiments in the Tarim Basin were dominated by reductions in dust AOTs. Although the simulated AOTs in the DA experiment were significantly lower than the ones in the FR experiment around the middle part of China, Korea, and the Sea of Japan, the simulated dust AOTs were slightly increased with GCOM-C AOT assimilation.
To quantitatively evaluate the influences of GCOM-C AOT assimilations on aerosol simulations, we first compared AOTs simulated by the FR and DA experiments with the AERONET-retrieved ones over the Fukuoka and Osaka sites in Figure 9. RMSEs over the Fukuoka and Osaka sites were reduced from 0.255 and 0.266 to 0.178 and 0.251, while CORRs over these two sites were increased from −0.120 and 0.055 to 0.567 and 0.605. On 2 April, simulated AOTs before assimilation were obviously lower than the observed AOTs over the Osaka site, whereas the ones after assimilation generally agreed with the AERONET observations.  [59,60] in the two CALIPSO orbit paths. One orbit path passed the west of Tarim Basin, while the other passed the downwind areas. According to CALIOP aerosol subtypes, the profiles of the aerosol extinctions in the two paths were both dominated by dust aerosols. Figure 10 shows that the increased AOTs in the DA experiment were probably due to increments of the eastward-transported aerosols. Although the FR experiment reproduced the vertical distributions of the aerosols, aerosol extinctions were underestimated. Assimilating GCOM-C AOTs improved the magnitude of the aerosol extinctions. As shown in Figure  11, the FR experiment tended to underestimate the aerosol extinctions, while the DA experiment showed more reasonable magnitude and vertical distributions.  [59,60] in the two CALIPSO orbit paths. One orbit path passed the west of Tarim Basin, while the other passed the downwind areas. According to CALIOP aerosol subtypes, the profiles of the aerosol extinctions in the two paths were both dominated by dust aerosols. Figure 10 shows that the increased AOTs in the DA experiment were probably due to increments of the eastward-transported aerosols. Although the FR experiment reproduced the vertical distributions of the aerosols, aerosol extinctions were underestimated. Assimilating GCOM-C AOTs improved the magnitude of the aerosol extinctions. As shown in Figure 11, the FR experiment tended to underestimate the aerosol extinctions, while the DA experiment showed more reasonable magnitude and vertical distributions.

Summary and Conclusions
As a new next-generation Earth-observing satellite, the Global Change Observation Mission-Climate (GCOM-C) satellite enables us to retrieve more precious optical properties of global aerosols, which is important for further aerosol research. Here, we showed the first utilization of GCOM-C AOTs in aerosol data assimilation. To evaluate its influence on the simulations of an Asian dust event from 28 to 31 March 2018, the aerosol-coupled NICAM with the four-dimensional local ensemble transform Kalman filter (4D-LETKF) was applied to assimilate the new aerosol observations acquired by the GCOM-C satellite.
Dust aerosols were the dominant aerosol type around Northeast China, the Mongolian Plateau, and Tianshan Mountain during this period. Compared with assimilated GCOM-C AOT observations, the FR experiment obviously underestimated the magnitude of this dust storm. Data assimilation with hourly GCOM-C AOT data significantly influenced the model simulation. The DA experiment generally reproduced spatial and temporal dust characteristics that were not captured without data assimilation. Merely 3.71% and 21.40% of the biases were within ±0.1 and ±0.5 in the FR experiment, whereas 20.68% and 59.34% of the biases were achieved within ±0.1 and ±0.5 in the DA experiment. The correlation between the simulated and observed fields was significantly increased from 0.137 to 0.687 by leveraging new GCOM-C satellite information.
Verification of the simulated results by comparison with independent ground-based measurements (AERONET and AD-Net) further confirmed the superiority of the data assimilation. The RMSE (bias) over the Beijing site was reduced from 0.744 (−0.613) to 0.466 (−0.085) by GCOM-C AOT assimilation. Significantly reduced differences between simulated and AERONET-retrieved AOTs indicated that the transport path of dust aerosols in the DA experiment was more authentic than that in the FR experiment. Simulated aerosol vertical extinctions were also improved by GCOM-C AOT assimilation, although the AOT assimilation included limited aerosol vertical information in correcting the aerosol profile. This illustrates that GCOM-C observations provide useful aerosol information to enhance model performances in dust simulation through data assimilation.
As a technical advance that leverages new satellite information, this study is promising for future aerosol simulations; however, the GCOM-C aerosol retrieval algorithm is still in development. We aim to further compare its product with that of other polar satellites and evaluate the possible benefits of data assimilation for forecasting dust storms.