放牧对退化草地近地面辐射的影响
1.
2.
3.
Effects of grazing on near surface radiation in degraded steppe in Inner Mongolia, China
1.
2.
3.
通讯作者:
收稿日期: 2021-08-17 修回日期: 2021-09-28
基金资助: |
|
Received: 2021-08-17 Revised: 2021-09-28
作者简介 About authors
李小院(1993—),男,内蒙古苏尼特右旗人,硕士研究生,主要从事水文过程与生态效应方面的研究E-mail:
关键词:
Keywords:
本文引用格式
李小院, 张圣微, 王帅, 李瑞燊, 赵星宇, 刘敏敏.
Li Xiaoyuan, Zhang Shengwei, Wang Shuai, Li Ruishen, Zhao Xingyu, Liu Minmin.
0 引言
陆地下垫面状况对气候影响的研究早在20世纪70年代末就已经开始,许多学者利用不同的大气环流模式(AGCM)模拟了陆面过程在气候变化中的作用[1],指出下垫面状况在很大程度上决定了陆表的能量和水分平衡,进而影响着局地、区域乃至全球大气环流和气候的基本特征。太阳辐射是地球能量的主要来源,是驱动地-气系统中各种物理和生命活动相互作用的原动力,也是影响其他分量变化的基本特征因子。因此,以地表辐射和热量平衡为主的能量交换过程便成为陆面过程的主要研究内容[2-4]。地气间能量交换包括发生在土壤-植被-大气之间的辐射过程,所以在研究陆气相互作用中的水、能平衡之前[5-6],对典型下垫面条件的地表辐射特征进行定位观测和深入分析非常必要[7-8]。
干旱、半干旱草原是全球陆地生态系统的重要类型,也是中国北方主要的地表类型[9],由于特殊的地表植被、土壤和反照率,其地表辐射平衡过程不同于湿润区[10-12]。放牧是草原土地利用方式中最为常见的一种,它通过减小草原地上生物量和叶面积指数,改变了植被覆盖状况,影响着地表温度、地表反射率和地表粗糙度等下垫面微气象特征[13],继而作用于地气间的能量分配、植被蒸散和CO2通量交换[14-15]。因此,禁牧、放牧对区域尺度上的碳、水循环都有着重要的影响。在生态条件脆弱的半干旱区,过度放牧是草原土壤荒漠化的重要原因[16]。而目前对于陆面过程的研究主要针对青藏高原特殊地形和西北干旱区[17],对中国北方半干旱草原区研究较少,并且已有的研究大多仅基于短时间的加密观测对草原生长期或典型天气条件下辐射特征的个例分析,利用长期的观测资料对半干旱草原下垫面辐射平衡特征的研究还稍显欠缺[18-20]。
鉴于此,本文利用内蒙古正镶白旗典型退化草地2020年生长季辐射通量观测资料,分析典型半干旱草原下垫面辐射平衡四分量(太阳总辐射Rsd,地表反射辐射Rsu,大气长波辐射Rld,地表长波辐射Rlu)、净辐射(Rn)及地表反照率(Albedo)的日变化和季节变化规律,并将禁牧、放牧草地作对比分析,探讨地表辐射收支的动态变化规律及其影响因素,有助于加深对中国北方半干旱草原下垫面辐射平衡过程的认识,为准确评估草原生态系统与大气之间物质和能量交换提供科学依据,并以期为改进北方半干旱草原区陆面过程中的辐射参数化和卫星遥感参数化过程中所需参数提供科学依据。
1 数据与方法
1.1 研究区及设备布设
研究区位于锡林郭勒草原西南部的正镶白旗,是内蒙古典型的退化草原(图1)。研究区属中温带半干旱大陆性草原气候,降水少而集中。年降水量314 mm,年蒸发量2 126 mm,年均气温1.9 ℃。生长季有天然植被覆盖,植被以禾本科和菊科植物为主,代表性植物有羊草(Leymus chinensis)、冷蒿(Artemisia frigida)、克氏针茅(Stipa krylovii)、糙隐子草(Cleistogenes squarrosa)等。
图1
试验区内分别设置了禁牧区和自由放牧区,每个区内设置3个处理和3个平行,共18个分区。在每个分区中心位置布设5、15、30、60、100 cm 5层土壤水分、温度和电导率传感器(CS650,Campbell,美国)。在禁牧和放牧区中部安设1套自动气象站进行空气温湿度(HMP155A,Campbell,美国)观测。2020年6月16日辐射计安装并成功运行,本文选择6—10月的辐射通量观测资料(W·m-2)进行分析。6月草地植被进入返青期,7—8月进入生长旺盛期,9—10月草地植被进入枯黄期。近地层定位辐射观测系统的辐射仪器(Apogee,美国)包括:在禁牧和放牧区分别安装了3套SL-610向上长波和SP-610向上短波辐射计,同时在禁牧区的对照3样地里安装1套SL-510向下长波和SP-510向下短波辐射计,用来测量不同处理下向上、向下长波、短波辐射。通过CR1000X(Campbell,美国)数据采集器对气温、湿度、降水、辐射、土壤水分、土壤温度、土壤电导率等进行每15 s一次自动观测,处理输出每1 min、5 min、10 min、15 min、0.5 h、1 h、1 d的气象、辐射和土壤数据。
1.2 气候与土壤特征
研究区的气候特征因其独特的地理位置和环境而异。图2描述了禁牧区和放牧区研究时期内(6—10月)逐日的降水量(P)、气温(Ta)、地表温度(Ts)和5 cm土层的土壤含水量(VWC)。由于研究区地处中纬度内陆地区,距海洋遥远,且海拔较高,气候变差大,是典型的干旱半干旱大陆性气候。该地区雨热同期,降水多在7、8、9月;秋季气温剧降,霜临较早。
图2
图2
禁牧、放牧草地6—10月气象因子和土壤水分特征
Fig.2
Meteorolgical factors and soil moisture change from June to October in the prohibited grazing and grazing areas
1.3 植被特征
对于现有基础上的NDVI数据集,用mice包进行数据填补。操作步骤按:观察缺失值,多重填补,分析结果优化模型,获得最终结果。首先用md.pattern( )函数来判断行的缺失值,并利用VIM包进行数据可视化。使用mice( )函数进行数据填充,填补用的方法为pmm预测均值匹配(predictive mean matching)法。多重填补法的填补矩阵数,默认为5次。获得填充结果。
对5组插补数据集进行多元线性回归分析,进行T检验和F检验,判断模型整体拟合的有效性。然后对插补的结果,再进行R2检验。最后,根据检验统计量,我们选择以第3次pmm的结果作为填充,获得填充后的完整数据集(图3)。
图3
图3
禁牧、放牧区NDVI缺失值分布散点图(A)和NDVI数据填充(B)
A图中:○.数据散点图,◆.禁牧列和放牧列的缺失值位置。B图中每一个面板为一个列的原始数据和5次填充的数据,○.原始数据取值,◆.填充的数据取值。x轴为禁牧NDVI,y轴为放牧NDVI
Fig.3
Scatter plot (A)of NDVI missing value distribution and NDVI data filling condition (B) in prohibited grazing and grazing areas
1.4 分段Piecewise结构方程模型
条件独立性检验的P值在单个Fisher’s C统计量中使用下式[25]:
这一统计量是自由度为2k的
C统计数据可用于计算SEM的AIC分数,以便在模型选择框架中进行嵌套比较[26]:
式中:K是自由度的可能性。
另一种变体AICc,可以根据样本大小进行选择:
本文采用基于R语言实现的Piecewise SEM[24],来研究辐射通量(Rn、Rlu、Rld、Rsu、Rsd、Albedo)与气候因子(P、Ta、RH、VPD)、土壤因子(Ts、VWC)以及植被因子(NDVI)之间的相互关系。
2 结果与分析
2.1 地表辐射收支的变化特征
2.1.1 禁牧和放牧草地辐射通量对比
内蒙古典型退化草原生长季地表辐射通量有显著单峰型日变化特征(图4A)。
图4
图4
禁牧、放牧草地辐射通量9—10月平均日变化(A)和逐日变化(B)
Fig.4
The average diurnal variation (A) and day-to-day variation (B) of radiation fluxes from September to October in prohibited grazing and grazing areas
太阳总辐射(Rsd)的变化特征主要受当地太阳高度角、大气透明度和云量等因素的影响[27]。研究区6—10月Rsd日平均最大值依次为779.63(13:00)、745.77(13:00)、726.03(12:00)、645.65(12:00)、635.53(12:00) W·m-2。Rsd随着时间推进逐月递减,6—8月是一年中光热条件最佳时段。
由图4A可以看出,禁牧区和放牧区各月地表反射辐射(Rsu)均呈现出典型的日变化规律。晚上由于没有Rsd致使Rsu很低,甚至为0,日出以后,Rsu开始慢慢提升,其日变化规律与Rsd相同,达到日最大值的时间也基本与Rsd一致。禁牧区6—10月日平均最大值依次为123.46(13:00)、110.42(14:00)、112.80 (12:00)、98.45(12:00)、105.30(12:00) W·m-2(图4B);而放牧区日平均最大值分别为127.73(13:00)、113.54(14:00)、126.41(12:00)、113.11(12:00)、121.02(13:00) W·m-2(图4B)。
2.1.2 地表净辐射
从昼夜变化来看,Rn呈现单峰变化趋势。日出之前净辐射通量是负值,日出之后,随着太阳高度角的增大,Rsd在逐渐升高,从而Rn也在逐渐上升,直到正午(12:00—13:00)达到一日中的最大值;而后持续下降,直到日落后Rn又变成负值。从整体变化来看,禁牧草地和放牧草地的Rn无显著差异,变化规律一致。
2.1.3 地表反照率
利用观测的Rsu与Rsd的比值得到禁牧、放牧草地的地表反照率(Albedo)。为了减小辐射较小时的测量误差,本文选取Rsd>50 W·m-2的辐射数据资料进行反照率的计算[30]。
图5
图5
生长季禁牧、放牧草地地表反照率平均日变化(A)和逐日变化(B)
Fig.5
The average diurnal variation (A) and the day-to-day variation (B) of albedo of prohibited grazing and grazing areas in growing season
总体而言,生长季放牧草地的地表反照率明显高于禁牧草地,放牧区由于自由放牧降低了植被覆盖率,使得地表反照率增大。
2.2 禁牧和放牧条件下辐射通量对气候、土壤和植被的响应
2.2.1 禁牧和放牧条件下的植被差异及其对地表反照率的影响
图6
图6
生长季禁牧、放牧草地地表反照率(Albedo)和植被覆盖度的逐月变化(A)和归一化差异植被指数(NDVI)逐日变化(B)
Fig.6
Month-to-month variation (A) of albedo and vegetation coverage and day-to-day variation (B) of NDVI of prohibited grazing and grazing areas in growing season
禁牧、放牧两种处理下的地表反照率变化相对平稳,而植被覆盖度有明显的先上升后下降的单峰变化趋势(图6A)。禁牧、放牧下总体变化趋势相一致。
研究期间,禁牧、放牧处理下的NDVI逐日变化特征以及变化率分布见图6B。在研究期间,禁牧、放牧草地NDVI整体上经历先上升后下降的变化趋势,于7月20日达到峰值。
2.2.2 禁牧和放牧条件下辐射通量对气候、土壤和植被因子的响应
由图7可知,禁牧区RH(
图7
图7
禁牧(A)、放牧(B)草地辐射通量与气候、土壤、植被因子的路径图
单箭头线称为直接路径,表示因果关系,方向由原因指向结果。实线表示正相关,虚线表示负相关,数字为路径系数
Fig.7
Path plot of radiation fluxes and climate, soil and vegetation factors in prohibited grazing (A) and grazing (B) areas. In the figure, the line with single arrow is called direct path, which indicates causality, and the direction is from cause to result. solid lines indicates positive correlation, dotted lines indicates negative correlation, and the number is the path coefficient β。The sign * indicates the significance level P value (***,0.001;**,0.01;*,0.05)
放牧区Rsu(
禁牧区Albedo和NDVI既作为内生变量受到辐射、气象和土壤因子的直接影响,又作为外生变量对Rn产生直接影响,同时Albedo对NDVI也产生直接的正向作用;放牧区Albedo作为中间变量,既受到外生变量的直接影响,又对NDVI(正向)和Rn(负向)产生直接作用。
3 讨论
3.1 生长季禁牧和放牧草地植被状态的差异性及其影响因素
6月草地植被刚进入返青期,此时研究区内植被覆盖低,NDVI处在研究期间的最低水平。随后因为雨热条件的齐备,植被迅速生长,NDVI快速上升并达到峰值[33]。这段时间内禁牧、放牧区的NDVI在6月26日出现最大差异,差异率达到0.71(图6B)。这说明在牧草萌芽期,禁牧、放牧区植被差异最为明显。此后,NDVI值处于逐渐下降趋势,直到生长季结束,其值相较于返青期略有降低。这一时期内,禁牧、放牧区的NDVI在8月5日出现第二次较大的差异,变化率达到仅次于第一次的0.69(图6B)。而此时的放牧区NDVI值也处于相对较低的水平。这说明植被生长状况良好,自由放牧导致放牧区的植被破坏较为严重,覆盖度较低,进而二者差异显著[34-35]。在8月28日至9月3日时间段,禁牧、放牧区NDVI出现一次时间短、差异小的负向变化,30日变化率达到负向最大(-0.16)。
3.2 生长季禁牧和放牧草地近地面辐射能量的差异性及其影响因素
6月由于草原植被刚进入返青期,禁牧区地表植被覆盖度略高于放牧区,从而导致Albedo略低于放牧区。10月牧草已经枯黄,但由于依然存在降水,土壤表层湿度较高,导致10月的Albedo仍然处于较低水平[38]。6—7月研究区内植物矮小,放牧区植被被牛羊啃食和践踏,地表与冬季裸露的土地无明显差异,植被盖度很低[39];而禁牧区由于围封作用,植被得以恢复,植被盖度较高,因此,这一时期内的禁牧、放牧草地植被覆盖度相差较大[40]。但此时的植物均低矮,因此Albedo相差无几。8—9月植被进入生长旺盛期,植被高而密且恢复较快,因此放牧区的植被盖度略低于禁牧区。10月植被进入枯黄期,盖度较8、9月相比显著降低。而8—10月期间,由于研究区植被以牧草为主,牛羊的啃食与破坏作用尤为明显,这一时期内的Albedo差异非常明显。
4 结论
内蒙古典型退化草原生长季地表辐射通量有显著单峰型日变化特征。研究区太阳总辐射随着时间(6—10月)推进逐月递减,6—8月是一年中光热条件最佳时段。禁牧区的地表反射辐射总体小于放牧区。研究区各月大气长波辐射日变化幅度很小,基本上没有极大值出现。地表长波辐射存在明显的日变化和季节变化特征,但是二者之间的差异甚微。
生长季放牧草地的地表反照率明显高于禁牧草地,放牧区植被覆盖率低,使得反照率增大。禁牧草地和放牧草地的净辐射无显著差异,变化规律一致。对于研究期间的植被而言,禁牧和放牧草地NDVI整体上经历先上升后下降的变化趋势,于7月20日达到峰值。期间出现两次较明显的差异时期(6月26日和8月5日)。
在禁牧区,气象因子(相对湿度和饱和水汽压差)、太阳总辐射与土壤条件(5 cm土层土壤含水量和土壤温度)均对地表反照率有重要的抑制作用。辐射分量(地表反射辐射、大气长波辐射)和植被指数(NDVI)对净辐射的影响是正向极显著的;而另外的两个辐射分量(地表反照率和地表长波辐射)对净辐射有显著的负向作用。土壤因子(5 cm土层土壤含水量)、气象因子(相对湿度)、辐射(太阳总辐射和地表反照率)均对NDVI有显著的促进作用。地表反照率(负向显著)和NDVI(正向显著)既作为内生变量受到辐射、气象和土壤因子的直接影响,又作为外生变量对净辐射产生直接影响,同时地表反照率对NDVI也产生直接显著的正向作用。
在放牧区,地表反射辐射对地表反照率的影响是正向显著的;而太阳总辐射、5 cm土层土壤含水量、相对湿度均对地表反照率的影响是负向显著的。地表反射辐射和大气长波辐射对净辐射有极显著的正向作用;而地表反照率和地表长波辐射对净辐射的作用则是负向显著的。大气长波辐射、地表反射辐射、相对湿度、5 cm土层土壤含水量和地表反照率均对NDVI有重要的促进作用;而地表长波辐射对NDVI却有着显著的抑制作用。地表反照率作为中间变量,既受到外生变量的直接影响,又对NDVI(正向显著)和净辐射(负向显著)产生直接作用。
参考文献
黑河地区绿州和沙漠地面辐射收支的若干特征
[J].
Review on estimation of land surface radiation and energy budgets from ground measurement,remote sensing and model simulations
[J].
云和降水扰动对黄土高原半干旱草地辐射收支及能量分配的影响
[J].
Grassland desertification by grazing and the resulting micrometeorological changes in Inner Mongolia
[J].
Evapotranspiration in a prairie ecosystem:effects of grazing by cattle
[J].
Effect of land management on ecosystem carbon fluxes at a subalpine grassland site in the Swiss Alps
[J].
Energy partitioning and its biophysical controls above a grazing steppe in central Mongolia
[J].
草原生长期地表辐射和能量通量月平均日变化特征
[J].
Mice:Multivariate imputation by chained equations in R
[J].
PiecewiseSEM:Piecewise structural equation modeling in R for ecology,evolution,and systematics
[J].
Confirmatory path analysis in a generalized multilevel context
[J].
The AIC model selection method applied to path analytic models compared using a d-separation test
[J].
Vegetation controls on surface energy partitioning and water budget over China
[J].
/
〈 |
|
〉 |
