Observational evidence for the Monin-Obukhov similarity under all stability conditions

Data collected in the surface layer in a northern suburban area of Nanjing from 15 November to 29 December 2007 were analyzed to examine the Monin-Obukhov similarity for describing the turbulent fluctu- ations of 3D winds under all stability conditions and to obtain the turbulence characteristics under different weather conditions. The results show that the dimensionless standard deviations of turbulent velocity com- ponents (σ u /u* , σ v /u* , σ w /u * ) and dimensionless turbulent kinetic energy (TKE) can be well described by "1/3" power law relationships under stable, neutral, and unstable conditions, with σ u /u * > σ v /u * > σ w /u* . Land use and land cover changes mainly impact dimensionless standard deviations of horizontal component fluctuations, but they have very little on those of the vertical component. The dimensionless standard devi- ations of wind components and dimensionless TKE are remarkably affected by different weather conditions; the deviations of horizontal wind component and dimensionless TKE present fog day > clear sky > overcast > cloudy; the trend of the vertical wind component is the reverse. The surface drag coefficient at a Nan- jing suburban measurement site during the observation period was obviously higher than at other reported plains and plateau areas, and was approximately one order larger in magnitude than the reported plains areas. Dimensionless standard deviation of temperature declined with increasing |z /L| with an approximate "-1/3" slope in unstable stratification and "-2/3" slope in stable stratification.     

Modeling Study of the Global Distribution of Radiative Forcing by Dust Aerosol

To quantitatively understand the dust aerosol effects on climate change,we calculated the global distribution of direct radiative forcing due to dust aerosol under clear and cloudy skies in both winter and summer,by using an improved radiative transfer model and the global distribution of dust mass concentration given by GADS (Global Aerosol Data Set).The results show that the global means of the solar forcing due to dust aerosol at the tropopause for winter and summer are –0.48 and –0.50 W m-2 ,respectivel...     

Characteristics of Solar Radiation and the Impact of Clouds at Yangbajing, Tibet

Yangbajing (YBJ) is located in the Tibetan Plateau, China. The characteristics of solar radiation and its relationship with clouds at YBJ from April 2009 to April 2010 were analyzed in this paper. The annual mean solar radiation was 478.4 W m 2 , and the annual mean transmittance was 0.713. The atmospheric mean trans- mittance of clear skies reaches 0.828 when the solar elevation angle (SEA) is greater than 10 degrees. Comparisons with numerical simulations show that the atmosphere of YBJ is clean. Impacts from atmospheric conditions on solar radiation are similar for clear skies during the year because the standard deviation of transmittance in clear skies was less than 0.05 when the SEA was greater than 10 degrees. It is important to understand the impact of clouds on solar radiation without considering other impact factors. In the last part of this article, the authors analyzed and established a statistical quantitative relationship between surface solar radiation and cloud fraction.     

Aircraft Measurements of Cloud–Aerosol Interaction over East Inner Mongolia

To investigate the potential effects of aerosols on the microphysical properties of warm clouds, airborne observational data collected from 2009 to 2011 in Tongliao, Inner Mongolia, China, were statistically analyzed in this study. The results demonstrated that the vertical distribution of the aerosol number concentration(N_a) was similar to that of the clean rural continent. The average aerosol effective diameter(D_e) was maintained at approximately 0.4 μm at all levels. The data obtained during cloud penetrations showed that there was a progressive increase in the cloud droplet concentration(N_c) and liquid water content(LWC) from outside to inside the clouds, while the Nawas negatively related to the Ncand LWC at the same height. The fluctuation of the N_a, Ncand LWC during cloud penetration was more obvious under polluted conditions(Type 1) than under clean conditions(Type 2). Moreover, the wet scavenging of cloud droplets had a significant impact on the accumulation mode of aerosols, especially on particles with diameters less than 0.4 μm. The minimum wet scavenging coefficient within the cloud was close to 0.02 under Type 1 conditions, while it increased to 0.1 under Type 2 conditions,which proved that the cloud wet scavenging effect under Type 1 conditions was stronger than that under Type 2 conditions.Additionally, cloud droplet spectra under Type 1 conditions were narrower, and their horizontal distributions were more homogeneous than those under Type 2 conditions.     

VARIATION CHARACTERISTICS OF THE PLANETARY BOUNDARY LAYER HEIGHT AND ITS RELATIONSHIP WITH PM2.5 CONCENTRATION OVER CHINA

The planetary boundary layer height (PBLH) was calculated using the radiosonde sounding data, including 120 L-band operational sites and 8 GPS sites in China. The diurnal and seasonal variations of PBLH were analyzed using radiosonde sounding (OBS-PBLH) and ERA data (ERA-PBLH). Based on comparison and error analyses, we discussed the main error sources in these data. The frequency distributions of PBLH variations under different regimes (the convective boundary layer, the neutral residual layer, and the stable boundary layer) can be well fitted by a Gamma distribution and the shape parameter k and scale parameter s values were obtained for different regions of China. The variation characteristics of PBLH were found in summer under these three regimes for different regions. The relationships between PBLH and PM2.5 concentration generally follow a power law under very low or no precipitation conditions in the region of Beijing, Tianjin and Hebei in summer. The results usually deviated from this power distribution only under strong precipitation or high relative humidity conditions because of the effects of hygroscopic growth of aerosols or wet deposition. The OBS-PBLH provided a reasonable spatial distribution relative to ERA-PBLH. This indicates that OBS-PBLH has the potential for identifying the variation of PM2.5 concentration.     

One-Dimensional Variational Retrieval of Temperature and Humidity Profiles from the FY4A GIIRS

A physical retrieval approach based on the one-dimensional variational(1 D-Var) algorithm is applied in this paper to simultaneously retrieve atmospheric temperature and humidity profiles under both clear-sky and partly cloudy conditions from FY-4 A GIIRS(geostationary interferometric infrared sounder) observations. Radiosonde observations from upper-air stations in China and level-2 operational products from the Chinese National Satellite Meteorological Center(NSMC)during the periods from December 2019 to January 2020(winter) and from July 2020 to August 2020(summer) are used to validate the accuracies of the retrieved temperature and humidity profiles. Comparing the 1 D-Var-retrieved profiles to radiosonde data, the accuracy of the temperature retrievals at each vertical level of the troposphere is characterized by a root mean square error(RMSE) within 2 K, except for at the bottom level of the atmosphere under clear conditions. The RMSE increases slightly for the higher atmospheric layers, owing to the lack of temperature sounding channels there.Under partly cloudy conditions, the temperature at each vertical level can be obtained, while the level-2 operational products obtain values only at altitudes above the cloud top. In addition, the accuracy of the retrieved temperature profiles is greatly improved compared with the accuracies of the operational products. For the humidity retrievals, the mean RMSEs in the troposphere in winter and summer are both within 2 g kg^–1. Moreover, the retrievals performed better compared with the ERA5 reanalysis data between 800 h Pa and 300 h Pa both in summer and winter in terms of RMSE.     

An Improved, Downscaled, Fine Model for Simulation of Daily Weather States

In this study,changes in daily weather states were treated as a complex Markov chain process,based on a continuous-time watershed model(soil water assessment tool,SWAT) developed by the Agricultural Research Service at the U.S.Department of Agriculture(USDA-ARS).A finer classification using total cloud amount for dry states was adopted,and dry days were classified into three states:clear,cloudy,and overcast(rain free).Multistate transition models for dry-and wet-day series were constructed to comprehensively downscale the simulation of regional daily climatic states.The results show that the finer,improved,downscaled model overcame the oversimplified treatment of a two-weather state model and is free of the shortcomings of a multistate model that neglects finer classification of dry days(i.e.,finer classification was applied only to wet days).As a result,overall simulation of weather states based on the SWAT greatly improved,and the improvement in simulating daily temperature and radiation was especially significant.     

Differences in Meteorological Conditions between Days with Persistent and Non-Persistent Pollution in Beijing,China

We compared the regional synoptic patterns and local meteorological conditions during persistent and non-persistent pollution events in Beijing using US NCEP–Department of Energy reanalysis outputs and observations from meteorological stations. The analysis focused on the impacts of high-frequency(period 90 days) variations in meteorological conditions on persistent pollution events(those lasting for at least 3 days). Persistent pollution events tended to occur in association with slow-moving weather systems producing stagnant weather conditions, whereas rapidly moving weather systems caused a dramatic change in the local weather conditions so that the pollution event was short-lived. Although Beijing was under the influence of anomalous southerly winds in all four seasons during pollution events, notable differences were identified in the regional patterns of sea-level pressure and local anomalies in relative humidity among persistent pollution events in different seasons. A region of lower pressure was present to the north of Beijing in spring, fall, and winter, whereas regions of lower and higher pressures were observed northwest and southeast of Beijing, respectively, in summer. The relative humidity near Beijing was higher in fall and winter, but lower in spring and summer. These differences may explain the seasonal dependence of the relationship between air pollution and the local meteorological variables. Our analysis showed that the temperature inversion in the lower troposphere played an important part in the occurrence of air pollution under stagnant weather conditions.Some results from this study are based on a limited number of events and thus require validation using more data.     

Albedo of coastal landfast sea ice in Prydz Bay,Antarctica: Observations and parameterization

The snow/sea-ice albedo was measured over coastal landfast sea ice in Prydz Bay, East Antarctica(off Zhongshan Station)during the austral spring and summer of 2010 and 2011. The variation of the observed albedo was a combination of a gradual seasonal transition from spring to summer and abrupt changes resulting from synoptic events, including snowfall, blowing snow, and overcast skies. The measured albedo ranged from 0.94 over thick fresh snow to 0.36 over melting sea ice. It was found that snow thickness was the most important factor influencing the albedo variation, while synoptic events and overcast skies could increase the albedo by about 0.18 and 0.06, respectively. The in-situ measured albedo and related physical parameters(e.g., snow thickness, ice thickness, surface temperature, and air temperature) were then used to evaluate four different snow/ice albedo parameterizations used in a variety of climate models. The parameterized albedos showed substantial discrepancies compared to the observed albedo, particularly during the summer melt period, even though more complex parameterizations yielded more realistic variations than simple ones. A modified parameterization was developed,which further considered synoptic events, cloud cover, and the local landfast sea-ice surface characteristics. The resulting parameterized albedo showed very good agreement with the observed albedo.     

An improved method for correction of air temperature measured using different radiation shields

The variation of air temperature measurement errors using two different radiation shields （DTR502B Vaisala,Finland,and HYTFZ01,Huayun Tongda Satcom,China） was studied.Datasets were collected in the field at the Daxing weather station in Beijing from June 2011 to May 2012.Most air temperature values obtained with these two commonly used radiation shields were lower than the reference records obtained with the new Fiber Reinforced Polymers （FRP） Stevenson screen.In most cases,the air temperature errors when using the two devices were smaller on overcast and rainy days than on sunny days; and smaller when using the imported rather than the Chinese shield.The measured errors changed sharply at sunrise and sunset,and reached maxima at noon.Their diurnal variation characteristics were,naturally,related to changes in solar radiation.The relationships between the record errors,global radiation,and wind speed were nonlinear.An improved correction method was proposed based on the approach described by Nakamura and Mahrt （2005） （NM05）,in which the impact of the solar zenith angle （SZA） on the temperature error is considered and extreme errors due to changes in SZA can be corrected effectively.Measurement errors were reduced significantly after correction by either method for both shields.The error reduction rate using the improved correction method for the Chinese and imported shields were 3.3％ and 40.4％ higher than those using the NM05 method,respectively.     

MODIS BRIGHTNESS TEMPERATURE DATA ASSIMILATION UNDER CLOUDY CONDITIONS: METHODS AND IDEAL TESTS

Clouds have important effects on the infrared radiances transmission in that the inclusion of cloud effects in data assimilation can not only improve the quality of the assimilated atmospheric parameters greatly, but also minimize the initial error of cloud parameters by adjusting part of the infrared radiances data. On the basis of the Grapes-3D-var (Global and Regional Assimilation and Prediction Enhanced System), cloud liquid water, cloud ice water and cloud cover are added as the governing variables in the assimilation. Under the conditions of clear sky, partly cloudy cover and totally cloudy cover, the brightness temperature of 16 MODIS channels are assimilated respectively in ideal tests. Results show that when the simulated background brightness temperatures are lower than the observation, the analyzed field will increase the simulated brightness temperature by increasing its temperature and reducing its moisture, cloud liquid water, cloud ice water, and cloud cover. The simulated brightness temperature can be reduced if adjustment is made in the contrary direction. The adjustment of the temperature and specific humidity under the clear sky conditions conforms well to the design of MODIS channels, but it is weakened for levels under cloud layers. The ideal tests demonstrate that by simultaneously adding both cloud parameters and atmospheric parameters as governing variables during the assimilation of infrared radiances, both the cloud parameters and atmospheric parameters can be adjusted using the observed infrared radiances and conventional meteorological elements to make full use of the infrared observations.     

Changes in Vegetation and Assessment of Meteorological Conditions in Ecologically Fragile Karst Areas

Meteorological conditions have an important impact on changes of vegetation in ecologically fragile karst areas.This study aims to explore a method for quantitative evaluation of these meteorological conditions. We analyzed the changing trend of vegetation during 2000–2018 and the correlations between vegetation changes and various meteorological factors in karst rocky areas of Guangxi Zhuang Autonomous Region, China. Key meteorological factors in vegetation areas with varying degrees of improvement were selected and evaluated at seasonal timescale. A quantitative evaluation model of comprehensive influences of meteorological factors on vegetation was built by using the partial least-square regression(PLS). About 91.45% of the vegetation tended to be improved, while only the rest 8.55% showed a trend of degradation from 2000 to 2018. Areas with evident vegetation improvement were mainly distributed in the middle and northeast, and those with obvious vegetation degradation were scattered. Meteorological factors affecting vegetation were significantly different among the four seasons. Overall, high air humidity, small temperature difference in spring and autumn, and low daily minimum temperature and air pressure were favorable conditions. Low temperature in winter as well as high temperature in summer and autumn were unfavorable conditions. The Climate Vegetation Index(CVI) model was established by PLS using the maximum, minimum, and average temperatures; vapor pressure; rainfall; and air pressure as key meteorological factors. The Enhanced Vegetation Index(EVI) was well fitted by the CVI model, with the average coefficient of determination(r2) and root mean square error(RMSE) of 0.856 and 0.042, respectively. Finally, an assessment model of comprehensive meteorological conditions was built based on the interannual differences in CVI. The meteorological conditions in the study area in 2014 were successfully evaluated by combining the model and selected seasonal key meteorological factors.     

Operational Implementation of the ATOVS Processing Procedure in KMA and Its Validation

The Korea Meteorological Administration (KMA) has processed the data from the advanced TOVS(ATOVS) onboard NOAA-16 satellite since May 2001. The operational production utilizes the AAPP(ATOVS and AVHRR Processing Package) of EUMETSAT and IAPP (International ATOVS ProcessingPackage) of the University of Wisconsin. For the initial guess profiles, the predicted fields (usually 6 to 12hour forecasted fields) from the global aviation model of NOAA/NCEP are used. The average number ofprofiles retrieved from the ATOVS data is about 1,300 for each morning and afternoon orbit at about 18 and06 UTC, respectively. The retrieved temperature and dew point temperatures are provided to forecastersin real time and used for initialization of prediction models. With the advanced microwave sensor (AMSU;Advanced Microwave Sounding Unit), accuracy of the ATOVS products is expected to be better than thatof the TOVS products, especially in cloudy conditions. Indeed, the preliminary results from a validationstudy with the collocated radiosonde data during a 8-month period, from May to December 2001, for theEast Asia region show an improved accuracy of the ATOVS products for cloudy skies versus the TOVS,especially for higher altitudes. The RMS (Root Mean Square) difference between the ATOVS productsand radiosonde data is about 1.3℃ for both clear and cloudy conditions, except for near the ground and athigher altitudes, at around 200 hPa. There is no significant temporal variation of the error statistics at allpressure levels. In case of the water vapor mixing ratio, the largest difference is shown at lower altitudes,while the accuracy is much better for the clear sky cases than the cloudy sky cases. The bias and RMSEat lower altitudes is about 0.557 g kg-1 and 2.5 g kg-1 and decrease significantly with increasing altitude.     

Diurnal and seasonal variations of CO2 fluxes and their climate controlling factors for a subtropical forest in Ningxiang

In this study, the diurnal and seasonal variations of CO2 fluxes in a subtropical mixed evergreen forest in Ningxiang of Hunan Province, part of the East Asian monsoon region, were quantified for the first time. The fluxes were based on eddy covariance measurements from a newly initiated flux tower. The relationship between the CO2 fluxes and climate factors was also analyzed. The results showed that the target ecosystem appeared to be a clear carbon sink in 2013, with integrated net ecosystem CO2exchange(NEE), ecosystem respiration(RE), and gross ecosystem productivity(GEP) of-428.8, 1534.8 and1963.6 g C m-2yr-1, respectively. The net carbon uptake(i.e. the-NEE), RE and GEP showed obvious seasonal variability,and were lower in winter and under drought conditions and higher in the growing season. The minimum NEE occurred on12 June(-7.4 g C m-2d-1), due mainly to strong radiation, adequate moisture, and moderate temperature; while a very low net CO2 uptake occurred in August(9 g C m-2month-1), attributable to extreme summer drought. In addition, the NEE and GEP showed obvious diurnal variability that changed with the seasons. In winter, solar radiation and temperature were the main controlling factors for GEP, while the soil water content and vapor pressure deficit were the controlling factors in summer. Furthermore, the daytime NEE was mainly limited by the water-stress effect under dry and warm atmospheric conditions, rather than by the direct temperature-stress effect.     

Growth and Interactions of Multi-Source Perturbations in Convection-Allowing Ensemble Forecasts

This study investigated the growth of forecast errors stemming from initial conditions (ICs), lateral boundary conditions (LBCs), and model (MO) perturbations, as well as their interactions, by conducting seven 36 h convection- allowing ensemble forecast (CAEF) experiments. Two cases, one with strong-forcing (SF) and the other with weak-forcing (WF), occurred over the Yangtze-Huai River basin (YHRB) in East China, were selected to examine the sources of uncertainties associated with perturbation growth under varying forcing backgrounds and the influence of these back- grounds on growth. The perturbations exhibited distinct characteristics in terms of temporal evolution, spatial propagation, and vertical distribution under different forcing backgrounds, indicating a dependence between perturbation growth and forcing background. A comparison of the perturbation growth in different precipitation areas revealed that IC and LBC perturbations were significantly influenced by the location of precipitation in the SF case, while MO perturbations were more responsive to convection triggering and dominated in the WF case. The vertical distribution of perturbations showed that the sources of uncertainties and the performance of perturbations varied between SF and WF cases, with LBC perturbations displaying notable case dependence. Furthermore, the interactions between perturbations were considered by exploring the added values of different source perturbations. For the SF case, the added values of IC, LBC, and MO perturbations were reflected in different forecast periods and different source uncertainties, suggesting that the combination of multi-source perturbations can yield positive interactions. In the WF case, MO perturbations provided a more accurate estimation of uncertainties downstream of the Dabie Mountain and need to be prioritized in the research on perturbation development.     

Numerical Study on the Impact of Ground Heating and Ambient Wind Speed on Flow Fields in Street Canyons

The impact of ground heating on flow fields in street canyons under different ambient wind speed conditions was studied based on numerical methods.A series of numerical tests were performed,and three factors including height-to-width(H/W) ratio,ambient wind speed and ground heating intensity were taken into account.Three types of street canyon with H/W ratios of 0.5,1.0 and 2.0,respectively,were used in the simulation and seven speed values ranging from 0.0 to 3.0 m s 1 were set for the ambient wind speed.The ground heating intensity,which was defined as the difference between the ground temperature and air temperature,ranged from 10 to 40 K with an increase of 10 K in the tests.The results showed that under calm conditions,ground heating could induce circulation with a wind speed of around 1.0 m s 1,which is enough to disperse pollutants in a street canyon.It was also found that an ambient wind speed threshold may exist for street canyons with a fixed H/W ratio.When ambient wind speed was lower than the threshold identified in this study,the impact of the thermal effect on the flow field was obvious,and there existed a multi-vortex flow pattern in the street canyon.When the ambient wind speed was higher than the threshold,the circulation pattern was basically determined by dynamic effects.The tests on the impact of heating intensity showed that a higher ground heating intensity could strengthen the vortical flow within the street canyon,which would help improve pollutant diffusion capability in street canyons.     

MODIS brightness temperature data assimilation under cloudy conditions Ⅱ: Impacts on rainstorm forecasting

Although satellite observations provide large amount of information of clouds and precipitation and play an important role in the forecast of heavy rainfall, they have not been fully taken advantage of in data assimilation of numerical weather predictions, especially those in infrared channels. Assimilating radiances is common under clear-sky conditions since it is extremely difficult to simulate infrared transmittance in cloudy sky. Based on the Global and Regional Assimilation and Prediction Enhanced System 3-dimensional variance (GRAPES-3DVar), cloud liquid water content, ice-water content and cloud cover are employed as governing variables in the assimilation system. This scheme can improve the simulation of infrared transmittance by a fast radiative transfer model for TOVS (RTTOV) and adjust the atmospheric and cloud parameters based on infrared radiance observations. In this paper, we investigate a heavy rainfall over Guangdong province on May 26, 2007, which is right after the onset of a South China Sea monsoon. In this case, channels of the Moderate Resolution Imaging Spectroradiometer (MODIS) for observing water vapor (Channel 27) and cloud top altitude (Channel 36) are selected for the assimilation. The process of heavy rainfall is simulated by the Weather Research and Forecasting (WRF) model. Our results show that the assimilated MODIS data can improve the distribution of water vapor and temperature in the first guess field and indirectly adjust the upper-level wind field. The tendency of adjustment agrees well with the satellite observations. The assimilation scheme has positive impacts on the short-range forecasting of rainstorm.     

Model Projections of East Asian Summer Climate under the ‘Free Arctic’ Scenario

This paper addresses the ‘ice-free Arctic’ issue under the future global warming scenario. Four coupled climate models used in the third phase of the Coupled Model Intercomparison Project (CMIP3) were selected to project summer climate conditions over East Asia once the Arctic becomes ice-free. The models project that an ice-free Arctic summer will begin in the 2060s under the SRESA1B (according to IPCC Special Reports on Emissions Scenarios) simulations. Our results show that the East Asian summer monsoons will tend to be stronger and that the water vapor transport to central northern China will be strengthened, leading to increased summer precipitation in central northern China. The models also project an intensified Antarctic Oscillation, a condition which favors increased precipitation in South China’s Yangtze River Valley. The overall precipitation in Northwest China is projected to increase under ice-free Arctic summer conditions.     

Airborne Observations of Cloud Condensation Nuclei Spectra and Aerosols over East Inner Mongolia

A set of vertical profiles of aerosol number concentrations, size distributions and cloud condensation nuclei(CCN)spectra was observed using a passive cloud and aerosol spectrometer(PCASP) and cloud condensation nuclei counter, over the Tongliao area, East Inner Mongolia, China. The results showed that the average aerosol number concentration in this region was much lower than that in heavily polluted areas. Monthly average aerosol number concentrations within the boundary layer reached a maximum in May and a minimum in September, and the variations in CCN number concentrations at different supersaturations showed the same trend. The parameters c and k of the empirical function N = c S~kwere 539 and1.477 under clean conditions, and their counterparts under polluted conditions were 1615 and 1.42. Measurements from the airborne probe mounted on a Yun-12(Y12) aircraft, together with Hybrid Single-Particle Lagrangian Integrated Trajectory model backward trajectories indicated that the air mass from the south of Tongliao contained a high concentration of aerosol particles(1000–2500 cm~(-3)) in the middle and lower parts of the troposphere. Moreover, detailed intercomparison of data obtained on two days in 2010 indicated that the activation efficiency in terms of the ratio of NCCNto N_a(aerosols measured from PCASP) was 0.74(0.4 supersaturations) when the air mass mainly came from south of Tongliao, and this value increased to 0.83 on the relatively cleaner day. Thus, long-range transport of anthropogenic pollutants from heavily polluted mega cities,such as Beijing and Tianjin, may result in slightly decreasing activation efficiencies.     

Below-Cloud Aerosol Scavenging by Different-Intensity Rains in Beijing City

Below-cloud aerosol scavenging process by precipitation is important for cleaning the polluted aerosols in the atmosphere, and is also a main process for acid rain formation. However, the related physical mechanism has not been well documented and clarified yet. In this paper, we investigated the below-cloud PM_(2.5)(particulate matter with aerodynamic diameter being 2.5 μm or less) scavenging by different-intensity rains under polluted conditions characterized by high PM_(2.5) concentrations, based on in-situ measurements from March 2014 to July 2016 in Beijing city. It was found that relatively more intense rainfall events were more efficient in removing the polluted aerosols in the atmosphere. The mean PM_(2.5) scavenging ratio and its standard deviation(SD) were 5.1% ± 25.7%, 38.5% ± 29.0%,and 50.6% ± 21.2% for light, moderate, and heavy rain events, respectively. We further found that the key impact factors on below-cloud PM_(2.5) scavenging ratio for light rain events were rain duration and wind speed rather than raindrop size distribution. However, the impacts of rain duration and wind speed on scavenging ratio were not important for moderate and heavy rain events. To our knowledge, this is the first statistical result about the effects of rain intensity, rain duration, and raindrop size distribution on below-cloud scavenging in China.