<p>Fog is a very complex phenomenon (Gultepe et al., 2007). In some areas it can contribute substantially to hydrological and chemical inputs and is therefore of high environmental relevance (Blas et al., 2010). Fog formation is affected by numerous factors, such as meteorology, air pollution, terrain (geomorphology), and land-use.</p><p>In our earlier studies we addressed the role of meteorology and air pollution on fog occurrence (H&#367;nov&#225; et al., 2018) and long-term trends in fog occurrence in Central Europe (H&#367;nov&#225; et al., 2020). This study builds on earlier model identification of year-to-year and seasonal components in fog occurrence and brings an analysis of the deformation of the above components due to the individual explanatory variables. The aim of this study was to indicate the geographical and environmental factors affecting the fog occurrence.</p><p>&#160;&#160;&#160;&#160;&#160;&#160; We have examined the data on fog occurrence from 56 meteorological stations of various types from Romania reflecting different environments and geographical areas. We used long-term records from the 1981&#8211;2017 period.&#160;</p><p>&#160;&#160;&#160;&#160;&#160;&#160; We considered both the individual explanatory variables and their interactions. With respect to geographical factors, we accounted for the altitude and landform. With respect to environmental factors,&#160;&#160; we accounted for proximity of large water bodies, and proximity of forests. Geographical data from Copernicus pan-European (e.g. CORINE land cover, high resolution layers) and local (e.g. Urban Atlas) projects were used. Elevation data from EU-DEM v1.1 were source for morphometric analysis (Copernicus, 2020).</p><p>&#160;&#160;&#160;&#160;&#160;&#160; &#160;We applied a generalized additive model, GAM (Wood, 2017; Hastie & Tibshirani, 1990) to address nonlinear trend shapes in a formalized and unified way. In particular, we employed penalized spline approach with cross-validated penalty coefficient estimation. To explore possible deformations of annual and seasonal components with various covariates of interest, we used (penalized) tensor product splines to model (two-way) interactions parsimoniously, Wood (2006).</p><p>&#160;&#160;&#160;&#160;&#160;&#160; The fog occurrence showed significant decrease over the period under review. In general the selected explanatory variables significantly affected the fog occurrence and their effect was non-linear. Our results indicated that, the geographical and environmental variables affected primarily the seasonal component of the model. Of the factors which were accounted for, it was mainly the altitude showing the clear effect on seasonal component deformation (H&#367;nov&#225; et al., in press).</p><p>&#160;&#160;&#160;&#160;&#160;&#160;</p><p>&#160;</p><p>References:</p><p>Blas, M, Polkowska, Z., Sobik, M., et al. (2010). Atmos. Res. 95, 455&#8211;469.</p><p>Copernicus Land Monitoring Service (2020). Accessed online at: https://land.copernicus.eu/.</p><p>Gultepe, I., Tardif, R., Michaelidis, S.C., Cermak, J., Bott, A. et al. (2007). Pure Appl Geophys, 164, 1121-1159.</p><p>Hastie, T.J., Tibshirani, R.J. (1990). Generalized Additive Models. Boca Raton, Chapman & Hall/CRC.</p><p>H&#367;nov&#225;, I., Brabec, M., Mal&#253;, M., Dumitrescu, A., Geleti&#269;, J. (in press) Sci. Total Environ. 144359.</p><p>H&#367;nov&#225;, I., Brabec, M., Mal&#253;, M., Valeri&#225;nov&#225;, A. (2018) Sci. Total Environ. 636, 1490&#8211;1499.</p><p>H&#367;nov&#225;, I., Brabec, M., Mal&#253;, M., Valeri&#225;nov&#225;, A. (2020) Sci. Total Environ. 711, 135018.</p><p>Wood, S.N. (2006) Low rank scale invariant tensor product smooths for generalized additive mixed models. Biometrics 62(4):1025-1036</p><p>Wood, S.N. (2017). Generalized Additive Models: An Introduction with R (2nd ed). Boca Raton, Chapman & Hall/CRC.</p><p>&#160;</p>