第五讲海洋模式中的参数化过程.ppt

1、第五讲 海洋模式中的参数化过程,刘海龙 2007年10月,2007年10月,2,模式动力框架和参数化过程,动力框架（数值方法） 方程、网格、差分格式、积分方案等等 参数化（物理过程） 湍流过程、中尺度涡、辐射传输等等,2007年10月,3,参数化对模式的重要性,有关近似和参数化的相关概念是贯穿我们整个阐述过程的主题。正如我们所强调的那样，按照惯例用于解决海洋行为的运动方程通过复杂的一系列运动学近似、物理的参数化和数值假定得到。任何一个或所有这些近似和参数化都可能对海洋模拟的质量产生重要影响。因此对于从事海洋环流模拟研究的新手来说，知道解的敏感性和潜在问题的源头是十分重要的。,（引自Numeri

2、cal ocean circulation modeling, by D. B. Haidvogel and A. Beckmann, 1999）,2007年10月,4,什么是参数化（Parameterization）？,to express in terms of parameters (Merriam-Websters Collegiate Dictionary) 在数值模式中，不考虑过程的细节，而是用其它一些确定的变量所表示的简化函数表示这个过程，这个方法称为“参数化”。 (Glossary of PO and Related Disciplines),2007年10月,5,例子：参数化

3、次网格湍流,表示密度，()表示扰动量，-表示Reynolds平均,Reynolds应力,单位时间内，单位面积，在z方向上输送的x方向的脉动动量通量的平均值。,K表示湍流粘性系数,海洋环流模式中哪些过程需要参数化？,2007年10月,7,纬圈和全球平均的海温分布,Thermocline 温跃层,Mixed layer 混合层,Tropical Thermocline,Permanent Thermocline,深对流,深对流,底边界层,短波辐射穿透,Temp,2007年10月,8,海水温度方程,深对流,水平混合,2007年10月,9,动量方程,2007年10月,10,海洋模式中的主要参数化过程,

4、动量方程中的水平和垂直粘性 温度和盐度方程中的水平和垂直混合 中尺度涡的混合和输运 深对流过程 短波辐射穿透 海底边界层 。,2007年10月,11,LASG 大洋环流模式,ML20 : 20 Levels, 4 5 (陈克明, 1994; 张学洪等,1996; 俞永强等, 1997) L30T63 : 30 Levels, 1.875 1.875 (金向泽等, 1999) LICOM 1.0 : LASG/IAP Climate System Ocean Model, 30 Levels, 0.5 0.5 (刘海龙, 2002；刘海龙等, 2004),2007年10月,12,2007年10月

5、,13,主要内容,1、次网格过程参数化 1.1 基本概念和理论 1.2 水平粘性方案 1.3 垂直混合方案 1.4 中尺度涡参数化 2、深对流 3、短波辐射穿透,1、次网格参数化 1.1基本概念和理论,2007年10月,15,海水中的分子运动,分子运动粘性系数,分子扩散,分子运动在距边界几毫米以内是重要的，对海洋内部的运动和示踪物扩散的影响可以忽略。,2007年10月,16,海洋中的湍流,流体微团 相对于分子足够大 相对于运动足够小 Reynolds数 海水是湍流流体 湍流运动在边界附近较强,By Leonardo da Vinci,2007年10月,17,海洋的中尺度运动,大气,海洋,100

6、0km,10100km Eddies,中尺度涡：10-100km 中尺度涡是海洋中最活跃的物理过程，在稳定状态下，海洋的混合主要是通过中尺度涡进行的。 中尺度涡的混合主要是沿着等密度面进行，穿越等面度面的分量较小。,2007年10月,18,海洋中的大尺度运动Ocean General Circulation,空间尺度100-1000 km，时间尺度月以上， 由风和浮力通量驱动，是气候研究和模拟的对象。,2007年10月,19,次网格尺度过程（1）,不能被模式的网格所分辨的过程，就称为“次网格”过程。次网格过程需要参数化。 “次网格”过程的空间尺度并是一个相对的概念。如涡分辨率模式。 提高模式的

7、分辨率可以减小模式对参数化的依赖程度。,2007年10月,20,次网格尺度过程（2）,海洋,10100km 中尺度涡,1cm 湍流,1mm 分子运动,1000km 大尺度环流,波长,能谱,次网格尺度,2007年10月,21,海洋模式的水平分辨率,Coarse: 2 Medium: 2/3 to 2 Eddy-permitting: 1/6 to 2/3 (涡相容的) Eddy-resolving: 1/6 (涡分辨的),目前用于气候研究的海洋模式大多为中等分辨率（1左右，约为100km），可以较好地分辨大尺度环流。,The classification of ocean horizontal

8、resolution（IPCC TAR, 2001）,2007年10月,22,Reynolds平均,Reynolds平均 对湍流足够大 对大尺度运动足够小,采样间隔2秒,Reynolds 平均的性质,2007年10月,23,Reynolds应力,已知,2007年10月,24,平均量方程组,X方向动量方程,Y方向动量方程,静力方程和连续方程,温度方程,5个方程9个未知数,2007年10月,25,扰动量方程,NS方程减平均量方程得,化简得,2007年10月,26,扰动协方差方程,并整理得,三阶项,2007年10月,27,湍流动能（TKE）方程,平均流的动能与TKE之间的转化,平均流的位能与TKE之

9、间的转化,TR表示输送和压力对TKE的再分配 表示TKE的耗散,2007年10月,28,方程组小结,扰动量u,全量方程,平均量方程,扰动量方程,扰动量协方差()方程,u,TKE方程,如何使平均量方程组闭合？,2007年10月,30,Prandtl的混合长假设,z,z,z-l,t0,t1,l即为混合长，可以理解为湍流涡旋所携带某种属性能够保持不变的距离，描述了湍流涡旋的空间尺度。,2007年10月,31,一阶闭合,Reynolds应力与平均速度梯度成正比， 引入涡动粘性系数,涡动粘性系数,2007年10月,32,一阶闭合的完整方程组,A为粘性/混合系数，各向异性。A可以随空间变化。 优点：形式简

10、单，易于实现、计算效率高， 尺度选择性，即有效地消除最小可分辨波数。,2007年10月,33,非局地混合,研究表明，水柱中特定层的湍流输运并不是仅仅依赖于局部的梯度，而是与整个水柱的状态有关，此即非局地（non-local）混合。,b为任意预报量，k为垂直混合系数 为非局地项，正比于表面通量， 反比于垂向摩擦速度和混合层深度,2007年10月,34,高阶闭合,采用扰动协方差方程、湍流动能方程， 利用观测或试验的假设，闭合方程。,2007年10月,35,小结,基本概念 分子运动、湍流、中尺度涡、大尺度环流、次网格尺度 湍流的基本理论 湍流的数学表示Reynolds应力 湍流方程的解法方程闭合问题

11、 一阶闭合（混合长理论，常用于水平混合的参数化） 高阶闭合（TKE方程，常用于垂直混合参数化） 非局地混合（应用于垂直混合参数化）,1.2 水平粘性方案,2007年10月,37,动量方程,2007年10月,38,调和形式方案（1）,Laplace方案,2007年10月,39,调和形式方案（2）,双调和方案（Biharmonic）,2007年10月,40,粘性系数,常系数方案 变系数方案 （1）在高纬度减小粘性系数 （2）与网格距的平方或3次方成比例,2007年10月,41,Smagorinsky方案（1963）,拉伸形变,切形变,分别为Karman常数和网格距,2007年10月,42,粘性系数

12、的约束,耗散方程的稳定性 平流和耗散的平衡 西边界流,网格距和时间步长,为特征平流速度,（1）稳定性条件是必须满足的。 （2）后两条约束如果不满足，虽然可以模式可以运行， 但是结果会变坏。,n为Munk边界层的格点数,2007年10月,43,例子,n=1,2007年10月,44,调和形式粘性方案的尺度选择性,为网格距，k为波数，A为扩散系数 当波长为2 时，即k= /，耗散最强,2007年10月,45,Laplacian和Biharmonic粘性方案比较,Laplacian,Biharmonic,时间,时间,波长,波长,0.9,0.9,0.1,0.1,2007年10月,46,小结,调和粘性方案的特点：尺度选择性、易用性。 粘性系数的约