Parameterization of Radiation Fog-Top Height and Methods Evaluation in Tianjin

: Di ﬀ erent methods have been developed to estimate the fog-top height of radiation fog and evaluated using the measurements obtained from a 255-m meteorological tower located in Tianjin in 2016. Di ﬀ erent indicators of turbulence intensity, friction velocity ( u * ), turbulence kinetic energy ( TKE ), and variance of vertical velocity ( σ w2 ) were used to estimate the fog-top height, respectively. Positive correlations between the fog-top height and u * , TKE , and σ w2 were observed, with empirical parameterization schemes H = 583.35 × u 1.12 ∗ , H = 205.4 × ( TKE ) 0.68 , and H = 420.10 × (cid:16) σ 2 w (cid:17) 0.51 being obtained. Among them, σ w2 is the most appropriate indicators of turbulence intensity to estimate the fog-top height. Compared with sensible ﬂux and condensation rate, the new form of convective velocity scale (w * ) was the most appropriate indicator of buoyancy induced by radiative cooling, and the relationship H = 328.33 × w 1.34 ∗ was obtained. σ w2 and with w * , which represents the intensity of turbulence and buoyancy, were used to estimate the fog-top height. The relationship H = 396.26 × ( σ w + 0.1 × w * ) − 16 was obtained, which can be used to accurately estimate the fog-top height. Moreover, the temperature convergence (TC) method was used to estimate the fog-top height; however, the results strongly rely on the threshold value.


Introduction
Fog is a boundary layer weather phenomenon composed of suspended water droplets or ice crystals, which reduces horizontal visibility to less than 1 km [1]. The low visibility due to fog events is a hazard to navigation, motor-vehicle accidents [2], plane delays [3], air crashes, and public health [4,5], and the impact of fog has significantly increased because of increasing air, marine, and road transportation [6,7]. Therefore, more attention needs to be paid on fog.
Previous studies have shown that fog and low stratus significantly influence the earth's radiation balance [8,9]. Fog-top height (fog thickness) is considered to be very useful information for aircraft maneuvers, especially important for aircrafts having to land in foggy conditions. Moreover, there is no doubt about the importance of an accurate information of fog thickness for data assimilation of Numerical Weather Prediction (NWP) models, due to the significant impact of this parameter on the radiation budget close to the surface [10]. The estimation of observed fog-top height is important and useful for validation of model simulations [11], since comparisons between observed and simulated fog thickness cannot be performed in many cases due to no supplementary observational evidence of fog-top height. Moreover, the height of the fog top at the end of the mature phase is useful information for estimating the beginning and duration of the fog dissipation phase [12,13], with the thicker the fog layer, the more time is needed for fog dissipation. Therefore, it is also crucial to improve the nowcasting of fog dissipation [14]. However, many studies cannot provide information about observed fog-top height due to a lack of vertical measurements, and the height of fog top is usually difficult to determinate, since the fog-layer top cannot be observed from the surface.
The vertical extent of fog (fog thickness) ranges between a few meters and several hundred meters. Previous studies showed that turbulent mixing and radiative cooling are thought to be responsible for the increase in the depth of the fog. Numerical results showed that the fog-top height depends on the intensity of turbulence; the stronger the turbulence, the higher the top of the fog. In addition, some researchers have pointed out that the depth of fog is mainly due to the intensity of the inversion layer and that the inversion layer plays an important role in suppressing the vertical extension of this fog event; that is, the stronger the inversion, the lower the height of the fog-top. Moreover, both the concentration and size of aerosol particles in the polluted fog can affect urban fog thickness [15]. Therefore, fog-top height is closely related to turbulence and the upper inversion. However, the quantitative relationships between fog-top height and intensity of turbulence and inversion are ambiguous in the present study [16].
Many instruments and methods have been developed to determine the height of fog; satellite data, ground remote sensing instruments, and atmospheric sounding are usually used to provide approximations of fog-top height. In many case, temperature and humidity data from atmospheric sounding are used to estimate fog thickness [17,18]. The relative humidity profile is the most common and simple method to estimate the height of fog top; however, it is still uncertain to precisely define the fog top using the relative humidity threshold because of measurement uncertainties [19]. Cai et al. [20] defined the height above with a relative humidity that was smaller than 80%, while Sorli et al. [21] set the threshold of relative humidity to 95%, and 100% was used in Huang et al. [22]. Moreover, a mixing ratio of cloud water has also been used to determinate the height of fog top, with a threshold value 0.05 g kg −1 [23]. By considering the inversion intensity of temperature profiles and the 100% threshold on relative humidity profiles, the fog height can be estimated relatively accurately [24]. Moreover, the height of the fog top can be determined using virtual temperature profiles [11]. However, to determine fog-top height accurately, the measure points under and above the fog top should be dense. Therefore, there is discrepancy between the height of fog top estimated by atmospheric sounding and the real value, due to the vertical resolution of atmospheric sounding data. Moreover, ground remote sensing instruments are used to study radiation fog [25]; as an active sensor, it has advantages of high vertical resolution and detection sensitivity, however, the results are not satisfactory. A method of fog detection has been proposed by Wu et al. [26] and this method seems to be more reliable than the passive satellite measurements; however, it fails to identify dense fog layers because of the high signal attenuation.
Data and products from satellite have been widely used to detect fog or low clouds in foggy conditions [27]. Most applications retrieve fog-top height by comparing the satellite-derived infrared (IR) (infrared) temperature with temperature of soundings nearby. The premise of the scheme is that fog is related to an uplifted temperature inversion with inversion base height delimiting the fog top; therefore, the inversion base height is the fog top height [28]. Yi et al. [29] presented a novel method to retrieve low stratus/fog top heights using the infrared (IR) water vapor and split-window bands from geostationary (GEO) systems. However, the equations to retrieve inversion and thickness under foggy conditions were taken from Liu and Key [30], which were originally developed for clear sky situation. Thus, in future work, the adjustment of the new scheme should be developed to improve the adaptation of the original scheme. Despite its merit, the study also showed that a proper discrimination between low stratus and fog from satellite data remain challenging. A discrimination between fog at the ground and other low stratus situations from satellite data requires information on cloud vertical geometry to establish whether the cloud touches the ground. Cermak and Bendix [31] introduced a technique that allows for the discrimination between low stratus and (ground) fog on the basis of geostationary satellite imagery. The cloud base height is derived using a sub-adiabatic model of cloud microphysics. In this model, the cloud base is varied until model liquid water path (LWP) matches that retrieved from satellite data. The performance of this technique was shown to be good in comparison with METeorological Aerodrome Reports (METAR) data comprising 1030 satellite scenes. The main challenge in model development is to accurately quantify the deviation from the adiabatic profile. However, there are still misclassifications between low stratus and fog in some situations. In conclusion, it is very different to differentiate between fog and low clouds using satellite data [29,32].
As Marchand et al. [28] stressed in their review, fog-top height retrievals are a hitherto not yet a solved problem. In addition, satellite data, ground remote sensing instruments, and atmospheric sounding are expensive and not always available. Román-Cascón et al. [33] proposed two methods based on the surface data, and line correlations were observed between the height of fog top and surface friction velocity, and surface sensible flux. However, the results showed a relatively large bias and their methods were thus not appropriate to estimate the height of shallow fog.
In this study, different methods for estimation of fog-top height were evaluated and validated using the profile and atmospheric turbulence data obtained from a 255-m meteorological tower located in Tianjin. Moreover, two new estimated methods using the surface measurements were developed and evaluated. In Section 2, the information of the observation site, the instruments and data set are described. Moreover, several important parameters used in the study are calculated. In Section 3, different methods are used to estimate the fog thickness and compare with the observation results. A few new estimated methods based on turbulence intensity and radiative cooling are developed and validated using radiation fog events. In Section 4, the conclusions are summarized.

Experiment Site and Data
Tianjin is located in the eastern of North China Plain (NCP), with rather flat terrain. Furthermore, it is the largest coastal city in North China, with the Bohai Sea to the east. The results of Quan et al. [34] have shown that the mean occurrence number of fog events over the NCP is 11.4 ± 6.9 days yr −1 . Long-term fog trends are observed at 13 weather stations in Tianjin from 1980 to 2010, showing that the total amount of fog has increased since 1980 [35]. The data used in this study were obtained from a meteorological tower at a height of 255 m, located at the atmospheric boundary layer Meteorological Observation Station in the south of Tianjin city (39.08 • N, 117.21 • E). There are no tall buildings around the tower in a radius of 50 m, and heights of the surrounding buildings in the radius of 300 m are within 30 m. Therefore, the observation site can represent the homogeneous urban underlying surface [36].
The data and instruments used in this study are presented in Table 1. The measurements of meteorological parameters (wind speed, wind direction, air temperature, etc.) were recorded automatically and continuously at a 10-min interval. The 10 min data were then converted to a 30 min moving average. The sampling frequency of atmospheric turbulence data is 10 Hz, and a preprocessing was performed using Eddy Pro software (Advanced 4.2.1, LI-COR Biosciences, Inc. LI-COR, Lincoln, NE, USA). The preprocessing includes spike removal [37], double coordinate rotation [38], and trend removal. In addition, the block time average method was used in the trend removal process [39], with an average time interval of 30 min. The results of Eddy Covariance (EC) depend on wind speed, stability parameter, and friction velocity; strict quality control was performed on the atmospheric turbulence data, for detailed criteria, please refer to Ye et al. [40] and Ren et al. [41]. Moreover, it should be noted that due to the limitation of the maximum height of the observation tower and because the fog-top heights of radiation fog were typically below 200 m [42], fog-top heights that were higher than 250 m are out of the scope of this paper.

Calculation Methods
Turbulence kinetic energy (TKE) and friction velocity (u * ) can be calculated according to Equations (1) and (2), respectively: where u , v , and w are fluctuating values respect to average values of the two horizontal wind speed u and v and vertical wind speed w (m s −1 ). The turbulent exchange coefficient K m (m 2 s −1 ) can be calculated according to Equation (3): where κ is the von Karman's constant and has a value of 0.4; z is the reference level and the observational height (m); L is the Obukhov length (m); z/L represents the stability factor; ϕ m is a dimensionless wind shear [43,44]. The sensible heat (H c , w m −2 ) is calculated as follows: where ρ is the air density (kg m −3 ), c p the specific heat of air (J kg −1 K −1 ), T is the fluctuation of air temperature (K), and T * is temperature scale (K). The parameter critical turbulent exchange coefficient (Kc) (m 2 s −1 ) can be calculated according to Equations (5) and (6), for a shallow fog and a deep fog, respectively: Atmosphere 2020, 11, 480 5 of 18 where α = 0.062 is the gravitational settling parameter for radiation fog [45], β(p,T)C 0 (z) is the condensation rate per unit mass due to cooling of the air, and C 0 (z) = −(∂T/∂t) is the total local cooling rate, hereinafter referred to as cooling rate. The slope β(p,T) can be expressed using the Clausius-Clapeyron equation as: where Lv is the latent heat of vaporization of water (J kg −1 ), e s represents the saturation vapor pressure (h pa), and R v the gas constant for vapor (J kg −1 K −1 ). C t is the cooling rate at the fog top (K h −1 ), H is the depth of the fog bank (m). Moreover, the parameter can be used as a threshold; when K m < K c , the fog can persist at a fog mature phase [46]. The convective velocity scale (w * ) represents the scaling velocity for the convective boundary layer, and can be calculated according to Equation (8) as follows: where z i is the convective boundary layer height and can be estimated using the revised Carson scheme, by setting the initial boundary layer height to 150 m according to the potential temperature profile. The term (g/θ)(w θ ) represents the buoyancy heat flux. The liner correlation between fog thickness and u * is obtained in Román-Cascón et al. [33] and the relationship is:

Results and Discussions
There are no unified fog criteria in previous literature, and different criteria have been used to identify fog events [47][48][49]. In this study, the criteria of fog were as follows [40]: (1) RH ≥ 90%, (2) precipitation less than 1 mm per hour to eliminate the effect of precipitation (precipitation fog is out of our scope), and (3) the 10-min visibility ≤1 km. Moreover, we only focus on radiation fog in this study, which is the primary and most important inland fog form. Therefore, the 'fog' in this study refers to radiation fog, unless otherwise stated. All radiation fog events observed in 2016 were selected to develop new parameterizations of the fog-top height and evaluate different methods. The selected radiation fog events all formed in a nocturnal boundary layer and mainly due to radiative cooling at the surface. The wind speed was always low, and the durations of fog were between 7 and 18 h, which is similar to the statistical result in Liu et al. [50]. In this study, the fog-top height was estimated using the relative humidity profile (not shown here), and the fog thickness was determined as the height above where the value of RH was smaller than 90%, and under where the value of RH was larger than 90%.

Fog-Top Height Estimations Using Turbulence Intensity
Previous studies [19] have shown that the whole fog life cycle is linked to the turbulence intensity. Though the role of turbulence in the formation of radiation fog remains controversial, turbulence within the fog layer is crucial for the exchange of material and energy, important for the development of fog. Therefore, turbulence mixing is considered to be the main mechanism that influences the fog-top height, and different described parameters related to turbulence intensity are used to estimate the fog-top height. In this study, friction velocity (u * ), TKE, and variance of vertical velocity (σ w 2 ) were used to estimate the fog-top height, and the results were evaluated. The relationships between the fog-top (H) height and u * , TKE, and σ w 2 are shown in Figure 1a,b,c, respectively. A positive correlation between H and u * is shown in Figure 1a, and the relationship Atmosphere 2020, 11, 480 6 of 18 can be obtained, with a coefficient of determination of R 2 = 0.41. The result confirmed the hypothesis that the turbulence intensity is crucial in the life cycle of fog and is closely linked to the fog-top height.
Although there is a low turbulence intensity during the fog formation phase, turbulence mixing leads above air to be saturated conditions at the higher levels, and hence, to the development of fog bank, while strong turbulence leads to the fog collapse [12]. Therefore, the range of u * in the scatterplot (Figure 1a) can be interpreted as the minimum and maximum turbulence intensity required for fog formation and dissipation, which are approximately 0.037 and 0.281 m s −1 , respectively. The minimum turbulence intensity was slightly larger than that in Román-Cascón et al. [33], which were 0.025 and 0.030 m s −1 for different sites. The larger values of turbulence in this study were due to the sonic instrument installation height being 40 m in this study, which was 3 m or 1.5 m in Román-Cascón et al. [33]. However, the result verified that there was a threshold value of turbulence that influenced the formation of fog [51]. The maximum turbulence intensity in this study was much larger than those in Román-Cascón et al. [33]; this discrepancy was due to the sonic instrument installation height and the thickness of fog, with the thicker the fog, the stronger the turbulence the fog can endure [46]. The estimations with H = 1369 × u * − 28 (Equation (9)) obtained in Román-Cascón et al. [33] and similar linear correlation H = 502 × u * − 4 obtained in this study are also shown in Figure 1a. The results show there are obvious overestimations with H = 1369 × u * − 28, compared with the observations in Tianjin. There were few discrepancies between the estimations with H = 502 × u * − 4 and Equation (10), though the coefficient of determination of estimations with Equation (10) was a little higher. Therefore, Equation (10) obtained in this study seemed to be more appropriate to estimate the fog-top height in Tianjin, compared with Equation (9). The main reason of the discrepancies between the result in this study and Román-Cascón et al. [33] was that the relationship between fog-top height and u * in Román-Cascón et al. [33] was obtained by fitting seven or four mean values of u * associated to each discrete value of fog thickness, while the relationship in this study was obtained from a large amount of original values. In addition, the vertical resolution in Román-Cascón et al. [33] was low, seven layers under 200 m or four layers under 100 m, which led to largely misestimate the actual fog-top height. The vertical resolution in this study was higher, 15 layers under 250 m and eight layers under 100 m, which is favorable for the accurate estimation of fog-top height. Therefore, the relationship (Equation (10)) in this study was more appropriate to estimate the fog-top height, which represents the effect of turbulence on fog-top height and can show the temporal variations in the fog life cycle. However, there are some underestimations in large u * and overestimations for weak turbulence, which may be due to the effect of other meteorological parameters and large-scale atmospheric background fields. Similar to the result of u * , the estimations of fog-top height with TKE are shown in Figure 1b, and the relationship H = 205.43 × (TKE) 0.68 (11) with a coefficient of determination of R 2 = 0.44 was obtained, which indicates that TKE can be used to estimate the fog-top height. Compared with the estimations of the fog-top height with u * and TKE, the estimations of the fog-top height with σ w 2 seem to be more agreed with the observations (Figure 1c).
An obvious correlation between H and σ w 2 is shown in Figure 1c, and the relationship H = 420.10 × σ 2 w 0.51 (12) can be obtained, with a coefficient of determination of R 2 = 0.53. The results suggest that the variance of vertical velocity (σ w 2 ) is a more appropriate turbulence intensity indicator than the first two parameters, which can be used to estimate the fog-top height. Moreover, the minimum value of σ w 2 is 0.0031 m 2 s −2 in this study, that is, in the threshold range for fog formation (0.002-0.005 m 2 s −2 ) obtained in Price [51], can prove the existence of a threshold value of turbulence for fog formation. Though there are some biases between estimations with turbulence and observations, the results in this section suggest the possibility of using turbulence to estimate the fog-top height.
Though there are some biases between estimations with turbulence and observations, the results in this section suggest the possibility of using turbulence to estimate the fog-top height.

Fog-Top Height Estimations using Radiative Cooling
Radiative cooling and heating impact the liquid water balance of fog, and therefore, play an important role in the evolution of fog. There are two main radiative processes that affect the evolution of fog: one is the longwave radiative cooling at the fog top, which produces liquid water by condensation; the other one is the heating of the ground by absorption of solar radiation, which can cause a sensible heat transfer to the fog, causing the fog dissipation [45]. Therefore, radiation fog often occurs at night, and dissipates a few hours after sunrise due to the solar radiation and turbulence. Therefore, condensation rate related to radiative cooling should to be closely related to fog-top height. However, the aim of this study was to estimate fog-top height only using the surface measurements, and previous result confirmed that an average or surface cooling rate can be used to replace the cooling at the fog top [46]. The condensation rate near the surface was used to estimate the fog-top height in this study; the result is shown in Figure 2a. The result shows that the fog-top height seems to be less related to the condensation rate near the surface, though fog optical depth is closely related to the condensation rate. Moreover, the relationship between the difference of fog-top height and condensation rate is shown in Figure 2b. The results show that the fog-top height increases (difference of fog-top height is larger than zero) when the condensation rate is larger than zero, while fog-top height decreases when condensation rate is smaller than zero. However, there are some cases where those are opposite to the conclusion that radiative cooling supports the development of fog. The main reason for this may be radiative heating and the difference of condensation rate between surface and fog top. Moreover, the result indicates that turbulence mixing seems to be more important for the development of fog-top height, though the condensation rate is crucial to the formation of fog and fog optical depth. Previous studies show that the fog-top height is influenced by the buoyancy generated by radiative cooling at the fog top, other than the turbulence. Therefore, sensible flux was chosen as the scaling parameter to represent the intensity of buoyancy and used to estimate the fog-top height. Compared with the estimations with turbulence intensity, the relationship between the fog-top height and sensible flux is complicated, and the scatter points are discrete, especially when the values of sensible flux are near zero (not shown here). It is obvious that the relationship between the fog-top height and sensible flux can be divided into two phases: the positive and negative sensible flux. Therefore, the relationship between the fog-top heights and absolute values of sensible flux is developed, and Equation (13) is obtained, with a coefficient of determination of R 2 = 0.40 (Figure 3a): Previous studies show that the fog-top height is influenced by the buoyancy generated by radiative cooling at the fog top, other than the turbulence. Therefore, sensible flux was chosen as the scaling parameter to represent the intensity of buoyancy and used to estimate the fog-top height. Compared with the estimations with turbulence intensity, the relationship between the fog-top height and sensible flux is complicated, and the scatter points are discrete, especially when the values of sensible flux are near zero (not shown here). It is obvious that the relationship between the fog-top height and sensible flux can be divided into two phases: the positive and negative sensible flux. Therefore, the relationship between the fog-top heights and absolute values of sensible flux is developed, and Equation (13) is obtained, with a coefficient of determination of R 2 = 0.40 (Figure 3a): Though the result indicates that the absolute values of sensible flux can be as a potential indicator of the fog-top height, there are some discrepancies between estimations and observations, especially for low sensible flux. The results are similar to the results in Román-Cascón et al. [33], with significantly dispersed distribution and larger error bars. Therefore, sensible flux may be not the appropriate indicator to represent the buoyancy and other described parameters should be sought out to represent the buoyancy induced by radiative cooling.
Once fog occur, fog stratification underwent a transition from thermally stable to weakly unstable, and the fog layer can be deemed to be well-mixed. To represent the intensity of buoyancy induced by radiative cooling at the fog top, a new form of the convective velocity scale (w * ) was introduced. It should to be noted that new w * in this study is calculated according to Equation (8) with the absolute value of sensible flux. The new w * represents the buoyancy and neglects the sign, which means whether the sensible flux is negative or positive, it contributes to the mixing of fog layer and dry air, and leads to the development of fog layer.
The relationships between the H and new w * are shown in Figure 3b and a positive correlation between H and w * is observed, with the relationship being obtained, with a coefficient of determination of R 2 = 0.44. The results confirm that buoyancy induced by radiative cooling is a contributor to fog development, with the mixing between a fog layer and dry air. There are some bias between estimations and observations, which is because the fog bank near the surface is not always well-mixed fog; therefore, the effect of radiative cooling at the fog top cannot be represented by the new w * . However, compared with sensible flux, the new w * is a more appropriate indicator of buoyancy induced by radiative cooling at the fog top, which can be used to calculate the fog-top height.
Atmosphere 2020, 11, x FOR PEER REVIEW 10 of 19 Once fog occur, fog stratification underwent a transition from thermally stable to weakly unstable, and the fog layer can be deemed to be well-mixed. To represent the intensity of buoyancy induced by radiative cooling at the fog top, a new form of the convective velocity scale (w * ) was introduced. It should to be noted that new w * in this study is calculated according to Equation (8) with the absolute value of sensible flux. The new w * represents the buoyancy and neglects the sign, which means whether the sensible flux is negative or positive, it contributes to the mixing of fog layer and dry air, and leads to the development of fog layer.
The relationships between the H and new w* are shown in Figure 3b and a positive correlation between H and * is observed, with the relationship H = 328.33 × w * 1.34 (14) being obtained, with a coefficient of determination of R 2 = 0.44. The results confirm that buoyancy induced by radiative cooling is a contributor to fog development, with the mixing between a fog layer and dry air. There are some bias between estimations and observations, which is because the fog bank near the surface is not always well-mixed fog; therefore, the effect of radiative cooling at the fog top cannot be represented by the new * . However, compared with sensible flux, the new * is a more appropriate indicator of buoyancy induced by radiative cooling at the fog top, which can be used to calculate the fog-top height.

Fog-Top Height Estimations Using Turbulence and Radiative Cooling
Based on the budget of liquid water content (LWC), the asymptotic LWC distribution was obtained, which is the result of the comprehensive effect among radiative cooling, droplet gravitational settling, and turbulence mixing in the liquid water budget of radiation fog. Based on the relationship between turbulence exchange coefficient and the fog thickness (Equation (6)), a new estimation method was developed: where 'a' represents the coefficient constant and the result is shown in Figure 4. The relationship H = 1.45 k 2 m β(p,T)C t 1 3 + 35 was obtained, with a coefficient of determination of R 2 = 0.26, and the scatter points were extremely discrete. Obvious underestimations were observed; thus, this method was not appropriate to estimate the fog thickness. The large bias between the estimations and observations may be because Equation (15) was obtained from the balance condition at the fog mature phase, which is not suitable for estimations of the whole fog life cycle. In conclusion, there are two main causes of the poor estimations, one is that the asymptotic formulation is only suitable for the beginning of dissipation phase, and the other cause is that the relationship between condensation rate and the fog-top height is complicated and influenced by other factors, such as radiative heating and droplet gravitational settling.

Fog-Top Height Estimations Using Turbulence and Radiative Cooling
Based on the budget of liquid water content (LWC), the asymptotic LWC distribution was obtained, which is the result of the comprehensive effect among radiative cooling, droplet gravitational settling, and turbulence mixing in the liquid water budget of radiation fog. Based on the relationship between turbulence exchange coefficient and the fog thickness (Equation (6) were extremely discrete. Obvious underestimations were observed; thus, this method was not appropriate to estimate the fog thickness. The large bias between the estimations and observations may be because Equation (15) was obtained from the balance condition at the fog mature phase, which is not suitable for estimations of the whole fog life cycle. In conclusion, there are two main causes of the poor estimations, one is that the asymptotic formulation is only suitable for the beginning of dissipation phase, and the other cause is that the relationship between condensation rate and the fog-top height is complicated and influenced by other factors, such as radiative heating and droplet gravitational settling. Previous studies have shown that the development of fog can be attributed to turbulence mixing and radiative cooling at the fog top. In brief, the life cycle of fog depends mainly on the balance between radiative cooling and turbulence. Moreover, the results in Section 3.1 and 3.2 suggest that σw 2 and * , with the same dimension, are the most appropriate parameters to represent the intensity of turbulence and buoyancy induced by radiative cooling, which are closely related to the fog-top Previous studies have shown that the development of fog can be attributed to turbulence mixing and radiative cooling at the fog top. In brief, the life cycle of fog depends mainly on the balance between radiative cooling and turbulence. Moreover, the results in Sections 3.1 and 3.2 suggest that σ w 2 and w * , with the same dimension, are the most appropriate parameters to represent the intensity of turbulence and buoyancy induced by radiative cooling, which are closely related to the fog-top height. Therefore, the fog-top height is estimated using the comprehensive parameters that include the two physics parameters. The relationship can be obtained, with a coefficient of determination of R 2 = 0.55. Some improvements can be observed ( Figure 5), which confirms the previous conclusion that the fog life cycle is closely related to the turbulence intensity and radiative cooling. In conclusion, Equation (16) can be used to quantitatively estimate the fog-top height, using only the surface measurements. However, the estimation with Equation (16)  can be parameters to estimate the fog-top height. More scaling parameters represent the turbulence and radiation and different combinations of these scaling parameters should be found to estimate the fog-top height. Moreover, more radiation fog events are required to be used to estimate the fog-top height and verify the results.

Fog-Top Height Estimation through Temperature Convergence (TC Method)
The radiative cooling near the surface can lead to saturated air conditions at the ground, and hence, to the formation of fog droplets, and it may promote condensation in a supersaturated surface layer of sufficient depth, leading to fog formation. Once the fog first forms at the ground, the fog stratification undergoes a transition from thermally stable to weakly unstable [52]. The turbulent mixing within the fog causes the homogenization of the layer, and the mixing causes the convergence of temperatures, which leads the temperature of different heights to be a constant [11]. Previous studies have shown that the fog thickness can be estimated through the vertical temperature profile [53,18]. Considering the fog to be present at a certain height, when the difference between the potential temperature (θz) at the height and surface level (θs) is less than a threshold value (θT): The fog-top height is estimated to be the maximum height under where this criterion (in Equation (17)) is always satisfied. With the use of potential temperature, the height-related influence can be avoided. The determination of threshold value is crucial to estimate the fog thickness accurately. The uncertainty between two level temperatures is 0.4 K, because the instrument-related uncertainty of temperature measurements is 0.2 K. Small differences in potential temperatures are allowed between different heights, and this strategy has been used in Price [52] and Román-Cascón

Fog-Top Height Estimation through Temperature Convergence (TC Method)
The radiative cooling near the surface can lead to saturated air conditions at the ground, and hence, to the formation of fog droplets, and it may promote condensation in a supersaturated surface layer of sufficient depth, leading to fog formation. Once the fog first forms at the ground, the fog stratification undergoes a transition from thermally stable to weakly unstable [52]. The turbulent mixing within the fog causes the homogenization of the layer, and the mixing causes the convergence of temperatures, which leads the temperature of different heights to be a constant [11]. Previous studies have shown that the fog thickness can be estimated through the vertical temperature profile [18,53]. Considering the fog to be present at a certain height, when the difference between the potential temperature (θ z ) at the height and surface level (θ s ) is less than a threshold value (θ T ): Atmosphere 2020, 11, 480 12 of 18 The fog-top height is estimated to be the maximum height under where this criterion (in Equation (17)) is always satisfied. With the use of potential temperature, the height-related influence can be avoided. The determination of threshold value is crucial to estimate the fog thickness accurately. The uncertainty between two level temperatures is 0.4 K, because the instrument-related uncertainty of temperature measurements is 0.2 K. Small differences in potential temperatures are allowed between different heights, and this strategy has been used in Price [52] and Román-Cascón et al. [33], who used threshold value settings of 0.8 K and 1.2 K, respectively. In this study, different threshold values were used in the estimations, and the threshold value was set to 1.2 K based on the estimated results of the fog thickness. The results of the fog thickness estimation versus the observations are shown in Figure 6a. The estimation results are consistent with the observations, with a coefficient of determination of R 2 = 0.61, which indicates that the TC method based on the differences in potential temperature is suitable for estimating the fog thickness, though there are some underestimations. Similar with the theory mentioned in Section 3.1, though the fog stratification is weakly unstable during the fog life cycle, atmospheric boundary structure may be different during the development phase and dissipation phase. Therefore, to estimate the fog thickness accurately, different threshold values were used, and the least square method was used to determine the θ T during the two phases. Ultimately, the threshold value during fog development phase and dissipation phase was set to 1.2 K and 1 K, respectively. The relationship between the fog thickness estimated by TC method and the observation are divided into two phases: the development phase and the dissipation phase; the results are shown in Figure 6b. There was only a little improvement with the segmentation of threshold values, with a coefficient of determination of R 2 = 0.64; thus, the same threshold value can be used for the whole fog life cycle. The TC method was based on the theory that turbulent mixing within fog causes the temperature convergence, which results in the same value of different heights. However, the physics mechanisms during fog development and dissipation phase are different, with radiative cooling of the fog top and turbulence mixing leading to the increase of fog top, while the fog dissipation is the comprehensive effect of turbulence, radiative cooling of fog top, and entrainment of dry air. There may be a little difference in the degree of temperature convergence during the two phases; however, this is not observed in this study, which suggests more fog events are required to improve the estimations. during the two phases. Ultimately, the threshold value during fog development phase and dissipation phase was set to 1.2 K and 1 K, respectively. The relationship between the fog thickness estimated by TC method and the observation are divided into two phases: the development phase and the dissipation phase; the results are shown in Figure 6b. There was only a little improvement with the segmentation of threshold values, with a coefficient of determination of R 2 = 0.64; thus, the same threshold value can be used for the whole fog life cycle. The TC method was based on the theory that turbulent mixing within fog causes the temperature convergence, which results in the same value of different heights. However, the physics mechanisms during fog development and dissipation phase are different, with radiative cooling of the fog top and turbulence mixing leading to the increase of fog top, while the fog dissipation is the comprehensive effect of turbulence, radiative cooling of fog top, and entrainment of dry air. There may be a little difference in the degree of temperature convergence during the two phases; however, this is not observed in this study, which suggests more fog events are required to improve the estimations. Figure 6. Scatterplots of observations of fog thickness against the estimation of fog thickness by temperature convergence (TC) method using in Equation (17), with (a) threshold value is set to 1.2 K, (b) threshold value during fog development phase and dissipation phase is set to 1.2 K and 1 K, respectively.
The virtual potential temperature seems to be more appropriate, because this parameter includes the effect of water vapor and liquid water content on temperature, which is more reasonable to describe the vertical distribution of temperature. Therefore, similar to in Equation (17), the difference between the virtual potential temperature at the height and surface level is used to determinate the fog top height. However, non-substantial differences are observed between using the potential temperature and virtual potential temperature (not shown here), which is consistent with the results of Román-Cascón et al. [33]. Moreover, previous studies also show that although some improvements have been made by using different thresholds for certain fog thickness estimations, it is always at the expense of a worsening of the results for other fog thickness; therefore, the determination of the threshold is still an intractable problem. Moreover, temperature convergence does not occur for Figure 6. Scatterplots of observations of fog thickness against the estimation of fog thickness by temperature convergence (TC) method using in Equation (17), with (a) threshold value is set to 1.2 K, (b) threshold value during fog development phase and dissipation phase is set to 1.2 K and 1 K, respectively.
The virtual potential temperature seems to be more appropriate, because this parameter includes the effect of water vapor and liquid water content on temperature, which is more reasonable to describe the vertical distribution of temperature. Therefore, similar to in Equation (17), the difference between the virtual potential temperature at the height and surface level is used to determinate the fog top height. However, non-substantial differences are observed between using the potential temperature and virtual potential temperature (not shown here), which is consistent with the results of Román-Cascón et al. [33]. Moreover, previous studies also show that although some improvements have been made by using different thresholds for certain fog thickness estimations, it is always at the expense of a worsening of the results for other fog thickness; therefore, the determination of the threshold is still an intractable problem. Moreover, temperature convergence does not occur for shallow fog (not shown here), which is linked to strong thermal inversions and weak turbulence; therefore, the TC method is not valid. The result suggests that although the TC method may be inadequate for estimating the fog thickness of shallow fogs, it is still an appropriate method to estimate the deepen fog thickness larger than 100 m, which is hazardous for transport.

Comparisons of Fog-Top Height Estimations with Different Methods
All parameterization schemes developed in this study and corresponding coefficients of determination are shown in Table 2. The results confirm that σ w 2 and w * , with the same dimension, are the most appropriate indicators to represent the intensity of turbulence and buoyancy induced by radiative cooling. Therefore, the combination of σ w 2 and w * (Equation (16)) is the best parameterization scheme developed in this study that can be used to estimate the fog-top height only using the surface measurement. Moreover, the TC method can be used to accurately estimate the fog-top height using the vertical profile data. The fog-top height estimations of a thick fog and shallow fog with H = 1369 ×u * − 28 obtained in Román-Cascón et al. [33] (Equation (9)), H =583.35 × u 1.12 * (Equation (10)), H = 396.26 × (σ w + 0.1 × w * ) − 16 (Equation (16)) developed in this study, and the TC method (in Equation (17)) were compared and are shown in Figure 7. The results show that there were obvious overestimations with the relationship (Equation (9)) obtained in Román-Cascón et al. [33], for both thick and shallow fog, while the other three estimation methods all perform better in thick fog than in shallow fog. The results in Section 3.4 has pointed out that the TC method is not suitable for estimating the fog-top height of shallow fog, in which temperature convergence does not occur. Comparing the estimations in shallow fog with only u * (Equation (10)) and the combination of u * and w * (Equation (16)), the estimations with the combination of u * and w * are closer to the observations, since the new parametrization scheme (Equation (16)) in this study considers the effect of radiative cooling and turbulence. Moreover, the new parametrization scheme (Equation (16)) can approximately simulate the temporal variation of fog-top height in the fog life cycle, though there are some underestimations. For thick fog, the TC method performed best, especially during fog formation and development phase. However, the TC method is expensive and not available in most areas. Moreover, to determine fog-top height with the TC method accurately, the measure points under and above the fog top should be dense. Therefore, the TC method can only be applied in a few sites, and the results strongly rely on the set of threshold value and vertical resolution of profile data. Comparing the estimations in thick fog with only u * and the combination of u * and w * , there are better agreements between the estimations of the new parametrization scheme (Equation (16)) and observations, especially during the fog development phase. The new parametrization scheme can also simulate the temporal variation of fog-top height during the fog dissipation phase, though there are some underestimations. In conclusion, if only surface measurements are used to estimate the fog-top height in Tianjin, the new parametrization scheme (Equation (16)) developed by combining the effect of turbulence and radiative cooling in this study is more suitable compared with other methods. However, the parametrization schemes, except the TC method, cannot used to simulate the fog-top height during fog formation phase, during which the effect of turbulence on fog is ambiguous. Therefore, more effort should be made to improve our understandings of the effect of turbulence on fog life cycle.  (17)) during a shallow and thick fog events in 2018.

Conclusions
In this study, different estimation methods for the fog-top height were developed and evaluated using all radiation fog events obtained from a 255-m meteorological tower located in Tianjin in 2016. The height of fog top was estimated using the relative humidity profile, and the threshold value of RH is 90%.
Different indicators of turbulence intensity, friction velocity (u*), TKE, and variance of vertical velocity (σw 2 ) were used to estimate the fog-top height, respectively, and the results were evaluated. The positive correlations between H and * and TKE were obtained, with empirical parameterization scheme H = 583.35 × u * 1. 12 and H = 205.43 × TKE 0.68 . Compared with * and TKE, the variance of vertical velocity (σw 2 ) was the most appropriate turbulence intensity indicator for estimating the fog-top height, with the relationship H = 420.10 × σ w 2 0.51 and a higher coefficient of determination being obtained. The results confirm that turbulence intensity is crucial to the life cycle of fog and is closely linked to the fog-top height. Moreover, the result confirmed the hypothesis that there are threshold values of turbulence that can influence the fog formation and dissipation.
Radiative cooling impacts the liquid water balance of fog and therefore plays an important role in the evolution of fog. Different indicators of radiative cooling, surface condensation rate, and sensible flux were used to estimate the fog-top height. The result in this study show that there is no obvious relationship between surface condensation rate and radiative cooling, which indicates that the condensation rate near the surface is not suitable to estimate the fog-top height, though it is positively related to fog optical depth. The estimation results with sensible flux, which was chosen as the described parameter representing the buoyancy generated by radiative cooling at the fog top, are discrete, with large biases. Consider the effect of buoyancy on vertical mixing, the absolute values of the sensible flux were used to associate with the fog-top height, and relationship H = 29.49 × |Hc| 0.45 was obtained. However, the coefficient of determination was low (R 2 = 0. 40 (9)), H =583.35 × u 1.12 * (Equation (10)), H = 396.26 × (σ w + 0.1 × w * ) − 16 (Equation (16)), and the TC method (in Equation (17)) during a shallow and thick fog events in 2018.

Conclusions
In this study, different estimation methods for the fog-top height were developed and evaluated using all radiation fog events obtained from a 255-m meteorological tower located in Tianjin in 2016. The height of fog top was estimated using the relative humidity profile, and the threshold value of RH is 90%.
Different indicators of turbulence intensity, friction velocity (u * ), TKE, and variance of vertical velocity (σ w 2 ) were used to estimate the fog-top height, respectively, and the results were evaluated. and a higher coefficient of determination being obtained. The results confirm that turbulence intensity is crucial to the life cycle of fog and is closely linked to the fog-top height. Moreover, the result confirmed the hypothesis that there are threshold values of turbulence that can influence the fog formation and dissipation.
Radiative cooling impacts the liquid water balance of fog and therefore plays an important role in the evolution of fog. Different indicators of radiative cooling, surface condensation rate, and sensible flux were used to estimate the fog-top height. The result in this study show that there is no obvious relationship between surface condensation rate and radiative cooling, which indicates that the condensation rate near the surface is not suitable to estimate the fog-top height, though it is positively related to fog optical depth. The estimation results with sensible flux, which was chosen as the described parameter representing the buoyancy generated by radiative cooling at the fog top, are discrete, with large biases. Consider the effect of buoyancy on vertical mixing, the absolute values of the sensible flux were used to associate with the fog-top height, and relationship H = 29.49×|Hc| 0.45 was obtained. However, the coefficient of determination was low (R 2 = 0.40), which indicates the effect of radiative cooling at the fog top on fog-top height cannot be represented by sensible flux. To represent the buoyancy generated by radiative cooling at the fog top, a new form of convective velocity scale (w * ), which was calculated with the absolute value of the sensible flux, was introduced and used to estimate the fog-top height. An obvious positive correlation between H and new w * was observed, and the relationship H =328.33 × w 1.34 * was obtained, with a higher coefficient of determination. The results show that the new w * is more appropriate to represent the effect of radiative cooling at the fog top and can be used to calculate the fog-top height.
To comprehensively consider the effect of turbulence and radiative cooling, different combinations of different parameters were used to estimate the fog-top height. Based on the asymptotic LWC distribution at the fog mature phase, which was the result of comprehensive effect among radiative cooling, droplet gravitational settling, and turbulence mixing in the liquid water budget of radiation fog, a new method (Equation (16)) was presented to estimate the fog thickness. However, the estimations were extremely discrete, which indicates that this method is not appropriate to estimate the fog-top height in the whole fog life cycle. σ w 2 and w * , which are the most appropriate indicators of the intensity of turbulence and buoyancy induced by radiative cooling, and are used to estimate the fog-top height. The relationship H = 396.26 ×(σ w + 0.1 ×w * ) − 16 was obtained, with a higher coefficient of determination, which can be used to accurately estimate the fog-top height.
The fog thickness was also estimated through the vertical temperature profile, based on the theory that the fog layer is homogeneous, and the temperature at different heights within the fog is approximately constant, due to the turbulent mixing within fog. The fog-top height was estimated based on Equation (17) and threshold value was set to 1.2 K. The results of fog thickness estimation using the TC method were consistent with the observations, with a coefficient of determination of R 2 = 0.61. Moreover, there was a little improvement with the segmentation of threshold values, with a coefficient of determination of R 2 = 0.64. The virtual potential temperature was also introduced to determine the fog-top height, and few differences were observed between using the potential temperature and virtual potential temperature. In addition, the determination of threshold value was still an intractable problem and the TC method is not valid for shallow fog due to strong thermal inversions and weak turbulence. In brief, the TC method based on the differences in the potential temperature is suitable for estimating the fog thickness of deep fog.
This study compared the estimations of fog-top height with H = 1369 u * − 28 obtained in Román-Cascón et al. [33], (Equation (9)), H =583.35 × u 1.12 * (Equation (10)), H = 396.26 × (σ w + 0.1 × w * ) − 16 (Equation (16)) developed in this study, and the TC method (in Equation (17)) in a thick fog and shallow fog, respectively. The results showed that TC method performs best in thick fog, though it is not valid in shallow fog, which is not well-mixed. However, the TC method is expensive and not available in most areas, and the results strongly rely on the set of threshold values, which are different in different sites. Moreover, to determine fog-top height accurately, the measure points under and above the fog top should be dense. Therefore, the TC method can only be applied in a few sites, and the results strongly rely on the set of threshold and vertical resolution of profile data. There are better agreements between the estimations of the new parameterization scheme (Equation (16)) developed in this study and observations in both thick and thin fog, compared with the estimations of fog-top height with relationship (Equation (9)) in Román-Cascón et al. [33] and with only u * (Equation (10)). The results confirmed that the new parametrization scheme (Equation (16)) developed in this study, considering the comprehensive effect of turbulence and radiative cooling, can be widely used to estimate fog-top height only using the surface measurements.
It is not possible to state definitive conclusions after the analysis of several radiation fog cases at one site; nevertheless, the main factors, turbulence intensity and radiative cooling, that influence the fog-top height were determined. However, due to the limitation of observation height and vertical resolution, more research is required to seek for more appropriate and common parameters that can be used to estimate the fog-top height. Funding: This research was jointly funded by the National Natural Science Foundation of China (41675018, 41675135, 41705045), the Natural Science Foundation of Tianjin (17JCYBJC23400), and the Bohai Rim Regional Fund (QYXM201801).

Conflicts of Interest:
The authors declare no conflict of interest.