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

Click to edit Master title style,Click to edit Master text styles,Second level,Third level,Fourth level,Fifth level,2007年10月,*,*,第五讲,海洋模式中的参数化过程,刘海龙,lhllasg,.,iap,.ac.,cn,2007,年,10,月,模式,动力框架和,参数化过程,动力框架（数值方法）,方程、网格、差分格式、积分方案等等,参数化（物理过程）,湍流过程、中尺度涡、辐射传输等等,2007年10月,2,参数化对模式的重要性,有关近似和参数化的相关概念是贯穿我们整个阐述过程的主题。正如我们所强调的那样，按照惯例用于解决海洋行为的运动方程通过复杂的一系列运动学近似、物理的参数化和数值假定得到。,任何一个或所有这些近似和参数化都可能对海洋模拟的质量产生重要影响。,因此对于从事海洋环流模拟研究的新手来说,，,知道解的敏感性和潜在问题的源头是十分重要的,。,（引自,Numerical ocean circulation modeling, by D. B.,Haidvogel,and A. Beckmann, 1999,）,2007年10月,3,什么是,参数化（,Parameterization,）,？,to express in terms of parameters,(Merriam-Websters Collegiate Dictionary),在,数值模式,中，不考虑过程的细节，而是用其它一些确定的变量所表示的,简化,函数表示这个过程，这个方法称为,“,参数化,”,。,(Glossary of PO and Related Disciplines),2007年10月,4,例子：,参数化次网格湍流,表示密度，,(),表示扰动量，,-,表示,Reynolds,平均,Reynolds,应力,单位时间内，单位面积，在,z,方向上输送的,x,方向的脉动动量通量的平均值。,K,表示湍流粘性系数,2007年10月,5,海,洋环流模式中哪些过程需要参数化？,纬圈和全球平均的海温分布,Thermocline,温跃层,Mixed layer,混合层,Tropical,Thermocline,Permanent,Thermocline,深对流,深对流,底边界层,短波辐射穿透,Temp,2007年10月,7,海水温度方程,深对流,短波辐射穿透,垂直混合,水平混合,中尺度涡,2007年10月,8,动量方程,垂直粘性,水平粘性,2007年10月,9,海洋模式中的主要参数化过程,动量方程中的水平,和垂直,粘性,温度和盐度方程中的,水平和,垂直混合,中尺度涡的混合和输运,深对流过程,短波辐射穿透,海底边界层,。,2007年10月,10,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月,11,http:/web.,lasg,.ac.,cn/FGCM/index,.,htm,2007年10月,12,主要内容,1,、次网格过程参数化,1.1,基本概念和理论,1.2,水平粘性方案,1.3,垂直混合方案,1.4,中尺度涡参数化,2,、深对流,3,、短波辐射穿透,2007年10月,13,1,、次网格参数化,1.1,基本概念和理论,海水中的分子运动,分子运动粘性系数,分子扩散,分子,运动,在距边界,几毫米,以内是重要的，对海洋内部的运动和示踪物扩散的影响可以忽略。,2007年10月,15,海洋中的湍流,流体微团,相对于分子足够大,相对于运动足够小,Reynolds,数,海水是湍流流体,湍流运动在边界附近较强,By Leonardo,da,Vinci,2007年10月,16,海洋的中尺度运动,大气,海洋,1000km,10100km,Eddies,中尺度,涡：,10-100km,中尺度涡是海洋中最活跃的物理过程，在稳定状态下，海洋的混合主要是通过中尺度涡进行的,。,中尺度涡的混合主要是沿着等密度面进行，穿越等面度面的分量较小。,2007年10月,17,海洋中的大尺度,运动,Ocean General Circulation,空间尺度,100-1000 km,，时间尺度月以上，,由风和浮力通量驱动,，是气候研究和模拟的对象。,2007年10月,18,次网格尺度过程（,1,）,不能被模式的网格所分辨的过程，就称为,“,次网格,”,过程,。次网格过程需要,参数化,。,“,次网格,”,过程的空间尺度并是一个,相对的,概念。如涡分辨率模式。,提高模式的分辨率可以减小模式对参数化的依赖程度。,2007年10月,19,次网格尺度过程（,2,）,海洋,10100km,中尺度涡,1cm,湍流,2,Medium: 2/3,to 2,Eddy-permitting: 1/6,to 2/3,(,涡相容的,),Eddy-resolving: 1/6,(,涡分辨的,),目前用于气候研究的海洋模式大多,为中等,分辨率（,1,左右,，,约为,100km,）,，,可以较好地分辨大尺度环流。,The classification of ocean horizontal resolution（IPCC TAR, 2001）,2007年10月,21,Reynolds,平均,Reynolds,平均,对湍流足够大,对大尺度运动足够小,采样间隔,2,秒,Reynolds,平均的性质,2007年10月,22,Reynolds,应力,已知,2007年10月,23,平均量方程组,X,方向动量方程,Y,方向动量方程,静力方程和连续方程,温度方程,5,个方程,9,个未知数,2007年10月,24,扰动量方程,NS,方程减平均量方程得,化简得,2007年10月,25,扰动协方差方程,并整理,得,三阶项,2007年10月,26,湍流动能（,TKE,）方程,平均流的动能与,TKE,之间的转化,平均流的位能与,TKE,之间的转化,TR,表示输送和压力对,TKE,的再分配,表示,TKE,的耗散,2007年10月,27,方程组小结,扰动量,u,全量方程,平均量方程,扰动量方程,扰动量协方差,(,),方程,u,TKE,方程,2007年10月,28,如何使平均量方程组闭合？,Prandtl,的混合长假设,z,z,z-l,t0,t1,l,即为混合长，可以理解为湍流涡旋所携带某种属性能够保持不变的距离，描述了湍流涡旋的空间尺度。,2007年10月,30,一阶闭合,Reynolds,应力与平均速度梯度成正比，,引入涡动粘性系数,涡动粘性系数,2007年10月,31,一阶闭合的完整方程组,A,为粘性,/,混合系数，各向异性。,A,可以随空间,变化,。,优点：形式简单，易于实现、,计算效率高，,尺度选择性，即有效地消除最小可分辨波数。,2007年10月,32,非局地混合,研究表明,，水柱中特定层的湍流输运并不是仅仅依赖于局部的梯度，而是与整个水柱的状态有关，此即非局地（,non-local,）混合。,b,为任意预报量，,k,为垂直混合系数,为非局地项，,正比于表面通量，,反比于垂向摩擦速度和混合层深度,2007年10月,33,高阶闭合,采用扰动协方差方程、湍流动能方程，,利用观测或试验的假设，闭合方程。,2007年10月,34,小结,基本概念,分子,运动,、,湍流、中尺度涡、大尺度环流、次网格尺度,湍流的基本理论,湍流的数学表示,Reynolds,应力,湍流方程的解法,方程闭合问题,一阶闭合,（混合长理论，常用于水平混合的参数化）,高阶闭合（,TKE,方程，,常用于垂直混合参数化）,非局地混合（应用,于,垂直混合参数化）,2007年10月,35,1.2 水平粘性方案,动量方程,垂直粘性,水平粘性,2007年10月,37,调和形式,方案,（,1,）,Laplace,方案,2007年10月,38,调和形式方案（,2,）,双调和方案（,Biharmonic,）,2007年10月,39,粘性系数,常系数方案,变系数方案,（,1,）在高纬度减小粘性系数,（,2,）与网格距的平方或,3,次方成比例,2007年10月,40,Smagorinsky,方案（,1963,）,拉伸形变,切形变,分别,为,Karman,常数和,网格距,2007年10月,41,粘性系数的约束,耗散方程的稳定性,平流和耗散的平衡,西边界流,网格距和时间步长,为特征平流,速度,（,1,）稳定性条件是必须满足的。,（,2,）后两条约束如果,不满足,，虽然可以模式可以,运行,，,但是结果会,变坏。,n,为,Munk,边界层的,格点数,2007年10月,42,例子,n=1,2007年10月,43,调和形式粘性方案的尺度选择性,为网格距，,k,为波数，,A,为扩散系数,当波长为,2,时，即,k=,/,，耗散最强,2007年10月,44,Laplacian,和,Biharmonic,粘性方案比较,Laplacian,Biharmonic,时间,时间,波长,波长,0.9,0.9,0.1,0.1,2007年10月,45,小结,调和粘性方案的特点：,尺度选择性、易用性。,粘性系,数的,约束,：,计算稳定性、耗散平流项导致的,扰动,、维持西边界流。,粘性系数和,分辨率：,在较低分辨率的海洋模式中,，使用调和方案可以得到较好的结果。在高分辨率模式中，可以较好体现部分中尺度涡谱。但是，当最小波长与,Rossby,变形半径接近时显得过于耗散,。,2007年10月,46,To be continuedThank you!,POP,，,0.1,，,Maltrud,and,McClean, 2005,2007年10月,48,