焦其顺(1988-),男,河北邢台人,硕士研究生,主要从事中尺度土壤水分观测研究.E-mail: dream2qsj@gmail.com
区域地表土壤水分的获取是长期以来的研究难点。近年来, 利用宇宙射线快中子法测量区域地表土壤含水量的研究得到了长足发展。作为一种新型的中尺度土壤水分无损测量方法, 宇宙射线快中子法能够获得半径300 m左右源区内的土壤水分, 较大的测量源区使其填补了传统点观测和遥感大面积观测尺度之间的空白。利用宇宙射线土壤水分观测系统(COSMOS)、无线传感器网络对张掖绿洲农田下垫面的土壤水分进行观测。研究结果表明COSMOS观测数据很好地反映了农田尺度下土壤水分的变化趋势。COSMOS土壤水分主要受源区内灌溉和降水的影响, 源区内不同地区灌溉顺序的不同使观测的COSMOS土壤水分呈现明显的双峰变化趋势。利用COSMOS源区内无线传感器网络数据对观测结果进行了验证, 灌溉期受地面积水的影响, COSMOS观测结果偏高, 剔除灌溉期数据后两者一致性较好, 均方根误差(RMSE)由0.037 m3/m3降低到0.028 m3/m3。宇宙射线土壤水分观测系统能够有效测量高异质性状况下的农田区域土壤水分, 为遥感反演的土壤水分提供真正意义上的面状地面验证数
It is very difficult to measure large region soil moisture. In recent years, measurement of surface soil moisture by the cosmicray fast neutron probe has gradually attracted more attention. By this passive, noninvasive and intermediate scale measurement, soil moisture at a horizontal scale of around 300 m can be observed, which makes this method available to fill the gap between little scale of traditional point measurement and large scale of remote sensing in the measurement of soil moisture. In this paper, Cosmic-ray Soil Moisture Observing System (COSMOS) and Wireless Sensor Network (WSN) were used to observe field soil moisture in Zhangye Oasis. The results of the COSMOS soil moisture well reflected the variation trend of soil moisture at the field scale. There were some regular changes with the cosmicray soil moisture during the irrigation period. Combined with irrigation data in cosmicray probe footprint, a bimodal change was showed in the cosmicray soil moisture figure during irrigation period. This was because that the order of irrigation of the three communities was different in cosmicray probe footprint. WSN data were used to validate this result and we found that rootmeansquare error between cosmicray soil moisture and SoilNET average soil moisture was very large during irrigation periods because of the impact of cosmicray measurement theory. While rootmeansquare error would decrease from 0.037 m3/m3 to 0.028 m3/m3 after we eliminated the data in irrigation period. So COSMOS can be used to measure field soil moisture under high heterogeneity condition and provide truly ground data for the validation of remote sensing
土壤水作为水资源的重要组成部分, 在多种尺度的水分运动过程中扮演着重要的角色。长期以来, 如何准确测量土壤含水量一直是相关学科研究的重点。到目前为止, 国内外科研人员提出的土壤水分测定方法多达几十种[ 1]。传统的测量方法包括:称重法(烘干法)、中子仪法、电容法、张力计法、时域反射计(Time-domain Reflectometry, TDR)[ 2]、频域反射计(Frequency-domain Reflectometry, FDR)等, 这些测量方法都属于有损、点测量方式, 其测量尺度多为小尺度, 空间代表性差, 难以得到土壤含水量在区域分布上的信息[ 3]。遥感监测土壤水分的可行性研究始于20世纪60年代, 应用研究也在70年代中期开展起来。利用遥感技术可以实现快速、大面积区域土壤水分监测[ 4][Vinnikov, 1999 #2960], 但遥感获取的土壤含水量是像元尺度的平均值[ 5], 其空间分辨率一般在数十米(航空遥感)到数十公里(航天遥感)之间, 且只是飞机或卫星过境时刻的瞬时值, 无法获取连续数据, 限制了遥感大尺度测量在实际土壤水分监测中的应用。
传统点测量和遥感大尺度监测之间存在一个缺口, 即无法获取农田或小流域尺度上的区域土壤水分信息。同样, 在验证遥感反演的土壤水分过程中无法获取所需的像元尺度土壤水分数据[ 6]。近几年来, 利用宇宙射线快中子测量地表土壤含水量的研究逐渐引起人们的重视, 地表以上宇宙射线快中子强度与土壤水分含量呈反比关系, 通过架设在地表上方的中子探头测量宇宙射线快中子的强度, 从而反演地表土壤水分含量[ 7, 8]。宇宙射线快中子法是一种被动、非侵入式的中尺度土壤水分测量方法, 其测量半径在300 m左右, 有效测量深度为12~76 cm[ 7]。本文利用宇宙射线仪器对黑河中游张掖地区农田土壤水分进行长时间序列连续监测, 并与土壤水分无线传感器数据(Wireless Sensor Network, WSN)[ 9]进行了对比分析, 以研究宇宙射线观测结果的可靠性。
本文使用的数据集均来自2012年“黑河流域生态-水文综合遥感观测试验”[ 10, 11], 主要包括宇宙射线土壤水分数据集、无线传感器网络数据集(SoilNET部分)[ 9]和降雨数据集(大满站[ 12])。
宇宙射线土壤水分数据集主要包括快中子原始数据、气压原始数据等, 采样频率为1小时, 利用快中子数据可以用来反演研究区的土壤水分信息;无线传感器网络数据集(SoilNET部分)主要包含位于大满超级站所在的962 m×962 m MODIS像元内的50个SoilNET土壤水分/温度节点观测数据, 每个节点可以连接4个SISOMOP(SImple SOil MOisture Probe, University of Karlsruhe, Germany)土壤水分/温度探头[ 13, 14], 分别探测4 cm、10 cm、20 cm和40 ㎝土壤水分和温度, 其中19个分布在宇宙射线源区内。SoilNET节点的采样频率为10分钟, 文中主要使用SoilNET节点土壤水分数据对宇宙射线观测结果进行验证, 使用时需将原始10分钟数据聚合为1小时数据, 以与宇宙射线数据同步;降雨数据集(大满超级站)主要包含原始降雨数据, 采样频率为10分钟, 使用过程中需将原始数据聚合为1小时数据, 以便于分析宇宙射线土壤水分变化与降雨之间的关系。
太空中存在大量的高能粒子流, 受地球磁场的影响, 这些高能粒子流被捕获进入大气层, 并与大气中存在的原子核发生碰撞产生次级宇宙射线(cascades of secondary cosmic rays), 次级宇宙射线穿过大气层进入土壤再与土壤中的原子核发生碰撞产生能量在1~2 MeV的快中子(fast neutrons)。这些快中子与土壤中的各种原子发生碰撞会失去部分能量而慢化, 成为热中子(能量为0.025 eV)或超热中子(能量1~10 eV), 因此土壤中的快中子一部分经过慢化而被土壤吸收, 而另一部分则会散射到空气中。
快中子被慢化的过程主要受3个方面影响:①物质基本的散射截面或散射能力;②每次碰撞中能量的对数衰减, 这直接说明每次碰撞的有效性;③单位质量中的原子数, 其与物质浓度成正比, 与质量数成反比。研究表明氢原子对快中子具有最高的减缓和慢化能力[ 15], 比如能量在2 MeV的快中子只需要与氢原子发生18次碰撞就能慢化到0.025 eV, 远远小于其他原子的碰撞次数, 碳原子要113次, 氧原子要149次[ 7]。
通常来讲, 扩散到空气中的快中子量与地面含水量呈反比例关系, 地表含水量越大, 含有的氢原子越多, 对快中子的慢化和吸收量越大, 从而在地面探测到的快中子数就越少。相对于土壤中的水分含量, 植被含水量、空气中水分含量[ 16]以及土壤中晶格水[ 17]的含量要相对较小, 虽然近期有学者研究它们对快中子的影响[ 18], 但研究还有待深入, 本文探讨中忽略其影响。快中子强度与土壤含水量的反比例关系是宇宙射线快中子法测量的物理基础, 宇宙射线快中子法通过监测土壤上方快中子的浓度从而实现陆地表面土壤含水量的反演[ 19]。
宇宙射线土壤水分观测系统(COsmic-ray Soil Moisture Observing System, COSMOS)是由美国Hydroinnova公司生产的土壤水分探测装置, 一套宇宙射线土壤水分探测仪器包括:采集器(内置气压传感器、湿度传感器和无线传输模块)、40 Ah电池、快中子探头、热中子探头和机箱, 其中快中子和热中子探头均为白色柱状物, 直径10 cm、高130 ㎝。在进行土壤水分观测时只需要安装快中子探头即可, 本试验初期即只安装一个快中子探头进行土壤水分的观测。
2.2.1探测源区和深度
宇宙射线快中子法测量的源区通常定义为86%仪器探测到的快中子所来源地区, 根据中子传输程序MCNPX[ 20](Monte Carlo N-Particle eXtended)模拟的结果, 在海拔为0 m的地区, 其源区半径约为300 m[ 7], 模拟的结果与早期中子传输理论相一致[ 21, 22]。源区大小主要跟空气的物理化学属性相关, 而与地表的土壤水分含量无关。这主要是因为中子的平均散射自由程与每立方米空气中的分子数呈反比, 空气密度越大, 单位体积空气中的分子数越多, 而中子传输的距离越小, 源区越小。反之, 源区越大。在高海拔地区, 气压较低, 测量源区会比低海拔地区大, 3 000 m海拔地区的源区会比海平面地区的源区大25%左右[ 18]。通常情况下源区半径与气压关系可用下式(根据文献[23]公式21改写)表示:
Rs= Rs[ Ps, 0/ Ps] (1)
式中: Rs, 0为参考气压 Ps, 0下的源区半径(通常使用1个标准气压下的半径值, 为300 m[ 23]); Rs为当前气压 Ps下的半径值。本试验研究区海拔为1 557 m, 经计算得试验探测源区在362 m左右。
宇宙射线快中子法的探测深度 z定义探测器记录的快中子86%来源于0~ z深度土壤中, z与源区内的土壤含水量呈反比。当土壤中的含水量高时, 其含有的氢原子数比较多, 使得土壤中的快中子快速被慢化和吸收而无法继续向下传播。根据MCNPX的模拟结果, 探测深度从12 cm(土壤湿润, 模拟含水量为0.40 m3/m3)增长到76 cm(土壤干燥, 模拟含水量为0)[ 7]。近年来, 有研究[ 24]指出反演深度可用下式计算得到:
Z=5.8/( ρbd× τ+ θ+0.0829) (2)
式中: ρbd为土壤的容重, τ为晶格含水量(lattice water, 定义为晶格水在矿物质和束缚水中的质量比[ 17], 通常比较小, 约为0~0.05, 本文计算中忽略), θ为土壤体积含水量, z为反演深度(单位:cm)。
2.2.2 反演函数
由于快中子传播过程中受气压影响较大, 在反演土壤水分时要去除气压随时间变化带来的影响, 根据仪器使用说明书中提供的公式对测量得到的快中子数进行气压订正, 见公式(3)。
N= Nraw×Exp( β×( P- P0)) (3)
式中: Nraw为原始快中子数, β为常数0.0077, P为测量气压值, P0为当地气压理论值。经过气压订正的快中子数与土壤质量含水量存在非线性关系, 本文使用仪器手册中的土壤含水量计算公式, 如下式:
![]() | (4) |
式中: θm为质量含水量, N为订正后快中子数, N0为干燥条件下的快中子数, 需要根据源区实测土壤含水量进行率定, a1=0.079, a2=0.64, a3=0.37和 a4=0.91为常数项。
研究区位于甘肃黑河流域中游张掖地区大满农业灌区[ 11], COSMOS仪器架设位置为100.37°E, 38.86°N。该地区处于温带干旱半干旱气候区, 多年平均降水量为121.5 mm。大满灌区属河西走廊的低平原部分, 为荒漠中的一片绿洲, 主要是由人工灌溉的农田组成。2012年5~9月试验开展期宇宙射线源区内农作物主要为制种玉米, 源区共包含五星村3个社区农田, 即四社农田(占源区面积51.6%)、五社农田(占源区面积42.6%)和六社农田(占源区面积0.04%), 源区东南方向分布有六社村落。图1为研究区示意图。试验中宇宙射线土壤水分探测仪器COSMOS架设在大满超级站西侧, 仪器装有一台快中子探测器, 探测器底部距离地面约50 cm, 采集器内置气压传感器, 采样频率为1小时。
为计算地表土壤含水量, 需要对宇宙射线仪器进行率定以确定公式(4)中 N0, 2012年7月11日下午和9月16日下午分别在COSMOS仪器东、南、西、北的50, 100, 200和300 m距离上进行土壤采样, 其中每个点分别在膜下和膜间各取一个样本。用称重法获取源区0~10 cm和10~20 cm深度上的土壤含水量, 并在源区内采集4个土壤容重样本, 测量土壤容重数据主要用于将宇宙射线测量得到的质量含水量转化为体积含水量。
对采集的土样进行烘干称重, 得到每个点位0~10 cm和10~20 cm的平均质量含水量, 通过对比发现, 膜下土壤水分和膜间土壤水分相差不大, 标准差(SD)比较小。采样点均匀分布于仪器源区内, 因此对所有测点进行算术平均得到0~10 cm和10~20 cm深度的区域土壤质量含水量[ 23, 25]。对采集的土壤容重样本进行烘干分析得到平均容重为1.49, 测量数据见表1。
![]() | 表1 土样采样数据 Table 1 data of soil samples |
根据测量得到的源区内土壤样本数据和公式(2), 宇宙射线的有效代表深度介于10~20 cm之间, 设定土样10~20 cm土壤水分线性分布, 根据公式(5)计算10+ x cm各深度土壤水分。
![]() | (5) |
式中: SWC0~10, SWC10~20分别为0~10 cm和10~20 cm平均土壤水分, x取值范围为1~10。
计算结果表明7月11日根据称重法测量的17 cm土壤含水量反求的宇宙射线对应深度与设定值最接近, 9月16日为18 cm深度值, 因此率定时分别选用7月11日17 cm深度土壤含水量和9月16日18 cm深度土壤含水量作为 θm。仪器探测到的快中子数进行气压订正后, 分别选取对应土壤样本采集时间(15~20时)的快中子订正数据进行平均, 得到对应土壤水分的快中子数 N, 最后根据公式(3)反求 N0, 经计算7月11日 N0为3 195, 9月16日 N0为3 364, 在计算土壤水分时选用两者的平均值3 280作为最终的 N0。
根据率定后 N0和公式(2)可以计算出6月1日至9月30日玉米生长季内的区域土壤水分变化, 如图2所示, 其中COSMOS绿色线表示未平滑处理前的土壤水分数据, 由于快中子传播过程中存在一定的随机性, 因此宇宙射线仪器记录的快中子数据存在一定波动性, 其测量得到的土壤水分数据在一小时数据中会出现较小范围的波动。分析中常使用12小时平滑后的数据, 图2中蓝色线COSMOS_12为经过12小时平滑后的数。
通过上图可以看出宇宙射线反演的土壤水分主要受灌溉和降水的影响, 并且在灌溉期源区内土壤含水量呈现双峰变化。其主要原因是源区内3个社的灌溉顺序不同, 五星四社农田(如图1所示)先进行灌溉, 灌溉完毕后出现第一个峰值;然后是五星六社农田灌溉, 由于五星六社田块占源区面积比重很小, 灌溉期约为1天, 源区内田块总含水量减少, 对快中子的慢化作用降低, 使探头测量到的快中子数有所上升, 从而使反演的土壤水分降低;在六社灌溉完毕, 五社开始灌溉时, 由于COSMOS源区42.6%位于五社农田内, 因此在其灌溉时源区内土壤含水量会明显升高, 此时出现第二个峰值。
灌溉期源区内土壤水分呈现双峰变化在前3轮灌溉(第一轮6月1~5日;第二轮灌溉7月1~5日, 第三轮灌溉7月27日至8月1日)过程中反应比较明显, 而第四轮灌溉(8月20~28日)时, 由于水流小, 灌溉持续时间较长, 且四社灌溉完后, 五社和六社同时灌溉, 并没有出现中间的低谷区, 双峰变化不明显。
从时间序列上, 源区内土壤水分可分为三个阶段, 6月26日以前土壤水分约为0.25 m3/m3, 在6月26日至9月1日期间, 平均土壤水分在0.35 m3/m3, 在9月1日之后, 土壤含水量持续减少, 约在0.25 m3/m3左右。第一阶段属于封垄前, 玉米植株较矮, 在强日照条件下, 地表蒸发比较大, 土壤水分含量减少很快, 整体含水量比较低;第二阶段为封垄后, 玉米植株高大茂盛, 减弱了地表土壤水分的蒸发, 虽然植物蒸腾量增大, 但考虑到该阶段内灌溉间隔短、降水较多, 土壤总体储水量大于其他两个阶段;第三阶段为玉米生长后期, 玉米接近成熟并逐渐枯萎, 地表蒸发较大, 期间无灌溉补充水分, 土壤含水量减少较快, 整体含水量较低。
在宇宙射线源区范围内共有19个SoilNET站点, 选取数据质量较好的6月22日-8月18日数据对宇宙射线观测值进行验证。由于在该时间段内宇宙射线反演的土壤水分平均值为0.31 m3/m3, 其探测深度在14 cm左右, 因此选取4 cm和10 cm数据进行算术平均作为验证值。图3表示6月22日-8月18日宇宙射线土壤水分和SoilNET平均土壤水分变化对比图, 可以看出两者的变化趋势比较一致, 相关系数为0.79, 均方根误差RMSE为0.037 m3/m3。
![]() | 图4 非灌溉期COSMOS与SoilNET土壤水分对比分析Fig. 4 comparison of the COSMOS and SoilNET soil moisture in non-irrigation time |
Zreda等[ 18]指出源区内的H原子是影响快中子强度的最重要元素。当地面存在积水时, 源区内H原子含量明显增高, 土壤中产生的快中子在穿透积水层过程中很容易被慢化, 因此宇宙射线探测到的快中子在灌溉期会明显减少, 从而使反演的土壤水分增高(如图3散点图中圆圈内数据)。而SoilNET探头埋设在地下, 并不受积水影响, 从而使得宇宙射线观测的值(包含地面积水)大于无线传感器网络观测的值(不含积水)。
剔除灌溉期(7月1~5日, 7月27日至8月1日)数据后相关系数明显升高, 达到0.86, RMSE降低为0.028 m3/m3 (如图4所示)。通过对比可以得出, 宇宙射线反映的是一种面上的土壤水分变化值, 能够很好的反映出农田下垫面下土壤水分的变化趋势, 但在积水情况下宇宙射线存在一定高估现象。不过, 不论在自然下垫面下, 还是农田下垫面下, 积水现象总是很短暂的, 对整体观测影响有限。
目前, 传统土壤水分测量方式多是有损、点测量方式, 在获取农田区域土壤水分数据时费时、费力, 并且在为遥感像元尺度土壤水分数据做验证时多采用以点代面或者插值的方式进行, 无法达到尺度间的匹配, 在遥感技术迅速发展的情况下, 如何获取遥感像元尺度上的土壤含水量信息, 以利于遥感反演精度的验证和评价, 是急需解决的问题。宇宙射线快中子法作为一种新型的区域土壤含水量测量方式, 其探测源区半径在300 m左右, 并且可以实时获取土壤水分数据。本文分析了张掖绿洲农田中开展的区域土壤含水量观测实验结果, 表明宇宙射线土壤水分观测系统能够有效地获取区域土壤含水量, 不仅可以准确获取长时间序列农田尺度土壤水分, 也将为遥感反演的土壤水分提供真正意义上的面状地面验证数据。
当然, 利用宇宙射线快中子法测量土壤水分还有许多问题要解决。首先, 宇宙射线源区内的氢原子不只存在于土壤水分中, 在晶格水、土壤有机质、植被以及水汽中均存在氢原子, 这些氢原子源的存在势必会对快中子的传播造成影响。Zreda等[ 18]人指出晶格水的存在会使得宇宙射线探测深度降低, 不过在认为晶格水、土壤有机质不变的情况下, 这些影响可以通过一次源区实地土壤水分校正进行消除[ 18];在稀疏植被区, 植被的影响也可以根据实地土壤水分校正进行消除, 但Zreda指出, 在快速生长植被下垫面下, 植被含水量对观测存在一定的影响, 目前还没有完善的方法对此进行校正;水汽的存在会影响宇宙射线的测量源区, MCNPX模拟表明饱和水汽下的探测半径会比干空气下的探测半径下降10%左右。水汽不仅仅影响宇宙射线测量源区, 还会对源区内快中子强度造成影响, Zreda等人建议增加水汽测量相关设备对快中子数据进行水汽校正, 以消除土壤水分反演过程中水汽变化带来的影响[ 18]。此外, 当地表被积水覆盖时, 宇宙射线反演土壤水分会受到很大的影响, 本文研究中指出在灌溉的情况下宇宙射线反演土壤水分较站点测量土壤水分高, 关于积水对宇宙射线快中子技术的影响还有待进一步研究。最后, 根据中子散射传播理论, 不同深度的土壤水分对宇宙射线的影响也不是完全相同, 深度加权也是目前研究的热点[ 24, 25]。
致谢:感谢北京师范大学地理学与遥感科学学院以及中国科学院寒区旱区环境与工程研究所各位老师给予的帮助和指导。
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|