Abstract

The lack of information on the vertical distribution of dust, in turn, results in large uncertainties when attempting to evaluate the impacts of dust on climate processes. We analyzed over two years of LIDAR measurements from NASA’s CALIPSO and CloudSat satellites to document the vertical pathways of transatlantic transport of Saharan dust. Our analysis overcomes the limitations of quantitative dust detections with passive satellite measurements over land and low clouds and provides the fine vertical resolved structures. The results show the strong seasonal shift in dust source regions and transportation pathways due to the meteorological and thermodynamical conditions, which also control the dust vertical distribution as well as the depth of the dust layer. The dust layer top descending rates of 35 m/degree in summer, 25 m/degree in autumn and spring, and 10 m/degree in winter are found, respectively, while the dust is being transported across the Atlantic. Comparison with the model simulation highlights the potentials of dust observations using CALIPSO LIDAR. The observed seasonal dependence of these pathways gives new insights into the transport of the Saharan dust and provides important guidance for simulations of the production and transport of the global dust aerosol.

1. Introduction

Mineral dust is mainly generated from the continents and can be transported both horizontally and vertically into the atmosphere and affect vast area on the earth. Mineral dust plays an important role in climate processes by affecting radiation [1, 2], modifying cloud and precipitation processes by acting as cloud condensation nuclei (CCN), giant CCN (GCCN), and effective ice nuclei (IN) [38]. Also mineral is an important source of nutrition which modulates a variety of marine biogeochemical processes [9]. However, among different types of aerosols, there are large uncertainties in dust sources, sinks, and their long-range transports [10, 11].

The Sahara is a primary source of dust, the majority of which is transported westward across the Atlantic to the Southeastern United States, the Caribbean Sea, and South America [1214]. African dust and its transport have been documented by several passive remote sensors, for example, the Advanced Vary High Resolution Radiometer (AVHRR), Meteosat, Total Ozone Mapping Spectrometer (TOMS), and Moderate Resolution Imaging Spectroradiometer (MODIS) [1518]. Modeling and observation results also suggest that Saharan dust and dust containing Sahara Air Layer (SAL) could impact the hurricane formation and evolution in the region of the North Atlantic [19, 20]. Although our knowledge of the properties of Saharan dust, its transport, and impact has improved dramatically over the last three decades, the transatlantic transport of dust remains poorly understood, particularly its vertical distribution. For example, among AeroCom models (http://nansen.ipsl.jussieu.fr/AEROCOM), global and annual average dust mass fractions of total mass above 5 km altitude range from less 5% to close 30% [10]. Our limited understanding results from the inherent inability of passive sensors to provide vertically resolved measurements over large scales. The lack of information on the vertical distribution of dust, in turn, results in large uncertainties when attempting to evaluate the impacts of dust on climate processes [1, 4].

Multiyear polarization LIDAR measurements from CALIPSO satellite provide new insights into the vertical distribution of dust on a global scale [14, 21]. An approach was developed to provide vertical extinction and mass mixing ratio distribution based on CALIPSO LIDAR measurements. With these data we show a 3D view of across Atlantic dust transport and discuss their associations with large-scale dynamics. We also compare with model simulations to illustrate major issues in simulated dust vertical distributions. Section 2 provides details of data analyses. Results and discussions are presented in Section 3. Finally, a brief summary is given in Section 4.

2. Dataset and Methodology

The CALIPSO and CloudSat satellites were launched together in April 2006 to join in the A-Train constellation and have been acquiring data since June 2006. The key instrument in the CALIPSO satellite payload is the Cloud-Aerosol LIDAR with Orthogonal Polarization (CALIOP) [21]. CALIOP provides high vertical resolution profiles of aerosols and clouds (especially optically thin clouds) globally. The vertical and horizontal resolutions of the CALIOP are 30 and 333 m below 8.2 km, and 60 and 1000 m between 8.2 and 20.2 km, respectively. CloudSat carried the first space borne cloud radar—a 94.05 GHz Cloud Profiling Radar (CPR) with sensitivity of −30 dBZ [22]. The CPR measurements provide radar reflectivity (Ze) profiles with a vertical resolution of 240 m and a footprint of 1.4 km × 1.8 km (cross and along track). Due to their different wavelengths, CPR and CALIOP complement with each other well to provide global three-dimensional aerosol and cloud structure [23].

Taking the advantage of the synergy of the CALIPSO and CloudSat satellites, the CALIOP L1B profiles are collocated to CloudSat radar footprints based on their surface location. A section of collocated CPR and CALIOP measurements with dust layer and clouds is presented in Figures 1(a) and 1(b). Although CALIOP measurement itself can be used to separate aerosol and clouds, optically thick aerosol layers, such as dust layers over the dust source regions and smoke layers, are often misclassified as clouds [14, 24, 25]. However, these geometrically and optically thick aerosol layers have Ze values much lower than the CPR detection limit due to their small particle sizes. For clouds with a similar geometrical thickness and cloud temperature, the layer maximum Ze is often over −10 dBZ [26]. With the additional Ze profile, we can reduce the occurrence of this misclassification. Figure 1(c) shows identified aerosol layer by removing cloudy profiles (indicated with black lines) with CALIOP and CPR measurements. Optically thick high and middle clouds can significantly attenuate CALIOP signals. Therefore, it is challenging to reliably identify aerosol layers below these clouds. In this study, aerosol layers indentified under cloud-free conditions, in conjunction with optically thin cirrus, or above low clouds (cloud top height 3 km, above ground level) are used for further analyses.

Dust aerosols have higher linear depolarization ratios (the ratio of perpendicular to parallel polarization components measured by a LIDAR) than the other types of aerosols due to dust nonspherical shape, which can be used to distinguish them from other types of aerosols [14, 2729]. In this study, dust layers are first identified with a threshold method based on 1-km layer-averaged volume depolarization ratio (VDR), defined as the ratio of layer integrated perpendicular to parallel components at 532 nm based on the horizontal 1.5-km averaged level 1 B data. As discussed in detail by Liu et al. [14], is used as a threshold for identifying dust layers and has been applied to the nighttime CALIPSO measurements during June 2006 through August 2008 to derive the vertical distribution of dust layers. With 1-km layer dust detection results, dust layer base and top heights are further refined [14]. The same definition of dust occurrence as Liu et al. [14] is used in this study.

Figure 1(d) gives the retrieved dust extinction coefficients based on Fernald backward integration method [30] combined with our proposed constrained-iteration scheme. Fernald method is often used to retrieve aerosol extinction profile which has both backward and forward integration algorithms. The backward integration is often adopted by the ground-based LIDAR while forward integration method for the space-borne one. For the atmospheric aerosol with the larger optical depth, the Fernald forward integration can produce larger errors because of the divergence. But here we apply backward integration algorithm to the CALIOP profile and constrain the retrieval iteration taking the advantage of the calibration of CALIOP using the averaged molecular backscatter between 30 and 34 km altitude on 532 nm wavelength [21]. Here is our proposed scheme: first, given an arbitrary extinction coefficient to the dust base, using the Fernald backward integration method, one can retrieve the extinction coefficients at 532 nm wavelength above the dust base, and the optical depth above the dust base can be obtained, recorded as ODBW. Second, since CALIOP is well calibrated at this wavelength, the difference between the initial arbitrary extinction coefficient and the attenuated backscatter coefficient at the dust base can also be used to calculate the attenuation above the dust base thus obtaining the optical depth, recorded as ODAT. Third, one can adjust the initial extinction coefficient at the dust base to start the iteration till these two optical depths are equal (in practice, assuming the absolute difference between ODBW and ODAT is smaller than 0.05). The LIDAR ratio of the dust aerosol uses the typical value of 40 Sr suggested by CALIOP [29].

A comparison of the retrieved dust optical depths with a method developed by Hu et al. [31] with CALIOP ocean surface returns is presented in Figure 2. It clearly shows [that dust optical depths from the approach we proposed are consistent Hu et al. with a mean bias about 0.05.] Considering uncertainties in Hu et al. due to errors in estimated near surface wind speeds used to estimate two-way transmittance and in neglecting dust in the marine boundary layer, the agreement between two is very good and we can use extinction coefficients confidently for further analyses.

One application of dust extinction coefficient is to estimate dust mass concentration, which provides a more direct way to compare with model dust simulation. The specific extinction coefficient (the ratio of the extinction coefficient to mass) is used to convert extinction coefficient to mass concentration of dust which has been reported in many literatures. Maring et al. [32] derived 0.5 m2/g near Africa and 0.6 m2/g in Barbados and showed a small dust size distribution change when transported across the Atlantic. Chin et al. [33] used 0.8 m2/g for dust in the model simulations. Moulin et al. [34] reviewed dust models giving 0.6–0.8 m2/g. Based on these published results and calculations with recent in situ dust size distributions, we decided to use 0.7 m2/g in this study with an uncertainty about 30%.

3. Results and Discussions

3.1. Horizontal Dust Distributions

The height resolved seasonal dust occurrence maps (in grid box) over 100°W to 30°E and 15°S to 45°N are presented in Figure 3, which clearly show seasonal variations in the horizontal pathways of cross-Atlantic African dust transport. CALIPSO dust samples over Africa are mainly observed in cloud-free conditions or in conjunction with optically thin cirrus, while dust samples over the ocean are observed primarily in conjunction with low clouds. Comparing results without considering dust over low clouds (see Figure 2 in Liu et al. [14]), significant dust aerosols are transported across the Atlantic above low-level clouds. Detection of dust layers above low clouds is extremely difficult for passive sensors. Thus, CALIOP measurements provide the most reliable dust detection over regions with high occurrences of low clouds [26].

Figure 3 clearly shows that the transatlantic African dust occurs all year long, but June-July-August (JJA) is the most active season. During JJA, over a large area of the Northern African continental has dust occurrence over 80%, and the majority of them have dust top heights reach above 4 km. Dense dust layers during JJA can reach above 6 km. March-April-May (MAM) is the second most active dust season with high occurrence areas concentrated over a smaller region than JJA. Compared with JJA, there are significant less dense dust layers reaching above 4 km during MAM. December-January-February (DJF) is the quietest period with a small high occurrence area close to equator with most dust layer concentrated below 2 km overland.

The high dust occurrence regions over the Atlantic indicate the dust transport pathway. During JJA, dust aerosols are transported the furthest west. Over the most area of the west Atlantic between 10–30°N, there is over 40% chance to find dust layer below 4 km. Compared with JJA and September-October-November (SON), dust-transporting pathways shift southwards 10° in DJF and MAM. During DJF, only limited dust storms are transported across the Atlantic, and dust aerosols are transported southward after reaching the west coast of the Atlantic while they are transported northward in the other seasons. It is clear that the vertical ranges of dust transported across the Atlantic depend on the dust vertical distributions of dust aerosol in the source regions. With the zonal dust cross-sections, this feature will be discussed further in the next subsection.

3.2. Vertical Dust Distributions

To illustrate the vertical pathway of African dust across the Atlantic Ocean, all dusty profiles are used to calculate the zonal cross-sections of vertical dust distributions as illustrated in Figure 3. The mean vertical dust occurrence is calculated as the ratio of total dust cases observed at a given height range to the total dust case over the region within a 2-degree longitude bin based on profiles containing dust layers below 8 km. The vertical resolution of the calculation is 60 m. Thus, a 100% occurrence in Figure 4(a) indicates that dust always occurs in this height range when a dust layer is transported into the region.

The vertical distribution of dust over source regions (east of −18° longitude) shows strong seasonal variations. During JJA, most dust layers have tops above 5 km, while the majority of dust layers lie below 4 km during other seasons. Except MAM, mean dust tops generally increase westward but with noticeably different slopes. As these layers are transported westward over the eastern edge of the Atlantic Ocean, their bases are gradually lifted due to strong inversions near the top of the marine boundary layer. The corresponding longitude of the maximum dust layer base height is about 25°W, 30°W, and 35°W in JJA, SON, and MAM, respectively. The net effect of the changes in layer base and top heights is to confine most dust layers to altitudes between 2 and 4 km as they are transported further west. Although most dust layers are lifted above the marine boundary layer, during winter a greater proportion of dust is transported into the marine boundary layer near the west coast of Africa.

The fates of dust aerosols are determined by dry deposition, wet removal, cloud microphysical processes, and dynamics [35, 36]. The mean dust layer top slowly descends across the Atlantic at rates of about 35 m/degree in JJA, 25 m/degree in SON and MAM, and 10 m/deg in DJF, mainly due to dry deposition and large-scale subsidence. Dust layers approaching the Caribbean Sea and the South American coast could either be lifted to higher altitudes or brought down into the boundary layer as a result of more vigorous convection. LIDAR data from the first two and half years shows large interannual variability in the vertical distribution of dust over these regions. Dust particles lifted to high altitudes and colder temperatures can act effective ice nuclei, thereby affecting cloud properties and precipitation generation, while those entrained into the boundary layer are more efficiently removed through deposition to the surface, attachment to cloud drops, or coalescence scavenging.

Along the pathway, the variations of vertical dust distributions reflect the dust layer thickness changes. Over the ocean, the dust occurrence magnitudes decrease significantly as dust is transported further west. This reflects that the dust layer thicknesses decrease along the transportation pathway, which could be caused by the horizontal SAL expansion or mixing across SAL vertical boundaries. Although dust layer thicknesses between JJA and MAM are similar over the land, thick dust layers are transported further west in JJA than in MAM.

Dust aerosol extinction coefficient and particle depolarization ratio show dust loading and particle shape variations along the transportation. Due to large dust particle size, high dust extinction (532 nm) could indicate regions with major dust storms. It clearly shows that these regions with high dust extinctions are not exactly correlated with high dust occurrence regions. Therefore, everyday boundary layer dust aerosols contribute significantly to the high occurrence of dust aerosols overland. This is consistent with year-round dusty boundary layer observed in West Africa with sophisticated ground-based instrumentation [37]. Overland, DJF has the highest mean extinction coefficient, although the dust layers are shallower during DJF than the other seasons. SON has the lowest mean extinction coefficient overland.

The particle depolarization ratios given in Figure 4(c) are controlled by two major factors: dust particle shape and size and mixing of dust with different aerosols. Over land, highest depolarization ratios are typically found near surface, which could be mainly due to large irregular dust particles. Other than DJF, vertical particle depolarization ratios have small variations with height, which indicates that the dust layer is well mixed. During DJF, particle depolarization ratios above 2 km are noticeably smaller than the aerosols below, which could be the results of different mixings of dust with other aerosols or vertical particle size and shape stratifications within the dust layer. The seasonal differences of dust transport through marine boundary layers are illustrated more clearly with dust occurrence and depolarization. The particle depolarization ratios revealed that dusty layers can push a few hundred kilometers west over the ocean before significantly mixed with sea salts (indicated with high extinction coefficients and low particle depolarization ratios) during DJF and MAM. Over ocean, high dust occurrence regions are coincident with high particle depolarization regions as well. Therefore, these regions are good indicators for the positions of unmixed SAL cores when the SALs propagate across the Atlantic. On the other hand, it suggests that strong mixings happen along the dust transportation, which could be very difficult for the model to simulate.

To see how the model captures these general features revealed with CALIOP measurements, dust extinction coefficients are converted to mass mixing ratio for model simulation comparisons. Figure 4(d) shows the dust mass mixing ratio estimated from retrieved extinctions. Although the uncertainties can be as high as 30% mainly due to the selection of specific extinction coefficient, the main height-longitude structure is reliable. It is well known that there are large differences in simulated dust mass distributions among state-of-art models [34, 36, 38, 39]. Thus, extended evaluations of these models are out of the scope of this paper, and we only use one published dust modeling result to highlight the potentials of using CALIOP observations to identify the model simulation weaknesses [36]. Figures 4(d)–4(f) clearly show that the DEAD dust model captures the general seasonal variations of dust vertical and horizontal distributions, but there are several noticeable differences. First, the DEAD model predicts more dust below 2 km than estimations from observations. Second, there are noticeable differences in longitude distributions of dust distributions overland. Although this can be potentially caused by different meteorology conditions, but CALIOP and MODIS observations indicated limited interannual variations [18]. Third, there are differences in how far dust aerosols are transported westerly. Model simulated dust aerosols are transported further west than the observations. Fourth, model simulations show that significant aerosols are transported above 6 km, but CALIOP measurements have difficulties to detect these dust aerosols due to low signal-to-noise ratio for weak aerosol layers or challenges in separating dust with optically thin ice clouds at colder temperatures. Due to the important roles of dust aerosols at these high elevations through both mixed-phase and ice cloud microphysical processes, improved dust detections for CALIOP measurements or new dust measurements over tropical upper troposphere are needed.

3.3. 3-Dimensional Dust Layer Top Variations

To further illustrate the advantage of active measurements, grid-box mean dust layer top distribution during JJA is illustrated in Figure 5. As the dust layer is transported west across the Atlantic Ocean, the dust layer top generally decreases further west, but there are interesting latitude variations at a given longitude. The center part of dust layer shows a smooth top variation with a higher top near the center and decreases to both sides. As illustrated in Figure 6, the southern side of dust layer is bounded by the Inter-Tropical Convergence Zone (ITCZ), and the northern side of the dust layer is bounded by strong westerly flow. Large top variations near the boundaries in Figure 6 clearly indicate that strong vertical mixings occur at these boundaries.

3.4. Dynamical Control of Dust Transport

Many studies explored the processes controlling dust emitted from Africa desert regions and their transportation in the atmosphere [18, 36, 40, 41]. However, these studies are limited with the available data. With passive sensor retrieved optical depths over water, Moulin et al. [34] showed the importance of large-scale climatic control on the dust export by changing precipitation and atmospheric circulations over the regions of dust mobilization and transport. Using a satellite-derived dust index, Engelstaedter and Washington [41] concluded that the dry convection associated with an increase in the occurrence of small-scale high-wind events and vertical velocity are the main processes controlling the annual dust cycle at the West African dust sources. With better dust detection and quantitatively estimation capabilities of CALIOP over land and clouds, long-term CALIOP measurements will offer a new source data to understand the processes controlling African dust.

Based on the first three-year CALIOP data, Figure 6 shows the importance of dynamics in African dust transport. In Figure 6(a), seasonal dust occurrences are overplotted with the 850 mb horizontal wind from NCEP/NCAR reanalysis data [42], which clearly show that the westward transport of African dust is strongly controlled by the prevailing easterlies. The southern bound of dust layer is limited to the positions of the ITCZ, which can be identified in Figure 6(d) as tropical regions with weak updraft velocities the over ocean and stronger updraft velocities over the land. The seasonal shifting of ITCZ and associated strong easterly wind zones determine the track of dust transportation across the Atlantic [40]. During DJF, strong updrafts over South America and high wind over the tropical east Atlantic bring significant African dust to the Amazon region, even though dust over the source regions is mainly confined below 2 km during this period.

Large-scale convergences and associated updrafts (see Figures 6(c) and 6(d)) overland are one of the important factors to control African dust aerosol production. As discussed by Engelstaedter and Washington [41], large-scale convergences favor dry convection developments and cause an increase in the occurrence of small-scale high-wind events to produce more dense dust storms. The updrafts associated with large-scale convergences not only support deep dust layer developments, but also increase dust lifetime by reducing gravity sedimentation [36]. During DJF, a narrow updraft band in the south of 10°N can be clearly identified in Figure 6(d), which collocates with high dust occurrence areas. Compared with DJF, there other three seasons have a more wide distributed dust overland, which could be attributed to multiple convergence zones developed in northern 10°N. Due to strong solar heating, JJA has wide and stronger updraft zones which produces the deepest dust layer in all season.

4. Summary

Mineral dust plays an important role in the climate system processes by affecting radiation, modifying cloud and precipitation processes, and modulating a variety of marine biogeochemical processes. The Sahara is a primary source of dust, the majority of which is transported westward across the Atlantic to the Southeastern United States, the Caribbean Sea, and South America. However, there are large uncertainties in dust sources, sinks, and their long-range transports. To overcome the limitations of quantitative dust detections with passive satellite measurements over the land and low clouds, we use multiyear combined CALIPSO and CloudSat data to show three-dimensional dust distributions over the source regions and their transport pathways across the Atlantic.

To improve dense dust detection, CloudSat radar detected that clouds are used to together with CALIOP measurements to provide more reliable vertical resolved dust occurrence maps. The results show strong seasonal shifts in dust source regions and transportation pathways. Other than DJF, dust aerosols are produced over several large areas scattered across Africa, which result in a wide spread dust overland. These results confirm that both large-scale and small-scale meteorological conditions are important for dust emissions over Africa. The seasonal differences in meteorological and thermodynamical conditions also result in significant differences in dust vertical distributions overland with the deepest dust layer in JJA and the shallowest dust layer in DJF. Although the seasonal north-south shifts of dust bands across the Atlantic are linked with the movements of the ITCZ, the seasonal vertical pathways are also linked to many other factors. When dust aerosols are transported across the Atlantic, the mean dust layer top descends at different rates (about 35 m/degree in JJA, 25 m/degree in SON and MAM, and 10 m/deg in DJF). The amount of dust transported below and above 2-km also changes with seasons. JJA has more dust transported above 2 km and reached further west to impact large areas including the east coast of North America.

The extinction coefficients and particle depolarization ratio of dusty aerosols are also derived to understand dust distribution more quantitatively. Other than DJF, dust aerosols are vertically well mixed overland; however, there are significant horizontal variations. The spatial patterns of dust extinction coefficients and depolarization ratios indicate that the strong mixings occur as dust is transported across the Atlantic. These mixing processes could be challenging for models to simulate. Further studies with CALIOP and airborne measurements are needed to better understand these mixing processes.

To illustrate the potentials to use CALIOP measurements for dust transport model improvements, we use a specific extinction coefficient that is 0.7 m2g−1 to convert dust extinction coefficients to dust mass mixing ratios, which are predicted by most models. As an example, observations are compared with DEAD simulated dust spatial distributions. It is clear that the DEAD dust model captures the general seasonal variations of dust distributions, but several noticeable model biases are identified for further improvements. Although there are up to 30% uncertainties in estimating dust mass loading from CALIOP measurements, it is very valuable to reduce larger intermodal spreads in dust simulations, especially on dust vertical distributions. With two-wavelength CALIOP measurements, further improvements of dust mass loading estimations are possible. Because dust aerosols are efficient IN, validating and improving model simulations above 6 km require better observation capabilities for optically thin dust detection.

The general features of the vertical pathways of dust transport across the Atlantic demonstrate that African dust observed at surface sites within and along the eastern edge of the Atlantic is strongly affected by meteorological conditions. Therefore, surface observations of dust are not, by themselves, sufficient for evaluating simulations of transatlantic transport of African dust. Long-term observations from space-borne LIDARs, together with radiometer measurements, are key for understanding the production and transport of African dust and for improving its representation in the climate and weather models that are needed to quantify the impact of African dust on climate change.

Acknowledgments

This research is supported by the National Natural Science Foundation of China (no. 41075016) and the Major State Basic Research Development Program of China (no. 2011CB403400). Z. Wang is supported by NASA Grant NNX07AQ83G and NNX10AN18G. The authors would like to thank CALIPSO PIs Dave Winkler and Charles Trepte, CloudSat PIs Graeme Stephens and Deborah Vane, and the CALIPSO/CloudSat data group. NCEP Reanalysis data was provided by the NOAA/OAR/ESRL PSD, Boulder, CO, USA, from their website at http://www.esrl.noaa.gov/psd/.