Fig.1 Average sand emission flux simulated by WRF-Chem in the Gobi Desert during April 14-19， 2020. The box indicates the main area where dust was transported to the Tibetan Plateau （35°-40°N， 92°-102°E） during this dust event. These include the main dust coverage area on the north slope （35°-40°N，92°-96°E） and the main dust coverage area on the east slope （35°-40°N， 99°-102°E）

Fig.2 Comparison of WRF-chem AOD （A）， MODIS_Aqua AOD （C）， and MODIS_Terra AOD （D） during April 14-19， 2020，and comparison between the daily average PM10 concentration and PM10 concentration simulated by WRF-chem at most stations in northern China （30°-50°N， 80°-120°E， B）

Fig.3 PM10 concentration time series for sites and model results during April 14-19， 2020

Fig.4 Spatial distribution of near-surface dust concentration and wind field during April 14-19， 2020

Fig.5 Spatial distribution of dust concentration and wind field at 700 hPa during April 14-19， 2020

Fig.6 Spatial distribution of dust concentration and wind field at 500 hPa during April 14-19， 2020

Fig.7 The meridional section of Gobi Desert dust transport to the Tibet Plateau over the north slope of the Tibetan Plateau on April 16th （A）， 17th （B）， and over the east slope of the Tibetan Plateau on April 15th （C）， 16th （D）， 17th （E）

Fig.8 The 96 h dust HYSPLIT backward trajectories at 00:00 on April 17， 2020 （UTC） of the east （36°N， 100°E， A） and north （37°N， 95°E， B） slopes of the Tibetan Plateau affected by dust this time， and the 120 h dust HYSPLIT forward trajectories of Daranzadegade station （43.57°N， 104.43°E） on April 14th， 2020

Fig.9 Meridional wind speed field at 12:00 （A） and 18:00 （B） on April 16th， and at 00:00 （C） and 06:00（D） on April 17th （blue indicates north wind， red indicates south wind， unit m·s-1， horizontal resolution 0.25°×0.25°. The blue contour line indicates the rising area with a vertical velocity greater than or equal to 0.02 m·s-1， the contour interval is 0.02 m·s-1， the slash is the sinking area， the purple indicates the vertical velocity from -0.04 m·s-1 to -0.02 m·s-1， and the yellow indicates vertical speed from -0.08 m·s-1 to -0.04 m·s-1

Fig.10 The 700 hPa geopotential height （purple line， unit dagpm） and zonal wind speed （blue indicates west wind， red indicates east wind， unit m·s-1， horizontal resolution 0.25°× 0.25°） at 00:00 （A）， 06:00 （B）， 12:00 （C）， 18:00 （D） on April 16th， 15th （E） and 17th （F）

Fig.11 The 500 hPa geopotential height （blue contour， unit is dagpm） and temperature advection （the red area is warm advection， and the blue area is cold advection， unit ℃/3h） at 00:00 （A）， 12:00 on April 15th （B） and 00:00 （C）， 06:00 （D）， 12:00 （E）， 18:00 （F） on April 16th

Fig.12 The 200 hPa jet stream （wind speed in the colored area is greater than 30 m·s-1） and geopotential height （blue dotted contour line， unit dagpm） at 00:00 （A）， 06:00 （B） on April 15th， 06:00 （C）， 18:00 （D） on April 16th

Fig.13 The vertical profiles of the average potential temperature （black lines）， divergence （red indicates the convergence area， unit is s-1）， the flow field （vector arrow） on April 15th （A） and 16th （B）， and the wind speed （the red area is easterly， unit is m·s-1） and the synthetic wind vector on April 15th （E） and 16th （F） over the Hexi Corridor-Qilian Mountains （38°-41°N）； The vertical profiles of the average potential temperature （black lines）， divergence （red indicates the convergence area， unit is s-1）， the flow field （vector arrow） on April 15th （C） and 16th （D）， and the wind speed （the red area is easterly， unit is m·s-1） and the synthetic wind vector on April 15th （G） and 16th （H） over the western Gobi Desert of Mongolia （41°-44°N）