基于二阶矩闭合模式平衡解的海浪生湍流混合系数解析估计

从次小尺度海洋运动方程和运动近各向同性假定出发,本文给出高确定性的海洋湍流二阶矩闭合方程和特征量变化方程.这组方程描述波动、涡旋和环流等较大尺度运动对海洋湍流的剪切不稳定性生成机制.同时,基于海浪破碎是湍流海面生成主要机制的认识,本文按破碎海浪统计理论给出湍流特征量变化方程的海面边界条件,包括由破碎能量损耗和卷入深度分别确定的海面动能通量和动能耗散率.根据上层海洋湍流特征量方程的平衡解和关于动能耗散率的观测资料,本文得到海浪生湍流混合系数的解析估计.所得到的海浪生湍流动能耗散率与观测资料进行了定性与定量对比,混合系数与先前Prandtl混合长度理论估计的一致性检验.     

地震海洋学方法在海洋混合参数提取中的研究与应用——以南海内波和地中海涡旋为例

混合是海洋中普遍存在的一种海水运动形式,对多个海洋学分支的研究具有重要的影响.随着物理海洋学的研究重心从大尺度向中小尺度现象过渡,近年来混合问题的研究重心也逐渐转向了中小尺度现象.内波与中尺度涡都是非常重要的中小尺度物理海洋学现象,对海洋能量在不同尺度中的级联发挥着重要的作用.本文基于地震海洋学研究了海洋混合参数的提取方法,并以南海内波和地中海涡旋为例进行了计算和分析.结果显示,南海内波在200~600 m深度范围内所引起的混合可达10-2.79 m2·s-1左右,比大洋的统计结果10-5 m2·s-1高出两个数量级以上.而地中海涡旋所引起的湍流混合率可达10-3.44m2·s-1左右,与大洋统计结果相比高出1.5个数量级左右,并且地中海涡旋下边界的混合要强于上边界,这一特征与前人的研究一致,另外涡旋上边界之上以及侧边界的外侧也具有非常高的混合率.     

南极冰盖极昼期间近地面湍流特征的实验观测

根据我国首次在东南极大陆冰盖上获得的超声风温仪直接测量资料,采用涡旋相关法计算和分析了极昼期间雪面上的湍流强度、湍流动能、湍流感热、粗糙度、拖曳系数和归一化方差.结果表明,南极冰盖表面的粗糙度为4.3×10~(-4)m,拖曳系数为1.8X10~(-3),湍流参数有明显的日变化.     

波浪对Ekman层的能量输入

基于海洋Ekman层流速剖面结构受波浪影响的3点观测事实,研究了波浪对Ekman层的机械能输入,在常涡粘系数的假定下得到了波浪对Ekman层能量输入的解析表达式．研究结果表明,海洋Ekman层的机械能输入包括两部分：一是风应力对Ekman层的直接能量输入,其受波浪的影响：二是波浪对Ekman层的能量输入,其与风应力、波浪参数和波向与风应力的相对方向有关．估计了波浪对Ekman层能量输入的大小,结果显示,在中高风速和中高纬度条件下,波浪对Ekman层的能量输入与不考虑波浪影响的Ekman层风能输入相比较可达10％以上,在绕极流区可达20％以上;在波浪方向与风应力方向不一致时波浪的能量输入更为显著．给出了波浪对非定常Ekman层能量输入的计算公式．研究结果在研究海浪的大尺度效应方面有启示意义．     

中层和热层大气的垂直动力传输及其影响

本文根据气体分子运动论导出了垂直动力传输函数,此函数可以代替Nicolet和Mangc的混合比,较好地描述高层大气中从完全混合状态向扩散分离状态的过渡过程。分析了高层大气中湍流运动的观测结果,给出了湍流系数随高度分布的几个模式,利用这些模式计算了惰性气体氩、氦与氮的数密度比的高度剖面,并与箭载质谱仪探空资料作了比较。根据光化-动力传输理论,计算了上中层及低热层大气主要化学成份的高度分布。结果表明:垂直动力传输明显地影响着大气成份的分布。文中还指出,氧分子的光离解系数与其数密度的高度剖面耦合在一起,因而应当由上向下数值求解有源有汇的连续方程组,这样可避免由于光离解系数的假设所造成的误差。     

强潮海湾中的湍流与混合:观测与数值模拟

对垂直剪切应力(雷诺应力)、湍动能剪切生成率与湍动能耗散率以及垂直涡黏性系数等湍流参数的现场观测与数值模拟,是理解强潮海湾中垂直混合过程及其动力与生物地球化学效应的基础.利用坐底ADCP流速观测资料、湍流剖面观测资料以及二阶湍流封闭模型,本文对一典型强潮海湾(厦门湾)中的湍流混合过程进行了观测与数值模拟研究.观测结果表明:在M2分潮占优的厦门湾中,湍动能剪切生成率与湍动能耗散率等湍流参数呈现明显的M4分潮周期,并具有显著的涨-落潮不对称性,涨潮时强,落潮时弱.在涨潮期间,整个水体在垂向上近似混合均匀,湍动能剪切生成率和垂直涡黏性系数的垂向平均值分别达到5×10-6 m2 s-3和2×10-2 m2 s-1;在落潮期间,强湍流则主要局限于厚度为5~8 m的底边界层内,且湍流强度随离底距离的增大而逐渐减小.基于重力位能异常控制方程的分析表明,潮流拉伸是造成湍流强度涨-落潮不对称性的主要原因,潮混合的贡献较小.基于二阶湍流封闭模型的模拟结果与实测结果的对比分析表明,在弱层化强潮海湾中,二阶湍流封闭模型能够较好地模拟湍动能耗散率与垂直涡黏性系等湍流参数的垂向结构与潮内变化;模拟结果与实测值之间的差异主要源于一些理想化的模型假设和数学近似.     

近地层能量闭合度对陆面过程模式影响

大量近地层观测试验表明,利用涡动相关法观测的湍流通量小于近地层可利用能量,即近地层能量是不闭合的,这种不闭合度一般为20%甚至更高.而陆面过程模式是基于地气间能量平衡建立,并且模式中的湍流边界层参数化方案通常根据实际观测的湍流通量来确定,因此能量不闭合必将对陆面过程模式造成一定的影响.本文利用2007年春季SACOL站的近地层观测资料,依据能量守恒将能量不闭合中的残余能量通过波文比分配到观测的湍流通量中,即修正涡动相关法观测的湍流通量使得近地层能量达到平衡;之后分别利用观测和修正的湍流通量,建立了能量不闭合和闭合情形下的湍流参数化方案,借助陆面过程模式SHAW,通过数值模拟和对比分析方法考察近地层能量闭合度对陆面过程模式的影响.研究结果表明近地层能量闭合对陆面过程模式有显著的影响:在陆面过程数值模拟中,当应用近地层能量不闭合的湍流通量形成的湍流参数化方案时,陆面过程模式会明显高估地表长波辐射及土壤温度;但当应用修正湍流通量使得近地层能量达到闭合形成的湍流参数化方案后,在不改变任何地表土壤物理生化属性的情况下,陆面过程模式能较好地模拟地表长波辐射和土壤温度.     

探知海洋七十载——2010年Crafoord奖获得者Walter H Munk教授成就解读

2010年Crafoord奖授予了著名物理海洋学家美国Scripps海洋研究所的Walter Heinrich Munk教授,以表彰他在海洋环流、海浪和潮汐及其在地球动力学中的作用等领域所做出的卓越贡献.本文通过介绍Munk教授七十多年的科研工作经历,解读他在物理海洋研究中所取得的部分卓越成就(如揭示了西边界流强化的物理机制,奠定了海浪和潮汐的预测技术,阐明了深海大洋混合的能量来源,提出了描述海洋内波特征的普适谱等),希望能够为中国科学家在科研求索之路上提供一些借鉴和启示.     

边界层参数化方案对边界层热力和动力结构特征影响的比较

本文利用高分辨率中尺度WRF模式,通过改变边界层参数化方案进行多组试验,评估该模式对美国北部森林地区边界层结构的模拟能力,同时比较了五种不同边界层参数化方案模拟得出的边界层热力和动力结构.结果表明：除个别方案外,配合不同边界层方案的WRF模式都能成功模拟出白天对流边界层强湍流混合特征和夜间稳定边界层内强逆温、逆湿和低空急流等热力和动力结构.非局地YSU、ACM2方案在白天表现出强的湍流混合和卷夹,相比于局地MYJ、UW方案,模拟的对流边界层温度更高、湿度更低、混合层高度更高、感热通量更大,更接近实际观测,这表明在不稳定层结下考虑非局地大涡输送更为合理,但局地方案在风速和风向的预报上存在一定优势.TEMF方案得到的白天局地湍流混合强度为所有方案中最弱,混合层难以发展,无法体现对流边界层内气象要素垂直分布均匀的特点.对于夜间稳定边界层的模拟,不同参数化方案之间的差异较小,但是YSU方案在一定程度上高估了机械湍流,导致局地湍流混合偏强,从而影响了其对稳定边界层的模拟能力.     

边界层参数化对海南岛山地环流结构和湍流特征模拟的影响

本文基于WRF-ARW (V4.0)中尺度数值模式, 选用耦合同一近地层参数化方案(Eta)的五种边界层参数化方案(MYJ、MYNN2、MYNN3、BouLac、UW), 对2020年5月1—2日海南岛一次典型山地环流个例进行模拟, 对比分析了这五种方案所模拟的山地环流结构和湍流特征的差异, 旨在为模式应用于山地环流研究和模式改进提供一定的科学依据.研究结果表明, 这五种湍流动能边界层参数化方案均能模拟出山谷风特征, 对环流结构和湍流特征的影响表现为谷风时段大于山风时段.对于山地环流水平结构的模拟, 因平原风推进距离的不同, 五种方案模拟的近地面风速差异可达4 m·s^-1以上, 其中MYJ方案模拟的谷风最强盛, 而MYNN3方案最弱, 山区多为静风或小风.对于山地环流垂直结构的模拟, MYNN2、UW方案模拟的谷风环流较强, 表现为谷风厚度较厚, 推进距离较远, 同时由于模拟的谷风环流可越过山顶, 从而模拟的高海拔地区上升区的覆盖范围和强度均较大; MYJ、BouLac方案模拟的谷风环流均未能越过山顶, 且其中BouLac方案的平原风环流未能与上坡风环流耦合; 而MYNN3方案模拟的环流结构最不明显.湍流特征对比分析表明, 五种方案模拟的湍流动能强度和分布特征存在显著差异, 其中MYNN2方案湍流混合最强, MYJ方案最弱.由于MYNN2方案下垂直湍流混合较强, 大气不稳定度增加, 因此其对应的谷风厚度较大, 同时该方案模拟的边界层结构与实况最为接近, 更适用于海南岛山地环流的模拟.

     

青藏高原上空UTLS区域一次地形重力波过程中的物质上传

利用美国航空航天局MERRA(Modern-Era Retrospective Analysis for Research and Applications)再分析资料和MODIS(Moderate-Resolution Imaging Spectroradiometer)卫星资料以及欧洲气象中心ECMWF-Interim(European Centre for Medium-Range Weather Forecasts)再分析资料,分析了发生于青藏高原北侧上空的一次地形重力波事件,并使用中尺度预报模式WRF-ARW.V3.0(Weather Research and Forecasting model,V3.0)对其进行了数值模拟.在此基础上,诊断分析了此次地形重力波在UTLS(Upper Troposphere and Lower Stratosphere)区域造成的物质和能量垂直传输特征.分析结果表明这一中尺度地形重力波信号的水平波长约为600km,与地形扰动水平尺度接近,重力波在对流层中传播的垂直波长约为3km,在垂直方向上随着高度的增加呈现出由东向西倾斜的结构特征.此次地形重力波上传进入平流层并在150hPa附近破碎,波破碎后动量通量在短时间内发生了强烈的衰减,重力波携带的能量在破碎高度附近释放.重力波破碎的同时垂直方向湍流混合变得异常强烈,湍流交换系数可在短时间内增加到背景值的8倍以上,剧烈湍流混合过程导致了对流层上层的空气进入平流层,使下平流层空气出现了位势涡度和臭氧的低值区,在浮力频率的垂直剖面中也可以看到由于地形重力波过程造成的平流层下层浮力频率异常低值区.     

2001年发生厄尔尼诺事件的天文条件

用潮汐产生的地球流体与固体的差异旋转,可以解释1500-1800a周期的气候变化.强潮汐加大垂直方向和水平方向海水的混合,将太平洋深层冷水翻上表面(或使太平洋西部暖水流向东部,北部暖水流向南部),使海洋上方空气变冷(或变暖),产生拉尼娜冷事件(或厄尔尼诺热事件).天文资料表明,强潮汐与厄尔尼诺事件有很好的对应关系,火山活动亦受强潮汐的控制.火山喷发使海洋底层暖水上升到海洋表面,火山灰遮蔽阳光使气候变冷,它们是控制厄尔尼诺事件发生的重要外强迫因子.当火山灰在低纬度地区上升到平流层,较小的气溶胶可在数月内传播至全球,并可在平流层内持续飘浮1-3a,最后降落在两极地区,减弱被遮蔽区的太阳辐射,与该区发生日食的效果相同,是厄尔尼诺事件的延迟因子.     

不同垂向分布特征的垂直黏性系数对经典定常Ekman螺旋结构的影响

经典定常Ekman螺旋结构的表面漂流应该在南北半球分别沿风向右偏和左偏45°,而已有观测和研究结果表明:在北极和高纬度海域,表面流相对海面风向的偏角小于45°;在大洋低纬度海域,该偏角大于45°.针对该现象,本文设计了理想化数值实验,研究了不同垂向分布特征的垂直黏性系数对经典定常Ekman螺旋结构的影响,结果显示:垂直黏性系数的垂向分布特征对Ekman螺旋结构影响显著,当垂直黏性系数随深度增加而减小时,表面流偏角大于45°;当垂直黏性系数随深度增加而增大时,表面流偏角小于45°;观测到的低纬度海域与极区表面流偏角的不同应该主要由海洋上层垂直黏性系数的不同垂向分布特征产生;螺旋扁平度不等于常数1,但无明显的变化规律;体积输运的大小和方向与经典Ekman理论结果一致,不受垂直黏性系数分布的影响.     

森林冠层上下湍流谱结构和耗散率

利用新型三维超声风速 /温度仪 ,根据在长白山原始森林冠层上下两个高度上测量的湍流资料 ,采用涡动相关法计算和分析了森林冠层上下湍流谱的结构、局地各向同性和耗散率 .结果表明 ,在惯性副区 ,归一化湍流谱遵从 - 23的指数规律 ;在森林冠层上 ,尽管谱的形状与均匀表面的一致 ,但是 ,大气湍流是非各向同性的 ;而在森林冠层内则是近似各向同性的 ;森林冠层上下湍流能量和热量耗散率比典型草原和牧场的结果大 ,揭示了森林粗糙表面在湍流输送过程中的动力扰动和对大涡的破碎作用     

风垂直切变对大气静力适应过程的影响——第三部分:能量转换

基于前两部分的研究,第三部分从能量转换的角度揭示了风垂直切变对大气静力适应过程的影响.研究显示,对非静态等温大气,初始时刻加以有效势能和有效弹性势能的强迫后,风垂直切变改变了适应过程中扰动能量在各能量形式中的分配比例,这种作用随扰动尺度和时间的变化有显著差异,系统尺度越小影响越显著.比较东风急流模型和西风急流模型显示,负的风切变应力使得四类扰动能量谱随系统尺度的变化趋于平缓,使得微尺度和小尺度(中尺度)系统中扰动垂直动能(扰动水平动能)的比例随时间减小,正的风切变使得扰动能量谱随系统尺度的变化显现跳跃的特征,使得中尺度(微尺度和小尺度)系统中扰动垂直动能(扰动水平动能)的比例随时间减小.风垂直切变引起的动量通量输送是扰动能量与基本气流能量之间交换的通道,当风切变应力和动量输送方向相同时,风垂直切变使得能量向基本态转移,维持基本气流,削弱扰动能量,缩短静力适应过程的周期;当风切变应力和动量输送方向相反时,风垂直切变作用相反;此作用随着扰动尺度的变化差异明显.     

斜压对流边界层中风廓线和湍流特征的研究

在美国Ｋａｎｓａｓ和Ｎｅｂｒａｓｋａ交界处于１９９２年进行的ＳＴＯＲＭ－ＦＥＳＴ（ＳＴｏｒｍ－ｓｃａｌｅ Ｏｐｅｒａｔｉｏｎａｌ　ａｎｄ　Ｒｅｓｅａｒｃｈ　Ｍｅｔｅｏｒｏｌｏｇｙ－Ｆｒｏｎｔｓ　Ｅｘｐｅｒｉｍｅｎｔ Ｓｙｓｔｅｍ　Ｔｅｓｔ）资料基础上,分析了 斜压对流边界层中温度、湿度和风速的垂直分布．分析表明,在对流边界层中温度和湿度常 显示很好的混合;但风速分量有时混合很好,有时则存在切变．计算了存在和不存在风切变 两种情况下的湍流动能收支．最后,讨论了对流边界层中风切变形成的可能原因．     

地球自转与El Ni 0--海气耦合理论

从低纬的海气耦合的浅水模式方程组出发,运用正交模和特殊函数的方法进一步讨论地球自转速率变化对海气耦合系统的影响.研究表明:地球自转速率的变化通过海气耦合一方面使大气和海洋的Kelvin波和Rossby波的移动及稳定性发生变化,另一方面使纬向风、洋流和海表温度发生变化.特别是在地球自转减慢时,通过海气耦合,出现纬向风和洋流异常和大洋东部海表温度增加,从而导致引起全球气候异常的El Ni o现象.     

地球自转与El Nio——海气耦合理论

从低纬的海气耦合的浅水模式方程组出发 ,运用正交模和特殊函数的方法进一步讨论地球自转速率变化对海气耦合系统的影响 .研究表明 :地球自转速率的变化通过海气耦合一方面使大气和海洋的Kelvin波和Rossby波的移动及稳定性发生变化 ,另一方面使纬向风、洋流和海表温度发生变化 .特别是在地球自转减慢时 ,通过海气耦合 ,出现纬向风和洋流异常和大洋东部海表温度增加 ,从而导致引起全球气候异常的ElNi no现象     

基于全球气象再分析资料的InSAR对流层延迟改正研究

气压、温度和水汽含量等大气物理参数的时空变化导致的对流层延迟是制约合成孔径雷达干涉测量(Interferometric Synthetic Aperture Radar,InSAR)高精度应用的重要因素之一.最新研究显示气象再分析资料在补偿对流层延迟影响方面具有巨大的应用潜力,这促使我们对其有效性和鲁棒性做进一步的研究和探索.本文首先推导了利用气象再分析资料对InSAR进行对流层延迟校正的算法;然后以美国南加州地区的ENVISAT ASAR数据为例,分析了基于两种气象再分析资料(ERA-Interim和North American Regional Reanalysis,NARR)校正InSAR对流层延迟改正的效果;通过与MERIS水汽延迟改正结果比较,验证了该方法的有效性.实验结果表明:(1)不能简单忽略干延迟,可通过气象再分析资料进行有效估计;(2)通过与MERIS水汽产品获得的对流层延迟比较发现,气象再分析资料能够取得接近于MERIS的改善效果;(3)对ERA-Interim和NARR两种气象再分析资料而言,虽然后者具有更高的时间和空间分辨率,但在改正InSAR对流层延迟方面并没有表现出比前者更明显的优势;(4)气象再分析资料可以很好地估计与地形强相关的垂直分层延迟,但对于小尺度的湍流混合延迟的捕捉能力有限.综合分析认为,气象再分析资料的优势在于其数据可随时获得、免费和全球覆盖,它可以显著减弱大尺度的垂直分层延迟对干涉图相位的影响,从而有助于InSAR获取更真实可靠的地形高程和地表形变信息.     

北京城市大气边界层低层垂直动力和热力特征及其与污染物浓度关系的研究

利用2001年3月19~29日和2003年8月11~25日中国科学院大气物理研究所325m大气边界层观测塔资料,分析研讨了北京城市大气边界层低层的垂直动力结构特征及其与污染物浓度分布的关系.对比分析了北京城市大气边界层低层不同高度的风、温度和湿度梯度资料、大气湍流和大气化学观测系统资料,综合分析获取了无因次速度、温度湍流方差和湍流通量、湍能分布特征及其与污染物浓度空间分布的关系,同时分析了北京地区沙尘天气过程中城市边界层低层垂直结构特征及其污染物浓度的分布与变化特征.分析结果表明,在不稳定层结条件下,47和120m高度上无因次速度湍流方差和温度湍流方差遵循莫宁-奥布克霍夫相似规律,并给出相应的拟合公式.稳定大气边界层可按层结参数z′/L分成二分区,z′/L<0.1为弱稳定区,此时相似规律可适用,z′/L>0.1为强稳定区,在此区内无因次速度方差随稳定度增大有增大的趋势,而无因次温度方差则保持不变.白天近地层包含了47和120m,而280m已在近地层之上.对2001年3月北京地区一次沙尘天气过程的城市边界层资料分析发现,320m高度上总悬浮颗粒物浓度最高达到913.3μg/m3,在边界层内大气颗粒物从上层向低层输送,这与锋面过境时低空急流从上层向下发展过程并伴随的强下沉运动有关.