1 引　言

云内微物理过程复杂，建立更好的云微物理参数化方案是提高数值模式预报准确性关键所在（胡志晋，1986；谭超等，2013）。数值模式中处理湿物理过程的方案主要包括陆面过程、边界层方案、积云对流参数化方案和云微物理方案（姜晓飞，2015）。目前天气预报模式的微物理过程多使用显式云方案（刘奇俊等，2003）描述。显式云方案的关键问题之一是如何解决格点尺度水汽凝结问题。传统的云微物理方案需要格点尺度水汽达到饱和情况下才启动凝结、凝华等微物理过程的计算（Zhao，et al，1997）。但全球模式水平分辨率较低，特别是在热带对流活跃区域，其对流的尺度很可能小于格点尺度，因此若继续沿用水汽达到饱和才启动凝结等过程，就会使得低纬度对流活跃的区域用微物理过程描述的水汽凝结过程偏弱。针对这一问题，目前多数全球模式云方案中均已引入大尺度宏观云方案，用以描述水汽的蒸发、凝结和凝华过程，并假定在格点相对湿度达到一定阈值的情况下便启动水汽凝结和凝华等过程（Lin，et al，1983；Sundqvist，1988；Reisner，et al，1998，Thompson，et al，2004），同时利用云微物理方案描述水成物的碰并、冻结、融化、沉降等其他微物理过程（Thompson，et al，2008）。另外，对流过程对大尺度云降水过程有重要影响，如何合理考虑次格点对流过程对云微物理过程的影响，增强云方案与对流参数化间的联系，也是全球模式云方案亟需解决的问题（Gettelman，et al，2008；Morrison，et al，2008）。

当前比较有代表性的全球模式云方案是欧洲中期天气预报中心（ECMWF）和美国环境预报中心（NCEP）的全球模式云方案以及GRAPES全球模式云方案。3种全球模式原理相似，均是利用大尺度宏观云方案描述水汽的凝结、凝华等过程，之后根据温度将水成物含量分成液水和冰水，用微物理方案描述碰并、凇附、自动转化等其他微物理过程。不同的是，ECMWF全球模式云方案和GRAPES全球模式云方案的云量采用的是云量显式预报方案，即云量作为模式预报量直接给出。而NCEP全球模式云方案的云量是诊断得到（Zhao，et al，1997），即根据格点平均水汽和格点预报的含水量拟合出的经验公式得到，该经验公式一般受模式分辨率影响较大，且缺乏足够的物理基础，对云量的预报性能稍差。此外，ECMWF全球模式云方案和GRAPES全球模式云方案均考虑了次格点对流卷出的影响，以期增加云方案与对流参数化的联系，而NCEP全球模式云方案并未考虑该影响。

GRAPES_GFS1.0版的云方案是借鉴NCEP全球模式云方案以及WSM6方案，2013年双参数云微物理方案（刘奇俊等，2003）加入到全球模式中，并且实现与大尺度宏观云方案、云量显式预报方案以及次格点对流卷出影响项的耦合（谭超等，2013），自此，GRAPES_GFS拥有了自主研发的云方案。之后，由于模式对低纬度地区云预报性能稍差，在低云预报方面做了改进（姜晓飞，2015），一定程度上改善了模式预报热带低层云水凝结量少的问题。

数值模式云方案改进的基础建立在对云方案的诊断评估上（吴伟，2011；马严枝等，2016；聂皓浩，2016），ECMWF曾经对自己的全球模式进行过诊断评估工作，并在此基础上对云方案做了改进优化。目前专门针对GRAPES_GFS云方案的诊断评估及改进工作开展相对较少，且GRAPES_GFS在云量预报及降水预报上仍有欠缺，热带地区格点降水量低，中低纬度地区云量预报量少等问题较为突出。因此，对GRAPES_GFS云方案进行诊断评估，并在此基础上对其进行改进是有科学意义和应用价值的。其中，云方案关于次格点对流卷出影响项的描述以及大尺度宏观云方案中关于水汽凝结等过程的描述，是改进中低纬度地区云量及降水预报的关键所在（姜晓飞，2015）。

文中在对GRAPES_GFS2.3版预报结果和ERA5分析评估的基础上，对GRAPES_GFS云方案进行优化改进。研究重点包括云量、云含水量、柱含水量等云宏、微观物理量以及降水、温度、向外长波辐射（OLR）、物理过程对温度和湿度的反馈倾向等物理量，通过与ERA5相应物理量对比分析，检验云方案改进是否合理。

2 GRAPES_GFS云微物理参数化方案的改进 2.1 GRAPES_GFS云方案简介

GRAPES全球模式云方案包括双参数云微物理方案、大尺度宏观云方案、云量显式预报方案和次格点对流卷出的影响过程4个主要部分，在此基础上，水成物含水量、粒子数浓度、云量以及云物理过程对温度影响的预报方程如下

$\begin{split}\frac{{\partial {q_x}}}{{\partial t}} =& A\left({{q_x}} \right) + S_{{\rm{micro}}}^{{q_x}} + S_{{\rm{macro}}}^{{q_x}} +\\& {S_{\rm{convd}}}\left({{q_x}} \right) + {S_{\rm{convs}}}\left({{q_x}} \right) + D\left({{q_x}} \right)\end{split}$

(1)

$\begin{split}\frac{{\partial {N_x}}}{{\partial t}} =& A\left({{N_x}} \right) + S_{{\rm{micro}}}^{{N_x}} + S_{{\rm{macro}}}^{{N_x}} +\\& {S_{\rm{convd}}}\left({{N_x}} \right) + {S_{\rm{convs}}}\left({{N_x}} \right) + D\left({{N_x}} \right)\end{split}$

(2)

$\frac{{\partial a}}{{\partial t}} = A\left(a \right) + S_{{\rm{strat}}}^a + S_{{\rm{convd}}}^a + S_{{\rm{convs}}}^a + D\left(a \right)$

(3)

$\frac{{\partial T}}{{\partial t}} = A\left(T \right) + S_{\rm{micro}}^{ T} + S_{\rm{macro}}^{ T} + S_{\rm{convd}}^{ T} + S_{\rm{convs}}^{ T} + D\left(T \right)$

(4)

式（1）为水成物含量（ $ q $ ）的预报方程，其中 $ {q}_{{{x}}} $ 包括水汽（ $ {q}_{{\rm{v}}} $ ）以及云水（ $ {q}_{{\rm{c}}} $ ）、雨水（ $ {q}_{{\rm{r}}} $ ）、云冰（ $ {q}_{{\rm{i}}} $ ）、雪（ $ {q}_{{\rm{s}}} $ ）四种水成物的含水量，等式右边各项分别表示平流过程、微物理过程、大尺度宏观云过程、次格点深对流过程、次格点浅对流过程和扩散过程引起的水成物含量变化。式（2）为水成物数浓度（ $ N $ ）预报方程，其中 $ {N}_{x} $ 包括雨水（ $ {N}_{{\rm{r}}} $ ）、云冰（ $ {N}_{{\rm{i}}} $ ）、雪（ $ {N}_{{\rm{s}}} $ ）三种水成物的数浓度，等式右边对应的各项与式（1）相同。式（3）为云量（ $ a $ ）的预报方程，等式右边分别为平流引起云量变化项、层积云云量变化、次格点深对流过程、次格点浅对流过程和扩散过程对云量的影响。式（4）为湿过程引起的温度变化方程，等式右边分别为平流过程、大尺度宏观云过程、次格点深对流过程、次格点浅对流过程及扩散过程对温度影响。

2.2 GRAPES_GFS云方案的优化讨论

在对GRAPES_GFS云方案原理解读、脚本分析以及预报结果诊断评估的基础上，对GRAPES_GFS云方案做了以下优化（表1）。

云方案具体改进方案包括：

（1）GRAPES_GFS原云方案中，计算水汽凝结、凝华时，需要用到权重饱和比湿（ $q_{\rm{vs\_weight}}$ ），来计算云区内云水和冰水凝结量。原方案中， $q_{\rm{vs\_weight}}$ 计算公式为

$ {q_{\rm{vs\_weight}}} = \left\{ {\begin{array}{*{20}{c}} {{q_{{\rm{vs}}}}}\quad\quad\quad\quad\quad\quad\quad\quad\\ {\left({1 - \mu } \right) \times {q_{{\rm{vs}}}} + \mu \times {q_{{\rm{vsi}}}}} \end{array}} \right.\begin{array}{*{20}{c}} {T {\text{≥}} {0^ \circ }{\rm{C}}}\\ {T {\text{＜}} {0^ \circ }{\rm{C}}} \end{array} $

(5)

$ \mu =\left({q}_{{\rm{i}}}+{q}_{{\rm{s}}}\right)/q $

(6)

考虑到模式计算云微物理项稳定性不高，容易造成水汽凝结、凝华量计算误差较大，故改用更稳定的温度作为约束条件，计算权重饱和比湿，公式为

$ {q_{\rm{vs\_weight}}} = \left\{ {\begin{array}{*{20}{c}} \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!{{q_{\rm{vs}}}}&{\!\!T {\text{≥}} 273.15\;{\rm{K}}\qquad}\\ {\alpha \times {q_{\rm{vs}}} + \left({1 - \alpha } \right) \times {q_{\rm{vsi}}}}&{{T_{\rm{ice}}} {\text{≤}} T {\text{＜}} 273.15\;{\rm{K}}}\\ \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!{{q_{\rm{vsi}}}}&{\!\!\!\!\!T {\text{＜}} {T_{\rm{ice}}}\qquad\quad\quad} \end{array}} \right. $

(7)

$ \alpha ={\left[\left(T-{T}_{\rm{ice}}\right)/\left(273.15-{T}_{\rm{ice}}\right)\right]}^{2} $

(8)

式中， $ {q}_{\rm{vs}} $ 为相对水面饱和比湿， $ {q}_{\rm{vsi}} $ 为相对冰面饱和比湿， $ q={q}_{{\rm{i}}}+{q}_{{\rm{s}}}+{q}_{{\rm{c}}}+{q}_{{\rm{r}}} $ ， $T_{\rm{ice}}=250.16 \;{\rm{K}}$

（2）GRAPES_GFS原云方案中，在描述与水汽有关的相变（水汽凝结、凝华等）时，考虑云量影响。为保证云方案改进前、后对云量影响的处理一致，在原方案的基础上，增加了云量对自动转化、碰并、融化等其他微物理过程的影响。

（3）云方案中，云量由大尺度云降水过程和对流过程决定，其中，格点尺度云量由水成物的含水量变率计算得到。

云方案中格点尺度云量计算，即在原云量计算公式（谭超等，2013）基础上，调整系数 $ \beta $ ，经多次敏感性试验后，确定 $ \beta $ 合适值，使得云量模拟更合理。

具体如下

$ {a_{\rm{strat}}} = \beta \cdot \frac{{{{\left({1 - a} \right)}^2}}}{{\left({{q_{\rm{v}}} - q_{\rm{vs}}} \right)}} \cdot \frac{{{\rm d}q}}{{{\rm d}t}} $

(9)

$ \beta = \left\{ {\begin{array}{*{20}{l}} {40}&{\qquad k {\text{＜}} 15} & {\left( {\text{低云}} \right)}\\ {10}&{15 {\text{≤}} k {\text{＜}} 29}&{\left( {\text{中云}} \right)}\\ {25}&{\qquad k {\text{≥}} 29}&{\left( {\text{高云}} \right)} \end{array}} \right. $

式中， $ {a}_{\rm{strat}} $ 为云量生成项， $q_{\rm{vs}}$ 为格点尺度水汽饱和比湿， $ {\rm d}q/{\rm d}t $ 为水成物的含水量变率， $ k $ 为模式层层数。

此外，优化了云量计算的部分约束条件，如原方案中，假定高、中、低云凝结生成量最高均不超过0.6，根据敏感性试验结果，新方案将其细化，分别为0.6、0.8、0.8。

（4）调整积云对流参数化对云方案的影响项计算公式（谭超等，2013）。

对流卷出对格点尺度云的影响分云量、液水和冰水3项，反馈到云方案时，上述3项又分为深对流项和浅对流项。计算公式为

${S_{{\rm{conv}}}} = {D_{{\rm{up}}}}{l_{{\rm{up}}}} + \frac{{{M_{\rm{up}}}}}{\rho }\frac{{\partial {q_x}}}{{\partial z}}$

(10)

$\delta {a_{\rm{conv}}} = {D_{{\rm{up}}}} + \frac{{{M_{\rm{up}}}}}{\rho }\frac{{\partial a}}{{\partial z}}$

(11)

${D_{{\rm{up}}}} = {\delta _{\rm{up}}}\frac{{{M_{\rm{up}}}}}{{\bar \rho }}$

(12)

${l_{{\rm{up}}}} = \frac{{{q_x}}}{a}$

(13)

式中， $ {D}_{\rm{up}} $ 为积云上升气流卷出率， $ {M}_{\rm{up}} $ 为上升气流质量通量， $ {l}_{\rm{up}} $ 为积云云水或云冰含水量， $ {q}_{x} $ 为格点平均含水量， $ {\delta }_{\rm{up}}=0.75\times 1{0}^{-4}\;{\rm m}^{-1} $ 。

原方案中，并未将式（10）中的 $ {D}_{\rm{up}}{l}_{x} $ 项及式（11）中的 $ {D}_{\rm{up}} $ 项加入到次格点对流对云方案的影响项中，新方案考虑了这两项的影响。

（5）优化和改进次格点总云水凝结率的计算。

原云方案中，宏观云水汽、温度及水成物控制方程（Zhao，et al，1997）为

$\frac{{\partial q}}{{\partial t}} = {q_{\rm{non}}} + {E_{\rm{c}}} + {E_{\rm{r}}} - {C_{\rm{b}}} - {C_{\rm{g}}}$

(14)

$\frac{{\partial T}}{{\partial t}} = {T_{\rm{non}}} - \frac{L}{{{c_{{p}}}}}\left({{E_{\rm{c}}} + {E_{\rm{r}}} - {C_{\rm{b}}} - {C_{\rm{g}}}} \right) - \frac{{{L_{\rm{f}}}}}{{{c_{{p}}}}}{P_{\rm{sm}}}$

(15)

$\frac{{\partial Q}}{{\partial t}} = {Q_{\rm{non}}} - {E_{\rm{c}}} + {C_{\rm{b}}} + {C_{\rm{g}}} - P$

(16)

式中，T为温度，q为水汽混合比，Q为云水/冰混合比。 ${q}_{\rm{non}}{\text{、}}\;{T}_{\rm{non}}{\text{、}}\;{Q}_{\rm{non}}$ 为非凝结项（平流和湍流项之和）， ${C}_{\rm{{g}}}{\text{、}}\;{C}_{\rm{{b}}}$ 为格点尺度和对流尺度的凝结率， ${E}_{{\rm{c}}}{\text{、}}\;{E}_{{\rm{r}}}$ 是格点平均总云水和总降水的蒸发率， $ {P}_{\rm{sm}} $ 是融化层以下雪的融化率， $ P $ 为来自云水（冰）混合比的降水生成率， ${c}_{{{p}}}$ 为定压空气比热容， $ L $ 是凝结和凝华的潜热系数， $ {L}_{{\rm{f}}} $ 是冻结潜热系数。

在格点尺度下，上述公式可简化为

$ \frac{\partial q}{\partial t}={A}_{{\rm{q}}}+{E}_{{\rm{c}}}-{C}_{{\rm{g}}} $

(17)

$ \frac{\partial T}{\partial t}={A}_{{\rm{t}}}-\frac{L}{{c}_{{{p}}}}\left({E}_{{\rm{c}}}-{C}_{{\rm{g}}}\right) $

(18)

$ \text{其中}\quad\quad\quad{A}_{{\rm{q}}}={q}_{\rm{non}}+{E}_{{\rm{r}}}-{C}_{{\rm{b}}}\quad\quad\quad\quad $

(19)

$ \quad\quad\quad\quad\quad{A}_{{\rm{t}}}={T}_{\rm{non}}-\frac{L}{{c}_{{{p}}}}\left({E}_{{\rm{r}}}-{C}_{\rm{{b}}}\right)-\frac{{L}_{{\rm{f}}}}{{c}_{{{p}}}}{P}_{\rm{sm}} $

(20)

凝结率 $ {C}_{g} $ 计算公式为

$ {C}_{{\rm{g}}}=\frac{M-{q}_{{\rm{s}}}{f}_{{\rm{t}}}}{1+f\varepsilon {L}^{2}{q}_{{\rm{s}}}/\left(R{c}_{p}{T}^{2}\right)} $

(21)

$ f_{\rm{t}}=\frac{{\rm{2}}\left({\rm{1}}-{{a}}\right)\left(f_{\rm{s}}-f_{\rm{0}}\right) \left[\left(1-a\right)M+E_{\rm{C}}\right]}{{\rm{2}}q_{\rm{s}}\left({\rm{1}}-{{a}}\right)\left(f_{\rm{s}}-f_{\rm{0}}\right)+Q/{{a}}} $

(22)

式中，f_t为相对湿度的时间变率， ${f_{\rm{0}}} $ 为凝结蒸发的相对湿度阈值， ${f_{\rm{s}}}=1$ ， $\varepsilon $ =0.622， $q_{\rm{s}}$ = $q_{\rm{vs\_weight}} $

$ M={A}_{{\rm{q}}}-\frac{f\varepsilon L{q}_{{\rm{s}}}}{R{T}^{2}}{A}_{{\rm{t}}}+\frac{f{q}_{{\rm{s}}}}{P}\cdot \frac{\partial P}{\partial t} $

(23)

式（23）中最后一项，NCEP全球模式云方案将其作为一个小项，模拟时忽略不计。GRAPES_GFS原云方案根据理想气体状态方程，将该项近似为 $\dfrac{f{q}_{{\rm{s}}}}{P}{\text •} \rho R\dfrac{\partial T}{\partial t}$ ，用微物理相变对温度的影响确定该项。考虑到温度变化不单是水汽相变决定，新方案改用A_t来表征该项，其中A_t代表平流、辐射、对流、边界层过程和扩散的总云水温度和的倾向。

新方案中，大尺度宏观云计算凝结等过程时，M计算公式简化为

$ M={A}_{{\rm{q}}}-\frac{f{q}_{{\rm{s}}}}{T}{\text •} \left(\frac{\varepsilon L}{RT}-1\right){A}_{{\rm{t}}} $

(24)

（6）考虑到云滴尺度小于雨滴，蒸发率远大于雨滴蒸发率，故新方案假定若格点内云水和雨水同时存在，雨水不蒸发，只有在格点内无云存在或者云水蒸发完后，雨水蒸发才启动。

3 试验结果与资料分析

本研究使用GRAPES_GFS2.3版，模式水平分辨率为0.25°×0.25°，垂直共分60层，积分步长取300 s，选取2009年6月1—5日的一次天气过程作为试验个例，预报时间共计5 d，起报时间为1日00时（世界时，下同），结束时间为6日00时，输出时间间隔为6 h，共计20个时次预报结果。同时选取2017年7月1—31日资料，进行为期1个月的批量试验。在对试验预报结果与ERA5再分析资料对比分析基础上，对云方案进行改进。

ERA5是ECMWF第5代全球气候监测数据集（张永林等，2020），水平分辨率为31 km，相较于ERA-Interim，ERA5时间分辨率大幅度提高，并提供了特定天气事件更为详细的演变过程。本研究主要使用ERA5的云量、降水及水成物含量等相关数据。

3.1 原方案与新方案云水和云冰含水量预报对比分析

GRAPES_GFS与云相关预报结果主要包括液相水成物和冰相水成物含水量。本研究通过对云方案改进前、后各水成物柱含量、垂直分布进行检验分析，参照ERA5结果，对云方案改进结果进行分析验证。

图1分别给出了ERA5以及GRAPES_GFS云方案改进前、后液相和冰相水成物柱含量预报为2009年6月1—5日20个预报时次的平均值。

对比图1a、b，可以看到GRAPES_GFS云方案改进前，中高纬度地区液水含量明显高于ERA5，低纬度地区大部分区域柱液水含量低于ERA5，ERA5柱液水含量在赤道太平洋洋面上有明显的带状分布，GRAPES_GFS柱液水含量在印度尼西亚地区明显少于ERA5。根据图1c，云方案改进后GRAPES_GFS赤道地区柱液水含量升高，特别是在赤道太平洋海域柱液水含量出现与ERA5类似的带状分布，大洋洲北部和印度尼西亚地区液水含量也明显增大，柱液水总体分布情况和ERA5分布更接近。对比图1d、e，可以看出云方案改进前GRAPES_GFS绝大部分区域冰水含量明显高于ERA5，只有赤道附近太平洋洋面上，ERA5柱冰水含量相对集中且高于GRAPES_GFS。根据图1f，云方案改进后赤道附近柱冰水含量略有增加，中高纬度地区改进不明显，需对冰水进一步优化。

为分析含水量垂直分布情况，图2给出个例试验期间ERA5以及GRAPES_GFS云方案改进前、后水成物总量经向平均后的垂直剖面。对比图2a、b、c可见，在纬度40°S—0°，高度900—700 hPa范围内，云方案改进前，总含水量明显低于ERA5，改进后该区域含水量增大，呈现连续分布。在中高纬度地区，尤其是南北纬60°附近，GRAPES_GFS总含水量高于ERA5，但在北半球热带区域，ERA5总含水量分布高度最高可达100 hPa，GRAPES_GFS水成物最高分布大概在130 hPa。总体来看，GRAPES_GFS云方案的改进提高了热带地区低层含水量预报性能，对低云预报准确度的提高有帮助。

3.2 原方案与新方案云量预报对比分析

云量预报主要包括高、中、低云量和总云量。云量与周围环境温度、湿度、含水量以及空气运动等有密切关系。

图3是ERA5及GRAPES-GFS云方案改进前、后个例试验期间总云量和低云量全球平均分布，对比图3a、b、c，ERA5和GRAPES_GFS的总云量全球分布相似，但ERA5总云量相对较大，北极极地和赤道附近更为明显，ERA5云量分布更为集中，小而散的云块较少，这和含水量预报结果类似。GRAPES_GFS云方案改进后，对低层含水量影响较大，热带地区极为明显，而低云量可以较好反映低层含水量变化。从图3d、e、f可以看出，低云量分布与总云量分布相似，ERA5热带地区云区成片分布，特别是太平洋上空，但云量较小，多数区域小于0.3。GRAPES_GFS云方案改进前，热带低云量少于0.1的区域多于ERA5，云量总体来看小于ERA5，云方案改进后，热带低云量略有增加，与低云含水量增加一致。

GRAPES_GFS预报低云时，在南极洲出现空白区域，该区域海拔较高，推测该现象可能是地形因素造成的。总体来看，调整浅对流卷入对云方案影响过程，增加云量对微物理过程的影响，对低云影响较大。这是由于次格点浅对流卷入在低层贡献较大，水成物含量变化集中在低层，使得对流活动旺盛的热带地区液水含量增大，对流卷出的云水蒸发导致格点湿度升高，对云量预报结果改善产生正效果。

3.3 原方案与新方案降水预报对比分析

GRAPES_GFS降水预报包括次格点尺度降水和格点尺度降水两部分。图4给出个例试验期间ERA5及GRAPES_GFS云方案改进前、后总降水量日平均分布及云方案改进前、后的差值，分析图4a、b、c，GRAPES_GFS和ERA5日降水区域分布相差不大，ERA5总降水量高于GRAPES_GFS，热带地区尤为明显，特别是靠近赤道地区沿纬向分布的降水带，在中纬度地区，ERA5日降水量略高于GRAPES_GFS。结合图4d总降水差值结果得到，云方案改进后，总降水量略有增大，热带地区增加明显，尤其是南美洲北部地区以及赤道纬向降水带区域，低纬度地区降水量增加1—4 mm。

图5是ERA5及GRAPES_GFS云方案改进前、后次格点尺度降水和格点尺度降水的日平均预报。对比图5a、b、c可见，ERA5次格点尺度降水明显高于GRAPES_GFS次格点尺度降水，特别是印度尼西亚地区以及中纬度区域尤为明显，云方案改进后，热带次格点尺度降水明显增大，中纬度降水量也稍有增加。根据图5d、e、f，GRAPES_GFS格点尺度降水集中在中纬度地区，热带很少，ERA5在热带地区格点尺度降水明显高于GRAPES_GFS格点尺度降水。云方案改进后，热带地区格点尺度降水有所增加。总体来看，GRAPES_GFS无论次格点尺度降水还是格点尺度降水都仍小于ERA5。

3.4 GRAPES_GFS云方案改进前、后辐射预报对比分析

辐射和云量密切相关，云预报结果好坏直接影响辐射预报质量。图6给出个例试验期间20个预报时次向外长波辐射和高云量小时平均预报全球分布，对比图6a、b、c可知，ERA5的向外长波辐射值明显高于GRAPRS_GFS，尤其是北半球高纬度地区和热带地区，约高15 W/m²，但总体分布相似。云方案改进后，热带地区云量增加，向外长波辐射减小，在赤道太平洋洋面及印度尼西亚地区表现较为明显。从图6d、e、f可以看出，ERA5在印度尼西亚区域高云云量较大，极区高云量稍小，向外长波辐射则在赤道附近小于GRAPRS_GFS。总体来看，云方案的改进，对于向外长波辐射的改进效果并不明显，仅在热带部分地区（印度尼西亚等地）略有效果。

3.5 GRAPES_GFS云方案改进前后温度、含水量倾向项分析

辐射过程、对流过程及微物理过程等都对温度或水汽含水量有重要反馈作用，在模式中，诊断GRAPES_GFS中辐射过程、对流参数化过程以及云微物理过程对温度和水汽的反馈倾向项，分析模式温度场、湿度场对这些物理过程的反馈结果，是优化和改进模式物理方案的重要手段。

图7a为GRAPES_GFS云方案改进前热带辐射过程、浅对流过程、深对流过程及云微物理过程温度倾向项的预报。对流过程对周围环境的影响主要是升温降湿，根据图7a，深、浅对流过程温度倾向项均为正值，深对流深厚，在100 hPa高度以下，对温度都有影响，在700—300 hPa高度最为明显；浅对流浅薄，对温度影响主要集中在1000—700 hPa。云微物理过程在700 hPa高度以下的低层，使温度降低，这和低层云滴、雨滴等的蒸发吸热有关，高层温度倾向项为正值，这和水汽凝结或凝华，云滴、雨滴冻结等物理过程放热有关。可见，温度场对对流和微物理过程都产生合理反馈。图7b、c表示深对流、浅对流过程的水成物含水量倾向项和云微物理过程的水汽、液水、冰水含量倾向项的经向平均预报结果。可以看出，对流过程对湿度场的影响主要表现为减湿效果。深对流使水汽减少，尤其是低层，浅对流在900 hPa高度以下其水汽含量倾向项为负值，900—700 hPa区域为正值，但很小，其上受浅对流尺度限制，浅对流的水汽含量倾向项趋近于0。从云微物理过程含水量倾向项变化可以看出，当液水倾向项、冰水倾向项为正值时，水汽倾向项为负值，表征水汽与液（冰）水的相互转化关系。水汽倾向项多为正值，峰值约出现在850 hPa和400 hPa高度，达8×10⁻⁴ g/（kg·d）。液水在高层碰并成冰，倾向项为负值，低层积聚成雨滴下落，倾向项为正值，冰水倾向项正值主要集中在150—500 hPa高度，这与成冰过程主要集中在高层有关。从GRAPES_GFS云方案改进前、后含水量倾向项差值（图7d）可以看出，云方案改进使深对流水汽含量倾向项增大，加大了模式的去湿过程。浅对流含水量倾向项差值变化较小。而云方案水汽倾向项略有减小，液水倾向项变化不大，冰水倾向项增大。

3.6 GRAPES_GFS云方案改进前、后批量试验预报

在个例试验的基础上，进行2017年7月1—31日为期1个月的批量试验。分析批量试验期间云方案改进前、后柱液水含量和柱冰水含量差值（图8）得出，云方案改进后，中低纬度柱液水含量略有增加，增值30—50 g/m²，柱冰水含量在低纬度和高纬度地区增加较为明显，中纬度地区有增有减。总体来看，云方案改进后，柱液水含量的改变较符合预期，柱冰水含量的改进效果不明显。

批量试验期间ERA5及GRAPES_GFS云方案改进前、改进后次格点尺度降水量、格点尺度降水量、总降水量日平均分布结果（图9）显示，GRAPES_GFS和ERA5降水区域分布差别不大，对比图9a、b、c可以看到，GRAPES_GFS日平均次格点尺度降水量略低于ERA5，特别是在赤道附近区域，次格点尺度降水低于ERA5，但ERA5日平均次格点尺度降水量0.1—1 mm的小降水量级区域分布多于GRAPES_GFS，这与GRAPES_GFS次格点尺度降水不考滤浅对流过程形成的降水有直接关系。GRAPES_GFS云方案改进后热带区域次格点尺度降水略有增加，特别是印度尼西亚附近增加更明显。GRAPES_GFS格点尺度降水低于ERA5，在热带地区表现更为明显，云方案改进后，格点尺度降水增加，中、低纬度区域都有明显改进。总降水量预报结果显示，云方案改进前，GRAPES_GFS降水量低于ERA5，格点尺度降水在热带区域偏小。云方案改进后，GRAPES_GFS总降水量增加，与EAR5总降水量分布更接近，特别是在热带区域格点尺度降水的预报，改进效果明显。

4 总结与讨论

基于GRAPES_GFS2.3版，对云方案中对流卷入影响过程、宏观云过程进行优化，同时增加云量对碰并、融化等微物理过程的影响，通过个例分析，对云方案改进前、后各水成物含水量、云量、降水、向外长波辐射等预报结果进行分析，以ERA5资料结果作为参照，得到以下结论：

（1）将云方案中计算水汽凝结时用到的权重饱和比湿改为由温度确定后，GRAPES_GFS液相水成物在热带含量增多较为明显，且低云预报性能有明显改进。

（2）增加云量对云微物理过程的影响，云量、含水量及降水预报准确度均有提高，特别是在热带地区，改进效果更明显。具体改进效果主要表现为云方案改进后，热带区域柱水含量出现“带状”分布。

（3）优化格点云量计算过程，调整高云、中云、低云生成过程中部分约束条件之后，GRAPES_GFS云量明显增多，特别是中、低纬度地区，云量预报准确度明显提高。

（4）调整宏观云方案中次格点总云水凝结率和雨水蒸发，在一定程度上增加了水成物含量预报量，使得GRAPES_GFS含水量预报结果与ERA5相近，含水量垂直分布特征与ERA5含水量垂直分布相似，特别是低层含水量的垂直分布，同时，缓解了GRAPES_GFS热带地区格点尺度降水过少现象。

（5）GRAPES_GFS双参数云微物理方案改进之后，含水量、云量、降水预报都有明显改进，这也有助于提高模式中其他与云有关过程预报质量，如向外长波辐射的预报。

本研究虽然对云方案做了一些改进，但是GRAPES_GFS冰水含量预报偏高，中纬度区域降水偏少等问题仍存在，未来需要进一步研究完善。且全球范围的云微物理观测资料相对缺乏，后期可增加GRAPES预报结果与实测资料的比较。