Regional Saharan dust modelling during the SAMUM 2006 campaign

The regional dust model system LM-MUSCAT-DES was developed in the framework of the SAMUM project. Using the unique comprehensive data set of near-source dust properties during the 2006 SAMUM field campaign, the performance of the model system is evaluated for two time periods in May and June 2006. Dust optical thicknesses, number size distributions and the position of the maximum dust extinction in the vertical profiles agree well with the observations. However, the spatio-temporal evolution of the dust plumes is not always reproduced due to inaccuracies in the dust source placement by the model. While simulated winds and dust distributions are well matched for dust events caused by dry synoptic-scale dynamics, they are often misrepresented when dust emissions are caused by moist convection or influenced by small-scale topography that is not resolved by the model. In contrast to long-range dust transport, in the vicinity of source regions the model performance strongly depends on the correct prediction of the exact location of sources. Insufficiently resolved vertical grid spacing causes the absence of inversions in the model vertical profiles and likely explains the absence of the observed sharply defined dust layers.     

Spatial distribution and optical properties of Saharan dust observed by airborne high spectral resolution lidar during SAMUM 2006

Airborne measurements of pure Saharan dust extinction and backscatter coefficients, the corresponding lidar ratio and the aerosol optical thickness (AOT) have been performed during the Saharan Mineral Dust Experiment 2006, with a high spectral resolution lidar. Dust layers were found to range from ground up to 4–6 km above sea level (asl). Maximum AOT values at 532 nm, encountered within these layers during the DLR Falcon research flights were 0.50–0.55. A significant horizontal variability of the AOT south of the High Atlas mountain range was observed even in cases of a well-mixed dust layer. High vertical variations of the dust lidar ratio of 38–50 sr were observed in cases of stratified dust layers. The variability of the lidar ratio was attributed to dust advection from different source regions. The aerosol depolarization ratio was about 30% at 532 nm during all measurements and showed only marginal vertical variations.     

Saharan dust absorption and refractive index from aircraft-based observations during SAMUM 2006

During the Saharan Mineral Dust Experiment (SAMUM) conducted in summer 2006 in southeast Morocco, the complex refractive index of desert dust was determined from airborne measurements of particle size distributions and aerosol absorption coefficients at three different wavelengths in the blue (467 nm), green (530 nm) and red (660 nm) spectral regions. The vertical structure of the dust layers was analysed by an airborne high spectral resolution lidar (HSRL). The origin of the investigated dust layers was estimated from trajectory analyses, combined with Meteosat 2nd Generation (MSG) scenes and wind field data analyses. The real part n of the dust refractive index was found almost constant with values between 1.55 and 1.56, independent of the wavelength. The values of the imaginary part k varied between the blue and red spectral regions by a factor of three to ten depending on the dust source region. Absolute values of k ranged from 3.1 × 10⁻³ to 5.2 × 10⁻³ at 450 nm and from 0.3 × 10⁻³ to 2.5 × 10⁻³ at 700 nm. Groupings of k values could be attributed to different source regions.     

Depolarization ratio profiling at several wavelengths in pure Saharan dust during SAMUM 2006

Vertical profiles of the linear particle depolarization ratio of pure dust clouds were measured during the Saharan Mineral Dust Experiment (SAMUM) at Ouarzazate, Morocco (30.9°N, –6.9°E), close to source regions in May–June 2006, with four lidar systems at four wavelengths (355, 532, 710 and 1064 nm). The intercomparison of the lidar systems is accompanied by a discussion of the different calibration methods, including a new, advanced method, and a detailed error analysis. Over the whole SAMUM periode pure dust layers show a mean linear particle depolarization ratio at 532 nm of 0.31, in the range between 0.27 and 0.35, with a mean Ångström exponent (AE, 440–870 nm) of 0.18 (range 0.04–0.34) and still high mean linear particle depolarization ratio between 0.21 and 0.25 during periods with aerosol optical thickness less than 0.1, with a mean AE of 0.76 (range 0.65–1.00), which represents a negative correlation of the linear particle depolarization ratio with the AE. A slight decrease of the linear particle depolarization ratio with wavelength was found between 532 and 1064 nm from 0.31 ± 0.03 to 0.27 ± 0.04.     

Airborne measurements of dust layer properties, particle size distribution and mixing state of Saharan dust during SAMUM 2006

The Saharan Mineral Dust Experiment (SAMUM) was conducted in May/June 2006 in southern Morocco. As part of SAMUM, airborne in situ measurements of the particle size distribution in the diameter range 4 nm < D _p < 100 μm were conducted. The aerosol mixing state was determined below D _p < 2.5 μm. Furthermore, the vertical structure of the dust layers was investigated with a nadir-looking high spectral resolution lidar (HSRL). The desert dust aerosol exhibited two size regimes of different mixing states: below 0.5 μm, the particles had a non-volatile core and a volatile coating; larger particles above 0.5 μm consisted of non-volatile components and contained light absorbing material. In all cases, particles larger than 10 μm were present, and in 80% of the measurements no particles larger than 40 μm were present. The abundance of large particles showed almost no height dependence. The effective diameter D _eff in the dust plumes investigated showed two main ranges: the first range of D _eff peaked around 5 μm and the second range of D _eff around 8 μm. The two ranges of D _eff suggest that it may be inadequate to use one average effective diameter or one parametrization for a typical dust size distribution.     

Vertical profiling of Saharan dust with Raman lidars and airborne HSRL in southern Morocco during SAMUM

Three ground-based Raman lidars and an airborne high-spectral-resolution lidar (HSRL) were operated during SAMUM 2006 in southern Morocco to measure height profiles of the volume extinction coefficient, the extinction-to-backscatter ratio and the depolarization ratio of dust particles in the Saharan dust layer at several wavelengths. Aerosol Robotic Network (AERONET) Sun photometer observations and radiosoundings of meteorological parameters complemented the ground-based activities at the SAMUM station of Ouarzazate. Four case studies are presented. Two case studies deal with the comparison of observations of the three ground-based lidars during a heavy dust outbreak and of the ground-based lidars with the airborne lidar. Two further cases show profile observations during satellite overpasses on 19 May and 4 June 2006. The height resolved statistical analysis reveals that the dust layer top typically reaches 4–6 km height above sea level (a.s.l.), sometimes even 7 km a.s.l.. Usually, a vertically inhomogeneous dust plume with internal dust layers was observed in the morning before the evolution of the boundary layer started. The Saharan dust layer was well mixed in the early evening. The 500 nm dust optical depth ranged from 0.2–0.8 at the field site south of the High Atlas mountains, Ångström exponents derived from photometer and lidar data were between 0–0.4. The volume extinction coefficients (355, 532 nm) varied from 30–300 Mm⁻¹ with a mean value of 100 Mm⁻¹ in the lowest 4 km a.s.l.. On average, extinction-to-backscatter ratios of 53–55 sr (±7–13 sr) were obtained at 355, 532 and 1064 nm.     

Vertical profiling of convective dust plumes in southern Morocco during SAMUM

Lifting of dust particles by dust devils and convective plumes may significantly contribute to the global mineral dust budget. During the Saharan Mineral Dust Experiment (SAMUM) in May–June 2006 vertical profiling of dusty plumes was performed for the first time. Polarization lidar observations taken at Ouarzazate (30.9°N, 6.9°W, 1133 m height above sea level) are analyzed. Two cases with typical and vigorous formation of convective plumes and statistical results of 5 d are discussed. The majority of observed convective plumes have diameters on order of 100–400 m. Most of the plumes (typically 50–95%) show top heights <1 km or 0.3DLH with the Saharan dust layer height DLH of typically 3–4 km. Height-to-diameter ratio is mostly 2–10. Maximum plume top height ranges from 1.1 to 2.9 km on the 5 d. 5–26 isolated plumes and clusters of plumes per hour were detected. A low dust optical depth (<0.3) favours plume evolution. Observed surface, 1 and 2–m air temperatures indicate that a difference of 17–20 K between surface and 2-m air temperature and of 0.9–1 K between the 1 and 2-m temperatures are required before convective plumes develop. Favourable horizontal wind speeds are 2–7 m s⁻¹.     

Spectral aerosol optical depth characterization of desert dust during SAMUM 2006

The aerosol optical depth (AOD) in the range 340–1550 nm was monitored at Ouarzazate (Morocco) during the Saharan Mineral Dust Experiment (SAMUM) experiment in May–June 2006. Two different sun photometers were used for this purpose. The mean AOD at 500 nm was 0.28, with a maximum of 0.83, and the mean Ångström exponent (AE) was 0.35. The aerosol content over the site changed alternatively from very low turbidity, associated to Atlantic air masses, to moderate dust load, associated to air masses arriving in the site from Algeria, Tunisia and Libya. The dusty conditions were predominant in the measurement period (78% of data), with AOD (500 nm) above 0.15 and AE below 0.4. The spectral features of the AOD under dusty conditions are discussed. Air mass back trajectory analysis is carried out to investigate the origin and height patterns of the dust loaded air masses. The advection of dust occurred mainly at atmospheric heights below 3000 m, where east flow is the predominant. At the 5000 m level, the air masses originate mainly over the Atlantic Ocean. Finally the Optical Properties of Aerosols and Clouds (OPAC) model is used to perform a set of simulations with different aerosol mixtures to illustrate the measured AOD and AE values under varying dust concentrations, and a brief comparison with other measurement sites is presented.     

Measurements of desert dust optical characteristics at Porte au Sahara during SAMUM

Main optical characteristics of desert dust, such as phase function and single scattering albedo, have been derived from combinations of sun-/sky-radiometer and satellite measurements during the SAMUM experiment (10 May–10 June 2006) at the site Porte au Sahara (30.237°N, 5.607°W) in South Morocco. Scattering phase functions have been retrieved using combined data of spectral aerosol optical thickness (AOT) and spectral sky brightness in the almucantar, considering non-spherical light scattering. Intercomparisons of modelled top-of-atmosphere (TOA) reflectance with satellite observations of the Medium Resolution Imaging Spectrometer (MERIS) and Scanning Imaging Absorption Spectrometer for Atmospheric Chartography () instrument have been used for the estimation of spectral single scattering albedo. For the radiative transfer calculations scattering phase functions and AOT from ground-based observations have been used. The spectral single scattering albedo ranges from 0.93 in the blue to 0.98 at 753 nm.     

Spectral surface albedo over Morocco and its impact on radiative forcing of Saharan dust

In May–June 2006, airborne and ground-based solar (0.3–2.2 μm) and thermal infrared (4–42 μm) radiation measurements have been performed in Morocco within the Saharan Mineral Dust Experiment (SAMUM). Upwelling and downwelling solar irradiances have been measured using the Spectral Modular Airborne Radiation Measurement System (SMART)-Albedometer. With these data, the areal spectral surface albedo for typical surface types in southeastern Morocco was derived from airborne measurements for the first time. The results are compared to the surface albedo retrieved from collocated satellite measurements, and partly considerable deviations are observed. Using measured surface and atmospheric properties, the spectral and broad-band dust radiative forcing at top-of-atmosphere (TOA) and at the surface has been estimated. The impact of the surface albedo on the solar radiative forcing of Saharan dust is quantified. In the SAMUM case of 19 May 2006, TOA solar radiative forcing varies by 12 W m⁻² per 0.1 surface-albedo change. For the thermal infrared component, values of up to +22 W m⁻² were derived. The net (solar plus thermal infrared) TOA radiative forcing varies between −19 and +24 W m⁻² for a broad-band solar surface albedo of 0.0 and 0.32, respectively. Over the bright surface of southeastern Morocco, the Saharan dust always has a net warming effect.     

In situ measurements of optical properties at Tinfou (Morocco) during the Saharan Mineral Dust Experiment SAMUM 2006

In situ measurements of optical and physical properties of mineral dust were performed at the outskirts of the Saharan Desert in the framework of the Saharan Mineral Dust Experiment part 1 (SAMUM-1). Goals of the field study were to achieve information on the extent and composition of the dust particle size distribution and the optical properties of dust at the ground. For the particle number size distribution, measured with a DMPS/APS, a size dependent dynamic shape factor was considered. The mean refractive index of the particles in this field study is  1.53–4.1 × 10⁻³ i   at 537 nm wavelength and  1.53–3.1 × 10⁻³ i   at 637 nm wavelength derived from measurements of scattering and absorption coefficients, as well as the particle size distribution. Whereas the real part of the refractive index is rather constant, the imaginary part varies depending on the mineral dust concentrations. For high dust concentration the single scattering albedo is primarily influenced by iron oxide and is  0.96 ± 0.02  and  0.98 ± 0.01  at 537 nm and 637 nm wavelength, respectively. During low dust concentration the single scattering albedo is more influenced by a soot-type absorber and is  0.89 ± 0.02  and  0.93 ± 0.01  for the same wavelengths.     

EARLINET observations of the 14–22-May long-range dust transport event during SAMUM 2006: validation of results from dust transport modelling

We observed a long-range transport event of mineral dust from North Africa to South Europe during the Saharan Mineral Dust Experiment (SAMUM) 2006. Geometrical and optical properties of that dust plume were determined with Sun photometer of the Aerosol Robotic Network (AERONET) and Raman lidar near the North African source region, and with Sun photometers of AERONET and lidars of the European Aerosol Research Lidar Network (EARLINET) in the far field in Europe. Extinction-to-backscatter ratios of the dust plume over Morocco and Southern Europe do not differ. Ångström exponents increase with distance from Morocco. We simulated the transport, and geometrical and optical properties of the dust plume with a dust transport model. The model results and the experimental data show similar times regarding the appearance of the dust plume over each EARLINET site. Dust optical depth from the model agrees in most cases to particle optical depth measured with the Sun photometers. The vertical distribution of the mineral dust could be satisfactorily reproduced, if we use as benchmark the extinction profiles measured with lidar. In some cases we find differences. We assume that insufficient vertical resolution of the dust plume in the model calculations is one reason for these deviations.     

Spectral absorption coefficients and imaginary parts of refractive indices of Saharan dust during SAMUM-1

During the SAMUM-1 experiment, absorption coefficients and imaginary parts of refractive indices of mineral dust particles were investigated in southern Morocco. Main absorbing constituents of airborne samples were identified to be iron oxide and soot. Spectral absorption coefficients were measured using a spectral optical absorption photometer (SOAP) in the wavelength range from 300 to 800 nm with a resolution of 50 nm. A new method that accounts for a loading-dependent correction of fibre filter based absorption photometers, was developed. The imaginary part of the refractive index was determined using Mie calculations from 350 to 800 nm. The spectral absorption coefficient allowed a separation between dust and soot absorption. A correlation analysis showed that the dust absorption coefficient is correlated ( R ² up to 0.55) with the particle number concentration for particle diameters larger than 0.5 μm, whereas the coefficient of determination R ² for smaller particles is below 0.1. Refractive indices were derived for both the total aerosol and a dust aerosol that was corrected for soot absorption. Average imaginary parts of refractive indices of the entire aerosol are 7.4 × 10⁻³, 3.4 × 10⁻³ and 2.0 × 10⁻³ at wavelengths of 450, 550 and 650 nm. After a correction for the soot absorption, imaginary parts of refractive indices are 5.1 × 10⁻³, 1.6 × 10⁻³ and 4.5 × 10⁻⁴.     

Relative effects of solar and infrared radiative forcing in a mesoscale model

A new efficient parameterization scheme for solar short-wave radiative heating, as a component of the net radiative effects in the atmosphere, is tested in a three-dimensional mesoscale model. This model is designed with moist convective processes in mind, so that the radiative parameterization (solar plus thermal infrared) are interactive with the cloud field. Previous work by the authors with only an infrared scheme has demonstrated that cloud-radiation interactions are characterized by strong cloud-top cooling, leading to upper cloud-layer destabilization. The effects of including solar heating are to modulate the strength of the strong infrared cooling, thereby leading to weaker interactions between clouds, radiation, and mesoscale fields. The present study shows that even on the mesoscale and for relatively short time-spans, radiative processes in the presence of clouds are not negligible.As a further step, a simple fractional cloud cover parameterization is introduced and the model response is compared with results omitting this parameterization.     

Microphysical and radiative properties of mixed-phase altocumulus: A model evaluation of glaciation effects

Using a cloud model with explicit microphysics and radiation, we evaluate the microphysical changes in a supercooled liquid altocumulus cloud with increasing ice content, until glaciation occurs. The properties of the ice and water particle constituents are resolved independently, revealing the relative radiative contributions from the water phase source cloud versus the ice phase virga. Cloud contents are also converted to millimeter-wave radar reflectivities to shed light on the ability of such radars to study these ubiquitous clouds. The results show that the radiative and radar backscattering properties of mixed-phase clouds are dominated by the cloud droplets and ice crystals they contain, respectively.     

The linearization and adjoint of radiation transfer processes in the NMC spectral model part I: Solar radiative transfer

Summary Linearization and adjoint-model derivation for the solar radiation transfer codes in the NMC spectral model have been carried out. Verification of the validity of resulting tangent linear model and the correctness of the corresponding adjoint have been performed. Applications of derived adjoint model are considered, including parameter estimation for inputs to solar radiation codes with aid of the physics (i.e., the solar radiation codes) and a sensitivity study of the downward solar radiation flux at the earth surface with respect to water vapor amount at various heights.With 5 Figures     

Direct effects and feedback of desert dust on the climate of the Arabian Peninsula during the wet season: a regional climate model study

We investigate the dust radiative forcing and its feedback on the Arabian Peninsula’s wet season climate using the International Centre for Theoretical Physics-Regional Climate Model (ICTP-RegCM4). We have found that the dust plumes exert a negative (positive) radiative forcing at the surface (top of the atmosphere) by reducing incoming solar radiation reaching the ground and locally heating up the atmosphere column. Consequently, the surface air temperature is cooler, hence indicating a decrease in the warm bias and an increase in the temperature gradient. This reduces the geopotential heights and enhances the low-level wind convergence, suggesting stronger upward motion. These changes increase evaporation, the difference between precipitation and evaporation in the atmosphere and rainfall over the Peninsula, indicating an intensification of the hydrologic cycle. The decrease in the precipitation dry bias and the large reduction in the temperature warm bias caused by the impact of dust over the entire Peninsula represent a significant success for the RegCM4 simulation. Therefore, the inclusion of dust in the simulation of the Arabian Peninsula’s climate for the wet season contributes to an improved performance of this regional climate model over the region.     

Incorporation of internal fluctuations in a meandering plume model of concentration fluctuations

A meandering plume model that explicitly incorporates the effects of small-scale structure in the instantaneous plume has been formulated. The model requires the specification of two physically based input parameters; namely, the meander ratio,M, which is dependent on the ratio of the meandering plume dispersion to the instantaneous relative plume dispersion and, a relative in-plume fluctuation measure,k, that is related inversely to the fluctuation intensity in relative coordinates. Simple analytical expressions for crosswind profiles of the higher moments (including the important shape parameters such as fluctuation intensity, skewness, and kurtosis) and for the concentration pdf have been derived from the model. The model has been tested against some field data sets, indicating that it can reproduce many key aspects of the observed behavior of concentration fluctuations, particularly with respect to modeling the change in shape of the concentration pdf in the crosswind direction.List of Symbols C Mean concentration in absolute coordinates - C _r Mean concentration in relative coordinates - C₀ Centerline mean concentration in absolute coordinates - C _r,0 Centerline mean concentration in relative coordinates - f Probability density function of concentration in absolute coordinates - f _c Probability density function of plume centroid position - f _r Probability density function of concentration in relative coordinates - i Absolute concentration fluctuation intensity (standard deviation to mean ratio) - i _r Relative concentration fluctuation intensity (standard deviation to mean ratio) - k Relative in-plume fluctuation measure:k=1/i _r ² - K Concentration fluctuation kurtosis - M Meander ratio of meandering plume variance to relative plume variance - S Concentration fluctuation skewness - x Downwind distance from source - y Crosswind distance from mean-plume centerline - z Vertical distance above ground - Instantaneous (random) concentration - Crosswind dispersion ofnth concentration moment about zero - _ny Mean-plume crosswind (absolute) dispersion - _y Plume centroid (meandering) dispersion in crosswind direction - _y,c Instantaneous plume crosswind (relative) dispersion - Normalized mean concentration in absolute coordinates:C/C ₀ - Particular value taken on by instantaneous concentration,      

A parameterization of radiative fluxes suitable for use in a statistical-dynamical model

Summary A parameterization of shortwave and longwave radiation fluxes derived from detailed radiative transfer models is included in a global primitive equation statistical-dynamical model (SDM) with two bulk atmospheric layers. The model is validated comparing the model simulations with the observed mean annual and seasonal zonally averaged climate. The results show that the simulation of the shortwave and longwave radiation fluxes matches well with the observations. The SDM variables such as surface and 500 hPa temperatures, zonal winds at 250 hPa and 750 hPa, vertical velocity at 500 hPa and precipitation are also in good agreement with the observations. A comparison between the results obtained with the present SDM and those with the previous version of the model indicates that the model results improved when the parameterization of the radiative fluxes based on detailed radiative transfer models are included into the SDM.The SDM is used to investigate its response to the greenhouse effect. Sensitivity experiments regarding the doubling of CO₂ and the changing of the cloud amount and height are performed. In the case 2×CO₂ the model results are consistent with those obtained from GCMs, showing a warming of the climate system. An enhancement of the greenhouse effect is also noted when the cloud layer is higher. However, an increase of the cloud amount in all the latitude belts provokes an increase of the surface temperature near poles and a decrease in all the other regions. This suggests that the greenhouse effect overcomes the albedo effect in the polar latitudes and the opposite occurs in other regions. In all the experiments the changes in the surface temperature are larger near poles, mainly in the Southern Hemisphere.With 8 Figures