Fog droplet size distributions (DSDs) evolve under the influence of many physical processes, yet their development through the fog lifecycle remains insufficiently understood and challenging to represent in numerical models, constraining the accuracy of fog forecasting. To improve understanding of the fog evolution, field observations under a polluted background were conducted during winters from 2006–2009 and 2017–2018 in Nanjing, China. Among the 27 observed fog events, microphysical properties including fog droplet number concentration (Nf), liquid water content (LWC), volume-mean radius (Rv), and effective radius (Reff) varied substantially. Unimodal, bimodal, and trimodal DSDs were observed, with mode separating diameters of 2 µm for unimodal; 2 and 6–18 µm for bimodal; and 2, 6–12, and 18–26 µm for trimodal DSDs. Both the number of modes and the mode separating diameters vary over the fog life cycle, with more frequent and pronounced changes occurring during fog formation and dissipation or during periods of strong fluctuations in Nf and LWC. Compared with unimodal DSDs, bimodal and trimodal DSDs exhibited broader PDF distributions of LWC, Rv and Reff. Based on these observational features, segmented gamma and lognormal fits were applied to mean DSDs using partition points at 10 and 20 µm. Comparisons between microphysical parameters derived from the fitted DSD and observations show that three-segment fitting improved estimates of Nf and LWC, while substantially enhanced the representation of Reff, absorption coefficient, and optical thickness, reducing deviations from up to 90 % to within 20 %.
National Natural Science Foundation of China420611340094257508541975176Introduction
Composed of suspended small water droplets or individual ice crystals in the air near the surface, fog has multiple impacts ranging from transportation, vegetation, air quality, human health and economy (Gultepe et al., 2014; Jia et al., 2019; Lakra and Avishek, 2022). Due the sharp decline in visibility, long duration, wide spatial coverage associated with fog, it is essential and urgent to improve fog modelling and forecasting. The formation and evolution of fog is driven by macro and microphysical processes including precipitation, radiation, advection, cloud-base lowering, turbulent, aerosol activation and condensation (Gultepe et al., 2007; Mazoyer et al., 2022; Shao et al., 2023; Wang et al., 2020). With different physical processes interact with each other nonlinearly, fog remains as a challenging problem for numerical weather prediction (NWP), even though progress have been made in recent years (Boutle et al., 2018; Martinet et al., 2020; Tudor, 2010).
In order to gain a better understanding in mechanisms of fog evolution, in situ observations have been conducted worldwide with different regions and aerosol backgrounds (Elias et al., 2009; Gultepe et al., 2007, 2009; Gultepe and Milbrandt, 2007; Haeffelin et al., 2010; Mazoyer et al., 2022). Marked variabilities of microphysical parameters such as fog number concentration (Nf) and liquid water content (LWC) have been found. Fogs in polluted area show higher Nf due to higher aerosol number concentration (Na), while in relative clean regions including mountains, rainforests and rural areas there are more big droplets, which contribute to LWC significantly (Gultepe and Milbrandt, 2010; Guo et al., 2015; Li et al., 2017; Nelli et al., 2024). Also, the fierce competition for water vapor associated with higher Na suppresses the condensation growth in urban areas, resulting a lower Nf of large droplets and smaller dispersion in urban fog compared to clean regions (Ge et al., 2024).
Meteorological variables and large-scale processes strongly influence fog formation. In general, fog can be classified into two categories: airmass fog and frontal fog, which can be further divided into cold-advection and warm-advection fog, radiation fog, and sea fog etc. (Willett, 1928). Another fog-forming mechanism is the overall lowering of a cloud layer, including its cloud top (Koračin et al., 2014). Radiation fog typically forms near the surface under clear skies and weak winds associated with anticyclonic conditions. Its primary mechanism is radiative cooling, while opposing effects include upward soil heat flux and the warming and moisture loss caused by turbulent mixing within the stable boundary layer (Brown, 1980; Roach, 1976; Turton and Brown, 1987). The advection fog is associated with the advection of a moist air mass with a temperature contrast relative to the underlying surface, which is mainly coastal but also can be observed over land (Friedlein, 2004). Advection-radiation fog is produced by the radiative cooling of moist air that has been advected inland from the ocean or another large water body (Ryznar, 1977). The C-FOG field campaign along the Atlantic Canada and northeastern US coastlines showed that coastal fog was influenced by multiple weather systems, including northeastern high pressure, west–northwest low pressure, and tropical cyclonic activity (Gultepe et al., 2021). Another study in the region found that fog associated with cyclonic systems was consistently produced by cloud-base lowering and subsequent downward extension to the surface, whereas anticyclonic fog developed either from surface radiative cooling or from the downward extension of low-level stratus to the surface (Dorman et al., 2021).
The fog droplet size distribution is a key characteristic of fog microphysical processes (Niu et al., 2012), which is influenced by aerosol chemical composition and number concentration as well as various environment factors such as temperature, humidity, wind speed and direction (Mazoyer et al., 2017; Price, 2019). Fog droplet size distributions (DSDs) often exhibit one or more distinct modes, referred to as unimodal, bimodal, or trimodal DSD, and can be attributed to different origins of the fog and processes within it (Elias et al., 2015; Hammer et al., 2014; Sampurno Bruijnzeel et al., 2006). Kunkel (1982) finds various shapes in DSDs measured in advection fogs. Many other studies have shown the existence of bimodal DSDs in mature radiation fogs (Meyer et al., 1980; Pinnick et al., 1978; Roach et al., 1976; Wendisch et al., 1998). Gultepe and Milbrandt (2007) reported DSD modes near 4 and 23 µm during winter fog events in the Toronto region. Boudala et al. (2022) investigated the seasonal and microphysical characteristics of fog at Cold Lake Airport in northern Alberta, Canada, and found that radiation fog exhibited a bimodal droplet spectrum with peaks at 4 and 17–25 µm. Mazoyer et al. (2022) observed that fog DSD in a semi-urban area of Paris exhibited both single mode (about 11 µm) and double mode (about 11 and 22 µm). When the fog DSD is bimodal, there is a mass transfer from smaller droplets to larger droplets which may due to collision-coalescence process, while sedimentation by gravity speeds up the removal of fog droplets (Mazoyer et al., 2022). The initial fog DSD is influenced by environmental supersaturation and background aerosol properties. As visibility decreases and fog develops, the DSD broadens, transitioning from unimodal to multimodal (Mazoyer et al., 2022). The characteristics of the DSD also strongly influence the optical properties of fog. Stewart and Essenwanger (1982) showed that the attenuation of electromagnetic radiation by fog depends sensitively on the shape of the droplet size distribution. The DSD and water vapor determine the overall optical properties of fog and its effects on visibility and radiative transfer together.
Integrated with in situ measurement, numerical experiment is a commonly used approach to gain a better understanding of the physical mechanism in fog. Over the past decade, numerous numerical experiments have been conducted to evaluate the fog forecasting capabilities and limitations of various mesoscale NWP models, leading to notable progress (Cui et al., 2019; Payra and Mohan, 2014). Despite WRF has made progress in forecasting certain variables such as temperature and wind, it often struggles to capture the accurate fog lifecycle (Peterka et al., 2024; Román-Cascón et al., 2016). The simulated evolution of fog exhibits a sensitivity to the shape of the DSD comparable to its sensitivity to aerosol loading or cloud droplet number concentration (CDNC), yet it remains one of the least investigated and rarely adjusted components of microphysical parameterization schemes (Boutle et al., 2022). A simulation of a heavy fog event in North China Plain found that effective radius of fog droplet decreases nonlinearly with aerosol number concentration (Jia et al., 2019). Since the effective radius was obtained under the assumption of a monodisperse DSD and the dispersion effect was neglected, it may have been overestimated or underestimated due to fog-aerosol interactions (Chen et al., 2016; Liu and Daum, 2002). Currently, fog DSDs are described using various spectral distribution functions such as exponential, gamma or lognormal functions in bulk and bin microphysical scheme (Kessler, 1969; Khain et al., 2015). However, the fog DSD exhibits strong spatial and temporal variability and evolves throughout the fog lifecycle, often displaying distinct features and deviating from idealized distributions due to turbulent mixing, radiative effects, and gravitational settling (Gultepe et al., 2007; Nelli et al., 2024; Tampieri and Tomasi, 1976; Wang et al., 2021). Such variabilities in DSD could cause substantial deviations from the predefined spectral distribution functions, further bringing challenges for fog parameterization (Khain et al., 2015; Lakra and Avishek, 2022). Therefore, a more physically consistent and adaptable representation of DSD is required to improve simulation reliability of fog evolvement.
In this study, based on the observation data of the 27 fog events obtained in Nanjing, China during the winters of 2006–2009, 2017–2018, we focus on the characteristics and evolution of DSDs, how they are associated with microphysical characteristics and how to improve the representation of multimodal size distributions using the gamma and lognormal function. The rest of the article is organized as follows. Section 2 describes the observation site, data, and methods used in this study. Section 3 presents the results, including an overview of microphysics across 27 fog events, an analysis the fog lifecycle under different modes, the correlations between microphysical characteristics and varying DSD modes, and a refinement of the gamma and lognormal fitting approach with an evaluation of its performance. The main conclusions are presented in Sect. 4.
Data set and methods
The field campaign was conducted during the winters of 2006–2009 and 2017–2018, with each campaign lasting approximately one month per year. The sampling site was located in the northwestern suburban area of Nanjing, Jiangsu Province, China (32.2° N, 118.7° E; 22 m above sea level), north of the Yangtze River and surrounded by industrial facilities, residential areas, and major roads (Niu et al., 2010, 2012). The DSD was measured with a fog monitor (FM-100) from Droplet Measurement Technologies (DMT, USA) with diameters ranging from 1 to 50 µm into 20 bins, at a sampling frequency of 1 Hz. The width of each bin is 2 µm for the first 10 bins and 3 µm for the last 10 bins. To exclude the influence of large unactivated aerosol particles, data from the first bin (1–2 µm) are omitted (Lu et al., 2013). Fog with Nf> 10 cm−3 and LWC > 10−3 g m−3 was identified (Lu et al., 2020; Wang et al., 2021). Microphysical characteristics including fog number concentration (Nf), liquid water content (LWC), volume-mean radius (Rv), effective radius (Reff), relative dispersion (ε), autoconversion threshold (T) and first bin strength (FBS) were calculated through following formulas:
1Nf=∑nr2LWC=1×10-6ρ∑4π3r3n(r)3Rv=∑nrr3Nf134Reff=∑n(r)r3∑n(r)r25σ=nrr-r‾2Nf126ε=σr‾7T=PP0=∫rc∞r6nrdr∫0∞r6nrdr∫rc∞r3nrdr∫0∞r3nrdr8FBS=N1stNf
where ρ is the density of water, r is the fog droplet radius of each bin, r‾ is the mean arithmetic radius defined with r‾=∑nrrNf, rc is defined with rc≈4.09×10-4βcon16N16LWC13, in which βcon=1.15×1023, N1st is the number concentration of the first bin (2–4 µm) following the exclusion of the 1–2 µm bin.
Because visibility observations were unavailable for fog cases 1–3 and 18–20, visibility for these cases was estimated using the observed fog droplet spectra and the microphysical parameterization scheme developed by Gultepe et al. (2006), based on the extinction theory of visible light in fog:
V=-lnαβext
in which
βext=π∑Qextnrr2dr
where α is the constant threshold, typically set to 0.02, βext represents the extinction coefficient, Qext is the Mie extinction efficiency, which depends on particle radius, number concentration, and the wavelength of visible light. When droplet size exceeds about 4 µm, Qext approaches a constant value of 2. For smaller droplets (less than 4 µm), Qext varies between 0.9 and 3.8 (Brenguier et al., 2000; Koenig, 1971).
Previous studies have shown that DSDs can exhibit one or multiple modes (Elias et al., 2015; Hammer et al., 2014; Sampurno Bruijnzeel et al., 2006). Therefore, in this study, two approaches were employed to identify the presence of multimodal DSDs and to determine the corresponding number of modes and their mode separating diameters: (1) identifying turning points where the number concentration transitions from decreasing to increasing with droplet size, and (2) determining modes using unimodal and multimodal gamma and lognormal functions. The detailed procedures are described below.
Identify turning points
For a DSD defined on discrete diameter bins, the change rate of number concentration was first computed as
Δni(D)=ni(D)-ni+1(D)ni(D),i=1,2,…,19
where D is the droplet diameter, nD is the number concentration for each bin. The sign of Δni(D) characterized the local trend of the DSD: a positive value indicates decreasing concentration with increasing diameter, whereas a negative value indicates increasing concentration with increasing diameter. The turning diameter Dturn was identified as the diameter where the sign of Δni(D) change from positive to negative, with the additional requirement that the number concentrations on its right side of are nonzero. In addition, if the number concentration in the first bin is the highest in the DSD, the Dturn of the first mode is assigned a value of 2 µm.
Unimodal and multimodal gamma and lognormal distribution
To determine the number of modes and mode separating positions, each fog DSD was sequentially fitted with unimodal (i=1), bimodal (i=2), and trimodal (i=3) gamma and lognormal distributions as show below:
For the gamma distribution:
nD=∑i=13niD=∑i=13N0,iDμie-λiD
where D is the droplet diameter, nD is the number concentration for each bin, N0, μ and λ are the intercept, shape and slope parameters, respectively. For the lognormal distribution:
nD=∑i=13niD=∑i=13N0,i2π1/2Dilnσg,iexp-(lnDi-lnDg,i)22(lnσg,i)2
where D is the droplet diameter, nD is the number concentration for each bin, N0 is the total number concentration, Dg is the geometric mean diameter and σg is the geometric standard deviation. Details of the upper and lower bounds used in the fitting of Eqs. (12) and (13) are provided in Appendix A1.
For each DSD, unimodal, bimodal, and trimodal gamma and lognormal fits were performed, yielding the corresponding sets of fitting parameters as well as the individual modal components. A fit was retained only if adjacent modal components intersected (mode 1–2 for bimodal fits; mode 1–2 and mode 2–3 for trimodal fits). The turning diameter Dturn was identified as the diameter of the first intersection of two adjacent modes, indicating that the mode begins at the next bin. If the first size bin of the DSD had the highest number concentration, it was identified as Dturn of the first mode. A fit was accepted only if all identified Dturns did not fall within the same or adjacent size bins. Finally, the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC) are used to evaluate the performance of the accepted fits.
AIC and BIC are two influential and widely used model selection criteria in machine learning, engineering, and related scientific fields (Akaike, 1974; Schwarz, 1978; Zhang et al., 2023). AIC provides a numerical basis for ranking competing models by their information loss in approximating the unknown true process, with the model yielding the lowest AIC considered the best approximating model (Symonds and Moussalli, 2011). BIC is consistent in the sense that it selects the true model with probability approaching one. A lower BIC corresponds to a higher posterior probability for the model and is therefore regarded as indicating a better model (Chakrabarti and Ghosh, 2011). The corresponding formulas are as follows:
14AIC=-2lnL^+2p15BIC=-2lnL^+pln(N)
where L^ is the maximum likelihood estimate, p is the number of independently adjusted parameters, N is the number of samples. Among the accepted fits, the one with the smallest combined AIC and BIC is used to determine the number of modes and the corresponding Dturn of the DSD.
Compared with the original unimodal (i=1) fit, the multimodal composite distributions yield lower AIC and BIC values (Fig. 1a). In addition, Nf, LWC, Rv and Reff retrieved from the optimal fit are closer to the observations (Fig. 1b). These results further confirm the validity of this approach and the necessity of representing DSDs with multimodal distributions in this study. A comparison of the optimal fit results based on gamma and lognormal distributions reveals that the lognormal distribution yield lower AIC and BIC values (Fig. 1a) as well as smaller absolute deviations of Nf, LWC, Rv and Reff (Fig. 1b). Therefore, in this study, the lognormal fitting results are more suitable than those based on the gamma distribution for determining the number of modes and the Dturn.
Distributions of AIC/BIC values for the unimodal and optimal gamma and lognormal distribution fits (a), and their absolute deviations between the retrieved and observed microphysical properties (b).
The Dturn identified from the turning point are mainly located at approximately 2, 6–14, and 16–20 µm, whereas those determined from lognormal distributions are primarily found at approximately 2, 6–10, and 20–23 µm. Overall, the two methods exhibit similar Dturn distribution characteristics. In contrast, the Dturn distribution obtained from the gamma distribution lacks a mode near 20 µm. This indicating that in this study, the gamma distribution may fail to capture modes larger than 20 µm. In addition, considering that the observational data may be affected by sampling uncertainty or small fluctuations in number concentration across size bins, the lognormal fitting results were adopted in this study to determine the number of modes and Dturn of the DSD. Possible modes in the 30–50 µm range are excluded when Dturn derived from lognormal fitting was adopted. The sensitivity of the gamma and lognormal representations of the DSD to large droplets is discussed in Appendix A6, with the corresponding results shown in Fig. A20.
Distribution of Dturn derived from the turning-point method (a), gamma fitting (b), and lognormal fitting (c).
Given that temporal resolution may influence the identification of DSD mode numbers and Dturns, we derived these quantities using both 1 and 5 min averaged DSDs. The results from the lognormal distribution are presented here, while those from the gamma distribution are provided in the Fig. A11 in the Appendix. As the two resolutions produced only minor differences (Fig. 3), a 5 min averaging interval was used in this work to reduce the influence of noise.
The Dturn distributions of all DSDs obtained using lognormal distribution at 1 min resolution (a), 5 min resolution (b), and the Dturn distributions by DSD type at both resolutions (c).
ResultsOverview of fog microphysics
The fog type, meteorological variables (temperature, wind speed), visibility, and summary statistics of Nf, LWC, Rv, Reff, ε, T and FBS for each event are provided in Table A1 of the Appendix. Among the 27 fog events, 14 were classified as radiation fog, 8 as radiation-advection fog, 4 as advection fog, while 1 event was likely associated with raindrop evaporation as it occurred after precipitation. Compared to other fog types, radiation-advection fog typically exhibits longer duration, slightly higher wind speed and temperature. The average Nf, LWC, Rv and Reff vary over the ranges of 25–586 cm−3, 0–0.27 g m−3, 1.6–6 µm, 1.9–8.2 µm, respectively, which shows greater Nf, lower LWC and smaller droplet sizes comparing to semi-urban area in Paris, France (Mazoyer et al., 2022) and rainforest area in Xishuangbanna, China (Wang et al., 2021). In the meantime, significant variability in the microphysical properties is observed between different events. Among all DSDs, unimodal, bimodal, and trimodal distributions account for 5 %, 49 %, and 46 %, respectively. Unimodal DSDs exhibit the highest mean Nf, LWC and T, whereas bimodal DSDs have the largest Rv, Reff and relative dispersion (ε), indicating larger droplet sizes and a broader size distribution. Trimodal DSDs show the highest FBS value.
Microphysical characteristics of single mode, double mode and triple mode DSDs, the first row shows the mean values, in the parentheses are the 25th and 75th percentiles of each characteristic.
Figure 4 presents the mean DSDs for different mode types averaged over all fog cases (panel a) and over cases excluding Fog Case 1–2 (panel b). The mean DSD spectrums are sensitive to cases with high Nf and long durations, such as Fog Case 1–2 (F1–2). In particular, 50 % of the unimodal DSDs originate from F1–2, leading to higher Nf in the mean unimodal DSD spectrum than in other DSD types when all cases are included. For all cases, the mean spectrum of unimodal DSDs exhibits the highest Nf, particularly for droplets with diameters below 10 µm. In contrast, bimodal and trimodal DSDs show Nf at diameters below 10 µm that are very close to those of the mean spectrum of all DSDs. As aerosol activation is governed by environmental supersaturation and aerosol hygroscopic properties (Shen et al., 2018; Wang et al., 2019), the similar Nf at small diameters may indicate that aerosol activation persists throughout fog development, despite changes in the number of DSD modes. Within the diameter range of 10–50 µm, bimodal DSDs exhibit higher Nf than trimodal DSDs, consistent with their larger Rv, Reff and ε.
Average spectrums of fog DSDs for different mode types over all fog cases (a), and average spectra of fog DSDs with different modes excluding Fog Case 1–2 (b).
In the following section, five representative cases are selected to examine the lifecycle characteristics of fog events with different modes and the evolution of DSDs throughout the stages of fog formation, development, and dissipation.
Mode transitions and possible mechanisms
Figures 5–9 presents the temporal evolution of meteorological variables, visibility, microphysical properties, and droplet spectrum for the five selected cases. Surface, 850, and 700 hPa synoptic conditions preceding each fog event are provided in the Appendix (Figs. A13–A17). MODIS 3.9 µm shortwave infrared and visible channel imagery for Fog Case 1–1 (F1–1), Fog Case 10 (F10), and Fog Case 20 (F20) are provided in the Appendix (Figs. A8–A10). For Fog Case 4 (F4) and Fog Case 18 (F18), no satellite imagery is available within the fog period because the overpass times of the polar-orbiting satellite did not coincide with the observations. The 1 s temporal evolution of the DSDs are shown in Fig. A18. Among the five fog cases, F4 and F20 have relatively short durations and distinct formation-dissipation lifecycle. Both F1–1 and F18 experience rapid growth, whereas F10 exhibits pronounced temporal fluctuations during its lifecycle evolution. The maximum LWC used as an indicator of mature phase is marked in Figs. 5–9.
Panel (a) is the temporal evolution of Nf, LWC, Rv, Reff, FBS and ε for fog case 1–1, the dark lines represent 5 min averaged values while the light lines are 1 min averaged values. Panel (b) is the 5 min average DSD. As the fog develops, the colors vary from blue to red, with the DSD at the time of maximum LWC marked in purple. From the time of maximum LWC to fog dissipation, the colors gradually shift from purple back to blue. Each DSD in (b) is marked in (a) using colors and symbols (circles, lines, and plus signs), where the colors denote time and the symbols indicate the modal type of the DSD. Panel (c) is the temporal evolution of visibility, temperature and wind speed, with wind direction represented by wind barbs. Panel (d) is the 1 min temporal evolution of DSDs.
Same as Fig. 5, but for fog case 4.
Same as Fig. 5, but for fog case 10.
Same as Fig. 5, but for fog case 18.
Same as Fig. 5, but for fog case 20.
Fog Case 1 (F1) was a radiation-advection fog event lasting 39 h. It formed under radiative cooling conditions, with sustained southwest warm and moist airflow supporting its long duration (Fig. A13). To enable a clearer and more detailed analysis of its microphysical characteristics, it was divided into two events at approximately the 14th hour after fog formation, based on the temporal evolution of Nf and LWC. The exact initial and end times of these two events, as well as their positions within the full fog lifecycle (F1), are provided in Table A1 and Fig. A12.
F1–1 experienced a rapid development after formation. During the first hour after fog formation, visibility decreased rapidly, and Nf increasing to approximately 600 cm−3. Following the drastic intensification there is a relatively stable phase lasting about 11 h, during which slowly decreasing temperature and steady wind direction created favorable conditions that maintain relatively high Nf and LWC. During fog formation and development, the DSD transitions among unimodal, bimodal, and trimodal. As the LWC approaches its maximum, Nf in the 5–20 µm range decreased, while those above 20 µm increased rapidly, suggesting a mass transfer from smaller to larger droplets as the DSD broadens toward larger droplets. During this stage, the DSD is predominantly bimodal. After 25 December 2006 10:00 LT, as Nf and LWC gradually decrease, the DSD again alternates between bimodal and trimodal (Fig. 5).
F4 is an advection fog event. Under the influence of a high-pressure ridge, the low-level atmosphere was stable with weak winds. Warm, moist air at 850 hPa was advected over a colder surface, cooled and condensed to form fog (Fig. A14). It shows a clear formation-dissipation evolution. About 40 min after fog formation, visibility gradually decreased as Nf and LWC increased. Temperature was relatively low during the formation and development stages, favoring fog persistence. Before LWC reaches its maximum, Nf and LWC continue to increase, followed by a gradual decrease after maximum LWC. Rv, Reff and ε are positively correlated with LWC, while FBS shows a negative correlation. At the beginning of fog formation, transitions between bimodal and trimodal are more frequent. As the fog develops, the Nf of the DSD gradually increases, with a more pronounced enhancement in droplets within the 10–20 µm size range. During fog dissipation, transitions also occur with Nf gradually decreasing across all size ranges (Fig. 6).
F10 formed as a radiation fog under high-pressure control and weak surface winds (Fig. A15) and remained stable with relatively low Nf and LWC for the first 6 h. Under these stable conditions, the DSDs persisted in trimodal with little variation. After 20 December 2007 06:00, rising temperature and variable wind direction enhanced turbulent mixing and promoted fog development. Nf and LWC increased with pronounced fluctuations, corresponding changes were observed in the DSDs, as transitions between bimodal and trimodal distributions occurring frequently. Both Rv and Reff increased with fluctuations, corresponding to the marked increase in Nf of droplets larger than 20 µm as shown by the DSD. During the dissipation stage, when Nf and LWC gradually decrease, trimodal DSDs are more prevalent with occasional transitions to bimodal (Fig. 7).
F18 is a radiation fog event formed under high-pressure conditions and near surface weak winds, driven by radiative cooling (Fig. A16). Blank areas in Fig. 8 indicate missing data. After fog formation, visibility decreases substantially, while Nf and LWC increase rapidly. Correspondingly, Nf across all size ranges rise quickly, and the DSD gradually transitions from a trimodal to a bimodal. After the LWC reaches its maximum, temperature decreases slowly and near-surface wind speeds remain weak, favouring fog maintenance. During this period and until gradual fog dissipation, the DSD alternates between bimodal and trimodal distributions.
F20 is an advection fog event. Under high-pressure control, easterly and southeasterly winds transported warm, moist marine air over a cold surface cooled by nocturnal radiative loss, providing favorable conditions for fog formation (Fig. A17). The stable temperature and wind direction as well as wind speeds below 3 m s−1 favor fog maintenance. And the FBS remained consistently high at above 60 %. For F20, the DSDs were predominantly trimodal, without the pronounced increase in the Nf of droplets larger than 10 µm as observed in F18 and F10. Although this fog case exhibits two pronounced formation and dissipation cycles, transitions between bimodal and trimodal DSDs occur infrequently, likely due to its relatively low Nf and LWC.
Across the five fog cases analyzed above, transitions of the DSD among different modal types occur in all cases and are closely linked to the characteristics of the fog life cycle. In cases with relatively low Nf and LWC or weak temporal variability, such as F4 and F20, these transitions occur less frequently. In F10, frequent changes in DSD modal types coincide with pronounced oscillations in Nf and LWC. Both F1–1 and F18 exhibit explosive increases in Nf and LWC during fog formation, and as LWC approaches its maximum, bimodal DSDs occur more frequently. However, LWC in F1–1 is slightly higher than in F18, while the mean Nf in F18 is approximately 1.5 times that in F1–1. The resulting stronger competition for water vapor may suppress the formation of droplets larger than 20 µm, preventing the sustained occurrence of a bimodal DSD.
PDF distributions of Nf(a), LWC (b), Rv(c) and Reff(d) for DSDs with different modes.
Average contributions of each DSD bin to Nf(a) and LWC (b) with different modes.
Correlation of modes and microphysical properties
Figure 10 presents the probability distributions of Nf, LWC, Rv and Reff across unimodal and multimodal DSDs. The Nf distribution of unimodal DSDs differs markedly from those of bimodal and trimodal DSDs, with its probability density function (PDF) distributions mainly concentrated below 200 and above 400 cm−3. In particular, the PDF distributions above 400 cm−3 is higher than that of bimodal and trimodal DSDs. For bimodal and trimodal DSDs, the PDF of Nf decreases with increasing Nf, although bimodal DSDs show slightly higher PDF values than trimodal DSDs in the range of 250–500 cm−3. The LWC distribution of unimodal DSDs shows characteristics similar to its Nf distribution, with the highest PDF occurring at values exceeding 0.4 g m−3. Bimodal and trimodal DSDs exhibit similar distribution when LWC is below 0.2 g m−3, whereas bimodal DSDs have higher PDF values and a broader distribution at LWC higher than 0.2 g m−3. For the Rv distribution, unimodal DSDs exhibit the narrowest range, mainly concentrated between 2 and 6 µm. Bimodal and trimodal DSDs have comparable ranges. However, trimodal DSDs are primarily concentrated between 2 and 4 µm, with the PDF decreasing rapidly beyond 4 µm, whereas bimodal DSDs show a more uniform distribution over the range of 2–7 µm. In terms of Reff, unimodal DSDs exhibit the highest PDF values in the 8–10 µm range. Bimodal DSDs show the widest Reff distribution, with PDF values above 8 µm exceeding those of trimodal DSDs. Trimodal DSDs are mainly distributed between 2 and 10 µm, with slightly higher PDF values than bimodal DSDs in the 2–5 µm range.
Figure 11 shows the average contributions of each droplet size bin to total Nf and LWC with different modes. The mean contribution of each size bin to the total Nf shows good consistency with the DSD modal distribution. For unimodal DSDs, the contribution of individual bins to Nf decreases with increasing droplet size. In contrast, bimodal DSDs exhibit an enhanced contribution to Nf in the 18–29 µm size range. Compared with trimodal DSDs, droplets larger than 10 µm contribute more to the Nf of bimodal DSDs. For LWC, in unimodal DSDs droplets in the 2–10, 16–18, and 20–29 µm diameter ranges all make noticeable contributions. In bimodal and trimodal DSDs, droplets at other sizes contribute relatively less to LWC compared with those in the 2–6 and 20–29 µm diameter ranges. However, compare to trimodal DSDs, droplets in the 12–20 µm range contribute slightly more to LWC in bimodal DSDs.
Performances and improvement of gamma and lognormal fitting
Bulk microphysical schemes commonly represent DSDs with gamma or lognormal distributions, making the accuracy of these representations critical to performance of numerical simulation. To evaluate the validity of the gamma and lognormal distribution for winter fog in Nanjing, mean DSDs of 27 observed fog events are fitted using the form of Eqs. (12) and (13) with i=1. For the fog events examined in this study, both gamma and lognormal distribution provide a good fit to the average DSD in the small-droplet range (2–10 µm), but significantly underestimates number concentrations as droplet size increases, especially the gamma distribution. In multimodal DSDs of the examined events, additional modes appear besides the mode near 2 µm, which likely contribute to the poor fit. The probability distribution of Dturn obtained from lognormal fitting across all DSDs (Fig. 2c) indicates that the other two modes occur primarily around 10 and 20 µm. Therefore, segmented gamma and lognormal fitting was conducted using 10 and 20 µm as partition points (Fig. 12b, c, d, f, g, h). When DSDs are segmented at 10 µm, the fit increasingly underestimates number concentrations for diameters above 20 µm. Segmentation at 20 µm produces good agreement in the 2–10 and 20–50 µm ranges, although substantial deviations remain in the intermediate 10–20 µm range. Based on these results, a three-segment gamma fitting approach was applied using 10 and 20 µm as partition points (Fig. 12d, h). This approach significantly improves the overall fit ranging from 2 to 50 µm, providing a more accurate representation of the whole DSDs.
Gamma and lognormal fitting of the mean spectrum: original fit (a, e), two-segment fitting with a breakpoint at 10 µm (b, f), two-segment fitting with a breakpoint at 20 µm (c, g), and three-segment fitting with breakpoints at 10 and 20 µm (d, h).
To demonstrate that the superior performance of the three-segment gamma and lognormal fitting is due to the physically meaningful segmentation based on the characteristics of DSDs, rather than merely the increased number of segments, we evaluated the performance of alternative segmented fitting. Since the gamma and lognormal distributions are nonlinear, two fitting points would fall on a straight line and cannot constrain the distribution, potentially leading to non-identifiable or ill-posed parameter estimates. Therefore, the segmentation points must satisfy two conditions: (a) the full spectrum must be divided into three segments, and (b) each segment must contain at least three bins. Under these constraints, 66 feasible segmentation combinations exist. Using each set of segmentation points, we performed gamma and lognormal fits for all DSDs and retrieved the corresponding Nf, LWC, Rv and Reff. The absolute deviations between the retrieved values and the observed ones were compared for both the alternative fittings and the fixed 10 and 20 µm segmentation (Fig. 13). For all four microphysical parameters, the deviations from the fixed-segmentation fitting are significantly smaller than those from alternative segmentation. A two-sided binomial test was conducted to evaluate the probability that the fixed segmentation outperforms the alternative segmentation. For both distributions, the 95 % confidence interval is [0.986, 1.000] with a p-value of 2.23×10-308. These p-values are far below the 0.05 significance threshold. These results confirm the effectiveness of the 10 and 20 µm segmentation for both gamma and lognormal distribution.
Boxplots of the mean absolute deviations from 66 alternative segmented fittings, with red and blue dots indicating the mean deviations from the fixed 10 and 20 µm segmentation.
To further evaluate the performance of the gamma and lognormal fitting with different breaking points, Nf, LWC, Rv and Reff were calculated based on both the original and segmented fit. These results were then compared with those derived from observations. The deviation distribution between fitted and observed results are analysed in Fig. 14, and the correlation between them are showed in Fig. A19. The non-segmented gamma fit significantly underestimates the Nf of droplets larger than 10 µm, leading to underestimation of all derived microphysical quantities. The Nf calculated from the two-segment fit with a breakpoint at 10 µm is close to observation, but the derived LWC are still underestimated. For the fit segmented at 20 µm all microphysical quantities are underestimated, likely due to underrepresentation of droplet concentrations in the 10–20 µm range. Compared to the non-segmented gamma fit approach, the three-segment fitting shows substantial improvement in the high Nf regime (Nf> 700 cm−3) and in LWC estimation as the results are tightly clustered around the zero-deviation line and the 1 : 1 line (Fig. A19). Also, the fitting accuracy for Reff and Rv is also improved for both gamma and lognormal fitting.
Deviation distribution of Nf, LWC, Rv and Reff between observation spectrum and the gamma and lognormal fitted spectrum. The number above each boxplot indicates the segmentation position (µm); “none” denotes no segmentation point applied.
Except for the non-segmented gamma fitting, the other three segmented approaches exhibit a clear pattern in estimating Reff: underestimation primarily occurs when the observed Reff< 6 µm, while overestimation tends to occur when the observed Reff> 6 µm (Fig. A19). This pattern is particularly pronounced in the fitting segmented at 20 µm, whereas the three-segment fitting shows notable improvement in reducing underestimation at lower Reff values.
Cloud optical thickness (τ) and single-scattering albedo (ω0) are key parameters for evaluating the Twomey effect (Stephens, 1984; Twomey and Bohren, 1980). τ can be calculated with
τ=∫3LWCDeffdz
where z is the thickness of the cloud or fog layer (Stephens, 1978). When assuming cloud or fog is vertically homogeneous, Eq. (16) can be simplified as
τdz=3LWCDeff
where τdz represents the average optical thickness per layer. The single-scattering albedo (ω0) can be expressed by
1-ω0=1.7kwReff
where kw is the complex part of the refractive index of water. Equation (18) indicates the critical role of Reff in the Twomey effect (Wang et al., 2019).
To more precisely assess the potential climate impact of inaccuracy in Reff estimates from gamma and lognormal fitting, Eqs. (17) and (18) were used to calculate absorption coefficient and optical thickness based on both observed and fitted Reff. This allowed evaluation of the extent to which gamma and lognormal fitting overestimates or underestimates these parameters. Results are showed in Fig. 15.
Correlation between absorption coefficient (1-ω0) and optical thickness (τ) derived from observed spectrum and those computed from the gamma and lognormal fitted spectrum with no breakpoint (a, b, i, j), breakpoint at 10 µm (c, d, k, l), breakpoint at 20 µm (e, f, m, n) and breakpoint at 10 and 20 µm (g, h, o, p).
Due to its significant underestimation of the Nf of droplets larger than 10 µm, the non-segmented gamma and lognormal fitting notably underestimate absorption coefficient (1-ω0) and optical thickness (τ) by up to 90 %. Compared to the fitting segmented at 20 µm, the fitting segmented at 10 µm more accurately captures absorption coefficient and optical thickness, yet still exhibits up to 50 % overestimation or underestimation of absorption coefficient. It is noteworthy that the 20 µm-segmented fitting generally underestimates optical thickness, likely due to the underestimation of the Nf for droplets in the 10–20 µm range.
In contrast, the three-segment fitting significantly improves the estimation of both absorption coefficient and optical thickness, with most deviations confined within ±20 %. The most notable improvements lie in reducing the underestimation of absorption coefficient and the overestimation of optical thickness, as the percentages of both underestimation and overestimation become very small.
Correlation between the N0, μ and λ derived from the three-segment gamma fitting.
The interrelationships among the three fitting parameters obtained from the three-segment gamma fitting are shown in Fig. 16. For each segment, N0 exhibits a negative correlation with both μ and λ, while μ and λ are positively correlated. In the D<10µm segment, larger N0 and λ are observed, indicating a narrower and more concentrated spectrum with a steep decline in number concentration as droplet size increases. For D>10µm, N0 shows a wider distribution with smaller λ, reflecting broader spectrum and higher number concentrations of larger droplets. Compared to the D>10µm segments, N0 in the D<10µm range is more tightly gathered in a range of 100–103 cm−3µm−(1+μ), suggesting that small droplet concentrations are higher and more consistent across different DSDs. Because the fitting parameters of lognormal distribution do not exhibit clear correlations, they are not shown here.
Conclusions
As a key parameter of fog microphysical processes, the droplet size distribution (DSD) is influenced by multiple macro- and micro-scale factors, exhibits significant temporal and spatial variability, and evolves throughout the fog lifecycle, thereby posing challenges for accurate fog prediction (Niu et al., 2012; Nelli et al., 2024). To gain a better understanding of how the DSD evolves over the fog lifecycle, this study investigates the microphysical characteristics of 27 winter fog events in Nanjing under polluted conditions, with a focus on the evolution of droplet size distributions (DSDs) throughout the fog lifecycle. Among the 27 fog cases, DSDs with single mode (2 µm), double mode (2, 6–18 µm) and triple mode (2, 6–12, 18–26 µm) were observed. The average Nf, LWC, Rv and Reff vary over the ranges of 25–586 cm−3, 0–0.27 g m−3, 1.6–6 µm, 1.9–8.2 µm, which shows greater Nf, lower LWC and smaller droplets comparing to other clean regions such as the tropical rainforests of southwestern China (Wang et al., 2021). The main findings are as follows:
Among all fog cases, radiation fog accounts for the largest proportion and radiation-advection fog tends to persist longer. Variations in the number of DSD modes are closely linked to fog lifecycle characteristics and to changes in physical variables such as Nf and LWC. For fog events with relatively low Nf and LWC or relatively less intense formation and dissipation processes, transitions in DSD modal type are infrequent and occur mainly during the formation and dissipation stages. When Nf and LWC increase sharply, the DSD shows frequent modal transitions with an increasing prevalence of bimodal DSD. Moreover, when Nf and LWC exhibit strong oscillations, transitions among different modal types occur simultaneously and at high frequency.
Comparison of the retrieved physical parameters from segmented gamma and lognormal fitting with observations indicates that the three-segment fitting yields the best performance, especially in improving Nf and LWC estimation. Meanwhile, the three-segment fitting reduces the estimation deviations in Reff, absorption coefficient and optical thickness from up to 90 % in the non-segmented fitting to below 20 %, demonstrating its effectiveness in improving fog DSD representation and microphysical characteristic retrieval.
These findings advance our understanding of fog droplet size distribution (DSD) evolution during fog lifecycles and the correlations between DSD modes and microphysical properties in polluted urban regions. The improved segmented gamma and lognormal fitting offer a new perspective for DSD parameterization and demonstrates strong potential for improving the representation of cloud/fog microphysical processes in weather prediction and climate models. It should be noted that this study is primarily based on observational data and focuses on analyzing the microphysical characteristics and their evolutions. The underlying physical processes and mechanisms governing fog evolution as well as the interactions and relative importance of different controlling factors, remain to be clarified through enhanced analyses of aerosol background conditions and sensitivity experiments on relevant physical processes using numerical simulations. In addition, only a three-parameter gamma and lognormal distribution was used to fit and refine the mean DSD. The comparative performance of alternative distribution and evaluate the influence of different parameterizations on fitting accuracy could be explored in future studies.
Details on the upper and lower bounds of parameters used for gamma and lognormal fitting
In this study, DSDs were fitted using gamma and lognormal distributions to determine the number of modes and their corresponding Dturn, which are also shown below.
For gamma distribution:
nD=∑i=13niD=∑i=13N0,iDμie-λiD
For the lognormal distribution:
nD=∑i=13niD=∑i=13N0,i2π1/2Dilnσg,iexp-(lnDi-lnDg,i)22(lnσg,i)2
Fits with different numbers of modes were performed, and the optimal solution was obtained for each distribution.
For the gamma distribution, no predefined size ranges were prescribed. Unimodal, bimodal, and trimodal gamma distributions were fitted separately, allowing the fitting procedure to automatically identify the optimal solution.
In the lognormal distribution, three possible modes (approximately 2, 10–18, and 18–23 µm) were identified based on the mean DSD of all cases (Fig. A1). To minimize subjective influence, relatively broad bounds of Dg were assigned to the parameters under different values of i. For unimodal distributions (i=1), Dg were constrained to the range 2–50 µm. In bimodal distributions (i=2), two candidate mode diameter combinations were considered: (2, 11 µm) and (2, 21 µm). Accordingly, different lower and upper bounds were specified for Dg. For the 2 and 11 µm combination, the bounds were set to 2–10 and 5–20 µm; for the 2 and 21 µm combination, the bounds were 2–10 and 15–50 µm. In trimodal distributions (i=3), the bounds for Dg were set to 2–50, 5–20, 15–50 µm. The discrete representation of the lognormal distribution parameter σg is given by
lnσg=∑ni(lnDi-lnDg)2∑ni-11/2
in which lnDg=∑nilnDi∑ni. In this study, the maximum calculated value of σg is approximately 2.34. Therefore, the range of σg was set to 0–2.5. In the meantime, N0 was set to 0–2000 cm−3.
Average spectrums of all fog DSDs.
Model selection criteria for unimodal, bimodal, and trimodal fits and representative examples
In this study, based on comparisons with the turning-point method and the gamma fitting results, the lognormal distribution was selected to determine the number of DSD modes and the corresponding Dturn. Because two sets of fitting parameters were specified for the bimodal distribution, two fitting results were obtained, resulting in four candidate fits for each DSD in total. Intersections between adjacent fitted components were required, and the Dturn derived from these fits were not allowed to occur in the same or adjacent bins. Among the fits that satisfy this criterion, the one with the lowest AIC and BIC was selected to determine the number of modes and the mode diameters of the DSD.
The following six representative examples are provided to illustrate the four candidate fits obtained for each individual DSD and the procedure by which the optimal fit is selected under the proposed criteria.
Single mode examples
Figure A2 shows a unimodal DSD during the fog development stage prior to the maximum LWC in F1–1. The Dturn of mode 2, defined by the intersection between mode 1 and mode 2, falls within the same bin as the Dturn of mode 1 and is therefore considered invalid and shown as NaN in the figure. Despite the trimodal fit yielding the lowest AIC and BIC values, this fit is discarded. Consequently, the unimodal fit (panel a) with the second-lowest AIC and BIC values was selected to determine the number of modes in this DSD.
Example of single mode DSD.
Double mode examples
Figure A3 shows the DSD during the fog development stage prior to the maximum LWC in F1–1, corresponding to the light-blue curve in Fig. 7b. For this DSD, the trimodal fit in panel (d) and the bimodal fit in panel (c) with lower AIC and BIC values were discarded because mode 1 and mode 2 do not intersect (appearing as NaN in the figure), making it impossible to determine the second mode. In the bimodal fit shown in panel (b), both modes 1 and 2 have corresponding Dturn that do not fall within the same or adjacent bins, and the sum of the AIC and BIC values is smaller than that of the unimodal fit. Therefore, this fit was used to determine the number of modes and the corresponding Dturn for this DSD.
Example 1 of double mode DSD.
Example 2 of double mode DSD.
Example 3 of double mode DSD.
Example 4 of double mode DSD.
Example of triple mode DSD.
Figure A4 shows the DSD at the time when LWC in F1–1 reaches its maximum, corresponding to the purple curve in Fig. 7b. For this DSD, the trimodal fit (panel d) is discarded because the Dturn of the first and second mode fall within the same bin. Accordingly, the Dturn of mode 2 is shown as NaN in the figure. Among the remaining fits, the fit in panel (c) yields the lowest AIC and BIC values. The intersection between mode 1 and mode 2 defines a corresponding Dturn that does not fall within the same or adjacent bins. Therefore, the bimodal fit in panel (c), which has lower AIC and BIC values, is selected to represent the number of modes and the mode diameters of this DSD.
Figure A5 corresponds to the orange DSD in F10 prior to the maximum LWC (Fig. 8b). Compared with Fig. A7, this DSD exhibits a broader mode in the 10–25 µm size range. In the trimodal fit (panel d), mode 2 and mode 3 do not intersect, so the modes cannot be determined. This fit is therefore discarded, with the Dturn of mode 3 shown as NaN in the figure. The other three fits are all valid. Therefore, the bimodal fit in panel (c) with the smaller combined AIC and BIC is selected to determine the number of modes and the Dturn of this DSD.
Figure A6 corresponds to the green DSD in F1–1 prior to the maximum LWC (Fig. 7b). In the fitting procedure, Dg was constrained to 5–20 µm for panel (b) and to 15–20 µm for panel (c). For this DSD, the retrieved Dg falls within the overlapping range of the two parameter bounds (15–20 µm), resulting in identical fitting results in panels (b) and (c). Because they yield the lowest AIC and BIC values, and the corresponding Dturns do not fall within the same or adjacent bins, this DSD is classified as bimodal.
Triple mode examples
Figure A7 shows the mode identification result for the yellow curve prior to the maximum LWC in F10 (Fig. 8b). In the trimodal fit (panel d), intersections exist between mode 1 and mode 2 as well as between mode 2 and mode 3, and the resulting Dturn satisfy the physical requirement that they do not fall within the same or adjacent bins. Moreover, this fit yields the minimum combined AIC and BIC among all fits. Therefore, it is used to define the number of modes and the Dturn of this DSD.
MODIS 3.9 <inline-formula><mml:math id="M399" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi></mml:mrow></mml:math></inline-formula>m shortwave infrared and visible channel imagery for F1–1, F10, and F20
For F4 and F18, no satellite imagery is available within the fog period because the overpass times of the polar-orbiting satellite did not coincide with the observations.
Aqua MODIS 3.9 µm shortwave infrared and visible channel imagery for Fog 1–1, with the observation site marked by a red cross. Satellite imagery from MODIS (Aqua), NASA LAADS DAAC. Administrative boundaries derived from the CN-border dataset, based on data from the National Catalogue Service For Geographic Information (https://www.webmap.cn, last access: 25 February 2026).
Terra MODIS 3.9 µm shortwave infrared and visible channel imagery for Fog 10, with the observation site marked by a red cross. Satellite imagery from MODIS (Terra), NASA LAADS DAAC; administrative boundaries from the CN-border dataset (National Catalogue Service For Geographic Information, https://www.webmap.cn).
Terra MODIS 3.9 µm shortwave infrared and visible channel imagery for Fog 20, with the observation site marked by a red cross. Satellite imagery from MODIS (Terra), NASA LAADS DAAC; administrative boundaries from the CN-border dataset (National Catalogue Service For Geographic Information, https://www.webmap.cn).
Summary table of the 27 fog events and details of fog type identification
In this study, some cases classified as radiation fog have initial times in the early morning rather than at night. For these less typical cases, additional details on the fog types identify are provided here. Fog case 6 formed during the night and dissipated around 06:00. According to the observation log, light mist persisted at the site until the formation of Fog case 7, but it did not meet the criteria used in this study to define fog (Nf>10 cm−3 and LWC > 10−3 g m−3). Therefore, Fog case 6 and Fog case 7 were treated separately, although their formation mechanisms were the same. For Fog case 8, the night preceding fog formation was clear, with strong radiative cooling and very weak winds, conditions favorable for fog development. Fog case 8 was therefore classified as a radiation fog. Fog case 26 developed in the morning. Light drizzle occurred the previous afternoon, followed by light mist during the night, indicating that precipitation played a role in its formation.
Initial and end times, classification, visibility, meteorological variables (temperature, wind speed) and microphysical properties of the 27 fog events, the first row shows the mean values, in the parentheses are the 25th and 75th percentiles of each characteristic. In the first row of the column titled Modes, the DSD types observed in each case are listed, where U, D, and T denote unimodal, bimodal, and trimodal DSD, respectively.
The Dturn distributions of all DSDs obtained using gamma distribution at 1 min resolution (a), 5 min resolution (b), and the Dturn distributions by DSD type at both resolutions (c).
The temporal evolution of Nf, LWC, Rv, Reff, FBS and ε for fog case 1, the dark lines represent 5 min averaged values while the light lines are 1 min averaged values, the black dashed line is the boundary between F1–1 and F1–2.
Synoptic conditions of Fog 1–1. Red dots indicate observation stations; blue lines represent isobars or geopotential height contours; red lines are isotherms, dashed when temperatures are below 0 °C; green lines denote constant specific humidity.
Synoptic conditions of Fog 4.
Synoptic conditions of Fog 10.
Synoptic conditions of Fog 18.
Synoptic conditions of Fog 20.
The time series of DSDs at 1s resolution of Fog Case 1–1, Fog Case 4, Fog Case 10, Fog Case 18 and Fog Case 20.
Correlation between Nf, LWC, Rv and Reff derived from observed spectrum and those computed from the gamma and lognormal fitted spectrum with no breakpoint (first column), breakpoint at 10 µm (second column), breakpoint at 20 µm (third column) and breakpoint at 10 and 20 µm (fourth column). Red dots indicate overestimation of the microphysical properties by the gamma fit, while blue dots indicate underestimation.
Discussion of the sensitivity of segmented gamma and lognormal fits to the exclusion of potential modes in 30–50 <inline-formula><mml:math id="M490" display="inline"><mml:mrow class="unit"><mml:mi mathvariant="normal">µ</mml:mi></mml:mrow></mml:math></inline-formula>m
To examine the impact of excluding potential modes in the 30–50 µm range on segmented gamma and lognormal fitting, an additional segmentation points at 40 µm was considered. Using the same approach as in Sect. 3.4, segmented fits were applied to all DSDs, and the fitting results were used to retrieve Nf, LWC, Rv, Reff, absorption coefficient and optical thickness. The deviations of these retrieved quantities from the observations were then calculated. Figure A19 compares the deviations obtained from fits with segmentation points at 10 and 20 µm and those including the additional 40 µm segmentation point. The results show that introducing the 40 µm segmentation point does not lead to a significant improvement in fitting performance, indicating that the omission of potential modes larger than 30 µm does not substantially affect the results of this study.
Comparison between segmented fitting results with breakpoints at 10 and 20 µm (a) and those with additional breakpoints at 10, 20, and 40 µm (b).
Data availability
Fog DSD data used in this study is available at Zenodo (Zhang, 2026, 10.5281/zenodo.18552099). Except for fog case 1–3 and 18–20, where visibility was calculated using Eqs. (9) and (10), visibility data for other cases were obtained from PWD measurements. Meteorological variables including temperature, wind speed, and wind direction were sourced from the Nanjing University of Information Science and Technology (NUIST) automatic weather station, except for wind speed in fog case 1–3 and temperature, wind speed, and wind direction in fog case 15–17, which were taken from 10.24381/cds.adbb2d47 (Copernicus Climate Change Service, 2023a). Synoptic conditions for all fog cases were derived from 10.24381/cds.bd0915c6 (Copernicus Climate Change Service, 2023b) and 10.24381/cds.adbb2d47 (Copernicus Climate Change Service, 2023a). The satellite imagery in the Appendix (Figs. A8–A10) is obtained from https://ladsweb.modaps.eosdis.nasa.gov/ (Aqua MODIS and Terra MODIS, last access: 25 February 2026).
Author contributions
JZ and XL shaped the concept of this study and performed data analysis. JZ prepared the figures and wrote the initial draft. XL supervised this work and revised the article. ZA gave suggestions on data processing and visualization. JL provided guidance for the synoptic analysis. DX provided guidance on locating and downloading the MODIS data.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Acknowledgements
We gratefully acknowledge Professor Shengjie Niu and his research team for their efforts in acquiring and sharing the fog microphysical observation data that supported this study. We would like to acknowledge the High Performance Computing Center of Nanjing University of Information Science and Technology and the National Key Scientific and Technological Infrastructure project “Earth System Numerical Simulation Facility” (EarthLab) for their support of this work.
Financial support
This study was supported by the National Natural Science Foundation of China (grant nos. 42575085, 42061134009 and 41975176) and the National Key Scientific and Technological Infrastructure project “Earth System Numerical Simulation Facility” (EarthLab).
Review statement
This paper was edited by Greg McFarquhar and reviewed by two anonymous referees.
ReferencesAkaike, H.: A new look at the statistical model identification, IEEE Trans. Automat. Contr., 19, 716–723, 10.1109/TAC.1974.1100705, 1974.Boudala, F. S., Wu, D., Isaac, G. A., and Gultepe, I.: Seasonal and Microphysical Characteristics of Fog at a Northern Airport in Alberta, Canada, Remote Sens., 14, 4865, 10.3390/rs14194865, 2022.Boutle, I., Price, J., Kudzotsa, I., Kokkola, H., and Romakkaniemi, S.: Aerosol–fog interaction and the transition to well-mixed radiation fog, Atmos. Chem. Phys., 18, 7827–7840, 10.5194/acp-18-7827-2018, 2018.Boutle, I., Angevine, W., Bao, J.-W., Bergot, T., Bhattacharya, R., Bott, A., Ducongé, L., Forbes, R., Goecke, T., Grell, E., Hill, A., Igel, A. L., Kudzotsa, I., Lac, C., Maronga, B., Romakkaniemi, S., Schmidli, J., Schwenkel, J., Steeneveld, G.-J., and Vié, B.: Demistify: a large-eddy simulation (LES) and single-column model (SCM) intercomparison of radiation fog, Atmos. Chem. Phys., 22, 319–333, 10.5194/acp-22-319-2022, 2022.Brenguier, J.-L., Pawlowska, H., Schüller, L., Preusker, R., Fischer, J., and Fouquart, Y.: Radiative Properties of Boundary Layer Clouds: Droplet Effective Radius versus Number Concentration, J. Atmos. Sci., 57, 803–821, 10.1175/1520-0469(2000)057<0803:RPOBLC>2.0.CO;2, 2000.Brown, R.: A numerical study of radiation fog with an explicit formulation of the microphysics, Q. J. Roy. Meteor. Soc., 106, 781–802, 10.1002/qj.49710645010, 1980.Chakrabarti, A. and Ghosh, J. K.: AIC, BIC and Recent Advances in Model Selection, in: Philosophy of Statistics, vol. 7, edited by: Bandyopadhyay, P. S. and Forster, M. R., North-Holland, Amsterdam, 583–605, 10.1016/B978-0-444-51862-0.50018-6, 2011.Chen, J., Liu, Y., Zhang, M., and Peng, Y.: New understanding and quantification of the regime dependence of aerosol-cloud interaction for studying aerosol indirect effects, Geophys. Res. Lett., 43, 1780–1787, 10.1002/2016GL067683, 2016.Copernicus Climate Change Service: ERA5 hourly data on single levels from 1940 to present, last access: 25 February 2026, Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [data set], 10.24381/cds.adbb2d47, 2023a.Copernicus Climate Change Service: ERA5 hourly data on pressure levels from 1940 to present, last access: 25 February 2026, Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [data set], 10.24381/cds.bd0915c6, 2023b.Cui, C., Bao, Y., Yuan, C., Li, Z., and Zong, C.: Comparison of the performances between the WRF and WRF-LES models in radiation fog – A case study, Atmos. Res., 226, 76–86, 10.1016/j.atmosres.2019.04.003, 2019.Dorman, C. E., Hoch, S. W., Gultepe, I., Wang, Q., Yamaguchi, R. T., Fernando, H. J. S., and Krishnamurthy, R.: Large-Scale Synoptic Systems and Fog During the C-FOG Field Experiment, Bound.-Lay. Meteorol., 181, 171–202, 10.1007/s10546-021-00641-1, 2021.Elias, T., Haeffelin, M., Drobinski, P., Gomes, L., Rangognio, J., Bergot, T., Chazette, P., Raut, J.-C., and Colomb, M.: Particulate contribution to extinction of visible radiation: Pollution, haze, and fog, Atmos. Res., 92, 443–454, 10.1016/j.atmosres.2009.01.006, 2009.Elias, T., Dupont, J.-C., Hammer, E., Hoyle, C. R., Haeffelin, M., Burnet, F., and Jolivet, D.: Enhanced extinction of visible radiation due to hydrated aerosols in mist and fog, Atmos. Chem. Phys., 15, 6605–6623, 10.5194/acp-15-6605-2015, 2015. Friedlein, M. T.: DENSE FOG CLIMATOLOGY, B. Am. Meteorol. Soc., 85, 515–517, 2004.Ge, P., Zhang, Y., Fan, S., Wang, Y., Wu, H., Wang, X., and Zhang, S.: Observational study of microphysical and chemical characteristics of size-resolved fog in different regional backgrounds in China, Sci. Total Environ., 950, 175329, 10.1016/j.scitotenv.2024.175329, 2024.Gultepe, I. and Milbrandt, J. A.: Microphysical Observations and Mesoscale Model Simulation of a Warm Fog Case during FRAM Project, in: Fog and Boundary Layer Clouds: Fog Visibility and Forecasting, Birkhäuser Basel, 1161–1178, 10.1007/978-3-7643-8419-7_4, 2007.Gultepe, I. and Milbrandt, J. A.: Probabilistic Parameterizations of Visibility Using Observations of Rain Precipitation Rate, Relative Humidity, and Visibility, J. Appl. Meteorol. Clim., 49, 36–46, 10.1175/2009JAMC1927.1, 2010.Gultepe, I., Müller, M. D., and Boybeyi, Z.: A New Visibility Parameterization for Warm-Fog Applications in Numerical Weather Prediction Models, J. Appl. Meteorol. Clim., 45, 1469–1480, 10.1175/JAM2423.1, 2006.Gultepe, I., Tardif, R., Michaelides, S. C., Cermak, J., Bott, A., Bendix, J., Müller, M. D., Pagowski, M., Hansen, B., Ellrod, G., Jacobs, W., Toth, G., and Cober, S. G.: Fog Research: A Review of Past Achievements and Future Perspectives, Pure Appl. Geophys., 164, 1121–1159, 10.1007/s00024-007-0211-x, 2007.Gultepe, I., Pearson, G., Milbrandt, J. A., Hansen, B., Platnick, S., Taylor, P., Gordon, M., Oakley, J. P., and Cober, S. G.: The Fog Remote Sensing and Modeling Field Project, B. Am. Meteorol. Soc., 90, 341–360, 10.1175/2008BAMS2354.1, 2009.Gultepe, I., Kuhn, T., Pavolonis, M., Calvert, C., Gurka, J., Heymsfield, A. J., Liu, P. S. K., Zhou, B., Ware, R., Ferrier, B., Milbrandt, J., and Bernstein, B.: Ice Fog in Arctic During FRAM–Ice Fog Project: Aviation and Nowcasting Applications, B. Am. Meteorol. Soc., 95, 211–226, 10.1175/BAMS-D-11-00071.1, 2014.Gultepe, I., Heymsfield, A. J., Fernando, H. J. S., Pardyjak, E., Dorman, C. E., Wang, Q., Creegan, E., Hoch, S. W., Flagg, D. D., Yamaguchi, R., Krishnamurthy, R., Gaberšek, S., Perrie, W., Perelet, A., Singh, D. K., Chang, R., Nagare, B., Wagh, S., and Wang, S.: A Review of Coastal Fog Microphysics During C-FOG, Bound.-Lay. Meteorol., 181, 227–265, 10.1007/s10546-021-00659-5, 2021.Guo, L., Guo, X., Fang, C., and Zhu, S.: Observation analysis on characteristics of formation, evolution and transition of a long-lasting severe fog and haze episode in North China, Sci. China Earth Sci., 58, 329–344, 10.1007/s11430-014-4924-2, 2015.Haeffelin, M., Bergot, T., Elias, T., Tardif, R., Carrer, D., Chazette, P., Colomb, M., Drobinski, P., Dupont, E., Dupont, J.-C., Gomes, L., Musson-Genon, L., Pietras, C., Plana-Fattori, A., Protat, A., Rangognio, J., Raut, J.-C., Rémy, S., Richard, D., Sciare, J., and Zhang, X.: Parisfog: Shedding new Light on Fog Physical Processes, B. Am. Meteorol. Soc., 91, 767–783, 10.1175/2009BAMS2671.1, 2010.Hammer, E., Gysel, M., Roberts, G. C., Elias, T., Hofer, J., Hoyle, C. R., Bukowiecki, N., Dupont, J.-C., Burnet, F., Baltensperger, U., and Weingartner, E.: Size-dependent particle activation properties in fog during the ParisFog 2012/13 field campaign, Atmos. Chem. Phys., 14, 10517–10533, 10.5194/acp-14-10517-2014, 2014.Jia, X., Quan, J., Zheng, Z., Liu, X., Liu, Q., He, H., and Liu, Y.: Impacts of Anthropogenic Aerosols on Fog in North China Plain, J. Geophys. Res.-Atmos., 124, 252–265, 10.1029/2018JD029437, 2019.Kessler, E.: On the Distribution and Continuity of Water Substance in Atmospheric Circulations, in: On the Distribution and Continuity of Water Substance in Atmospheric Circulations, edited by: Kessler, E., American Meteorological Society, Boston, MA, 1–84, 10.1007/978-1-935704-36-2_1, 1969.Koenig, L. R.: Numerical Experiments Pertaining to Warm-Fog Clearing, Mon. Weather Rev., 99, 227–241, 10.1175/1520-0493(1971)099<0227:NEPTWC>2.3.CO;2, 1971.Koračin, D., Dorman, C. E., Lewis, J. M., Hudson, J. G., Wilcox, E. M., and Torregrosa, A.: Marine fog: A review, Atmos. Res., 143, 142–175, 10.1016/j.atmosres.2013.12.012, 2014.Khain, A. P., Beheng, K. D., Heymsfield, A., Korolev, A., Krichak, S. O., Levin, Z., Pinsky, M., Phillips, V., Prabhakaran, T., Teller, A., van den Heever, S. C., and Yano, J.-I.: Representation of microphysical processes in cloud-resolving models: Spectral (bin) microphysics versus bulk parameterization, Rev. Geophys., 53, 247–322, 10.1002/2014RG000468, 2015. Kunkel, B. A.: Microphysical properties of fog at Otis AFB, Air Force Geophysics Laboratories, Air Force Systems Command, United States, Report no. AFGL-TR-82-0026, 1982.Lakra, K. and Avishek, K.: A review on factors influencing fog formation, classification, forecasting, detection and impacts, Rend. Lincei-Sci. Fis., 33, 319–353, 10.1007/s12210-022-01060-1, 2022.Li, J., Wang, X., Chen, J., Zhu, C., Li, W., Li, C., Liu, L., Xu, C., Wen, L., Xue, L., Wang, W., Ding, A., and Herrmann, H.: Chemical composition and droplet size distribution of cloud at the summit of Mount Tai, China, Atmos. Chem. Phys., 17, 9885–9896, 10.5194/acp-17-9885-2017, 2017.Liu, Y. and Daum, P. H.: Indirect warming effect from dispersion forcing, Nature, 419, 580–581, 10.1038/419580a, 2002.Lu, C., Liu, Y., Niu, S., Zhao, L., Yu, H., and Cheng, M.: Examination of microphysical relationships and corresponding microphysical processes in warm fogs, Acta Meteorol. Sin., 27, 832–848, 10.1007/s13351-013-0610-0, 2013.Lu, C., Liu, Y., Yum, S. S., Chen, J., Zhu, L., Gao, S., Yin, Y., Jia, X., and Wang, Y.: Reconciling Contrasting Relationships Between Relative Dispersion and Volume-Mean Radius of Cloud Droplet Size Distributions, J. Geophys. Res.-Atmos., 125, e2019JD031868, 10.1029/2019JD031868, 2020.Martinet, P., Cimini, D., Burnet, F., Ménétrier, B., Michel, Y., and Unger, V.: Improvement of numerical weather prediction model analysis during fog conditions through the assimilation of ground-based microwave radiometer observations: a 1D-Var study, Atmos. Meas. Tech., 13, 6593–6611, 10.5194/amt-13-6593-2020, 2020.Mazoyer, M., Lac, C., Thouron, O., Bergot, T., Masson, V., and Musson-Genon, L.: Large eddy simulation of radiation fog: impact of dynamics on the fog life cycle, Atmos. Chem. Phys., 17, 13017–13035, 10.5194/acp-17-13017-2017, 2017.Mazoyer, M., Burnet, F., and Denjean, C.: Experimental study on the evolution of droplet size distribution during the fog life cycle, Atmos. Chem. Phys., 22, 11305–11321, 10.5194/acp-22-11305-2022, 2022. Meyer, M. B., Jiusto, J. E., and Lala, G. G.: Measurements of visual range and radiation-fog (haze) microphysics, J. Atmos. Sci., 37, 622–629, 1980.Nelli, N., Francis, D., Abida, R., Fonseca, R., Masson, O., and Bosc, E.: In-situ measurements of fog microphysics: Visibility parameterization and estimation of fog droplet sedimentation velocity, Atmos. Res., 309, 107570, 10.1016/j.atmosres.2024.107570, 2024.Niu, S., Lu, C., Liu, Y., Zhao, L., Lü, J., and Yang, J.: Analysis of the microphysical structure of heavy fog using a droplet spectrometer: A case study, Adv. Atmos. Sci., 27, 1259–1275, 10.1007/s00376-010-8192-6, 2010.Niu, S. J., Liu, D. Y., Zhao, L. J., Lu, C. S., Lü, J. J., and Yang, J.: Summary of a 4-Year Fog Field Study in Northern Nanjing, Part 2: Fog Microphysics, Pure Appl. Geophys., 169, 1137–1155, 10.1007/s00024-011-0344-9, 2012.Payra, S. and Mohan, M.: Multirule Based Diagnostic Approach for the Fog Predictions Using WRF Modelling Tool, Adv. Meteorol., 2014, 456065, 10.1155/2014/456065, 2014.Peterka, A., Thompson, G., and Geresdi, I.: Numerical prediction of fog: A novel parameterization for droplet formation, Q. J. Roy. Meteor. Soc., 150, 2203–2222, 10.1002/qj.4704, 2024.Pinnick, R. G., Hoihjelle, D. L., Fernandez, G., Stenmark, E. B., Lindberg, J. D., Hoidale, G. B., and Jennings, S. G.: Vertical Structure in Atmospheric Fog and Haze and Its Effects on Visible and Infrared Extinction, Journal of Atmospheric Sciences, 35, 2020–2032, 10.1175/1520-0469(1978)035<2020:VSIAFA>2.0.CO;2, 1978.Price, J. D.: On the Formation and Development of Radiation Fog: An Observational Study, Bound.-Lay. Meteorol., 172, 167–197, 10.1007/s10546-019-00444-5, 2019.Roach, W. T.: On the effect of radiative exchange on the growth by condensation of a cloud or fog droplet, Q. J. Roy. Meteor. Soc., 102, 361–372, 10.1002/qj.49710243207, 1976.Roach, W. T., Brown, R., Caughey, S. J., Garland, J. A., and Readings, C. J.: The physics of radiation fog: I – a field study, Q. J. Roy. Meteor. Soc., 102, 313–333, 10.1002/qj.49710243204, 1976.Román-Cascón, C., Steeneveld, G. J., Yagüe, C., Sastre, M., Arrillaga, J. A., and Maqueda, G.: Forecasting radiation fog at climatologically contrasting sites: evaluation of statistical methods and WRF, Q. J. Roy. Meteor. Soc., 142, 1048–1063, 10.1002/qj.2708, 2016.Ryznar, E.: Advection-radiation fog near lake Michigan, Atmos. Environ., 11, 427–430, 10.1016/0004-6981(77)90004-X, 1977.Sampurno Bruijnzeel, L., Eugster, W., and Burkard, R.: Fog as a Hydrologic Input, in: Encyclopedia of Hydrological Sciences, John Wiley & Sons, 10.1002/0470848944.hsa041, 2006. Schwarz, G.: Estimating the Dimension of a Model, Ann. Stat., 6, 461–464, 1978.Shao, N., Lu, C., Jia, X., Wang, Y., Li, Y., Yin, Y., Zhu, B., Zhao, T., Liu, D., Niu, S., Fan, S., Yan, S., and Lv, J.: Radiation fog properties in two consecutive events under polluted and clean conditions in the Yangtze River Delta, China: a simulation study, Atmos. Chem. Phys., 23, 9873–9890, 10.5194/acp-23-9873-2023, 2023.Shen, C., Zhao, C., Ma, N., Tao, J., Zhao, G., Yu, Y., and Kuang, Y.: Method to Estimate Water Vapor Supersaturation in the Ambient Activation Process Using Aerosol and Droplet Measurement Data, J. Geophys. Res.-Atmos., 123, 10606–10619, 10.1029/2018JD028315, 2018.Stephens, G. L.: Radiation Profiles in Extended Water Clouds. I: Theory, J. Atmos. Sci., 35, 2111–2122, 10.1175/1520-0469(1978)035<2111:RPIEWC>2.0.CO;2, 1978.Stephens, G. L.: The Parameterization of Radiation for Numerical Weather Prediction and Climate Models, Mon. Weather Rev., 112, 826–867, 10.1175/1520-0493(1984)112<0826:TPORFN>2.0.CO;2, 1984.Stewart, D. A. and Essenwanger, O. M.: A survey of fog and related optical propagation characteristics, Rev. Geophys., 20, 481–495, 10.1029/RG020i003p00481, 1982. Symonds, M. R. E. and Moussalli, A.: A brief guide to model selection, multimodel inference and model averaging in behavioural ecology using Akaike's information criterion, Behav. Ecol. Sociobiol., 65, 13–21, 10.1007/s00265-010-1037-6, 2011.Tampieri, F. and Tomasi, C.: Size distribution models of fog and cloud droplets in terms of the modified gamma function, Tellus, 28, 333–347, 10.3402/tellusa.v28i4.10300, 1976.Tudor, M.: Impact of horizontal diffusion, radiation and cloudiness parameterization schemes on fog forecasting in valleys, Meteorol. Atmos. Phys., 108, 57–70, 10.1007/s00703-010-0084-x, 2010.Turton, J. D. and Brown, R.: A comparison of a numerical model of radiation fog with detailed observations, Q. J. Roy. Meteor. Soc., 113, 37–54, 10.1002/qj.49711347504, 1987.Twomey, S. and Bohren, C. F.: Simple Approximations for Calculations of Absorption in Clouds, J. Atmos. Sci., 37, 2086–2095, 10.1175/1520-0469(1980)037<2086:SAFCOA>2.0.CO;2, 1980.Wang, H., Zhang, Z., Liu, D., Zhu, Y., Zhang, X., and Yuan, C.: Study on a Large-Scale Persistent Strong Dense Fog Event in Central and Eastern China, Adv. Meteorol., 2020, 8872334, 10.1155/2020/8872334, 2020.Wang, Y., Niu, S., Lv, J., Lu, C., Xu, X., Wang, Y., Ding, J., Zhang, H., Wang, T., and Kang, B.: A New Method for Distinguishing Unactivated Particles in Cloud Condensation Nuclei Measurements: Implications for Aerosol Indirect Effect Evaluation, Geophys. Res. Lett., 46, 14185–14194, 10.1029/2019GL085379, 2019.Wang, Y., Niu, S., Lu, C., Lv, J., Zhang, J., Zhang, H., Zhang, S., Shao, N., Sun, W., Jin, Y., and Song, Q.: Observational study of the physical and chemical characteristics of the winter radiation fog in the tropical rainforest in Xishuangbanna, China, Sci. China Earth Sci., 64, 1982–1995, 10.1007/s11430-020-9766-4, 2021. Wendisch, M., Mertes, S., Heintzenberg, J., Wiedensohler, A., Schell, D., Wobrock, W., Frank, G., Martinsson, B. G., Fuzzi, S., and Orsi, G.: Drop size distribution and LWC in Po Valley fog, Contributions to Atmospheric Physics, 71, 87–100, 1998.Willett, H. C.: Fog and Haze, their Causes, Distribution, and Forecasting, Mon. Weather Rev., 56, 435–468, 10.1175/1520-0493(1928)56<435:FAHTCD>2.0.CO;2, 1928.Zhang, J.: Fog droplet size distribution for Fog Case 1-1, Fog Case 4, Fog Case 10, Fog Case 18, and Fog Case 20, Zenodo [data set], 10.5281/zenodo.18552099, 2026.Zhang, J., Yang, Y., and Ding, J.: Information criteria for model selection, WIREs Comput. Stat., 15, e1607, 10.1002/wics.1607, 2023.