According to statistics, in the 12 years from 2008 to 2019, the spatially average number of total fog days in Shandong Province was 218 (the spatially annual average was 18 days). The number of total fog days observed at most stations concentrated between 90 and 270 days. Eighty-five percent of the stations reported fewer than 300 fog days (Fig. 2a). However, some special stations observed more than 600 fog days: Taishan station at a high altitude of 1533.7 m reported 1054 days, and Chengshantou station reported 693 days. The spatially averaged number of light fog days was 74 (the spatially annual average was 6 days), accounting for 33.9% of the total fog days. Most stations reported between 35 and 110 light fog days, while a few stations experienced more than 600 fog days (Fig. 2b). The spatially averaged number of dense fog days was 53 (the spatially annual average was 4 days). The number of dense fog days observed at most stations concentrated between 25 and 70, and nearly 100% of all stations reported fewer than 150 dense fog days (Fig. 2c). The spatially averaged number of extreme fog days was 90 (the spatially annual average was 8 days). Most stations reported between 25 and 160 extreme fog days (Fig. 2d). Compared with light fog and dense fog days, extreme fog days were more widely distributed and more skewed with a larger standard deviation and a higher degree of dispersion.
3.1. The spatial distribution characteristics and the mechanisms
To analyze the spatial distribution of fog days in Shandong Province during 2008–2019, this study used the Cressman interpolation method to interpolate the statistical station data to a latitude and longitude grid with a spatial resolution of 0.1°×0.1°; accordingly, the distributions of the numbers of total fog, light fog, dense fog, and extreme fog days were obtained (Fig. 3). Because the altitude of Taishan station is 1533.7 m, almost within the cloud layer, most of the fog observed at this station may be low clouds. The observation data showed a high frequency of extreme fog days at this station; however, the formation mechanism of this low visibility differed from other regions. Therefore, the Taishan station was excluded from the analyses of the spatiotemporal distribution of fog. For the spatial distribution of fog in Shandong province, a few fog days region appeared a northeast-southwest trending zone in the central region, where was higher altitude topography, and more fog days on both sides can be found at lower altitude areas (Fig. 1 and Fig. 3). A few total fog days region located in the central region in Shandong province were about 130 days in the 12 years from 2008 to 2019, including the light fog of about 50 days (Fig. 3b). The numbers of total fog located in the west of Shandong province reached about 340 days, which included about 80 dense fog days and 140 extreme fog days (Figs. 3c and d). The coastal stations in eastern Shandong (such as Chengshantou and Shidao) reported high numbers of fog days (Fig. 3a), especially extreme fog days (Fig. 3d), which might be related to the adjacent ocean (Huang et al. 2011; Choi 2013; Pithani et al. 2019; Tian et al. 2020).
Shandong features unique terrain, and the elevation of central Shandong (Taishan Mountains) is higher than those of the surrounding areas, and the spatial distribution characteristics of fog in Shandong Province showed that the higher-elevation area in the center of the study area corresponded to low numbers of total fog, light fog, dense fog, and extreme fog days, and the lower-elevation areas on both sides exhibited high numbers (Fig. 1 and Fig. 3), which indicates that the formation of fog in Shandong Province might be related to the terrain (Fu et al. 2014).
On the other hand, the formation of fog requires a nearly saturated or supersaturated water vapor environment, which necessitates favorable meteorological conditions. Figure 4 shows the box-whisker plot of the average daily average RH and DTR at low-altitude (< 100 m) stations and high-altitude (> 100 m) stations (excluding Taishan station). The daily average RH was present on 65.0% at low-altitude stations and 62.5% at high-altitude stations (Fig. 4a). DTR was present on 9.7°C at low-altitude stations and 9.4°C at high-altitude stations (Fig. 4b). Compared with high altitude areas, low altitude areas had more abundant water vapor conditions and larger DTR, which was more conducive to the formation of fog in low altitude areas.
To quantify the respective effects of altitude, background RH, and DTR on the spatial distribution of fog frequency, the method from Liu et al. (2017) and Sharma et al. (2019) was used. Therefore, multivariate linear regression was established:
where a, b, and c are regression coefficients, and d is constant. F is fog days, and x1, x2, and x3 indicate altitude, background RH, and DTR of meteorological stations, respectively.
Then, the relative effects of the three factors are as follows:
where Δ is the standard deviation and In this way, the relative effects of three factors on the spatial distribution of fog frequency can be obtained.
The relative effects of the altitude, daily average RH, and DTR on the spatial distribution of fog were calculated in Table 1. The altitude had negative relative effects on the spatial distribution of fog days, while daily average RH and DTR had positive relative effects. Specifically, daily average RH was the dominant factor among the three factors on the spatial distribution of fog, which accounted for approximately 61.95% in total fog, 44.24% in light fog, 69.07% in dense fog, and 72.57% in extreme fog. And the relative effect of RH on its spatial distribution became more important with increasing fog levels. Altitude was the factor that accounted for approximately − 19.68% in total fog, -15.89% in light fog, -19.89% in dense fog, and − 22.64% in extreme fog, which means that high level fog corresponded to a relatively larger negative relative effect of altitude. With the enhancement of the level of fog, the effect of the DTR on the spatial distribution of fog days decreased, which accounted for approximately 39.87% in light fog, 11.04% in dense fog, and 4.79% in extreme fog. These indicated that the relatively low altitude, high daily average RH, and high DTR correspond to the high-frequency area of fog days, and these factors had different relative effects on the spatial distributions of different levels of fog.
Table 1
Relative effects (%) of altitude (F1), RH (F2), and DTR (F3) on the spatial distribution of fog (total, light, dense, and extreme fog) days in Shandong Province.
Classification
|
Light fog
|
Dense fog
|
Extreme fog
|
Total fog
|
F1
|
-15.89
|
-19.89
|
-22.64
|
-19.68
|
F2
|
44.24
|
69.07
|
72.57
|
61.95
|
F3
|
39.87
|
11.04
|
4.79
|
18.37
|
Under the influence of the uneven terrain, the development of the wind and water vapor field often determines the location and time of fog formation (Golding 1993; Walmsley et al. 1996). A comparison between fog days and non-fog days reveals that the wind field in central Shandong was a divergent anticyclone on non-fog days that passed through southern and western Shandong and then flowed from northern Shandong to the Bohai Sea; namely, the airflow flowed from the mainland to the Yellow Sea (Fig. 5b). Figure 5d shows that the water vapor flux was relatively small on non-fog days in Shandong Province, and there was no water vapor transport channel, so no water vapor convergence occurred in the whole area. In addition, for the entire study area, the wind speed on non-fog days was greater than that on fog days, which was not conducive to the formation of stable atmospheric stratification. Under these conditions, the water vapor in the whole region on non-fog days was far less than that on fog days (Figs. 5a and b), and these meteorological conditions were not conducive to fog formation in this area.
However, on foggy days, the northerly wind originating from the Bohai Sea was divided into two branches under the influence of the terrain (Fig. 5c). The first airflow flowed to the west of Shandong, whereas the second airflow flowed to the eastern sea area of the Shandong Peninsula, then flowed through the Yellow Sea to southern Shandong in the form of an anticyclone and finally turned southeasterly. Figure 5c shows that southern Shandong Province exhibited high water vapor flux, namely, a strong water vapor transport belt. Sufficient water vapor conditions were conducive to the formation of fog in southern Shandong and coastal areas. The second airflow carrying warm and humid air from the Yellow Sea along this water vapor transport channel, northeasterly cold air (the first airflow) from the Bohai Sea, and westerly (or northwesterly) wind from the foot of the Taihang Mountains converged along a line in western Shandong Province. The area corresponding to this convergence line was characterized by negative water vapor flux divergence (Fig. 5e); namely, water vapor converged in this area. Therefore, the second airflow carried water vapor to eastern Shandong Province, and water vapor accumulated under the convergence of the wind field. In addition, many studies have shown that the convergence of cold and warm winds easily forms frontal fog (Tardif and Rasmussen 2008, 2010; Guo and Guo 2016). Therefore, such wind field and water vapor conditions under terrain conditions were conducive to the formation of fog in western Shandong Province, constituting one of the reasons for the frequent occurrence of fog in this region. In the Taishan Mountains in central Shandong Province, the wind field was divergent, the water vapor flux was small, and the divergence of the water vapor flux was positive; namely, water vapor divergence occurred (Fig. 5e). Hence, the low water vapor caused by the wind field under the terrain was beneficial to form low-frequency fog in this area.
3.2. Interannual variation and distribution of fog and its impact factors
The numbers of fog days at each station in 2008–2019 were linearly fitted, and the annual change rates of total fog, light fog, dense fog, and extreme fog days were obtained. Overall, 88.5% of all stations in Shandong Province showed an increasing trend of total fog days; a greater increase was in the western region, where Chiping station had the fastest growth (6.1 d a− 1 (d a− 1 stands day year− 1)). The number of fog days in most central stations experienced slow growth, and only a few stations in the eastern coastal region had slight decreasing trends (Fig. 6a). Light fog was the main contributor to fog growth. Approximately 75.2% of all stations showed an increasing trend of light fog days. The growth rate in the western region was larger than that in the eastern region overall, and the fastest growth rate was observed at Guangrao station (3.1 d a− 1) (Fig. 6b). For dense fog, 68.0% of the stations displayed an increasing trend, and the fastest growth occurred at Linqing station (1.5 d a− 1). Thirty-two percent of the stations had a decreasing trend, and most of them were distributed in the coastal and central areas; among them, Shidao station (-1.4 d a− 1) had the fastest decline (Fig. 6c). For extreme fog, 67.9% of the stations reported an increasing trend; almost all the stations in the western region experienced increasing trends, some of which were growing rapidly, and Chiping station (3.2 d a− 1) reported the fastest growth rate. In contrast, the stations with decreasing trends were distributed in the central, southern, and coastal areas, and the station with the fastest decline was Chengshantou station (-3.0 d a− 1) (Fig. 6d). It should be pointed out that more than half of the stations passed the significance test of 5%. In addition, the haze days presented an increasing trend, which indicated that the air pollution in Shandong has been aggravating in recent years. The haze days increased at most sites, and the distribution of the growth rate of haze days was basically consistent with that of fog days (Fig. 6e). The growth rate of haze days at the western stations was also greater than that at the eastern stations. This shows that the background of aggravating air pollution may be an important factor for the gradual increase in the frequency of fog.
From 2008 to 2019, the trends of the spatially averaged annual numbers of total fog, light fog, dense fog, and extreme fog days in Shandong Province were similar and could be divided into three stages (Figs. 7a and b). The first stage was 2008–2012, during which the spatially averaged number of fog days remained relatively unchanged with a slight decrease over time; this stage had the fewest annual numbers of fog days during the whole research period. The spatially averaged annual occurrences of light fog and dense fog were maintained at approximately 4 days, that of extreme fog was maintained at approximately 6 days and slightly decreased, and that of total fog occurred approximately 13 days with a slightly decreasing trend (Fig. 7a). The slope of the linear fit of the corresponding spatially averaged 3-year moving average in this stage was approximately 0 (Fig. 7b). The second stage was 2012–2016, during which the spatially averaged annual number of fog days increased rapidly. The number of total fog days increased by 125% from 12 days in 2012 to 27 days in 2016, and the number of light fog days increased by 125% from 4 days in 2012 to 9 days in 2016; in contrast, the numbers of dense and extreme fog days increased relatively slowly by 50% and 100%, respectively, from 4 to 6 days and 5 to 10 days. The slope corresponding to the 3-year moving average linear fitting in this stage far exceeded 0. The third stage was 2016–2019, during which the spatially averaged annual number of fog days remained constant with a slightly decreasing trend.
Corresponding to the annual variation of fog days, the spatially averaged annual specific humidity at 1000 hPa over Shandong Province also showed an increasing trend from 2008 to 2019 (Figs. 7c and d), while the spatially averaged annual wind speed showed a decreasing trend (Figs. 7e and f). Specific humidity exhibited a strong positive correlation with the number of fog days: the correlation coefficients with the numbers of total fog, light fog, dense fog, and extreme fog days were 0.591*, 0.594*, 0.484, and 0.579*, respectively (the symbol (*) indicates significance at the 5% significance level (two-tailed)), which indicated that the occurrence frequency of fog increased with increasing specific humidity. The specific humidity in 2013–2018 increased significantly compared with that before 2012; this abundance of water vapor during 2013–2018 enabled the increase in fog frequency after 2012. The numbers of total fog, light fog, dense fog, and extreme fog days exhibited negative correlations with the wind speed with coefficients of -0.487, -0.426*, -0.574, and − 0.483, respectively. In particular, the fog frequency increased rapidly with the abrupt decrease of wind speed after 2012, indicating that a lower wind speed was more conducive to the formation of fog. Moreover, the spatially averaged annual variation of haze days showed an increasing trend (Figs. 7g and h), which had a strong positive correlation with total fog, light fog, dense fog, and extreme fog days, and the correlation coefficients were 0.867**, 0.885**, 0.717**, and 0.826**, respectively (the symbol (**) indicates significance at the 1% significance level (two-tailed)). Therefore, the increase in water vapor, the aggravation in air pollution, and the decrease in wind speed may have been important factors governing the increased fog frequency in Shandong Province during 2008–2019.
Moreover, Figs. 6 and 7 indicate that the growth rate of light fog days was much larger than that of dense fog and extreme fog days, in which the growth rates were produced by the interaction of the specific humidity, wind speed, and air pollution. The meteorological conditions that are conducive to the formation of various types of fog are different (Tang et al. 2009; Wei et al. 2010). Therefore, the responses of variations of different levels of fog days to the changes of meteorological conditions are also different. This needs further detailed study.
3.3. Monthly and seasonal interannual and intra-annual variations in fog
Table 2 demonstrates that in the 12 years from 2008 to 2019, the spatially average numbers of total fog days in Shandong Province in spring (Mar., Apr., and May), summer (Jun., Jul., and Aug.), autumn (Sept., Oct., and Nov.), and winter (Dec., Jan., and Feb.) were 37, 35, 67, and 79, respectively (the spatially annual average in four seasons were 3.1, 2.9, 5.6, and 6.6 days, respectively). More fog days occurred in winter and autumn than in spring and summer. The standard deviations of fog days in summer and autumn were larger than those in spring and winter, and the degree of dispersion in the former two seasons was higher than that in the latter, suggesting a high and concentrated frequency of fog days at the meteorological stations in summer and autumn.
Table 2
Spatial average numbers and standard deviations of total fog, light fog, dense fog, and extreme fog in each season from 2008–2019.
Classification
|
Spring
(Std. Dev.)
|
Summer
(Std. Dev.)
|
Autumn
(Std. Dev.)
|
Winter
(Std. Dev.)
|
Light fog
|
11 (7)
|
13 (11)
|
24 (16)
|
27 (12)
|
Dense fog
|
10 (10)
|
10 (12)
|
16 (10)
|
18 (9)
|
Extreme fog
|
16 (20)
|
11 (30)
|
28 (20)
|
35 (16)
|
Total fog
|
37 (31)
|
35 (44)
|
67 (41)
|
79 (31)
|
The monthly and seasonal interannual variations in fog in Shandong Province during 2008–2019 were further analyzed. In general, the overall frequency of fog in Shandong Province in autumn and winter was higher than that in spring and summer during 2008–2019 (Fig. 8a). Specifically, the distribution of light fog in each season was almost consistent with the distribution of total fog (Fig. 8b). The seasonal frequency distributions of dense fog and extreme fog were slightly different from those of total fog and light fog but were also similarly higher in autumn and winter than in spring and summer (Figs. 8c and d). Except for summer, the frequency of dense fog and extreme fog increased in other seasons over the study period. The frequency increased the fastest in winter for different levels of fog. The frequency of light fog remained nearly unchanged in all seasons from 2008 to 2012 and began to increase rapidly after 2012, but the frequency fell rapidly in spring, summer, and autumn after 2017. In contrast, the frequency of light fog in winter increased rapidly after 2017 (Fig. 8b). Dense fog and extreme fog displayed similar trends. The difference was that the frequency of dense fog in autumn was the highest during 2015–2018, and in winter, it fell rapidly after 2014 (Figs. 8c and d). All levels of fog had high frequencies in Nov., Dec., and Jan., followed by Jul.-Oct. The frequency of fog increased nearly every month over the study period; among the three levels of fog, light fog was the main contributor to the increase in total fog frequency. The months with the highest frequency of total fog, light fog, dense fog, and extreme fog were Dec. 2015, Dec. 2015, Dec. 2015, and Jan. 2013, respectively (Figs. 8e-h).
In order to further analyze the reason for the higher fog frequency in autumn and winter, we calculated the spatial average lower-tropospheric stability (LTS). The LTS reflects the strength of the inversion, which is the difference in potential temperature at 700 hPa and the surface (LTS = θ700 - θ0, Klein and Hartmann 1993; Wood and Hartmann 2006; Zhang et al. 2010). In this way, the strength of the inversion of each month and season was captured.
Figure 9 shows the spatial average LTS of spring, summer, autumn, winter, and each month in Shandong Province from 2008 to 2019. It can be seen that LTS had a distinct seasonal cycle. The LTS in autumn and winter was greater than that in spring and summer, which indicated that during autumn and winter, the boundary layer tends to form a stronger low-level temperature inversion than that in spring and summer. It was obvious that the seasonal cycles of fog frequency and LTS were generally in phase, which illustrated that the fog frequency was high when LTS was large, while low when LTS was small (Figs. 8 and 9). Inversion is conducive to the formation of fog because it stabilizes the lower atmosphere (Liu et al. 2012; Liu et al. 2016). Therefore, the stronger inversion in autumn and winter formed a favorable low-level atmospheric stability condition for frequent fog. Moreover, the statistical results also show that the temporal-spatial average wind speeds at 1000 hPa in spring, summer, autumn, and winter in Shandong from 2008 to 2019 were 1.35 m/s, 1.27 m/s, 1.22 m/s, and 1.18 m/s, respectively, revealing that the wind speed in autumn and winter was smaller than that in spring and summer.
3.4. Diurnal variation and duration of fog
Statistical analyses were conducted on the frequency of total fog, light fog, dense fog, and extreme fog at each time point, the start time point of fog, and the end time point of fog over the day at the 122 meteorological stations in Shandong Province during 2016–2019. It was found that foggy weather in Shandong Province can occur at any time of the day, with fog frequently occurring from 00:00 to 08:00 (Fig. 10a) (Note: All the time in this paper is provided in China Standard Time, CST). The frequencies of light fog, dense fog, and extreme fog in this period accounted for 66.8%, 64.8%, and 70.8% in one day, respectively, all peaking at 06:00 (accounting for 9.1%, 9.7%, and 10.1%, respectively). The minimum frequencies of different levels of fog were reached at 15:00 (accounting for 0.9%, 0.9%, and 0.6%, respectively). Near this time, the temperature was typically very high (Fig. 10a), and the supersaturation required for fog formation was difficult to achieve, hindering the formation of fog. Figure 10a shows that the temperature gradually decreased from 15:00 to approximately 06:00 the next day. Low temperature at night creates the condition for the surface air to become saturated, allowing water vapor to condense and form fog (Belorid et al. 2015; Bari et al. 2016); therefore, the proportion of fog gradually increased during this period. After sunrise, the temperature gradually rose, and the fog began to dissipate slowly. From 08:00 to 14:00, the temperature increased sharply and then reached the maximum temperature (Fig. 10a), so most of the fog tended to dissipate during this period. The daily cycles of fog frequency and temperature were generally reversed, which indicated that the diurnal variation in temperature may be responsible for the diurnal variation in fog frequency.
In Shandong Province, the start time concentrated mainly from night to the early morning, peaking at 05:00 (Fig. 10b), and the end time also concentrated primarily from night to morning, peaking at 06:00–07:00 (Fig. 10c). The diurnal variations of the start and end times of the different levels of fog were similar, and the end time was pushed forward by 1–3 h compared to the start time, indicating that the duration of fog was mainly 1–3 h.
According to the statistics of the fog duration, the duration of most fog events was 1 h in Shandong Province, with the frequency of this period accounting for 50.9%. The proportions of fog duration of 1–3 h, 4–6 h, 7–9 h, 10–12 h, 13–24 h, and 25 h and above were 69.0%, 13.7%, 7.0%, 4.0%, 4.8%, and 1.5%, respectively. The fog duration in Shandong Province was mainly between 1–3 h, but fog days still occurred with durations exceeding 12 h or even 24 h (Fig. 11).