The Microphysics1M.jl
module describes a 1-moment bulk parameterization of
cloud microphysical processes.
The module is based on the ideas of
Kessler1995,
Grabowski1998
and Kaul2015.
The cloud microphysics variables are expressed as specific humidities:
q_tot
- total water specific humidity,q_vap
- water vapor specific humidity,q_liq
- cloud water specific humidity,q_ice
- cloud ice specific humidity,q_rai
- rain specific humidity,q_sno
- snow specific humidity.
Particles are assumed to follow power-law relationships involving the mass(radius),
denoted by m(r)
, the cross section(radius), denoted by a(r)
, and the
terminal velocity(radius), denoted by v_{term}(r)
, respectively.
The coefficients are defined in the
ClimaParams.jl package
and are shown in the table below.
For rain and ice they correspond to spherical liquid water drops
and ice particles, respectively.
There is no assumed particle shape for snow, and the relationships are
empirical.
where:
r
is the particle radius,r_0
is the typical particle radius used to nondimensionalize,m_0, \, m_e, \, a_0, \, a_e, \, v_0, \, v_e \,
are the default coefficients,\chi_m
,\Delta_m
,\chi_a
,\Delta_a
,\chi_v
,\Delta_v
are the coefficients that can be used during model calibration to expand around the default values. Without calibration all\chi
parameters are set to 1 and all\Delta
parameters are set to 0.
The above coefficients, similarly to all other microphysics parameters, are not hardcoded in the final microphysics parameterizations. The goal is to allow easy flexibility when calibrating the scheme. With that said, the assumption about the shape of the particles is used three times when deriving the microphysics formulae:
- The mass transfer equation (\ref{eq:mass_rate}) used in snow autoconversion, rain evaporation, snow sublimation and snow melt rates is derived assuming spherical particle shapes. To correct for non-spherical shape it should be multiplied by a function of the particle aspect ratio.
- The geometric collision kernel used for deriving rain-snow accretion rate assumes that both colliding particles are spherical. It does not take into account the reduced cross-section of snow particles that is used when modeling snow - cloud liquid water and snow - cloud ice accretion.
- In the definition of the Reynolds number that is used when computing ventilation factors.
symbol | definition | units | default value | reference |
---|---|---|---|---|
r_0^{rai} |
typical rain drop radius | m |
10^{-3} |
|
m_0^{rai} |
coefficient in m(r) for rain |
kg |
\frac{4}{3} \, \pi \, \rho_{water} \, r_0^3 |
|
m_e^{rai} |
exponent in m(r) for rain |
- | 3 |
|
a_0^{rai} |
coefficient in a(r) for rain |
m^2 |
\pi \, r_0^2 |
|
a_e^{rai} |
exponent in a(r) for rain |
- | 2 |
|
v_e^{rai} |
exponent in v_{term}(r) for rain |
- | 0.5 |
|
r_0^{ice} |
typical ice crystal radius | m |
10^{-5} |
|
m_0^{ice} |
coefficient in m(r) for ice |
kg |
\frac{4}{3} \, \pi \, \rho_{ice} \, r_0^3 |
|
m_e^{ice} |
exponent in m(r) for ice |
- | 3 |
|
r_0^{sno} |
typical snow crystal radius | m |
10^{-3} |
|
m_0^{sno} |
coefficient in m(r) for snow |
kg |
0.1 \, r_0^2 |
eq (6b) Grabowski1998 |
m_e^{sno} |
exponent in m(r) for snow |
- | 2 |
eq (6b) Grabowski1998 |
a_0^{sno} |
coefficient in a(r) for snow |
m^2 |
0.3 \pi \, r_0^2 |
\alpha in eq(16b) Grabowski1998. |
a_e^{sno} |
exponent in a(r) for snow |
- | 2 |
|
v_0^{sno} |
coefficient in v_{term}(r) for snow |
\frac{m}{s} |
2^{9/4} r_0^{1/4} |
eq (6b) Grabowski1998 |
v_e^{sno} |
exponent in v_{term}(r) for snow |
- | 0.25 |
eq (6b) Grabowski1998 |
where:
\rho_{water}
is the density of water,\rho_{ice}
is the density of ice.
The terminal velocity of an individual rain drop is defined by the balance
between the gravitational acceleration (taking into account
the density difference between water and air) and the drag force.
Therefore the v_0^{rai}
is defined as
where:
g
is the gravitational acceleration,C_{drag}
is the drag coefficient,\rho
is the density of air.
!!! note Assuming a constant drag coefficient is an approximation and it should be size and flow dependent, see for example here. We are in the process of updating the 1-moment microphysics scheme to formulae from Chen2022. Other possibilities: Khvorostyanov2002 or Karrer2020
Chen2022 provides a terminal velocity parameterisation
based on an empirical fit to a high accuracy model.
The terminal velocity depends on particle shape, size and density,
consideres the deformation effects of large rain drops,
as well as size-specific air density dependence.
The fall speed of individual particles
where:
D
is the particle diameter,a_i
,b_i
,c_i
are the free parameters,\phi
is the aspect ratio, and\kappa
is a parameter that depends on the particle shape (\kappa=0
for spheres,\kappa=-1/3
for oblate and\kappa=1/6
for prolate spheroids).
For ice and snow j=2
and for rain j=3
, to account for deformation at larger sizes.
For rain and ice we assume \phi=1
(spherical).
For snow we assume \kappa = -1/3
and
find the aspect ratio that is consistent with the assumed m(r)
and a(r)
relationships.
The aspect ratio is defined as:
where:
a
is the basal plane axial half-length, andc
is perpendicular to the basal plane.
The volume of a spheroid can be represented as V_p = 4\pi/3 \; a^2 c
and the area can be represented as A_p = \pi a c
.
It follows that
c = (4A_p^2) / (3 \pi V_p)
,
a = (3V_p) / (4A_p)
, and
\phi = (16 A_p^3) / (9 \pi V_p^2)
.
The volume and area are defined by the assumed power-law size relations
V_p = m(r) / (\rho_{ice})
, A_p = a(r)
.
As a result the terminal velocity of individual snow particles as:
where
Here we plot the terminal velocity formulae from the current default 1-moment scheme and Chen2022. We also show the aspect ratio of snow particles.
include("plots/TerminalVelocityComparisons.jl")
The rain a_i
, b_i
, and c_i
are listed in the table below.
The formula is applicable when D > 0.1 mm
,
q = e^{0.115231 \; \rho_a}
, where \rho_a
is air density [kg/m3].
The units are: [v] = m/s, [D]=mm, [a_i
] = mm^{-b_i} m/s
, [b_i
] is dimensionless, [c_i
] = 1/mm.
i |
a_i |
b_i |
c_i |
---|---|---|---|
1 | 0.044612 \; q |
2.2955 \; -0.038465 \; \rho_a |
0 |
2 | -0.263166 \; q |
2.2955 \; -0.038465 \; \rho_a |
0.184325 |
3 | 4.7178 \; q \; (\rho_a)^{-0.47335} |
1.1451 \; -0.038465 \; \rho_a |
0.184325 |
The ice and snow a_i
, b_i
, and c_i
are listed in the table below.
The formula is applicable when D < 0.625 mm
.
i |
a_i |
b_i |
c_i |
---|---|---|---|
1 | E_s (\rho_a)^{A_s} |
B_s + C_s \rho_a |
0 |
2 | F_s (\rho_a)^{A_s} |
B_s + C_s \rho_a |
G_s |
Coefficient | Formula |
---|---|
A_s |
0.00174079 \log{(\rho_{ice})}^2 − 0.0378769 \log{(\rho_{ice})} - 0.263503 |
B_s |
(0.575231 + 0.0909307 \log{(\rho_{ice})} + 0.515579 / \sqrt{\rho_{ice}})^{-1} |
C_s |
-0.345387 + 0.177362 \, \exp{(-0.000427794 \rho_{ice})} + 0.00419647 \sqrt{\rho_{ice}} |
E_s |
-0.156593 - 0.0189334 \log{(\rho_{ice})}^2 + 0.1377817 \sqrt{\rho_{ice}} |
F_s |
- \exp{[-3.35641 - 0.0156199 \log{\rho_{ice}}^2 + 0.765337 \log{\rho_{ice}}]} |
G_s |
(-0.0309715 + 1.55054 / \log{(\rho_{ice})} - 0.518349 log{(\rho_{ice})} / \rho_{ice})^{-1} |
The particle size distributions are assumed to follow Marshall-Palmer distribution Marshall1948 eq. 1:
where:
n_{0}
and\lambda
are the Marshall-Palmer distribution parameters.
The n_0
for rain and ice is assumed constant.
The n_0
for snow is defined as
where:
\mu^{sno}
and\nu^{sno}
are the coefficients\rho_0
is the typical air density used to nondimensionalize the equation and is equal to1 \, kg/m^3
The coefficients are defined in ClimaParams.jl package and are shown in the table below.
symbol | definition | units | default value | reference |
---|---|---|---|---|
n_{0}^{rai} |
rain drop size distribution parameter | \frac{1}{m^4} |
16 \cdot 10^6 |
eq (2) Marshall1948 |
n_{0}^{ice} |
cloud ice size distribution parameter | \frac{1}{m^4} |
2 \cdot 10^7 |
bottom of page 4396 Kaul2015 |
\mu^{sno} |
snow size distribution parameter coefficient | \frac{1}{m^4} |
4.36 \cdot 10^9 \, \rho_0^{\nu^{sno}} |
eq (A1) Kaul2015 |
\nu^{sno} |
snow size distribution parameter exponent | - |
0.63 |
eq (A1) Kaul2015 |
The \lambda
parameter is defined as
where:
q
is rain, snow or ice specific humidity\chi_m
,m_0
,m_e
,\Delta_m
,r_0
, andn_0
are the corresponding mass(radius) and size distribution parameters\Gamma()
is the gamma function
The cloud-ice size distribution is used when computing snow autoconversion rate and rain sink due to accretion. In other derivations cloud ice, similar to cloud liquid water, is treated as continuous.
!!! note
- Do we want to keep the n_0
for rain constant
and n_0
for snow empirical?
- Do we want to test different size distributions?
Here we plot the Marshall-Palmer particle size distribution for 4 different values for the rain specific humidity (q_rai).
include("plots/MarshallPalmer_distribution.jl")
Parameterized processes include:
- autoconversion of rain and snow,
- accretion,
- evaporation of rain water,
- sublimation, vapor deposition and melting of snow,
- sedimentation of rain and snow with mass weighted average terminal velocity (cloud water and cloud ice are part of the working fluid and do not sediment).
Parameters used in the parameterization are defined in ClimaParams.jl package. They consist of:
symbol | definition | units | default value | reference |
---|---|---|---|---|
C_{drag} |
rain drop drag coefficient | - | 0.55 |
C_{drag} is such that the mass averaged terminal velocity is close to Grabowski1996 |
\tau_{acnv\_rain} |
cloud liquid to rain water autoconversion timescale | s |
10^3 |
eq (5a) Grabowski1996 |
\tau_{acnv\_snow} |
cloud ice to snow autoconversion timescale | s |
10^2 |
|
q_{liq\_threshold} |
cloud liquid to rain water autoconversion threshold | - | 5 \cdot 10^{-4} |
eq (5a) Grabowski1996 |
q_{ice\_threshold} |
cloud ice snow autoconversion threshold | - | 1 \cdot 10^{-6} |
|
r_{is} |
threshold particle radius between ice and snow | m |
62.5 \cdot 10^{-6} |
abstract Harrington1995 |
E_{lr} |
collision efficiency between rain drops and cloud droplets | - | 0.8 |
eq (16a) Grabowski1998 |
E_{ls} |
collision efficiency between snow and cloud droplets | - | 0.1 |
Appendix B Rutledge1983 |
E_{ir} |
collision efficiency between rain drops and cloud ice | - | 1 |
Appendix B Rutledge1984 |
E_{is} |
collision efficiency between snow and cloud ice | - | 0.1 |
bottom page 3649 Morrison2008 |
E_{rs} |
collision efficiency between rain drops and snow | - | 1 |
top page 3650 Morrison2008 |
a_{vent}^{rai}, b_{vent}^{rai} |
rain drop ventilation factor coefficients | - | 1.5 \; ,\; 0.53 |
chosen such that at q_{tot}=15 g/kg and T=288K the evap. rate is close to empirical evap. rate in Grabowski1996 |
a_{vent}^{sno}, b_{vent}^{sno} |
snow ventilation factor coefficients | - | 0.65 \; ,\; 0.44 |
eq (A19) Kaul2015 |
K_{therm} |
thermal conductivity of air | \frac{J}{m \; s \; K} |
2.4 \cdot 10^{-2} |
|
\nu_{air} |
kinematic viscosity of air | \frac{m^2}{s} |
1.6 \cdot 10^{-5} |
|
D_{vapor} |
diffusivity of water vapor | \frac{m^2}{s} |
2.26 \cdot 10^{-5} |
The ventilation factor parameterizes the increase in the mass and heat exchange for falling particles. Following SeifertBeheng2006 eq. 24 the ventilation factor is defined as:
where:
a_{vent}
,b_{vent}
are coefficients,N_{Sc}
is the Schmidt number,N_{Re}
is the Reynolds number of a falling particle.
The Schmidt number is assumed constant:
where:
\nu_{air}
is the kinematic viscosity of air,D_{vapor}
is the diffusivity of water.
The Reynolds number of a spherical drop is defined as:
Applying the terminal velocity(radius) relationship results in
The mass weighted terminal velocity v_t
(following Ogura1971) is:
Integrating the default 1-moment m(r)
and v_{term}(r)
relationships
over the assumed Marshall-Palmer distribution results in group terminal velocity:
Integrating Chen2022 formulae for rain and ice over the assumed Marshall-Palmer size distribution, results in group terminal velocity:
Finally, integrating Chen2022 formulae for snow over the assumed Marshall-Palmer distribution, results in group terminal velocity:
where:
and
Rain autoconversion defines the rate of conversion form cloud liquid water to rain water due to collisions between cloud droplets. It is parameterized following Kessler1995:
where:
q_{liq}
- liquid water specific humidity,\tau_{acnv\_rain}
- timescale,q_{liq\_threshold}
- autoconversion threshold.
!!! note This is the simplest possible autoconversion parameterization. It would be great to implement others and test the impact on precipitation. See for example Wood2005 Table 1 for other simple choices.
Snow autoconversion defines the rate of conversion form cloud ice to snow due
the growth of cloud ice by water vapor deposition.
It is defined as the change of mass of cloud ice for cloud ice particles
larger than threshold r_{is}
.
See Harrington1995
for derivation.
The \frac{dm}{dt}
is obtained by solving the water vapor diffusion equation
in spherical coordinates and linking the changes in temperature at the drop
surface to the changes in saturated vapor pressure via the Clausius-Clapeyron
equation, following
Mason2010.
For the simplest case of spherical particles and not taking into account ventilation effects:
where:
S(q_{vap}, q_{vap}^{sat}) = \frac{q_{vap}}{q_{vap}^{sat}}
is saturation,q_{vap}^{sat}
is the saturation vapor specific humidity (over ice in this case),G(T) = \left(\frac{L_s}{KT} \left(\frac{L_s}{R_v T} - 1 \right) + \frac{R_v T}{p_{vap}^{sat} D} \right)^{-1}
combines the effects of thermal conductivity and water diffusivity.L_s
is the latent heat of sublimation,K_{thermo}
is the thermal conductivity of air,R_v
is the gas constant of water vapor,D_{vapor}
is the diffusivity of water vapor
Using eq. (\ref{eq:mass_rate}) and the assumed m(r)
relationship we obtain
Finally the snow autoconversion rate is computed as
!!! note We should include ventilation effects.
For non-spherical particles the mass rate of growth
should be multiplied by a function depending on the particle aspect ratio.
For functions proposed for different crystal habitats see
[Harrington1995](@cite) Appendix B.
We also have a simplified version of snow autoconversion rate, to be used in modeling configurations that don't allow supersaturation to be present in the computational domain. It is formulated similarly to the rain autoconversion:
where:
q_{liq}
- liquid water specific humidity,\tau_{acnv\_rain}
- timescale,q_{liq\_threshold}
- autoconversion threshold.
Accretion defines the rates of conversion between different categories due to collisions between particles.
For the case of collisions between cloud water (liquid water or ice) and precipitation (rain or snow) the sink of cloud water is defined as:
where:
c
subscript indicates cloud water category (cloud liquid water or ice)p
subscript indicates precipitation category (rain or snow)E_{cp}
is the collision efficiency.
Integrating over the distribution yields:
where:
\Pi_{a, v}^p = a_0^p \, v_0^p \, \chi_a^p \, \chi_v^p
\Sigma_{a, v}^p = a_e^p + v_e^p + \Delta_a^p + \Delta_v^p
For the case of cloud liquid water and rain and cloud ice and snow collisions,
the sink of cloud water becomes simply the source for precipitation.
For the case of cloud liquid water and snow collisions
for temperatures below freezing, the sink of cloud liquid water is
a source for snow.
For temperatures above freezing, the accreted cloud droplets
along with some melted snow are converted to rain.
In this case eq. (\ref{eq:accrfin}) describes the sink of cloud liquid water.
The sink of snow is proportional to the sink of cloud liquid water with
the coefficient \frac{c_{vl}}{L_f}(T - T_{freeze})
,
where c_{vl}
is the isochoric specific heat of liquid water,
L_f
is the latent heat of freezing,
and T_{freeze}
is the freezing temperature.
The collisions between cloud ice and rain create snow. The source of snow in this case is a sum of sinks from cloud ice and rain. The sink of cloud ice is defined by eq. (\ref{eq:accrfin}). The sink of rain is defined as:
where:
E_{ir}
is the collision efficiency between rain and cloud icen_i
andn_r
are the cloud ice and rain size distributionsm_r
,a_r
andv_{term}
are the mass(radius), cross section(radius) and terminal velocity(radius) relations for rainr_i
andr_r
mark integration over cloud ice and rain size distributions
Integrating eq.(\ref{eq:accr_ir}) yields:
where:
\Pi_{m, a, v}^{rai} = m_0^{rai} \, a_0^{rai} \, v_0^{rai} \, \chi_m^{rai} \, \chi_a^{rai} \, \chi_v^{rai}
\Sigma_{m, a, v}^{rai} = m_e^{rai} + a_e^{rai} + v_e^{rai} + \Delta_m^{rai} + \Delta_a^{rai} + \Delta_v^{rai}
Collisions between rain and snow result in snow in temperatures below freezing and in rain in temperatures above freezing. The source term is defined as:
where
i
stands for rain (T>T_{freezing}
) or snow (T<T_{freezing}
)j
stands for snow (T>T_{freezing}
) or rain (T<T_{freezing}
)a(r_i, r_j)
is the crossection of the two colliding particles
There are two additional assumptions that we make to integrate eq.(\ref{eq:accr_sr1}):
-
\left|v_{term}(r_i) - v_{term}(r_j)\right| \approx \left| v_{ti} - v_{tj} \right|
We approximate the terminal velocity difference for each particle pair with a difference between mass-weighted mean terminal velocities and move it outside of the integral. See the discussion in Ikawa_and_Saito_1991 page 88. -
We assume that
a(r_i, r_j) = \pi (r_i + r_j)^2
. This corresponds to a geometric formulation of the collision kernel, aka cylindrical formulation, see Wang2006 for discussion.
The eq.(\ref{eq:accr_sr1}) can then be integrated as:
!!! note Both of the assumptions needed to integrate the snow-rain accretion rate could be revisited:
The discussion on page 88 in
[Ikawa\_and\_Saito\_1991](https://www.mri-jma.go.jp/Publish/Technical/DATA/VOL_28/28_005.pdf)
suggests an alternative approximation of the velocity difference.
The ``(r_i + r_j)^2`` assumption for the crossection is inconsistent
with the snow crossection used when modelling collisions with cloud water
and cloud ice.
We start from a similar equation as when computing snow autoconversion rate
but integrate it from 0
to \infty
.
In contrast to eq.(\ref{eq:mass_rate}), now we are taking into account ventilation effects:
where:
F(r)
is the rain drop ventilation factor defined in (\ref{eq:ventil_factor})- saturation S is computed over water or ice
The final integral is:
For the case of rain we only consider evaporation (S - 1 < 0
).
For the case of snow we consider both the source term due to vapor deposition
on snow (S - 1 > 0
) and the sink due to vapor sublimation (S - 1 < 0
).
!!! note We should take into account the non-spherical snow shape. Modify the Reynolds number and growth equation.
If snow occurs in temperatures above freezing it is melting into rain. The sink for snow is parameterized again as
For snow melt
where:
F(r)
is the ventilation factor defined in (\ref{eq:ventil_factor})L_f
is the latent heat of freezing.
If T > T_{freeze}
:
The rain radar reflectivity factor (Z
) is used to measure the power returned by a radar signal when it encounters rain particles, and it is defined as the sixth moment of the rain particles distribution:
Integrating over the assumed Marshall-Palmer distribution (eq. 6) leads to
where:
n_{0}^{rai}
- rain drop size distribution parameter,\lambda
- as defined in eq. 7
By dividing Z
with the equivalent return of a 1 mm
drop in a volume of a meter cube (Z_0
) and applying the decimal logarithm to the result, we obtains the logarithmic rain radar reflectivity L_Z
, which is the variable that is commonly used to refer to the radar reflectivity values:
The resulting logarithmic dimensionless unit is decibel relative to Z
, or dBZ
.
include("plots/Microphysics1M_plots.jl")