冰川冻土, 2023, 45(6): 1887-1896 doi: 10.7522/j.issn.1000-0240.2023.0144

冰冻圈水文与水资源

基于SHAW模型的祁连山高寒草甸土壤水分储量及通量的模拟研究

张法伟,1,2, 仪律北3, 郭小伟2, 杨永胜2, 李杰霞2, 曹广民2, 李英年1,2

1.中国科学院 三江源国家公园研究院,青海 西宁 810008

2.中国科学院 西北高原生物研究所 高原生物适应与进化重点实验室,青海 西宁 810008

3.青海省林业和草原局 林业碳汇服务中心,青海 西宁 810008

Simulation of soil water storage and flux of alpine meadow using SHAW model on the Qilian Mountains

ZHANG Fawei,1,2, YI Lübei3, GUO Xiaowei2, YANG Yongsheng2, LI Jiexia2, CAO Guangmin2, LI Yingnian1,2

1.Institute of Sanjiangyuan National Park, Chinese Academy of Sciences, Xining 810008, China

2.Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 810008, China

3.Forestry Carbon Sequestration Service Center, Qinghai Forestry and Grassland Bureau, Xining 810008, China

收稿日期: 2022-03-25   修回日期: 2022-06-01  

基金资助: 中国科学院-青海省人民政府三江源国家公园联合研究专项.  LHZX-2020-07
国家重点研发计划项目.  2017YFA0604802
国家林业和草原局经济发展研究中心业务委托项目.  JYFL-2021-00020
国家自然科学基金项目.  41730752.  41877547
青海省科技成果转化专项.  2020-SF-145

Received: 2022-03-25   Revised: 2022-06-01  

作者简介 About authors

张法伟,副研究员,主要从事高寒草地碳水过程研究.E-mail:mywing963@126.com , E-mail:mywing963@126.com

摘要

高寒草甸是祁连山国家公园青海片区的重要植被类型之一,量化该下垫面土壤水分储量及交换过程是评估区域水源涵养功能的关键科学基础。基于2017年8月1日至2018年7月31日的祁连山南麓高寒草甸生态系统水热特征的连续观测数据,对SHAW(simultaneous heat and water)模型进行了参数优化,分析了0~100 cm土壤水分储量及通量的变化特征及环境影响。结果表明:SHAW模型可以相对准确地模拟高寒草甸土壤温、湿的季节变化特征,土壤水分的模拟效果略好于土壤温度的模拟效果。日均0~100 cm土壤水分储量(SWS0-100)为(274.99±19.57,平均值±标准差)mm,5—10月植被生长季的平均SWS0-100较非生长季低21.92 mm。SWS0-100的季节变异主要受控于群落叶面积指数的正效应和蒸散发的负效应,二者通过调控浅层(0~20 cm)和中层(20~60 cm)土壤水分储量间接影响SWS0-100。日均0~100 cm土壤水分通量(SWF0-100)为(0.16±9.52) mm·d-1,表现为向下传输。生长季和非生长季的SWF0-100均值分别为3.27和-3.23 mm·d-1。降水和土壤温度梯度通过驱动浅层、中层土壤水分通量和深层(60~100 cm)土壤水分通量间接影响SWF0-100。研究结果可为科学评估祁连山高寒草甸水源涵养功能提供数据支撑和理论依据。

关键词: 高寒草甸 ; 土壤水分储量 ; 土壤水分通量 ; SHAW模型 ; 结构方程模型

Abstract

Alpine meadows are one of the most important vegetation types in Qilian Mountain National Park. Quantifying soil water storage and exchange is essential to evaluating the function of regional water conservation. Based on the continuous observational hydrothermal data from August 1, 2017, to July 31, 2018, and the parameterized SHAW (simultaneous heat and water) model, this paper aimed to explore the seasonal variations of 0-100 cm soil water storage and soil water flux and their environmental controls. The results showed that the SHAW model could simulate soil hydrothermal processes reasonably. The simulation performance of soil water content was better than that of soil temperature. The daily 0-100 cm soil water storage (SWS0-100) was 274.99±19.57 mm (Mean±S.D.), which was 21.92 mm lower in the growing season from May to October than that of the nongrowing season. The seasonal variations of SWS0-100 were mainly controlled negatively by evapotranspiration and positively by leaf area index, through the effects of shallow (0-20 cm) and middle (20-60 cm) soil moisture storage. The daily 0-100 cm soil water flux (SWF0-100) exhibited downward transmission and averaged 0.16±9.52 mm·d-1. SWF0-100 migrated downward at a rate of 3.27 mm·d-1 during the growing season and converged upward at a rate of -3.23 mm·d-1 during the nongrowing season. The seasonal variations of SWF0-100 were indirectly determined by precipitation and deep soil temperature gradients, via the effects of the shallow, middle soil water flux and the deep (60-100 cm) soil water fluxes. The research results can provide data support and a theoretical basis for the scientific evaluation of the water conservation function of alpine meadows in the Qilian Mountains.

Keywords: alpine meadow ; soil water storage ; soil water flux ; SHAW model ; structural equation model

PDF (2783KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

张法伟, 仪律北, 郭小伟, 杨永胜, 李杰霞, 曹广民, 李英年. 基于SHAW模型的祁连山高寒草甸土壤水分储量及通量的模拟研究[J]. 冰川冻土, 2023, 45(6): 1887-1896 doi:10.7522/j.issn.1000-0240.2023.0144

ZHANG Fawei, YI Lübei, GUO Xiaowei, YANG Yongsheng, LI Jiexia, CAO Guangmin, LI Yingnian. Simulation of soil water storage and flux of alpine meadow using SHAW model on the Qilian Mountains[J]. Journal of Glaciology and Geocryology, 2023, 45(6): 1887-1896 doi:10.7522/j.issn.1000-0240.2023.0144

0 引言

土壤水是土壤-植被-大气连续体的主要组成,也是四水(地表水、地下水、大气水和土壤水)转换的纽带及土壤系统内物质和能量交换的载体1,对区域水资源的形成、保持与转换具有重要意义2。祁连山国家公园地处青藏高原东北部,总面积5.02万km2,是河西走廊主要内陆河的发源地和黄河流域重要的水源产流地3。高寒草甸是祁连山国家公园青海片区重要的植被类型和区域水文调节功能发挥的载体4-6,频繁发生的夏季融化、冬季冻结的冻融循环是其地表显著的物理特征之一7-8。冻融过程导致土壤水热特性的变化会直接影响到下垫面水源涵养功能的质量和状况9-10,是高寒地区地表水循环的关键过程之一11-12,对植被群落演替及生态系统稳定性维持具有重要的影响13-14。因此,深入研究冻融过程对土壤水分储量及交换的影响机制是量化高寒草甸水源涵养功能的基础和认知区域土壤水分植被承载力的关键15-17,也是评估中下游的水资源安全的科学支撑18-19,具有重要的理论价值和现实意义。

土壤水分储量及通量是土壤水势对水分运移的综合作用结果2,受到降水、蒸散发和土壤渗漏及冻融过程的强烈影响20-21。基于涡度相关技术的研究表明高寒草甸蒸散发的季节变异主要受辐射有效能控制,和降水量基本持平22-23。深层(100 cm)土壤渗漏量约相当于年降水的5%24。高寒草甸植被生长季土壤水分储量呈现出双峰型季节变化,密切关联于土壤温度25-26,与降水的关系存在时间尺度效应1221。由于自动观测技术无法测定土壤中固态水含量,加之冻融过程的复杂性,冻融期土壤水分交换过程是目前高寒草地生态水文学研究的难点和热点之一1319。由于冻土具有较高的导热性和较低的导水性2124,在温度梯度及水分相变所导致的土壤水势变化的驱动下,土壤水分向冻结峰或融化峰汇聚27。同时,冻融期中土壤水势的温度效应十分显著28,土壤水热传输强烈耦合,进一步增加了高寒草甸土壤水分储量及通量对冻融过程响应研究的难度1927-29。基于水热传输的模型模拟成为研究冻融期水热交换过程的重要手段82730。相关研究表明,SHAW模型综合考虑了冻融期土壤水热物理传输过程,对冻融阶段土壤水分的模拟结果更接近于实际观测30,能够较好地模拟高寒草地土壤水热交换82629。但早期研究多偏重于模型的评估和改进,少有应用模型探讨高寒草地土壤水分储量及通量的变化特征。因此,以祁连山东段冷龙岭南麓的高寒草甸为对象,在SHAW模型参数本地化的基础上,研究冻融期0~100 cm土壤水分储量及通量的季节变化特征和环境调控过程,以期为祁连山国家公园高寒草甸土壤水源涵养功能评估提供数据支撑和理论依据。

1 数据与方法

1.1 研究区域和数据来源

研究区位于青海海北高寒草地生态系统国家野外科学观测研究站(简称海北站,37°37′ N,101°19′ E,海拔3 200 m)。海北站地处青藏高原东北隅祁连山东支冷龙岭南麓,为典型的高原大陆性气候,年均气温和降水分别为-1.7 ℃和580 mm,雨热同季31。植被类型为禾草-矮嵩草草甸,优势种为矮生嵩草(Kobresia humilis)、垂穗披碱草(Elymus nutans)、异针茅(Stipa aliena)、麻花艽(Gentiana straminea)和黄花棘豆(Oxytropis ochrocephala)等。植被地上净初级生产力和0~40 cm地下根系现存量分别约为400 g·m-2和3 200 g·m-2,其中约90%的根系集中在0~20 cm土层中,形成草毡层31。土壤类型为高山草甸土,发育相对年轻,粗骨性强。由于地表草毡层具有较高的导水性和较低的导热性,土壤水热特性具有显著的垂向分层特异性12

研究区地势平缓开阔,系统水热状况由相关仪器自动监测(表1)。土壤水热数据观测频率为30 min(CR800,Campbell,美国)。同时,基于涡度相关技术对下垫面与大气间的水热交换进行连续观测。涡度相关系统距离土壤水热观测点约100 m,采样频率为10 Hz。利用EddyPro 7.06(Li-Cor,美国)对10 Hz高频数据进行二次坐标选择、除趋势和WPL密度校正,输出潜热通量、显热通量、风速、风向及空气温湿等30 min平均值。同时利用CM11观测1.5 m高度处的太阳短波辐射。群落叶面积指数(LAI)来源于MODIS发布的陆地植被产品(MOD15A2)。LAI以通量塔为中心,空间和时间分辨率分别为250 m×250 m和8 d,来自于美国橡树岭国家实验室的分布式主动存档中心(Distributed Active Archive Center,Oak Ridge National Laboratory,http://daac.ornl.gov/MODIS/modis.html),用于提供SHAW模型中植被生长的LAI信息。土壤水热、潜热通量、显热通量及其他气象因子也集成为8 d的时间序列以进行后续分析。

表1   研究区系统水热状况的观测信息

Table 1  The observation information of hydrothermal feature in the study site

观测项目观测指标观测仪器埋设深度/观测高度/cm
土壤水热观测土壤温度Hydra II,Stevens,美国5、10、15、20、30、40、50、60、80和100
土壤体积含水量
涡度相关系统空气温度LI-7500A,Li-COR,美国220
空气湿度LI-7500A,Li-COR,美国220
风速、风向CSAT-3,Campbell,美国220
潜热通量LI-7500A + CSAT-3220
显热通量220
太阳短波辐射CM11, Kipp&Zonen, 荷兰150

新窗口打开| 下载CSV


1.2 模型介绍和参数率定

SHAW模型为Flerchinger和Saxton于1989年开发的一维水热耦合模型,该模型基于表面边界层的能量平衡原理及在土壤-植被-大气连续体中的水热物理传输过程,采用隐式有限差分公式表述土壤的水热交换,然后通过Newton-Raphson迭代求解20,已成为冻土区陆面过程研究的常用手段之一82629。由于观测仪器无法测定土壤中固态水含量,只有生长季土壤体积含水量的观测值才能代表土壤全部含水量,方可满足SHAW模型初始和结束时段约束数据的模拟要求,因此本研究利用2017年8月1日至9月30日和2018年6月1日至7月31日的土壤水热观测数据作为模型的初始和结束数据,以2017年8月1日至2018年7月31日的逐时气象数据[气温、相对湿度、风速、降水、新雪密度(设定为0,模型依据气温计算)和太阳辐射]作为强迫资料,进行SHAW模型模拟研究。根据研究区特征,将坡度、坡向、地表粗糙度分别设置为0°、0°和1.0 cm。0~100 cm土壤被分为11层,分别为0、5、10、15、20、30、40、50、60、80 cm和100 cm,土壤水热参数设置详见表2。其中砂粒、粉粒、黏粒、有机质含量和容重等指标来自于站点的长期监测数据32,饱和导水率、空气进入势、饱和含水量和孔径指数首先根据文献资料设定初始值829-30,然后根据SHAW模型土壤水热的拟合值与观测值的接近程度进行参数优化。干土和植被反照率分别为0.35和0.2522。植被最小气孔阻抗设置为380 s·m-1[23。由于是冬季牧场,没有考虑植被残留层,生长季的植被高度、地上生物量和LAI等参数来自观测数据,叶片宽度、丛生参数和有效根深则为经验值,分别设置为1.0、0.9 cm和45 cm。其余部分参数选取模型提供的建议值。采用SHAW 3.0.2(https://www.ars.usda.gov/pacific-west-area/boise-id/northwest-watershed-research-center/docs/shaw-model/)进行模拟研究。土壤水分储量(SWS)和水分通量(SWF)分别表示土壤中液态水与固态水的水分总量和液态水与气态水的累计交换量20,为参数优化后SHAW模型的输出结果。其中土壤固态水含量是土壤水势和温度的函数,通过设定冻结温度依据能量守恒计算而来。SWS0-100为0~100 cm土壤剖面SWS与土壤深度的积分,SWF0-100为土壤剖面SWF的累计值,以表征水分的交换强度20。依据高寒草甸植物根系和土壤水热分层特性,将0~20 cm、20~60 cm和60~100 cm依次划分为浅层、中层和深层土壤25,分别用0、5、10、15、20 cm和30、40、50、60 cm及80、100 cm的平均值表示浅层和中层及深层土壤温湿性状。浅层土壤水分储量(SWS0-20)、中层土壤水分储量(SWS20-60)和深层土壤水分储量(SWS60-100)分别是浅层、中层、深层土壤水分储量之和,浅层土壤水分通量(SWF0-20)、中层土壤水分通量(SWF20-60)及深层土壤水分通量(SWF60-100)分别浅层、中层、深层土壤的土壤水分通量之和。

表2   SHAW模型各层土壤节点参数

Table 2  The soil parameters of the SHAW model at different soil nodes

土壤深度/m砂粒含量/%粉粒含量/%黏粒含量/%有机质含量/%容重/(kg·m-3饱和导水率/(cm·h-1空气进入势/m饱和含水量/(m3·m-3孔径指数
0.0039.6055.504.9018.29750.009.50-0.030.503.50
0.0539.6055.504.9013.29850.007.80-0.550.503.50
0.1044.4549.306.259.31950.005.50-0.600.454.50
0.1549.3043.107.605.321 110.002.50-0.120.454.00
0.2034.0056.359.654.341 110.001.50-0.920.425.50
0.3018.7069.6011.703.351 120.001.00-0.420.425.20
0.4018.6368.1013.272.801 130.000.50-0.850.405.80
0.5018.5766.6014.832.251 130.000.50-1.750.357.10
0.6018.5065.1016.401.701 130.000.50-0.700.355.70
0.8022.2061.0316.771.441 130.000.50-0.850.356.50
1.0025.9056.9717.131.191 130.000.50-0.050.206.50

新窗口打开| 下载CSV


1.3 模型评估和统计方法

利用统计学参数,如水热观测值与模拟值的决定系数(R2)、均方根误差(RMSE)和平均绝对误差(MAE)评价SHAW模型模拟效果729。在8 d时间尺度上,利用结构方程模型评估环境因子对0~100 cm土壤水分储量(SWS0-100)及水分通量(SWF0-100)的驱动过程及相对强度。结构方程模型是在对整体模型拟合和判断的基础上,阐释多因子之间相互关系的一种多元统计分析方法,在近年来的相关研究中应用广泛2333。根据相关研究结果25,SWS0-100的主要影响因子包括气温(Ta)、降水(Rain)、蒸散发(LE)、LAI、SWS0-20、SWS20-60及土壤温度(Ts20-60)、SWS60-100及土壤温度(Ts60-100);SWF0-100的主要影响因子包括Rain、Ta与浅层土壤的温度梯度[Δ(Ta-Ts0)]、浅层与中层土壤温度梯度[Δ(Ts0-Ts20)]、中层与深层土壤温度梯度[Δ(Ts20-Ts60)]、SWF0-20、SWF20-60及SWF60-100。基于结构方程模型中的标准作用系数和作用路径表征环境因子对SWS0-100和SWF0-100的相对强度和驱动过程。结构方程模型在R 4.0.334平台上利用piecewiseSEM35软件包实现。缩写SWS0-20、SWS20-60、SWS60-100和SWS0-100分别是0~20、20~60、60~100 cm和0~100 cm土壤水分储量;LAI为群落叶面积指数;Rain为降水量;Ta为气温;LE为蒸散发;Ts20-60为20~60 cm土壤温度;Ts60-100为60~100 cm土壤温度;SWF0-20、SWF20-60、SWF60-100和SWF0-100分别是0~20、20~60、60~100 cm和0~100 cm土壤水分通量;Δ(Ta-Ts0)为TaTs0-20差值;Δ(Ts0-20-Ts20-60)为Ts0-20Ts20-60差值;Δ(Ts20-Ts60)为Ts20-60Ts60-100差值。

2 结果与分析

2.1 强迫资料概况

研究时段的日均Ta和累计Rain分别为-1.30 ℃和494.4 mm。日均太阳短波辐射(Swin)和相对空气湿度(RH)分别为194.24 W·m-2和63.6%,均表现出植被生长季(5—10月)较高而非生长季(11月—翌年4月)较低的季节动态(图1)。浅层、中层和深层的日均Ts分别为2.85、2.61 ℃和2.38 ℃,Ts表现出随土壤深度(x)增加而指数下降的特征(Ts=3.76×(1-x-0.10R2=0.63,P<0.01)。浅层、中层和深层的日均土壤液态水含量(SWC)分别为0.22、0.25和0.17 cm3·cm-3,除在15 cm[(0.15±0.08) cm3·cm-3]和100 cm [(0.09±0.03) cm3·cm-3]处较低外,其余各层稳定在0.22~0.27 cm3·cm-3之间。Ts和SWC的变异系数(CV)分别平均为40.1%和200.0%,均随着土壤深度增加而指数下降(R2>0.79,P<0.001)。值得说明的是,相邻深度的Ts和SWC的相关系数均大于0.95,表明土壤水热具有显著的层次相关性。

图1

图1   2017年8月1日—2018年7月31日空气温度和降水(a)、相对湿度和太阳短波辐射(b)、土壤体积含水量(c)和土壤温度(d)的季节特征

Fig. 1   Seasonal variations of air temperature (Ta) and precipitation (Rain) (a), relative humidity (RH) and solar radiation (Swin) (b), volumetric soil water content (SWC) (c), and soil temperature (Ts) (d) from August 1, 2017 to July 31, 2018


2.2 SHAW模型效果评价

逐日SWC和逐日Ts的模拟值和实测值的对比结果表明,模型对SWC的模拟效果较好,R2在0.80以上,RMSE和MAE分别小于0.05 cm3·cm-3和0.03 cm3·cm-3表3)。模型对Ts的模拟效果随土壤深度增加而降低,R2平均为0.71,RMSE和MAE分别平均为2.78 ℃和2.11 ℃。SHAW模型模拟的日均显热通量和潜热通量分别为(29.8±17.3) W·m-2和(45.2±42.9) W·m-2,较涡度相关技术观测值分别高约4.2%和1.3%。模型对潜热通量的模拟效果好于显热通量,R2分别为0.81和0.14(图2),RMSE分别为18.67 W·m-2和15.97 W·m-2

表3   SHAW模型对各节点土壤温度和土壤水分的模拟效果的统计参数

Table 3  Statistical parameters of the simulation performance of the estimated soil temperature and volumetric soil water content at each node with the SHAW model

统计参数土壤深度/m
0.000.050.100.150.200.300.400.500.600.801.00
决定系数(R2土壤温度/℃0.920.850.830.840.760.740.710.630.550.450.49
土壤水分/(cm3·cm-30.850.900.900.900.870.840.850.820.800.860.85
均方根误差(RMSE)土壤温度/℃2.882.782.742.582.912.842.792.912.892.832.44
土壤水分/(cm3·cm-30.050.030.020.040.030.040.040.040.040.010.05
平均绝对误差(MAE)土壤温度/℃2.262.132.121.942.272.102.102.192.222.131.80
土壤水分/(cm3·cm-30.030.020.010.020.020.020.020.020.020.010.03

新窗口打开| 下载CSV


图2

图2   热量通量的涡度相关技术观测值与SHAW模型模拟值的关系:潜热通量(a),显热通量(b)

Fig. 2   Relationships between the simulated heat flux with SHAW model and the observed heat flux with eddy covariance techniques: Latet heat flux (a), sensible heat flux (b)


2.3 土壤水分储量的季节特征及环境影响

日均SWS0-100为(274.99±19.57) mm,表现出生长季低(264.51 mm)而非生长季高(286.42 mm)的季节特征,最低值(224.20 mm)和最高值(298.67 mm)分别出现在6月初和1月初(图3),变异系数(CV)为7.1%。日均SWS0-20和SWS20-60分别为(59.39±9.07) mm和(134.25±12.37) mm,也表现出生长季低而非生长季高的趋势。SWS60-100平均为(81.35±5.72) mm,但表现出生长季高(84.22 mm)而非生长季低(78.22 mm)的特征。SWS0-20、SWS20-60和SWS60-100的CV分别为15.3%、9.2%和7.0%,随土壤深度增加而降低。SWS0-20和SWS20-60的CV在生长季和非生长季差别较小,分别平均为10.1%和7.4%,但生长季SWS60-100的CV为7.3%,远高于非生长季的3.8%。

图3

图3   高寒草甸不同层次土壤水分储量的季节变化

Fig. 3   Seasonal variations of soil water storage of different layers in the alpine meadow


SWS0-100的结构方程模型表明,其季节变异主要受SWS20-60和SWS0-20调控,其标准作用系数分别为0.63和0.46[图4(a)],二者可解释SWS0-100季节变异的91.4%。SWS0-20与LAI正相关,与Ta和LE负相关,但LAI、Ta和LE的作用强度基本相当。SWS20-60主要受SWS0-20TaTs20-60显著影响,其中Ta的正效应(0.80)大于Ts20-60的负效应(-0.66)。Ts60-100和SWS20-60为SWS60-100的主要影响因素,其中Ts60-100正效应的标准作用系数为0.93,远大于SWS20-60的正效应。整体上看,LE和LAI对SWS0-100的综合效应最大,分别为-0.67和0.59,LE主要通过对SWS0-20(-0.29)和SWS20-60(-0.38)的负效应而LAI则通过对SWS0-20(0.25)和SWS20-60的正效应(0.34)影响SWS0-100Ta对SWS0-100的整体效应为-0.12。值得注意的是,在8 d时间尺度上,研究时段内降水对SWS的季节变异无显著影响。

图4

图4   高寒草甸土壤水分储量(a)和土壤水分通量(b)的结构方程模型

注:实线表示显著正效应,虚线表示显著负效应,点线表示不显著。线段上的数字为标准作用系数,线段的粗细为标准作用系数的10倍。显著性水平用ns表示(P>0.05),*(P<0.05)、**(P<0.01)和***(P<0.001)。

Fig. 4   Piecewise SEM models for soil water storage (SWS) (a) and soil water flux (SWF) (b) in the alpine meadow

Note:Solid (dashed) lines showed positive (negative) relationships;dotted lines showed that relationships were not statistically significant. The numbers that annotate the arrows were the standardized path coefficients,the weight for each line was ten times the standardized coefficient. With significance levels indicated by ns (P>0.05),*(P<0.05),**(P<0.01) and ***(P<0.001).


2.4 土壤水分通量的季节特征及环境调控

日均SWF0-20、SWF20-60和SWF60-100分别为(0.92±4.88)、(-0.34±4.09) mm·d-1和(-0.19±1.81) mm·d-1,都表现出生长季水分通量为正而非生长季为负的季节特征,即土壤水分在生长季向下运移,非生长季向上汇聚(图5)。生长季中,SWF0-20 (2.85±5.88 mm·d-1)远高于SWF20-60 (0.48±5.08 mm·d-1)和SWF60-100 (0.30±4.34 mm·d-1)。非生长季中,SWF0-20、SWF20-60和SWF60-100差异较小,分别为(-1.18±2.10)、(-1.24±2.44) mm·d-1和(-0.73±2.24) mm·d-1,三者平均为(-1.05±2.26) mm·d-1。冻结期各层土壤水分通量基本为0 mm·d-1,但不同层次存在时间滞后,滞后速率约为0.56 d·cm-1。日均SWF0-100为(0.16±9.52) mm·d-1,在生长季和非生长季中分别为(3.27±11.74) mm·d-1和(-3.23±4.50) mm·d-1

图5

图5   高寒草甸不同层次土壤水分通量的季节变化

Fig. 5   Seasonal variations of soil water flux of different layers in the alpine meadow


SWF0-100的结构方程模型表明,SWF0-20、SWF2060-和SWF60-100对SWF0-100的标准作用系数分别为0.43、0.36和0.32,三者作用强度基本相当[图4(b)]。SWF0-20对SWF20-60、SWF20-60对SWF60-100具有显著的正效应。SWF0-20的季节变异主要受Rain驱动,其标准作用系数高达0.95。SWF20-60主要受浅层和深层温度梯度影响,二者标准作用系数分别为0.56和-0.46。SWF60-100主要受中层、深层土壤温度梯度影响,二者标准作用系数分别为0.62和-0.40。降水对SWF0-100的总体作用系数为0.99,主要通过对SWF0-20 (0.40)和SWF20-60 (0.35)的间接效应来影响SWF0-100。浅层、中层和深层温度梯度的标准作用系数分别为0.38、0.20和-0.44。即SWF0-100季节变异的主要影响因子是降水,其次是深层土壤温度梯度。

3 讨论

3.1 SHAW模型模拟效果评价

本研究SHAW模型对高寒草甸土壤水热变化的模拟效果相对较好,与该模型在其他相关研究中表现(R2>0.70)较为一致72629,表明SHAW模型适合于高寒草甸0~100 cm土壤水热传输的模拟,模拟结果能够用以评估土壤水分状况。然而,本研究SHAW模型对SWC的模拟效果略好于Ts的模拟效果,这与郭林茂等8研究结果不同。一方面,土壤水分参数的不确定性是陆面过程模型的重要限制之一2629,本研究的一部分土壤水分参数(如容重、机械组成和有机质)是基于地面实测数据而来,另一部分难以直接测定的参数(如孔径指数和空气进入势等)是通过降低模型模拟结果和观测数据之间的残差优化获取的,有效地提高了SHAW模型对SWC的模拟效果7。另一方面,Ts模拟效果略差主要是由于地表(0 cm)Ts是按照气象观测规范所测定的裸露地表Ts,并非有植被覆盖的真实地表Ts,尤其在植被生长季中二者的差异较大36。而SHAW模型是按照土壤节点顺序进行土壤水热模拟,在保障表层Ts模拟效果的基础上,对深层Ts进行线性插补20,导致深层Ts的模拟值与观测值差别较大,即Ts的模拟效果随着土壤深度增加而下降,这和唐古拉地区SHAW模型对土壤剖面Ts的模拟效果相似26Ts的模拟误差也影响了地气温度梯度及显热通量的模拟效果,这与西藏那曲29和唐古拉地区7的相关研究结论一致。

3.2 土壤水分储量

由于研究区100 cm下土壤渗漏量很少24,SWS0-100主要取决于Rain和LE的平衡。SWS0-100在生长季低而非生长季较高,表明土壤冻结有利于土壤系统水分的保持1019,这由于冻土具有较高的隔水特性,限制了非生长季LE的过度损失8,加之冻结过程对下层土壤水分具有汇聚作用12。因为研究区在6月初的太阳辐射较强而LAI较小22,系统获取的辐射有效能较大,加之冻土消融提高了表层土壤水分含量,刺激了LE的迅速增加23,导致SWS0-100出现低值。

SWS0-100的季节变异主要受SWS0-20和SWS20-60的影响,这主要与浅层和中层水分变化较为活跃,深层水分含量相对稳定有关25。SWS0-20受LE、Ta和LAI的综合影响[图4(a)]。LE对SWS0-20为负效应,表明LE是SWS0-20散失的主要途径22-23Ta的负效应主要由于Ta可以表征大气蒸发需求,且与LE在季节变化上也具有显著的正相关22。LAI越大,植被蒸腾在LE中占比越高,系统可通过调节植被气孔行为控制水分散失22,加之高寒植被反照率较湿润裸地高,可以降低系统的辐射有效能及水分的过度消耗23,有利于浅层土壤水分保持。SWS20-60主要受SWS0-20的正效应影响,因为二者水分储量的季节变化具有一致性(R2=0.50,P<0.01)。由于SWS60-100的季节变异主要受控于生长季SWS60-100的变异,而后者与生长季Ts60-100显著正相关(R2=0.73,P<0.001),进而导致SWS60-100的季节变异与Ts60-100正相关,这也与相邻区域高寒灌丛草甸的研究结果一致21。降水对SWS0-100在8 d时间尺度上的季节变异无显著作用,一方面由于年降水的85%集中在植被生长季,且以小降水(<3.0 mm·d-1)事件为主,加之高原太阳辐射较强,极易被LE迅速消耗21-22。另一方面,强降水事件[(31.4 mm)·(48 h)-1]最多只能作用到高寒草甸40 cm的SWC2124,也一定程度上削弱了降水对SWS0-100的影响强度。

3.3 土壤水分通量

SWF的实质是各种势能综合作用下的水分质量的转移,主要包括温度梯度驱动的未冻水迁移、重力作用下的自由水下渗和土壤毛细作用下的毛细水迁移19。生长季中,SWF主要驱动力是由降水、蒸散发等过程引起的土壤含水量增减而导致的土壤剖面基质势梯度的变化,土壤水迁移以自由水和毛细水为主27,降水对SWF的影响也随着土壤深度增加而下降8,本研究也发现SWF20-60和SWF60-100的变化主要来源于降水的间接作用[图4(b)]。这主要由于高寒草甸植物根系集中在浅层土壤中形成了一系列大孔隙网络,提高了土壤的导水率12。非生长季中,土壤与大气间的水气交换被冻结的表层土壤所阻隔,土壤水分迁移与冻土融冻导致的温度变化密切相关1027,即温度梯度是SWF的主要影响因素[图4(b)]。一方面,冻结过程中,土壤水分在温度梯度的作用下向冻结峰汇集,土壤水分从深层向浅层聚集27,从而导致了Δ(Ts20-60-Ts60-100)对SWF20-60和SWF60-100的负效应。消融过程中,由于季节冻土的双向消融,加之浅层融化速率较快1012,在浅层土壤中形成一个水分的高值区,在温度梯度、重力势和基质势的综合作用下,土壤水分向下运移10,导致Δ(Ta-Ts0-20)和Δ(Ts0-20-Ts20-60)对SWF20-60和SWF60-100产生了正效应[图4(b)]。另一方面,温度是非生长季中土壤水分相变、未冻水含量及土水势的重要影响因素,还可通过改变土壤的导水率和水黏滞系数28来影响土壤水分通量。

4 结论

利用祁连山南麓高寒草甸2017年8月1日—2018年7月31观测资料作为SHAW模型的强迫资料,进行了冻融过程中土壤水热交换的模拟研究,并应用模型结果分析了季节冻结层SWS及SWF的变化特征及环境驱动,得到如下主要结论:

(1)日均SWS0-100为274.9 mm,约50%集中在20~60 cm中层土壤,生长季和非生长季的日均SWS0-100分别为264.5和286.4 mm。SWS0-100最低值和最高值分别出现在6月初和1月初,季节变异较小。蒸散发和群落叶面积指数通过影响SWS20-60和SWS0-20驱动SWS0-100的季节变化。

(2)日均SWF0-100为0.16 mm·d-1,生长季和非生长季的分别为3.27和-3.23 mm·d-1,表现为土壤水分在生长季向下传输而非生长季向上汇聚。生长季浅层水分交换活跃,SWF0-20为2.85 mm·d-1,而非生长季中浅层、中层和深层差异较小,SWF平均为-1.05 mm·d-1。降水和深层土壤温度梯度通过影响SWF0-20、SWF20-60和SWF60-100驱动SWF0-100的季节变化。

参考文献

Wang GenxuZhang ZhiqiangLi Xiaoyanet al. Introduction to ecohydrology[M]. BeijingScience Press2020.

[本文引用: 1]

王根绪张志强李小雁. 生态水文学概论[M]. 北京科学出版社2020.

[本文引用: 1]

Chapin FMatson PMooney H. Principles of terrestrial ecosystem ecology[M]. 2nd ed. New York, USASpringer-Verlag2011142-145.

[本文引用: 2]

Fu JianxinCao GuangchaoGuo Wenjiong.

Spatial-temporal differentiation of mountain-water-forest-farmland-lake-grass system in Qinghai area of the Qilian Mountain National Park, China

[J]. Chinese Journal of Applied Ecology, 2021328): 2866-2874.

[本文引用: 1]

付建新曹广超郭文炯.

祁连山国家公园青海片区山水林田湖草的时空分异

[J]. 应用生态学报, 2021328): 2866-2874.

[本文引用: 1]

Zhou Xingmin. Kobresia meadow in China[M]. BeijingScience Press2001.

[本文引用: 1]

周兴民. 中国嵩草草甸[M]. 北京科学出版社2001.

[本文引用: 1]

Li HongqinWu XialuZhang Faweiet al.

The response of soil water holding capacity of alpine meadow in winter pasture to enclosure

[J]. Chinese Journal of Grassland, 2021438): 1-8.

李红琴吴夏璐张法伟.

高寒草甸冬季牧场土壤持水能力对草地封育的响应

[J]. 中国草地学报, 2021438): 1-8.

Xue JianLi ZongxingFeng Qiet al.

Spatiotemporal variation characteristics of water conservation amount in the Qilian Mountains from 1980 to 2017

[J]. Journal of Glaciology and Geocryology, 2022441): 1-13.

[本文引用: 1]

薛健李宗省冯起.

1980—2017年祁连山水源涵养量时空变化特征

[J]. 冰川冻土, 2022441): 1-13.

[本文引用: 1]

Liu YangZhao LinLi Ren.

Simulation of the soil water-thermal features within the active layer in Tanggula region, Tibetan Plateau, by using SHAW model

[J]. Journal of Glaciology and Geocryology, 2013352): 280-290.

[本文引用: 5]

刘杨赵林李韧.

基于SHAW模型的青藏高原唐古拉地区活动层土壤水热特征模拟

[J]. 冰川冻土, 2013352): 280-290.

[本文引用: 5]

Guo LinmaoChang JuanZhou Jianet al.

Analysis of thermal-moisture conditions of active layer and energy-water balance of land-atmosphere system in Tanggula area

[J]. Plateau Meteorology, 2020392): 254-265.

[本文引用: 8]

郭林茂常娟周剑.

唐古拉地区活动层水热状况及地气系统能水平衡分析

[J]. 高原气象, 2020392): 254-265.

[本文引用: 8]

Yang MeixueYao TandongGou Xiaohuaet al.

The soil moisture distribution, thawing-freezing processes and their effects on the seasonal transition on the Qinghai-Xizang (Tibetan) Plateau

[J]. Journal of Asian Earth Sciences, 2003215): 457-465.

[本文引用: 1]

Dai LicongKe XunZhang Faweiet al.

Characteristics of hydro-thermal coupling during soil freezing-thawing process in seasonally frozen soil regions on the Tibetan Plateau

[J]. Journal of Glaciology and Geocryology, 2020422): 390-398.

[本文引用: 5]

戴黎聪柯浔张法伟.

青藏高原季节冻土区土壤冻融过程水热耦合特征

[J]. 冰川冻土, 2020422): 390-398.

[本文引用: 5]

Gu SongTang YanhongCui Xiaoyonget al.

Energy exchange between the atmosphere and a meadow ecosystem on the Qinghai-Tibetan Plateau

[J]. Agricultural and Forest Meteorology, 2005129175-185.

[本文引用: 1]

Zhang FaweiLi HongqinLi Yikanget al.

Strong seasonal connectivity between shallow groundwater and soil frost in a humid alpine meadow, northeastern Qinghai-Tibetan Plateau

[J]. Journal of Hydrology, 2019574926-935.

[本文引用: 6]

Wang GenxuYang YanZhang Guangtaoet al.

Cryosphere ecosystems: outpost and barrier in global change

[J]. Bulletin of Chinese Academy of Sciences, 2020354): 425-433.

[本文引用: 2]

王根绪杨燕张光涛.

冰冻圈生态系统:全球变化的前哨与屏障

[J]. 中国科学院院刊, 2020354): 425-433.

[本文引用: 2]

Zhou JianWang GenxuLi Xinet al.

Energy-water balance of meadow ecosystem in cold frozen soil areas

[J]. Journal of Glaciology and Geocryology, 2008303): 398-407.

[本文引用: 1]

周剑王根绪李新.

高寒冻土地区草甸草地生态系统的能量-水分平衡分析

[J]. 冰川冻土, 2008303): 398-407.

[本文引用: 1]

Liu BingxiaShao Ming’an.

Modeling soil-water dynamics and soil-water carrying capacity for vegetation on the Loess Plateau, China

[J]. Agricultural Water Management, 2015159176-184.

[本文引用: 1]

Wang GenxuXia JunLi Xiaoyanet al.

Critical advances in understanding ecohydrological processes of terrestrial vegetation: from leaf to watershed scale

[J]. Chinese Science Bulletin, 202166(): 3667-3683.

王根绪夏军李小雁.

陆地植被生态水文过程前沿进展: 从植物叶片到流域

[J]. 科学通报, 202166(): 3667-3683.

Zhou GuoyiXia JunZhou Pinget al.

Not vegetation itself but mis-revegetation reduces water resources

[J]. Scientia Sinica: Terrae, 2021512): 175-182.

[本文引用: 1]

周国逸夏军周平.

不恰当的植被恢复导致水资源减少

[J]. 中国科学: 地球科学, 2021512): 175-182.

[本文引用: 1]

Fu BojieOuyang ZhiyunShi Penget al.

Current condition and protection strategies of Qinghai-Tibet Plateau ecological security barrier

[J]. Bulletin of Chinese Academy of Sciences, 20213611): 1298-1306.

[本文引用: 1]

傅伯杰欧阳志云施鹏.

青藏高原生态安全屏障状况与保护对策

[J]. 中国科学院院刊, 20213611): 1298-1306.

[本文引用: 1]

Zhao LinHu GuojieZou Defuet al.

Permafrost changes and its effects on hydrological processes on Qinghai-Tibet Plateau

[J]. Bulletin of Chinese Academy of Sciences, 20193411): 1233-1246.

[本文引用: 5]

赵林胡国杰邹德富.

青藏高原多年冻土变化对水文过程的影响

[J]. 中国科学院院刊, 20193411): 1233-1246.

[本文引用: 5]

Flerchinger G NCaldwell T GCho Jet al.

Simultaneous heat and water (SHAW) model: model use, calibration, and validation

[J]. Transactions of the ASABE, 2012554): 1395-1411.

[本文引用: 5]

Li JingZhang FaweiSi Mengkeet al.

Response of soil water storage to meteorological factors in alpine shrub meadow on Northeastern Qinghai-Tibetan Plateau

[J]. Diversity, 2022143): 185.

[本文引用: 6]

Zhang FaweiLi HongqinWang Wenyinget al.

Net radiation rather than surface moisture limits evapotranspiration over a humid alpine meadow on the northeastern Qinghai-Tibetan Plateau

[J]. Ecohydrology, 2018112): e1925.

[本文引用: 7]

Li HongqinZhang FaweiZhu Jingbinet al.

Precipitation rather than evapotranspiration determines the warm-season water supply in an alpine shrub and an alpine meadow

[J]. Agricultural and Forest Meteorology, 2021300108318.

[本文引用: 6]

Dai LicongGuo XiaoweiZhang Faweiet al.

Seasonal dynamics and controls of deep soil water infiltration in the seasonally-frozen region of the Qinghai-Tibet Plateau

[J]. Journal of Hydrology, 2019571740-748.

[本文引用: 4]

Li HongqinWang ZhuoquanZhang Faweiet al.

Soil water conservation function of an alpine graminoid-Kobresia meadow on the southern slope of the Qilian Mountains

[J]. Research of Soil and Water Conservation, 2022293): 135-141 146.

[本文引用: 4]

李红琴王卓权张法伟.

祁连山南麓高寒禾草-矮嵩草草甸土壤水源涵养功能的特征

[J]. 水土保持研究, 2022293): 135-141 146.

[本文引用: 4]

Zhao LinLi RenDing Yongjian.

Simulation on the soil water-thermal characteristics of the active layer in Tanggula range

[J]. Journal of Glaciology and Geocryology, 2008306): 930-937.

[本文引用: 6]

赵林李韧丁永建.

唐古拉地区活动层土壤水热特征的模拟研究

[J]. 冰川冻土, 2008306): 930-937.

[本文引用: 6]

Zhao LinCheng GuodongLi Shuxunet al.

Freezing and melting process of permafrost active layer near Wudaoliang in Qinghai-Tibet Plateau

[J]. Chinese Science Bulletin, 20004511): 1205-1211.

[本文引用: 6]

赵林程国栋李述训.

青藏高原五道梁附近多年冻土活动层冻结和融化过程

[J]. 科学通报, 20004511): 1205-1211.

[本文引用: 6]

Zhang FaweiWang JunbangLin Liet al.

Temporal variations of soil surface resistance to vapor transfer and its quantitative relationship between soil temperature and soil moisture during non-growing season on an alpine meadow

[J]. Chinese Journal of Agrometeorology, 2017382): 96-103.

[本文引用: 2]

张法伟王军邦林丽.

高寒草甸非生长季土壤表层水汽传输阻抗的变化特征和水热驱动

[J]. 中国农业气象, 2017382): 96-103.

[本文引用: 2]

Guo DonglinYang Meixue.

Simulation of soil temperature and moisture in seasonally frozen ground of central Tibetan Plateau by SHAW model

[J]. Plateau Meteorology, 2010296): 1369-1377.

[本文引用: 8]

郭东林杨梅学.

SHAW模式对青藏高原中部季节冻土区土壤温、湿度的模拟

[J]. 高原气象, 2010296): 1369-1377.

[本文引用: 8]

Yang YongChen RenshengYe Baishenget al.

Heat and water transfer processes on the typical underlying surfaces of frozen soil in cold regions(Ⅰ): model comparison

[J]. Journal of Glaciology and Geocryology, 2013356): 1545-1554.

[本文引用: 3]

阳勇陈仁升叶柏生.

寒区典型下垫面冻土水热过程对比研究(Ⅰ): 模型对比

[J]. 冰川冻土, 2013356): 1545-1554.

[本文引用: 3]

Li HongqinZhang FaweiLi Yingnianet al.

Thirty-year variations of above-ground net primary production and precipitation-use efficiency of an alpine meadow in the north-eastern Qinghai-Tibetan Plateau

[J]. Grass and Forage Science, 2015712): 208-218.

[本文引用: 2]

Cao GuangminZhang Fawei. Chinese ecosystem observation and research dataset: grassland and desert ecosystem volume, Qinghai Haibei (2004—2015)[M]. BeijingChina Agriculture Press2021.

[本文引用: 1]

曹广民张法伟. 中国生态系统定位观测与研究数据集: 草地与荒漠生态系统卷-青海海北站(2004—2015)[M]. 北京中国农业出版社2021.

[本文引用: 1]

Li HongqinZhu JingbinZhang Faweiet al.

Growth stage-dependant variability in water vapor and CO2 exchanges over a humid alpine shrubland on the northeastern Qinghai-Tibetan Plateau

[J]. Agricultural and Forest Meteorology, 201926855-62.

[本文引用: 1]

Development Core Team R. R: A language and environment for statistical computing[DB/CD]. ViennaR Foundation for Statistical Computing2006.

[本文引用: 1]

Lefcheck J S.

piecewiseSEM: Piecewise structural equationmodelling in R for ecology, evolution, and systematics

[J]. Methods in Ecology and Evolution, 201675): 573-579.

[本文引用: 1]

Zhang FaweiGuo ZhujunLi Yikanget al.

Effect of soil temperature and volumetric moisture content on vegetation degeneration over the Achnatherum Splendens steppe around the Qinghai Lake, China

[J]. Arid Zone Research, 2013302): 219-225.

[本文引用: 1]

张法伟郭竹筠李以康.

青海湖芨芨草干草原植被退化对土壤温湿特征的影响

[J]. 干旱区研究, 2013302): 219-225.

[本文引用: 1]

/