2.31. Dust Model¶
Atmospheric dust is mobilized from the land by wind in the CLM. The most important factors determining soil erodibility and dust emission include the wind friction speed, the vegetation cover, and the soil moisture. The CLM dust mobilization scheme (Mahowald et al. 2006) accounts for these factors based on the DEAD (Dust Entrainment and Deposition) model of Zender et al. (2003). Please refer to the Zender et al. (2003) article for additional information regarding the equations presented in this section.
The total vertical mass flux of dust, (kg m^{2} s^{1}), from the ground into transport bin is given by
(2.31.1)¶
where is a global factor that compensates for the DEAD model’s sensitivity to horizontal and temporal resolution and equals 5 x 10^{4} in the CLM instead of 7 x 10^{4} in Zender et al. (2003). is the source erodibility factor set to 1 in the CLM and serves as a place holder at this time.
The grid cell fraction of exposed bare soil suitable for dust mobilization is given by
(2.31.2)¶
where and are the CLM grid cell fractions of lake (section 2.2.3.3) and snow cover (section 2.8.1), all ranging from zero to one. Not mentioned by Zender et al. (2003), and are the CLM top soil layer liquid water and ice contents (mm) entered as a ratio expressing the decreasing ability of dust to mobilize from increasingly frozen soil. The grid cell fraction of vegetation cover,, is defined as
(2.31.3)¶
where equation applies only for dust mobilization and is not related to the plant functional type fractions prescribed from the CLM input data or simulated by the CLM dynamic vegetation model (Chapter 22). and are the CLM leaf and stem area index values (m ^{2} m^{2}) averaged at the land unit level so as to include all the pfts and the bare ground present in a vegetated land unit. and may be prescribed from the CLM input data (section 2.2.1.4) or simulated by the CLM biogeochemistry model (Chapter 2.20).
The sandblasting mass efficiency (m ^{1}) is calculated as
(2.31.4)¶
where is the mass fraction of clay particles in the soil and %clay is determined from the surface dataset (section 2.2.3.3). corresponds to sand and to sandy loam.
is the total horizontally saltating mass flux (kg m^{1} s^{1}) of “large” particles (Table 2.31.1), also referred to as the vertically integrated streamwise mass flux
(2.31.5)¶
where the saltation constant equals 2.61 and is the atmospheric density (kg m^{3}) (Table 2.2.4), the acceleration of gravity (m s^{2}) (Table 2.2.7). The threshold wind friction speed for saltation (m s^{1}) is
(2.31.6)¶
where is a factor dependent on surface roughness but set to 1 as a place holder for now, and are the density (2650 kg m^{3}) and diameter (75 x 10 m) of optimal saltation particles, and is a factor dependent on soil moisture:
(2.31.7)¶
where
(2.31.8)¶
and
(2.31.9)¶
where for tuning purposes, is the volumetric soil moisture in the top soil layer (m m^{3}) (section 2.7.3), is the density of liquid water (kg m^{3}) (Table 2.2.7), and is the bulk density of soil in the top soil layer (kg m^{3}) defined as in section 2.6.3 rather than as in Zender et al. (2003). from equation is the threshold friction Reynolds factor
(2.31.10)¶
and is the threshold friction Reynolds number approximation for optimally sized particles
(2.31.11)¶
In (2.31.5) , is defined as the wind friction speed (m s^{1}) accounting for the Owen effect (Owen 1964)
(2.31.12)¶
where is the CLM wind friction speed (m s^{1}), also known as friction velocity (section 2.5.1), is the 10m wind speed (m s^{1}) calculated as the wind speed at the top of the canopy in section 4.3 of Bonan (1996) but here for 10 m above the ground, and is the threshold wind speed at 10 m (m s^{1})
(2.31.13)¶
In equation we sum over source modes where is the mass fraction of each source mode carried in each of :math:`J=4` transport bins
(2.31.14)¶
where , , and are the mass fraction, mass median diameter, and geometric standard deviation assigned to each particle source mode (Table 2.31.1), while and are the minimum and maximum diameters (m) in each transport bin (Table 2.31.2).
(fraction) 
(m) 


1 
0.036 
0.832 x 10 
2.1 
2 
0.957 
4.820 x 10 
1.9 
3 
0.007 
19.38 x 10 
1.6 
(m) 
(m) 


1 
0.1 x 10 
1.0 x 10 
2 
1.0 x 10 
2.5 x 10 
3 
2.5 x 10 
5.0 x 10 
4 
5.0 x 10 
10.0 x 10 