Abstract. The Gridpoint Statistical Interpolation data assimilation (DA) system was developed for the four-size bin sectional Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) aerosol mechanism in the Weather Research and Forecasting-Chemistry (WRF-Chem) model. The forward and adjoint operators for the aerosol optical depth (AOD) analysis were derived from WRF-Chem aerosol optical code. We applied three-dimensional variational DA to assimilate the multi-wavelength AOD, ambient aerosol scattering coefficient, and aerosol absorption coefficient, measured by the sun-sky photometer, nephelometer, and aethalometer, respectively. These were undertaken during a dust observation field campaign at Kashi in northwestern China in April 2019. The results showed that the DA analyses decreased the low biases in the model aerosols; however, it had had some deficiencies. Assimilating the surface particle concentration increased the coarse particles in the dust episodes, but AOD, and the coefficients for aerosol scattering and absorption, were still lower than observed values. Assimilating aerosol scattering coefficient separately from AOD improved the two optical quantities. However, it caused an overestimation of the particle concentrations at the surface. Assimilating the aerosol absorption coefficient yielded the highest positive bias in the surface particle concentration, aerosol scattering coefficient, and AOD. The positive biases in the DA analysis were caused by the forward operator underestimating particle scattering and absorption efficiency. As a compensation, the DA system increased particle concentrations excessively so as to fit the observed optical values. The best overall improvements were obtained from the simultaneous assimilation of the surface particle concentration and AOD. The assimilation did not substantially change the aerosol chemical fractions. After DA, the clear-sky aerosol radiative forcing at Kashi was −10.5 W m−2 at the top of the atmosphere, which was 46 % higher than the background radiative forcing value.