Synoptic-scale controls of fog and low-cloud variability in the Namib Desert
- 1Institute of Meteorology and Climate Research, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
- 2Institute of Photogrammetry and Remote Sensing, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
- 3IPSL LISA, CNRS, Université Paris-Est Créteil, Université Paris, Créteil, France
- 4IPSL LATMOS, CNRS, Sorbonne Université, Université Paris-Saclay, Paris, France
- 5Institute for Atmospheric and Climate Science, ETH Zurich, Zurich, Switzerland
- 6Norwegian Institute of Bioeconomy Research, Ås, Norway
- 7Department of Environmental Sciences, University of Basel, Basel, Switzerland
Correspondence: Hendrik Andersen (firstname.lastname@example.org)
Fog is a defining characteristic of the climate of the Namib Desert, and its water and nutrient input are important for local ecosystems. In part due to sparse observation data, the local mechanisms that lead to fog occurrence in the Namib are not yet fully understood, and to date, potential synoptic-scale controls have not been investigated. In this study, a recently established 14-year data set of satellite observations of fog and low clouds in the central Namib is analyzed in conjunction with reanalysis data in order to identify synoptic-scale patterns associated with fog and low-cloud variability in the central Namib during two seasons with different spatial fog occurrence patterns. It is found that during both seasons, mean sea level pressure and geopotential height at 500 hPa differ markedly between fog/low-cloud and clear days, with patterns indicating the presence of synoptic-scale disturbances on fog and low-cloud days. These regularly occurring disturbances increase the probability of fog and low-cloud occurrence in the central Namib in two main ways: (1) an anomalously dry free troposphere in the coastal region of the Namib leads to stronger longwave cooling of the marine boundary layer, increasing low-cloud cover, especially over the ocean where the anomaly is strongest; (2) local wind systems are modulated, leading to an onshore anomaly of marine boundary-layer air masses. This is consistent with air mass back trajectories and a principal component analysis of spatial wind patterns that point to advected marine boundary-layer air masses on fog and low-cloud days, whereas subsiding continental air masses dominate on clear days. Large-scale free-tropospheric moisture transport into southern Africa seems to be a key factor modulating the onshore advection of marine boundary-layer air masses during April, May, and June, as the associated increase in greenhouse gas warming and thus surface heating are observed to contribute to a continental heat low anomaly. A statistical model is trained to discriminate between fog/low-cloud and clear days based on information on large-scale dynamics. The model accurately predicts fog and low-cloud days, illustrating the importance of large-scale pressure modulation and advective processes. It can be concluded that regional fog in the Namib is predominantly of an advective nature and that fog and low-cloud cover is effectively maintained by increased cloud-top radiative cooling. Seasonally different manifestations of synoptic-scale disturbances act to modify its day-to-day variability and the balance of mechanisms leading to its formation and maintenance. The results are the basis for a new conceptual model of the synoptic-scale mechanisms that control fog and low-cloud variability in the Namib Desert and will guide future studies of coastal fog regimes.
In moist climates, fog is typically viewed as an atmospheric phenomenon that disturbs traffic systems and negatively affects physical and psychological health (e.g., Bendix et al., 2011). In the hyperarid Namib Desert, however, the water input of fog is key to the survival of many species (e.g., Seely et al., 1977; Seely, 1979; Ebner et al., 2011; Roth-Nebelsick et al., 2012; Warren-Rhodes et al., 2013; Henschel et al., 2018; Gottlieb et al., 2019). Despite this ecological significance, the local mechanisms that control the formation and spatiotemporal patterns of regional fog in the Namib are not yet fully understood, and potential linkages to synoptic-scale variability have yet to be explored. With regional climate simulations suggesting a warmer and even dryer climate (James and Washington, 2013; Maúre et al., 2018), fog could become an even more essential water source for regional ecosystems in the future. However, the lack of understanding concerning fog and low-cloud (FLC) processes and their interactions with dynamics, thermodynamics, aerosols, and radiation in this region (Zuidema et al., 2016; Formenti et al., 2019) limits the accuracy of and confidence in projected changes in fog patterns (e.g., Haensler et al., 2011).
Field observations of local meteorological parameters and fog have led to the distinction between two main fog types occurring in the region: advection fog and high fog. Advection fog can form when a moist warm air mass is transported over a cool ocean (Gultepe et al., 2007) and has been reported to occur mainly during austral winter, affecting a coastal strip of <30–40 km (Seely and Henschel, 1998). High fog is described as a low stratus that frequently reaches more than 60 km inland between September and March and leads to fog where the advected stratus base intercepts the terrain (Seely and Henschel, 1998). While the two fog types are reported to be transported inland with different wind systems (for a review of local wind systems see Lindesay and Tyson, 1990), they are both described to be of an advective nature. In Olivier and Stockton (1989), a coastal low is described as the mechanism that, in case of a narrow coastal upwelling region, drives the onshore advection of foggy air masses into the region of Lüderitz in southern Namibia during austral summer, while during winter they find fog to be associated with cold fronts. However, they assume that, while undetected, coastal lows were also present in these cases, as they typically precede the passage of a cold front (Olivier and Stockton, 1989; Reason and Jury, 1990). Recent analyses of diurnal FLC characteristics have shown that the timing of FLC occurrence depends on the distance to the coastline, with FLCs occurring significantly earlier at the coast than further inland, which is an indication of the dominance of advective processes (Andersen and Cermak, 2018; Andersen et al., 2019). Also, measurements of fog microphysics during the AEROCLO-sA field campaign in the Namib suggest that the observed fog events were advected cloudy air masses from the ocean (Formenti et al., 2019). While it has long been acknowledged that other fog types (e.g., radiation fog and frontal fog) can occur in the Namib as well (e.g., Jackson, 1941; Nagel, 1959), many statements regarding fog formation mechanisms in the historical literature do not seem to be founded on extensive and coherent observational evidence. Until recently, the occurrence of radiation fog, i.e., fog formation near the surface due to local radiative cooling under clear-sky conditions and without advective influence (Gultepe et al., 2007), was seen as a comparably rare situation (e.g., Seely and Henschel, 1998; Eckardt et al., 2013). This was questioned when, based on analyses of stable isotopes of fog water samples, Kaseke et al. (2017, 2018) found that the majority of their collected fog water samples stemmed from sweet water sources and interpreted this as evidence of predominant occurrence of radiation fog. Based on these findings, they postulated a potential shift from advection-dominated fog to radiation-dominated fog in the Namib Desert (Kaseke et al., 2017). Thus, the importance of the various fog formation mechanisms is currently a subject of scientific debate.
The goals of this study are to better understand the synoptic-scale conditions under which regional FLCs in the Namib occur, to understand how synoptic-scale variability changes local conditions, and thereby to assess the relevance of different potential fog formation mechanisms. To these ends, a 14-year time series of geostationary satellite observations of FLCs in the central Namib is combined with reanalysis data and air-mass back trajectories to systematically analyze the large-scale dynamic conditions and air-mass characteristics that are associated with FLC occurrence in the Namib. The guiding hypothesis for this study is given in the paragraph below.
Fog and low clouds in the central Namib are primarily of an advective nature and therefore associated with distinct synoptic-scale patterns of atmospheric dynamics and air-mass history. Thus, they can be statistically predicted with information on atmospheric circulation.
2.1 Satellite observations of FLCs
The Spinning Enhanced Visible and Infrared Imager (SEVIRI) sensor, mounted on the geostationary Meteosat Second Generation (MSG) satellites, is ideally suited to provide spatiotemporally coherent observations of clouds. It features a spatial resolution of 3 km at nadir and scans its full disk every 15 min (96 hemispheric scans per day; Schmetz et al., 2002). In the context of this study, 14 years (2004–2017) of SEVIRI data are used to continuously detect FLCs with the algorithm developed by Andersen and Cermak (2018). The algorithm relies mostly on a channel difference in the thermal infrared (12.0–8.7 µm), and in an extensive validation against surface observations this technique has shown a good performance (97 % overall correctness of the classification). The 14-year FLC data set used here has already been applied to study spatial and temporal patterns of FLC occurrence along the southwestern African coast in Andersen et al. (2019). It should be noted that this satellite technique does not discriminate between fog and other low clouds.
The focus of this study is on FLCs in the central Namib, from which the majority of historical and present-day station measurements stem (e.g., Nagel, 1959; Nieman et al., 1978; Lancaster et al., 1984; Seely and Henschel, 1998; Kaseke et al., 2017; Spirig et al., 2019). To provide a representative measure of the overall FLC cover in the central Namib on a daily basis, FLC occurrence is averaged between 03:00 and 09:00 UTC (local time is UTC +2 h) in the region between 22 and 24∘ S and up to 100 km inland. Only pixels with at least a 5 % FLC occurrence frequency in the climatology (as in Andersen et al., 2019) are used. A specified averaging time period is needed to avoid statistically mixing two separate FLC events occurring on successive nights, which would be the case in a daily average FLC occurrence data set. The specific time period is chosen to include all periods of the diurnal cycle, with FLC occurrence rising, peaking, and starting to dissipate (Andersen and Cermak, 2018) during this time. The spatial and temporal averaging is illustrated for an exemplary day in Fig. 1a. While the specific day shown here is arbitrary, the general feature of maximum FLC cover in the early morning hours and rapid decline shortly after sunrise is typical of the region (Andersen and Cermak, 2018). For further analyses, the data set is divided into “FLC days” with mean regional FLC cover exceeding 50 % between 03:00 and 09:00 UTC and “clear days” with mean FLC cover below 3 %. These thresholds are chosen to represent two clearly separated parts of the FLC cover distribution that occur with similar frequencies. The resulting distribution of daily average FLC cover and the number of cases in each class are shown in Fig. 1b. As the time of sunrise varies by season, the constructed data set is likely to feature a seasonal bias in FLC occurrence. It should be noted that this has no effect on the separation of FLC days and clear days within seasons, the analysis of which is the main purpose of this data set. The resulting monthly average FLC cover in the central Namib (Fig. 1c) should not be used in a quantitative sense but rather to illustrate the general seasonal cycle of FLCs in this region. It is interesting to note that the seasonal cycle of FLCs is not necessarily coupled to the seasonal cycle of fog occurrence due to the seasonal cycle in the vertical position of the low-cloud layer. For example, at coastal locations fog peaks in the central Namib between April and August (Andersen et al., 2019), while marine fog over the adjacent Atlantic has been found to peak between March and May, with a minimum occurrence between June and August (Dorman et al., 2020).
2.2 ERA5 reanalysis
To investigate the large-scale meteorological conditions associated with FLCs in the central Namib, ERA5 reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF) are used. ERA5 is the new generation of reanalysis and follow-up of ERA-Interim (Dee et al., 2011). In comparison to ERA-Interim, it features higher spatial (0.25∘) and temporal (hourly) resolutions, along with other improvements (Hersbach, 2016).
In the context of this study, 14 years (2004–2017) of meteorological fields are analyzed. To characterize large-scale dynamic and thermodynamic conditions, fields of mean sea level pressure (MSLP), geopotential height at 500, 700, 850, and 925 hPa (Z500, Z700, Z850, and Z925), 2 m air temperature (T2m), sea surface temperature (SST), total columnar water vapor (TCWV), specific humidity (Q), as well as winds at 10 m and at all ERA5 pressure levels between 1000 and 500 hPa, and lower tropospheric stability (LTS; computed as the difference between potential temperature at 700 hPa and T2m, Klein and Hartmann, 1993) are used. To represent the morning conditions for which FLC is averaged, 06:00 UTC fields of ERA5 data are selected. While, for additional analysis, T2m fields are also used at nighttime (01:00 and 03:00 UTC), the 06:00 UTC fields are used if no specific information on time is given.
2.3 Trajectory analysis
The 24 h backward trajectories are calculated using the Lagrangian analysis tool (LAGRANTO; Sprenger and Wernli, 2015). The three wind components needed for the trajectory calculations are taken from ERA5 on a regular 0.5∘ latitude–longitude grid, at 137 model levels in the vertical, and at a 3-hourly temporal resolution. The spatial resolution is used to reduce the data volume and computational cost. While the native resolution would be preferable, the general patterns of the trajectories are not expected to be affected, as tests with lower-resolution ERA-Interim data showed comparable results. The trajectories are started daily at 06:00 UTC for the periods April, May, June, September, October, and November 2004–2017. Their starting points are located in the central Namib, close to Gobabeb at 23∘ S, 15∘ E, 25 hPa above the surface (at ≈940 hPa), which corresponds roughly to 200 m above ground level. By doing so, the back trajectories represent air masses for the levels where fog and low clouds in the region are typically observed (Andersen et al., 2019). In order to obtain insights about the physical properties of the air masses, the temperature, potential temperature, specific humidity, and relative humidity are tracked along the trajectories. The location is chosen, as it is the main site of both historic and present-day scientific activity in the region (Lancaster et al., 1984; Seely and Henschel, 1998; Kaseke et al., 2017; Spirig et al., 2019).
2.4 Principal component analysis
The atmospheric variability in the South Atlantic Ocean and southern African region is characterized by means of a principal component analysis (PCA; see Storch and Zwiers, 1999). PCA solves the eigenvalues of the data covariance matrix and projects data variability onto an orthogonal basis, i.e., decomposes data variability into independent variability modes. Each mode explains a fraction of the total variance and is represented by a spatial anomaly pattern and a standardized time series (namely, the principal components, PCs) accounting for the amplitude of the anomaly pattern. Here, the PCA is used to analyze daily fields of the zonal and meridional components of 10 m wind at 06:00 UTC in a domain centered on the Namib (0∘–40∘ E, 40∘ S–0∘ N). In the context of this study, the main modes of the wind variability are used to understand possible linkages between atmospheric circulation at the synoptic scale and the daily occurrence of FLCs in the Namib region. The wind fields are first remapped onto a 1∘ regular grid (PCAs are computationally expensive; Pham-Thanh et al., 2019). Then daily 06:00 UTC anomalies are computed by subtracting the 14-year climatological average wind components at each grid point. The PCA is applied to the covariance matrix of both components in the domain. Remapping to 1∘ resolution allows not only an accurate description of the atmospheric variability at synoptic scale but also a smoothing out of the variability associated with small-scale effects. The sensitivity of the PCA to the spatial resolution is tested by conducting the analysis based on wind fields remapped to a 2∘ resolution. The results of the two PCAs at the different resolutions are very similar, demonstrating their robustness. Daily anomalies are computed with respect to the 14-year sampling of the FLC data set, in order to compare wind and FLC variability over a homogeneous climatology.
2.5 Statistical prediction of FLCs
Statistical modeling of fog or low clouds is typically done by using local fields of a set of predictors, i.e., relevant meteorological fields and aerosol properties (e.g., Andersen et al., 2017; Adebiyi and Zuidema, 2018; Fuchs et al., 2018). The circulation-induced variability can be captured by spatial patterns of pressure fields (Deloncle et al., 2007; Yu and Kim, 2010; Sippel et al., 2019). A major challenge when using pressure fields (denoted X, as an n×p matrix of n samples and p predictors located on a grid) to predict a target variable is, however, that the number of (strongly correlated) predictors can quickly outgrow the number of observations. This typically leads to high-variance problems (overfitting) in classical statistical models. The issue can be overcome with shrinkage methods, as, e.g., regularized linear models (Hastie et al., 2001). These provide an extension of linear regression techniques that shrink the regression coefficients of a model by penalizing their size, thereby addressing the aforementioned high-variance issues (Hastie et al., 2001). Ridge regression is a specific example of a regularized linear model where the shrinkage is controlled by a value λ that shrinks the coefficients of the model towards zero using the L2 penalty (the squared magnitude of the coefficient value is added as a penalty term to the loss function). This method is well suited for cases with a large number of correlated predictors that are all relevant (coefficients >0) (Friedman et al., 2010). The method can be used for classification and regression (Friedman, 2012).
Here, the statistical learning method is used in a classification setting. That is, a binary response variable (FLC day or clear day) is modeled using logistic regression regularized with the ridge penalty. In logistic regression with a binary response variable, the “odds ratio”, , is estimated as a linear function of the predictors for any given day as follows:
with β0 the intercept and βT the model coefficients. From the odds ratio, the estimated probabilities and the corresponding class (FLC day or clear day) are determined for each sample. The ridge regression penalty based on the L2 norm, i.e., , is then incorporated as a constraint on the size of the regression coefficients in the objective function that is minimized to fit the model. The tuning parameter λ directly trades off between a more flexible regression model (small penalty, i.e., low λ value) that, however, possibly suffers from high-variance issues and a less flexible regression model. Accordingly, a larger value of λ enforces smaller (but nonzero) regression coefficients, and a smoother spatial map of regression coefficients is obtained as a result. The optimal λ value is derived through 10-fold cross-validation. For a more complete description of regularized (logistic) regression, the reader is referred to Hastie et al. (2001) and the ElasticNet vignette for a hands-on tutorial (https://web.stanford.edu/~hastie/glmnet/glmnet_alpha.html, last access: 5 October 2019). Model estimation and cross-validation was performed using the scikit-learn package in Python (Pedregosa et al., 2011).
The ridge regression method is used to predict FLC and clear days over the complete 14-year time series, using spatial patterns of 06:00 UTC (representative of averaging time of FLC cover; see Sect. 2.1) ERA5 MSLP fields in a large spatial domain centered on the central Namib (0–45∘ S, 8∘ W–38∘ E; shown in Fig. 2). The ERA5 pressure fields feature a spatial resolution of and as such lead to 33 485 predictor fields.
3.1 Dynamic and thermodynamic conditions
Figure 2 shows a climatology of the dynamic and thermodynamic characteristics of the southeastern Atlantic and southern African region based on 14 years (2004–2017) of ERA5 data. Two seasons are shown in the figure that are representative of two different fog regimes (described in the next paragraph; Seely and Henschel, 1998; Andersen et al., 2019): September, October, and November (SON) in Fig. 2a and c and April, May, and June (AMJ) in Fig. 2b and d. At this spatial scale, the South Atlantic High and the continental high control the characteristic near-surface flow patterns during both seasons. During SON, the South Atlantic High is more prominent and, in combination with the thermal contrast between land and ocean, results in the formation of a low-level jet during this time (Nicholson, 2010). This alongshore coastal jet intensifies the upwelling of cold water, which feeds back to amplify the jet by increasing the thermal land–ocean contrast (Nicholson, 2010). On a local scale, the near-coastal winds that drive the upwelling are additionally modulated by the coastal topography (Koračin et al., 2004). While more prominent during SON, when the more pronounced South Atlantic High produces stronger winds, coastal upwelling water of the Benguela Current is apparent in the relatively low SSTs along the southwestern African coastline during both seasons (Fig. 2c and d) and throughout the year (Nelson and Hutchings, 1983). During AMJ, continental high-pressure situations are the most prominent circulation pattern (Tyson et al., 1996; Garstang et al., 1996). This is visible in the more pronounced continental high-pressure system and leads to a marked amplification of the easterly flow over the southern African continent. In the Namib Desert, thermally and topographically induced local wind systems within the boundary layer modulate these synoptic air-flow patterns, and the significance of the induced diurnal oscillations can exceed that of the synoptic scale (Goldreich and Tyson, 1988; Lindesay and Tyson, 1990). The combination of large-scale subsidence and low SSTs along the coastline produces high-LTS conditions in the coastal marine regions adjacent to the Namib, specifically during SON (LTS contours in Fig. 2c and d). In the adjacent marine regions downwind of the central Namibian coast, these stable conditions promote the formation of the southeastern Atlantic stratocumulus cloud deck and controls its seasonal cycle (Klein and Hartmann, 1993; Andersen et al., 2017), where the SST component is responsible for most of the LTS seasonality. One should note that MSLP and LTS are both affected by the high elevation of the central plateau in southern Africa (see Fig. 6) and are not likely to be a perfect representation of near-surface pressure conditions and lower-tropospheric stability in this region. However, due to the joint consideration of regions in southern Africa with high topography, the low-lying central Namib, and marine regions, no one specific pressure level of geopotential height can adequately summarize near-surface conditions throughout this large domain. Additional analyses show that patterns obtained from MSLP fields in southern Africa are similar to those at 925 and 850 hPa (not shown).
As outlined in the introduction, distinct seasonal fog and FLC patterns have been identified in the central Namib (Lancaster et al., 1984; Seely and Henschel, 1998; Cermak, 2012; Andersen et al., 2019). During SON, described as “high FLC season” in Andersen et al. (2019), FLCs frequently occur in the central Namib as a low stratus or high fog (cloud base height on average ≈400 m a.s.l. – above sea level) that touches the ground inland, whereas during the “low FLC season” in AMJ, FLCs occur less frequently, do not extend as far inland, and are typically lower, at ≈200 m a.s.l., and thus register as fog (termed “advection fog” in Seely and Henschel, 1998) at locations closer to the coastline (Andersen et al., 2019). While the FLC occurrence in the central Namib peaks in austral summer and is lowest during winter, fog peaks at coastal locations in AMJ and at inland locations during SON (Seely and Henschel, 1998; Andersen et al., 2019) due to the seasonal cycle in the vertical position of the cloud layer (Andersen et al., 2019). For these reasons, this study focuses on mechanisms determining FLC variability within these two characteristic fog seasons.
It has been assumed that the occurrence of FLCs in the Namib region and their variability on diurnal to seasonal scales is driven by the position and strength of large-scale pressure systems, as this would affect occurrence and advection of low-level clouds, atmospheric stability, and SSTs (Lancaster et al., 1984; Cermak, 2012; Andersen and Cermak, 2018; Andersen et al., 2019). Coastal upwelling, which has been shown to determine marine sea fog patterns along the Namibian coastline (Dorman et al., 2020), in combination with the presence of a coastal low that drives the onshore advection of foggy air masses has been found to be a major driver of fog occurrence in southern Namibia during austral summer (Olivier and Stockton, 1989). One should note though that the relationship between SSTs and fog in the Namib region is complex; Olivier and Stockton (1989) point out that too large an upwelling extent can also lead to less fog in southern Namibia. Based on these insights and also on knowledge from related coastal upwelling systems (Cereceda et al., 2008; Johnstone and Dawson, 2010; Del Río et al., 2018; Dorman et al., 2020), it is clear that the Atlantic anticyclone, the SSTs, and the large-scale subsidence are main drivers of this coastal FLC system. While all of these links play a role in FLCs in the Namib, the influence of synoptic-scale variability has not been explored, and a more in-depth analysis is needed to estimate the importance of the different mechanisms for the day-to-day variability in FLCs in the Namib region.
3.2 Differences in meteorological conditions on FLC days and clear days
Figure 3 shows large-scale patterns of averaged monthly mean differences in (a) MSLP and 10 m winds, (b) Z500 and winds at the same pressure level, (c) T2m, (d) LTS, (e) SST, and (f) TCWV on FLC versus clear days (as defined in Sect. 2.1) in the central Namib (marked with a star) during the investigated 14-year period (all months are considered here). The average of monthly mean differences is chosen rather than the overall mean differences to account for the distinct seasonal cycle of FLC occurrence in the Namib (Fig. 1c). In each pixel, an independent two-sided t test is computed to identify significant differences between the two classes (contours show p values <0.01). It is apparent that the dynamical conditions (Fig. 3a and b) on FLC days differ significantly on the synoptic scale. On FLC days, MSLP over continental southern Africa is systematically lower by about 3–5 hPa. This anomaly of lower MSLP extends over the southeastern Atlantic ocean at about 30∘ S. In a smaller oceanic region along the coastline north of 23∘ S, MSLP is significantly higher, leading to an overall anomalously high land–sea pressure gradient in this region and an onshore flow anomaly of near-surface winds in the central Namib on FLC days. The land–sea contrast in MSLP indicates a heat low over land, where the heat anomaly (Fig. 3c) could be driven by northerly advection ahead of the trough or enhanced surface warming. As discussed in Sect. 3.1, MSLP and 10 m winds may not be a good representation of near-surface level characteristics where topography is high; however, additional analyses of geopotential height at 850 and 925 hPa corroborate observed MSLP patterns. Differences exist in winds north of the central Namib, where at 925 and 850 hPa (not shown), a stronger onshore flow anomaly is observed than at 10 m, possibly indicating a topographical blocking of the onshore flow below the inversion. Z500 on FLC days (Fig. 3b) is significantly lower over the southeastern Atlantic between 30 and 40∘ S. This pattern is an indication of upper-level waves disturbing the mean tropospheric circulation of the southeastern Atlantic and southern Africa (Tyson et al., 1996; Fuchs et al., 2017). In combination, MSLP and Z500 show a weakly baroclinic structure with the mid-level trough shifted to the west (see Fig. A1). While a coastal low, which has been described in Olivier and Stockton (1989) as a local feature that can determine onshore flow, may still be present on FLC days, the composite differences between FLC days and clear days do not provide a clear indication of an increase in its presence on FLC days on average. However, as Reason and Jury (1990) describe, the coastal low is frequently followed by a frontal passage, which is a synoptic-scale signal observed here (Fig. A1).
There is a coherent pattern of slightly lower SSTs (≈0.5 K; Fig. 3e) along the coastline on FLC days; however, the difference between SSTs on FLC and clear days is not significant at the 0.01 level (and also not at the 0.05 level). It is interesting to note that SSTs tend to be lower on FLC days, although the coast-parallel, near-surface wind that partly governs the upwelling is slightly weaker in these cases (Fig. 3a), potentially hinting at a time-lag response of SSTs. This is to be expected, as Ekman transport produces a steady-state situation only after a few pendulum days (Pond and Pickard, 2013), although an initial upwelling response can be expected earlier (Lentz, 1992). It appears likely that effects of SST patterns on FLC variability are most pronounced on longer timescales (i.e., seasonal to interannual) that feature higher SST variability (Hutchings et al., 2009; Goubanova et al., 2013; Tim et al., 2015), as also observed in the Chilean Atacama Desert (Del Río et al., 2018). Differences in TCWV on FLC and clear days are pronounced (Fig. 3f). A coherent region of a significantly dryer column stretches from the central Namib over the coastal Atlantic, where the anomaly is strongest. This is likely the dry slot (Browning, 1997) or dry air intrusion of the synoptic-scale disturbance, which leads to increased longwave cooling at cloud top in case of FLC presence and has been shown to be a main determinant of cooling within the marine boundary layer (Koračin et al., 2005). This enhanced cooling can increase FLC cover, which has been observed to be a significant mechanism for stratocumulus clouds over the southeastern Atlantic (Adebiyi et al., 2015; Adebiyi and Zuidema, 2018). A substantial moist anomaly is visible over the southern African continent, likely driven by large-scale free-tropospheric moisture transport from the northwest (Fig. 3b). These moist air masses may contribute to the observed T2m heat anomaly via greenhouse warming (Fig. 3c). This effect of free-tropospheric moisture on surface temperatures has been observed in the Kalahari Desert (Manatsa and Reason, 2017) and other arid or semi-arid regions before (Evan et al., 2015; Oueslati et al., 2017; Alamirew et al., 2018). Along the coastal strip that is typically overcast with FLCs (Olivier, 1995; Cermak, 2012; Andersen and Cermak, 2018; Andersen et al., 2019), T2m is significantly lower by about 4 K, which is likely a feedback of FLCs reflecting solar radiation and slowing down the surface heating in the early morning (Iacobellis and Cayan, 2013) or due to air-mass differences. The observed difference patterns in LTS (Fig. 3d) between FLC and clear days match those of T2m so that they can be assumed to be mostly driven by its surface component (Pearson correlation coefficient is −0.90 for land pixels).
The observed anomaly patterns indicate that different mechanisms are triggered by the observed synoptic-scale disturbances and may contribute to FLC occurrence in the central Namib in two main ways:
increased FLC cover due to increased longwave cooling under the dry anomaly close to the coast;
onshore flow anomaly of marine boundary-layer air masses due to (a) a modulation of coastal winds and (b) a formation of a southern African heat low due to greenhouse warming by moist air masses and northerly warm air advection.
As both synoptic and FLC characteristics differ substantially between the SON and AMJ, the following section focuses on specific characteristics and differences in these mechanisms during these seasons.
3.3 Seasonal differences in synoptic-scale mechanisms
Figures 4 and 5 show seasonally averaged differences between FLC and clear days of all analyzed parameters during the two seasons SON and AMJ. During both seasons, MSLP (Fig. 4a and b) and Z500 (Fig. 4c and d) indicate synoptic-scale disturbances on FLC days. However, seasonal differences exist, as the disturbance is more pronounced during AMJ. The negative continental MSLP anomalies on FLC days are larger during AMJ, likely amplified by the more pronounced T2m anomalies and subsequent effects on a continental heat low during this time (see Fig. 5c and d). As noted above, the continental heat anomaly can be caused by northerly warm air advection or enhanced warming due to changes in the radiative balance. The observed seasonal MSLP and 10 m wind anomalies (Fig. 4a and b) that result in a transport of warm air from the northeast into the anomaly region (see Fig. A2), as well as the TCWV anomalies (Fig. 5a and b), suggest that during SON, the heat anomaly on FLC days is mostly due to northerly advection of warm air, whereas during AMJ, TCWV is significantly increased over the southern African continent. Here, the T2m anomalies closely follow those of TCWV (Pearson correlation coefficient of 0.75 in continental regions, with significantly higher T2m on FLC days than clear days), suggesting that the increased moisture causes an additional surface heating due to greenhouse warming as discussed in Sect. 3.2. It is likely that the TCWV anomaly is caused by a large-scale free-tropospheric moisture transport from the tropics, which is supported by the marked wind anomalies at 500 hPa (Fig. 4d) that show a northwesterly anomaly and the absolute wind and moisture fields at 700 hPa during this time (Fig. A1). It should be noted that a Lagrangian transport of moisture at this scale takes time and as such is likely to occur when the disturbance is relatively stationary or if two consecutive systems pass within a short timeframe (Knippertz and Martin, 2005).
While the yearly averaged composites show that over land, LTS is driven to a large extent by T2m (Fig. 3c and d), this is not quite as pronounced during SON (correlation coefficient ; Fig. 5c and e). Over continental southern Africa, the differences in T2m (Fig. 5c) are frequently compensated for by similar differences in potential temperature at 700 hPa (not shown). The most pronounced LTS feature during both seasons, however, is the coastal anomaly of increased LTS (over land and weaker over the adjacent ocean), which is driven by T2m. As this anomaly is also apparent during nighttime (01:00 and 03:00 UTC, not shown), it is likely that this pattern is mainly due to the relatively warm subsiding continental outflow that is apparent on clear days, rather than a radiative effect of FLCs as found in California (Iacobellis and Cayan, 2013). During AMJ, LTS is significantly lower over a large marine region south of 25∘ S, which is likely caused by the synoptic-scale disturbance.
During both seasons, SSTs in the coastal upwelling region are slightly lower on FLC days than on clear days, although these differences are not significant at the 0.01 level for the most part (very localized regions at ≈28∘ S are significantly lower during AMJ). In isolated patches further south, upwind of the study area, SSTs tend to be significantly higher on FLC days. This could lead to increased surface latent heat fluxes, increasing the moisture content of the marine boundary layer, particularly during AMJ when stronger near-surface winds are also apparent. A few 100 km to the west and south of the Namibian coastline, SSTs could similarly add to the increased moisture within the marine boundary layer. It is not clear yet, however, what exactly drives the observed anomaly patterns of SSTs. As upwelling reacts to the time-integrated wind field forcing over longer timescales than analyzed here (Pond and Pickard, 2013), the SST response to the instantaneous winds that are considered here is expected to be relatively weak. However, in the case of a relatively stationary disturbance as discussed above, the upwelling patterns could indeed reflect an SST response to a synoptic forcing. While the seasonally varying TCWV and SST anomalies (Figs. 4e, f and 5a, b, respectively) illustrate the seasonal variability in the mechanisms that can contribute to FLC occurrence in the central Namib, during all months, the outlined systematic patterns of significant negative MSLP anomalies over continental southern Africa and the localized coastal high-pressure anomaly are apparent. It can be concluded that a low-pressure anomaly in continental southern Africa and the associated onshore advection of marine boundary-layer air masses facilitate FLC occurrence in the central Namib during the entire year.
To better understand the characteristics of the observed moisture transport and its relevance for FLC occurrence in the central Namib, information on the vertical patterns of moisture and wind anomalies is needed. Figure 6 shows average seasonal differences in Q and winds on FLC versus clear days at different pressure levels during (a) SON and (b) AMJ (averaged between 20 and 25∘ S). During both seasons, a complex vertical structure of Q anomalies is apparent, which is assumed to be disturbance-induced. During both seasons, the marine boundary layer features an onshore flow anomaly and is more humid on FLC than on clear days, especially during AMJ, when this is a synoptic-scale feature, likely related to the cold front of the disturbance. These differences are caused by the subsiding dry continental easterly air masses that dominate on clear days, whereas on FLC days, a slight onshore flow of the more humid marine boundary-layer air is observed in the central Namib. Over land, these marine air masses flow against the dominant continental easterly winds (Lindesay and Tyson, 1990), producing a northerly wind flow at N (not shown) that has been found to be associated with fog occurrence in the central Namib (Seely and Henschel, 1998; Spirig et al., 2019). Above the moist marine boundary layer, the free troposphere is relatively dry on FLC days during both seasons, a feature which is not as clearly visible in the columnar TCWV composites during AMJ, as it is masked by the moist anomaly in the marine boundary layer (Fig. 5a). It is interesting to note that the marine dry anomaly peaks between December and February (not shown), the season with maximum FLC cover in the central Namib, with TCWV anomalies exceeding 10 kg m−2. The seasonal difference in the free-tropospheric Q anomalies over the continent is clear, and the vertical distribution of Q anomalies during AMJ corroborates the assumption that the observed positive TCWV anomalies are due to free-tropospheric moisture transport (Fig. 5b). Expressed in relative terms, Q is approximately halved within the dry anomaly region on FLC days during both seasons, suggesting that radiative cooling is an important factor for FLC cover, especially over marine regions where the dry anomaly is most pronounced. During AMJ, the free-tropospheric relative moisture difference between FLC days and clear days is observed to be as high as 220 %. This substantial increase in free-tropospheric moisture in this otherwise dry central plateau region induces a substantial surface heating, contributing to the formation of the observed heat low, which modulates regional wind systems and leads to the onshore flow anomaly.
It should be noted that in a comparable upwelling system (coastal California), Clemesha et al. (2017) also find a positive relationship between T2m over land and coastal low-level cloudiness, with the T2m anomaly shifted poleward by about 5∘ latitude from the cloud field. They propose that the T2m–cloud relationship is due to spatially offset associations between coastal low-level cloudiness and stability (potential temperature at 700 hPa), which is strongly correlated to T2m over land, thereby resulting in the T2m anomaly, rather than T2m driving the onshore advection. While in the central Namib the anomaly patterns between potential temperature at 700 hPa and T2m are similar in that they are also positively correlated during SON (and therefore compensate each other in terms of LTS; Fig. 5c and e), they are uncorrelated during AMJ (and also in the annual averages), when T2m over land is strongly correlated to TCWV. Also, during all times of year, the T2m and MSLP anomalies are directly inland from the cloud field, suggesting an influence on onshore advection.
3.4 The role of air-mass history and dynamical regimes
Air-mass back trajectories, initiated in the central Namib, close to Gobabeb at 23∘ S, 15∘ E (indicated by the star in Fig. 7), 06:00 UTC, and 25 hPa above ground level (approximates 200 m above ground level), are computed for the 14-year observational period. Figure 7 shows the back trajectories for FLC days (a, b) and clear days (c, d) for the two seasons SON (a, c) and AMJ (b, d). During both seasons, air masses on FLC days nearly exclusively stem from the marine boundary layer and have traversed over the cool upwelling ocean water along the coastline for the time span of 24 h. This is in agreement with findings from Koračin et al. (2005), which noted that a marine origin of air masses is critical, and potential mixing with continental air masses along the trajectory that would lead to a warming and drying. While the number of FLC days during SON is higher than during AMJ, following the general seasonality of FLCs in the region (see Fig. 1c), no clear seasonal differences in air-mass dynamics can be observed in such situations. This suggests that during both seasons, similar local dynamic conditions drive FLCs or air masses that develop into FLCs inland into the Namib Desert but that due to seasonal differences in large-scale dynamics, these situations occur with varying frequency during different seasons.
On clear days, air-mass histories are more diverse and show distinct seasonal differences but are frequently characterized by subsiding continental air masses. While on clear days during SON a considerable fraction of the air masses is still transported from the marine boundary layer, during AMJ, subsiding northeasterly continental air masses dominate. This seasonal shift in air-mass dynamics is likely driven by the seasonality of the two dominating high-pressure systems of the region that is shown in Fig. 2. During AMJ, the continental high-pressure system is enhanced and leads to the stronger easterly flow. These observations support the hypothesis by Lancaster et al. (1984) that the seasonality of fog in the central Namib is to some extent controlled by the southern African high-pressure system, as the associated easterly winds are likely to inhibit large-scale onshore advection of cloudy marine boundary-layer air masses. The results also suggest that aerosols from the biomass burning season in continental southern Africa (Swap et al., 2003) are unlikely to play a large role for fog formation by acting as cloud condensation nuclei, as biomass burning aerosols within the boundary layer are mostly associated with continental air masses in this region (Formenti et al., 2018). However, biomass burning aerosols may influence FLCs in the Namib region by absorbing solar radiation and modifying the thermodynamic conditions, which has been modeled and observed to influence the Namibian stratocumulus deck (Zhou and Penner, 2017; Deaconu et al., 2019).
While systematic differences in air masses exist between FLC days and clear days, clear days may still feature air masses that are advected from the marine boundary layer (see Fig.7c). To understand the differences between the FLC days and clear days in such situations, these are isolated and analyzed in the following. Figure 8 shows the average Q, relative humidity (RH), air temperature (T), potential temperature (pot. T), and pressure (P) along all of the back trajectories that are advected from the marine boundary layer (here, P>900 hPa over ocean). It is apparent that these air masses contain significantly more moisture and feature significantly lower pot. T on FLC days than on clear days, which explains most of the difference in RH. The back trajectories of FLC days feature a stronger cooling during the last 10 h of advection (hours 0–10), resulting in an additional increase in RH. The deviation in T between FLC and clear days seems to be driven by the vertical movement of the air masses, rather than differences in radiative cooling, as no changes in pot. T are apparent. At 10 h before initialization, air masses on clear days are located ≈20 hPa higher than on FLC days, not cooling off as they are advected due to their simultaneous subsidence. Other potential factors that may drive the observed deviation in T, such as the free-tropospheric moisture content and the surface temperature along the back trajectories, were not found to be systematically different on FLC and clear days (not shown). These findings highlight that FLCs Namib in the region are dependent on dynamics and also that marine boundary-layer moisture content and temperature changes during advection are important controls as well.
It is likely that the computed air-mass back trajectories do not fully capture thermally and topographically induced local air-flow patterns (see Lindesay and Tyson, 1990, for a review) that contribute to local FLC occurrence patterns and possibly formation. However, the larger-scale patterns of air-mass history of marine boundary-layer air masses versus the subsiding continental air masses from the free troposphere are clearly evident from the analysis presented and offer a consistent physical explanation of the large-scale FLC occurrence patterns. The observations suggest that FLCs in the Namib region are either advected after forming over the cool adjacent ocean or condensation takes place during advection of the marine boundary-layer air masses over land due to higher humidity levels, lower temperatures or radiative cooling, though a mix of these processes is likely.
The analysis of air-mass back trajectories shows that the discrimination between FLC and clear days is not possible using dynamics alone and that seasonal differences exist in the link between the probability of FLC days and advection patterns. To further investigate the role of different dynamical regimes for FLC occurrence, a PCA is conducted on spatial patterns of synoptic-scale near-surface winds (see Sect. 2.4 for details on the method). Figure 9a shows correlations between daily FLC cover in the central Namib and the PCs associated with the first six modes of variability in near-surface winds during all months of the year. All PCs are significantly correlated to FLC cover during some months of the year. Clear correlation patterns are evident; PCs 1, 2, 4, and 5 show negative correlations with FLC cover, while PCs 3 and 6 feature positive correlations. These PCs that facilitate FLC occurrence (3 and 6) show westerly or northwesterly wind anomalies in the central Namib, while PCs that are negatively associated with FLC cover in the Namib region feature anomalously strong continental easterly winds, consistent with results presented in Sect. 3.2 and 3.3. Figure 9c and d show the spatial patterns of near-surface wind anomalies of PCs 3 (explained variance: 11 %) and 4 (explained variance: 6 %) as examples for PCs that promote and impede FLC occurrence, respectively. A seasonal dependence of the correlations on PCs and FLC cover is apparent and seems to be related to the seasonality of FLC cover (Fig. 9b); PCs associated with onshore circulation in the central Namib feature the strongest positive correlations during winter when FLC cover is generally lowest over the Namib, especially evident for PC 3. This appears plausible, as during winter, the typical dynamical setting is less conducive to FLCs (see Fig. 2 for fall/early winter conditions during AMJ), and consequently, FLC occurrence is dependent on a stronger dynamical disruption during this time. During summer, when FLCs frequently occur in the central Namib, dynamical conditions associated with PCs 4 and 5 (dominance of continental easterlies) seem to impede the occurrence of FLCs. The results underscore that the advection of marine air masses is crucial for the occurrence of FLCs in the central Namib.
3.5 Statistical fog and low-cloud prediction with pressure fields
Based on the evidence presented above, showing that FLC occurrence is tightly connected to synoptic-scale patterns, it can be assumed that FLC occurrence can be predicted to some extent with a statistical learning technique that utilizes spatial patterns of dynamical information. Here, a ridge regression is applied to classify FLC days and clear days based on MSLP fields in a region spanning that is centered on the central Namib (see Sect. 2.4). MSLP fields are used, as their anomaly patterns on FLC days are similar during the different analyzed seasons and thus summarize the controlling mechanisms of onshore advection of marine boundary-layer air masses. Figure 10a shows the resulting coefficients, i.e., regression slopes, of the statistical model. The sign and spatial patterns of the coefficients are similar to the observed MSLP anomalies shown in Fig. 4, where coefficients (and anomalies) are negative in the inland region of Namibia and positive along the northern part of the Namibian coastline. It should be noted that the statistical model seems to mostly rely on regional MSLP fields, resulting in low coefficients at the synoptic scale, e.g., the Atlantic high-pressure system. It can be concluded that the synoptic-scale pressure patterns set the stage for more localized pressure and wind modulations that determine FLC occurrence and that regional MSLP fields contain information on both.
Figure 10b gives a summary of statistical measures of the skill of the model to classify between FLC and clear days in the central Namib. Using MSLP fields at 06:00 UTC on the day of the FLC cover information, the ridge regression model has a probability to correctly detect FLC days of 94 %, with 17 % of the reported FLC days being false alarms, leading to an overall correctness of the model of 86 % and a positive bias of 14 %. The critical success index (CSI: 0.79) and the Heidke skill score (HSS: 0.72) combine these scores and show that the model is skillful in distinguishing between the defined fog and clear days. As MSLP fields in southern Africa may not be representative due to the high topography, the model was additionally run based on Z850 and Z925 fields. The model performances were nearly the same (overall PC of 84 % in both cases), suggesting that it is adequate to use MSLP in this context. The colored dots in Fig. 10b illustrate the progression of the model skill when the training is carried out based on MSLP fields of 1, 2, 3, or 4 d prior to the FLC observation. While, as expected, the model skill deteriorates with an increasing temporal gap between the MSLP predictors and the time of FLC occurrence, the model is capable of predicting fog occurrence fairly well 1 d in advance, as the time series of day-to-day FLC occurrence features a significant autocorrelation of some days. To some extent, this may be connected to the strong persistence of synoptic-scale dynamics in the subtropics. Even though the model only uses MSLP fields, ignoring, e.g., effects of radiative cooling due to moisture anomalies, surface temperatures, and seasonal characteristics, which have been shown to modify FLC occurrence, the results still illustrate the potential of a dynamics-based statistical fog forecast in this region. It should be noted that changes in circulation additionally influence upwelling intensity, (e.g., Hutchings et al., 2009) such that some of the explained variability may also be attributed to factors influencing FLC formation rather than advection. However, due to the longer time scale of SST responses, and due to the marked contrasting differences in air mass history on FLC and clear days, the latter is thought to be the first-order mechanism leading to the high model skill.
It should be noted that the distinction between FLC and clear days is based on spatially and temporally averaged FLC occurrence (see Sect. 2.1) and that days are omitted that feature an FLC cover between 3 % and 50 %. Also, the exact location and time of FLC occurrence is likely to be dependent on local temperature gradients and topography that lead to local modulation of winds (Lindesay and Tyson, 1990). Still, the model produces promising results that may be built upon in future studies by testing a similar model setup to predict the timing and duration of FLCs at specific locations.
In this study, the occurrence of FLCs in the Namib Desert, derived from 14 years of satellite observations, is systematically analyzed within the context of the regional climate and related to large-scale patterns of MSLP, Z500, T2m, LTS, SST, and TCWV, as well as Q and winds at different pressure levels from ERA5 reanalyses. The satellite data set of FLC occurrence is separated into FLC days and clear days, which are further investigated in terms of their meteorological conditions, air-mass histories, and statistical predictability during two seasons (AMJ and SON).
It is found that MSLP and Z500 patterns on FLC days are systematically and significantly different from clear days on synoptic scales. On FLC days, a systematic pattern of significantly lower MSLP over continental southern Africa is observed, which, in combination with higher pressure over a marine coastal region at about 20∘ S, leads to an onshore flow anomaly of marine boundary-layer air. Together with significantly lower Z500 in the southeastern Atlantic region on FLC days, these dynamic patterns are an indication of synoptical-scale disturbances. These modify circulation systems, which in turn alter moisture transport, resulting in characteristic moisture patterns on FLC and clear days. Over the coastal boundary layer, the free troposphere is observed to be significantly drier on FLC days during both seasons, increasing radiative cooling, which likely increases FLC coverage, especially over the ocean where the dry anomaly is observed to be most pronounced. During AMJ, free-tropospheric moisture over the southern African continent is substantially increased, leading to greenhouse warming at the surface. While northerly warm air advection also contributes to the observed positive T2m anomalies on FLC days (during both seasons), the additional increase in T2m on FLC days during AMJ clearly corresponds to regions of increased free-tropospheric moisture content (correlation = 0.75). The increase in T2m leads to the development of a heat low that amplifies the upper-level, disturbance-induced low-MSLP anomaly, thereby contributing to the onshore flow anomaly of marine boundary-layer air masses. In the localized coastal region where FLCs typically occur, T2m is found to be significantly lower on FLC days, likely a combination of a local feedback of FLCs that slow down surface heating in the morning hours and air mass differences. A significant pattern of SST anomalies is found only in AMJ, with anomalously high SSTs off the coast possibly acting together with increased near-surface winds to enhance surface latent heat fluxes, which may contribute to the observed higher levels of specific humidity in the marine boundary layer.
The analysis of back trajectories initialized in the central Namib at typical cloud level shows systematic differences in air-mass dynamics on FLC days and clear days. Air masses on FLC days are nearly exclusively transported within the marine boundary layer over the cool upwelling waters along the coastline, whereas clear days are frequently associated with subsiding northeasterly air masses, especially during AMJ. During SON, when advection of marine boundary-layer air masses can also occur on clear days, air masses on clear days feature significantly less moisture and tend to be advected from higher altitudes than on FLC days. The findings clearly demonstrate the strong dependence of FLC occurrence in the central Namib on the advection of moist marine boundary-layer air masses, contrasting the notion of predominant radiation fog (Kaseke et al., 2017), but in agreement with many other studies (e.g., Olivier and Stockton, 1989; Seely and Henschel, 1998; Formenti et al., 2018; Andersen et al., 2019; Spirig et al., 2019). These results are supported by a principal component analysis of near-surface winds that show a clear connection of FLC cover to synoptic-scale dynamics. Principal components of spatial wind patterns that feature positive onshore flow anomalies are positively related to FLC cover. This relationship is especially strong during winter, when FLC occurrence is at its minimum, as then the dominant continental easterly flow typically inhibits inland advection of FLCs or locally developing FLCs. This suggests that during this time, a stronger dynamical forcing is needed to overcome this characteristic flow that is unfavorable for inland advection of cloudy marine boundary-layer air masses.
As the results show that spatial pressure patterns are connected to FLC occurrence, a ridge regression model is used to classify FLC days versus clear days based on regional MSLP fields. The resulting spatial pattern of model coefficients is similar to the observed MSLP anomaly patterns within the region of Namibia and the adjacent ocean areas. The spatial domain of relevant model coefficients seems to be smaller than the spatial extent of the pressure anomalies, probably because the regional fields contain information on synoptic-scale disturbance as well as local modulation. On this basis, the model is capable of skillfully delineating FLC days from clear days. The model is trained with MSLP fields with different temporal offsets and found to be capable of skillfully predicting FLC occurrence 1 d in advance, highlighting the potential of a statistical forecast of FLCs in this region. Future work should focus, however, on the development of a statistical model that links information on, e.g., MSLP, free-tropospheric moisture, SSTs, Z500, and aerosol loading with FLC occurrence in order to quantify the effects of the different processes and mechanisms outlined in this study.
The findings of this study suggest that FLCs in the central Namib are facilitated by synoptic-scale disturbances in two main ways:
increased longwave cooling due to an anomalously dry free troposphere, especially over the ocean that increases low-cloud cover;
onshore flow anomaly of these cloudy marine boundary-layer air masses due to
disturbance-induced modulation of local winds, and
a heat low over continental southern Africa.
The magnitude and characteristics of the disturbance and the related mechanisms depend on season, with a more pronounced disturbance during AMJ, when the typical dynamic setting is less conducive to FLC occurrence. Figure 11 is a schematic illustration that summarizes these seasonally varying mechanisms.
While a 14-year sample is not optimal to capture climatological variability, the mechanisms documented here for the first time are unlikely to be fundamentally different in other climatological periods. While it seems settled that, at least at the scales considered in this study, FLC occurrence is mostly driven by advective processes, the quantitative contributions of humidity and temperature changes and radiative cooling for low-cloud formation in the Namib during the advection of marine boundary-layer air masses are still unclear. A heat budget analysis as in, e.g., Adler et al. (2019) or Babić et al. (2019), based on ground-based measurements conducted during the field campaign of the Namib Fog Life Cycle Analysis (NaFoLiCA) project (Spirig et al., 2019), is necessary to better understand the origin, development, and life cycle of FLCs within the advected marine boundary-layer air masses. Future work should also focus on understanding the local and possibly synoptic-scale drivers of the vertical structure of FLCs in the Namib region on diurnal to seasonal scales and the day-to-day variability in (marine) boundary-layer humidity. As FLCs in the Namib are clearly connected to marine stratus/stratocumulus clouds, findings of recent and ongoing field campaigns over the southeastern Atlantic (Zuidema et al., 2016; Formenti et al., 2019) and related insights concerning the aerosol–cloud–meteorology system of the Namibian stratocumulus cloud field (e.g., Adebiyi and Zuidema, 2018; Andersen and Cermak, 2015; Diamond et al., 2018; Formenti et al., 2019; Fuchs et al., 2017, 2018; Gordon et al., 2018) are relevant to fully understand FLCs in the Namib Desert.
with a as the number of hits, b as the number of false alarms, c as the number of misses, and d as the number of correct negatives.
ERA5 data were generated using Copernicus Climate Change Service Information (2019). Satellite data and code for data processing are available from the corresponding author upon reasonable request.
HA and JC had the idea for the analysis. HA obtained and analyzed most of the data sets, conducted the original research, and wrote the article. JC and JF contributed to the study design, and JF computed initial back trajectories. PK helped to develop a conceptual understanding of the synoptic-scale patterns and physical mechanisms. JQ computed the back trajectories with LAGRANTO. MG conducted the PCA analysis. SS contributed to the design of the statistical model, and RV contributed insights to local-scale processes. JC, JF, PK, JQ, MG, SS, and RV contributed to article preparation and the interpretation of findings.
The authors declare that they have no conflict of interest.
This article is part of the special issue “New observations and related modelling studies of the aerosol–cloud–climate system in the Southeast Atlantic and southern Africa regions (ACP/AMT inter-journal SI)”. It is not associated with a conference.
Funding for this study was provided by Deutsche Forschungsgemeinschaft (DFG) in the project Namib Fog Life Cycle Analysis (NaFoLiCA), CE 163/7-1. HA acknowledges receiving a Research Travel Grant from the Karlsruhe House of Young Scientists that supported a stay at ETH Zürich and thus facilitated the collaboration with Sebastian Sippel. The contribution of Julian Quinting was funded by the Helmholtz Association (grant VH-NG-1243). Marco Gaetani was supported by the French National Research Agency under grant agreement no. ANR-15-CE01-0014-01 (AEROCLO-sA). We acknowledge support by the KIT-Publication Fund of the Karlsruhe Institute of Technology. We thank the three anonymous reviewers for their valuable and constructive comments.
This research has been supported by the Deutsche Forschungsgemeinschaft (grant no. CE 163/7-1). The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.
This paper was edited by Paquita Zuidema and reviewed by Ian Faloona and two anonymous referees.
Adebiyi, A. A. and Zuidema, P.: Low cloud cover sensitivity to biomass-burning aerosols and meteorology over the southeast Atlantic, J. Climate, 2, 4329–4346, https://doi.org/10.1175/JCLI-D-17-0406.1, 2018. a, b, c
Adebiyi, A. A., Zuidema, P., and Abel, S. J.: The Convolution of Dynamics and Moisture with the Presence of Shortwave Absorbing Aerosols over the Southeast Atlantic, J. Climate, 28, 1997–2024, https://doi.org/10.1175/JCLI-D-14-00352.1, 2015. a
Adler, B., Babić, K., Kalthoff, N., Lohou, F., Lothon, M., Dione, C., Pedruzo-Bagazgoitia, X., and Andersen, H.: Nocturnal low-level clouds in the atmospheric boundary layer over southern West Africa: an observation-based analysis of conditions and processes, Atmos. Chem. Phys., 19, 663–681, https://doi.org/10.5194/acp-19-663-2019, 2019. a
Alamirew, N. K., Todd, M. C., Ryder, C. L., Marsham, J. H., and Wang, Y.: The early summertime Saharan heat low: sensitivity of the radiation budget and atmospheric heating to water vapour and dust aerosol, Atmos. Chem. Phys., 18, 1241–1262, https://doi.org/10.5194/acp-18-1241-2018, 2018. a
Andersen, H. and Cermak, J.: How thermodynamic environments control stratocumulus microphysics and interactions with aerosols, Environ. Res. Lett., 10, 24004, https://doi.org/10.1088/1748-9326/10/2/024004, 2015. a
Andersen, H. and Cermak, J.: First fully diurnal fog and low cloud satellite detection reveals life cycle in the Namib, Atmos. Meas. Tech., 11, 5461–5470, https://doi.org/10.5194/amt-11-5461-2018, 2018. a, b, c, d, e, f
Andersen, H., Cermak, J., Fuchs, J., Knutti, R., and Lohmann, U.: Understanding the drivers of marine liquid-water cloud occurrence and properties with global observations using neural networks, Atmos. Chem. Phys., 17, 9535–9546, https://doi.org/10.5194/acp-17-9535-2017, 2017. a, b
Andersen, H., Cermak, J., Solodovnik, I., Lelli, L., and Vogt, R.: Spatiotemporal dynamics of fog and low clouds in the Namib unveiled with ground- and space-based observations, Atmos. Chem. Phys., 19, 4383–4392, https://doi.org/10.5194/acp-19-4383-2019, 2019. a, b, c, d, e, f, g, h, i, j, k, l, m, n
Babić, K., Adler, B., Kalthoff, N., Andersen, H., Dione, C., Lohou, F., Lothon, M., and Pedruzo-Bagazgoitia, X.: The observed diurnal cycle of low-level stratus clouds over southern West Africa: a case study, Atmos. Chem. Phys., 19, 1281–1299, https://doi.org/10.5194/acp-19-1281-2019, 2019. a
Bendix, J., Eugster, W., and Klemm, O.: FOG – BOON OR BANE?, Erdkunde, 65, 229–232, 2011. a
Cereceda, P., Larrain, H., Osses, P., Farías, M., and Egaña, I.: The spatial and temporal variability of fog and its relation to fog oases in the Atacama Desert, Chile, Atmos. Res., 87, 312–323, https://doi.org/10.1016/j.atmosres.2007.11.012, 2008. a
Cermak, J.: Low clouds and fog along the South-Western African coast – Satellite-based retrieval and spatial patterns, Atmos. Res., 116, 15–21, https://doi.org/10.1016/j.atmosres.2011.02.012, 2012. a, b, c
Clemesha, R. E., Gershunov, A., Iacobellis, S. F., and Cayan, D. R.: Daily variability of California coastal low cloudiness: A balancing act between stability and subsidence, Geophys. Res. Lett., 44, 3330–3338, https://doi.org/10.1002/2017GL073075, 2017. a
Deaconu, L. T., Ferlay, N., Waquet, F., Peers, F., Thieuleux, F., and Goloub, P.: Satellite inference of water vapour and above-cloud aerosol combined effect on radiative budget and cloud-top processes in the southeastern Atlantic Ocean, Atmos. Chem. Phys., 19, 11613–11634, https://doi.org/10.5194/acp-19-11613-2019, 2019. a
Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011. a
Del Río, C., Garcia, J. L., Osses, P., Zanetta, N., Lambert, F., Rivera, D., Siegmund, A., Wolf, N., Cereceda, P., Larraín, H., and Lobos, F.: ENSO influence on coastal fog-water yield in the Atacama Desert, Chile, Aerosol Air Qual. Res., 18, 127–144, https://doi.org/10.4209/aaqr.2017.01.0022, 2018. a, b
Diamond, M. S., Dobracki, A., Freitag, S., Small Griswold, J. D., Heikkila, A., Howell, S. G., Kacarab, M. E., Podolske, J. R., Saide, P. E., and Wood, R.: Time-dependent entrainment of smoke presents an observational challenge for assessing aerosol–cloud interactions over the southeast Atlantic Ocean, Atmos. Chem. Phys., 18, 14623–14636, https://doi.org/10.5194/acp-18-14623-2018, 2018. a
Dorman, C. E., Mejia, J., Koračin, D., and McEvoy, D.: World marine fog analysis based on 58-years of ship observations, Int. J. Climatol., 40, 145–168, https://doi.org/10.1002/joc.6200, 2020. a, b, c
Ebner, M., Miranda, T., and Roth-Nebelsick, A.: Efficient fog harvesting by Stipagrostis sabulicola (Namib dune bushman grass), J. Arid Environ., 75, 524–531, https://doi.org/10.1016/j.jaridenv.2011.01.004, 2011. a
Eckardt, F., Soderberg, K., Coop, L., Muller, A., Vickery, K., Grandin, R., Jack, C., Kapalanga, T., and Henschel, J.: The nature of moisture at Gobabeb, in the central Namib Desert, J. Arid Environ., 93, 7–19, https://doi.org/10.1016/j.jaridenv.2012.01.011, 2013. a
Evan, A. T., Flamant, C., Lavaysse, C., Kocha, C., and Saci, A.: Water vapor-forced greenhouse warming over the Sahara desert and the recent recovery from the Sahelian drought, J. Climate, 28, 108–123, https://doi.org/10.1175/JCLI-D-14-00039.1, 2015. a
Formenti, P., Piketh, S. J., Namwoonde, A., Klopper, D., Burger, R., Cazaunau, M., Feron, A., Gaimoz, C., Broccardo, S., Walton, N., Desboeufs, K., Siour, G., Hanghome, M., Mafwila, S., Omoregie, E., Junkermann, W., and Maenhaut, W.: Three years of measurements of light-absorbing aerosols over coastal Namibia: seasonality, origin, and transport, Atmos. Chem. Phys., 18, 17003–17016, https://doi.org/10.5194/acp-18-17003-2018, 2018. a, b
Formenti, P., D'Anna, B., Flamant, C., Mallet, M., Piketh, S. J., Schepanski, K., Waquet, F., Auriol, F., Brogniez, G., Burnet, F., Chaboureau, J.-P., Chauvigné, A., Chazette, P., Denjean, C., Desboeufs, K., Doussin, J.-F., Elguindi, N., Feuerstein, S., Gaetani, M., Giorio, C., Klopper, D., Mallet, M. D., Nabat, P., Monod, A., Solmon, F., Namwoonde, A., Chikwililwa, C., Mushi, R., Welton, E. J., and Holben, B.: The Aerosols, Radiation and Clouds in southern Africa (AEROCLO-sA) field campaign in Namibia: overview, illustrative observations and way forward, B. Am. Meteorol. Soc., 100, 1277–1298, https://doi.org/10.1175/BAMS-D-17-0278.1, 2019. a, b, c, d
Friedman, J., Hastie, T., and Tibshirani, R.: Regularization Paths for Generalized Linear Models via Coordinate Descent, J. Stat. Softw., 33, 1–22, 2010. a
Fuchs, J., Cermak, J., Andersen, H., Hollmann, R., and Schwarz, K.: On the Influence of Air Mass Origin on Low-Cloud Properties in the Southeast Atlantic, J. Geophys. Res.-Atmos., 122, 11076–11091, https://doi.org/10.1002/2017JD027184, 2017. a, b
Fuchs, J., Cermak, J., and Andersen, H.: Building a cloud in the southeast Atlantic: understanding low-cloud controls based on satellite observations with machine learning, Atmos. Chem. Phys., 18, 16537–16552, https://doi.org/10.5194/acp-18-16537-2018, 2018. a, b
Garstang, M., Tyson, P. D., Swap, R., Edwards, M., Kållberg, P., and Lindesay, J. A.: Horizontal and vertical transport of air over southern Africa, J. Geophys. Res.-Atmos., 101, 23721–23736, https://doi.org/10.1029/95JD00844, 1996. a
Goldreich, Y. and Tyson, P. D.: Diurnal and Inter-Diurnal Variations in Large-Scale Atmospheric Turbulence Over Southern Africa, S. Afr. Geogr. J., 70, 48–56, https://doi.org/10.1080/03736245.1988.10559755, 1988. a
Gordon, H., Field, P. R., Abel, S. J., Dalvi, M., Grosvenor, D. P., Hill, A. A., Johnson, B. T., Miltenberger, A. K., Yoshioka, M., and Carslaw, K. S.: Large simulated radiative effects of smoke in the south-east Atlantic, Atmos. Chem. Phys., 18, 15261–15289, https://doi.org/10.5194/acp-18-15261-2018, 2018. a
Gottlieb, T. R., Eckardt, F. D., Venter, Z. S., and Cramer, M. D.: The contribution of fog to water and nutrient supply to Arthraerua leubnitziae in the central Namib Desert, Namibia, J. Arid Environ., 161, 35–46, https://doi.org/10.1016/j.jaridenv.2018.11.002, 2019. a
Goubanova, K., Illig, S., Machu, E., Garçon, V., and Dewitte, B.: SST subseasonal variability in the central Benguela upwelling system as inferred from satellite observations (1999-2009), J. Geophys. Res.-Oceans, 118, 4092–4110, https://doi.org/10.1002/jgrc.20287, 2013. a
Gultepe, I., Tardif, R., Michaelides, 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, https://doi.org/10.1007/s00024-007-0211-x, 2007. a, b
Haensler, A., Cermak, J., Hagemann, S., and Jacob, D.: Will the southern african west coast fog be affected by future climate change? Results of an initial fog projection using a regional climate model, Erdkunde, 65, 261–275, https://doi.org/10.3112/erdkunde.2011.03.04, 2011. a
Henschel, J. R., Wassenaar, T. D., Kanandjembo, A., Louw, M. K., Neef, G., Shuuya, T., and Soderberg, K.: Roots point to water sources of Welwitschia mirabilis in a hyperarid desert, Ecohydrology, 12, e2039, https://doi.org/10.1002/eco.2039, 2018. a
Hersbach, H.: The ERA5 Atmospheric Reanalysis, AGU Fall Meeting Abstracts, 12–16 December 2016, San Fransisco, 2016. a
Hutchings, L., van der Lingen, C. D., Shannon, L. J., Crawford, R. J., Verheye, H. M., Bartholomae, C. H., van der Plas, A. K., Louw, D., Kreiner, A., Ostrowski, M., Fidel, Q., Barlow, R. G., Lamont, T., Coetzee, J., Shillington, F., Veitch, J., Currie, J. C., and Monteiro, P. M.: The Benguela Current: An ecosystem of four components, Prog. Oceanogr., 83, 15–32, https://doi.org/10.1016/j.pocean.2009.07.046, 2009. a, b
Iacobellis, S. F. and Cayan, D. R.: The variability of California summertime marine stratus: Impacts on surface air temperatures, J. Geophys. Res.-Atmos., 118, 9105–9122, https://doi.org/10.1002/jgrd.50652, 2013. a, b
James, R. and Washington, R.: Changes in African temperature and precipitation associated with degrees of global warming, Climatic Change, 117, 859–872, https://doi.org/10.1007/s10584-012-0581-7, 2013. a
Johnstone, J. A. and Dawson, T. E.: Climatic context and ecological implications of summer fog decline in the coast redwood region, P. Natl. Acad. Sci. USA, 107, 4533–4538, https://doi.org/10.1073/pnas.0915062107, 2010. a
Kaseke, K. F., Tian, C., Wang, L., Seely, M., Vogt, R., Wassenaar, T., and Mushi, R.: Fog spatial distributions over the central Namib Desert – An isotope approach, Aerosol Air Qual. Res., 18, 49–61, https://doi.org/10.4209/aaqr.2017.01.0062, 2018. a
Klein, S. A. and Hartmann, D. L.: The Seasonal Cycle of Low Stratiform Clouds, J. Climate, 6, 1587–1606, https://doi.org/10.1175/1520-0442(1993)006<1587:TSCOLS>2.0.CO;2, 1993. a, b
Knippertz, P. and Martin, J. E.: Tropical plumes and extreme precipitation in subtropical and tropical West Africa, Q. J. Roy. Meteorol. Soc., 131, 2337–2365, https://doi.org/10.1256/qj.04.148, 2005. a
Koračin, D., Dorman, C. E., and Dever, E. P.: Coastal perturbations of marine-layer winds, wind stress, and wind stress curl along California and Baja California in June 1999, J. Phys. Oceanogr., 34, 1152–1173, https://doi.org/10.1175/1520-0485(2004)034<1152:CPOMWW>2.0.CO;2, 2004. a
Koračin, D., Businger, J. A., Dorman, C. E., and Lewis, J. M.: Formation, evolution, and dissipation of coastal sea fog, Bound.-Lay. Meteorol., 117, 447–478, https://doi.org/10.1007/s10546-005-2772-5, 2005. a, b
Lentz, S. J.: The Surface Boundary Layer in Coastal Upwelling Regions, J. Am. Meteorol. Soc., 22, 1517–1539, https://doi.org/10.1175/1520-0485(1992)022<1517:tsblic>2.0.co;2, 1992. a
Maúre, G., Pinti, I., Ndebele-Murisa, M., Muthige, M., Lennard, C., Nikulin, G., Dosio, A., and Meque, A.: The southern African climate under 1.5 ∘C and 2 ∘C of global warming as simulated by CORDEX regional climate models, Environ. Res. Lett., 13, 065002, https://doi.org/10.1088/1748-9326/aab190, 2018. a
Nieman, W. A., Heyns, C., and Seely, M. K.: A note on precipitation at Swakopmund, Madoqua, 11, 69–73, 1978. a
Olivier, J. and Stockton, P. L.: The influence of upwelling extent upon fog incidence at Lüderitz, southern Africa, Int. J. Climatol., 9, 69–75, https://doi.org/10.1002/joc.3370090106, 1989. a, b, c, d, e, f
Oueslati, B., Pohl, B., Moron, V., Rome, S., and Janicot, S.: Characterization of heat waves in the Sahel and associated physical mechanisms, J. Climate, 30, 3095–3115, https://doi.org/10.1175/JCLI-D-16-0432.1, 2017. a
Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., Blondel, M., Louppe, G., Prettenhofer, P., Weiss, R., Dubourg, V., Vanderplas, J., Passos, A., Cournapeau, D., Brucher, M., Perrot, M., and Duchesnay, É.: Scikit-learn: Machine Learning in Python, J. Mach. Learn. Res., 12, 2825–2830, 2011. a
Pham-Thanh, H., Linden, R., Ngo-Duc, T., Nguyen-Dang, Q., Fink, A. H., and Phan-Van, T.: Predictability of the rainy season onset date in Central Highlands of Vietnam, Int. J. Climatol., 6383, https://doi.org/10.1002/joc.6383, 2019. a
Roth-Nebelsick, A., Ebner, M., Miranda, T., Gottschalk, V., Voigt, D., Gorb, S., Stegmaier, T., Sarsour, J., Linke, M., and Konrad, W.: Leaf surface structures enable the endemic Namib desert grass Stipagrostis sabulicola to irrigate itself with fog water, J. Roy. Soc. Interface, 9, 1965–74, https://doi.org/10.1098/rsif.2011.0847, 2012. a
Schmetz, J., Pili, P., Tjemkes, S., Just, D., Kerkmann, J., Rota, S., and Ratier, A.: An Introduction to Meteosat Second Generation (MSG), B. Am. Meteorol. Soc., 83, 977–992, 2002. a
Seely, M. K. and Henschel, J. R.: The Climatology of Namib Fog, Proceedings, First International Conference on Fog and Fog Collection, Vancouver, Canada, 353–356, 1998. a, b, c, d, e, f, g, h, i, j, k
Seely, M. K., De-Vos, M. P., and Louw, G. N.: Fog imbibition satellite fauna and unusual leaf structure in a namib desert dune plant trianthema hereroensis, S. Afr. J. Sci., 73, 169–172, 1977. a
Sippel, S., Meinshausen, N., Merrifield, A., Lehner, F., Pendergrass, A. G., Fischer, E., and Knutti, R.: Uncovering the Forced Climate Response from a Single Ensemble Member Using Statistical Learning, J. Climate, 32, 5677–5699, https://doi.org/10.1175/jcli-d-18-0882.1, 2019. a
Spirig, R., Vogt, R., Larsen, J. A., Feigenwinter, C., Wicki, A., Parlow, E., Adler, B., Kalthoff, N., Cermak, J., Andersen, H., Fuchs, J., Bott, A., Hacker, M., Wagner, N., Maggs-Kölling, G., Wassenaar, T., and Seely, M.: Probing the fog life-cycles in the Namib desert, B. Am. Meteorol. Soc., https://doi.org/10.1175/BAMS-D-18-0142.1, 2019. a, b, c, d, e
Swap, R. J., Annegarn, H. J., Suttles, J. T., King, M. D., Platnick, S., Privette, J. L., and Scholes, R. J.: Africa burning: A thematic analysis of the Southern African Regional Science Initiative (SAFARI 2000), J. Geophys. Res.-Atmos., 108, 1–15, https://doi.org/10.1029/2003JD003747, 2003. a
Tim, N., Zorita, E., and Hünicke, B.: Decadal variability and trends of the Benguela upwelling system as simulated in a high-resolution ocean simulation, Ocean Sci., 11, 483–502, https://doi.org/10.5194/os-11-483-2015, 2015. a
Warren-Rhodes, K. A., Mckay, C. P., Boyle, L. N., Wing, M. R., Kiekebusch, E. M., Cowan, D. A., Stomeo, F., Pointing, S. B., Kaseke, K. F., Eckardt, F., Henschel, J. R., Anisfeld, A., Seely, M., and Rhodes, K. L.: Physical ecology of hypolithic communities in the central Namib Desert: The role of fog, rain, rock habitat, and light, J. Geophys. Res.-Biogeo., 118, 1451–1460, https://doi.org/10.1002/jgrg.20117, 2013. a
Yu, J.-Y. and Kim, S. T.: Relationships between Extratropical Sea Level Pressure Variations and the Central Pacific and Eastern Pacific Types of ENSO, J. Climate, 24, 708–720, https://doi.org/10.1175/2010JCLI3688.1, 2010. a
Zhou, C. and Penner, J. E.: Why do general circulation models overestimate the aerosol cloud lifetime effect? A case study comparing CAM5 and a CRM, Atmos. Chem. Phys., 17, 21–29, https://doi.org/10.5194/acp-17-21-2017, 2017. a
Zuidema, P., Redemann, J., Haywood, J., Wood, R., Piketh, S., Hipondoka, M., Formenti, P., Zuidema, P., Redemann, J., Haywood, J., Wood, R., Piketh, S., Hipondoka, M., and Formenti, P.: Smoke and Clouds above the Southeast Atlantic: Upcoming Field Campaigns Probe Absorbing Aerosol's Impact on Climate, B. Am. Meteorol. Soc., 97, 1131–1135, https://doi.org/10.1175/BAMS-D-15-00082.1, 2016. a, b