气候系统模式会_南京_李煜斌研究.pptx

数值模式中表面层通量非迭代方案的改进 李煜斌1，高志球1,2 1 南京信息工程大学 2 中国科学院大气物理研究所 江苏 南京 2014年12月8日 2014年度气候模式发展研讨会 Outline 背景 目前常用的方案 新方案 结果对比 结论 背景：表面层通量的计算 Businger-Dyer 方程 (Businger, 1966; Dyer, 1967). BD 方程是继 之后的又一开创性工作，但是这一方程在稳定情况下低估通量，当理查森数超过临界值（约0.2）时BD方程不适用 (Louis, 1979)。 Holtslag and De Bruin (1988) 以及 Beljaars and Holtslag (1991)提出了新的适用于稳定情况下的方程。 Cheng and Brutsaert (2005, CB05) 进一步修正了稳定情况下的稳定度修正方程，后人的研究证实其计算结果更准确(Guo and Zhang, 2007; Jimenez et al, 2012). 但是这些方案都需要迭代 Loius (1979), Garratt (1992), Launiainen (1995), Song (1998), De Bruinet al. (2000), Yang et al. (2001); Li et al. (2010) 提出了非迭代的方案，但是这些方程不适用于完整的粗糙度区间 -0.5≤log(z0/z0h)≤30 和10 ≤z/z0≤105 Wouters et al. (2012, WRL12) 是目前较新的一个非迭代方案，它适用于整个粗糙度区间 0.5≤log(z0/z0h)≤30 and 10 ≤z/z0≤105 , 然而它的计算结果与迭代方案有所偏差. Sharan and Srivastava (2014, SS14) 在 WRL12的基础上，对不稳定情况下的非迭代方案进行了改进， 目标：非迭代，整个粗糙度区间，结果准确 Monin-Obukhov相似理论 背景 表面层通量计算： Momentum flux: Sensible heat flux： 目前常用的迭代方案: L (i.e., ζ) 是它自己的函数，需要迭代求解 稳定情况下：CB05 不稳定情况下：B-D方程 CB05方案迭代至5%误差范围内所需的步数 在天气和气候模式中，一般采取只迭代5步即跳出循环 kB-1 = log(z0/z0h) 稳定情况下： 目前稳定层结下的非迭代方案: WRL12 用RiB计算 ζ ，无需迭代 但是对于整个z0 和 z0h 区间只用一个方程来计算，导致误差 某些z0 和 z0h 情况下 5步迭代法的误差 某些z0 和 z0h 情况下 WRL12方案的误差 WRL12 通过 ζt 将 ζ 分成两段，然后分段回归。 我们是否可以同样将z0 和 z0h分段？ CB05 ζ 精确解 与 RiB, log(z/z0) and kB-1 之间的关系 黑色粗线是三次方程的拟合线 Region z/z0 z0/z0h 1 10~160 0.607~100 2 160~105 0.607~100 3 10~80 100~107 4 80~105 100~107 5 10~40 107~1011 6 40~105 107~1011 7 10~40 1011~1.07×1013 8 40~105 1011~1.07×1013 对于每个区间，通过下式进行回归： 通过不同的尝试 ,我们将 z/z0和 z0/z0h 分成8个区间: 新方案 1) 通过z/z0 and z0/z0h; 找到所属区间 2) 通过 RiB 和方程 找到所属RiB分段 3) 使用  计算 ζ . Cmn and Cijk 为通过多元回归分析得到的系数. 新方案 已发表于GMD 某些z0 和 z0h 情况下 新方案的误差 某些z0 和 z0h 情况下 5步迭代法的误差 某些z0 和 z0h 情况下 WRL12方案的误差 结果对比 结果对比（稳定情况下） 计算方案 遍历整个区间所需计算时间 最大 ζ 误差 平均 ζ 误差 特性和建议 CB05 方案 6260 s N/A N/A 目前的最优解，但是耗费大量计算时间。如果不考虑计算时间那么该方案最优。 CB05五步迭代法 3960 s 超过 50% 超过15% 花费相对较少的计算时间，但是增大了计算结果的不确定性。 Wouters12方案 261 s 超过50% 超过15% 花费极少的计算时间，但是增大了计算结果的不