巨菌草地表能量交换特征及其影响因素_刘凤山.pdf

收稿日期:20220616修回日期:20221014基金项目:福建省科技厅重大专项(2021NZ029009);福建省自然资源厅项目(KKY22003XA);福建省农业农村厅项目(KKY22001XA);福建农林大学项目(XKJC-712021030)作者简介:刘凤山(1986)，男，副研究员研究方向:地表水热收支动态Email:liufengshan0225 163com巨菌草地表能量交换特征及其影响因素刘凤山1，2，白妮妮1，3，林占熺1，林冬梅1，王子怡1，3，丁凰文2，吴宗楠2，张雅婷2，杨谨宁2，王培丹4(1福建农林大学国家菌草工程技术研究中心，福建 福州 350002;2福建农林大学生命科学学院，福建 福州 350002;3福建农林大学林学院，福建 福州 350002;4顺昌县第一中学，福建 南平 353200)摘要:基于波文比能量平衡法，研究了巨菌草生长过程入射太阳辐射(d)、反射太阳辐射(u)、净辐射(n)、潜热(LH)、显热(SH)、土壤热通量(G)和波文比()等的变化动态;通过耦合系数()评价地表能量通量的影响因素，提出一种地表通量异常值插值方法的判断标准，确定 异常值为200050结果表明:采用随机森林插值方法可获得质量良好的通量数据生育期 d、u、n、LH、SH、G、和 的平均值分别为 18599 Wm2、7256 Wm2、11135 Wm2、8294Wm2、3126 Wm2、272 Wm2、038 和 06，d是季节尺度上其他参数变化的主要影响因子，LH 和 SH 同时受到冠层导度(gc)和饱和水汽压差(VPD)等因素的影响晴天 n、LH 和 SH 较高，阴天 n、LH 和 SH 较低晴天和阴天的 由 gc和n决定，生育早期具有低 gc和高 n的特征，晴天的 大于阴天;生育盛期具有高 gc和低 n的特征，晴天的 低于阴天开放科学(资源服务)标识码(OSID)关键词:巨菌草;地表能量平衡;波文比;潜热比例中图分类号:S181文献标识码:A文章编号:1671-5470(2023)03-0383-08DOI:1013323/jcnkijfafu(natsci)20230 014Patterns and influencing factors of surface energy exchangesin Giant Juncao ecosystemLIU Fengshan1，2，BAI Nini1，3，LIN Zhanxi1，LIN Dongmei1，WANG Ziyi1，3，DING Huangwen2，WU Zongnan2，ZHANG Yating2，YANG Jinning2，WANG Peidan4(1National Engineering esearch Center of JUNCAO Technology，Fujian Agriculture and Forestry University，Fuzhou，Fujian 350002，China;2College of Life Sciences，Fujian Agriculture and Forestry University，Fuzhou，Fujian 350002，China;3College of Forestry，Fujian Agriculture and Forestry University，Fuzhou，Fujian 350002，China;4Shunchang No1 Middle School，Nanping，Fujian 353200，China)Abstract:Based on Bowen ratio-energy balance method，the incident solar radiation(d)，reflect solar radiation(u)，net radia-tion(n)，latent heat(LH)，sensible heat(SH)，soil heat flux(G)and Bowen ratio(=SH/LH)were analyzed during thegrowth of Giant Juncao(Cenchrus fungigraminus)，and the factors affecting surface energy flux were evaluated by coupling coeffi-cient()The results provided a judgment standard for the interpolation method of surface flux outlier The outliers were esti-mated to range from 200 to 050，indicating that good flux data can be obtained by random forest interpolation method The aver-ages of d，u，n，LH，SH and G during the entire growth stage were 18599，7256，11135，8294，3126，272 Wm2，with and being 038 and 06，respectively dwas the main factor affecting other parameters at the seasonal scale，and LH andSH were restricted by canopy conductance(gc)and saturated vapor pressure deficit(VPD)n，LH and SH values were relativelyhigher in the sunny days while lower in the cloudy days Moreover，the relationship of between sunny and cloudy days depended ongcand n:the early growth stage was featured by lower gcand higher nvalues，with higher in the sunny days than that in thecloudy days;at the vigorous growth stage，gcwitnessed an uprising trend while nfollowed a downward trend，with lower in thesunny daysKey words:Cenchrus fungigraminus;surface energy balance;Bowen ratio;latent heat proportion福建农林大学学报(自然科学版)第 52 卷 第 3 期Journal of Fujian Agriculture and Forestry University(Natural Science Edition)2023 年 5 月3地表能量收支特征决定了植被生态系统的生产力和功能1 不同植被具有不同的地表能量收支特征，深入了解植物地表能量收支各组分的动态具有重要意义地表能量收支是大气中热量、水分、气体等要素的主要来源，通过改变地表能量收支来改变植被类型是区域气候变化的重要驱动因素2 地表辐射管理是减缓气候变暖的重要途径之一35 气溶胶污染严重影响植物光合作用和生态系统生产力6 地表能量收支各组分包括入射太阳辐射(downward radiation，d)、反射太阳辐射(upward radiation，u)、净辐射(net radiation，n)、潜热通量(latent heat flux，LH)、显热通量(sensible heat flux，SH)、土壤热通量(ground heat flux，G)和波文比(bowen ratio，)等d是陆地生态系统的唯一能量来源;u反映植被对太阳辐射的吸收能力，受太阳高度角、叶片光谱特征、冠层结构等因素的影响78 n的主要分配形式为LH、SH 和 G9 在郁闭冠层中，n分配与植物的气孔行为紧密相关，改变植物生理过程会对 n分配产生影响1011 反映了 n在 LH 和 SH 之间的分配情况于贵瑞等1216 对陆地生态系统如森林、草地、农田和苔原等开展地表能量收支研究，对膜下滴灌14、沙漠腹地15 和绿洲16 等进行通量监测试验巨菌草(Cen-chrus fungigraminus Z X Lin D M Lin S Lan sp nov)是一种蒺藜草属植物17，具有光合效率高、生产力高和生态功能完善等优点，在食用菌培养料、牲畜饲料、生物燃料及生态治理等领域得到广泛应用1719 巨菌草的生产效率和形态学建构与地表能量收支特征密切相关9 本研究基于波文比能量平衡法研究了巨菌草地表能量交换特征及其影响因素，对于揭示巨菌草与大气之间交互作用的机理具有重要意义1材料与方法11研究区概况研究区位于宁夏石嘴山市平罗县宝丰镇鸿盛菌草科技创新产业园(3944836N，106432796E)，海拔 1 068 m该区域年平均温度 108，年平均降水量 1917 mm，年平均相对湿度 444%，属于典型的温带大陆性半干旱气候12巨菌草的生长与管理巨菌草种植在当地耕种多年的熟地，土壤类型为壤土，060 cm 土层的氮、磷、钾的平均质量分数分别为(725050)、(105058)和(2675222)gkg1，土壤含水量为 15%25%巨菌草株行距为 06 m06 m2021 年 5 月 1 日利用种茎种植，10 月 15 日收割，8 月 13 日到 10 月 15 日收集数据其间巨菌草从07 m 长到 27 m，温度 1025，平均温度 18;8 月 25 日采用漫灌灌溉 1 次，40 cm 土层含水量从2212%提高到 2873%;发生 3 次较大的降水，总降水量 264 mm13波文比通量塔的设计波文比通量塔安装在研究区中心，塔高 10 m(图 1)通量塔 5 和 10 m 高的横臂上搭载空气温度、空气湿度、风速、风向等传感器;通过 7 m 高横臂搭载的辐射传感器测量净辐射，上行、下行短波辐射等数据土壤温度和湿度探头分别埋设在 2、4、10、20 和 40 cm 深度的土壤里，土壤热通量板分别埋设在东、西方向 5cm 深度的土壤里上述仪器的型号和精度参考文献 20 14地表能量收支的测定141数据获取方式地表能量收支各组分包括 d、u、n、LH、SH、G 和 等其中，d、u、n和 G 通过安装在波文比通量塔上的仪器直接测量;、LH 和 SH 采用波文比能量平衡法获得=SH LH=T e(1)LH=nG()1+()(2)SH=nG()1+()(3)式中:为干湿表常数(hPa1);T 为空气温度();T 为两层空气的温度差();e 为两层空气的水汽压差(hPa)483福建农林大学学报(自然科学版)第 52 卷图 1波文比通量塔Fig1Bowen ratio flux tower142异常值判定规则根据波文比能量平衡法计算 LH 和 SH，趋于1 时，LH 和 SH 的计算结果趋于无穷大，从而出现通量异常值(1)0(1)(4)异常值范围可由 值决定(式 4)，位于该范围内的 LH 和 SH 都为异常值随着 取值的不同，异常值发生变化 越小，异常值范围越大异常值判断的另一个标准是 LH 和 SH 都不超过 n2021 随着 的增大，异常值的取值范围变小，异常值比例的平均值从5774%降到 1179%，且降低的速率减慢，如表 1 所示143异常值插值根据异常值判定规则确定的通量异常值，需要进行插值填补，使之成为正常值采用预测均值匹配(PMM)和随机森林(F)两种方法进行插值PMM 方法类似于回归分析，利用没有缺失的数据构建变量之间的对应关系，根据变量之间的最近原则确定缺失值，即发现缺失值;根据缺失值对应的其他变量，在其他完整数据集中寻找与该变量最接近的数据进行填补，通常选择最接近的几个数据求平均值后进行填补F 方法是一种基于决策树的机器学习算法，在训练时构建多个决策树并输出，以类(或平均预测的类)的模式进行操作22 表 1根据 决定的 异常值的范围和比例Table 1ange and ratio of Bowen ratio outiers determined by 异常值范围 异常值比例/%1 057742200 05029573150 06719384133 07514495125 0801179选择以下参数构建对应关系:d、u、n、n-G、5m 温度、10 m 温度、5 m 与 10 m 的温度差、5 m 湿度、10 m 湿度、5 m 与 10 m 的湿度差、巨菌草高度等根据 PMM 和 F 与 LH 和 SH 的对应关系利用studio 中 MICE(版本 3140)软件包对缺失值进行插补，插值次数为 3 次，迭代次数为 5 次结果表明，采用两种插值方法得到的 3 次插值数据集之间的通量平均值差异较小，2 时通量插值结果的变化率很小20 通过几个插值数据集平均值的微小差异，难以判断插值方式的优劣其他研究者采用的插值方法缺少对插值结果的定量评价14，2124 波文比能量平衡法是一种通过地表能量平衡法间接计算湍流通量，插值后的数据集最大程度保持能量之间的平衡从图 2 可知:异常值对 n与 LH、SH 和 G 之间的残差(S)有较大影响，=1 时，S 较大;2 时，S 较小(小于 5 Wm2)，且呈现随机波动的变化特征采用这 2 种插值方法得到的 3 次插值结果之间无显著差异PMM 方法的 S 平均值为(070151)Wm2，F 方法的 S 平均值为(013203)Wm2本研究选择 F 方法进行数据集插值(=2)，异常值为200 050，异常值比例为 2957%25，结果如图 2 所示残差计算公式为:S=nLHSHG(5)15数据分析为了评价地表能量通量的影响因子，使用以下公式计算空气动力学导度(ga)、植被冠层导度(gc)和解耦系数():d=(2/3)H(6)z0m=0123H(7)z0h=01z0m(8)583第 3 期刘凤山等:巨菌草地表能量交换特征及其影响因素ga=k2uzlnzmd()z0m lnzhd()z0h(9)1gc=CpVPDLH+1ga(10)=+1+gagc()(11)式中:d 为零平面位移(m);H 为植物高度(m);z0m为动量传输粗糙度(m);z0h为水热传输粗糙度(m);k 为卡曼常数(041);uz为 z 高度的风速(ms1);zm为风速的测量高度(m);zh为湿度的测量高度(m);为平均空气密度(103 kgm329);Cp为平均定压比热(1 01339 Jkg1129);VPD 为大气饱和水汽压差(kPa);为干湿表常数(0059 kPa1);为饱和水汽压斜率(kPa1)图 2不同异常值处理下插值方法对地表能量平衡的影响Fig2Influence of different gap-filling methods on surface energy balance2结果与分析d为入射短波辐射;u为反射短波辐射;n为净辐射;LH 为潜热通量;SH 为显热通量;G 为土壤热通量图 3地表能量收支的季节动态Fig3Seasonal dynamics of surface energy budget21巨菌草地表能量收支的季节动态巨菌草生育期中的地表能量收支的季节动态如图 3 所示d从 250 Wm2降到 100 Wm2，意味着天气由晴天变为阴天地表能量收支的其他分量与太阳辐射密切相关，u、n、LH 和 SH 等随着d的变化而变化例如 n从 200 Wm2下降到 50Wm2，LH 在33715058 Wm2波动，生育中期呈现相对稳定的变化趋势;SH 在11519145Wm2波动，生育后期显著增大;G 在净辐射中的总比重较低，在1017149 Wm2波动季节动态上，监测期间太阳逐渐向南半球移动，日照时数和强度均减少、降低，使得各能量收支分量呈波动下降的趋势生育期 d、u、n、LH、SH、G 和 的平均值分别为 18599 Wm2、7256 Wm2、11135Wm2、8294 Wm2、3126 Wm2、272 Wm2和 038宁夏 8、9、10 月份 d的平均值分别为 600、500 和 400 MJm2，8、9 和 10 月份的日辐射量从 2025MJm2下降到 5 MJm2左右27 该地区玉米在 59 月份的 LH 为 103118 Wm2 28，巨菌草的 LH波动范围更大，平均值略低巨菌草的 低于草地和灌木生态系统，与农田和森林接近29 683福建农林大学学报(自然科学版)第 52 卷22不同天气条件下巨菌草地表能量平衡的日动态d下降说明到达地表的太阳辐射量减少，出现阴天晴天的 d200 Wm2，阴天的 d110 Wm2图 4晴天和阴天条件下地表能量通量的密度分布Fig4Density map of surface energy fluxes onsunny and cloudy days分析不同天气条件下 10:0014:00 地表能量通量的密度分布，结果如图 4 所示从图 4 可知:晴天 n、LH 和 SH 均较高;n为 400700 Wm2，存在两个峰值，分别为 500 和 600Wm2;LH 在 300550 Wm2存在峰值;SH 在 50250 Wm2存在峰值阴天地表能量分配的各分量急剧降低，n分布在 0 150 Wm2;LH 分布在50 150 Wm2;SH 在 0 左右波动G 在晴天和阴天没有明显差异，均在 0 左右波动地表能量分配的各分量存在多个峰值，与生育期长短密切相关如晴天的 n200 Wm2，8 月的 n大于 9 月(图 3)能量通量的频率分布图包含多个峰值，每个峰值对应特定时间段内通量的分布频率其他通量存在多个峰值的原因与 n类似天气影响 n的分配过程(图 5)与阴天相比，晴天具有更高的潜热比例(LH/n)、更低的显热比例(SH/n)和土壤热通量比例(G/n)晴天，LH/n的平均值为 7297%，SH/n的平均值为 2724%，G/n的平均值为018%;阴天，LH/n的平均值为 6054%，SH/n的平均值为 3894%，G/n的平均值为412%早晨的潜热比例较低，中午达到峰值后保持不变显热比例的峰值出现在上午，下午 n更多分配到 LH阴天具有类似的规律，但是峰值时间推迟了 1 h图 5晴天和阴天条件下地表能量的分配比例Fig5Distribution of surface energy on sunny and cloudy days从能量闭合来看，晴天条件下 LH/n、SH/n和 G/n3 个比例之和为 10003%，地表能量各分量略大于 n，意味着 F 插值方法对晴天的插补效果较好;阴天这 3 个比例之和为 9536%，地表能量各分量低于nF 方法在阴天的 S 平均值为(013203)Wm2(图 5)23 的变化规律从图 6 可看出 在不同天气条件和月份的日变化过程总体上 在55 波动，夜间 在5 2 波动，白天 多在 01 波动晴天和阴天的 存在差异，8 月晴天的 大于阴天;9、10 月晴天的 小于阴天由正转负或由负转正的时间转折点也有差异晴天，8 月的 转折点为9:00 和19:00，9 月的转折点为8:00和 16:00，10 月的转折点为 9:00 和 16:00;阴天，8 月的 转折点为 8:00 和 15:00，9 月的转折点为 12:00和 18:00，10 月 多为正值，转折点不明显8 月和 9 月的 均小于 10 月783第 3 期刘凤山等:巨菌草地表能量交换特征及其影响因素夜间大气处于逆温状态24，n和 G 均较小，能量通量从地面向大气传输;日夜转换时间的 接近于1，表明此时地表能量通量的异常值较多巨菌草的 平均值为 038(图 6)，该值与水肥条件良好的草地和农田等植被2931 接近24地表能量平衡的影响因素 反映地表大气之间水热交换的耦合状况，可以解释地表能量平衡的主要影响因素3233 巨菌草 的季节特征为单峰变化过程，从 02 升高到 08 后再降到 0，峰值对应的儒略日为 260 d(图 7)在儒略日为240275 d 时 大于 05，说明地表能量通量主要受 n影响32，3536;儒略日在 240275 d 之外的 均小于 05，说明地表能量通量主要受 VPD 影响36 巨菌草 平均值为 06，该值与草地36 接近图 6 在不同月份和天气条件下的变化动态Fig6Hourly variation of Bowen ratio in different months and weather conditions巨菌草 gc的季节特征为单峰变化过程，从0005 ms1增加到0015 ms1后再下降到0，峰值的时间出现在儒略日为 255 d(图 7)巨菌草 gc的变化动态与玉米等作物37 相似晴天具有较强的太阳辐射，但巨菌草具有较低的 gc(图 7)，巨菌草将 n分配到 LH 的能力减弱，更多的 n是以 SH 形式消耗掉;而阴天巨菌草具有较低的 n，巨菌草将 n分配到 LH 的比例较高，从而造成晴天的 大于阴天(图 6)9、10 月巨菌草具有较高的 gc，巨菌草将 n分配到 LH 的能力增强，更多的 n是以 LH 形式消耗掉，从而造成晴天 小于阴天(图 6)图 7解耦系数和冠层导度的季节动态Fig7Seasonal dynamics of decoupling coefficient and canopy conductance3小结地表能量交换是植被与大气之间能量和物质传输的重要过程巨菌草独特的形态学特性改变了地表能量交换过程本研究利用波文比通量塔得到巨菌草在温带大陆性半干旱气候条件下的地表能量交换特883福建农林大学学报(自然科学版)第 52 卷征，结果表示如下(1)利用残差法可以客观评价地表能量通量插值方法的优劣，采用随机森林方法插补异常值数据可以获得质量良好的通量数据(2)受太阳辐射、巨菌草生长情况、gc和 VPD 的影响，地表能量收支各分量(d、u、n、LH、SH、G、)呈现明显的动态变化(3)天气影响地表能量收支特征，晴天的 n、LH 和 SH 均较高，阴天的 n、LH 和 SH 均较低在生育早、晚期晴天的 大于阴天，在生育盛期晴天的 小于阴天(4)巨菌草 和 gc呈单峰变化趋势，表明在生育早、晚期地表能量分配受大气因素调控，在生育盛期主要受辐射因素调控参考文献 1龚婷婷，雷慧闽，焦阳，等黄土高原农牧交错带稀疏自然植被生态系统的地表能量特征J 应用生态学报，2015，26(6):16251633 2U X，SONG H，XIA H，et al Effects of land use and land cover change on temperature in summer over the Yellow iver Ba-sin，China J emote Sensing，2022，14(17):4352 3QIAN B，JING Q，ZHANG X，et al Adapting estimation methods of daily solar radiation for crop modelling applications inCanada J Canadian Journal of Soil Science，2019，99:115 4SENEVIATNE S I，PHIPPS S J，PITMAN A J，et al Land radiative management as contributor to regional-scale climate ad-aptation and mitigation J Nature Geoscience，2018，11(2):8896 5曹龙CMIP6 地球工程模式比较计划(GeoMIP)概况与评述 J 气候变化研究进展，2019，15(5):487492 6高源，律江，李春征，等大气气溶胶对北京杨树人工林生态系统生产力的影响 J 生态学报，2022，42(12):48924902 7PANG G，CHEN D，WANG X，et al Spatiotemporal variations of land surface albedo and associated influencing factors on theTibetan Plateau J Science of the Total Environment，2022，804:150100 8DONG N，LUO M，LIU Z，et al The roles of leaf area index and albedo in vegetation induced temperature changes across Chi-na using modelling and observations J Climate Dynamics，2022，58(9):25572573 9方成圆，江洪，牛晓栋，等天目山常绿、落叶阔叶混交林生长季能量通量及平衡分析 J 福建农林大学学报(自然科学版)，2016，45(4):391397 10LIU F，TAO F，LIU J，et al Effects of land use/cover change on land surface energy partitioning and climate in NortheastChina J Theoretical and Applied Climatology，2016，123(1):141150 11刘凤山，林兴生，林辉，等巨菌草蒸散发影响机制的研究 J 福建农林大学学报(自然科学版)，2018，47(6):748754 12于贵瑞，张雷明，孙晓敏中国陆地生态系统通量观测研究网络(ChinaFLUX)的主要进展及发展展望 J 地理科学进展，2014，33(7):903917 13尹茜茜，乐旭，周浩，等全球 FLUXNET 站点总初级生产力的年际变化及其主导因子解析 J 大气科学学报，2020，43(6):11061114 14张振宇，刘丽娟，李小玉干旱区膜下滴灌向日葵农田蒸散发特征 J 中国生态农业学报(中英文)，2019，27(8):11951204 15曹寰琦，何清，金莉莉，等塔克拉玛干沙漠腹地地表能量平衡闭合特征 J 干旱气象，2020，38(6):955964 16严坤，李生宇，雷加强，等沙漠腹地人工绿地地表能量交换特征 J 干旱区地理，2013，36(3):433440 17林占熺，林冬梅，刘仲健，等基于形态和基因组证据的禾本科新种 巨菌草 J 森林与环境学报，2022，42(5):514520 18郑华坤，林雄杰，林辉，等巨菌草(Pennisetum giganteum)研究进展 J 福建农林大学学报(自然科学版)，2019，48(6):681687 19周晶，林兴生，林辉，等菌草研究与应用进展 J 福建农林大学学报(自然科学版)，2020，49(2):145152 20刘凤山，白妮妮，林占熺基于波文比系统的巨菌草蒸散发过程及影响因素J 福建农林大学学报(自然科学版)，2023，51(2):822828983第 3 期刘凤山等:巨菌草地表能量交换特征及其影响因素 21向皎，李程，张清涛，等绿洲荒漠过渡带风况对波文比和蒸散发的影响 J 生态学报，2016，36(3):705720 22NIELSEN D，BIGGA J W Miscible displacement in soils:experimental information J Soil Science Society of Amer-ica Journal，1961，25(1):15 23李传金，胡顺军，郑博文古尔班通古特沙漠南缘梭梭(Haloxylon ammodendron)群落能量平衡及蒸散特征 J 生态学报，2021，41(1):92100 24张鑫，佟玲，李思恩，等干旱区两种微气象学法测定农田蒸散发的比较研究 J 水利学报，2011，42(12):14701478 25UNLAND H E，HOUSE P，SHUTTLEWOTH W J，et al Surface flux measurement and modeling at a semi-arid SonoranDesert site J Agricultural and Forest Meteorology，1996，82(1):119153 26ALLAN，PEEIA L，SMITH M Crop evapotranspiration-guidelines for computing crop water requirements FAO Irri-gation and drainage paper 56 M ome:FAO-Food and Agriculture Organization of the United Nations，1998 27杨勤，梁旭，赵光平，等宁夏太阳辐射逐日、月、年总量的变化特征 J 干旱区研究，2009，26(3):413423 28杨志朋基于不同模型的宁夏玉米、枸杞农田蒸散发估算研究 D 石家庄:河北农业大学，2021 29黄松宇，贾昕，郑甲佳，等中国典型陆地生态系统波文比特征及影响因素 J 植物生态学报，2021，45(2):119130 30张翔，刘晓琴，张立锋，等青藏高原三江源区人工草地能量平衡的变化特征 J 生态学报，2017，37(15):49734983 31LI S G，EUGSTE W，ASANUMA J，et al Energy partitioning and its biophysical controls above a grazing steppe in centralMongolia J Agricultural and Forest Meteorology，2006，137(1):89106 32ALVES J D N，IBEIO A，ODY Y P，et al Energy balance and surface decoupling factor of a pasture in the BrazilianCerrado J Agricultural and Forest Meteorology，2022，319:108912 33JAVIS P G Coupling of transpiration to the atmosphere in horticultural crops:the omega factorJ Acta Horticulturae，1985，171:187203 34SIDHA V，BILLAH M M，VALAYAMKUNNATH P Field-scale intercomparison analysis of ecosystems in partitioning sur-face energy balance components in a semi-arid environment J Ecohydrology Hydrobiology，2019，19(1):2437 35CHANG Y，DING Y，ZHANG S，et al Dynamics and environmental controls of evapotranspiration for typical alpine meadowin the northeastern Tibetan Plateau J Journal of Hydrology，2022，612:128282 36GU S，TANG Y，CUI X，et al Energy exchange between the atmosphere and a meadow ecosystem on the Qinghai-TibetanPlateau J Agricultural and Forest Meteorology，2005，129:175185 37XU J，WU B，YU D，et al A canopy conductance model with temporal physiological and environmental factors J Scienceof the Total Environment，2021，791:148283(责任编辑:叶济蓉)093福建农林大学学报(自然科学版)第 52 卷