大洋环流和海气相互作用的数值模拟：ch8_海洋模式中的物理参数化过程.pptVIP

* * * * * 冬季深、夏季浅 * * * * * * * * * * * * * * * * * * 小结（2） 叶绿素对于表层海温的影响是热量收支中多个过程平衡的结果。在赤道太平洋减少穿透深度，表层海温却变冷，主要是垂直平流冷平流超过其加热作用所致。 叶绿素浓度与短波穿透的关系是水圈和生物圈相互作用的一种重要机制。 The End of Chapter 8 Appendix A1 深对流 A2 KrausTurner的整体混合层方案 A1 深对流 海水温度方程 表示深对流， 没有具体的表达式 对流区和热盐环流输送带 对流区 深对流的特点 对流主要是由静力不稳定造成的，导致静力不稳定有表面的冷却以及海冰的盐析过程。 对流主要发生在高纬度的几个特定的区域。对流发生的空间尺度小，是现有的模式分辨率分辨不出来的。 现有的模式方程采用静力近似，也就是当有静力不稳定发生时不能产生垂直速度，故方程无法描述这种过程。 深对流参数化方案(1) 对流调整方案：在每一步积分之后，判断是否有不稳定发生。如果有就对相邻格点的温度和盐度进行混合，并在整个水柱中进行调整。 问题：不能保证整个水柱稳定，计算量大。 T1 T0 T T 深对流参数化方案(2) 增加垂直混合系数K。 优点：便于计算，渐进变化。 问题：不稳定存在时间较长，大垂直扩散系数会限制时间步长。 K~1-100 m2/s 小结 对流发生范围小，过程是非静力的。 两种主要的参数化方案：对流调整方案和增加垂直扩散系数方案。 * 2007-11 A4 KrausTurner的整体混合层方案 整体混合层方案（1）——描述海洋混合层基本方程 温度方程 盐度方程 We0时， 整体混合层方案（2） ——Kraus-Turner积分方程组 We为挾卷速度， 当混合强、混合层加深时，We0 当混合弱、混合层变浅时，We＝0 * 2007-11 整体混合层方案（3） ——挾卷速度的计算 上边界 底边界 湍流动能（TKE）方程 垂直积分TKE方程得: 温度方程, 湍流扩散, 短波穿透 K为垂直扩散系数 * * * * * * * * * * * * Transpire出现 dispersion色散、分散 * * * * * 3 中尺度涡和 GM90 方案 * 海水温度方程 （）*代表涡致速度 纬向平均的年平均位密度分布 AI = 103m2s-1 AD=10-4m2s-1 海洋内部沿等密度面的混合远大于穿过等密度面的混合 Veronis效应 夸大穿越等密度面的扩散 造成主温跃层偏深、大西洋翻转环流和向极热输送偏弱 坐标面 等密度面 Observed T(y,z) L204X5, z-混合 IPCC第一次评估报告： ·海洋模式模拟长时间尺度现象的能力在很大程度上依赖于次网格尺度混合参数化的精度； ·由于混合过程失真，许多粗网格模式模拟的永久性温跃层比观测结果偏暖、温度梯度明显偏弱. (Gates et al., 1990). 问题：如何更好地描述沿着等密度面的扩散？ 海洋 中尺度涡 湍流 分子运动 大尺度环流 长度尺度 能谱 次网格尺度 海洋动能谱, 次网格尺度过程 10km 1cm 海洋模式分辨率100km 中尺度涡是次网格尺度的主要运动形式，也是参数化的主要对象 Redi（1982）旋转坐标系 等z面 等密度面 等密度面上的对称混合张量： Redi（1982）旋转坐标系 h?—两个等位密度面之间的厚度， 采用雷诺平均的方法，带入方程得： ? -坐标下的连续方程和位温方程 使用GM90混合方案的温度方程 R是一个旋转算子，将等密度坐标下的混合转换到z坐标中 GM90中的涡致速度 其中k 为 “厚度扩散系数”，k 通常也取为103 m2 s-1 左右。这一方案是根据等密度面坐标中等密度面厚度沿梯度扩散的过程，这也是“厚度扩散系数”这一名词的由来。 由3X3 的矩阵给定了一个初始密度场来刻画y-z 平面上二维的分布 水平混合 3维混合 GM90 GM90方案的物理图像 GM90方案的理想试验结果 Gent et al.(1995) T=0 T=20 T=1000 Observed T(y,z) L30T63, GM90 IPCC第二次评估报告： Gent and McWilliams 等发展的一种同中尺度涡旋有关的混合参数化方法(即“GM90”)，取得了很有希望的结果: 模拟的温跃层变陡、深海变冷、经向热输送增加. (Danabasoglu et al., 1994) 小结 主温跃层混合中沿着等密度面的混合要远大于穿过等密度面的混合。 Z坐标模式中的Veronis效应。 中尺度涡混合参数化须在等位密度坐标下进行，而后再逐点转换到z坐标方程中。 GM90引入了涡致速度 u*,v*,