均匀裸土地表土壤热通量计算方法对比及对能量闭合的影响

利用"内蒙古微气象观测蒸发试验"的观测资料,对6种地表土壤热通量计算方法(Plate Cal法、TDEC法、谐波法、热传导对流法、振幅法和相位法)进行比较,检验了6种方法在不同干湿地表状况下的适用性,并研究了6种方法计算地表土壤热通量的差异以及对地表能量闭合度的影响。结果表明:一般情况下,Plate Cal法计算的2 cm土壤热通量与观测值最接近,计算结果的均方差为6.9 W/m2。在不同干湿地表状况下,干燥和降水条件下适合使用Plate Cal法,计算结果的均方差分别为14.0 W/m2和30.1 W/m2;湿润条件下适合使用谐波法,计算结果的均方差为21.4 W/m2。6种方法计算的地表土壤热通量存在明显差别,最大相差178.6 W/m2,不同方法计算地表土壤热通量的最大差值超过25 W/m2的时次占样本的96.3%。不同方法计算地表土壤热通量的差异对地表能量闭合度的大小有明显影响,但不影响近地层能量闭合度随湍流混合增强而增大的规律。     

青藏高原典型下垫面地表能量通量的模型估算与验证

青藏高原地区地表能量通量的估算与验证对高原及其周边地区能量和水循环研究具有重要意义,地表能量平衡系统SEBS(Surface Energy Balance System)模型为研究高原非均匀地表区域地表能量通量提供了一种行之有效的方法。基于中国科学院那曲高寒气候环境观测研究站(简称那曲站)、中国科学院纳木错多圈层综合观测研究站(简称纳木错站)和中国科学院珠穆朗玛大气与环境综合观测研究站(简称珠峰站) 2008年辐射资料、大气边界层塔站观测资料,结合MODIS卫星数据,利用SEBS模型估算地表能量通量,并用站点地表能量通量观测资料进行验证。结果表明,模型估算的感热通量和土壤热通量与站点实测值具有较好的一致性,且感热通量和土壤热通量的估算精度明显优于潜热通量;感热通量的估算精度最高,那曲站、纳木错站和珠峰站的均方根误差分别为54. 98,37. 37和27. 10 W·m~(-2);而模型估算的潜热通量验证结果偏差较大和站点实测数据存在"能量不闭合"问题相关。鉴于在地表能量通量观测中广泛存在"能量不闭合"的问题,利用波文比校正方法校正站点实测潜热通量。研究表明波文比校正方法可以明显改善地表通量观测数据"能量不闭合"的问题,那曲站、纳木错站和珠峰站的能量闭合率分别提高了19. 4%,21. 4%和19. 1%;与原始站点实测潜热通量相比,校正后的潜热通量与SEBS模型估算结果一致性较好,3个站点潜热通量的均方根误差分别减少了6. 78,33. 48和29. 30 W·m~(-2)。     

陇中黄土高原地表能量不平衡特征及其影响机制研究

利用兰州大学半干旱气候与环境观测站(SACOL)2008年夏季的观测资料,分析了陇中黄土高原地表能量不平衡特征及其影响机制,计算了地表和5 cm深处土壤热通量板之间的热储存,并对地表能量不平衡特征重新进行了计算和评价。当把土壤热储存作为地表能量收支的一部分考虑后,不平衡差额绝对值平均降低了23 W.m-2,能量不闭合度平均减小了0.104,利用最小二乘法(OLS)进行线性回归得到的平均不闭合度减小了0.085。从白天和晚上的能量不闭合度频率分布也能够看出,土壤热储存对地表能量收支平衡有很大改善。但即使考虑了土壤热储存项,地表能量不闭合仍然很明显。除了土壤热储存,边界层大气不同形式的输送能力作为对地气能量通量交换过程中的重要因素,对能量不平衡也会产生深远影响。结果显示,地表能量不闭合度分别与垂直速度w、水平风速u和对流速度尺度w*有较密切的关系。涡动相关法通量计算中常通过坐标旋转强迫垂直速度变为零,然而非零的垂直速度和垂直平流是真实存在的,忽略垂直平流热量输送是产生地表能量不平衡的重要原因。     

半干旱草原下垫面能量平衡特征及土壤热通量对能量闭合率的影响

土壤热通量在半干旱草原下垫面能量平衡研究中极为重要,土壤热通量估计不够准确是导致地表能量不平衡的一个重要原因。利用2008年6—9月锡林郭勒草原主生长期地表辐射、通量和土壤温度梯度观测资料,研究中纬度半干旱草原下垫面地表能量平衡特征。首先,在分析能量平衡各分量月平均日变化特征的基础上,通过对土壤热流量板观测的5 cm深度土壤热通量(G)的相位前移,研究了土壤热通量相位滞后对地表能量平衡产生的影响;其次,利用谐波分析方法,通过计算地表土壤热通量(Gs),分析了地表到热流量板之间的土壤热量储存对地表能量平衡的影响。结果表明:(1)将土壤热通量相位前移30 min,湍流通量与可利用能量(Rn-G)线性回归的斜率从0.835增加到0.842,地表能量闭合率提高了0.7%,但仍有15.8%的能量不闭合;(2)考虑了地表到热流量板之间的土壤热量储存之后,湍流通量与可利用能量之间的回归斜率达到0.979,能量不闭合程度仅为2.1%。     

塔克拉玛干沙漠腹地地表辐射和能量平衡及小气候特征

利用塔克拉玛干沙漠大气环境综合观测试验站——塔中站2014年30 min步长实测数据,系统分析塔克拉玛干沙漠腹地地表辐射和能量平衡特征及小气候特征。结果表明:塔中地区辐射平衡和能量平衡各分量均表现为春、夏季较大,秋、冬季较小,并具有明显的日变化特征,峰值出现在12:00(地方时,下同)左右。该地区地表能量通量(净辐射通量、感热通量、潜热通量和土壤热通量)中,净辐射、感热通量和土壤热通量均在08:00开始上升,至14:00左右达到峰值,而潜热通量变化很小。感热通量为能量的主要消耗形式,春、夏、秋、冬四季分别占净辐射的百分比为44. 6%、68. 1%、55. 2%和55. 3%,其次为土壤热通量,仅有很少能量用于水的相变。云和沙尘对能量通量各分量有明显削弱作用。近地层空气温度白天随高度的增加降低,夜晚则随高度的增加升高,空气相对湿度也呈明显变化规律。近地层水平风速随高度增加明显增大,风向以东北风为主。     

半干旱区陆面过程观测试验的仪器精度和观测误差分析

利用定西试验基地2017年5月的平行对比观测资料,分析陆面过程观测试验的仪器精度和观测误差。结果表明:HMP45D型空气温湿度传感器、107-L型土壤温度传感器、WAA151型风速仪、HFP01型土壤热通量板和涡动系统(CSAT3+Li-7500)观测数据的偏差较小,均在允许范围内,说明相同仪器之间的观测数据基本一致,但观测土壤湿度的CS616型仪器在使用时需要订正;不同型号仪器观测长短波辐射的一致性很高,其相关系数均近乎1;观测数据的偏差也直接影响能量闭合率,试验中湍流通量观测偏差造成地表能量闭合率0.61%~2.20%的不确定性,可利用能量观测偏差造成地表能量闭合率0.08%~1.67%的不确定性。     

青藏高原东缘峨眉山地区冬季地表能量交换特征研究

利用青藏高原东缘峨眉山站2018年11月至2019年2月的观测数据,分析了峨眉山地区近地层气象要素的变化特征,并运用涡动相关法、土壤温度预报校正法(TDEC)和最小二乘法等,讨论了地表能量交换特征,并与藏东南丹卡站和排龙站进行对比分析。结果表明:在峨眉山地区冬季感热通量占主导。地表辐射各分量日变化均呈单峰结构,短波辐射的峰值在峨眉山地区于13:00(北京时,下同)左右出现,而长波辐射到达峰值的时间段要晚于短波辐射,在14:00左右出现。地表反照率月变化明显,日变化呈"U"型,即日出后和日落前较大,日间较小的趋势,但在日出后和日落前的值并不相等。峨眉山站冬季地表反照率均值为0.29。峨眉山区域地表能量不闭合现象十分显著,日间和夜间的能量闭合程度差异很大。在考虑地表0～5 cm处的热量储存的条件下,白天闭合程度最好为67.22%,夜间闭合程度最差为65.09%。与藏东南丹卡站和排龙站对比分析表明:丹卡、排龙站地表辐射各分量达日峰值时间晚于峨眉山站。峨眉山站冬季的地表反照率高于丹卡、排龙站,月变化更加显著。峨眉山属于青藏高原东缘地区,其地表能量闭合程度大于青藏高原上部分站点。     

基于土壤温湿度观测资料估算藏北高原地区土壤热通量

基于"全球协调加强观测计划(CEOP)亚澳季风之青藏高原试验研究"(CAMP/Tibet)土壤温湿度观测资料(2004年6月8日至7月30日),利用TDEC(土壤温度预报校正)法估算并分析了藏北高原地区土壤热通量的变化特征。通过与热通量板测量结果比较发现,它们在热流方向和相位上基本一致,且在量值上具有很好的线性关系(R~20.74),但热通量板测量值普遍偏小。通过对比不同天气条件下TDEC法的估算结果表明,云的存在以及降水发生都可导致土壤热通量的波动,且表层土壤热通量波动更加强烈。对于TDEC法估算的表层土壤热通量,发现它与净辐射通量线性关系密切(R20.78),而在土壤热传导率较大的地区,其日变化相位明显先于其他地区。对于土质均匀的砂质土壤,其土壤热通量日变化振幅最大,土壤的垂直热量交换最剧烈。     

黄土高原半干旱区异常能量闭合率特征分析

以兰州大学半干旱气候与环境观测站（简称SACOL站）4 a的陆面通量数据为基础,利用普通最小二乘法（OLS）和能量平衡比率（EBR）方法,对能量平衡的异常闭合特征及其与相对垂直湍强（RIw）的关系进行了研究,并进行了能量滞后分析。结果表明,能量过闭合和负闭合现象分别主要发生在白天和夜间,大小遵从正态分布;较大异常闭合产生的原因主要是日出日落时净辐射与地表热通量接近以及降水影响造成湍流通量出现异常大值。一般来说,垂直湍流运动越强,异常闭合越少,闭合度越向1收敛,反之亦然。强湍流或极弱湍流都不利于产生异常闭合,过闭合、负闭合的最适相对垂直湍强RIw约为0.11、0.14。另外,能量支出项的相对滞后也是造成包括负闭合在内的异常闭合的原因之一。剔除湍强较弱的点或将地表热通量G0、感热H、潜热LE相位相对净辐射Rn提前30 min后,异常闭合所占比重减少;月平均EBR法过闭合度降低,OLS法闭合度提高。     

北京郊区草地夏季能量收支平衡的数值模拟

北京郊区地表能量分配可能影响北京地区的天气和气候。为了进一步检验陆面过程模式对北京郊区具有代表性的稀疏草地地表能量分配的模拟能力,利用原版和改进版简化生物圈模式(SiB2,Simple Biosphere Model 2)模拟了2010年7月22日-8月5日期间北京郊区阳坊镇坦克打靶场草地的辐射平衡、能量收支以及地表热通量。并将模拟结果与实际测量的数据进行对比,结果表明:1)原版SiB2低估净辐射11.32%,改进版SiB2则低估净辐射5.81%,主要原因是改进版SiB2更新了土壤热传导率计算方法,从而提高了土壤温度(包括地表温度)模拟结果的精度,进而改善了地表向上的长波辐射模拟结果的准确性;2)改进版SiB2同时改善了感热通量和潜热通量的模拟结果,但是原版SiB2和改进版SiB2均低估了土壤热通量。     

Estimation of ground heat flux from soil temperature over a bare soil

Ground soil heat flux, G ₀, is a difficult-to-measure but important component of the surface energy budget. Over the past years, many methods were proposed to estimate G ₀; however, the application of these methods was seldom validated and assessed under different weather conditions. In this study, three popular models (force-restore, conduction-convection, and harmonic) and one widely used method (plate calorimetric), which had well performance in publications, were investigated using field data to estimate daily G ₀ on clear, cloudy, and rainy days, while the gradient calorimetric method was regarded as the reference for assessing the accuracy. The results showed that harmonic model was well reproducing the G ₀ curve for clear days, but it yielded large errors on cloudy and rainy days. The force-restore model worked well only under rainfall condition, but it was poor to estimate G ₀ under rain-free conditions. On the contrary, the conduction-convection model was acceptable to determine G ₀ under rain-free conditions, but it generated large errors on rainfall days. More importantly, the plate calorimetric method was the best to estimate G ₀ under different weather conditions compared with the three models, but the performance of this method is affected by the placement depth of the heat flux plate. As a result, the heat flux plate was recommended to be buried as close as possible to the surface under clear condition. But under cloudy and rainy conditions, the plate placed at depth of around 0.075 m yielded G ₀ well. Overall, the findings of this paper provide guidelines to acquire more accurate estimation of G ₀ under different weather conditions, which could improve the surface energy balance in field.     

Accuracy of soil heat flux plate measurements in coarse substrates – Field measurements versus a laboratory test

Summary The in-situ performance of heat flux plates within coarse porous substrates might be limited due to poor contact between plate and substrate. We tested this behaviour with a simple laboratory set-up. Two test substrates were placed above a reference material of known thermal conductivity between a warm and a cold plate to establish a vertical heat flux. The temperature gradients and the response of a soil heat flux plate were measured. By means of the Fourier law of heat conduction the thermal conductivity of each test substrate was calculated, thus incorporating all heat transfer within the volume and representing the “effective” conductivity. The laboratory method had an accuracy of up to ±7% (±13% for a smaller set-up). In comparison, heat flux plate-derived heat fluxes showed errors of up to 26%. Use of heat flux plates in coarse substrates is not recommended without additional measurements.     

Comparison of in-situ measured ground heat fluxes within a heteorogeneous urban ballast layer

Summary This study concentrates on measurements of ground heat fluxes within a porous urban ballast layer that were conducted from June to September 2002 at the goods station in Osnabrück, Germany. To account for the limitation of accurately installing sensors within the heterogeneous and porous ballast bulk, the heat fluxes were calculated from four different methods to compare their variability, dynamics and shortcomings. Ground heat fluxes were gathered from 1) a heat flux plate with the inclusion of heat storage between the soil surface and the heat flux plate, 2) temperature gradient measurements with correction for heat storage, 3) temperature gradient measurements with a modelled surface temperature and a laboratory derived thermal conductivity, 4) as residual from the surface energy balance equation. The results show a distinct deviation of the four methods for absolute values of the ground heat flux as well as for temporal dynamics on the diurnal cycle. As indicated by the temporal dynamics of the ground heat flux times series and a simple error analysis of the four methods, the most plausible estimates for an urban application in a heterogeneous ballast layer were obtained by temperature gradient measurements between the surface and −0.05 m. Overall, the results indicate that accurate ground heat flux measurements in urban applications still prove difficult to acquire.     

麻多高寒湿地冻结过程中土壤热通量变化特征分析

准确量化高寒湿地下垫面冻结过程中土壤热通量的变化特征,对认识高寒湿地—大气间水热交换过程有重要的科学意义。本文利用中国科学院麻多气候与环境综合观测站2014年5月至2015年5月的观测资料,分析了下垫面冻结过程中土壤热通量变化特征,探讨了冻结潜热对土壤热通量的贡献。基于温度积分计算土壤热通量的算法,指出在计算冻结过程中的土壤热通量时,需要同时考虑土壤热通量板以上的土壤热贮存及热通量板以上的冻结潜热。研究表明:（1）冻结锋面形成后,锋面所在深度土壤体积含水量迅速降低,锋面以下土壤热通量接近于零,土壤液态水开始冻结,冻结潜热向上穿过热通量板所在土壤层;降水下渗土壤后冻结所释放的潜热能使次日凌晨5 cm深度土壤热通量接近于零。（2）季节性冻结期,凌晨气温较高时穿过5 cm土壤层的向上土壤热通量很小,可能是由表层土壤发生了日冻融循环所致。土壤水释放的冻结潜热使土壤温度波动减弱并维持在冰点附近。高寒湿地下垫面仅在很浅的表层发生日冻融循环,无法通过5 cm土壤温度资料判断下垫面循环出现日期。（3）加入冻结潜热项,土壤热通量的计算值与实测值之间的均方根误差将会从11.5 W m-2下降到6.2 W m-2。以上研究结果对认识寒区陆面过程有重要的贡献。     

Evaluation of Noah-LSM for soil hydrology parameters in the Indian summer monsoon conditions

The micrometeorological observations, collected over a station in Ranchi (23°45′N, 85°30′E) which is under the monsoon trough region of India, were used in the Noah-LSM (NCEP, OSU, Air Force and Office of Hydrology Land Surface Model) to investigate the model performance in wet (2009 and 2011) and dry (2010) conditions during the south-west summer monsoon season. With this analysis, it is seen that the Noah-LSM has simulated the diurnal cycle of heat fluxes (sensible and ground) reasonably. The simulated heat fluxes were compared with its direct measurements by sonic anemometer and soil heat flux plate. The net radiation and sensible heat flux are simulated well by the model, but the simulation of ground heat flux was found to be poor in both dry as well as wet conditions. The soil temperature simulations were also found to be poor in 0–5- and 5–10-cm layers compared to other deeper layers. The observations were also correlated with the Modern Era Retrospective-analysis for Research and Applications (MERRA) data. The correlation between the observations and ground heat flux was better in MERRA dataset than that of the Noah-LSM simulation.     

白洋淀水陆不均匀地区能量平衡特征分析

应用2005年9月在河北白洋淀地区进行的大气边界层综合观测实验资料, 对水陆不均匀地表条件下的白洋淀地区陆地的能量平衡特征进行了分析。结果表明: (1) 该地区存在能量不闭合现象。涡动相关法得到的感热、 潜热之和仅为有效能的75%, 其中涡动相关法得到的潜热通量为Bowen比法得到的潜热通量的70%, 而涡动相关法得到的感热通量为Bowen比法得到的感热通量的77% 。 (2) 地表潜热通量和感热通量随着净辐射的变化而变化。但潜热通量明显比感热通量大, 净辐射主要消耗于地表的水汽蒸发。 (3) 该地区白天的Bowen比平均在-0.4~0.4之间, 总体平均为0.131。受天气条件影响较大, 有明显的日变化, 午后15:00以后近地面层会出现逆温, Bowen比变为负值。 (4) 能量闭合程度有一定的日变化, 随着太阳高度角的增大而增大。     

广州番禺地区草地陆气相互作用观测研究

介绍了2004年在广州番禺进行的陆气相互作用观测试验。观测研究表明:新型的超声风温仪虽然带有雨滴防护沙网, 但雨滴对超声观测的误差影响显著, 尤其是对u, w方向。雨天情况下与z/L相关系数很低, 湍流强度与稳定度不存在1/3次方关系。涡动相关法和能量平衡法计算的感热通量及潜热通量比较接近, 感热、潜热通量相关系数分别达到0.8699和0.8633, 两种方法带来的误差主要发生在近地层稳定度发生明显变化的时刻, 即在午间热通量的峰值或傍晚或晚间热通量的低值, 其中能量平衡法计算的感热、潜热普遍具有较大的正负峰值。涡动相关法计算的Q_h+Q_e普遍偏小, 与可用能量R_n-Q_g多数情况下存在能量不平衡, 说明了忽略热存储项的地表能量平衡方程的局限性。番禺夏、秋季近地层各能量具有与太阳辐射相似的日变化特征, 但夏季的潜热大于感热, 而秋季则相反。近地面二氧化碳从5—8月是一个减低过程, 尔后上升到12月份浓度最高, 总体浓度值在350×10-6~400×10-6之间变化。     

利用简化动力学近似分析黑河实验区的感热通量

本文利用简化动力学近似．运用统计和参数化两种方法分析了黑河实验区的感热通量．结果表明：不稳定条件下绿洲区的感热通量与近地面辐射温度(T3)和空气温度(Ta)之差符合幂定律方程；沙漠地区符合形式为H=f[u(T2-Ta)]的幂定律方程。利用从实验得到的热传输附加阻尼表示式kB^-1=0．2u(T2-Ta)用于对空气动力学阻尼(ra)的参数化，可适用于不同下垫面．并得到了比较满意的感热通量值．最后利用以上结果计算了实验区的潜热通量并做了比较分析。     

An analysis of errors in the calculation of energy flux densities above vegetation by a Bowen-ratio profile method

The Bowen-ratio profile method for calculating total carbon dioxide, latent heat, and sensible heat flux density above a corn crop was used by measuring temperature, water vapor, and CO₂ concentrations at several heights in the aerodynamic boundary layer of the crop. The ratio (α) of sensible heat flux density to carbon dioxide flux density as well as the Bowen ratio (β) were used in the computations. The two ratios, α and β were determined graphically from the slopes of linear plots of temperature vs CO₂ concentration and vs water vapor. Each of the energy flux densities was computed from the two ratios and net radiation minus soil heat flux density. An analysis of probable error was performed on the Bowen-ratio profile method to evaluate the accuracy of the flux density estimates. Less than 10% error was found for latent heat flux density and less than 15% for carbon dioxide flux density under normal midday conditions for the instrumentation used. However, the carbon dioxide flux density error increased to over 40% when the sensible heat flux was small.