Atmospheric boundary layer in the Atlantic: the desert dust impact

Atmos. Chem. Phys., 25, 16491–16510, 2025 https://doi.org/10.5194/acp-25-16491-2025 © Author(s) 2025. This work is distributed under the Creative Commons Attribution 4.0 License. Research article Atmospheric boundary layer in the Atlantic: the desert dust impact Ioanna Tsikoudi1,2, Eleni Marinou1, Maria Tombrou2, Eleni Giannakaki2, Emmanouil Proestakis1, Konstantinos Rizos1, Ville Vakkari3,4, Holger Baars5, Annett Skupin5, Ronny Engelmann5, Zhenping Yin5, and Vassilis Amiridis1 1National Observatory of Athens (NOA), Institute for Astronomy, Astrophysics, Space Applications and Remote Sensing (IAASARS), Athens, Greece 2Department of Physics, National and Kapodistrian University of Athens, Athens, Greece 3Finnish Meteorological Institute (FMI), Helsinki, Finland 4Atmospheric Chemistry Research Group, Chemical Resource Beneficiation, North-West University, Potchefstroom, South Africa 5Leibniz Institute for Tropospheric Resear (TROPOS), Leipzig, Germany Correspondence: Ioanna Tsikoudi ([email protected]) Received: 7 March 2025 – Discussion started: 14 March 2025 Revised: 15 October 2025 – Accepted: 16 October 2025 – Published: 24 November 2025 Abstract. We investigate the dynamics of the atmospheric Boundary Layer (BL) over the Atlantic Ocean, with a focus on the region surrounding Cabo Verde during the Joint Aeolus Tropical Atlantic Campaign (JATAC) and the ASKOS experiment, using a combination of ground-based PollyXT and Doppler lidars, satellite lidar data from Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), radiosondes, and the model outputs of the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). The comparison of CALIPSO lidar results with ECMWF/IFS reanalysis for 2012–2022, revealed good agreement for BL top over open ocean regions but weaker relation over dust-affected areas of the African continent. In these regions, daytime CALIPSO retrievals typically indicate lower BL tops than ECMWF, while at night CALIPSO often detects aerosols within the residual layer, leading to higher estimates than the model. Observations in Cabo Verde highlight distinctive Marine Atmospheric Boundary Layer (MABL) characteristics, such as limited diurnal evolution, but also show the potential for BL heights to reach up to 1 km, driven by factors like strong winds that increase mechanical turbulence. Additionally, the technical and physical challenges in estimating the BL height using different datasets and methods are discussed, examining cases with different thermodynamical conditions and aerosol load that directly affect the dynamics of the BL. The findings underline the strengths and limitations of different observational and modeling approaches, and emphasizes on the importance of considering local meteorology and aerosol conditions when interpreting BL height. 1 Introduction The atmospheric Boundary Layer (BL) is characterized by complex interactions between surface-driven forces and meteorological conditions, which determine its height, structure, and the degree of turbulent mixing within (Stull, 1988). BL dynamics vary considerably across different environments, presenting challenges for weather modeling and prediction, especially in transitional zones like those between deserts and oceans (Seibert et al., 2000; Li et al., 2017). Monitoring the BL top reliably is a challenge, particularly in heterogeneous environments where traditional observation methods may fall short. Lidar systems have proven valuable for continuous profiling of aerosol and atmospheric structures, as their high vertical resolution enables detailed monitoring of BL height (Wiegner et al., 2006; Baars et al., Published by Copernicus Publications on behalf of the European Geosciences Union. 16492 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact 2008). Yet, automatic identification of the BL top from lidar data is challenging in complex areas, because BL structures can be influenced by surface type, time of day, and atmospheric stability (Tsikoudi et al., 2022). Up to now, lidarbased BL retrievals showed very good performance on relatively predictable areas with known BL patterns, such as open land surfaces or stable atmospheric conditions (Tsaknakis et al., 2011; Seidel et al., 2012). Expanding lidar BL retrievals to more complex environments, is an ongoing challenge especially when it comes to oceanic and coastal BLs where ground-based observation sites are limited. Over the open Atlantic, the Marine Atmospheric Boundary Layer (MABL) is typically shallow and influenced by the relatively constant sea surface temperature, while boundary layers in coastal and island regions experience terrestrialmarine interactions that increase their variability (Garratt, 1994; Wood, 2012). A limited number of studies over years have addressed the detection and analysis of MABL using lidar data, primarily due to practical and observational challenges over the ocean (e.g. Atlas et al., 1986; Flamant et al., 1997; Pena et al., 2015). Given these constraints, satellite observations can provide an important means of obtaining information in remote regions lacking in-situ and ground-based remote sensing data, while also enabling the development of global climatologies (Teixeira et al., 2025). Although the BL is a near-surface phenomenon, several satellite measurements can indirectly infer its properties, particularly its depth and spatial or temporal variability. The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, has been widely used to derive global BL height climatologies over ocean and land and is therefore essential for studying lower troposphere characteristics (e.g. Liu et al., 2024). Nevertheless, when interpreting satellite-derived BL characteristics, it is crucial to decode the measurements appropriately, as the definition and identification of the BL can vary depending on the chosen approach and physical parameter. The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) of the CALIPSO satellite, can measure, among others, backscattered light from aerosols and clouds. Hence, in this case, the top of the lowest aerosol layer often coincides with the BL top, since aerosols are typically well mixed within the BL and drop sharply above it (Li et al., 2017). The general circulation over the tropical Atlantic is dominated by the Inter-Tropical Convergence Zone (ITCZ) and affected by the presence of the Saharan Air Layer (SAL). The SAL is a typically warm and dry air layer that frequently occurs at large scales in the tropical North Atlantic Ocean and can reside up to 5 km in altitude, often accompanied by dust aerosols (Carlson and Prospero, 1972; Dunion and Velden, 2004; Wu, 2007). The ITCZ, migrates seasonally between the northern and southern tropics, influencing rainfall and convective activity, creating conditions conducive to both the formation of clouds and the aerosol convection over the Atlantic (Zhou et al., 2020). In tandem, the SAL, comprising of hot, dry air laden with desert dust from the Sahara, moves westward across the Atlantic Ocean, especially in summer, driven by the prevailing trade winds (Prospero and MayolBracero, 2013) and has consequences on the surface radiation budget (Evan et al., 2009; Yu et al., 2006). These circulation patterns are key in transporting dust from Africa to the Atlantic, affecting the radiative balance and potentially impacting cloud formation, atmospheric stability, and therefore BL behavior in the region (Sun and Zhao, 2020). A typical characteristic of the eastern sides of the Atlantic, is that the air subsiding into the subtropical north-east Atlantic is warmer and drier than the air that has been in contact with the relatively cold ocean surface influenced by upwelling, and a strong inversion forms at the interface of the two air masses (Hanson, 1991). As such, transported desert dust from Africa introduce another layer of complexity in tropospheric dynamics and clouds activity by altering radiation budget, atmospheric stability, and moisture distribution (e.g. Marinou et al., 2021; Ansmann et al., 2017; Marsham et al., 2008). This dual effect of dust – scattering and absorbing solar radiation while in the same time serving as cloud condensation and ice nucleation nuclei (CCN/IN) – leads to competing influences on the BL (e.g. radiative cooling can suppress turbulent mixing, yet CCN activation can lead to increased cloud cover and associated feedback on surface radiation). These processes have been observed to influence the vertical structure and stability of the BL, but their overall impact on BL dynamics is still not fully understood. Accurately representing BL–aerosol interactions in climate and chemical transport models is crucial because these processes affect surface conditions and large-scale atmospheric circulation (Menut et al., 2009; Pérez et al., 2006; Tombrou et al., 2015, 2007). Gaps in observational data over complex environments, such as the dust-laden, desertocean transition zone in the Atlantic, limit the model’s ability to accurately capture BL evolution and aerosol influences (Rémy et al., 2019; Rémy et al., 2022; Kallos et al., 2007). The need for observational data to validate and refine these models is high, especially given the impacts on cloud formation, energy distribution, and surface-air interactions. Addressing these gaps through both ground-based experimental campaigns such as Joint Aeolus Tropical Atlantic Campaign (JATAC) and satellite sensors such as space Lidars can significantly enhance understanding and modeling of BL processes in regions of critical climatic importance. In addition to investigating BL–aerosol interactions, this study aims to improve BL top detection methods in diverse and complex environments. By addressing challenges inherent to automated BL detection, particularly in areas affected by aerosols and variable atmospheric conditions, this work contributes to the development of more robust methods for BL identification. The structure of this paper is as follows: Sect. 2 provides an overview of the datasets and methods used, including ground-based lidar, space lidar, radiosonde data, and model Atmos. Chem. Phys., 25, 16491–16510, 2025 https://doi.org/10.5194/acp-25-16491-2025 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact 16493 outputs. Section 3 examines the BL characteristics across different environments, beginning with the Atlantic Ocean (Area 1) and the ocean-desert transition zone (Area 2), before focusing on Cabo Verde, where dust interactions with the BL are investigated. Finally, Sect. 4 presents the main conclusions of this study. 2 Data sources and analysis This study uses data from the ASKOS Campaign (Marinou et al., 2023), the ground-based component of the JATAC organised by the European Space Agency (ESA). The campaign was conducted at the Ocean Science Centre Mindelo (OSCM), at the island of São Vicente, Cabo Verde, during 2021–2022. In addition, CALIPSO observations and ECMWF model data are employed. The BL height is derived using the gradient method and the Wavelet Covariance Transform method on satellite and ground-based lidar data respectively. 2.1 Datasets The comprehensive ASKOS dataset includes active remotesensing observations and radiosonde profiles, both essential for characterizing atmospheric dynamics in the studied region. Specifically, ground-based PollyXT lidar and Wind Doppler lidar observations are examined, together with the LIVAS (LIdar climatology of Vertical Aerosol Structure for space-based lidar simulation studies) Climate Data Record (Amiridis et al., 2015) derived from CALIPSO. Furthermore, the BL heights obtained from measurements are compared against the ERA5 reanalysis dataset, produced with the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF), at 0.25° ×0.25° horizontal resolution. 2.1.1 Groundbased lidars The ground-based PollyXT Raman Lidar (Engelmann et al., 2016), consists of a compact, pulsed Nd:YAG laser, emitting at 355, 532, and 1064 nm at a 20Hz repetition rate, with the laser beam pointed into the atmosphere at an off-zenith angle of 5°. The backscattered signal is collected by a Newtonian telescope with a 0.9 m focal length, acquiring profiles with a vertical resolution of 7.5m, and a temporal resolution of 30s. The system was operated by the Leibniz Institute for Tropospheric Research (TROPOS) during the ASKOS Campaign, providing data coverage for the entire campaign. Figure 1 presents some indicative PollyXT measurements. Specifically, the attenuated backscatter coefficient of the 1064nm channel (Att BSC, Fig. 1, left panel) is examined to derive the BL top, and the volume linear depolarization ratio (VLDR, Fig. 1, right panel) is complementary investigated to infer the aerosol shape. The white points in the attenuated backscatter indicate the presence of clouds and were not included in the BL analysis. Additionally, complementary data from a Halo Photonics Stream Line scanning Doppler lidar were used to examine the horizontal wind speed and direction, as well as the vertical wind component. This lidar is a 1.5µm pulsed Doppler lidar with a heterodyne detector (Pearson et al., 2009). The Doppler lidar has a range resolution of 48 m and measures the attenuated aerosol backscatter and Doppler velocity along the beam direction. Horizontal wind profiles were retrieved from a velocity azimuth display (VAD) scan with 12 azimuthal angles at 60° elevation angle every 15min. Otherwise, the Doppler lidar operated in vertical stare mode, retrieving vertical wind profile time series. The Doppler lidar data was post-processed according to Vakkari et al. (2019) and a signal-to-noise ratio (SNR) threshold of 0.005 was applied to the vertically-pointing measurements. Turbulent kinetic energy (TKE) dissipation rate profiles were calculated from the vertically-pointing data using the method by O’Connor et al. (2010). Instrumental noise was calculated from signal-to-noise ratio according to Pearson et al. (2009) and subtracted from the vertical wind variance time series before the TKE dissipation rate calculation. To estimate mixed layer height (MLH) from the TKE dissipation rate profiles a threshold of 10−4m2s−3was applied, similar to previous studies (e.g. Vakkari et al., 2015). 2.1.2 Space lidar: CALIPSO–CALIOP Towards investigating the dynamics of the BL over the Atlantic Ocean and parts of West Africa, observations of the Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP; Hunt et al., 2009), the primary instrument on board the joint National Aeronautics and Space Administration (NASA) and Centre National D’Études Spatiales (CNES) Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) mission (Winker et al., 2010), are extensively used. CALIOP provided as integrated component of the Afternoon-Train constellation of polar-orbit sun-synchronous satellites (Stephens et al., 2018), profiles of aerosols and clouds along the CALIPSO orbitpath between June 2006 and August 2023. In the framework of the study, CALIOP Level 2 (L2) Version 4 (V4) aerosol profiles (APro) of backscatter coefficient at 532 nm and particulate depolarization ratio at 532 nm are used, provided at uniform 5 km horizontal resolution and 60 m vertical resolution for the altitudinal range between −0.5 and 20.2 km a.m.s.l. (above mean sea level), for the domain encompassing the broader North Atlantic Ocean–Western Saharan Desert. Prior implementation of CALIOP optical products, rigorous quality assurance procedures are applied (Marinou et al., 2017; Proestakis et al., 2024), following also the quality controls adopted towards the generalization of the official CALIPSO Level 3 (L3) products (Winker et al., 2013; Tackett et al., 2018). Towards this objective, the most aghttps://doi.org/10.5194/acp-25-16491-2025 Atmos. Chem. Phys., 25, 16491–16510, 2025 16494 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact Figure 1. Ground-based PollyXT Lidar at Mindelo (16.87° N, 24.99° W), Cabo Verde, on 10 September 2021, depicting the attenuated backscatter coefficient (Att Bsc) at 1064 nm (left panel), and volume linear depolarization ratio (VLDR) at 532 nm (right panel). gressive quality control procedure applied in the framework of the study is the cloud-free condition, removing the entire L2 profiles when detected atmospheric layers (Vaughan et al., 2009) along the CALIPSO orbit-path are classified as clouds in the feature-type classification algorithm (Liu et al., 2009; Zeng et al., 2019). Figure 2 provides an indicative example of the considered CALIOP observations and products, and more specifically the Feature Type (Fig. 2, top left panel) product and the profiles of particulate depolarization ratio at 532 nm (Fig. 2, top right panel), total backscatter coefficient at 532 nm (Fig. 2, bottom left panel), and quality-assured total backscatter coefficient at 532nm (Fig. 2 bottom right), along the CALIPSO overpass on the 10 September 2021. For this analysis, CALIPSO data from September 2012–2022 were employed. 2.1.3 Radiosondes and models Radiosonde profiles were analyzed to examine the dynamic structure of the lower troposphere and to evaluate the remote sensing measurements conducted during the intensive phase of ASKOS Campaign (June and September 2022). The GRAW DFM-09 radiosondes were launched to provide real-time, high-resolution measurements of temperature, humidity, and wind, which are essential for identifying the BL characteristics, such as height, stability, and thermodynamic properties. The sensors were equipped with a GPS receiver and transmit data via a radiofrequency link to the ground station. The measurements-derived BL height results are compared to values obtained from the ERA5 Reanalysis dataset, produced by the ECMWF/IFS. The ERA5 data, available at a horizontal resolution of 0.25° ×0.25° with 137 vertical levels (Vogelezang and Holtslag, 1996), offers a consistent representation of atmospheric conditions. The BL height in the ERA5 reanalysis dataset, cannot be explicitly resolved, as turbulent processes occur on scales smaller than the model grid. Instead, it is diagnosed from boundary-layer theory based on the critical Richardson number (Ric), and the parametrization of the mixed layer in the model uses a BL height from an entraining parcel model. Though, in order to get a continuous field, also in neutral and stable situations, a bulk Richardson method is used as a diagnostic, independent of the turbulence, parametrization. This method follows the conclusions of the study by Seidel et al. (2012), who showed that this algorithm is suitable for both convective and stable boundary layers, identifying a nonnegative height in all cases, and is not strongly dependent on the sounding vertical resolution. Several approximations are applied to the original algorithm – such as ignoring surface friction and setting winds near the surface (2m) to zero for radiosonde data – so that the bulk Richardson number can be consistently used to define the BL height as the lowest level where it reaches the critical value (Ric) of 0.25. The detailed description can be found at ECMWF (2017) (Chap. 3). Additionally, Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) is also considered to analyze the backward trajectories of air masses arriving at the site of the ASKOS Campaign, Mindelo, Cabo Verde. This model estimates the tracking of air parcels over time, providing valuable information about the origins of the air parcels and their potential interactions with dust and other atmospheric constituents (Rolph et al., 2017). By identifying these pathways, a clearer understanding of the sources and transport mechAtmos. Chem. Phys., 25, 16491–16510, 2025 https://doi.org/10.5194/acp-25-16491-2025 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact 16495 Figure 2. CALIPSO nighttime overpass in the ESA-ASKOS campaign region of interest in the proximity of Cabo Verde on 10 September 2021, depicting the Feature Type (top left panel), particulate depolarization ratio at 532 nm (top right panel), total backscatter coefficient at 532 nm (bottom left panel), and the quality-assured total backscatter coefficient at 532 nm (bottom right panel). anisms of the atmospheric conditions at Cabo Verde can be established. 2.2 Boundary Layer top retrieval methods As mentioned above, only cloud-free lidar profiles are used in this study. The CALIPSO backscatter coefficient profiles at 532nm (Fig. 2, bottom right) were horizontally averaged over ±100m along the satellite trajectory around each point of interest. For instance, if the satellite crosses latitude 16.87° N, all available profiles within 200m of the trajectory are averaged. This approach improves the signal-to-noise ratio (SNR) of the measurements and facilitates comparison with coarser-resolution datasets, such as model reanalysis outputs. The BL height is then retrieved from the CALIPSO profiles using the gradient method (Li et al., 2021). The ground-based PollyXT lidar backscatter profiles at 1064 nm (Fig. 1, left panel) are averaged in time, since the instrument is stationary. For this study, profiles were averaged over ±15 or ±30 min around the time of interest, depending on the scene homogeneity. The BL top is retrieved using the Haar Wavelet Covariance Transform (WCT) method (Brooks, 2003). The dilation factor α, which defines the effective window size of the wavelet, was empirically set to 100 m, corresponding to approximately 13 vertical bins, given the PollyXT vertical resolution of 7.5 m. In most cases this choice provided a good balance between sensitivity to sharp gradients and noise reduction, although in some situations the dilation factor was adjusted to better capture the layering. The Haar integration is performed symmetrically above and below each altitude level over the chosen dillation window, with integration limits of ±α/2 around the center altitude. For the radiosonde data, layer detection is performed using the gradient method applied to virtual potential temperature (θV) and relative humidity (RH) (Seidel et al., 2010). The θV, which accounts for moisture effects on air density, provides a reliable representation of buoyancy and atmospheric stability, while RH typically exhibits a sharp gradient near the BL top in the persistently humid environment of São Vicente Island. In some cases, particularly under high aerosol loading, identification of the BL top also requires visual inspection to accurately locate the inversion cap. Figure 3 shows profiles of the backscatter coefficient at 1064 nm from ground-based PollyXT (left panel), of the RH from radiosonde (middle panel) and backscatter coefficient at 532nm from CALIPSO satellite lidar (right panel) from 23 September 2022, around 19:30 UTC. The dotted grey lines represent the method applied for detecting BL top, namely WCT method for PollyXT Lidar and Gradient method for the rest two. A local maximum of the wavelet profile for WCT method, and a local minimum of the gradihttps://doi.org/10.5194/acp-25-16491-2025 Atmos. Chem. Phys., 25, 16491–16510, 2025 16496 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact ent for the gradient method, represent steep reduction in the investigated variable (orange dashed lines). Several significant challenges arise when identifying the BL top with lidars (both ground-based and satellite), particularly in aerosol-complex environments. For a satellite-based lidar like CALIOP, the signal can become highly attenuated as it approaches the Earth’s surface, due to the existence of clouds above the BL. The weakened return signals also result from longer travel distances from the satellite platform to the earth’s surface, which lead to a lowered SNR. This can compromise the reliability of detecting lower tropospheric features and lead to inaccurate identification of the BL top. To mitigate this, (i) only cloud-free profiles were selected to ensure data quality, though this restriction reduces the dataset and introduces observational limitations, and (ii) averaged profiles were considered to increase the SNR. Additionally, in marine environments, cumulus clouds frequently form at the BL top, which can serve as a useful, albeit indirect, marker for BL height for ground-based lidars that can detect the cloud base. Moreover, if a thin cumulus cloud is present above the BL top and allows partial laser penetration, the WCT may incorrectly identify the cloud’s upper boundary as the BL top instead of the actual BL height. A similar issue occurs in the presence of dust layers, as the WCT detects reductions in the lidar signal caused by these layers. This can lead to misclassification of the dust layer boundaries as the BL top, complicating the accurate identification of the atmospheric structure. These limitations underscore the need for visual inspection to ensure accuracy in identifying the BL top in such settings, as automated methods may struggle to locate the correct layering. 3 Boundary Layer characteristics in diverse environments The characteristics of the BL during JATAC Campaign are examined across the contrasting environments depicted in Fig. 4: over the Atlantic Ocean (blue rectangle – Area 1), within the ocean-desert transition zone (orange rectangle – Area 2), and at the area of São Vicente in Cabo Verde (red circle). The Sahara Desert and the Atlantic Ocean are characterized by distinct conditions in terms of weather, aerosol concentrations, and therefore atmospheric dynamics. These variations are anticipated to influence respectively the structure and evolution of BL in the Areas of Fig. 4. The lower troposphere above the Atlantic Ocean is rich in marine aerosols, and presents relatively stable meteorological conditions, typical for open-ocean broad-scale circulations (Croft et al., 2021). In contrast, the lower troposphere over the desert is characterized by high dust aerosol concentrations, intense solar heating, and variable atmospheric stability (Giménez et al., 2010). The border region between ocean and desert introduces an interaction zone where different aerosols co-exist in big concentrations, producing unique BL characteristics due to the convergence of these differing air masses. Moreover, the existence of SAL has an impact on the on the surface radiation budget (Evan et al., 2009) and hence on the sea surface temperature (SST). Foltz and McPhaden (2008) found that Saharan dust outflows at the Tropical North Atlantic, were consistently associated with a reduction in solar radiation, with approximately 35 % of SST variability attributed to dust outbreaks, while other SST cooling anomalies were linked to wind stress. The dust aerosol effect on SST depends on several factors, such as the temperature contrast between the dust layer and SST, the characteristics of the dust layer, concentration and altitude (Luo et al., 2021). 3.1 Analysis of Area 1: the BL in the Atlantic Ocean The Atlantic Ocean is characterized by dynamic weather systems and cyclonic activity, incorporating continuous exchange of heat and moisture between the sea surface and the adjacent air parcel (Schnitker, 1982). In open ocean areas such as Area 1, there is no direct interaction of the lower troposphere and the land, allowing for the development of a MABL. The MABL contains higher humidity levels and the airflow is smoother due to reduced friction from the water surface, comparing to land. Wind and temperature profiles in the MABL are mainly influenced by sea surface temperature, oceanic currents and large-scale atmospheric circulation. In this section, we focus on the MABL characteristics within the blue rectangle of Area 1 (Fig. 4). 10 years of CALIOP data (2012–2022) are examined, using only the profiles recorded in month September. By limiting the data to one month, we aim to achieve more homogeneous conditions to better capture the prevailing environmental characteristics (e.g. relatively consistent sea surface temperatures). Figure 5 (left panel) illustrates the conceptual trajectories of the CALIPSO satellite across the study area. The analysis investigated cloud-free averaged profiles measured within 100 km around latitude 16.87° N, corresponding to the latitude of ground-based measuring site in Cabo Verde, as represented by the red points in Fig. 5 (left panel). A total of 449 profiles from nighttime and daytime CALIPSO trajectories (conceptually indicated in green and purple, respectively) are analyzed across longitudes from 60 to 25°W. The spatial range of 100km is suitable for capturing representative MABL characteristics in the study area because the selected profiles are cloud-free and measured over the ocean surface, maintaining generally homogeneous conditions of temperature, and humidity. For each profile, the derivative of the backscatter-coefficient profile at 532 nm is calculated (as in Fig. 3, right panel) and the minima are constrained at the lower 3km. The results of the MABL analysis from the space lidar data are compared with BL heights derived from the ECMWF dataset. To account for longitudinal time differences, each profile’s measurement time is converted to local Atmos. Chem. Phys., 25, 16491–16510, 2025 https://doi.org/10.5194/acp-25-16491-2025 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact 16497 Figure 3. 23 September 2022, around 19:30 UTC: profiles of atmospheric variables and their corresponding detection methods for determining the boundary layer (BL) top. The solid blue lines represent the products of measurements, while the dotted grey lines correspond to the applied methods for BL top detection. Left panel: backscatter coefficient at 1064 nm from the ground-based PollyXT lidar, middle panel: Relative Humidity (RH) from radiosonde, and right panel: backscatter coefficient at 532 nm from the CALIPSO satellite lidar. The grey shading in the lidar profiles represent the standard deviation resulting from the averaging. The selected BL top is highlighted by the orange dashed lines. The shading around these lines corresponds to ±50 m. Figure 4. Map displaying the study areas for BL analysis: the blue rectangle (Area 1) represents the open-ocean Marine Atmospheric Boundary Layer (MABL) discussed in Sect. 3.1. The orange rectangle (Area 2) marks a transition zone at the ocean-desert interface, analysed in Sect. 3.2. The red circle is the ground-based measurements site at the Ocean Science Center Mindelo (OSCM) in Cabo Verde (Sect. 3.3). time based on its longitude. For each lidar profile, a temporally and spatially matched ECMWF point at the same local time is selected for direct comparison. The findings are presented in Fig. 5 (right panel). The blue circles display the MABL top heights derived from CALIPSO profiles, averaged hourly in local time. The orange points represent the corresponding hourly-averaged BL top heights from ECMWF. The data points are clustered within the 00:00– 04:00 and 12:00–16:00 LT (local time) windows, because they correspond to CALIPSO’s nighttime and daytime overpasses in the studied region for the month of September. The BL top in Area 1 under cloud-free conditions consistently ranges between 600 and 800m a.s.l. (above sea level) in both datasets. There is a strong agreement in the mean BL heights between the two datasets, each exhibiting uncertainties of approximately 20 %, indicating that both provide comparable estimates of the boundary layer top. This agreement suggests that CALIPSO and ECMWF are consistent in representing the overall distribution of BL heights; however, as discussed in Appendix A1, their agreement at the level of individual profiles remains limited. While uncertainties associated with BL retrievals and time averaging may broaden the range of 600–800m for BL top, these results are consistent with the expected behavior of the MABL, which typically exhibits limited diurnal variation. The time-averaging uncertainties shown in the figure arise from the methods used to capture the BL in the two datasets. For CALIOP profiles, lidar-based retrievals inherently carry significant uncertainty and sensitivity due to measurement https://doi.org/10.5194/acp-25-16491-2025 Atmos. Chem. Phys., 25, 16491–16510, 2025 16498 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact Figure 5. Left panel: conceptual illustration of the trajectories of the CALIPSO satellite across the study area. Right panel: comparison of BL top derived from CALIPSO (blue points) and ECMWF (orange points) for 10 years (2012–2022) in Area 1. noise. Here, the BL top is derived using the gradient method and aerosol layers as discussed in Sect. 2; however, this method can occasionally detect layers that do not correspond to the actual PBL, introducing additional variability. In contrast, the model provides an averaged representation over a relatively large grid (0.25°, or approximately 27.8 km around 16° N), which may introduce variability but is less sensitive to small-scale fluctuations compared to CALIOP. Consequently, the model typically exhibits slightly lower standard deviations. 3.2 Analysis of Area 2: the BL in the Ocean-desert transition zone Area 2, highlighted by the orange rectangle in Fig. 4, spans within longitudes of 35–0° W: from the eastern Atlantic Ocean to the Western Africa, including the region around Cabo Verde. This area lies at the interface of two significantly different environments, as land and water interact differently with solar radiation due to their distinct heat capacities and reflective properties. On the West Africa land side, the lower troposphere directly interacts with the continental surface and the air is enriched with desert dust aerosols originating from the Sahara, where high temperatures, dry conditions, and strong winds are dominant. In contrast, the Eastern Atlantic Ocean side is predominantly influenced by marine aerosols within the lower troposphere, reflecting the ocean’s stable, moisture-laden environment. In terms of heat capacity, land absorbs and releases heat quickly, leading to larger temperature fluctuations, while water absorbs energy more gradually, storing and slowly releasing it. These sharp contrasts in meteorological conditions and aerosol composition across the Area 2, are expected to have a notable impact on the BL structure. For this analysis, similarly to Sect. 3.1, cloud-free averaged profiles were selected from the CALIPSO satellite lidar for September 2012–2022 to derive the BL top and are compared with the corresponding ECMWF data. Figure 6 presents the BL top results obtained from CALIPSO lidar measurements (blue points), and from the corresponding ECMWF points (orange) along the cross-section at latitude 16.87° N (the latitude of the Mindelo observatory). The CALIPSO trajectories are divided into daytime (Fig. 6, left panel) and nighttime (Fig. 6, right panel) intervals after converting to local time, to highlight the distinct patterns of BL during different phases of the diurnal cycle. In the daytime plot (Fig. 6, left panel), the two datasets show better agreement over the ocean compared to over land. Over land, the variability increases significantly for both CALIPSO and ECMWF, sometimes reaching up to 40% (e.g., at lon = −8°), particularly for the ECMWF dataset. This increased variability can be attributed to the diurnal evolution of the boundary layer: the data include all BL tops from 06:00 to 18:00 LT. Since the boundary layer over land grows and decays throughout these hours, typical for continental and desert areas (Garcia-Carreras et al., 2015), averaging over this period naturally results in large standard deviations. A similar behavior is observed in the CALIPSO retrievals, which also show substantial variability above land. It is also worth noting that CALIPSO tends to detect lower BL tops than ECMWF over land. This difference likely arises from the way to define the BL top: ECMWF relies on thermodynamic criteria, while CALIPSO identifies a decrease in aerosol concentration. Consequently, aerosols detected by CALIPSO are mostly confined within the mixed layer (Liu et al., 2018), whereas ECMWF’s BL height may include the residual layer or even the entrainment zone above it. In the nighttime plot (Fig. 6, right panel), the retrieved BL tops are as expected significantly lower over land for both Atmos. Chem. Phys., 25, 16491–16510, 2025 https://doi.org/10.5194/acp-25-16491-2025 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact 16499 Figure 6. BL height along the latitude of 16.84° N for September 2012–2022 (Area 2), derived from CALIPSO lidar (blue points) and ECMWF model data (orange points). CALIPSO trajectories were collocated with the nearest ECMWF model grid, and data were averaged over 2° longitudinal intervals. The error bars represent the variability in the BL height. The brown shaded region represents the topography of West Africa, indicating landmass and orographic features influence on the BL structure (sourced from Google Earth). The left figure illustrates results from daytime and the right illustrates nighttime CALIPSO trajectories. datasets. Over the ocean, the agreement between ECMWF and CALIPSO remains good. Over land, however, a different pattern emerges: the ECMWF dataset shows little variability but reports lower BL heights than CALIPSO, particularly further inland (lon >−10°). This again can be explained by the use of thermodynamic criteria to identify the BL top in ECMWF. In contrast, CALIPSO often detects aerosols residing in the residual layer or within the stable nocturnal boundary layer, resulting in systematically higher BL than ECMWF. An additional factor to consider is the quality of the CALIPSO nighttime profiles. The CALIOP instrument has different signal-to-noise characteristics during day and night: while solar background noise degrades daytime profiles, nighttime profiles suffer from lower photon count rates, which makes them noisier, especially over land (Hunt et al., 2009). This effect is consistent with our findings in Appendix A2, where the correlation between ECMWF and CALIPSO is low. The results suggest that the underrepresentation of aerosols in the ECMWF model (Morcrette et al., 2008; Bozzo et al., 2020; Rémy et al., 2024) may also contribute to the observed differences. This effect is relevant both during daytime over land, where dust reduces incoming radiation to the surface and thus influences BL evolution, and during nighttime, when aerosols become trapped in the residual layer and can be detected by CALIPSO. Overall, the two datasets show generally good agreement over the ocean, where both daytime and nighttime results are consistent. This aligns with the findings from Sect. 3.1 (Area 1). The agreement is also stronger during the daytime compared to the nighttime, reflecting the limitations of the satellite nighttime measurements. Over land, however, discrepancies emerge due to the strong diurnal cycle and the different methodologies used to define the BL top. 3.3 Focusing on Cabo Verde and JATAC/ASKOS Cabo Verde is an archipelago in the eastern tropical Atlantic, with distinctive BL dynamics shaped by both the insular geography and the influence of surrounding mountains on airflow patterns. Specifically, the highest point is the Monte Verde (744 m) on the eastern side, but there are also Caixa (535 m) and Madeiral (680 m) on the southern part, as well as Monte Cara (490 m) on the western part. Another geographical characteristic, is that Cabo Verde is situated directly in the path of frequent Saharan dust transport, so the region is often impacted by large dust plumes originating from the African continent and crossing over the islands. These dust events vary significantly in intensity, sometimes accumulating right above the BL or penetrating into it, while at other times showing minimal impact due to lower dust loads. The islands of Cabo Verde, are located nearly 1000 km from the West African coast. The region of São Vicente spans approximately 227 km2, while the neighbouring (northern) island of Santo Antão covers around 785 km2, creating an interface where land and sea effects influence local atmospheric conditions. The origins of air drawn in to the trade winds arriving at Cape Verde are diverse depending on the season; from North America, the Atlantic, Arctic, European and African regions. During autumn, Cape Verde is situated in the direct transport pathway of easterly dust from Africa to the North Atlantic (Carpenter et al., 2010). These sea-air temperature contrasts, rough land surfaces, and fluctuating humidity contribute to a dynamic environment that reflects both marine and coastal BL characteristics. https://doi.org/10.5194/acp-25-16491-2025 Atmos. Chem. Phys., 25, 16491–16510, 2025 16506 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact Figure A2. Distribution of BL heights from CALIPSO (blue) and ECMWF (orange) for daytime (left panel) and nighttime (right panel) conditions. Solid histograms show the frequency of occurrence, with dashed lines indicating the corresponding kernel density estimates (KDEs). The shaded areas mark the overlap between the two datasets: 72 % during daytime and 64 % during nighttime. Figure A3. Scatter plots of BL heights from CALIPSO versus ECMWF for daytime (left panel) and nighttime (right panel). Black points indicate land observations, fitted with the maroon line; white points indicate ocean observations, fitted with the blue line; all data together are fitted with the solid red line. Correlation coefficients (r) are indicated for each subset. Atmos. Chem. Phys., 25, 16491–16510, 2025 https://doi.org/10.5194/acp-25-16491-2025 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact 16507 Figure A4. Correlation between boundary layer heights retrieved from radiosondes and the PollyXT Lidar for June and September 2022 (N=40). The grey dotted line represents the 1 :1 line (y=x), while the grey dashed lines indicate ±150 m from the 1:1 line. The blue dashed line shows the linear regression fit (y=0.91 +4.84) and correlation coefficient r=0.87 highlighting the strong correlation between the two datasets. Data availability. The ASKOS Campaign dataset is available from the ESA Atmospheric Validation Data Centre (EVDC) at https://evdc.esa.int (last access: 19 November 2025) at https://doi.org/10.60621/jatac.campaign.2021.2022.caboverde (Marinou, 2025). Visualized datasets of the ASKOS Campaign and additional information are also available at https://askos.space.noa.gr/data (last access: 19 November 2025); the ERA5 reanalysis dataset used in this study are available from the Copernicus Climate Change Service (C3S) Climate Data Store (https://doi.org/10.24381/cds.adbb2d47, Hersbach et al., 2023). The livas dataset is available upon request. Author contributions. IT and EM conducted the analysis and drafted the manuscript; MT, EG and VA provided methodological guidance and contributed to the interpretation of the data; EM, MT and VA designed the study framework and defined the research objectives; EP, KR and VV contributed to the data curation, processing, visualization and revisions of the results; HB, RE, and AS contributed to the PollyXT data acquisition; HB and ZY assured the Pollynet data processing; EM and VA provided funding acquisition and project administration; all authors edited and reviewed the original draft, provided critical feedback and helped shape the research, analysis and manuscript. Competing interests. The contact author has declared that none of the authors has any competing interests. Disclaimer. Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher. Special issue statement. This article is part of the special issue “The Joint Aeolus Tropical Atlantic Campaign (JATAC) (AMT/ACP inter-journal SI)”. It is not associated with a conference. Acknowledgements. This article is based upon work from COST Action EARLICOST, supported by COST (European Cooperation in Science and Technology). Emmanouil Proestakis acknowledges support from the AXA Research Fund for postdoctoral researchers under the project entitled “Earth Observation for Air-Quality – Dust Fine-Mode (EO4AQ-DustFM)”. Parts of this research have been supported by the German Federal Ministry for Economic Affairs and Energy (BMWi) (grant no. 50EE1721C) and by the German Federal Ministry of Education and Research (BMBF) under the FONA” (grant no. 01LK2001A). Finally, the authors would like to acknowledge the project of ASKOS (grant agreement no. 4000131861/20/NL/IA) from the European Space Agency. The authors would also like to thank the anonymous referees for their reviews, which improved the paper significantly. Financial support. This research has been supported by the Hellenic Foundation for Research and Innovation (project acronym: REVEAL, grant no. 07222), the HORIZON EUROPE Widening participation and spreading excellence (project acronym: PANGEA4CALVAL, grant no. 101079201) funded by the European Union , and the HORIZON EUROPE Climate, Energy and Mobility (project acronym: CERTAINTY, grant no. 101137680). Review statement. This paper was edited by Michael Tjernström and reviewed by three anonymous referees. https://doi.org/10.5194/acp-25-16491-2025 Atmos. Chem. Phys., 25, 16491–16510, 2025 16508 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact References Amiridis, V., Marinou, E., Tsekeri, A., Wandinger, U., Schwarz, A., Giannakaki, E., Mamouri, R., Kokkalis, P., Binietoglou, I., Solomos, S., Herekakis, T., Kazadzis, S., Gerasopoulos, E., Proestakis, E., Kottas, M., Balis, D., Papayannis, A., Kontoes, C., Kourtidis, K., Papagiannopoulos, N., Mona, L., Pappalardo, G., Le Rille, O., and Ansmann, A.: LIVAS: a 3-D multi-wavelength aerosol/cloud database based on CALIPSO and EARLINET, Atmos. Chem. Phys., 15, 7127–7153, https://doi.org/10.5194/acp15-7127-2015, 2015. Ansmann, A., Rittmeister, F., Engelmann, R., Basart, S., Jorba, O., Spyrou, C., Remy, S., Skupin, A., Baars, H., Seifert, P., Senf, F., and Kanitz, T.: Profiling of Saharan dust from the Caribbean to western Africa – Part 2: Shipborne lidar measurements versus forecasts, Atmos. Chem. Phys., 17, 14987–15006, https://doi.org/10.5194/acp-17-14987-2017, 2017. Atlas, D., Walter, B., Chou, S.-H., and Sheu, P.: The structure of the unstable marine boundary layer viewed by lidar and aircraft observations, J. Atmos. Sci., 43, 1301–1318, 1986. Baars, H., Ansmann, A., Engelmann, R., and Althausen, D.: Continuous monitoring of the boundary-layer top with lidar, Atmos. Chem. Phys., 8, 7281–7296, https://doi.org/10.5194/acp-8-72812008, 2008. Bozzo, A., Benedetti, A., Flemming, J., Kipling, Z., and Rémy, S.: An aerosol climatology for global models based on the tropospheric aerosol scheme in the Integrated Forecasting System of ECMWF, Geosci. Model Dev., 13, 1007–1034, https://doi.org/10.5194/gmd-13-1007-2020, 2020. Brooks, I. M.: Finding boundary layer top: Application of a wavelet covariance transform to lidar backscatter profiles, J. Atmos. Ocean. Tech., 20, 1092–1105, 2003. Brooks, I. M., Tjernström, M., Persson, P. O. G., Shupe, M. D., Atkinson, R. A., Canut, G., Birch, C. E., Mauritsen, T., Sedlar, J., and Brooks, B. J.: The turbulent structure of the Arctic summer boundary layer during the Arctic summer cloud-ocean study, J. Geophys. Res.-Atmos., 122, 9685–9704, 2017. Carlson, T. N. and Prospero, J. M.: The large-scale movement of Saharan air outbreaks over the northern equatorial Atlantic, J. Appl. Meteorol. Clim., 11, 283–297, 1972. Carpenter, L. J., Fleming, Z. L., Read, K. A., Lee, J. D., Møller, S. J., Hopkins, J. R., Purvis, R. M., Lewis, A. C., Müller, K., Heinold, B., Herrmann, H., Fomba, K. W., van Pinxteren, D., Müller, C., Tegen, I., Wiedensohler, A., Müller, T., Niedermeier, N., Achterberg, E. P., Patey, M. D., Kozlova, E. A., Heimann, M., Heard, D. E., Plane, J. M. C., Mahajan, A., Oetjen, H., Ingham, T., Stone, D., Whalley, L. K., Evans, M. J., Pilling, M. J., Leigh, R. J., Monks, P. S., Karunaharan, A., Vaughan, S., Arnold, S. R., Tschritter, J., Pöhler, D., Frieß, U., Holla, R., Mendes, L. M., Lopez, H., Faria, B., Manning, A. J., and Wallace, D. W. R.: Seasonal characteristics of tropical marine boundary layer air measured at the Cape Verde Atmospheric Observatory, J. Atmos. Chem., 67, 87–140, 2010. Croft, B., Martin, R. V., Moore, R. H., Ziemba, L. D., Crosbie, E. C., Liu, H., Russell, L. M., Saliba, G., Wisthaler, A., Müller, M., Schiller, A., Galí, M., Chang, R. Y.-W., McDuffie, E. E., Bilsback, K. R., and Pierce, J. R.: Factors controlling marine aerosol size distributions and their climate effects over the northwest Atlantic Ocean region, Atmos. Chem. Phys., 21, 1889–1916, https://doi.org/10.5194/acp-21-1889-2021, 2021. Dunion, J. P. and Velden, C. S.: The impact of the Saharan air layer on Atlantic tropical cyclone activity, B. Am. Meteorol. Soc., 85, 353–366, 2004. ECMWF: IFS Documentation CY43R3 – Part IV: Physical processes, 4, ECMWF, https://doi.org/10.21957/efyk72kl, 2017. Engelmann, R., Kanitz, T., Baars, H., Heese, B., Althausen, D., Skupin, A., Wandinger, U., Komppula, M., Stachlewska, I. S., Amiridis, V., Marinou, E., Mattis, I., Linn´ r, H., and Ansmann, A.: The automated multiwavelength Raman polarization and watervapor lidar PollyXT: the neXT generation, Atmos. Meas. Tech., 9, 1767–1784, https://doi.org/10.5194/amt-9-1767-2016, 2016. Evan, A. T., Vimont, D. J., Heidinger, A. K., Kossin, J. P., and Bennartz, R.: The role of aerosols in the evolution of tropical North Atlantic Ocean temperature anomalies, Science, 324, 778–781, 2009. Flamant, C., Pelon, J., Flamant, P. H., and Durand, P.: Lidar determination of the entrainment zone thickness at the top of the unstable marine atmospheric boundary layer, Bound.-Lay. Meteorol., 83, 247–284, 1997. Foltz, G. R. and McPhaden, M. J.: Impact of Saharan dust on tropical North Atlantic SST, J. Climate, 21, 5048–5060, 2008. Garcia-Carreras, L., Parker, D., Marsham, J., Rosenberg, P., Brooks, I., Lock, A., Marenco, F., McQuaid, J., and Hobby, M.: The turbulent structure and diurnal growth of the Saharan atmospheric boundary layer, J. Atmos. Sci., 72, 693–713, 2015. Garratt, J. R.: The atmospheric boundary layer, Earth-Sci. Rev., 37, 89–134, 1994. Giménez, J., Pastor, C., Castañer, R., Nicolás, J., Crespo, J., and Carratalá, A.: Influence of Saharan dust outbreaks and atmospheric stability upon vertical profiles of size-segregated aerosols and water vapor, Atmos. Environ., 44, 338–346, 2010. Hanson, H. P.: Marine stratocumulus climatologies, Int. J. Climatol., 11, 147–164, 1991. Hersbach, H., Bell, B., Berrisford, P., Biavati, G., Horányi, A., Muñoz Sabater, J., Nicolas, J., Peubey, C., Radu, R., Rozum, I., Schepers, D., Simmons, A., Soci, C., Dee, D., and Thépaut, J.-N.: ERA5 hourly data on single levels from 1940 to present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [data set], https://doi.org/10.24381/cds.adbb2d47, 2023. Hunt, W. H., Winker, D. M., Vaughan, M. A., Powell, K. A., Lucker, P. L., and Weimer, C.: CALIPSO lidar description and performance assessment, J. Atmos. Ocean. Tech., 26, 1214–1228, 2009. Kallos, G., Astitha, M., Katsafados, P., and Spyrou, C.: Long-range transport of anthropogenically and naturally produced particulate matter in the Mediterranean and North Atlantic: Current state of knowledge, J. Appl. Meteorol. Clim., 46, 1230–1251, 2007. Li, H., Liu, B., Ma, X., Jin, S., Ma, Y., Zhao, Y., and Gong, W.: Evaluation of retrieval methods for planetary boundary layer height based on radiosonde data, Atmos. Meas. Tech., 14, 5977–5986, https://doi.org/10.5194/amt-14-5977-2021, 2021. Li, Z., Guo, J., Ding, A., Liao, H., Liu, J., Sun, Y., Wang, T., Xue, H., Zhang, H., and Zhu, B.: Aerosol and boundary-layer interactions and impact on air quality, Natl. Sci. Rev., 4, 810–833, 2017. Liu, B., Ma, Y., Liu, J., Gong, W., Wang, W., and Zhang, M.: Graphics algorithm for deriving atmospheric boundary layer Atmos. Chem. Phys., 25, 16491–16510, 2025 https://doi.org/10.5194/acp-25-16491-2025 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact 16509 heights from CALIPSO data, Atmos. Meas. Tech., 11, 5075– 5085, https://doi.org/10.5194/amt-11-5075-2018, 2018. Liu, J.-W., Yang, L., and Xie, S.-P.: CALIPSO Observed low clouds over the subtropical northeast Pacific: Summer climatology and interannual variability, J. Geophys. Res.-Atmos., 129, e2023JD039168, https://doi.org/10.1029/2023JD039168, 2024. Liu, Z., Vaughan, M., Winker, D., Kittaka, C., Getzewich, B., Kuehn, R., Omar, A., Powell, K., Trepte, C., and Hostetler, C.: The CALIPSO lidar cloud and aerosol discrimination: Version 2 algorithm and initial assessment of performance, J. Atmos. Ocean. Tech., 26, 1198–1213, 2009. Luo, B., Minnett, P. J., Zuidema, P., Nalli, N. R., and Akella, S.: Saharan Dust Effects on North Atlantic Sea-Surface Skin Temperatures, J. Geophys. Res.-Oceans, 126, e2021JC017282, https://doi.org/10.1029/2021JC017282, 2021. Marinou, E., Amiridis, V., Binietoglou, I., Tsikerdekis, A., Solomos, S., Proestakis, E., Konsta, D., Papagiannopoulos, N., Tsekeri, A., Vlastou, G., Zanis, P., Balis, D., Wandinger, U., and Ansmann, A.: Three-dimensional evolution of Saharan dust transport towards Europe based on a 9-year EARLINEToptimized CALIPSO dataset, Atmos. Chem. Phys., 17, 5893– 5919, https://doi.org/10.5194/acp-17-5893-2017, 2017. Marinou, E., Voudouri, K. A., Tsikoudi, I., Drakaki, E., Tsekeri, A., Rosoldi, M., Ene, D., Baars, H., O’Connor, E., Amiridis, V., and Meleti, C.: Geometrical and microphysical properties of clouds formed in the presence of dust above the Eastern Mediterranean, Remote Sens., 13, 5001, https://doi.org/10.3390/rs13245001, 2021. Marinou, E., Paschou, P., Tsikoudi, I., Tsekeri, A., Daskalopoulou, V., Kouklaki, D., Siomos, N., Spanakis-Misirlis, V., Voudouri, K. A., Georgiou, T., Drakaki, E., Kampouri, A., Papachristopoulou, K., Mavropoulou, I., Mallios, S., Proestakis, E., Gkikas, A., Koutsoupi, I., Raptis, I. P., Kazadzis, S., Baars, H., Floutsi, A., Pirloaga, R., Nemuc, A., Marenco, F., Kezoudi, M., Papetta, A., Moˇ cnik, G., Díez, J. Y., Ryder, C. L., Ratcliffe, N., Kandler, K., Sudharaj, A., and Amiridis, V.: An overview of the ASKOS campaign in Cabo Verde, Environ. Sci. Proc., 26, 200, https://doi.org/10.3390/environsciproc2023026200, 2023. Marinou, E.: ASKOS Campaign Dataset (JATAC 2021–2022, Cabo Verde), ESA Atmospheric Validation Data Centre [data set], https://doi.org/10.60621/jatac.campaign.2021.2022.caboverde, 2023. Marsham, J. H., Parker, D. J., Grams, C. M., Johnson, B. T., Grey, W. M. F., and Ross, A. N.: Observations of mesoscale and boundary-layer scale circulations affecting dust transport and uplift over the Sahara, Atmos. Chem. Phys., 8, 6979–6993, https://doi.org/10.5194/acp-8-6979-2008, 2008. Menut, L., Masson, O., and Bessagnet, B.: Contribution of Saharan dust on radionuclide aerosol activity levels in Europe? The 21– 22 February 2004 case study, J. Geophys. Res.-Atmos., 114, 10 pp., https://doi.org/10.1029/2009JD011767, 2009. Morcrette, J.-J., Beljaars, A., Benedetti, A., Jones, L., and Boucher, O.: Sea-salt and dust aerosols in the ECMWF IFS model, Geophys. Res. Lett., 35, 5 pp., https://doi.org/10.1029/2008GL036041, 2008. O’Connor, E. J., Illingworth, A. J., Brooks, I. M., Westbrook, C. D., Hogan, R. J., Davies, F., and Brooks, B. J.: A method for estimating the turbulent kinetic energy dissipation rate from a vertically pointing Doppler lidar, and indep. Tech., 27, 1652–1664, 2010. Pearson, G., Davies, F., and Collier, C.: An analysis of the performance of the UFAM pulsed Doppler lidar for observing the boundary layer, J. Atmos. Ocean. Tech., 26, 240–250, 2009. Pena, A., Gryning, S.-E., and Floors, R. R.: Lidar observations of marine boundary-layer winds and heights: a preliminary study, Meteorol. Z., 24, 581–589, 2015. Pérez, C., Nickovic, S., Baldasano, J., Sicard, M., Rocadenbosch, F., and Cachorro, V.: A long Saharan dust event over the western Mediterranean: Lidar, Sun photometer observations, and regional dust modeling, J. Geophys. Res.-Atmos., 111, 16 pp., https://doi.org/10.1029/2005JD006579, 2006. Proestakis, E., Gkikas, A., Georgiou, T., Kampouri, A., Drakaki, E., Ryder, C. L., Marenco, F., Marinou, E., and Amiridis, V.: A near-global multiyear climate data record of the fine-mode and coarse-mode components of atmospheric pure dust, Atmos. Meas. Tech., 17, 3625–3667, https://doi.org/10.5194/amt17-3625-2024, 2024. Prospero, J. M. and Mayol-Bracero, O. L.: Understanding the transport and impact of African dust on the Caribbean Basin, B. Am. Meteorol. Soc., 94, 1329–1337, 2013. Rémy, S., Kipling, Z., Flemming, J., Boucher, O., Nabat, P., Michou, M., Bozzo, A., Ades, M., Huijnen, V., Benedetti, A., Engelen, R., Peuch, V.-H., and Morcrette, J.-J.: Description and evaluation of the tropospheric aerosol scheme in the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS-AER, cycle 45R1), Geosci. Model Dev., 12, 4627–4659, https://doi.org/10.5194/gmd-124627-2019, 2019. Rémy, S., Kipling, Z., Huijnen, V., Flemming, J., Nabat, P., Michou, M., Ades, M., Engelen, R., and Peuch, V.-H.: Description and evaluation of the tropospheric aerosol scheme in the Integrated Forecasting System (IFS-AER, cycle 47R1) of ECMWF, Geosci. Model Dev., 15, 4881–4912, https://doi.org/10.5194/gmd-154881-2022, 2022. Rémy, S., Metzger, S., Huijnen, V., Williams, J. E., and Flemming, J.: An improved representation of aerosol in the ECMWF IFS-COMPO 49R1 through the integration of EQSAM4Climv12 – a first attempt at simulating aerosol acidity, Geosci. Model Dev., 17, 7539–7567, https://doi.org/10.5194/gmd-177539-2024, 2024. Rolph, G., Stein, A., and Stunder, B.: Real-time environmental applications and display system: READY, Environ. Model. Softw., 95, 210–228, 2017. Schnitker, D.: Climatic variability and deep ocean circulation: evidence from the North Atlantic, Palaeogeogr. Palaeocl. Palaeoecol., 40, 213–234, 1982. Seibert, P., Beyrich, F., Gryning, S.-E., Joffre, S., Rasmussen, A., and Tercier, P.: Review and intercomparison of operational methods for the determination of the mixing height, Atmos. Environ., 34, 1001–1027, 2000. Seidel, D. J., Ao, C. O., and Li, K.: Estimating climatological planetary boundary layer heights from radiosonde observations: Comparison of methods and uncertainty analysis, J. Geophys. Res.- Atmos., 115, 195 pp., https://doi.org/10.1029/2009JD013680, 2010. Seidel, D. J., Zhang, Y., Beljaars, A., Golaz, J.-C., Jacobson, A. R., and Medeiros, B.: Climatology of the planetary boundary layer over the continental United States and Europe, J. Geophys. Res.- https://doi.org/10.5194/acp-25-16491-2025 Atmos. Chem. Phys., 25, 16491–16510, 2025 16510 I. Tsikoudi et al.: Atmospheric boundary layer in the Atlantic: the desert dust impact Atmos., 117, 15 pp., https://doi.org/10.1029/2012JD018143, 2012. Stephens, G., Winker, D., Pelon, J., Trepte, C., Vane, D., Yuhas, C., L’ecuyer, T., and Lebsock, M.: CloudSat and CALIPSO within the A-Train: Ten years of actively observing the Earth system, B. Am. Meteorol. Soc., 99, 569–581, 2018. Stull, R. B.: Mean boundary layer characteristics, in: An introduction to boundary layer meteorology, Springer, 1–27, https://doi.org/10.1007/978-94-009-3027-8_1, 1988. Sun, Y. and Zhao, C.: Influence of Saharan dust on the large-scale meteorological environment for development of tropical cyclone over North Atlantic Ocean Basin, J. Geophys. Res.-Atmos., 125, e2020JD033454, https://doi.org/10.1029/2020JD033454, 2020. Tackett, J. L., Winker, D. M., Getzewich, B. J., Vaughan, M. A., Young, S. A., and Kar, J.: CALIPSO lidar level 3 aerosol profile product: version 3 algorithm design, Atmos. Meas. Tech., 11, 4129–4152, https://doi.org/10.5194/amt-11-4129-2018, 2018. Teixeira, J., Piepmeier, J. R., Nehrir, A. R., Ao, C. O., Chen, S. S., Clayson, C. A., Fridlind, A. M., Lebsock, M., McCarty, W., Salmun, H., Santanello, J. A., Turner, D. D., Wang, Z., and Zeng, X.: Toward a Global Planetary Boundary Layer Observing System: A Summary, B. Am. Meteorol. Soc., E1566–E1579, https://doi.org/10.1175/BAMS-D-23-0228.1, 2025. Tombrou, M., Dandou, A., Helmis, C., Akylas, E., Angelopoulos, G., Flocas, H., Assimakopoulos, V., and Soulakellis, N.: Model evaluation of the atmospheric boundary layer and mixed-layer evolution, Bound.-Lay. Meteorol., 124, 61–79, 2007. Tombrou, M., Bossioli, E., Kalogiros, J., Allan, J., Bacak, A., Biskos, G., Coe, H., Dandou, A., Kouvarakis, G., Mihalopoulos, N., Percival, C., Protonotariou, A., and Szabó-Takács, B.: Physical and chemical processes of air masses in the Aegean Sea during Etesians: Aegean-GAME airborne campaign, Sci. Total Environ., 506, 201–216, 2015. Tsaknakis, G., Papayannis, A., Kokkalis, P., Amiridis, V., Kambezidis, H. D., Mamouri, R. E., Georgoussis, G., and Avdikos, G.: Inter-comparison of lidar and ceilometer retrievals for aerosol and Planetary Boundary Layer profiling over Athens, Greece, Atmos. Meas. Tech., 4, 1261–1273, https://doi.org/10.5194/amt-41261-2011, 2011. Tsikoudi, I., Marinou, E., Vakkari, V., Gialitaki, A., Tsichla, M., Amiridis, V., Komppula, M., Raptis, I. P., Panagiotis, I., Kampouri, A., Daskalopoulou, V., Mihalopoulos, N., Giannakaki, E., Tombrou, M., and Flocas, H.: PBL Height Retrievals at a Coastal Site Using Multi-Instrument Profiling Methods, Remote Sens., 14, 4057, https://doi.org/10.3390/rs14164057, 2022. Tsikoudi, I., Marinou, E., Voudouri, K., Koutsoupi, I., Drakaki, E., Kampouri, A., Vakkari, V., Baars, H., Giannakaki, E., Tombrou, M., and Amiridis, V.: PBL Height Retrievals during ASKOS Campaign, Environ. Sci. Proc., 26, 23, https://doi.org/10.3390/environsciproc2023026023, 2023. Vakkari, V., O’Connor, E. J., Nisantzi, A., Mamouri, R. E., and Hadjimitsis, D. G.: Low-level mixing height detection in coastal locations with a scanning Doppler lidar, Atmos. Meas. Tech., 8, 1875–1885, https://doi.org/10.5194/amt-8-1875-2015, 2015. Vakkari, V., Manninen, A. J., O’Connor, E. J., Schween, J. H., van Zyl, P. G., and Marinou, E.: A novel post-processing algorithm for Halo Doppler lidars, Atmos. Meas. Tech., 12, 839–852, https://doi.org/10.5194/amt-12-839-2019, 2019. Vaughan, M. A., Powell, K. A., Winker, D. M., Hostetler, C. A., Kuehn, R. E., Hunt, W. H., Getzewich, B. J., Young, S. A., Liu, Z., and McGill, M. J.: Fully automated detection of cloud and aerosol layers in the CALIPSO lidar measurements, J. Atmos. Ocean. Tech., 26, 2034–2050, 2009. Vogelezang, D. and Holtslag, A.: Evaluation and model impacts of alternative boundary-layer height formulations, Bound.-Lay. Meteorol., 81, 245–269, 1996. Wiegner, M., Emeis, S., Freudenthaler, V., Heese, B., Junkermann, W., Münkel, C., Schäfer, K., Seefeldner, M., and Vogt, S.: Mixing layer height over Munich, Germany: Variability and comparisons of different methodologies, J. Geophys. Res.-Atmos., 111, 17 pp., https://doi.org/10.1029/2005JD006593, 2006. Winker, D. M., Pelon, J., Coakley, J. A., Ackerman, S. A., Charlson, R. J., Colarco, P. R., Flamant, P., Fu, Q., Hoff, R. M., Kittaka, C., Kubar, T. L., Treut, H. L., Mccormick, M. P., Mégie, G., Poole, L., Powell, K., Trepte, C., Vaughan, M. A., and Wielicki, B. A.: The CALIPSO mission: A global 3D view of aerosols and clouds, B. Am. Meteorol. Soc., 91, 1211–1230, 2010. Winker, D. M., Tackett, J. L., Getzewich, B. J., Liu, Z., Vaughan, M. A., and Rogers, R. R.: The global 3-D distribution of tropospheric aerosols as characterized by CALIOP, Atmos. Chem. Phys., 13, 3345–3361, https://doi.org/10.5194/acp-133345-2013, 2013. Wood, R.: Stratocumulus clouds, Mon. Weather Rev., 140, 2373– 2423, 2012. Wu, L.: Impact of Saharan air layer on hurricane peak intensity, Geophys. Res. Lett., 34, 5 pp., https://doi.org/10.1029/2007GL029564, 2007. Yu, H., Kaufman, Y. J., Chin, M., Feingold, G., Remer, L. A., Anderson, T. L., Balkanski, Y., Bellouin, N., Boucher, O., Christopher, S., DeCola, P., Kahn, R., Koch, D., Loeb, N., Reddy, M. S., Schulz, M., Takemura, T., and Zhou, M.: A review of measurement-based assessments of the aerosol direct radiative effect and forcing, Atmos. Chem. Phys., 6, 613–666, https://doi.org/10.5194/acp-6-613-2006, 2006. Zeng, S., Vaughan, M., Liu, Z., Trepte, C., Kar, J., Omar, A., Winker, D., Lucker, P., Hu, Y., Getzewich, B., and Avery, M.: Application of high-dimensional fuzzy k-means cluster analysis to CALIOP/CALIPSO version 4.1 cloud– aerosol discrimination, Atmos. Meas. Tech., 12, 2261–2285, https://doi.org/10.5194/amt-12-2261-2019, 2019. Zhou, W., Leung, L. R., Lu, J., Yang, D., and Song, F.: Contrasting recent and future ITCZ changes from distinct tropical warming patterns, Geophys. Res. Lett., 47, e2020GL089846, https://doi.org/10.1029/2020GL089846, 2020. Atmos. Chem. Phys., 25, 16491–16510, 2025 https://doi.org/10.5194/acp-25-16491-2025