MOM4_L40模式对全球海洋CFC-11分布的模拟及其通风能力评估

CFC-11是评估全球海洋环流模式的一个重要工具,海水中溶解的CFC-11被用来分析全球海洋模式的通风模拟.本文在中国气象局国家气候中心发展的40层全球海洋环流模式(MOM4_L40)增加了示踪物CFC-11模块,然后利用该模式研究了CFC-11在全球海洋中的分布,并评估了模式的通风能力.对CFC-11的海表浓度、柱总含量以及大洋剖面的垂直浓度分布和渗透深度进行了分析,结果表明,与观测相比,模式较好地再现了CFC-11在海洋表面的水平分布特征,CFC-11主要储存区位于西北大西洋、副热带北太平洋及南大洋,其浓度分布与温度分布梯度相反.沿三个大洋的5个剖面的CFC-11垂直分布模拟也与观测基本吻合.模式模拟的CFC-11分布情况与全球平均经向流函数吻合,在南大洋模拟效果更加接近观测值,深海模拟效果较好,渗透深度接近观测.同时,模拟与观测相比也存在偏差.比如在北大西洋主要的存储区域,模式低估了CFC-11的吸收,这与高纬的CFC-11向低纬过度输送有关,可能是受温盐环流和强迫资料的影响.总体来说,MOM_L40模式模拟大洋吸收的CFC-11总量是理想的,通过模拟被动示踪物CFC-11很好地再现了海洋的通风能力.     

热带太平洋-印度洋温跃层海温异常联合模及其演变

利用SODA次表层海温再分析资料和卫星遥感海面高度异常数据,分析了热带太平洋和印度洋温跃层海温之间的联系,提出了太平洋-印度洋温跃层海温异常联合模(PITM)的概念、并定义了该联合模指数.结果表明,联合模指数具有准两年和3~5年的年际变化周期以及2011-2012年的年际变化周期,并具有季节锁相和振幅不对称等特征.联合模的演变过程与温跃层海温异常(TOTA)的发展和传播过程紧密相联:在太平洋,TOTA一般从西太平洋出发沿赤道(5°S-5°N)向东传播,到达东太平洋之后折向北,再沿10°N-14°N纬度带向西传播到达太平洋西岸并向赤道西太平洋扩展,形成一条回路;南太平洋也有类似回路但信号较弱;在印度洋,则主要沿8°S-12°S纬度带向西传播,到达西岸后折向北,然后迅速沿赤道(1.25°S-1.25°N)向东扩展,也形成一条回路.对NCEP/NCAR再分析风场资料的合成分析则表明,联合模的演变过程与大气环流尤其是纬向垂直环流(Walker环流)的变化密切相关,联合模的正位相对应着赤道印度洋区域顺时针的Walker环流以及赤道太平洋区域逆时针的Walker环流;而联合模的负相位则有相反的情况.此外,联合模演变过程中,TOTA的传播发展与850 hPa异常纬向风的传播发展有很好的相关.     

北太平洋海表面高度的年际变化及其机制

利用15年（1993～2007年）月平均的海表面高度（SSH）异常资料,分析了北太平洋海表面高度的年际变化的时空结构,并研究了热通量和风应力两个因子对其的强迫作用.结果表明,北太平洋年际时间尺度SSH变化的大值区在黑潮延伸区和西太平洋暖池区.EOF分解第一模态的空间结构沿纬向呈带状分布,第二模态为沿经向呈带状分布.热通量强迫作用在中纬度的东北太平洋可以解释SSH年际变化40%以上.风应力对SSH的作用包括正压和斜压两个方面.正压Sverdrup平衡模型模拟的SSH年际变化较弱,仅能解释高纬度副极地环流西部的20%~40%.由大尺度风应力强迫的第一阶斜压Rossby波模型可以解释热带地区的20%~60%,中纬度中部的20%~40%,以及阿拉斯加环流东部和副极地环流西部的20%~60%.风应力强迫的一阶斜压Rossby波模型对SSH的强迫机理又可分为局地风应力强迫和西传Rossby波作用.其中,风应力的局地强迫作用（Ekman抽吸）在东北太平洋、白令海以及热带中部有显著的预报技巧,可以解释SSH年际变异的40%以上.Rossby波的传播作用在中纬度海域的副热带环流中西部和夏威夷岛以东起着重要作用,可解释20%~60%.     

交互集合耦合模式系统对北太平洋SST变率和ENSO的模拟检验

清华大学地球系统科学研究中心在一个标准耦合模式（SC）的基础上建立了交互集合耦合模式系统（IE）,该系统可以实现多个不同大气模式或者同一大气模式采取不同初值组成的多个分量集合之后与海、陆、冰模式进行耦合.本文利用同一大气模式七个不同初值分量与其它模式分量开展在线集合耦合试验,利用积分稳定之后100年的试验结果,分析了IE在减小海-气界面大气噪音的情况下,对北太平洋海表面温度（SST）变率和ENSO的模拟,并与SC模拟结果进行了对比.分析表明,IE减小了北太平洋中高纬度SST方差的85%以上,表明该区域SST变率主要受大气的影响,且主要是通过改变海表湍流热通量实现的.黑潮延伸体区和北太平洋中部副热带涡旋区域平均SST 8年左右的低频周期主要受来自大气内部动力过程的驱动.在集合耦合模拟中,无论是副热带涡旋区SST与ENSO的联系,还是ENSO与北太平洋中高纬度SST的联系都能模拟出来,而标准模式未能模拟出这些现象,意味着大气噪音过强将掩盖ENSO与太平洋热带外SST的联系.IE对与ENSO关联的“太平洋-北美”（PNA）遥相关型的合理模拟,并通过湍流热通量对海表温度的影响,是其能够更好模拟ENSO与北太平洋中高纬度SST关系的重要原因.本文通过分析验证了所建立的交互集合耦合模式系统的合理性,揭示了该系统在海-气相互作用研究领域方面具有一定应用前景.     

华南前汛期降水异常与太平洋海表温度异常的关系

利用近50年华南地区站点逐日降水观测资料和全球大气、海洋分析资料,分析了华南前汛期降水异常的变化特征及其与太平洋海温异常的联系.结果表明,近50年来华南前汛期降水总体呈现减少趋势.影响华南前汛期降水异常的太平洋海温异常型是一个类似于ENSO的西太平洋暖池模态,即显著海温异常区域位于西太平洋暖池.西太平洋暖池区域(120°E-180°E,20°S-20°N)前期冬季海温异常同华南前汛期降水存在显著的负相关关系,是具有预报意义的海温关键区.该关键区海温异常影响华南前汛期降水的可能物理过程是:当前期冬季暖池异常偏暖时,菲律宾周围地区对流活动加强,导致Walker环流及东亚太平洋中低纬局地Hadley环流增强;该异常通过影响东亚－太平洋遥相关波列,使前汛期期间西太平洋副高加强西伸,脊线位置偏北,同时副热带西风急流减弱北退.随着Hadley环流上升支的增强,东亚副热带地区下沉运动也增强了,华南地区对流活动受到抑制.而且由于副高的增强,经过其北侧向华南地区的西南水汽输送辐合也减弱了,因此前汛期降水偏少.冷海温年的情形则相反,华南前汛期降水偏多.近50年来华南前汛期降水总体呈现趋势性减少正是由于前冬西太平洋暖池趋势性增暖所致.     

CMIP5模式对西北太平洋SST的模拟偏差及可能成因分析

利用参加第五次国际耦合模式比较计划(Coupled Model Inter-comparison Project Phase 5,CMIP5)的21个全球耦合气候模式的输出结果,利用海表热收支诊断方法,讨论了CMIP5模式对西北太平洋海表面温度(SST)的模拟偏差及其可能成因.结果表明,大部分CMIP5模式模拟的西北太平洋SST在全年都有持续的冷偏差,且在夏季最大,在冬季最小.CMIP5模式模拟的西北太平洋SST冷偏差是海盆尺度的,并且主要分布在中低纬度海洋内区.混合层热收支诊断分析研究表明,大气过程所致的海表净热通量模拟偏少是造成西北太平洋SST模拟冷偏差的主要原因,其中大气过程导致的潜热通量模拟偏多是造成冬、春、秋三季SST模拟冷偏差的最主要因素,而在夏季海表短波辐射模拟偏少和潜热模拟偏多对SST模拟冷偏差均有重要贡献.进一步分析表明,冬、春季西北太平洋低层风速偏大可能与海洋大陆区域的正降水偏差有关,而夏、秋季风速偏差则可能与热带辐合带的模拟偏差有关.冬、春季海洋大陆区域的正降水异常及夏、秋季ITCZ的模拟偏差使得北太平洋地区产生东北风异常,该东北风异常叠加在气候态风场上,使得北太平洋地区风速增强,造成海表潜热通量模拟偏多,海温模拟偏冷.     

北太平洋对人为二氧化碳吸收的数值模拟

采用一个开边界海盆尺度环流模式研究人为CO2在北太平洋的吸收和分布，并与闭边界模式的模拟结果进行了比较.模拟结果表明，西北太平洋和赤道东太平洋是两个重要的人为CO2汇.使用较大的等密度面扩散系数使得西北太平洋人为CO2通量增大，而赤道海区通量减小（RUN2）.与闭边界模式相比，开边界模式中该两个区域的人为CO2通量都增加了.1800～1997年间，北太平洋共吸收人为CO2 2375GtC(1Gt=1×1015 g, RUN1).副极地海区是人为CO2的一个重要输出区，能输出人为CO2吸收量的38％～54％，而20°N～30°N海区是人为CO2的一个重要贮存区，占整个北太平洋的24％.开边界对于10°N以南太平洋吸收和贮存人为CO2有很大影响.与基于观测资料的估计相比，虽然模式低估了西北太平洋的人为CO2的穿透，高估了东北太平洋的人为CO2的穿透，但总的说来，模式比较好地估计了人为CO2在北太平洋的贮存.     

边界层参数化方案及海气耦合对WRF模拟东亚夏季风的影响

区域气候模式的边界层参数化方案很大程度上影响着陆地-海洋-大气间水汽、动量及热量的交换,该方案的不确定性会给模式结果带来明显误差.本文基于WRF区域气候模式中四种常用的边界层参数化方案(YSU,ACM2,BouLac和MYJ)分别对东亚夏季风进行模拟研究,分析了不同的边界层方案对东亚夏季风环流及降水模拟的影响.结果表明,局地湍流动能方案BouLac和MYJ对东亚夏季风的模拟结果相对于非局地闭合方案YSU和ACM2更接近于观测,前者能更好的模拟出中国东部中低空西南风气流和西太平洋副热带高压.对于东亚夏季风降水,无论是空间分布还是季节内演变,BouLac和MYJ方案都要明显优于YSU和ACM2.此外,通过对比YSU和BouLac两种方案的模拟结果,发现边界层方案对东亚夏季风的模拟在海洋区域的影响更为显著.造成不同方案模拟差异的主要原因是非局地方案YSU和ACM2的边界层垂直混合偏强,使得海表向上输送的潜热通量明显偏强,对流更活跃,导致降水偏多以及相应季风环流的异常偏差.进一步研究指出缺少海气反馈过程使得WRF模式由边界层方案引起的模拟误差在海洋区域更为突出,引入海气耦合可以减小海表热通量误差并明显改善东亚夏季风的模拟结果.     

开边界海盆尺度环流模式模拟北太平洋CFCs分布

利用一个开边界北太平洋海盆尺度环流模式模拟了北太平洋CFCs的吸收和分布.模式计算了CFC_11海气界面通量及其在海洋中的输运和再分布，考察了开边界对模式模拟结果的影响以及开边界模式中不同沿等密度面扩散系数的影响.打开南边界可以很好地改善赤道以南海区的CFC_11模拟.增加沿等密度面扩散系数使得137°E和165°E断面模拟的高CFC_11水向南、向下输送的强度和影响范围都有明显增大，能很好地改善西北太平洋的模拟结果.西北太平洋是CFC_11的一个重要的汇区，此外在赤道东太平洋上升流区还存在一个小的汇区.     

变网格大气模式对1998年东亚夏季风异常的模拟研究

本文利用法国国家科研中心(CNRS)动力气象实验室(LMD)发展的可变网格的格点大气环流模式LMDZ4,对1998年东亚夏季降水进行了模拟,考查了变网格模式对东亚夏季降水的模拟性能.结果表明,模式在一定程度上能模拟出东亚夏季降水的极大值中心、夏季风雨带以及降水由东南向西北递减的空间分布特征.模式基本再现了1998年夏季两次雨带的进退特征,包括降水强度、雨带范围等,从而合理再现了1998年夏季江淮地区的"二度梅"现象.与观测相比,模拟的不足在于:在陡峭地形区附近存在虚假降水;江淮和华北地区以及四川盆地存在水汽输送的气旋式辐合偏差,同时高层环流辐散偏强,使得下层暖湿空气辐合上升、降水偏多;在东南地区存在反气旋式的水汽输送偏差,30°N以南地区降水偏少.对于1998年的"二度梅"现象,模拟偏差主要表现为长江中下游地区两次(特别是第二次)较强降水持续时间偏短,强降水范围偏小,而黄淮和华南地区却降水偏多.分析表明,模式对两次梅雨期降水的模拟偏差直接受环流形势模拟偏差的影响.LMDZ4区域模式版本的特点一是区域加密,二是加密区内预报场每10天向再分析资料恢复一次.敏感试验结果表明,LMDZ4加密区向强迫场的10天尺度恢复总体上有利于提高模式对华北降水的模拟能力,而对长江流域和华南降水的模拟具有不利影响.较之均匀网格模拟试验,加密试验由于在东亚的分辨率大大提高,对东亚夏季降水模拟效果更好.     

利用回归修正方法改善区域耦合模式中的ENSO模拟

通过对数值实验的比较和分析, 提出了一种旨在改善区域耦合模式中ENSO模拟的回归修正方法. 该方法主要用于修正耦合模式中海气间交换的通量. 具体步骤如下: 首先, 利用多年的观测资料计算得到驱动海洋模式所需的动量及热量通量, 驱动海洋模式进行长期积分; 其次, 用海洋模式模拟的SST作为大气模式的边界条件, 相应积分大气模式; 再利用大气模式模拟变量和相应观测资料建立线性关系, 通过线性拟合得到修正系数; 最后, 利用随时间和空间变化的回归修正系数修正计算动量及热量通量的变量, 并用修正后的变量计算海气交换通量, 进行耦合模式积分. 同时利用一个热带太平洋-全球大气耦合模式对该方案及常用的“距平耦合”方案进行了检验. 结果表明, 该方案优于“距平耦合”方案, 不仅可以更好的控制气候“漂移”现象, 而且, 能够改善区域耦合模式在热带太平洋区域的ENSO模拟.     

FGOALS－g快速耦合模式模拟的北太平洋年代际变率

本文分析了由中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室（LASG/IAP）最新发展的FGOALS_－g快速耦合模式300 a积分模拟结果,通过与多种观测资料的对比分析,讨论了北太平洋年代际变率的时空结构、主要年代际模态的演变特征以及与ENSO的联系等研究内容. 结果表明：该模式能成功模拟出北太平洋年代际变率的主要空间分布特征;模拟的年代际模态具有多时间尺度性,其中最显著的是周期约为10～20 a左右的准20年振荡模态,该模态上层海洋热容量异常的演变过程主要表现为大致沿副热带海洋涡旋做海盆尺度顺时针旋转的特征,相应的大气异常不仅与阿留申低压的变异有关,而且与太平洋-北美PNA）遥相关型以及上游的欧亚大气环流异常有密切关系;模拟的北太平洋年代际变率对年际ENSO循环的发生频率和强度有明显的调制作用. 但模拟的KOE区和阿拉斯加湾SST异常振幅比观测偏强,这与模式海冰偏多、高纬度SST偏冷的误差有关.     

海洋表面盐度与海洋环流和气候变化的关系

依据Argo的海洋盐度观测及降水、蒸发、海面高度等数据分析了全球海洋表层盐度的平均态和低频变化特征及其与海洋环流和气候变化的关系,并结合前人成果对海洋盐度的研究进行了回顾和展望.蒸发减降水在最低阶上决定了海洋表层盐度的分布,海洋环流动力过程在大部分海域调整了海洋表层盐度的空间结构,并影响了海洋表层盐度的低频变率.全球变暖的情况下,海洋表层盐度表现为盐度高的地方会变得更高、盐度低的地方会变得更低,表征为整个水循环系统的加强.在最近20年全球变暖停滞期,海表盐度变化与全球变暖的趋势不同,更多地体现了对气候系统的年代际及多年代际振荡的响应,此阶段对应着Walker环流加强;西太平洋和东南印度洋海表盐度的反向变化揭示了跨海盆尺度海洋动力过程的重要性,海洋Rossby波和印尼贯穿流的共同作用使得东南印度洋海表盐度显著降低;在东南太平洋和北大西洋更多地体现了全球变暖下的海表盐度单调升高或降低的响应.期待未来的盐度卫星能够提供更高时空分辨率的盐度数据,这将和Argo观测一起,拓展我们在中尺度动力学、区域海洋学以及气候变化等方面的认知.     

海温异常对西太平洋副热带高压脊面演变影响的机制研究

根据欧洲中期数值预报中心提供的大气环流资料和NOAA提供的海温资料, 应用改进的低阶谱模式方法研究了海温异常对西太平洋副热带高压脊面演变影响的物理机制. 结果表明, 在El Niño型海温强迫下, 大气环流内部动力过程中的大尺度波波和波流相互作用较弱. 随着外部热源强迫从冬季型向夏季型推进, 西太平洋副热带高压脊面北进不明显, 使得夏季西太平洋副热带高压脊面相对偏南. 在La Niña型海温强迫下, 大气环流内部动力过程中的大尺度波波和波流相互作用明显. 随着外源强迫从冬季型向夏季型转变, 西太平洋副热带高压脊面将随之向北移动19个纬度左右, 使得夏季西太平洋副热带高压脊面偏北; 且北移到达一定纬度后, 西太平洋副热带高压脊面表现出较为显著的30~60 d季节内南北振荡, 振幅在4°~7°之间.不同海温异常型所导致的大气内部动力过程差异是造成西太平洋副热带高压脊面年际变化的重要原因之一.     

FGOALS－g快速耦合模式模拟的北太平洋年代际变率

本文分析了由中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室（LASG/IAP）最新发展的FGOALS_－g快速耦合模式300 a积分模拟结果,通过与多种观测资料的对比分析,讨论了北太平洋年代际变率的时空结构、主要年代际模态的演变特征以及与ENSO的联系等研究内容. 结果表明：该模式能成功模拟出北太平洋年代际变率的主要空间分布特征;模拟的年代际模态具有多时间尺度性,其中最显著的是周期约为10～20 a左右的准20年振荡模态,该模态上层海洋热容量异常的演变过程主要表现为大致沿副热带海洋涡旋做海盆尺度顺时针旋转的特征,相应的大气异常不仅与阿留申低压的变异有关,而且与太平洋-北美PNA）遥相关型以及上游的欧亚大气环流异常有密切关系;模拟的北太平洋年代际变率对年际ENSO循环的发生频率和强度有明显的调制作用. 但模拟的KOE区和阿拉斯加湾SST异常振幅比观测偏强,这与模式海冰偏多、高纬度SST偏冷的误差有关.     

热带太平洋沃克环流变化的数值模拟

针对LASGBAP发展的耦合系统模式FGOALS-g2和FGOALS—s2,评估了其对热带太平洋沃克环流气候态的模拟能力,在此基础上,分析了沃克环流的变化特征,讨论了沃克环流变化的机理。对20世纪历史模拟结果的分析表明,两个模式均能够合理再现热带太平洋沃克环流气候态分布特征。观测中,过去百年（1900—2004年）和过去55年（1950～2004年）,沃克环流减弱,而近23年（1982～2004年）,沃克环流增强。在三个时间段内,FGOALS-g2模拟的沃克环流均减弱。FGOALS-s2中,过去百年赤道太平洋大气环流减弱,整个热带太平洋大气环流变化不明显;而在过去55年和近23年,模拟的沃克环流均增强。沃克环流变化模拟偏差与模式模拟的内部变率与观测不一致有关。降水与边界层向对流层输送的水汽相平衡的水循环约束关系可以很好地解释沃克环流的变化。FGOALS—g2中,在过去百年、55年和23年间,热带西太平洋降水相对变率（△P／P）增幅小于水汽相对变率（△q/q）增幅、东太平洋冷舌区△P/P增幅大于△q／q增幅,造成沃克环流减弱。FGOALS—s2中,在过去55年及23年间,热带西太平洋对流质量交换增强,东太平洋对流质量交换减弱,使得沃克环流增强。热带太平洋SST变化趋势分布型主导着沃克环流的变化。在过去55年和23年间,FGOALS—g2（FGOALS—s2）中,热带太平洋海表面温度（SST）变化表现为类El Nino（La Nina）型分布,对应沃克环流减弱（增强）。因此,气候系统模式合理模拟沃克环流变化的前提是对热带太平洋SST变化空间型的成功模拟。这一结论得到AMIP试验结果的支持。     

南海海面高度和输运流函数: 全球 变网格模式结果

为了研究南海环流结构和变异及其与外部水域的关系, 我们建立了一个嵌套于全球大洋环流模式中的高分辨率中国近海环流数值模式. 给出模拟所得南海月平均以及年平均的海面高度和流函数分布, 与TOPEX/POSEIDON资料比较表明, 所得海面高度距平与观测十分一致. 基于这些结果, 讨论了南海的环流结构, 尤其是上层环流结构. 结果表明: 对于表层海水来说, 黑潮在冬、春和秋季均通过吕宋海峡入侵南海, 夏季则表层没有入侵. 但对于整个海水而言, 全年均有海水从太平洋通过吕宋海峡进入南海. 这一差异表明, 在夏季, 太平洋的海水是在次表层和中层入侵南海的.南海北部陆坡附近全年受气旋式环流控制. 夏季的南海南部反气旋流圈、越南东南离岸流和冬季的南海南部气旋流圈都得到了很好的再现. 南海海面高度和海面高度距平之间的差异明显. 表明, 在利用卫星高度计资料研究南海的上层环流时, 长期平均海面高度的空间分布有重要意义.     

热带太平洋外海气耦合作用对东亚夏季风年代际变化的可能影响

利用1950—2000年逐月观测的热带太平洋海表温度分别驱动NCAR CCM3全球大气环流模式以及该模式耦合SOM(Slab Ocean Model)混合层海洋模式进行长时间积分试验,将两组试验结果相减来研究热带太平洋外海气耦合作用对20世纪70年代中后期东亚夏季风年代际变化的可能影响.模拟结果表明:热带太平洋外的海气耦合作用使得孟加拉湾、南海附近存在反气旋环流异常,导致我国东部存在偏南风异常,从而引起东亚夏季风年代际增强.其中热带印度洋的年代际增暖对东亚夏季风年代际增强有一定影响.     

厄尔尼诺现象的构造基础与激发因素

北太平洋对北极的半封闭状态和南太平洋对南极的开放状态是厄尔尼诺事件发生的构造基础，它导致北太平洋海表热能的积累和周期性向南太平洋输送，强潮汐振荡和火山喷发是其激发因素。     

伴随冬季北太平洋副热带海洋锋强度变化的大气扰动异常及对中纬度大气平均场的影响

中纬度海洋热力状况异常影响大气主要通过两种途径:非绝热加热的直接强迫作用和大气瞬变涡旋反馈的间接强迫作用,而后者的作用并没有被很好地认识.为了进一步理解间接强迫作用的物理机制,本文利用观测资料分析和区域大气模式模拟,研究了伴随冬季北太平洋副热带海洋锋强度变化的中纬度大气场异常,特别是对流层中高层不同频率的涡旋扰动活动的异常.实际观测和数值试验结果均表明,当北太平洋副热带海表面温度锋偏强时,其上空的中纬度大气经向温度梯度增强;对应此时的大气斜压性增强,且中纬度大气西风急流整层加速;然而增强的大气斜压性并不对应一致性增强的大气涡旋扰动活动.中纬度大气的涡旋扰动根据其生命周期,进一步划分为高频(2~7天)和低频(10~90天)涡旋扰动.研究结果表明偏强的北太平洋副热带海洋锋对应着增强的中纬度大气高频涡旋扰动和减弱的低频涡旋扰动;其中,中纬度大气高频扰动活动的增强,将有利于削弱中纬度大气经向温度梯度,从而减弱中纬度大气斜压性;而高频扰动对纬向风倾向项的正贡献,有利于中纬度急流中心北侧及下游区域的西风加速,形成中纬度西风相当正压结构的增强;大气低频扰动的减弱,对中纬度大气纬向风倾向项产生负贡献,不利于急流的纬向均匀化,而其热力强迫异常则有利于维持中纬度对流层中层大气的经向温度梯度.