冬季北太平洋海表热通量异常和海气相互作用——基于一个全球海气耦合模式长期积分的诊断分析

利用一个全球海气耦合模式长期积分所给出的资料,分析了冬季北太平洋海表湍流热通量（潜热和感热）异常及其对海表温度（SST）异常的影响,并比较了海表热通量诸分量和海洋内部的动力学过程对SST变化的相对重要性。结果表明,冬季热带外海洋上的湍流热通量是影响SST的主要因子,但在北太平洋中部海水的平流作用也不可忽视。冬季热带外海洋向大气释放的潜热和感热通量与SST倾向（而不是SST本身）之间存在着显著的相关,这同Cayan和Reynolds等利用COADS资料和NCEP资料同化模式分析的结果是一致的。模式诊断的结果支持这样一种看法：和热带海洋不同,冬季热带外海洋上的海气相互作用主要地表现为大气对海洋的强迫作用,而不是相反。模式给出的SST倾向的第一个EOF分量及其与海平面气压场的相关特征同Wallace等从观测资料分析所得到的结果是一致的;进一步的分析表明：在冬季北太平洋的大部分区域（特别是西太平洋）,大尺度大气环流异常在很大程度上决定着SST的异常,而这种决定作用正是通过它对湍流热通量的强烈影响来实现的。     

西太平洋-南海地区潜热通量长期变化趋势的南北差异及成因分析

利用1958~2014年美国伍兹霍尔海洋研究所客观分析海气通量项目(OAFlux)的月平均潜热通量和相关气象要素数据,以及NCEP/NCAR再分析表面气压数据,通过Trend-EOF分析方法,本文研究了西太平洋—南海地区潜热通量的长期变化趋势。发现西太平洋—南海地区潜热通量整体呈上升的趋势,其中冬季上升趋势最强。冬季潜热通量趋势存在明显的南北差异,特别是在南海地区,南海北部为上升趋势而南部为下降趋势。南海北部以及菲律宾海地区冬季潜热通量上升的主要原因是海气比湿差的增大,而南海南部潜热通量呈下降趋势,在东侧主要原因是风速减小,在西侧主要原因是海气比湿差减小。南海潜热通量长期趋势的南北差异是风速和海气比湿差的共同作用造成的。另外,研究发现风速变化趋势受到局地环流变化的影响,在表面气压下降中心线以北地区为上升趋势,在其以南为下降趋势,而海气比湿差的变化趋势则主要取决于海表温度的变化趋势。     

印度洋海气热通量交换研究

基于综合海洋大气资料集(COADS)资料的研究表明,热带印度洋的海气热通量交换具有明显的区域性特征,在部分海域,如冬季热带印度洋的中东部、夏季的热带西印度洋和北印度洋,它主要表现为海洋对大气的强迫.海洋对大气的这种强迫,主要是通过潜热加热实现的.与潜热加热相比,感热加热尽管是一个小量,但感热异常与表层海温的显著相关,较之潜热明显超前.无论冬季还是夏季,热带印度洋都存在大面积海域,其SST变化难以通过海气热通量交换来解释.     

南海地区热通量的时空变化特征

本文利用美国 NCEP1958-1998 年高斯网格月平均再分析资料,分析了南海及周边地区(0～20°N,100～125°E)热通量(包括潜热通量和感热通量)的时空变化,结果表明该潜热通量、感热通量具有明显的季节转换特征.南海中部海区是潜热通量、感热通量季节变化最剧烈的关键区,南海季风对潜热、感热输送均有影响,并且蒸发潜热输送大于感热输送.EOF 分析表明,风速对潜热、感热输送贡献较大,另外气温和相对湿度对潜热输送有贡献,而水温与气温差对感热输送有贡献.整个南海地区潜热通量、感热通量具有明显的年际变化特征,潜热通量存在5～8 a 以及 2 a 左右的周期振动,而感热通量只有 5～8 a 的周期振动.     

阵风效应对冬季北太平洋涛动年际变化模拟效果的影响

利用观测海温驱动NCAR CAM5模式进行了两组数值试验,对比分析了在CAM5模式海气湍流热通量参数化方案中引入阵风效应前后,模式对冬季NPO(North Pacific Oscillation,北太平洋涛动)年际变化模拟效果的差异。结果表明,考虑了阵风效应后,CAM5模式能较好地模拟冬季NPO的年际变化特征,模拟的1979—1999年冬季NPO指数序列与观测结果的相关系数由0.09提高到0.57;热带东太平洋海表温度异常通过影响海表湍流热通量异常,对冬季NPO年际变化产生影响;引入阵风效应后,模式对热带东太平洋海表湍流热通量异常的模拟结果更趋合理,从而使得模式可以较好地模拟出冬季NPO的年际变化。     

北太平洋冬季年代际尺度海气相互作用

为了更深了解北太平洋年代际尺度海气相互作用,本文通过对湍流热通量和SST的相关关系进行了研究。当热通量和SSTA正相关时,代表海洋影响大气,而当二者相关系数为负时,表示大气强迫海洋。将整个北太平洋SSTA求平均和把北太平洋分成3个区域求SSTA平均,分别对热通量和SST的关系进行研究。结果表明,北太平洋冬季年代际尺度上的海气相互作用主要表现为海洋强迫大气,而在日本附近海域以及北太平洋中东部,海洋对大气的强迫作用最明显,因此这两个区域是北太平洋冬季年代际尺度上海气相互作用的关键区。     

太平洋湍流热通量异常与中国夏季降水年代际变化

采用NCEP/NCAR再分析感热、潜热资料,利用EOF方法,分析了夏季太平洋湍流热通量年代际变化的时空特征及其与中国夏季降水的关系。分析结果显示,1951—2003年夏季太平洋湍流热通量距平的年代际分量的EOF前3个模态的方差贡献率分别为32.27%、12.52%、11.10%,且这3个模态分别与中国华北、华南和东北地区夏季降水的年代际关系密切;太平洋湍流热通量年代际异常EOF第1模态的时间系数与东亚夏季风指数年代际变化有显著的负相关。20世纪70年代中期以前,太平洋湍流热通量异常偏多,东亚夏季风偏强,对应华北降水偏多;20世纪70年代中期以后,太平洋湍流热通量异常偏少,东亚夏季风偏弱,对应华北夏季降水减少,干旱化趋势明显。     

东中国海夏季潜热通量的时空特征及其与中国东部降水的联系

利用1985—2006年OAflux、NCEP/NCAR再分析资料和中国大陆东部的120个测站降水资料,应用EOF分解和线性回归等方法分析了东中国海夏季海气热通量时空特征及其与中国东部夏季降水的关系。结果表明,东中国海夏季潜热通量EOF1模态空间场的强信号区为黄海、东海至台湾海峡和南海北部局部海域,时间系数(PC1)表现出显著的年际变化和1997年前后的突变。降水和经向环流的PC1回归分析还表明,夏季潜热通量的异常变化与中国东部降水和经向环流的变化密切相关,潜热通量正异常对应黄淮地区和华南地区的上升气流正异常,降水偏多,而东北与华北地区为下沉气流正异常,降水偏少。1997年前后东中国海潜热通量发生突变之后呈明显的上升增强趋势,使得上述现象更为显著,即:对应东中国海夏季潜热通量的正异常,黄淮地区和华南地区为上升气流正异常,降水偏多;东北地区和华北地区为下沉气流正异常,降水偏少。这表明夏季东中国海及邻近海域潜热通量的异常变化是导致中国东部汛期降水年际和年代际异常的重要原因之一。     

加利福尼亚附近的海温强迫及其与北太平洋年代际振荡的可能联系

本文利用1958～2018年期间海表面温度（SST）异常和湍流热通量异常变化的关系,探讨了其与北太平洋年代际振荡（PDO）相关的年际和年代际时间尺度上在不同海域的海气相互作用特征。结果表明:在年际尺度上,黑潮—亲潮延伸区（KOE）表现为显著大气强迫海洋,赤道中东太平洋表现为显著海洋强迫大气;在年代际尺度上,PDO北中心表现为大气强迫海洋,加利福尼亚附近则表现为显著海洋强迫大气。进一步分析表明:加利福尼亚附近区域是北太平洋准12年振荡的关键区域之一,与PDO准十年的周期类似,加利福尼亚附近的冷（暖）海温对应其上有反气旋（气旋）型环流,赤道中太平洋海水上翻和北太平洋东部副热带区域经向风应力的变化是北太平洋准12年振荡的另外两个重要环节。     

春夏季青藏高原与伊朗高原地表热通量的时空分布特征及相互联系

伊朗高原和青藏高原热力作用对东亚区域气候具有重要影响。基于1979—2014年欧洲中心ERA-interim月平均再分析地表热通量资料,分析了春、夏季青藏高原与伊朗高原地表热通量的时、空分布特征以及春、夏季青藏高原与伊朗高原地表热通量的关系。结果表明,春、夏季青藏高原与伊朗高原地表热通量在季节、年际和年代际尺度上具有不同的时、空分布特征。对于青藏高原,春、夏季地表感热呈西部大东部小、地表潜热呈东部大西部小;地表感热在春季最大且大于地表潜热,地表潜热在夏季最大且大于地表感热。在年际时间尺度上,春、夏季青藏高原地表热通量异常的年际变化在东、西部不一致,青藏高原西部,地表感热与地表潜热有较强的负相关关系。青藏高原地表感热异常具有很强的持续性,当春季地表感热较强（弱）时,夏季高原地表感热同样较强（弱）。青藏高原东部与西部地表热通量的年代际变化有明显差异,春（夏）季青藏高原东部地表感热呈显著的年代际减弱趋势,1998（2001）年发生年代际转折,由正异常转为负异常;而青藏高原西部地表感热在春季则有显著的增大趋势,2003年发生年代际转折,由负异常转为正异常。青藏高原东部地表潜热仅在春季为显著减弱趋势,2003年出现年代际转折,由正异常转为负异常;青藏高原西部地表潜热在春、夏季都有显著减弱趋势,年代际转折出现在21世纪初,由正异常转为负异常。对于伊朗高原,春、夏季地表热通量的空间分布在整个区域较一致,地表感热在夏季最大,地表潜热在春季大、夏季小,但各季节地表感热都大于地表潜热。相对于青藏高原地表感热,伊朗高原地表感热在各月都更大。在年际时间尺度上,春、夏季伊朗高原各区域地表热通量异常的年际变化较一致;地表感热与潜热有很强的负相关关系;伊朗高原地表感热、潜热异常都具有持续性,当春季地表感热（潜热）通量较强（弱）时,夏季地表感热（潜热）通量同样较强（弱）。伊朗高原北部与南部地表热通量的年代际变化存在差异。其中,春、夏季伊朗高原北部地表感热（潜热）呈显著增强（减弱）趋势,在20世纪末发生了年代际转折,春、夏季北部地表感热（潜热）由负（正）异常转为正（负）异常。而伊朗高原南部春、夏季地表热通量无显著变化趋势,但春季地表感热、潜热与夏季地表感热同样在20世纪末存在年代际转折,地表感热（潜热）由负（正）异常转为正（负）异常。春、夏季两个高原地区地表热通量的关系主要表现为:就春季同期变化而言,伊朗高原地表感热与青藏高原西部地表感热具有同相变化关系,与青藏高原东部地表感热具有反相变化关系,伊朗高原地表潜热与青藏高原东部地表潜热具有同相变化关系;就非同期变化而言,春季伊朗高原地表感热与夏季青藏高原东部地表感热存在反相变化关系。      

Turbulent fluxes over inhomogeneous landscape

Mesoscale surface turbulent fluxes over a complex terrain surrounded by oceans have been investigated using a 3-D numerical mesoscale model, under conditions with and without synoptic flows. The study indicated that under synoptically calm condition, the allocation and intensity of mesoscale surface turbulent fluxes (MSTFs) were greatly impacted by the thermally forced mesoscale circulation (TFMC) over mesoscale heterogeneous landscape. The max-imum values of sensible (Hs) and latent (LE) heat fluxes were located over the convergent zones and considerably im-pacted by the soil wetness (M), but did not depend strongly on the atmospheric background thermal stability (β0). The simulated results suggested that the sensible heat flux was closely proportional to the square of wind speed in the surface layer. By the action of synoptic flow, the allocation of LE was shifted to downwind, its intensity increased.     

Nonstationarity of turbulent heat fluxes at Summit,Greenland

Turbulence data collected over a total of 25 days during two summers are used to describe processes responsible for the nonstationarity of turbulent sensible heat fluxes at Summit, Greenland. A stationarity test shows that about 40% of the data are classified as nonstationary. Three main factors are explored to account for the large fraction of nonstationary runs: (1) intermittency of turbulence in stable conditions, (2) changes in net all-wave radiation in response to cloud forcing, and (3) diurnal trends in stability. A classification procedure that accounts for the intermittent nature of turbulence shows that during stable, nonstationary conditions 50% of the total sensible heat flux is realized in 22% of the sampling time. Intermittency often occurs at Summit during periods characterized by weak and irregular horizontal winds in combination with strong stability. Rapid changes in net all-wave radiation in response to cloud forcing results in nonstationarity during unstable conditions. Between 0930–1130 and 1900–1930 UTC turbulent heat fluxes are not only small in magnitude but also typically change sign, with nonstationarity during these periods often as high as 65%. These results should help resolve some of the present uncertainties in obtaining reliable fluxes at this site, in particular under stable atmospheric conditions.     

Estimations of Land Surface Characteristic Parameters and TurbulentHeat Fluxes over the Tibetan Plateau Based on FY-4A/AGRI Data

正Nan GE1, Lei ZHONG*1,2,3, Yaoming MA4,5,6, Yunfei FU1, Mijun ZOU1,Meilin CHENG1, Xian WANG1, and Ziyu HUANG1     

Evaluation of the ERA-40 air–sea surface heat flux spin-up

Ocean models depend strongly on surface fluxes. When computed from atmospheric models, fluxes are affected by spin-up, i.e. they increase (or decrease) with forecast length. Such behavior may bias ocean models. The European Centre for Medium Range Forecasts (ECMWF) 40-year re-analysis (ERA-40) has been used to quantify short-range spin-ups of radiative and turbulent heat fluxes. Fluxes are compared as differences between two runs with the same initialization time. This method allows flux analysis over short-range forecasts as a function of distance from the initialization time. Results indicate that (i) latent heat flux spin-up increases with time but levels off after 24 h; (ii) sensible and radiative flux spin-ups remain constant after 6 h; (iii) regional spin-up of turbulent fluxes are systematic and can be larger than 30% for sensible heat but never exceeds 15% for latent heat; (iv) spin-up depends upon the season. The same analysis has been carried out with the ECMWF 15-year re-analysis (ERA-15); spin-ups in ERA-40 have been generally smaller than those in ERA-15.     

Measured And Simulated Latent And Sensible Heat Fluxes At Two Marine Sites In The Baltic Sea

In this study, turbulent heat flux data from two sites within the Baltic Sea are compared with estimates from two models. The main focus is on the latent heat flux. The measuring sites are located on small islands close to the islands of Bornholm and Gotland. Both sites have a wide wind direction sector with undisturbed over-water fetch. Mean parameters and direct fluxes were measured on masts during May to December 1998.The two models used in this study are the regional-scale atmospheric model HIRLAM and the ocean model PROBE-Baltic. It is shown that both models overestimate the sensible and latent heat fluxes. The overestimation can, to a large extent, be explained by errors in the air-water temperature and humidity differences. From comparing observed and modelled data, the estimated 8-month mean errors in temperature and humidity are up to 1 °C and 1 g kg^-1, respectively. The mean errors in the sensible and latent heat fluxes for the same period are approximately 15 and 30 W m^-2, respectively.Bulk transfer coefficients used for calculating heat and humidity fluxes at the surface were shown to agree rather well with the measurements, at least for the unstable data. For stable stratification, the scatter in data is generally large, and it appears that the bulk formulation chosen overestimates turbulent heat fluxes.     

Turbulent carbon exchange in very stable conditions

Turbulent fluxes obtained using the conventional eddy covariance approach result in erratic results with large time fluctuations in extremely stable conditions. This can limit efforts to estimate components of the nocturnal energy budget and respiratory CO₂ fluxes. Well-organized fluxes that show a clear dependence on turbulent intensity were obtained when multiresolution decomposition was used to estimate turbulent exchanges. CO₂, heat and water vapour fluxes were observed at a site in the eastern Amazon basin that had been cleared for agricultural purposes. Temporal scales of the carbon transfer were determined and shown to be similar to those of latent heat, but as much as three times larger than those of sensible heat. CO₂ eddy diffusivities at the temporal scales on which most of the vertical CO₂ exchange occurs are shown to be 50 times larger than the eddy diffusivity for heat. A process associated with the vertical scale of the scalar accumulation layer is suggested to explain these different scales and turbulent diffusivities of carbon and sensible heat transfer. For an appreciable range of turbulence intensities, the observed vertical turbulent carbon exchange is insufficient to account for the locally respired CO₂ estimated independently. Evidence that shallow drainage currents may account for this is given.     

Surface turbulent flux measurements over the Loess Plateau for a semi-arid climate change study

In order to provide high quality data for climate change studies, the data quality of turbulent flux measurements at the station of SACOL (Semi-Arid Climate & Environment Observatory of Lanzhou University), which is located on a semi-arid grassland over the Loess Plateau in China, has been analyzed in detail. The effects of different procedures of the flux corrections on CO₂, momentum, and latent and sensible heat fluxes were assessed. The result showed that coordinate rotation has a great influence on the momentum flux but little on scalar fluxes. For coordinate rotation using the planar fit method, different regression planes should be determined for different wind direction sectors due to the heterogeneous nature of the ground surface. Sonic temperature correction decreased the sensible heat flux by about 9%, while WPL correction (correction for density fluctuations) increased the latent heat flux by about 10%. WPL correction is also particularly important for CO₂ fluxes. Other procedures of flux corrections, such as the time delay correction and frequency response correction, do not significantly influence the turbulent fluxes. Furthermore, quality tests on stationarity and turbulence development conditions were discussed. Parameterizations of integral turbulent characteristics (ITC) were tested and a specific parameterization scheme was provided for SACOL. The ITC test on turbulence development conditions was suggested to be applied only for the vertical velocity. The combined results of the quality tests showed that about 62%–65% of the total data were of high quality for the latent heat flux and CO₂ flux, and as much as about 76% for the sensible heat flux. For the momentum flux, however, only about 35% of the data were of high quality.     

SENSITIVITY TESTS OF INTERACTION BETWEEN LAND SURFACE PHYSICAL PROCESS AND ATMOSPHERIC BOUNDARY LAYER*

In this paper,an interactive model between land surface physical process and atmosphere boundary layer is established,and is used to simulate the features of soil environmental physics,surface heat fluxes,evaporation from soil and evapotranspiration from vegetation and structures of atmosphere boundary layer over grassland underlying.The sensitivity experiments are engaged in primary physics parameters.The results show that this model can obtain reasonable simulation for diurnal variations of heat balance,soil volumetric water content,resistance of vegetation evaporation,flux of surface moisture,and profiles of turbulent exchange coefficient,turbulent momentum,potential temperature,and specific humidity.The model developed can be used to study the interaction between land surface processes and atmospheric boundary layer in city regions,and can also be used in the simulation of regional climate incorporating a mesoscale model.